consulta con ecuaciones y variables dependientes de otra funcion
Publicado por Camilo García (52 intervenciones) el 26/10/2011 06:20:05
Estimados señores:
Junto con saludarlos, les presento la siguiente duda:
tengo que encontrar la solucion a una ecuacion, pero la variable a encontrar depende de otra funcion.
En primera instancia se me ocurrio trabajar la variable como una syms, pero el programa me deja un error.
les dejo el codigo para ver si pueden ayudarme en mi problema.
Si les sirve de algo tengo la resolución de este problema en mathcadV14. Si les sirve de apoyo dejenme su correo y se los envio.
les dejo el codigo que intente realizar a ver si me pueden apoyar. De antemano muchas gracias, Saludos
programa principal
%% problema con flexión pura.
b=30;
h=30;
fc=300;
fy=4200;
Es=2100000;
r=4;
d=[r
h-r];
as=[2
2];
P1=0;
syms cfp
for i=1:2
es(i,1)=findes(cfp,d(i,1));
fs(i,1)=findfs(findes(cfp,d(1,1),fy,Es);
end
cf=fsolve(0.85*fc*cfp*b+as(1,1)*(h/2-d(1,1))*fs(1,1)+as(2,1)*(h/2-d(2,1))*fs(2,1)-P1,1)
funcion findes
function [es] = findes(cc,dd)
es = 0.003*(cc-dd)/(cc);
end
funcion findfs
function [es] = findfs(ep,fy,Es)
if ep>=fy/Es
es=fy
elseif ep>=-fy/Es && ep<=fy/Es
es=ep*Es
elseif ep<=-fy/Es
es=-fy
end
%-----------------------------------------------------------------------
Junto con saludarlos, les presento la siguiente duda:
tengo que encontrar la solucion a una ecuacion, pero la variable a encontrar depende de otra funcion.
En primera instancia se me ocurrio trabajar la variable como una syms, pero el programa me deja un error.
les dejo el codigo para ver si pueden ayudarme en mi problema.
Si les sirve de algo tengo la resolución de este problema en mathcadV14. Si les sirve de apoyo dejenme su correo y se los envio.
les dejo el codigo que intente realizar a ver si me pueden apoyar. De antemano muchas gracias, Saludos
programa principal
%% problema con flexión pura.
b=30;
h=30;
fc=300;
fy=4200;
Es=2100000;
r=4;
d=[r
h-r];
as=[2
2];
P1=0;
syms cfp
for i=1:2
es(i,1)=findes(cfp,d(i,1));
fs(i,1)=findfs(findes(cfp,d(1,1),fy,Es);
end
cf=fsolve(0.85*fc*cfp*b+as(1,1)*(h/2-d(1,1))*fs(1,1)+as(2,1)*(h/2-d(2,1))*fs(2,1)-P1,1)
funcion findes
function [es] = findes(cc,dd)
es = 0.003*(cc-dd)/(cc);
end
funcion findfs
function [es] = findfs(ep,fy,Es)
if ep>=fy/Es
es=fy
elseif ep>=-fy/Es && ep<=fy/Es
es=ep*Es
elseif ep<=-fy/Es
es=-fy
end
%-----------------------------------------------------------------------
Valora esta pregunta
0