Télécharger aactum.eso

Retour à la liste

Numérotation des lignes :

  1. C AACTUM SOURCE CB215821 16/04/21 21:15:00 8920
  2. SUBROUTINE AACTUM(SIGT,EPST,ENDT,XLIMT,SIGD,EPSD,ENDD,XLIMD,
  3. * VDEP,EPSC,EPSPLS,XLMD,PPLUS,
  4. * VDEP1,SIFF,VDEP2,FDDD,DEPS,DSIG,EDED,FD,
  5. * PROD,ENDOPL,EPL,E,EP,D,CONT3,CONT4,EPLUS,
  6. * EDPLUS,U,V,R,W1,W2)
  7. C===================================================================
  8. C
  9. C CALCUL DES GRANDEURS A LA FIN DU PAS
  10. C
  11. C====================================================================
  12. C
  13. C CREATION : F.CORMERY
  14. C E.N.S.M.A - LMPM
  15. C JUIN 1993
  16. C
  17. C====================================================================
  18. IMPLICIT INTEGER(I-N)
  19. IMPLICIT REAL*8 (A-H,O-Z)
  20. REAL*8 LAMB,MUT
  21. C-----------------------------------------------------------------------
  22. C COMMON SPECIFIQUE
  23. C-----------------------------------------------------------------------
  24. COMMON/APENTE/ LAMB,MUT,G,ALP1,BETA
  25. COMMON/MAED/ B,C0,CPRIME
  26. COMMON/DSIGEP/ H
  27. C-----------------------------------------------------------------------
  28. C DATA ET DIMENSION
  29. C-----------------------------------------------------------------------
  30. C N6
  31. CC DIMENSION SIGT(6),EPST(6),ENDT(6),SIFF(6),DEPS(6),
  32. C N12 N18
  33. CC * SIGD(6),EPSD(6),ENDD(6),DSIG(6),EDED(6),
  34. C N24
  35. CC * VDEP(6,6),VDEP1(6,6),FD(6),
  36. C N30 N36 N42
  37. CC * EPSC(6),EPSPLS(6),PROD(6),ENDOPL(6),EPL(6)
  38. C N51 N60 N69 N78 N87 N96
  39. CC DIMENSION E(3,3),EP(3,3),D(3,3),CONT3(3,3),CONT4(3,3),EPLUS(3,3)
  40. C N105
  41. CC * ,EDPLUS(3,3),VDEP2(6,6)
  42. C N186 N267 N348
  43. CC DIMENSION U(3,3,3,3),V(3,3,3,3),R(3,3,3,3),PPLUS(3,3,3,3),
  44. C N429 N510
  45. CC * W1(3,3,3,3),W2(3,3,3,3)
  46. C
  47. DIMENSION SIGT(*),EPST(*),ENDT(*),SIFF(*),DEPS(*),
  48. * SIGD(*),EPSD(*),ENDD(*),DSIG(*),EDED(*),
  49. * VDEP(6,*),VDEP1(6,*),FD(*),
  50. * EPSC(*),EPSPLS(*),PROD(*),ENDOPL(*),EPL(*)
  51. DIMENSION E(3,*),EP(3,*),D(3,*),CONT3(3,*),CONT4(3,*),EPLUS(3,*)
  52. * ,EDPLUS(3,*),VDEP2(6,*)
  53. DIMENSION U(3,3,3,*),V(3,3,3,*),R(3,3,3,*),PPLUS(3,3,3,*),
  54. * W1(3,3,3,*),W2(3,3,3,*)
  55. C--------------
  56. DATA ZERO/0.D0/,UN/1.D0/,DEUX/2.D0/
  57. DATA UD/0.5D0/
  58. C-----------------------------------------------------------------------
  59. C CALCUL DES TRACES
  60. C-----------------------------------------------------------------------
  61. TRAC1=EPSPLS(1)**2+EPSPLS(2)**2+EPSPLS(3)**2
  62. * +UD*(EPSPLS(4)**2+EPSPLS(5)**2+EPSPLS(6)**2)
  63. C----------------
  64. TRAC3=EPSD(1)+EPSD(2)+EPSD(3)
  65. C----------------
  66. TRAC4=EPSPLS(1)*EPSD(1)+EPSPLS(2)*EPSD(2)+EPSPLS(3)*EPSD(3)
  67. * +UD*(EPSPLS(4)*EPSD(4)+EPSPLS(5)*EPSD(5)
  68. * +EPSPLS(6)*EPSD(6))
  69. C----------------
  70. TRAC5=EPSPLS(1)+EPSPLS(2)+EPSPLS(3)
  71. C************************************************************************
  72. C ACTUALISATION DES GRANDEURS
  73. C************************************************************************
  74. C
  75. C------------------------------------------------------------------------
  76. C CALCUL DE D n+1
  77. C------------------------------------------------------------------------
  78. DO 10 I=1,6
  79. IF(I.LE.3) ENDD(I)=ENDT(I)+(XLMD/(UN-B*XLMD))*
  80. * (EPSPLS(I)/SQRT(DEUX*TRAC1)+B*ENDT(I))
  81. IF(I.GE.4) ENDD(I)=ENDT(I)+(XLMD/(UN-B*XLMD))*
  82. * ((EPSPLS(I)/DEUX)/SQRT(DEUX*TRAC1)+B*ENDT(I))
  83. IF(ABS(ENDD(I)).LE.1E-15) ENDD(I)=0.D0
  84. 10 CONTINUE
  85. C------------------------------------------------------------------------
  86. C CALCUL DES FORCES THERMO
  87. C------------------------------------------------------------------------
  88. C---------CALCUL DU PRODUIT EPS.EPS
  89. EDED(1)=EPSD(1)**2+EPSD(6)**2+EPSD(5)**2
  90. EDED(2)=EPSD(2)**2+EPSD(6)**2+EPSD(4)**2
  91. EDED(1)=EPSD(3)**2+EPSD(4)**2+EPSD(5)**2
  92. EDED(4)=EPSD(6)*EPSD(5)+EPSD(2)*EPSD(4)+EPSD(4)*EPSD(3)
  93. EDED(5)=EPSD(1)*EPSD(5)+EPSD(6)*EPSD(4)+EPSD(5)*EPSD(3)
  94. EDED(6)=EPSD(6)*EPSD(1)+EPSD(2)*EPSD(6)+EPSD(4)*EPSD(5)
  95. DO 11 I=1,6
  96. FD(I)=-G*EPSD(I)-ALP1*TRAC3*EPSD(I)-2*BETA*EDED(I)
  97. 11 CONTINUE
  98. C------------------------------------------------------------------------
  99. C CALCUL DE FD:DD
  100. C------------------------------------------------------------------------
  101. FDDD=ZERO
  102. DO 12 I=1,6
  103. FDDD=FDDD+FD(I)*(ENDD(I)-ENDT(I))
  104. 12 CONTINUE
  105. C------------------------------------------------------------------------
  106. C CALCUL DES TRACES
  107. C------------------------------------------------------------------------
  108. TRAC2=EPSPLS(1)*ENDD(1)+EPSPLS(2)*ENDD(2)+EPSPLS(3)*ENDD(3)
  109. * +(EPSPLS(4)*ENDD(4)+EPSPLS(5)*ENDD(5)
  110. * +EPSPLS(6)*ENDD(6))
  111. TRAC6=ENDD(1)+ENDD(2)+ENDD(3)
  112. C------------------------------------------------------------------------
  113. C CALCUL DE C0+C1 tr(Dn+1)
  114. C------------------------------------------------------------------------
  115. XLIMD=C0+CPRIME*(ENDD(1)+ENDD(2)+ENDD(3))
  116. C------------------------------------------------------------------------
  117. C CALCUL DE D+ n+1
  118. C------------------------------------------------------------------------
  119. MOD1=0
  120. IF(MOD1.EQ.1)THEN
  121. CALL T4CT2(PPLUS,ENDD,ENDOPL,6)
  122. GOTO 61
  123. ENDIF
  124. DO 62 I=1,6
  125. ENDOPL(I)=ENDD(I)
  126. 62 CONTINUE
  127. 61 CONTINUE
  128. C------------------------------------------------------------------------
  129. C
  130. C-------------------------------------------------------------------------
  131. DO 13 I=4,6
  132. EPSPLS(I)=EPSPLS(I)/2
  133. 13 CONTINUE
  134. CALL T4CT2(PPLUS,EPSPLS,EPL,6)
  135. DO 14 I=4,6
  136. EPSPLS(I)=EPSPLS(I)*2
  137. 14 CONTINUE
  138. C-------------------------------------------------------------------------
  139. C CALCUL DE SIG n+1
  140. C-------------------------------------------------------------------------
  141. TRACE2=EPSD(1)*ENDD(1)+EPSD(2)*ENDD(2)+EPSD(3)*ENDD(3)+
  142. * EPSD(4)*ENDD(4)+EPSD(5)*ENDD(5)+EPSD(6)*ENDD(6)
  143. C-------
  144. PROD(1)=DEUX*(EPSD(1)*ENDD(1))+EPSD(5)*ENDD(5)+EPSD(6)*ENDD(6)
  145. PROD(2)=DEUX*(EPSD(2)*ENDD(2))+EPSD(4)*ENDD(4)+EPSD(6)*ENDD(6)
  146. PROD(3)=DEUX*(EPSD(3)*ENDD(3))+EPSD(5)*ENDD(5)+EPSD(4)*ENDD(4)
  147. PROD(4)=UD*EPSD(6)*ENDD(5)+EPSD(2)*ENDD(4)+UD*EPSD(4)*ENDD(3)+
  148. * UD*EPSD(5)*ENDD(6)+UD*EPSD(4)*ENDD(2)+EPSD(3)*ENDD(4)
  149. PROD(5)=UD*EPSD(5)*ENDD(1)+UD*EPSD(4)*ENDD(6)+EPSD(3)*ENDD(5)+
  150. * EPSD(1)*ENDD(5)+UD*EPSD(6)*ENDD(4)+UD*EPSD(5)*ENDD(3)
  151. PROD(6)=UD*EPSD(6)*ENDD(1)+EPSD(2)*ENDD(6)+UD*EPSD(4)*ENDD(5)+
  152. * EPSD(1)*ENDD(6)+UD*EPSD(6)*ENDD(2)+UD*EPSD(5)*ENDD(4)
  153. C-------
  154. NP=1
  155. IF(NP.EQ.1)THEN
  156. DO 30 I=1,3
  157. SIGD(I)=G*ENDOPL(I)+LAMB*TRAC3+ALP1*(TRACE2+TRAC3*ENDD(I))
  158. * +DEUX*BETA*PROD(I)+DEUX*MUT*EPSD(I)
  159. 30 CONTINUE
  160. DO 35 I=4,6
  161. SIGD(I)=G*ENDOPL(I)+ALP1*TRAC3*ENDD(I)+MUT*EPSD(I)
  162. * +DEUX*BETA*PROD(I)
  163. 35 CONTINUE
  164. GOTO 21
  165. ENDIF
  166. C------------------------------------------------------------------
  167. C CALCUL MATRICE TANGENTE continue
  168. C ET MATRICE TANGENTE consistante
  169. C------------------------------------------------------------------
  170. C
  171. C----------------------------------Calcul des coef.----------------
  172. 21 XA=DEUX*TRAC1
  173. XM=B*ABS(G)*SQRT(TRAC1/DEUX)+(B**2)*ABS(G)*TRAC2
  174. * -CPRIME*TRAC5/SQRT(DEUX*TRAC1)-CPRIME*B*TRAC6
  175. V2=-XLMD/(SQRT(XA)*TRAC5)
  176. V3=XLMD/SQRT(XA)
  177. V5=ABS(G)/(SQRT(XA)*CPRIME*TRAC5)
  178. C---------------------------------Passage en notation indicielle---
  179. DO 42 I=1,3
  180. D(I,I)=ENDD(I)
  181. EP(I,I)=EPSPLS(I)
  182. E(I,I)=EPSD(I)
  183. EDPLUS(I,I)=ENDOPL(I)
  184. EPLUS(I,I)=EPL(I)
  185. DO 43 J=I+1,3
  186. IF ((I.EQ.1).AND.(J.EQ.2)) K=6
  187. IF ((I.EQ.1).AND.(J.EQ.3)) K=5
  188. IF ((I.EQ.2).AND.(J.EQ.3)) K=4
  189. D(I,J)=ENDD(K)
  190. D(J,I)=ENDD(K)
  191. E(I,J)=UD*EPSD(K)
  192. E(J,I)=UD*EPSD(K)
  193. EP(I,J)=UD*EPSPLS(K)
  194. EP(J,I)=UD*EPSPLS(K)
  195. EDPLUS(I,J)=ENDOPL(K)
  196. EDPLUS(J,I)=ENDOPL(K)
  197. EPLUS(I,J)=EPL(K)
  198. EPLUS(J,I)=EPL(K)
  199. 43 CONTINUE
  200. 42 CONTINUE
  201.  
  202. C--------------------- Calcul des produit contractes -----------------
  203. CONT1=ZERO
  204. CONT2=ZERO
  205. DO 44 I=1,3
  206. DO 45 J=1,3
  207. CONT1=CONT1+E(I,J)*EP(I,J)
  208. CONT2=CONT2+E(I,J)*D(I,J)
  209. CONT3(I,J)=ZERO
  210. CONT4(I,J)=ZERO
  211. DO 46 K=1,3
  212. CONT3(I,J)=CONT3(I,J)+E(I,K)*EP(K,J)+EP(I,K)*E(K,J)
  213. CONT4(I,J)=CONT4(I,J)+E(I,K)*D(K,J)+D(I,K)*E(K,J)
  214. 46 CONTINUE
  215. 45 CONTINUE
  216. 44 CONTINUE
  217.  
  218. C----------------------------Calcul de la matrice tangente---------
  219. DO 50 I=1,3
  220. DO 60 J=1,3
  221. DO 70 K=1,3
  222. DO 80 L=1,3
  223. C-------------------------
  224. U(I,J,K,L)=LAMB*DELT(I,J)*DELT(K,L)+MUT*(DELT(I,K)*DELT(J,L)
  225. * +DELT(I,L)*DELT(J,K))+ALP1*(DELT(I,J)*D(K,L)+D(I,J)*DELT(K,L))
  226. * +BETA*(D(I,K)*DELT(J,L)+D(L,J)*DELT(I,K)+DELT(I,L)
  227. * *D(J,K)+D(I,L)*DELT(J,K))
  228. C-------------------------
  229. V(I,J,K,L)=(1/XM)*(ALP1*G*CONT1*DELT(I,J)/XA+ALP1*B*G
  230. * *CONT2*DELT(I,J)/SQRT(XA)+(ALP1*G/XA)*EP(I,J)
  231. * *TRAC3+(ALP1*B*G/SQRT(XA))*D(I,J)*TRAC3+(DEUX*BETA*G/XA)
  232. * *CONT3(I,J)+DEUX*BETA*(B*G/SQRT(XA))*CONT4(I,J)
  233. * +(G**2/XA)*EPLUS(I,J)+(B*(G**2)/SQRT(XA))
  234. * *EDPLUS(I,J))*EPLUS(K,L)
  235. C-------------------------
  236. R(I,J,K,L)=(1/XM)*((ALP1*B*G/SQRT(XA))*CONT1*DELT(I,J)
  237. * +ALP1*(B**2)*G*CONT2*DELT(I,J)+(ALP1*B*G/SQRT(XA))
  238. * *EP(I,J)*TRAC3+ALP1*(B**2)*G*D(I,J)*TRAC3+DEUX*BETA
  239. * *(B*G/SQRT(XA))*CONT3(I,J)+DEUX*BETA*(B**2)*G*
  240. * CONT4(I,J)+(B*(G**2)/SQRT(XA))*EPLUS(I,J)
  241. * +B**2*G**2*EDPLUS(I,J))*EDPLUS(K,L)
  242. C-------------------------
  243. W1(I,J,K,L)=0.D0
  244. W2(I,J,K,L)=0.D0
  245. DO 81 M=1,3
  246. DO 81 N=1,3
  247. W1(I,J,K,L)=(G*V5*EPLUS(I,J)*E(M,N)+G*V2*EPLUS(I,J)*DELT(M,N)
  248. * +G*V3*PPLUS(I,J,M,N)+DEUX*BETA*V5*CONT3(I,J)*E(M,N)
  249. * +DEUX*BETA*V2*CONT3(I,J)*DELT(M,N)+BETA*V3*
  250. * (E(I,M)*DELT(J,N)+DELT(J,M)*E(I,N)+DELT(I,M)*E(J,N)+
  251. * E(J,M)*DELT(I,N))+ALP1*V5*TRAC4*DELT(I,J)*EP(M,N)
  252. * +ALP1*V2*TRAC4*DELT(I,J)*DELT(M,N)+ALP1*V3*
  253. * DELT(I,J)*E(M,N)+ALP1*V5*TRAC3*EP(I,J)*EP(M,N)
  254. * )*PPLUS(M,N,K,L)+W1(I,J,K,L)
  255. W2(I,J,K,L)=(ALP1*V2*TRAC3*EP(I,J)*DELT(M,N)
  256. * +ALP1*V3*TRAC3*(DELT(I,M)*DELT(N,J)+DELT(J,M)*DELT(N,I))
  257. * /2.D0)*PPLUS(M,N,K,L)+W2(I,J,K,L)
  258. 81 CONTINUE
  259. C-------------------------
  260. 80 CONTINUE
  261. 70 CONTINUE
  262. 60 CONTINUE
  263. 50 CONTINUE
  264. C--------------------------
  265. DO 110 I=1,6
  266. DO 120 J=1,6
  267. IF (I.EQ.1) THEN
  268. M=1
  269. N=1
  270. ENDIF
  271. IF (I.EQ.2) THEN
  272. M=2
  273. N=2
  274. ENDIF
  275. IF (I.EQ.3) THEN
  276. M=3
  277. N=3
  278. ENDIF
  279. IF (I.EQ.4) THEN
  280. M=1
  281. N=2
  282. ENDIF
  283. IF (I.EQ.5) THEN
  284. M=2
  285. N=3
  286. ENDIF
  287. IF (I.EQ.6) THEN
  288. M=1
  289. N=3
  290. ENDIF
  291. IF (J.EQ.1) THEN
  292. K=1
  293. L=1
  294. ENDIF
  295. IF (J.EQ.2) THEN
  296. K=2
  297. L=2
  298. ENDIF
  299. IF (J.EQ.3) THEN
  300. K=3
  301. L=3
  302. ENDIF
  303. IF (J.EQ.4) THEN
  304. K=1
  305. L=2
  306. ENDIF
  307. IF (J.EQ.5) THEN
  308. K=2
  309. L=3
  310. ENDIF
  311. IF (J.EQ.6) THEN
  312. K=1
  313. L=3
  314. ENDIF
  315. VDEP2(I,J)=U(M,N,K,L)
  316. VDEP(I,J)=U(M,N,K,L)+V(M,N,K,L)+R(M,N,K,L)
  317. VDEP1(I,J)=U(M,N,K,L)+W1(M,N,K,L)+W2(M,N,K,L)
  318. IF(ABS(VDEP2(I,J)).LE.1E-15)VDEP2(I,J)=0.D0
  319. IF(ABS(VDEP1(I,J)).LE.1E-15)VDEP1(I,J)=0.D0
  320. IF(ABS(VDEP(I,J)).LE.1E-15)VDEP(I,J)=0.D0
  321. 120 CONTINUE
  322. 110 CONTINUE
  323. C---------------------------------------------------------------------
  324. C CALCUL DE SIG n+1
  325. C---------------------------------------------------------------------
  326. DO 111 I=1,6
  327. DO 121 J=1,6
  328. IF (I.EQ.1) THEN
  329. M=1
  330. N=1
  331. ENDIF
  332. IF (I.EQ.2) THEN
  333. M=2
  334. N=2
  335. ENDIF
  336. IF (I.EQ.3) THEN
  337. M=3
  338. N=3
  339. ENDIF
  340. IF (I.EQ.4) THEN
  341. M=2
  342. N=3
  343. ENDIF
  344. IF (I.EQ.5) THEN
  345. M=1
  346. N=3
  347. ENDIF
  348. IF (I.EQ.6) THEN
  349. M=1
  350. N=2
  351. ENDIF
  352. IF (J.EQ.1) THEN
  353. K=1
  354. L=1
  355. ENDIF
  356. IF (J.EQ.2) THEN
  357. K=2
  358. L=2
  359. ENDIF
  360. IF (J.EQ.3) THEN
  361. K=3
  362. L=3
  363. ENDIF
  364. IF (J.EQ.4) THEN
  365. K=2
  366. L=3
  367. ENDIF
  368. IF (J.EQ.5) THEN
  369. K=1
  370. L=3
  371. ENDIF
  372. IF (J.EQ.6) THEN
  373. K=1
  374. L=2
  375. ENDIF
  376. VDEP2(I,J)=U(M,N,K,L)
  377. 121 CONTINUE
  378. 111 CONTINUE
  379. C-------------------------------------------------------------
  380. NP=2
  381. IF(NP.EQ.2)THEN
  382. DO 302 I=1,6
  383. SIFF(I)=0.D0
  384. DO 301 J=1,6
  385. SIFF(I)=VDEP2(I,J)*EPSD(J)+SIFF(I)
  386. 301 CONTINUE
  387. SIFF(I)=G*ENDOPL(I)+SIFF(I)
  388. 302 CONTINUE
  389. ENDIF
  390. C-------------------------------------------------------------
  391. H=0
  392. DO 307 I=1,6
  393. DEPS(I)=EPSD(I)-EPST(I)
  394. DSIG(I)=SIGD(I)-SIGT(I)
  395. H=H+DEPS(I)*DSIG(I)
  396. 307 CONTINUE
  397. C WRITE(10,1203)H
  398. C1203 FORMAT(/5X,'DSIG : DEPS = ',E12.5)
  399. RETURN
  400. END
  401.  
  402.  
  403.  
  404.  
  405.  
  406.  

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