DESARROLLO DE HERRAMIENTAS - FI-UAEMex

136 HERRAMIENTAS DE CÓMPUTO
c
10 FORMAT(/,1X,'ARCHIVO DE DATOS: ')
READ(*,42)INPUT
42 FORMAT(A20)
WRITE(*,11)
READ(*,113)OUTPUT
113 FORMAT (A20)
11 FORMAT(/,1X,'ARCHIVO DE SALIDA: ')
OPEN(1,FILE=INPUT,STATUS='OLD')
OPEN(2,FILE=OUTPUT,STATUS='unknown')
WRITE(*,*)
WRITE(2,*)'******************************************************'
WRITE(2,*)'* *'
WRITE(2,*)'* ANALISIS DE ARMADURAS PLANAS *'
WRITE(2,*)'* ( A R M A 2 D ) *'
WRITE(2,*)'* *'
WRITE(2,*)'******************************************************'
C
C LECTURA DE DATOS GENERALES
C
READ (1,*)NB,NU,NAP
C
nnu=NU+NAP
nuu=2*NU
nun=nuu+1
C
C LECTUA DE COORDENADAS DE NUDOS Y FUERZAS EN LOS MISMOS
C
DO 200 I=1,NNU
READ(1,*)X(I),Y(I),K(2*i-1,nun),K(2*i,nun)
200 CONTINUE
C
C GENERACION DE LA MATRIZ DE CONTINUIDAD [ A ]
C
DO 250 I=1,NB
C LECTURA DE LA RIGIDEZ AXIAL, EL NUDO INICIAL Y EL NUDO FINAL D LAS BARRAS
C
READ(1,*)E(I),AR(I),IIN(i),IFI(i)
c
L(i)=((x(Ifi(i))-x(iin(i)))**2+(y(Ifi(i))-y(iin(i)))**2)**.5
ux=(x(Ifi(i))-x(iin(i)))/l(i)
uy=(y(Ifi(i))-y(iin(i)))/l(i)
IF (iin(I).lE.NU) THEN
a(i,2*iin(i)-1)=-ux
a(i,2*iin(i))=-uy
ENDIF
2 IF ( Ifi(i).lE.NU) THEN
a(i,2*IFI(i)-1)=ux
a(i,2*IFI(i))=uy
ENDIF
250 CONTINUE
write(2,22)
22 format(//'Matriz de Continuidad [A]'//)
WRITE(2,39)((a(I,J),j=1,2*nu),i=1,Nb)
39 FORMAT(4F10.4)
C [AT][K][A]
C
c nuu : dimension de la matriz de rigideces [K]
C EAL(I)=(E(I)*AR(I))/(L(I))
nuu=2*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(5f10.4)
C
C SOLUCION DEL SISTEMA POR GAUSS-JORDAN
C
DESARROLLO DE HERRAMIENTAS DE ANÁLISIS ESTRUCTURAL
PARA SU USO DESDE LA INTERNET