C Rutina 052
C-------------------------------------------------------------CIJ_1---
C
C SUBRUTINA CIJ_1
C
C CALCULA EL TERMINO LIBRE (ELASTICIDAD AXISIMETRICA)
C
SUBROUTINE CIJ_1(CIJ,NU,Na,Nd,T,COD)
C
IMPLICIT NONE
INCLUDE '_param.for'
INCLUDE '_struct.for'
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES EXTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C CIJ(,) : Matriz del termino libre
REAL*8 CIJ(MAX_DIM,MAX_DIM)
C Na,Nd : Normal correspondiente a los elementos anterior y posterior
REAL*8 Na(2),Nd(2)
C T : Angulo interior en el nodo
REAL*8 T
C NU : Coef. Poisson
REAL*8 NU
C SIG : Variable auxiliar
REAL*8 SIG
C COD : COD=1,4 : Nodo normal
C COD=2 : Nodo perteneciente al eje (R=0). Nodo inicial.
C COD=3 : Nodo perteneciente al eje (R=0). Nodo final.
INTEGER COD
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES INTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C A,AS,PI : Variables auxiliares
REAL*8 A,AS,PI
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VALORES INICIALES
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
PI=DACOS(-1.d0)
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C CALCULO TERMINO LIBRE
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
A =1.d0/(4.d0*Pi*(1.d0-NU))
AS=1.d0/(4.d0*(1.d0-NU))
C
SELECT CASE(COD)
CASE(1,4)
CIJ(1,1)=T/2.d0/Pi + A*(Nd(1)*Nd(2)-Na(1)*Na(2))
CIJ(2,1)= A*(Nd(2)**2.d0-Na(2)**2.d0)
CIJ(1,2)= A*(Nd(2)**2.d0-Na(2)**2.d0)
CIJ(2,2)=T/2.d0/Pi + A*(Na(1)*Na(2)-Nd(1)*Nd(2))
CASE(2)
SIG=DSIGN(1.d0,-Nd(2))
CIJ(1,1)= 0.d0
CIJ(1,2)= 0.d0
CIJ(2,1)=SIG*(-AS*Nd(2)**3.d0)
CIJ(2,2)=SIG*(1.d0/2.d0 - Nd(1)/2.d0 + AS*Nd(1)*Nd(2)**2.d0)
CASE(3)
SIG=DSIGN(1.d0,Nd(2))
CIJ(1,1)= 0.d0
CIJ(1,2)= 0.d0
CIJ(2,1)=SIG*(-AS*Na(2)**3.d0)
CIJ(2,2)=SIG*( 1.d0/2.d0 - Na(1)/2.d0 + AS*Na(1)*Na(2)**2.d0)
CASE DEFAULT
WRITE(F_W,08110) 1
WRITE (*,08110) 1
STOP
END SELECT
C
C
C FIN DE LA RUTINA
C
RETURN
C
C FORMATOS DE ESCRITURA
C
08110 FORMAT(/,2x,'ERROR 08110 - CIJ_1 - Error Programa (#',i1,'#)')
C
END
C
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
