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

C TES - matriz que armazena os pivots auxiliares para a condensação de SF nos
resíduos
SUBROUTINE COND(NTF,NEF,PIVOTS,TES,ASALL)
IMPLICIT REAL*4 (A-H,O-Z)
INTEGER I,J,K,K1,L,N,NIF,NTF,NEF
REAL*4 PIVOTS(14,3),TES(12,3),ASALL(3)
REAL*8 T,C
COMMON/EM/ LM(24),ND,NS,PF(24),RF(24,4),XM(24),SA(12,24),SF(12,4),
* ASA(24,24),LMI(4),IX(7),IY(7),EL(4),TLO(5,4),QQQ(696)
NIF = NTF - NEF
DO 200 N = 1,NIF
K = NTF - N
L = K + 1
IF (ASA(L,L).LE.0.)GO TO 200
DO 165 I = 1,K
IF (ASA(L,I).EQ.0.) GO TO 165
T = ASA(L,I)/ASA(L,L)
ASALL(N)=ASA(L,L)
C
C Matriz de rigidez [ASA]
DO 160 J = I,K
ASA(J,I) = ASA(J,I) - ASA(J,L)*T
ASA(I,J) = ASA(J,I)
160 END DO
165 END DO
C=RF(L,1)/ASA(L,L)
C
C Vetor de cargas [RF] com uma coluna
DO 180 I=1,K
PIVOTS(I,N)=ASA(I,L)/ASA(L,L)
RF(I,1)=RF(I,1)-ASA(I,L)*C
180 CONTINUE
C
C Matriz Momento x Deslocamento [SA] e Momentos Iniciais [SF]
IF(ASA(L,L).LT.1.0D-8) GO TO 200
DO 195 K1=1,12
T=Sa(K1,L)
TES(K1,N)=T
Sa(K1,L)=Sa(K1,L)/ASA(L,L)
DO 190 I=1,K
Sa(K1,I) = Sa(K1,I) - (ASA(I,L)/ASA(L,L))*T
190 CONTINUE
SF(K1,1)=SF(K1,1)+T*C
195 CONTINUE
200 CONTINUE
RETURN
END
Apêndice 2 - Subrotina para a determinação da raiz de f( γ ) = 0