Iteracion matlab elementos finitos
Publicado por Marimar (4 intervenciones) el 30/03/2019 19:34:29
hola tengo un problema de fractura de elementos finitos (aunque no sepais del tema podeis ayudar). En este tengo dos ecuaciones.
>> z(1) = f_t/y(2)*(h-a_0-y(1))^2*b/2-y(1)*f_t/2*b-f_t/2*max(1-y(1)/y(2)*eps_t/eps_max,0)*(y(1)-y(2))*b;
>>z(2) = f_t/y(2)*(h-a_0-y(1))^3*b/3+f_t*y(2)^2*b/3+f_t*max(1-y(1)/y(2)*eps_t/eps_max,0)*(y(1)^2-y(2)^2)*b/2 +f_t/2*(1-max(1-y(1)/y(2)*eps_t/eps_max,0))*(y(1)^2 +y(1)*y(2)-2*y(2)^2)*b/3-P*L/4;
Estas dos las defino en una funcion y realizo un fsolve (z=0) y hallo mi incognita y y el error cometido (FVAL de fsolve). Mi problema es que para realizar esto debo de probar diferentes valores de P hasta que me salga un error 0. Me preguntaba si habria una manera de resolver el problema para no ir yo manualmente o con un for, variando el valor de P. Es decir, si habria una manera de obtener "y" y "P" minimizando el error.
NOTA: de las ecuaciones todos son valores conocidos menos P e y(1), y(2).
Muchas gracias !!!!!1
>> z(1) = f_t/y(2)*(h-a_0-y(1))^2*b/2-y(1)*f_t/2*b-f_t/2*max(1-y(1)/y(2)*eps_t/eps_max,0)*(y(1)-y(2))*b;
>>z(2) = f_t/y(2)*(h-a_0-y(1))^3*b/3+f_t*y(2)^2*b/3+f_t*max(1-y(1)/y(2)*eps_t/eps_max,0)*(y(1)^2-y(2)^2)*b/2 +f_t/2*(1-max(1-y(1)/y(2)*eps_t/eps_max,0))*(y(1)^2 +y(1)*y(2)-2*y(2)^2)*b/3-P*L/4;
Estas dos las defino en una funcion y realizo un fsolve (z=0) y hallo mi incognita y y el error cometido (FVAL de fsolve). Mi problema es que para realizar esto debo de probar diferentes valores de P hasta que me salga un error 0. Me preguntaba si habria una manera de resolver el problema para no ir yo manualmente o con un for, variando el valor de P. Es decir, si habria una manera de obtener "y" y "P" minimizando el error.
NOTA: de las ecuaciones todos son valores conocidos menos P e y(1), y(2).
Muchas gracias !!!!!1
Valora esta pregunta
0
