Télécharger hookp.eso

Retour à la liste

Numérotation des lignes :

  1. C HOOKP SOURCE PV 11/02/18 21:15:00 6875
  2. C
  3. SUBROUTINE HOOKP(VAR,VARD,NVARI,XMAT,NCOMAT,MFR,NSTRS,DINV)
  4. C
  5. C----------------------------------------------------------------------------
  6. C Objet: Calcul de la variation de la matrice de Hooke endommagee inverse
  7. C dans le cas de la viscoplasticite avec endommagement
  8. C anisotrope(materiau 142) en regime isotherme
  9. C----------------------------------------------------------------------------
  10. C
  11. C----------------------------------------------------------------------------
  12. C Entree: NSTRS nombre de composantes des contraintes ou
  13. C des deformations
  14. C NVARI nombre de variables internes pilotant les
  15. C equations du modele
  16. C VAR(NVARI) tableau des variables internes
  17. C VARD(NVARI) tableau des variations des variables internes
  18. C XMAT(NCOMAT) tableau des parametres scalaires du materiau
  19. C a une temperature T donnee
  20. C dont XMAT(1) module d'Young et XMAT(2) coeff de Poisson
  21. C MFR indice de la formulation mecanique(seulement massif
  22. C pour les materiaux endommageables)
  23. C--------------------------------------------------------------------------------
  24. C
  25. C--------------------------------------------------------------------------------
  26. C Sortie: DINV(NSTRS,NSTRS) variation de l'inverse de la matrice de Hooke endommagee
  27. C _______________________________________________________________________________
  28. C
  29. IMPLICIT INTEGER(I-N)
  30. IMPLICIT REAL*8(A-H,O-Z)
  31.  
  32. -INC PPARAM
  33. -INC CCOPTIO
  34. DIMENSION VAR(*),VARD(*),XMAT(*),DINV(NSTRS,NSTRS)
  35. C
  36. C --- Termes d'endommagement ---
  37. C
  38. D11=VAR(2)
  39. D22=VAR(3)
  40. D33=VAR(4)
  41. D12=VAR(5)
  42. D13=VAR(6)
  43. D23=VAR(7)
  44.  
  45. DE1 = 1.D0-D11
  46. DE2 = 1.D0-D22
  47. DE3 = 1.D0-D33
  48. IF (D11.GE.1.D0) DE1 = 1.D-6
  49. IF (D22.GE.1.D0) DE2 = 1.D-6
  50. IF (D33.GE.1.D0) DE3 = 1.D-6
  51. C write(6,*) DE1,DE2,DE3
  52. C
  53. C --- Coefficients de la matrive de Hooke inversée ---
  54.  
  55. YUNG=XMAT(1)
  56. XNU=XMAT(2)
  57. COEFF1=(1.D0+XNU)/YUNG
  58. COEFF2=XNU/YUNG
  59. COEFF3=1.D0/YUNG
  60. C
  61. C --- Calcul du déterminant de I - D ---
  62.  
  63. det=(DE1)*(DE2)*(DE3)
  64. & - 2.D0*D12*D13*D23-D12*D12*(DE3)
  65. & - D23*D23*(DE1)-D13*D13*(DE2)
  66.  
  67. det2=det*det
  68. det4=det2*det2
  69. C
  70. C --- Calcul des dérivées partielles de det par rapport aux termes du tenseur d'endommagement
  71.  
  72. detp11=D23*D23-(DE2)*(DE3)
  73. detp22=D13*D13-(DE1)*(DE3)
  74. detp33=D12*D12-(DE1)*(DE2)
  75. detp12=-2.D0*(D23*D13+(DE3)*D12)
  76. detp13=-2.D0*(D12*D23+(DE2)*D13)
  77. detp23=-2.D0*(D12*D13+(DE1)*D23)
  78. C
  79. C --- Calcul des termes de la matrice de l'inverse de I - D (multipliés par det)
  80.  
  81. A11=(DE2)*(DE3)-D23*D23
  82. A22=(DE1)*(DE3)-D13*D13
  83. A33=(DE2)*(DE1)-D12*D12
  84. A12=D12*(DE3)+D23*D13
  85. A13=D13*(DE2)+D23*D12
  86. A23=D23*(DE1)+D13*D12
  87. C
  88. C --- Calcul des termes de la variation de la matrice de Hooke inverse
  89.  
  90. DINV(1,1)=COEFF3*(vard(2)*(-2.D0*detp11*det*A11*A11)
  91. & +vard(3)*(2.D0*(D33-1.D0)*A11*det2-2.D0*detp22*det*A11*A11)
  92. & +vard(4)*(2.D0*(D22-1.D0)*A11*det2-2.D0*detp33*det*A11*A11)
  93. & +vard(5)*(-2.D0*detp12*det*A11*A11)
  94. & +vard(6)*(-2.D0*detp13*det*A11*A11)
  95. & +vard(7)*(-4.D0*D23*A11*det2-2.D0*detp23*det*A11*A11))/det4
  96.  
  97. DINV(1,2)=(COEFF1*(vard(2)*(-2.D0*detp11*det*A12*A12)
  98. & +vard(3)*(-2.D0*detp22*det*A12*A12)
  99. & +vard(4)*(-2.D0*D12*A12*det2-2.D0*detp33*det*A12*A12)
  100. & +vard(5)*(2.D0*(DE3)*A12*det2-2.D0*detp12*det*A12*A12)
  101. & +vard(6)*(2.D0*D23*A12*det2-2.D0*detp13*det*A12*A12)
  102. & +vard(7)*(2.D0*D13*A12*det2-2.D0*detp23*det*A12*A12))
  103. & -COEFF2*(vard(2)*((D33-1.D0)*A11*det2
  104. & -2.D0*detp11*det*A11*A22)
  105. & +vard(3)*((D33-1.D0)*A22*det2-2.D0*detp22*det*A11*A22)
  106. & +vard(4)*(((D22-1.D0)*A22+(D11-1.D0)*A11)*det2
  107. & -2.D0*detp33*det*A11*A22)
  108. & +vard(5)*(-2.D0*detp12*det*A11*A22)
  109. & +vard(6)*(-2.D0*D13*A11*det2-2.D0*detp13*det*A11*A22)
  110. & +vard(7)*(-2.D0*D23*A22*det2-2.D0*detp23*det*A11*A22)))/det4
  111.  
  112. DINV(1,3)=(COEFF1*(vard(2)*(-2.D0*detp11*det*A13*A13)
  113. & +vard(3)*(-2.D0*D13*A13*det2-2.D0*detp22*det*A13*A13)
  114. & +vard(4)*(-2.D0*detp33*det*A13*A13)
  115. & +vard(5)*(2.D0*D23*A13*det2-2.D0*detp12*det*A13*A13)
  116. & +vard(6)*(2.D0*(DE2)*A13*det2-2.D0*detp13*det*A13*A13)
  117. & +vard(7)*(2.D0*D12*A13*det2-2.D0*detp23*det*A13*A13))
  118. & -COEFF2*(vard(2)*((D22-1.D0)*A11*det2
  119. & -2.D0*detp11*det*A11*A33)
  120. & +vard(3)*(((D33-1.D0)*A33+(D11-1.D0)*A11)*det2
  121. & -2.D0*detp22*det*A11*A33)
  122. & +vard(4)*((D22-1.D0)*A33*det2-2.D0*detp33*det*A11*A33)
  123. & +vard(5)*(-2.D0*D12*A11*det2-2.D0*detp12*det*A11*A33)
  124. & +vard(6)*(-2.D0*detp13*det*A11*A33)
  125. & +vard(7)*(-2.D0*D23*A33*det2-2.D0*detp23*det*A11*A33)))/det4
  126.  
  127. DINV(1,4)=(2.D0*COEFF3/det4)*(vard(2)*(-2.D0*detp11*det*A11*A12)
  128. & +vard(3)*((D33-1.D0)*A12*det2-2.D0*detp22*det*A11*A12)
  129. & +vard(4)*(((D22-1.D0)*A12-D12*A11)*det2
  130. & -2.D0*detp33*det*A11*A12)
  131. & +vard(5)*((DE3)*A11*det2-2.D0*detp12*det*A11*A12)
  132. & +vard(6)*(D23*A11*det2-2.D0*detp13*det*A11*A12)
  133. & +vard(7)*((-2.D0*D23*A12+D13*A11)*det2
  134. & -2.D0*detp23*det*A11*A12))
  135.  
  136. IF (IFOUR.GT.0) THEN
  137. DINV(1,5)=(2.D0*COEFF3/det4)*(vard(2)*(-2.D0*detp11*det*A11*A13)
  138. & +vard(3)*(((D33-1.D0)*A13-D13*A11)*det2
  139. & -2.D0*detp22*det*A11*A13)
  140. & +vard(4)*((D22-1.D0)*A13*det2-2.D0*detp33*det*A11*A13)
  141. & +vard(5)*(D23*A11*det2-2.D0*detp12*det*A11*A13)
  142. & +vard(6)*((DE2)*A11*det2-2.D0*detp13*det*A11*A13)
  143. & +vard(7)*((-2.D0*D23*A13+D12*A11)*det2
  144. & -2.D0*detp23*det*A11*A13))
  145.  
  146. DINV(1,6)=(COEFF1*(vard(2)*(-2.D0*detp11*det*A12*A13)
  147. & +vard(3)*(-D13*A12*det2-2.D0*detp22*det*A12*A13)
  148. & +vard(4)*(-D12*A13*det2-2.D0*detp33*det*A12*A13)
  149. & +vard(5)*(((DE3)*A13+D23*A12)*det2
  150. & -2.D0*detp12*det*A12*A13)
  151. & +vard(6)*((D23*A13+(DE2)*A12)*det2
  152. & -2.D0*detp13*det*A12*A13)
  153. & +vard(7)*((D13*A13+D12*A12)*det2-2.D0*detp23*det*A12*A13))
  154. & -COEFF2*(vard(2)*(-D23*A11*det2-2.D0*detp11*det*A11*A23)
  155. & +vard(3)*((D33-1.D0)*A23*det2-2.D0*detp22*det*A11*A23)
  156. & +vard(4)*((D22-1.D0)*A23*det2-2.D0*detp33*det*A11*A23)
  157. & +vard(5)*(D13*A11*det2-2.D0*detp12*det*A11*A23)
  158. & +vard(6)*(D12*A11*det2-2.D0*detp13*det*A11*A23)
  159. & +vard(7)*(((DE1)*A11-2.D0*D23*A23)*det2
  160. & -2.D0*detp23*det*A11*A23)))*2.D0/det4
  161. ENDIF
  162.  
  163. DINV(2,2)=COEFF3*(vard(2)*(2.D0*(D33-1.D0)*A22*det2
  164. & -2.D0*detp11*det*A22*A22)
  165. & +vard(3)*(-2.D0*detp22*det*A22*A22)
  166. & +vard(4)*(2.D0*(D11-1.D0)*A22*det2-2.D0*detp33*det*A22*A22)
  167. & +vard(5)*(-2.D0*detp12*det*A22*A22)
  168. & +vard(6)*(-4.D0*D13*A22*det2-2.D0*detp13*det*A22*A22))
  169. & +vard(7)*(-2.D0*detp23*det*A22*A22)/det4
  170.  
  171. DINV(2,3)=(COEFF1*(vard(2)*(-2.D0*D23*A23*det2
  172. & -2.D0*detp11*det*A23*A23)
  173. & +vard(3)*(-2.D0*detp22*det*A23*A23)
  174. & +vard(4)*(-2.D0*detp33*det*A23*A23)
  175. & +vard(5)*(2.D0*D13*A23*det2-2.D0*detp12*det*A23*A23)
  176. & +vard(6)*(2.D0*D12*A23*det2-2.D0*detp13*det*A23*A23)
  177. & +vard(7)*(2.D0*(DE1)*A23*det2-2.D0*detp23*det*A23*A23))
  178. & -COEFF2*(vard(2)*(((D33-1.D0)*A33+(D22-1.D0)*A22)*det2
  179. & -2.D0*detp11*det*A22*A33)
  180. & +vard(3)*((D11-1.D0)*A22*det2-2.D0*detp22*det*A22*A33)
  181. & +vard(4)*((D11-1.D0)*A33*det2-2.D0*detp33*det*A22*A33)
  182. & +vard(5)*(-2.D0*D12*A22*det2-2.D0*detp12*det*A22*A33)
  183. & +vard(6)*(-2.D0*D13*A33*det2-2.D0*detp13*det*A22*A33)
  184. & +vard(7)*(-2.D0*detp23*det*A22*A33)))/det4
  185.  
  186. DINV(2,4)=(vard(2)*((D33-1.D0)*A12*det2-2.D0*detp11*det*A22*A12)
  187. & +vard(3)*(-2.D0*detp22*det*A22*A12)
  188. & +vard(4)*(((D11-1.D0)*A12-D12*A22)*det2
  189. & -2.D0*detp33*det*A22*A12)
  190. & +vard(5)*((DE3)*A22*det2-2.D0*detp12*det*A22*A12)
  191. & +vard(6)*((-2.D0*D13*A12+D23*A22)*det2
  192. & -2.D0*detp13*det*A22*A12)
  193. & +vard(7)*(D13*A22*det2-2.D0*detp23*det*A22*A12))
  194. & *(2.D0*COEFF3/det4)
  195.  
  196. IF (IFOUR.GT.0) THEN
  197. DINV(2,5)=(COEFF1*(vard(2)*(-D23*A12*det2
  198. & -2.D0*detp11*det*A12*A23)
  199. & +vard(3)*(-2.D0*detp22*det*A12*A23)
  200. & +vard(4)*(-D12*A23*det2-2.D0*detp33*det*A12*A23)
  201. & +vard(5)*(((DE3)*A23+D13*A12)*det2
  202. & -2.D0*detp12*det*A12*A23)
  203. & +vard(6)*((D23*A23+D12*A12)*det2-2.D0*detp13*det*A12*A23)
  204. & +vard(7)*((D13*A23+(DE1)*A12)*det2
  205. & -2.D0*detp23*det*A12*A23))
  206. & -COEFF2*(vard(2)*((D33-1.D0)*A13*det2
  207. & -2.D0*detp11*det*A22*A13)
  208. & +vard(3)*(-D13*A22*det2-2.D0*detp22*det*A22*A13)
  209. & +vard(4)*((D11-1.D0)*A13*det2-2.D0*detp33*det*A22*A13)
  210. & +vard(5)*(D23*A22*det2-2.D0*detp12*det*A22*A13)
  211. & +vard(6)*(((DE2)*A22-2.D0*D13*A13)*det2
  212. & -2.D0*detp13*det*A22*A13)
  213. & +vard(7)*(D12*A22*det2-2.D0*detp23*det*A22*A13)))*2.D0/det4
  214.  
  215. DINV(2,6)=2.D0*COEFF3*(vard(2)*(((D33-1.D0)*A23-D23*A22)*det2
  216. & -2.D0*detp11*det*A22*A23)
  217. & +vard(3)*(-2.D0*detp22*det*A22*A23)
  218. & +vard(4)*((D11-1.D0)*A23*det2-2.D0*detp33*det*A22*A23)
  219. & +vard(5)*(D13*A22*det2-2.D0*detp12*det*A22*A23)
  220. & +vard(6)*((-2.D0*D13*A23+D12*A22)*det2
  221. & -2.D0*detp13*det*A22*A23)
  222. & +vard(7)*((DE1)*A22*det2-2.D0*detp23*det*A22*A23))/det4
  223. ENDIF
  224.  
  225. DINV(3,3)=COEFF3*(vard(2)*(2.D0*(D22-1.D0)*A33*det2
  226. & -2.D0*detp11*det*A33*A33)
  227. & +vard(3)*(2.D0*(D11-1.D0)*A33*det2-2.D0*detp22*det*A33*A33)
  228. & +vard(4)*(-2.D0*detp33*det*A33*A33)
  229. & +vard(5)*(-4.D0*D12*A33*det2-2.D0*detp12*det*A33*A33)
  230. & +vard(6)*(-2.D0*detp13*det*A33*A33))
  231. & +vard(7)*(-2.D0*detp23*det*A33*A33)/det4
  232.  
  233. DINV(3,4)=(COEFF1*(vard(2)*(-D23*A13*det2-2.D0*detp11*det*A13*A23)
  234. & +vard(3)*(-D13*A23*det2-2.D0*detp22*det*A13*A23)
  235. & +vard(4)*(-2.D0*detp33*det*A13*A23)
  236. & +vard(5)*((D23*A23+D13*A13)*det2-2.D0*detp12*det*A13*A23)
  237. & +vard(6)*(((DE2)*A23+D12*A13)*det2
  238. & -2.D0*detp13*det*A13*A23)
  239. & +vard(7)*((D12*A23+(DE1)*A13)*det2
  240. & -2.D0*detp23*det*A13*A23))
  241. & -COEFF2*(vard(2)*((D22-1.D0)*A12*det2
  242. & -2.D0*detp11*det*A33*A12)
  243. & +vard(3)*((D11-1.D0)*A12*det2-2.D0*detp22*det*A33*A12)
  244. & +vard(4)*(-D12*A33*det2-2.D0*detp33*det*A33*A12)
  245. & +vard(5)*(((DE3)*A33-2.D0*D12*A12)*det2
  246. & -2.D0*detp12*det*A33*A12)
  247. & +vard(6)*(D23*A33*det2-2.D0*detp13*det*A33*A12)
  248. & +vard(7)*(D13*A33*det2-2.D0*detp23*det*A33*A12)))*2.D0/det4
  249.  
  250. IF (IFOUR.GT.0) THEN
  251. DINV(3,5)=2.D0*COEFF3*(vard(2)*((D22-1.D0)*A13*det2
  252. & -2.D0*detp11*det*A33*A13)
  253. & +vard(3)*(((D11-1.D0)*A13-D13*A33)*det2
  254. & -2.D0*detp22*det*A33*A13)
  255. & +vard(4)*(-2.D0*detp33*det*A33*A13)
  256. & +vard(5)*((-2.D0*D12*A13+D23*A33)*det2
  257. & -2.D0*detp12*det*A33*A13)
  258. & +vard(6)*((DE2)*A33*det2-2.D0*detp13*det*A33*A13)
  259. & +vard(7)*(D12*A33*det2-2.D0*detp23*det*A33*A13))/det4
  260.  
  261. DINV(3,6)=2.D0*COEFF3*(vard(2)*(((D22-1.D0)*A23-D23*A33)*det2
  262. & -2.D0*detp11*det*A33*A23)
  263. & +vard(3)*((D11-1.D0)*A23*det2-2.D0*detp22*det*A33*A23)
  264. & +vard(4)*(-2.D0*detp33*det*A33*A23)
  265. & +vard(5)*((-2.D0*D12*A23+D13*A33)*det2
  266. & -2.D0*detp12*det*A33*A23)
  267. & +vard(6)*(D12*A33*det2-2.D0*detp13*det*A33*A23)
  268. & +vard(7)*((DE1)*A33*det2
  269. & -2.D0*detp23*det*A33*A23))/det4
  270. ENDIF
  271.  
  272. DINV(4,4)=(2.D0/det4)*(COEFF1*(vard(2)*((D33-1.D0)*A11*det2
  273. & -2.D0*detp11*det*(A12*A12+A11*A22))
  274. & +vard(3)*((D33-1.D0)*A22*det2
  275. & -2.D0*detp22*det*(A12*A12+A11*A22))
  276. & +vard(4)*((-2.D0*D12*A12+(D22-1.D0)*A22+(D11-1.D0)*A11)
  277. & *det2-2.D0*detp33*det*(A12*A12+A11*A22))
  278. & +vard(5)*(2.D0*(DE3)*A12*det2
  279. & -2.D0*detp12*det*(A12*A12+A11*A22))
  280. & +vard(6)*((2.D0*D23*A12-2.D0*D13*A11)*det2
  281. & -2.D0*detp13*det*(A12*A12+A11*A22))
  282. & +vard(7)*((2.D0*D13*A12-2.D0*D23*A22)*det2
  283. & -2.D0*detp23*det*(A12*A12+A11*A22)))
  284. & -2.D0*COEFF2*(vard(2)*(-2.D0*detp11*det*A12*A12)
  285. & +vard(3)*(-2.D0*detp22*det*A12*A12)
  286. & +vard(4)*(-2.D0*D12*A12*det2-2.D0*detp33*det*A12*A12)
  287. & +vard(5)*(2.D0*(DE3)*A12*det2-2.D0*detp12*det*A12*A12)
  288. & +vard(6)*(2.D0*D23*A12*det2-2.D0*detp13*det*A12*A12)
  289. & +vard(7)*(2.D0*D13*A12*det2-2.D0*detp23*det*A12*A12)))
  290.  
  291. IF (IFOUR.GT.0) THEN
  292. DINV(4,5)=(2.D0/det4)*(COEFF1*(vard(2)*(-D23*A11*det2
  293. & -2.D0*detp11*det*(A13*A12+A11*A23))
  294. & +vard(3)*((-D13*A12+(D33-1.D0)*A23)*det2
  295. & -2.D0*detp22*det*(A13*A12+A11*A23))
  296. & +vard(4)*((-D12*A13+(D22-1.D0)*A23)*det2
  297. & -2.D0*detp33*det*(A13*A12+A11*A23))
  298. & +vard(5)*((D23*A12+(DE3)*A13+D13*A11)*det2
  299. & -2.D0*detp12*det*(A13*A12+A11*A23))
  300. & +vard(6)*(((DE2)*A12+D23*A13+D12*A11)*det2
  301. & -2.D0*detp13*det*(A13*A12+A11*A23))
  302. & +vard(7)*((D12*A12+D13*A13-2.D0*D23*A23+(DE1)*A11)
  303. & *det2-2.D0*detp23*det*(A13*A12+A11*A23)))
  304. & -2.D0*COEFF2*(vard(2)*(-2.D0*detp11*det*A12*A13)
  305. & +vard(3)*(-D13*A12*det2-2.D0*detp22*det*A12*A13)
  306. & +vard(4)*(-D12*A13*det2-2.D0*detp33*det*A12*A13)
  307. & +vard(5)*(((DE3)*A13+D23*A12)*det2
  308. & -2.D0*detp12*det*A12*A13) +vard(6)*((D23*A13+(DE2)
  309. & *A12)*det2-2.D0*detp13*det*A12*A13)
  310. & +vard(7)*((D13*A13+D12*A12)*det2-2.D0*detp23*det*A12*A13)))
  311.  
  312. DINV(4,6)=(2.D0/det4)*(COEFF1*(vard(2)*((-D23*A12
  313. & +(D33-1.D0)*A13)*det2-2.D0*detp11*det*(A23*A12+A22*A13))
  314. & +vard(3)*(-D13*A22*det2-2.D0*detp22*det*(A23*A12+A22*A13))
  315. & +vard(4)*((-D12*A23+(D11-1.D0)*A13)*det2
  316. & -2.D0*detp33*det*(A23*A12+A22*A13))
  317. & +vard(5)*((D13*A12+(DE3)*A23+D23*A22)*det2
  318. & -2.D0*detp12*det*(A23*A12+A22*A13))
  319. & +vard(6)*((D12*A12+D23*A23-2.D0*D13*A13+(DE2)*A22)
  320. & *det2-2.D0*detp13*det*(A23*A12+A22*A13))
  321. & +vard(7)*(((DE1)*A12+D13*A23+D12*A22)*det2
  322. & -2.D0*detp23*det*(A23*A12+A22*A13)))
  323. & -2.D0*COEFF2*(vard(2)*(-D23*A12*det2-2.D0*detp11*det*A12
  324. & *A23) + vard(3)*(-2.D0*detp22*det*A12*A23)
  325. & +vard(4)*(-D12*A23*det2-2.D0*detp33*det*A12*A23)
  326. & +vard(5)*(((DE3)*A23+D13*A12)*det2
  327. & -2.D0*detp12*det*A12*A23)
  328. & +vard(6)*((D23*A23+D12*A12)*det2-2.D0*detp13*det*A12*A23)
  329. & +vard(7)*((D13*A23+(DE1)*A12)*det2
  330. & -2.D0*detp23*det*A12*A23)))
  331.  
  332. DINV(5,5)=(2.D0/det4)*(COEFF1*(vard(2)*((D22-1.D0)*A11*det2
  333. & -2.D0*detp11*det*(A13*A13+A11*A33))
  334. & +vard(3)*((-2.D0*D13*A13+(D33-1.D0)*A33+(D11-1.D0)*A11)*det2
  335. & -2.D0*detp22*det*(A13*A13+A11*A33))
  336. & +vard(4)*((D22-1.D0)*A33*det2
  337. & -2.D0*detp33*det*(A13*A13+A11*A33))
  338. & +vard(5)*((2.D0*D23*A13-2.D0*D12*A11)*det2
  339. & -2.D0*detp12*det*(A13*A13+A11*A33))
  340. & +vard(6)*(2.D0*(DE2)*A13*det2
  341. & -2.D0*detp13*det*(A13*A13+A11*A33))
  342. & +vard(7)*((2.D0*D12*A13-2.D0*D23*A33)*det2
  343. & -2.D0*detp23*det*(A13*A13+A11*A33)))
  344. & -2.D0*COEFF2*(vard(2)*(-2.D0*detp11*det*A13*A13)
  345. & +vard(3)*(-2.D0*D13*A13*det2-2.D0*detp22*det*A13*A13)
  346. & +vard(4)*(-2.D0*detp33*det*A13*A13)
  347. & +vard(5)*(2.D0*D23*A13*det2-2.D0*detp12*det*A13*A13)
  348. & +vard(6)*(2.D0*(DE2)*A13*det2-2.D0*detp13*det*A13*A13)
  349. & +vard(7)*(2.D0*D12*A13*det2-2.D0*detp23*det*A13*A13)))
  350.  
  351. DINV(5,6)=(2.D0/det4)*(COEFF1*(vard(2)*((-D23*A13
  352. & +(D22-1.D0)*A12)*det2-2.D0*detp11*det*(A23*A13+A33*A12))
  353. & +vard(3)*((-D13*A23+(D11-1.D0)*A12)*det2
  354. & -2.D0*detp22*det*(A23*A13+A33*A12))
  355. & +vard(4)*(-D12*A33*det2-2.D0*detp33*det*(A23*A13+A33*A12))
  356. & +vard(5)*((D13*A13+D23*A23-2.D0*D12*A12+(DE3)*A33)*det2
  357. & -2.D0*detp12*det*(A23*A13+A33*A12))
  358. & +vard(6)*((D12*A13+(DE2)*A23+D23*A33)*det2
  359. & -2.D0*detp13*det*(A23*A13+A33*A12))
  360. & +vard(7)*(((DE1)*A13+D12*A23+D13*A33)*det2
  361. & -2.D0*detp23*det*(A23*A13+A33*A12)))
  362. & -2.D0*COEFF2*(vard(2)*(-D23*A13*det2-2.D0*detp11*det*A13
  363. & *A23) +vard(3)*(-D13*A23*det2-2.D0*detp22*det*A12*A23)
  364. & +vard(4)*(-2.D0*detp33*det*A13*A23)
  365. & +vard(5)*((D23*A23+D13*A13)*det2-2.D0*detp12*det*A13*A23)
  366. & +vard(6)*(((DE2)*A23+D12*A13)*det2
  367. & -2.D0*detp13*det*A13*A23)
  368. & +vard(7)*((D12*A23+(DE1)*A13)*det2
  369. & -2.D0*detp23*det*A13*A23)))
  370.  
  371. DINV(6,6)=(2.D0/det4)*(COEFF1*(vard(2)*((-2.D0*D23*A23
  372. & +(D33-1.D0)*A33+(D22-1.D0)*A22)*det2
  373. & -2.D0*detp11*det*(A23*A23+A22*A33))
  374. & +vard(3)*((D11-1.D0)*A22*det2
  375. & -2.D0*detp22*det*(A23*A23+A22*A33))
  376. & +vard(4)*((D11-1.D0)*A33*det2
  377. & -2.D0*detp33*det*(A23*A23+A22*A33))
  378. & +vard(5)*((2.D0*D13*A23-2.D0*D12*A22)*det2
  379. & -2.D0*detp12*det*(A23*A23+A22*A33))
  380. & +vard(6)*((2.D0*D12*A23-2.D0*D13*A33)*det2
  381. & -2.D0*detp13*det*(A23*A23+A22*A33))
  382. & +vard(7)*(2.D0*(DE1)*A23*det2
  383. & -2.D0*detp23*det*(A23*A23+A22*A33)))
  384. & -2.D0*COEFF2*(vard(2)*(-2.D0*D23*A23*det2-2.D0
  385. & *detp11*det*A23*A23) +vard(3)*(-2.D0*detp22
  386. & *det*A23*A23) +vard(4)*(-2.D0*detp33*det*A23*A23)
  387. & +vard(5)*(2.D0*D13*A23*det2-2.D0*detp12*det*A23*A23)
  388. & +vard(6)*(2.D0*D12*A23*det2-2.D0*detp13*det*A23*A23)
  389. & +vard(7)*(2.D0*(DE1)*A23*det2
  390. & -2.D0*detp23*det*A23*A23)))
  391. ENDIF
  392.  
  393. DINV(2,1)=DINV(1,2)
  394.  
  395. DINV(3,1)=DINV(1,3)
  396. DINV(3,2)=DINV(2,3)
  397.  
  398. DINV(4,1)=DINV(1,4)
  399. DINV(4,2)=DINV(2,4)
  400. DINV(4,3)=DINV(3,4)
  401.  
  402. IF (IFOUR.GT.0) THEN
  403. DINV(5,1)=DINV(1,5)
  404. DINV(5,2)=DINV(2,5)
  405. DINV(5,3)=DINV(3,5)
  406. DINV(5,4)=DINV(4,5)
  407.  
  408. DINV(6,1)=DINV(1,6)
  409. DINV(6,2)=DINV(2,6)
  410. DINV(6,3)=DINV(3,6)
  411. DINV(6,4)=DINV(4,6)
  412. DINV(6,5)=DINV(5,6)
  413. ENDIF
  414.  
  415. RETURN
  416. END
  417.  
  418.  
  419.  
  420.  
  421.  
  422.  

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