Télécharger elast.eso

Retour à la liste

Numérotation des lignes :

  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
  27. CALL CELAST (TT, PCOEL, E1,E2,E3, E1P,E2P,E3P,
  28. & XNU12,XNU23,XNU31, XNU12P,XNU23P,XNU31P,
  29. & XMU12,XMU13,XMU23, XMU12P,XMU13P,XMU23P)
  30. SOEL(1,1) = 1./E1
  31. SOEL(2,2) = 1./E2
  32. SOEL(3,3) = 1./E3
  33. SKL = -XNU12/E1
  34. SOEL(1,2) = SKL
  35. SOEL(2,1) = SKL
  36. SKL = -XNU23/E2
  37. SOEL(2,3) = SKL
  38. SOEL(3,2) = SKL
  39. SKL = -XNU31/E3
  40. SOEL(3,1) = SKL
  41. SOEL(1,3) = SKL
  42. SOEL(4,4) = 0.5/XMU12
  43. SOEL(5,5) = 0.5/XMU13
  44. SOEL(6,6) = 0.5/XMU23
  45. S11P = -E1P/E1**2
  46. S22P = -E2P/E2**2
  47. S33P = -E3P/E3**2
  48. S12P = -XNU12P/E1-XNU12*S11P
  49. S23P = -XNU23P/E2-XNU23*S22P
  50. S31P = -XNU31P/E3-XNU31*S33P
  51. S44P = -0.5*XMU12P/XMU12**2
  52. S55P = -0.5*XMU13P/XMU13**2
  53. S66P = -0.5*XMU23P/XMU23**2
  54. VEPEL0(1) = (S11P*SIG(1)+S12P*SIG(2)+S31P*SIG(3))*TTP
  55. VEPEL0(2) = (S12P*SIG(1)+S22P*SIG(2)+S23P*SIG(3))*TTP
  56. VEPEL0(3) = (S31P*SIG(1)+S23P*SIG(2)+S33P*SIG(3))*TTP
  57. VEPEL0(4) = S44P*SIG(4)*TTP
  58. VEPEL0(5) = S55P*SIG(5)*TTP
  59. VEPEL0(6) = S66P*SIG(6)*TTP
  60. RETURN
  61. END
  62.  
  63.  

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