Télécharger la plaquette Cast3M


fr en

Télécharger celast.eso

Retour à la liste

Numérotation des lignes : 

celast

  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.  
  71.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales