Esope source files

   1 : C CELAST    SOURCE    STRU      08/12/22    21:15:01     6246           
   2 : c     Sous-Programme du module MISTRAL0
   3 : C     --------------------------------------------------------------------------   
   4 :       SUBROUTINE CELAST (TT, PCOEL, E1,E2,E3, E1P,E2P,E3P,
   5 :      &                   XNU12,XNU23,XNU31, XNU12P,XNU23P,XNU31P,
   6 :      &                   XMU12,XMU13,XMU23, XMU12P,XMU13P,XMU23P)           
   7 : C     --------------------------------------------------------------------------
   8 : C     Calcul des :
   9 : C     - modules d'Young :          E1, E2, E3,
  10 : C     - coefficients de Poisson :  XNU12, XNU23, XNU31,
  11 : C     - modules de cisaillement :  XMU12, XMU13, XMU23
  12 : C     et de leurs derivees en fonction de la temperature TT :
  13 : C     E1P, E2P, E3P, XNU12P, XNU23P, XNU31P, XMU12P, XMU13P, XMU23P.
  14 : C     --------------------------------------------------------------------------
  15 :       IMPLICIT INTEGER (I-N)
  16 :       IMPLICIT REAL*8 (A-H, O-Z)
  17 :         DIMENSION PCOEL(1:*)
  18 :       ISOEL = NINT(PCOEL(1))
  19 :       IF (ISOEL.EQ.1) THEN
  20 : C     -----------------
  21 : C     Materiau isotrope
  22 : C       
  23 :          YOUNG = PCOEL(2)+PCOEL(3)*TT
  24 :          POIS = PCOEL(4)+PCOEL(5)*TT
  25 :          CISAIL = YOUNG/(2.*(1.+POIS))
  26 :          CISAIP = (PCOEL(3)/(1.+POIS)-YOUNG*PCOEL(5)/(1.+POIS)**2)/2.
  27 :          E1 = YOUNG   
  28 :          E1P = PCOEL(3)
  29 :          E2 = YOUNG
  30 :          E2P = PCOEL(3)
  31 :          E3 = YOUNG
  32 :          E3P = PCOEL(3)
  33 :          XNU12 = POIS
  34 :          XNU12P = PCOEL(5)
  35 :          XNU23 = POIS     
  36 :          XNU23P = PCOEL(5)
  37 :          XNU31 = POIS     
  38 :          XNU31P = PCOEL(5)
  39 :          XMU12 = CISAIL   
  40 :          XMU12P = CISAIP
  41 :          XMU13 = CISAIL   
  42 :          XMU13P = CISAIP
  43 :          XMU23 = CISAIL   
  44 :          XMU23P = CISAIP
  45 :       ELSE
  46 : C     -------------------
  47 : C     Materiau orthotrope
  48 : C       
  49 :          E1 = PCOEL(2)+PCOEL(3)*TT
  50 :          E1P = PCOEL(3)
  51 :          E2 = PCOEL(4)+PCOEL(5)*TT         
  52 :          E2P = PCOEL(5)
  53 :          E3 = PCOEL(6)+PCOEL(7)*TT         
  54 :          E3P = PCOEL(7)
  55 :          XNU12 = PCOEL(8)+PCOEL(9)*TT       
  56 :          XNU12P = PCOEL(9)
  57 :          XNU23 = PCOEL(10)+PCOEL(11)*TT     
  58 :          XNU23P = PCOEL(11)
  59 :          XNU31 = PCOEL(12)+PCOEL(13)*TT     
  60 :          XNU31P = PCOEL(13)
  61 :          XMU12 = PCOEL(14)+PCOEL(15)*TT     
  62 :          XMU12P = PCOEL(15)
  63 :          XMU13 = PCOEL(16)+PCOEL(17)*TT     
  64 :          XMU13P = PCOEL(17)
  65 :          XMU23 = PCOEL(18)+PCOEL(19)*TT     
  66 :          XMU23P = PCOEL(19)
  67 :         END IF
  68 :       RETURN
  69 :         END
  70 :