Télécharger celast.eso

Retour à la liste

Numérotation des lignes :

  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