Télécharger elast4.eso

Retour à la liste

Numérotation des lignes :

  1. C ELAST4 SOURCE AM 08/12/19 21:16:35 6237
  2. SUBROUTINE ELAST4(ISENS,IFOUR,VAR,NNVARI,XMATT,NCOMAT,
  3. & YUNGV,XNUV,XCAR,ICARA,MFR,NSTRS,DD,DDV,KERRE,INDIC,ITHHER)
  4. C
  5. C----------------------------------------------------------------------------
  6. C Objet: Calcul de la matrice de Hooke endommagee(ISENS=1)
  7. C ou son inverse(ISENS=2) dans le cas de la viscoplasticite avec
  8. C endommagement anisotrope (materiau 142) en regime isotherme.
  9. C La matrice est calculée à partir de l'expression de la loi de Hooke
  10. C----------------------------------------------------------------------------
  11. C
  12. C----------------------------------------------------------------------------
  13. C Entree: ISENS= 1 DD en sortie est la matrice de Hooke endommagee
  14. C d'ou SIG=DD*EPSELAS
  15. C = 2 DD en sortie est la matrice de Hooke inverse
  16. C endommagee d'ou EPSELAS=DD*SIG
  17. C avec SIG contraintes et EPSELAS deformations elastiques
  18. C IFOUR= -2 EN CONTR. PLANES
  19. C -1 EN DEFORM. PLANES
  20. C 0 EN AXISYMETRIE
  21. C 1 EN SERIE DE FOURIER
  22. C 2 EN TRIDIM
  23. C NSTRS nombre de composantes des contraintes ou
  24. C des deformations
  25. C NNVARI nombre de variables internes pilotant les
  26. C equations du modele
  27. C VAR(NVARI) tableau des variables internes
  28. C XMATT(NCOMAT) tableau des parametres scalaires du materiau
  29. C a une temperature T donnee
  30. C dont XMATT(1) module d'Young et XMATT(2) coeff de Poisson
  31. C YUNGV derivee /T du module d'Young a T
  32. C XNUV derivee /T du coeff de Poisson a T
  33. C MFR indice de la formulation mecanique(seulement massif
  34. C pour les materiaux endommageables)
  35. C ICARA nombre de caracteristiques geometriques des elements
  36. C finis
  37. C XCAR(ICARA) tableau des caracteristiques geometriques des
  38. C elements finis
  39. C INDIC=0, 1 OU -1 pour plasticite avec endommagement
  40. C =2 OU -2 pour viscoplasticite avec endommagement
  41. C ITHHER = 0 pas de chargement thermique et materiau constant
  42. C = 1 chargement thermique et materiau constant
  43. C = 2 chargement thermique et materiau(T)
  44. C--------------------------------------------------------------------------------
  45. C
  46. C--------------------------------------------------------------------------------
  47. C Sortie: DD(NSTRS,NSTRS) matrice de Hooke
  48. C endommagee si ISENS=1
  49. C son inverse si ISENS=2
  50. C DDV(NSTRS,NSTRS) derivee de la matrice de Hooke
  51. C endommagee si ISENS=1
  52. C son inverse si ISENS=2
  53. C KERRE indice qui regit les erreurs
  54. C Il vaut 99 si la formulation mecanique n'est pas disponible ou
  55. C s'il y a incompatibilite entre MFR et IFOUR
  56. C _______________________________________________________________________________
  57. C
  58. IMPLICIT INTEGER(I-N)
  59. IMPLICIT REAL*8 (A-H,O-Z)
  60. DIMENSION VAR(*),XCAR(*),XMATT(*)
  61. DIMENSION DD(NSTRS,NSTRS),DDV(NSTRS,NSTRS)
  62. C
  63. CALL ZERO(DD,NSTRS,NSTRS)
  64. CALL ZERO(DDV,NSTRS,NSTRS)
  65. KERRE=0
  66. C
  67. C--------------------------------------------------------------------------------
  68. C TERMES DU TENSEUR D'ENDOMMAGEMENT
  69. C--------------------------------------------------------------------------------
  70. C
  71. YUNG=XMATT(1)
  72. XNU=XMATT(2)
  73. D11=VAR(2)
  74. D22=VAR(3)
  75. D33=VAR(4)
  76. D12=VAR(5)
  77. D13=VAR(6)
  78. D23=VAR(7)
  79.  
  80. IF (D11.LE.1.D0) THEN
  81. DEN1 = 1.D0-D11
  82. ELSE
  83. DEN1 = 1.D-6
  84. ENDIF
  85. IF (D22.LE.1.D0) THEN
  86. DEN2 = 1.D0-D22
  87. ELSE
  88. DEN2 = 1.D-6
  89. ENDIF
  90. IF (D33.LE.1.D0) THEN
  91. DEN3 = 1.D0-D33
  92. ELSE
  93. DEN3 = 1.D-6
  94. ENDIF
  95.  
  96. C --- Déterminant de la matrice I - D ---
  97. det=(DEN1)*(DEN2)*(DEN3)
  98. & - 2.D0*D12*D13*D23-D12*D12*(DEN3)
  99. & - D23*D23*(DEN1)-D13*D13*(DEN2)
  100. C
  101. det2=det*det
  102. det4=det2*det2
  103. C
  104. C --- Expression des cofacteurs de la matrice I - D ---
  105. A11=(DEN2)*(DEN3)-D23*D23
  106. A22=(DEN1)*(DEN3)-D13*D13
  107. A33=(DEN2)*(DEN1)-D12*D12
  108. A12=D12*(DEN3)+D23*D13
  109. A13=D13*(DEN2)+D23*D12
  110. A23=D23*(DEN1)+D13*D12
  111. C
  112. COEFF1=(1.D0+XNU)/YUNG
  113. COEFF2=XNU/YUNG
  114. COEFF3=1.D0/YUNG
  115. C
  116. IF (MFR.NE.1) GOTO 40
  117. C
  118. C------------
  119. C MASSIFS
  120. C------------
  121. IF (IFOUR.EQ.-2) GOTO 40
  122. C
  123. C-------------------
  124. C EPSELAS=DD*SIG
  125. C-------------------
  126.  
  127. DD(1,1)=COEFF3*(A11*A11/det2)
  128. DD(1,2)=(COEFF1*A12*A12-COEFF2*A11*A22)/det2
  129. DD(1,3)=(COEFF1*A13*A13-COEFF2*A11*A33)/det2
  130. DD(1,4)=2.D0*COEFF3*A11*A12/det2
  131. IF (IFOUR.GT.0) THEN
  132. DD(1,5)=2.D0*COEFF3*A11*A13/det2
  133. DD(1,6)=2.D0*(COEFF1*A12*A13-COEFF2*A11*A23)/det2
  134. ENDIF
  135.  
  136. DD(2,1)=DD(1,2)
  137. DD(2,2)=COEFF3*A22*A22/det2
  138. DD(2,3)=(COEFF1*A23*A23-COEFF2*A22*A33)/det2
  139. DD(2,4)=2.D0*COEFF3*A12*A22/det2
  140. IF (IFOUR.GT.0) THEN
  141. DD(2,5)=2.D0*(COEFF1*A12*A23-COEFF2*A13*A22)/det2
  142. DD(2,6)=2.D0*COEFF3*A22*A23/det2
  143. ENDIF
  144.  
  145. DD(3,1)=DD(1,3)
  146. DD(3,2)=DD(2,3)
  147. DD(3,3)=COEFF3*A33*A33/det2
  148. DD(3,4)=2.D0*(COEFF1*A13*A23-COEFF2*A33*A12)/det2
  149. IF (IFOUR.GT.0) THEN
  150. DD(3,5)=2.D0*COEFF3*A13*A33/det2
  151. DD(3,6)=2.D0*COEFF3*A33*A23/det2
  152. ENDIF
  153.  
  154. DD(4,1)=DD(1,4)
  155. DD(4,2)=DD(2,4)
  156. DD(4,3)=DD(3,4)
  157. DD(4,4)=2.D0*(COEFF1*(A12*A12+A11*A22)-COEFF2*2.D0*A12*A12)/det2
  158. IF (IFOUR.GT.0) THEN
  159. DD(4,5)=2.D0*(COEFF1*(A12*A13+A11*A23)-COEFF2*2.D0*A13*A12)/det2
  160. DD(4,6)=2.D0*(COEFF1*(A12*A23+A22*A13)-COEFF2*2.D0*A23*A12)/det2
  161. ENDIF
  162.  
  163. IF (IFOUR.GT.0) THEN
  164. DD(5,1)=DD(1,5)
  165. DD(5,2)=DD(2,5)
  166. DD(5,3)=DD(3,5)
  167. DD(5,4)=DD(4,5)
  168. DD(5,5)=2.D0*(COEFF1*(A13*A13+A11*A33)-COEFF2*2.D0*A13*A13)/det2
  169. DD(5,6)=2.D0*(COEFF1*(A13*A23+A33*A12)-COEFF2*2.D0*A23*A13)/det2
  170. ENDIF
  171.  
  172.  
  173. IF (IFOUR.GT.0) THEN
  174. DD(6,1)=DD(1,6)
  175. DD(6,2)=DD(2,6)
  176. DD(6,3)=DD(3,6)
  177. DD(6,4)=DD(4,6)
  178. DD(6,5)=DD(5,6)
  179. DD(6,6)=2.D0*(COEFF1*(A23*A23+A33*A22)-COEFF2*2.D0*A23*A23)/det2
  180. ENDIF
  181. C
  182. IF (ISENS.EQ.2) GOTO 500
  183. C
  184. C----------------
  185. C SIG=DD*EPSELAS
  186. C----------------
  187. CALL INVALM (DD, NSTRS , NSTRS , KER, 1.D-12)
  188. GOTO 500
  189. C
  190. C FORMULATION NON DISPONIBLE
  191. 40 KERRE=99
  192. 500 CONTINUE
  193. RETURN
  194. END
  195.  
  196.  
  197.  
  198.  
  199.  
  200.  
  201.  
  202.  
  203.  
  204.  
  205.  
  206.  
  207.  
  208.  
  209.  

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