DESARROLLO DE HERRAMIENTAS - FI-UAEMex

HERRAMIENTAS DE CÓMPUTO 143
2615 CONTINUE
2515 CONTINUE
C
ENDIF
c Imprime A
write(2,22)
22 format('Matriz de Continuidad [A]'/)
WRITE(2,39)((a(i,j),j=1,6),i=1,4)
39 FORMAT(6F10.4)
c
C Termina el ciclo del barrido y limpia la matriz de continuidad [ A ]
DO 320 I=1,4
DO 340 J=1,3*nu
a(I,J)=0
340 CONTINUE
320 CONTINUE
C
250 CONTINUE
c
c Hace simetrica la matriz [ K ]
c
DO 267 I=1,6
DO 287 J=1,6
K(j,i)=K(i,j)
287 CONTINUE
267 CONTINUE
c
c Impresion de la matriz de rigidez global [ K ]
c
write(2,23)
23 format(//'Matriz Golbal de Rigideces [ K ]'//)
WRITE(2,37)((K(I,J),j=1,NUU),i=1,NUU)
37 format(6f10.4)
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
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)
DESARROLLO DE HERRAMIENTAS DE ANÁLISIS ESTRUCTURAL
PARA SU USO DESDE LA INTERNET