DESARROLLO DE HERRAMIENTAS - FI-UAEMex

142 HERRAMIENTAS DE CÓMPUTO
R(1)=2*E*IZ/(L)
R(2)=2*E*IZ/(L)
R(3)=2*E*IZ/(L)
R(4)=E*AA/L
c
c El extremo A de la barra, es nudo !
c
IF(IIN.LE.NU)THEN
a(1,1)=-uy/L
a(1,2)=ux/L
a(1,3)=1
a(2,1)=-(2*uy)/L
a(2,2)=(2*ux)/L
a(2,3)=1
a(3,1)=-uy/L
a(3,2)=ux/L
a(3,3)=0
a(4,1)=-ux
a(4,2)=-uy
a(4,3)=0
c
c Producto [AT][K][A]
c
DO 565 I=1,3
DO 585 J=1,3
DO 505 M=1,4
K(3*iin-3+i,3*iin-3+j)=K(3*iin-3+i,3*iin-3+j)+a(M,i)*a(M,j)*R(M)
505 CONTINUE
585 CONTINUE
565 CONTINUE
ENDIF
c
c El extremo B de la barra, es nudo !
c
IF(IFI.LE.NU)THEN
a(1,4)=uy/L
a(1,5)=-ux/L
a(1,6)=0
a(2,4)=(2*uy)/L
a(2,5)=-(2*ux)/L
a(2,6)=1
a(3,4)=uy/L
a(3,5)=-ux/L
a(3,6)=1
a(4,4)=ux
a(4,5)=uy
a(4,6)=0
c
c Producto [AT][K][A]
c
DO 555 I=4,6
DO 655 J=4,6
DO 755 M=1,4
K(3*ifi-6+i,3*ifi-6+j)=K(3*ifi-6+i,3*ifi-6+j)+a(M,i)*a(M,j)*R(M)
755 CONTINUE
655 CONTINUE
555 CONTINUE
ENDIF
c
c A y B son nudos
c Producto [AT][K][A]
IF((IIN.LE.NU).and.(ifi.le.nu))THEN
DO 515 I=1,3
DO 615 J=4,6
DO 715 M=1,4
K(3*iin-3+i,3*ifi-6+j)=K(3*iin-3+i,3*ifi-6+j)+a(M,i)*a(M,j)*R(M)
715 CONTINUE
615 CONTINUE
515 CONTINUE
c
DO 2515 I=1,3
DO 2615 J=4,6
DO 2715 M=1,4
K(3*ifi-6+j,3*iin-3+i)=K(3*ifi-6+j,3*iin-3+i)+a(M,i)*a(M,j)*R(M)
2715 CONTINUE
DESARROLLO DE HERRAMIENTAS DE ANÁLISIS ESTRUCTURAL
PARA SU USO DESDE LA INTERNET