DESARROLLO DE HERRAMIENTAS - FI-UAEMex

HERRAMIENTAS DE CÓMPUTO 163
c ******************************************************************
c GENERACION DE LA MATRIZ DE CONTINUIDAD [ A ] DE CADA BARRA Y
c ENSAMBLE DE SU PARTICIPACION A LA MATRIZ DE RIGIDEZ GLOBAL
c
c Barrido de elementos
c
DO 1250 IB=1,NB
Write(2,782)ib
782 format(/' barra ',i5/)
C
C Lectura de propiedades ( E,I,G,J) y conectividades ( A,B ) de la barra
C
READ(1,*)E,IY,G,JJ,IIN,IFI
C
c Longitud y cosenos directores de la barra
c
L=((x(ifi)-x(iin))**2+(y(ifi)-y(iin))**2)**.5
ux=(x(ifi)-x(iin))/L
uy=(y(ifi)-y(iin))/L
c
C Matriz de rigidez diagonal de la barra
C
R(1)=2*E*Iy/(L)
R(2)=2*E*Iy/(L)
R(3)=2*E*Iy/(L)
R(4)=G*JJ/L
c
c El extremo A de la barra, es nudo !
c
IF(IIN.LE.NU)THEN
a(1,1)=-uy
a(1,2)=ux
a(1,3)=-1/L
a(2,1)=-uy
a(2,2)=ux
a(2,3)=-2/L
a(3,1)=0
a(3,2)=0
a(3,3)=-1/L
a(4,1)=-ux
a(4,2)=-uy
a(4,3)=0
C Identifica los desplazamientos en el nudo A de la barra
DO 1650 I=1,3
DAB(I)=K(3*IIN-3+I,NUN)
1650 CONTINUE
ENDIF
c
c El extremo B de la barra, es nudo !
c
IF(IFI.LE.NU)THEN
a(1,4)=0
a(1,5)=0
a(1,6)=1/L
a(2,4)=-uy
a(2,5)=ux
a(2,6)=2/L
a(3,4)=-uy
a(3,5)=ux
a(3,6)=1/L
a(4,4)=ux
a(4,5)=uy
a(4,6)=0
C Identifica los desplazamientos en el nudo B de la barra
DO 1550 I=4,6
DAB(I)=K(3*IFI-6+I,NUN)
1550 CONTINUE
ENDIF
c
WRITE(2,192)(DAB(j),j=1,6)
192 FORMAT(' DESPLAZAMIENTOS EN A Y B ', 6F10.4)
c Producto {e}= [A]{d}
c
DO 8000 I=1,4
DO 8001 J=1,6
DESARROLLO DE HERRAMIENTAS DE ANÁLISIS ESTRUCTURAL
PARA SU USO DESDE LA INTERNET