análise de pavimentos de edifícios de concreto armado com a ...

C
C Determinacao do escoamento ( My,(1/r)y )
alfae = Es/Ec
xis(1) = (-alfae*(As+Als)+(SQRT(alfae*(alfae*((As+Als)*(As+Als))+
* 2*b*(Als*dl+As*d)))))/b
ecc = fyk*As/( (Ec*xis(1)*b/2.)+(Als*Es*(xis(1)-dl)/xis(1)) )
scc = ecc*Ec
IF (scc.GT.fcm) THEN
scc = fcm
ecc = scc/Ec
END IF
esc = (ecc*(xis(1)-dl))/xis(1)
Rcc = (scc*xis(1)*b)/2.
Rsc = esc*Es*Als
M(2) = Rcc*(d-xis(1)/3.)+Rsc*(d-dl)
est = ecc*(d-xis(1))/xis(1)
curv(2) = est/(d-xis(1))
C
C Determinacao da situação última ( Mu,(1/r)u )
aa = 0.8*b*fcm
bb = -((0.8*b*fcm*d)+(0.01*Als*Es)+As*fyk)
cc = (0.01*Als*Es*dl)+As*fyk*d
xis(2) = (-bb-SQRT((bb*bb)-(4*aa*cc)))/(2.*aa)
ecc = 0.01*(xis(2))/(d-xis(2))
IF (ecc.LT.0.02)THEN
aa = Ec*b
bb = (2*Es*Als)+(200*fyk*As)
cc = -( (2*Es*Als*dl)+(200*fyk*As*d) )
xis(2) = (-bb+SQRT((bb*bb)-(4*aa*cc)))/(2.*aa)
ecc = 0.01*(xis(2))/(d-xis(2))
END IF
esc = 0.01*(xis(2)-dl)/(d-xis(2))
curv(3) = 0.010/(d-xis(2))
IF (ecc.GT.0.0035)THEN
ecc = 0.0035
esc = 0.0035*(xis(2)-dl)/xis(2)
curv(3) = 0.0035/xis(2)
END IF
Rcc = 0.8*xis(2)*b*fcm
Rsc = esc*Es*Als
M(3) = Rcc*(d-(0.4*xis(2)))+Rsc*(d-dl)
IF (ecc.LT.0.02)THEN
M(3) = ( As*fyk*(d-xis(2)/3.) )-( Rsc*(dl-xis(2)) )
END IF
C
C Determinacao dos parametros de encruamento
Dr = M(1)/curv(1)
Dy = (M(2)-M(1))/(curv(2)-curv(1))
Du = (M(3)-M(2))/(curv(3)-curv(2))
ky = Dy/(1-Dy/Dr)
ku = Du/(1-Du/Dr)
C
C Gravação dos resultados
WRITE (5,10)
10 FORMAT (/3X,'********************************',/5X,
* 'DIAGRAMA MOMENTO x CURVATURA',/3X,'*********************************',/5X,
* 'Mr',8X,'(1/r)r',7X,'My',8X,'(1/r)y',7X,'Mu',7X,'(1/r)u',/)
WRITE (5,15)M(1),curv(1),M(2),curv(2),M(3),curv(3)
15 FORMAT(3(F9.2,2X,E10.2,2X))