C Rutina 065
C-------------------------------------------------------------AIJ_BIJ---
C
C SUBRUTINA AIJ_BIJ
C
C Incrementa las matrices Aij y Bij en la integracion numerica
C
SUBROUTINE AIJ_BIJ(Aij,Bij,W,SF_U,SF_T,FF,Lk)
C
IMPLICIT NONE
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES EXTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C Aij,Bij : Matrices clasicas de constantes de integracion
REAL*8 Aij(2,4),Bij(2,4)
C Lk : Longitud del elmento
REAL*8 Lk
C SF_U(,),SF_T(,) : Nucleos de la integral de contorno
C SF_U(i,j) : Uij
C SF_T(i,j) : Tij
REAL*8 SF_U(2,2),SF_T(2,2)
C FF() : FUNCIONES DE FORMA
REAL*8 FF(2)
C W : Coficiente de peso de la integral de Gauss
REAL*8 W
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C CALC. MATRICES Aij y Eij
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
Aij(1,1)=Aij(1,1)+(W*SF_T(1,1)*FF(1))*Lk/2.d0
Aij(1,2)=Aij(1,2)+(W*SF_T(1,1)*FF(2))*Lk/2.d0
Aij(1,3)=Aij(1,3)+(W*SF_T(1,2)*FF(1))*Lk/2.d0
Aij(1,4)=Aij(1,4)+(W*SF_T(1,2)*FF(2))*Lk/2.d0
Aij(2,1)=Aij(2,1)+(W*SF_T(2,1)*FF(1))*Lk/2.d0
Aij(2,2)=Aij(2,2)+(W*SF_T(2,1)*FF(2))*Lk/2.d0
Aij(2,3)=Aij(2,3)+(W*SF_T(2,2)*FF(1))*Lk/2.d0
Aij(2,4)=Aij(2,4)+(W*SF_T(2,2)*FF(2))*Lk/2.d0
C
Bij(1,1)=Bij(1,1)+(W*SF_U(1,1)*FF(1))*Lk/2.d0
Bij(1,2)=Bij(1,2)+(W*SF_U(1,1)*FF(2))*Lk/2.d0
Bij(1,3)=Bij(1,3)+(W*SF_U(1,2)*FF(1))*Lk/2.d0
Bij(1,4)=Bij(1,4)+(W*SF_U(1,2)*FF(2))*Lk/2.d0
Bij(2,1)=Bij(2,1)+(W*SF_U(2,1)*FF(1))*Lk/2.d0
Bij(2,2)=Bij(2,2)+(W*SF_U(2,1)*FF(2))*Lk/2.d0
Bij(2,3)=Bij(2,3)+(W*SF_U(2,2)*FF(1))*Lk/2.d0
Bij(2,4)=Bij(2,4)+(W*SF_U(2,2)*FF(2))*Lk/2.d0
C
C FIN DE LA RUTINA
C
RETURN
09210 FORMAT(/,2x,'ERROR 09210 - AIJBIJ - Error Programa (#',i1,'#)')
C
END
C
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
