DESARROLLO DE HERRAMIENTAS - FI-UAEMex

uz=(z(Ifi(i))-z(iin(i)))/l(i)
IF (iin(I).lE.NU) THEN
a(i,3*iin(i)-2)=-ux
a(i,3*iin(i)-1)=-uy
a(i,3*iin(i))=-uz
ENDIF
2 IF ( Ifi(i).lE.NU) THEN
a(i,3*IFI(i)-2)=ux
a(i,3*IFI(i)-1)=uy
a(i,3*IFI(i))=uz
ENDIF
250 CONTINUE
write(2,22)
22 format(//'Matriz de Continuidad [A]'//)
WRITE(2,39)((a(I,J),j=1,3*nu),i=1,Nb)
39 FORMAT(6F10.4)
C [AT][K][A]
C
c nuu : dimension de la matriz de rigideces [K]
nuu=3*nu
DO 260 I=1,NUU
DO 280 J=1,NUU
DO 300 M=1,NB
EAL(M)=(E(M)*AR(M))/(L(M))
K(i,j)=K(i,j)+a(M,i)*a(M,j)*EAL(M)
300 CONTINUE
280 CONTINUE
260 CONTINUE
write(2,23)
23 format(//'Matriz de Rigideces [K]'//)
WRITE(2,37)((K(I,J),j=1,NUN),i=1,NUU)
37 format(7f12.2)
C
C SOLUCION DEL SISTEMA POR GAUSS-JORDAN
C
N=NUU
c
DO 146 ZZ=1,N
DO 144 I=1,N
DO 150 J=N+1,ZZ,-1
IF (I.EQ.ZZ) GOTO 144
IF (K(ZZ,ZZ).EQ.0) THEN
C
DO 132 Q=ZZ+1,N
IF (K(Q,ZZ).NE.0) THEN
DO 134 QQ=1,N+1
W=K(Q,QQ)
K(Q,QQ)=K(ZZ,QQ)
K(ZZ,QQ)=W
134 CONTINUE
GOTO 142
ENDIF
132 CONTINUE
C
WRITE(*,*)'EL SISTEMA ES INDETERMINADO'
STOP
ENDIF
C
142 K(I,J)=K(I,J)+K(ZZ,J)*(-K(I,ZZ))/K(ZZ,ZZ)
C
150 CONTINUE
144 CONTINUE
146 CONTINUE
C
DO 128 I=1,N
K(I,N+1)=K(I,N+1)/K(I,I)
128 CONTINUE
C
C IMPRIME LOS DESPLAZAMIENTOS DE LOS NUDOS
C
write(2,47)
47 format(//'Desplazamientos de los nudos :'//)
DO 600 I=1,NU
WRITE(2,*)I,'DX',K(3*I-2,N+1)
WRITE(2,*)I,'DY',K(3*I-1,N+1)
HERRAMIENTAS DE CÓMPUTO 139
DESARROLLO DE HERRAMIENTAS DE ANÁLISIS ESTRUCTURAL
PARA SU USO DESDE LA INTERNET