C Rutina 052
C-------------------------------------------------------------CIJ_2---
C
C SUBRUTINA CIJ_2
C
C PROBLEMA DE ELASTICIDAD PLANA
C CALCULO DEL TERMINO LIBRE
C
SUBROUTINE CIJ_2(C,POIS,Na,Nd,T,TP)
C
IMPLICIT NONE
INCLUDE '_param.for'
INCLUDE '_struct.for'
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES EXTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C C(,) : Termino libre
REAL*8 C(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 POIS : Coef. Poisson 3-D
REAL*8 POIS
C TP : Tipo de problema
INTEGER TP
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES INTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C NU : Coef. Poisson 2-D
REAL*8 NU
C PI : Numero Pi
REAL*8 PI
REAL*8 A
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VALORES INICIALES
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
PI=DACOS(-1.d0)
C
SELECT CASE (TP)
C Tension plana
CASE(21)
nu=POIS/(1.d0+POIS)
C Deformacion plana
CASE(22)
nu=POIS
CASE DEFAULT
WRITE(F_W,08120) 1
WRITE (*,08120) 1
STOP
END SELECT
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C CALCULO TERMINO LIBRE
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
A=1.d0/(4.d0*Pi*(1.d0-NU))
C(1,1)=T/2.d0/Pi + A*(Nd(1)*Nd(2)-Na(1)*Na(2))
C(2,1)= A*(Nd(2)**2.d0-Na(2)**2.d0)
C(1,2)= A*(Nd(2)**2.d0-Na(2)**2.d0)
C(2,2)=T/2.d0/Pi + A*(Na(1)*Na(2)-Nd(1)*Nd(2))
C
C FIN DE LA RUTINA
C
RETURN
C
C FORMATOS DE ESCRITURA
C
08120 FORMAT(/,2x,'ERROR 08120 - CIJ_2 - Error Programa (#',i1,'#)')
C
END
C
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
