Esope source files

   1 : C ELAST     SOURCE    STRU      08/12/22    21:15:15     6246           
   2 : c     Sous-Programme du module MISTRAL0
   3 : C     --------------------------------------------------------------------------   
   4 :       SUBROUTINE ELAST (TT,TTP,SIG, PCOEL, SOEL,VEPEL0)
   5 : C     --------------------------------------------------------------------------
   6 : C     Calcul :
   7 : C     - des coefficients d'elasticite skl (de souplesse) :
   8 : C             s11 = 1/E1          s22 = 1/E2          s33 = 1/E3 
   9 : C             s12 = -XNU12/E1     s23 = -XNU23/E2     s31 = -XNU31/E3 
  10 : C             s44 = 1/(2.XMU12)   s55 = 1/(2.XMU13)   s66 = 1/(2.XMU23)
  11 : C     stockes dans SOEL(K,L),
  12 : C     - de la matrice VEPEL0 definie par :
  13 : C                           VEPSEL = VEPEL0 + SOEL*SIGP
  14 : C     donnant les vitesses de deformations elastiques VEPSEL(IJ) en fonction
  15 : C     des derivees par rapport au temps des contraintes SIGP(IPJP). 
  16 : C     --------------------------------------------------------------------------
  17 :       IMPLICIT INTEGER (I-N)
  18 :       IMPLICIT REAL*8 (A-H, O-Z)
  19 :         DIMENSION SIG(1:*),PCOEL(1:*)
  20 :         DIMENSION SOEL(1:6,1:*),VEPEL0(1:*)
  21 :       DO K = 4,6
  22 :          DO L = 1,K-1
  23 :           SOEL(K,L) = 0.
  24 :           SOEL(L,K) = 0.
  25 :          END DO
  26 :         END DO
  26 :       CALL CELAST (TT, PCOEL, E1,E2,E3, E1P,E2P,E3P,
  27 :      &             XNU12,XNU23,XNU31, XNU12P,XNU23P,XNU31P,
  28 :      &             XMU12,XMU13,XMU23, XMU12P,XMU13P,XMU23P)
  29 :       SOEL(1,1) = 1./E1
  30 :         SOEL(2,2) = 1./E2
  31 :         SOEL(3,3) = 1./E3
  32 :         SKL = -XNU12/E1
  33 :         SOEL(1,2) = SKL
  34 :         SOEL(2,1) = SKL
  35 :         SKL = -XNU23/E2
  36 :         SOEL(2,3) = SKL
  37 :         SOEL(3,2) = SKL
  38 :         SKL = -XNU31/E3
  39 :         SOEL(3,1) = SKL
  40 :         SOEL(1,3) = SKL
  41 :       SOEL(4,4) = 0.5/XMU12
  42 :       SOEL(5,5) = 0.5/XMU13
  43 :       SOEL(6,6) = 0.5/XMU23
  44 :       S11P = -E1P/E1**2
  45 :       S22P = -E2P/E2**2
  46 :       S33P = -E3P/E3**2
  47 :       S12P = -XNU12P/E1-XNU12*S11P
  48 :       S23P = -XNU23P/E2-XNU23*S22P
  49 :       S31P = -XNU31P/E3-XNU31*S33P
  50 :       S44P = -0.5*XMU12P/XMU12**2
  51 :       S55P = -0.5*XMU13P/XMU13**2
  52 :       S66P = -0.5*XMU23P/XMU23**2
  53 :         VEPEL0(1) = (S11P*SIG(1)+S12P*SIG(2)+S31P*SIG(3))*TTP
  54 :         VEPEL0(2) = (S12P*SIG(1)+S22P*SIG(2)+S23P*SIG(3))*TTP
  55 :         VEPEL0(3) = (S31P*SIG(1)+S23P*SIG(2)+S33P*SIG(3))*TTP
  56 :         VEPEL0(4) = S44P*SIG(4)*TTP      
  57 :         VEPEL0(5) = S55P*SIG(5)*TTP      
  58 :         VEPEL0(6) = S66P*SIG(6)*TTP      
  59 :         RETURN
  60 :         END
  61 :