% 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=50; %input('Ingrese el Ancho Total de la Viga en [cm]: ');
b=100; %input('Ingrese el Alto Total de la Viga en [cm]: ');
t1=14; %input('Ingrese el espesor del ala en [cm]: ');
t2=25; %input('Ingrese el espesor del alma en [cm]: ');
num_barras=3; %input('Ingrese numero de barras de acero: ');
diametro=26; %input('Ingrese el diametro de las barras en [mm]: ');
d=97.5; %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; %segun bibliografia
yt=b-y1;
Mcr1=fr*I1/yt;
Curv1=Mcr1/(Ec*I1);
%display(Mcr1);
%display(Curv1);
 
%despues de la rotura (igualacion areas) T=C
tolerancia=9999;
c=1;

while tolerancia > 150
    aux2=a*t1*(c-0.5*t1)+0.5*t2*(c-t1)^2;
    aux3=n*As*(d-c);
    tolerancia=abs(aux2-aux3);
    c=c+0.1;
end
%display(c);
I2=(((a/12)*t1^3)+(a*t1*(c-0.5*t1)^2))+(((t2/12)*(c-t1)^3)+(0.5*t2*(c-t1)^3))+(n*As*(d-c)^2);
% falta calcular Mcr2 y Curv 2 ¿?¿?¿


%Empieza a ser no lineal en
fc=0.4*fc_prima; %segun bibliografia
e_cm=fc/Ec;
Curv3=e_cm/c;
e_s=e_cm*(d-c)/c;
Mcr3=Curv3*Ec*I2;


%hormigon no lineal

e_0=2*fc_prima/Ec;
if e_0 >= e_cm
    fc=fc_prima*((/)-(/)^2);
else

end