C Rutina 066
C----------------------------------------------------------PF_AIJBIJ2---
C
C SUBRUTINA PF_AIJBIJ2
C
C PROBLEMA 2-D Plano
C Calcula Aij y Bij analiticamente
C
SUBROUTINE PF_AIJBIJ2(Aij,Bij,PMAT,Lk,TP,N1,N2,COD)
C
IMPLICIT NONE
C
INCLUDE '_param.for'
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 PMAT() : Propiedades del material
REAL*8 PMAT(2)
C Lk : Longitud del elmento y semilongitudes (sing. en el int.)
REAL*8 Lk
C TP : Tipo de problema
INTEGER TP
C N1,N2 : Normal al elemento
REAL*8 N1,N2
C COD = 1 : P.C. no pertenece al elemento
C 2 : P.C. en el origen del elemento
C 3 : P.C. en el interior del elemento
C 4 : P.C. en el extremo del elemento
INTEGER COD
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES INTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C nu,PI,G
REAL*8 nu,PI,G
C FCsAa,FCsAb,FCsBa,Ca,Cb1,Cb2 : Factor comun del termino singular
REAL*8 Ca,C1b,C2b
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C CALCULO DE Aij y Bij
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
Pi=DACOS(-1.d0)
C
SELECT CASE (TP)
C tension plana
CASE(21)
G=PMAT(1)/2.d0/(1+PMAT(2)) !EG000605
nu=PMAT(2)/(1.d0+PMAT(2))
c deformacion plana
CASE(22)
G=PMAT(1)/2.d0/(1+PMAT(2)) !EG000605
nu=PMAT(2)
CASE DEFAULT
WRITE(F_W,09230) 3
WRITE (*,09230) 3
STOP
END SELECT
C
SELECT CASE (COD)
C
C Singularidad en el origen
C
CASE(2)
C1b=(3.d0-4.d0*nu)/8.d0/Pi/(1.d0-nu)/G !EG000605
C2b= 1.d0/8.d0/Pi/(1.d0-nu)/G !EG000605
Ca =-(1.d0-2.d0*nu)/4.d0/Pi/(1.d0-nu)
C
Bij(1,1)=C2b*Lk/2.d0*n2**2.d0-C1b*Lk/2.d0*(DLOG(Lk)-3.d0/2.d0)
Bij(1,2)=C2b*Lk/2.d0*n2**2.d0-C1b*Lk/2.d0*(DLOG(Lk)-1.d0/2.d0)
Bij(1,3)=-C2b*Lk/2.d0*n1*n2
Bij(1,4)=Bij(1,3)
Bij(2,1)=Bij(1,4)
Bij(2,2)=Bij(1,3)
Bij(2,3)=C2b*Lk/2.d0*n1**2.d0-C1b*Lk/2.d0*(DLOG(Lk)-3.d0/2.d0)
Bij(2,4)=C2b*Lk/2.d0*n1**2.d0-C1b*Lk/2.d0*(DLOG(Lk)-1.d0/2.d0)
C
Aij(1,3)=-Ca*(1.d0-DLOG(Lk))
Aij(1,4)= Ca
Aij(2,1)= Ca*(1.d0-DLOG(Lk))
Aij(2,2)=-Ca
C
C Singularidad en el extremo
C
CASE(4)
C1b=(3.d0-4.d0*nu)/8.d0/Pi/(1.d0-nu)/G !EG000605
C2b= 1.d0/8.d0/Pi/(1.d0-nu)/G !EG000605
Ca=-(1.d0-2.d0*nu)/4.d0/Pi/(1.d0-nu)
C
Bij(1,1)=C2b*Lk/2.d0*n2**2.d0-C1b*Lk/2.d0*(DLOG(Lk)-1.d0/2.d0)
Bij(1,2)=C2b*Lk/2.d0*n2**2.d0-C1b*Lk/2.d0*(DLOG(Lk)-3.d0/2.d0)
Bij(1,3)=-C2b*Lk/2.d0*n1*n2
Bij(1,4)=Bij(1,3)
Bij(2,1)=Bij(1,4)
Bij(2,2)=Bij(1,3)
Bij(2,3)=C2b*Lk/2.d0*n1**2.d0-C1b*Lk/2.d0*(DLOG(Lk)-1.d0/2.d0)
Bij(2,4)=C2b*Lk/2.d0*n1**2.d0-C1b*Lk/2.d0*(DLOG(Lk)-3.d0/2.d0)
Aij(1,3)=-Ca
Aij(1,4)= Ca*(1.d0-DLOG(Lk))
Aij(2,1)= Ca
Aij(2,2)=-Ca*(1.d0-DLOG(Lk))
C
C
C Singularidad en el interior
C
CASE(3)
WRITE(F_W,09230) 1
WRITE (*,09230) 1
STOP
CASE DEFAULT
WRITE(F_W,09230) 2
WRITE (*,09230) 2
STOP
END SELECT
C
C FIN DE LA RUTINA
C
RETURN
09230 FORMAT(/,2x,'ERROR 09230 - PF_AIJBIJ2 - Error Programa (#',i1,
@ '#)')
C
END
C
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
