Untitled

% momento v/s curvatura, viga H.A. "T" simple

fc_prima=200; %kgf/cm2
Fy=2800; %tipo de acero, A44-28H
Es=2100000;  %modulo elasticidad acero
Ec=15100*(fc_prima)^0.5;
n=Es/Ec;
a=input('Ingrese el Ancho Total de la Viga en [cm]: ');
b=input('Ingrese el Alto Total de la Viga en [cm]: ');
t1=input('Ingrese el espesor del ala en [cm]: ');
t2=input('Ingrese el espesor del alma en [cm]: ');
num_barras=input('Ingrese numero de barras de acero: ');
diametro=input('Ingrese el diametro de las barras en [mm]: ');
d=input('Ingrese distancia (altura) del centro del fierro a cota superior de la Viga [cm]: ');
while d >= b
    d=input('Ingrese distancia (altura) del centro del fierro a cota superior de la Viga [cm]: ');
end
As=num_barras*3.14*(0.5*diametro/10)^2; % en cm2
aux1=b-t1;
Aeq=a*t1+t2*aux1+(n-1)*As;
%calculo centroide seccion transformada
y1=(0.5*a*t1*t1+aux1*t2*(0.5*aux1+t1)+(n-1)*As*d)/(a*t1+aux1*t2+(n-1)*As);
%calculo inercia seccion transformada
I1=((a*t1^3)/12)+(a*t1*(0.5*t1-y1)^2)+((t2*aux1^3)/12)+(t2*aux1*(0.5*aux1-y1)^2)+(n-1)*As*(d-y1)^2;

%calculo rotura
fr=2*fc_prima^0.5;
yt=b-y1;
Mcr=fr*I1/yt;
Curvatura=Mcr/(Ec*I1);
display(Mcr);
display(Curvatura);