% 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;
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);