Télécharger elast3.eso

Retour à la liste

Numérotation des lignes :

  1. C ELAST3 SOURCE CHAT 05/01/12 23:29:35 5004
  2. SUBROUTINE ELAST3(ISENS,IFOUR,XMATT,MFR,NSTRS,DD,KERRE,INDIC)
  3. C
  4. C----------------------------------------------------------------------------
  5. C Objet: Calcul de la matrice de Hooke endommagee(ISENS=1)
  6. C ou son inverse(ISENS=2)
  7. C dans le cadre du modele CERAMIQUE (materiau 65)
  8. C----------------------------------------------------------------------------
  9. C
  10. C----------------------------------------------------------------------------
  11. C Entree: ISENS= 1 DD en sortie est la matrice de Hooke endommagee
  12. C d'ou SIG=DD*EPSELAS
  13. C = 2 DD en sortie est la matrice de Hooke inverse
  14. C endommagee d'ou EPSELAS=DD*SIG
  15. C avec SIG contraintes et EPSELAS deformations elastiques
  16. C IFOUR= -2 EN CONTR. PLANES
  17. C -1 EN DEFORM. PLANES
  18. C 0 EN AXISYMETRIE
  19. C 1 EN SERIE DE FOURIER
  20. C 2 EN TRIDIM
  21. C NSTRS nombre de composantes des contraintes ou
  22. C des deformations
  23. C XMATT(NCOMAT) tableau des parametres scalaires du materiau
  24. C a une temperature T donnee
  25. C dont XMATT(1) module d'Young et XMATT(2) coeff de Poisson
  26. C MFR indice de la formulation mecanique (seulement massif
  27. C pour les materiaux endommageables)
  28. C INDIC = 0 ====> le point de Gauss n'est pas cassé
  29. C INDIC = 1 ====> le point de Gauss est cassé
  30. C--------------------------------------------------------------------------------
  31. C
  32. C--------------------------------------------------------------------------------
  33. C Sortie: DD(NSTRS,NSTRS) matrice de Hooke endommagee si ISENS=1
  34. C son inverse si ISENS=2
  35. C
  36. C KERRE indice qui regit les erreurs
  37. C Il vaut 99 si la formulation mecanique n'est pas disponible ou
  38. C s'il y a incompatibilite entre MFR et IFOUR
  39. C _______________________________________________________________________________
  40. C
  41. IMPLICIT INTEGER(I-N)
  42. IMPLICIT REAL*8 (A-H,O-Z)
  43. DIMENSION XMATT(*)
  44. DIMENSION DD(NSTRS,NSTRS)
  45. C
  46. CALL ZERO(DD,NSTRS,NSTRS)
  47. KERRE=0
  48. C
  49. YUNG=XMATT(1)
  50. XNU=XMATT(2)
  51. C Modifs Cmau
  52. C D1=1.D0-XD
  53. D1=1.D0
  54. C
  55. C ====================================================================
  56. C SI LE PT. DE GAUSS EST ROMPU, ON DIV. SON MOD. D'YOUNG PAR 1 MILLION
  57. C ====================================================================
  58. IF (INDIC.EQ.1) THEN
  59. D1=1.D-6
  60. ENDIF
  61. C
  62. COEFF1=(D1*YUNG)/(1.D0-XNU*XNU)
  63. COEFF2=1.D0/(D1*YUNG)
  64. COEFF3=(D1*YUNG)/((1.D0+XNU)*(1.D0-2.D0*XNU))
  65. COEFF4=(1.D0-2.D0*XNU)/2.D0
  66. C
  67. IF (MFR.EQ.1) THEN
  68. C
  69. C------------
  70. C MASSIFS
  71. C------------
  72. IF (IFOUR.NE.-2) THEN
  73. C
  74. C ----------------------------------
  75. C DEFORMATIONS PLANES ET AXISYMETRIE
  76. C ----------------------------------
  77. IF (ISENS.EQ.1) THEN
  78. C
  79. C----------------
  80. C SIG=DD*EPSELAS
  81. C----------------
  82. DD(1,1)=COEFF3*(1.D0-XNU)
  83. DD(1,2)=COEFF3*XNU
  84. DD(1,3)=DD(1,2)
  85. DD(2,2)=DD(1,1)
  86. DD(2,3)=DD(1,2)
  87. DD(3,3)=DD(1,1)
  88. DD(2,1)=DD(1,2)
  89. DD(3,1)=DD(1,2)
  90. DD(3,2)=DD(1,2)
  91. C
  92. LA=4
  93. LB=4
  94. C
  95. C ---------------------------
  96. C LB=6 POUR TRIDIM ET FOURIER
  97. C ---------------------------
  98. IF (IFOUR.GT.0) LB=6
  99. DO 5 IA=LA,LB
  100. DD(IA,IA)=COEFF3*COEFF4
  101. 5 CONTINUE
  102. ELSE
  103. C
  104. C-------------------
  105. C EPSELAS=DD*SIG
  106. C-------------------
  107. DD(1,1)=COEFF2
  108. DD(1,2)=-XNU*COEFF2
  109. DD(1,3)=DD(1,2)
  110. DD(2,2)=DD(1,1)
  111. DD(2,3)=DD(1,2)
  112. DD(3,3)=DD(1,1)
  113. DD(2,1)=DD(1,2)
  114. DD(3,1)=DD(1,3)
  115. DD(3,2)=DD(2,3)
  116. LA=4
  117. LB=4
  118. C
  119. C ---------------------------
  120. C LB=6 POUR TRIDIM ET FOURIER
  121. C ---------------------------
  122. IF (IFOUR.GT.0) LB=6
  123. DO 10 IA=LA,LB
  124. DD(IA,IA)=(1.D0+XNU)*DD(1,1)*2.D0
  125. 10 CONTINUE
  126. C
  127. ENDIF
  128. C
  129. C -------------
  130. C CONTR. PLANES
  131. C -------------
  132. ELSE
  133. IF (ISENS.EQ.1) THEN
  134. C
  135. C-----------------
  136. C SIG=DD*EPSELAS
  137. C-----------------
  138. DD(1,1)=COEFF1
  139. DD(1,2)=COEFF1*XNU
  140. DD(2,2)=COEFF1
  141. DD(4,4)=COEFF1*((1.D0-XNU)/2.D0)
  142. DD(2,1)=DD(1,2)
  143. ELSE
  144. C
  145. C------------------
  146. C EPSELAS=DD*SIG
  147. C------------------
  148. DD(1,1)=COEFF2
  149. DD(1,2)=-COEFF2*XNU
  150. DD(2,2)=COEFF2
  151. DD(4,4)=COEFF2*2.D0*(1.D0+XNU)
  152. DD(2,1)=DD(1,2)
  153. DD(3,1)=DD(1,2)
  154. DD(3,2)=DD(1,2)
  155. C
  156. ENDIF
  157. ENDIF
  158. C
  159. ELSE
  160. IF (MFR.EQ.5) THEN
  161. C
  162. C------------------------------------------
  163. C COQUES EPAISSES (NON ENCORE DISPONIBLE)
  164. C------------------------------------------
  165. C
  166. C FORMULATION NON DISPONIBLE
  167. ENDIF
  168. KERRE=99
  169. ENDIF
  170. C
  171. RETURN
  172. END
  173.  
  174.  
  175.  
  176.  
  177.  

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