C Rutina 053
C------------------------------------------------------------ORTOMAT---
C
SUBROUTINE ortoMat(TP,MAT) !MT000225-1
C
C Esta rutina calcula las constantes que aparecen en la ley de com_
C portamiento y en la solucion fundamental de materiales ortotropos.
C
IMPLICIT NONE
INCLUDE '_PARAM.FOR'
INCLUDE '_STRUCT.FOR'
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES EXTERNAS (ENTRADA)
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C TP : tipo de problema:
INTEGER TP
C MAT : Material
RECORD /S_material/ MAT
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES EXTERNAS (SALIDA)
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C CTESORTO : Ctes del material ortotropo.
C RECORD /S_orto/ MO !MT000225-1
C !MT000225-1
C RECORD /S_material/ MAT !MT000225-1
C !MT000225-1
CCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES INTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C NTIP : Tipo de problema
C TP=31 : Tension plana
C TP=32 : Deformacion plana
INTEGER NTIP
C E(3) : modulos de elasticidad E11, E22, E33
C POIS(3) : coeficiente de Poison POIS12, POIS13, POIS23
C G12 : modulo de cortadura G12
C TETA : Angulo entre los ejes ppales del mat y los globales
REAL*8 E(3),POIS(3),G12,TETA
C PI : 3.14...
REAL*8 PI
C
CCCCCCCCCCCCCCCCCCCCCCCCCC
C VALORES INICIALES
CCCCCCCCCCCCCCCCCCCCCCCCCC
C
PI=DACOS(-1.D0)
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C CONSTANTES ELASTICAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
SELECT CASE (MAT.TIP)
CASE(1)
C Material isotropo
E(1)=MAT.PRP(1)
E(2)=MAT.PRP(1)
E(3)=MAT.PRP(1)
POIS(1)=MAT.PRP(2)
POIS(2)=MAT.PRP(2)
POIS(3)=MAT.PRP(2)
G12=MAT.PRP(1)/2.d0/(1.d0+MAT.PRP(2))
TETA=0.D0
CASE(2)
C Material ortotropo
E(1)=MAT.PRP(1)
E(2)=MAT.PRP(2)
E(3)=MAT.PRP(3)
POIS(1)=MAT.PRP(4)
POIS(2)=MAT.PRP(5)
POIS(3)=MAT.PRP(6)
G12=MAT.PRP(7)
TETA=MAT.PRP(8)*2.D0*PI/360.D0
CASE DEFAULT
WRITE(F_W,09111) MAT.TIP
WRITE( *,09111) MAT.TIP
STOP
END SELECT
C
SELECT CASE (TP)
CASE(31)
NTIP=1
CASE(32)
NTIP=0
CASE DEFAULT
WRITE(F_W,09112) TP
WRITE( *,09112) TP
STOP
END SELECT
C
C Calculo del factor de escala de las tensiones !MT00025-1
C !MT00025-1
SELECT CASE (MAT.TIP) !MT00025-1
CASE(1) !MT00025-1
MAT.E=MAT.PRP(1) !MT00025-1
MAT.D=MAT.PRP(3) !MT00025-1
CASE(2) !MT00025-1
MAT.E=MAT.PRP(1) !MT00025-1
MAT.D=MAT.PRP(9) !MT00025-1
CASE DEFAULT !MT00025-1
WRITE(F_W,04113) MAT.TIP !MT00025-1
WRITE (*,04113) MAT.TIP !MT00025-1
STOP !MT00025-1
END SELECT !MT00025-1
C !MT00025-1
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C CONSTANTES DEL MATERIAL
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
CALL CALCONS(NTIP,E,POIS,G12,MAT.MO.B11,MAT.MO.B22,MAT.MO.B12, !MT00025-1
@ MAT.MO.B66,MAT.MO.GAMMA,MAT.MO.ACU,MAT.MO.BMAS, !MT00025-1
@ MAT.MO.K11,MAT.MO.K22,MAT.MO.K12,MAT.MO.C11, !MT00025-1
@ MAT.MO.C22,MAT.MO.A1,MAT.MO.A2,MAT.MO.B1, !MT00025-1
@ MAT.MO.B2,MAT.MO.C1,MAT.MO.C2,MAT.MO.D1, !MT00025-1
@ MAT.MO.D2,MAT.MO.BMEN,MAT.MO.AMAS,MAT.MO.AMEN, !MT00025-1
@ MAT.MO.MU1,MAT.MO.MU2,MAT.MO.P,MAT.MO.Q, !MT00025-1
@ MAT.MO.CUK,MAT.MO.QP,MAT.MO.PP,MAT.MO.QQP, !MT00025-1
@ MAT.MO.QQQ,MAT.MO.U1,MAT.MO.U2,MAT.MO.T1, !MT00025-1
@ MAT.MO.T2,MAT.MO.DD1,MAT.MO.DD2,MAT.MO.SS1, !MT00025-1
@ MAT.MO.SS2) !MT00025-1
C
RETURN
C
C FORMATOS DE ESCRITURA
C
09111 FORMAT(/,2x,'ERROR 09111 - ORTOMAT - Tipo de material (',i2,
@ ') no implementado.')
09112 FORMAT(/,2x,'ERROR 09112 - ORTOMAT - Tipo de problema (',i2,
@ ') no implementado.')
04113 FORMAT(/,2x,'ERROR 04113 - ORTOMAT',/, !MT00025-1
@ ' Tipo de material ',i2,' no implementado.') !MT00025-1
C
END
C
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
