C Rutina 052
C---------------------------------------------------------------CIJ---
C
C SUBRUTINA CIJ
C
C CALCULA EL TERMINO LIBRE
C
SUBROUTINE CIJ(PC,MAT,NrmA,NrmD,DT,TP,COD)
C
C
IMPLICIT NONE
INCLUDE '_param.for'
INCLUDE '_struct.for'
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES EXTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C PC : Puntos de colocacion
RECORD /S_coloc/ PC
C MAT : Constantes del material
RECORD /S_material/ MAT
C NrmA,NrmD : Normal correspondiente a los elementos anterior y posterior
REAL*8 NrmA(2),NrmD(2)
C DT : Angulo interior en el nodo
REAL*8 DT
C TM : Tipo de material
INTEGER TM
C TP : Tipo de problema
INTEGER TP
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
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VALORES INICIALES
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C CALCULO TERMINO LIBRE
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C
SELECT CASE(TP)
CASE(1)
C Elasticidad axisimétrica
CALL CIJ_1(PC.CIJ,MAT.PRP(2),NrmA,NrmD,DT,COD)
CASE(21,22)
C Elasticidad plana (clásica)
CALL CIJ_2(PC.CIJ,MAT.PRP(2),NrmA,NrmD,DT,TP)
CASE(31,32)
C Elasticidad plana (comleja) ortótropos
CALL CIJ_3(PC,MAT,NrmA,NrmD)
CASE(51,52)
C Elasticidad plana (comleja) anisótropos !MT000705
CALL CIJ_4(PC,MAT,NrmA,NrmD,TM)
CASE DEFAULT
WRITE(F_W,08101) TP
WRITE (*,08101) TP
STOP
END SELECT
C
C FIN DE LA RUTINA
C
RETURN
C
C FORMATOS DE ESCRITURA
C
08100 FORMAT(/,2x,'ERROR 08100 - CIJ - Error Programa (#',i1,'#)')
08101 FORMAT(/,2x,'ERROR 08101 - CIJ - Tipo de problema (',i2,
@ ') no implementado.')
C
END
C
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
