Télécharger ricbet2d.eso

Retour à la liste

Numérotation des lignes :

  1. C RICBET2D SOURCE BP208322 17/03/01 21:18:07 9325
  2. SUBROUTINE RICBET2D(XMAT,SIG0,SIGF,VAR0,VARF,DEPST,ISTEP)
  3. C
  4. C====&===1=========2=========3=========4=========5=========6=========7==
  5. C Commentaires : Subroutine permettant de mettre en oeuvre le
  6. C modele RICBET (B. Richard) pour representer
  7. C le comportement 2D d'un béton ordinaire sous
  8. C sollicitations cycliques
  9. C
  10. C Traits : - Endommagement scalaire
  11. C - Glissement frottant
  12. C - 2D
  13. C - Monotone et cyclique
  14. C - Effet unilatteral
  15. C - Comportement linŽaire en compression
  16. C - Non localite portant sur le seuil d'endommagement (Fd)
  17. C
  18. C Auteurs : B. Richard (Dr - Ing.) & F. Ragueneau (Pr) - LMT/ENS
  19. C Date : Mars 2011
  20. C====&===1=========2=========3=========4=========5=========6=========7==
  21. C
  22. C-----DECLARATION GENERALE----------------------------------------------
  23. C
  24. IMPLICIT REAL*8(A-H,O-Z)
  25. C
  26. C-----LISTE DES INCLUDES A CHARGER--------------------------------------
  27. C
  28. C -INC CCOPTIO
  29. C -INC DECHE
  30. C
  31. C-----DECLARATION GENERALE----------------------------------------------
  32. C
  33. REAL*8 X(3,3),UNIT(3,3),EPSPI(3,3),EPSI(3,3),SIGMA(3,3)
  34. REAL*8 EPSDEV(3,3),SIGMA1(3,3),EPSPIDEV(3,3)
  35. REAL*8 EPSIPOS(3,3),EPSINEG(3,3),SIGPI(3,3),DEPSN(3,3)
  36. REAL*8 TERME1(3,3),TERME2(3,3),TERME3(3,3),TERME4(3,3)
  37. REAL*8 SIGMOX(3,3),TRAV(3,3),SIG1DEV(3,3),TERME5(3,3)
  38. REAL*8 XDT,XTRAV2(3,3),XTRAV3(3,3),DEFN(3,3),EPSIT(3,3)
  39. REAL*8 XMAT(*),SIG0(*),SIGF(*),VAR0(*),VARF(*),DEPST(*)
  40. REAL*8 XDFP(3,3),XDGP(3,3),XDJ2P(3,3),DEVGP(3,3)
  41. REAL*8 DEPSI(3,3)
  42. C
  43. C-----DEFINITION DE PARAMETRES------------------------------------------
  44. C
  45. PARAMETER (IMAX = 1000)
  46. PARAMETER (IMAX2 = 10)
  47. C
  48. C-----MISE EN DONNEE----------------------------------------------------
  49. C
  50. XE = XMAT(1)
  51. XNU = XMAT(2)
  52. XFT = XMAT(5)
  53. ALDIR = XMAT(6)
  54. GAMMA1 = XMAT(7)
  55. A1 = XMAT(8)
  56. SIGRF = XMAT(9)
  57.  
  58. XAF = XMAT(10)
  59. XBF = XMAT(11)
  60. XAG = XMAT(12)
  61. XBG = XMAT(13)
  62. XAC = XMAT(14)
  63. XBC = XMAT(15)
  64. XSIGU = XMAT(16)
  65. XFC = XMAT(17)
  66.  
  67. XK = XE/((1.0D0-2.D0*XNU)*3.0D0)
  68. XG = XE/(2.0D0*(1.D0+XNU))
  69. C
  70. C-----DEFINITION DE LA MATRICE UNITE D ORDRE 2--------------------------
  71. C
  72. DO I=1,3
  73. DO J=1,3
  74. IF (I.EQ.J) THEN
  75. UNIT(I,J) = 1.0D0
  76. ELSE
  77. UNIT(I,J) = 0.0D0
  78. ENDIF
  79. ENDDO
  80. ENDDO
  81. C
  82. C-----ACTUALISATION DES DEFORMATIONS------------------------------------
  83. C
  84.  
  85. C Indicateur
  86. XBEN = VAR0(22)+VAR0(23)+VAR0(24)*0.0D0
  87.  
  88. C Boucle pour les contraintes planes
  89. DO NIDX=1,IMAX2
  90.  
  91. IF (NIDX.EQ.IMAX) THEN
  92. C print*,'non converge'
  93. C stop
  94. ENDIF
  95.  
  96. IF (ABS(DEPST(3)).GE.1.0D10) THEN
  97. C print*,'pb depst'
  98. ENDIF
  99.  
  100. C Contraintes tests
  101. SIGMA1(1,1) = VAR0(22)
  102. SIGMA1(2,2) = VAR0(23)
  103. SIGMA1(3,3) = VAR0(24)
  104. SIGMA1(1,2) = VAR0(25)
  105. SIGMA1(1,3) = VAR0(26)
  106. SIGMA1(2,3) = VAR0(27)
  107. SIGMA1(2,1) = VAR0(25)
  108. SIGMA1(3,1) = VAR0(26)
  109. SIGMA1(3,2) = VAR0(27)
  110.  
  111. C Tenseur de glissement interne
  112. EPSPI(1,1) = VAR0(2)
  113. EPSPI(2,2) = VAR0(3)
  114. EPSPI(3,3) = VAR0(4)
  115. EPSPI(1,2) = VAR0(5)
  116. EPSPI(1,3) = VAR0(6)
  117. EPSPI(2,3) = VAR0(7)
  118. EPSPI(2,1) = VAR0(5)
  119. EPSPI(3,1) = VAR0(6)
  120. EPSPI(3,2) = VAR0(7)
  121.  
  122. C Tenseur des contraintes de rappel
  123. X(1,1) = VAR0(8)
  124. X(2,2) = VAR0(9)
  125. X(3,3) = VAR0(10)
  126. X(1,2) = VAR0(11)
  127. X(1,3) = VAR0(12)
  128. X(2,3) = VAR0(13)
  129. X(2,1) = VAR0(11)
  130. X(3,1) = VAR0(12)
  131. X(3,2) = VAR0(13)
  132.  
  133. C Tenseur des deformations plastiques
  134. DEFN(1,1) = VAR0(35)
  135. DEFN(2,2) = VAR0(36)
  136. DEFN(3,3) = VAR0(37)
  137. DEFN(1,2) = VAR0(38)
  138. DEFN(1,3) = VAR0(39)
  139. DEFN(2,3) = VAR0(40)
  140. DEFN(2,1) = VAR0(38)
  141. DEFN(3,1) = VAR0(39)
  142. DEFN(3,2) = VAR0(40)
  143.  
  144. C Ecrouissage plastique
  145. XECO = VAR0(41)
  146.  
  147. C Deformation plastique cumulee
  148. XDPP = VAR0(42)
  149.  
  150. C Endommagement
  151. XDT = VAR0(14)
  152.  
  153. C Force associee a l ecrouissage isotrope
  154. XFFT = VAR0(15)
  155.  
  156. C Deformations totales
  157. EPSI(1,1) = VAR0(16)
  158. EPSI(2,2) = VAR0(17)
  159. EPSI(3,3) = VAR0(18)
  160. EPSI(1,2) = VAR0(19)
  161. EPSI(1,3) = VAR0(20)
  162. EPSI(2,3) = VAR0(21)
  163. EPSI(2,1) = VAR0(19)
  164. EPSI(3,1) = VAR0(20)
  165. EPSI(3,2) = VAR0(21)
  166.  
  167. EPSI(1,1) = EPSI(1,1) + DEPST(1)
  168. EPSI(2,2) = EPSI(2,2) + DEPST(2)
  169. EPSI(3,3) = EPSI(3,3) + DEPST(3)
  170. EPSI(1,2) = EPSI(1,2) + (DEPST(4)/2.0D0)
  171. EPSI(1,3) = 0.0D0
  172. EPSI(2,3) = 0.0D0
  173. EPSI(2,1) = EPSI(2,1) + (DEPST(4)/2.0D0)
  174. EPSI(3,1) = 0.0D0
  175. EPSI(3,2) = 0.0D0
  176.  
  177. DEPSI(1,1) = DEPST(1)
  178. DEPSI(2,2) = DEPST(2)
  179. DEPSI(3,3) = DEPST(3)
  180. DEPSI(1,2) = (DEPST(4)/2.0D0)
  181. DEPSI(1,3) = 0.0D0
  182. DEPSI(2,3) = 0.0D0
  183. DEPSI(2,1) = (DEPST(4)/2.0D0)
  184. DEPSI(3,1) = 0.0D0
  185. DEPSI(3,2) = 0.0D0
  186. C
  187. C-----CALCUL DE QUELQUES QUANTITES--------------------------------------
  188. C
  189.  
  190. C Calcul de la trace de EPSI
  191. TRA = EPSI(1,1)+EPSI(2,2)+EPSI(3,3)
  192.  
  193. C Calcul du deviateur de EPSI
  194. CALL DEVIAT(EPSI,EPSDEV)
  195.  
  196. C Calcul de la trace de EPSPI
  197. TRAPI = EPSPI(1,1)+EPSPI(2,2)+EPSPI(3,3)
  198.  
  199. C Calcul du deviateur de EPSPI
  200. CALL DEVIAT(EPSPI,EPSPIDEV)
  201.  
  202. C Calcul de la trace de SIGMA1
  203. TRAS1 = SIGMA1(1,1)+SIGMA1(2,2)+SIGMA1(3,3)
  204.  
  205. C Calcul du deviateur de DEFN
  206. CALL DEVIAT(DEFN,DEPSN)
  207.  
  208. C Calcul de la trace de DEFN
  209. TRADN = DEFN(1,1)+DEFN(2,2)+DEFN(3,3)
  210.  
  211. C
  212. C-----PARTIE POSITIVES ET NEGATIVES DE (EPS-EPSP)------------------------
  213. C
  214. DO I=1,3
  215. DO J=1,3
  216. EPSIT(I,J) = EPSI(I,J)-DEFN(I,J)
  217. ENDDO
  218. ENDDO
  219.  
  220. CALL POSNEG(EPSIT,EPSIPOS,EPSINEG)
  221. C
  222. C-----SEUIL ENDOMMAGEMENT------------------------------------------------
  223. C
  224.  
  225. C Calcul de l enerigie Y+
  226. YP = 0.0D0
  227.  
  228. DO I=1,3
  229. DO J=1,3
  230. YP = YP + EPSIPOS(I,J)*EPSIPOS(I,J)
  231. ENDDO
  232. ENDDO
  233.  
  234. YP = 0.5D0*YP*XE
  235. C
  236. C-----NON LOCAL----------------------------------------------------------
  237. C
  238. IF (ISTEP.EQ.0) THEN
  239. SEUIL1=YP
  240. VARF(1)=YP
  241. ELSE IF (ISTEP.EQ.1) THEN
  242. VARF(1)=YP
  243. DO I=2,42
  244. VARF(I)=VAR0(I)
  245. ENDDO
  246. GOTO 2000
  247. ELSE IF (ISTEP.EQ.2) THEN
  248. SEUIL1=VAR0(1)
  249. VARF(1)=SEUIL1
  250. ENDIF
  251. C
  252. C-----EVOLUTION DE L ENDOMMAGEMENT---------------------------------------
  253. C
  254.  
  255. C Calcul du seuil
  256. XFD = SEUIL1 - ((XFT**2)/(2.0D0*XE) + XFFT)
  257.  
  258. TRAS1 = 1.0D0
  259.  
  260. C Evolution explicite de l'endommagement
  261. IF ((XFD.GT.0.0D0).AND.(XBEN.GE.0.0D0)) THEN
  262. XDT = 1.00D0-1.00D0/
  263. & (1.0D0+ALDIR*(SEUIL1-((XFT**2)/(2.0D0*XE))))
  264. XFFT = SEUIL1 - ((XFT**2)/(2.0D0*XE))
  265. ENDIF
  266.  
  267. C
  268. C-----SEUIL GLISSEMENT INTERNE-------------------------------------------
  269. C
  270.  
  271. C Calcul de sigma_pi
  272. DO I=1,3
  273. DO J=1,3
  274. SIGPI(I,J) = XDT*XK*(TRA-TRAPI-TRADN)*UNIT(I,J)+
  275. & 2.0D0*XG*XDT*(EPSDEV(I,J)-EPSPIDEV(I,J)-DEPSN(I,J))
  276. ENDDO
  277. ENDDO
  278.  
  279. C Calcul de la quantite (sigma_pi - X)
  280. DO I=1,3
  281. DO J=1,3
  282. SIGMOX(I,J) = SIGPI(I,J)-X(I,J)
  283. ENDDO
  284. ENDDO
  285.  
  286. C Calcul de I2(sigma_pi-X)
  287. CALL I2(SIGMOX,XFPI)
  288.  
  289. C Mise a 1 du multiplicateur plastique
  290. XLAM = 1.0D0
  291.  
  292. C CAlcul de la trace de Sigma_finale
  293. TRASIG1 = SIGMA1(1,1)+SIGMA1(2,2)+SIGMA1(3,3)
  294.  
  295. C Test sur la positivite du seuil de glissement interne
  296. IF ((XFPI.GT.0.0D0).AND.(XBEN.GE.0.D0)) THEN
  297.  
  298. C On retient la valeur initiale du seuil
  299. XFPI0 = XFPI
  300.  
  301. C On itere en bloquant le nombre d iteration max a IMAX
  302. DO NI=1,IMAX
  303.  
  304. C On realise un test de sortie
  305. IF ((ABS(XFPI/XFPI0 ).LE.(1.0D-8)).OR.
  306. & (ABS(XLAM).LE.(1.0D-20))) THEN
  307.  
  308. GOTO 400
  309.  
  310. ELSE
  311.  
  312. C Calcul des derivees dfpi/dsigpi et dphipi/dsigpi
  313. CALL DI2(SIGMOX,TERME1)
  314. CALL DI2(SIGMOX,TERME2)
  315.  
  316. C Calcul des derivees dfpi/dX et dphipi/dX
  317. DO I=1,3
  318. DO J=1,3
  319. TERME3(I,J) = 0.0D0
  320. TERME4(I,J) = 0.0D0
  321.  
  322. TERME3(I,J) = -1.0D0*TERME2(I,J)
  323. TERME4(I,J) = -1.0D0*TERME2(I,J)+A1*X(I,J)
  324. ENDDO
  325. ENDDO
  326.  
  327. C On calcul le terme C:dfpi/dsigpi
  328. TRA_SE = TERME2(1,1)+TERME2(2,2)+TERME2(3,3)
  329.  
  330. CALL DEVIAT(TERME2,TRAV)
  331.  
  332. DO I=1,3
  333. DO J=1,3
  334. TERME5(I,J) = 0.0D0
  335. TERME5(I,J) = XK*TRA_SE*UNIT(I,J)+
  336. & 2.0D0*XG*TRAV(I,J)
  337. ENDDO
  338. ENDDO
  339.  
  340. C On calcul le denominateur du multiplicateur plastique
  341. XT1 = 0.0D0
  342. DO I=1,3
  343. DO J=1,3
  344. XT1 = XT1 + TERME1(I,J)*TERME5(I,J)
  345. ENDDO
  346. ENDDO
  347.  
  348. XT2 = 0.0D0
  349. DO I=1,3
  350. DO J=1,3
  351. XT2 = XT2 + TERME3(I,J)*TERME4(I,J)
  352. ENDDO
  353. ENDDO
  354.  
  355. C Calcul du multiplicateur plastique
  356. XLAM = XFPI/(XDT*XT1 + GAMMA1*XT2)
  357.  
  358. C On actualisatise X et sigpi
  359. DO I=1,3
  360. DO J=1,3
  361. SIGPI(I,J)=SIGPI(I,J)-XDT*XLAM*TERME5(I,J)
  362. X(I,J)=X(I,J)-GAMMA1*XLAM*TERME4(I,J)
  363. ENDDO
  364. ENDDO
  365.  
  366. C On recalcule le nouveau seuil
  367. DO I=1,3
  368. DO J=1,3
  369. SIGMOX(I,J) = 0.0D0
  370. SIGMOX(I,J) = SIGPI(I,J)-X(I,J)
  371. ENDDO
  372. ENDDO
  373.  
  374. CALL I2(SIGMOX,XFPI)
  375.  
  376. ENDIF
  377.  
  378. ENDDO
  379.  
  380. 400 CONTINUE
  381.  
  382.  
  383. C On calcul epspi
  384. TRA_SE = SIGPI(1,1)+SIGPI(2,2)+SIGPI(3,3)
  385. CALL DEVIAT(SIGPI,TRAV)
  386.  
  387. TRAPI = (TRA-TRADN) - (TRA_SE/(XDT*3.0D0*XK))
  388.  
  389. DO I=1,3
  390. DO J=1,3
  391. EPSPIDEV(I,J) = (EPSDEV(I,J)-DEPSN(I,J))-
  392. & (1.0D0/(2.0D0*XG*XDT))*TRAV(I,J)
  393. ENDDO
  394. ENDDO
  395.  
  396. DO I=1,3
  397. DO J=1,3
  398. EPSPI(I,J) = (TRAPI/3.0D0)*UNIT(I,J)+EPSPIDEV(I,J)
  399. ENDDO
  400. ENDDO
  401.  
  402. ENDIF
  403. C
  404. C-----CONTRAINTES ACTUALISEE
  405. C
  406. TRA1 = SIGMA1(1,1)+SIGMA1(2,2)+SIGMA1(3,3)
  407.  
  408. IF (XBEN.GE.0.0D0) THEN
  409.  
  410. 752 CONTINUE
  411.  
  412. DO I=1,3
  413. DO J=1,3
  414. SIGMA1(I,J)=XK*(TRA-XDT*TRAPI-TRADN)*UNIT(I,J)+
  415. & 2.0D0*XG*
  416. & (EPSDEV(I,J)-XDT*EPSPIDEV(I,J)-DEPSN(I,J))
  417. ENDDO
  418. ENDDO
  419.  
  420. ELSEIF (XBEN.LE.SIGRF) THEN
  421.  
  422. C La contrainte effective
  423. DO I=1,3
  424. DO J=1,3
  425. SIGMA1(I,J)=XK*(TRA-TRADN)*UNIT(I,J)+
  426. & 2.0D0*XG*(EPSDEV(I,J)-DEPSN(I,J))
  427. ENDDO
  428. ENDDO
  429.  
  430. C On verifie le seuile de plasticite
  431. CALL J2(SIGMA1,TEMP1)
  432.  
  433. TEMP2 = SIGMA1(1,1)+SIGMA1(2,2)+SIGMA1(3,3)
  434.  
  435. TEMP2 = TEMP2/3.00
  436.  
  437. XFP = XAF*TEMP1+XBF*TEMP2 - (XFC + XECO)
  438.  
  439. C On verifie la positivite du seuil
  440. IF (XFP.GT.0.0D0) THEN
  441.  
  442. C On retient la valeur initiale pour le critere relatif
  443. XFP0 = XFP
  444.  
  445. C On commence les choses serieuses
  446. DO NP=1,IMAX
  447.  
  448. C La derivee du critere
  449. CALL DI2(SIGMA1,XDJ2P)
  450.  
  451. DO I=1,3
  452. DO J=1,3
  453. XDFP(I,J) = XAF*XDJ2P(I,J)+XBF*UNIT(I,J)/3.0D0
  454. XDGP(I,J) = XAG*XDJ2P(I,J)+XBG*UNIT(I,J)/3.0D0
  455. ENDDO
  456. ENDDO
  457.  
  458. C Le hessien elasique
  459. CALL DEVIAT(XDGP,DEVGP)
  460.  
  461. XTGP = XDGP(1,1)+XDGP(2,2)+XDGP(3,3)
  462.  
  463. DO I=1,3
  464. DO J=1,3
  465. TRAV(I,J)=0.0D0
  466. TRAV(I,J)=XK*XTGP*UNIT(I,J)+2.0D0*XG*DEVGP(I,J)
  467. ENDDO
  468. ENDDO
  469.  
  470. C La derivee de l ecrouissage
  471. XTEMP= EXP(-1.0D0*XBC*XDPP)*(XAC-(XAC*XDPP+XFC+XSIGU)*XBC)
  472.  
  473. C Le critere actuel
  474.  
  475. CALL J2(SIGMA1,TEMP1)
  476.  
  477. TEMP2 = SIGMA1(1,1)+SIGMA1(2,2)+SIGMA1(3,3)
  478.  
  479. TEMP2 = TEMP2/3.00
  480.  
  481. XFP = XAF*TEMP1+XBF*TEMP2 - (XFC + XECO)
  482.  
  483. C Le multiplicateur plastique
  484. CALL DBLECONT(XDFP,TRAV,TEMP9)
  485.  
  486. XLAMP = XFP/(TEMP9+XTEMP)
  487.  
  488. C Actualisation des variables internes
  489. XECO = XECO + XLAMP*XTEMP
  490.  
  491. XDPP = XDPP + XLAMP
  492.  
  493. DO I=1,3
  494. DO J=1,3
  495. DEFN(I,J)=DEFN(I,J) + XLAMP*XDGP(I,J)
  496. ENDDO
  497. ENDDO
  498.  
  499. DO I=1,3
  500. DO J=1,3
  501. SIGMA1(I,J)=SIGMA1(I,J) - XLAMP*TRAV(I,J)
  502. ENDDO
  503. ENDDO
  504.  
  505. C On calcul le seuil pour un eventuel arret
  506. CALL J2(SIGMA1,TEMP1)
  507. TEMP2 = SIGMA1(1,1)+SIGMA1(2,2)+SIGMA1(3,3)
  508. TEMP2 = TEMP2/3.00
  509.  
  510. XFP = XAF*TEMP1+XBF*TEMP2 - (XFC + XECO)
  511.  
  512. C Le critere d arret des iterations internes
  513. IF ((ABS(XFP/XFP0).LE.(1.0D-8)).OR.(NP.GE.IMAX)) THEN
  514.  
  515. GOTO 294
  516.  
  517. ENDIF
  518.  
  519.  
  520. ENDDO
  521.  
  522. 294 CONTINUE
  523.  
  524. ENDIF
  525.  
  526. ELSE
  527.  
  528. C Calcul du deviateur de DEFN
  529. CALL DEVIAT(DEFN,DEPSN)
  530.  
  531. C Calcul de la trace de DEFN
  532. TRADN = DEFN(1,1)+DEFN(2,2)+DEFN(3,3)
  533.  
  534. TRASIG1 = 3.0D0*XK*(TRA-XDT*TRAPI-TRADN)/
  535. & (1.0D0-(3.0D0*XK/SIGRF)*XDT*TRAPI)
  536.  
  537. DO I=1,3
  538. DO J=1,3
  539. SIG1DEV(I,J)=2.00D0*XG*(EPSDEV(I,J)-XDT*EPSPIDEV(I,J)*
  540. & (1.0D0-(TRASIG1/SIGRF))-DEPSN(I,J))
  541. ENDDO
  542. ENDDO
  543.  
  544. DO I=1,3
  545. DO J=1,3
  546. SIGMA1(I,J)=TRASIG1*UNIT(I,J)/3.0D0+SIG1DEV(I,J)
  547. ENDDO
  548. ENDDO
  549.  
  550. ENDIF
  551.  
  552. 101 CONTINUE
  553.  
  554.  
  555. IF (ABS(SIGMA1(1,1)).GE.1.0D50) THEN
  556. print*,'TRASIG1',TRASIG1
  557. print*,'QUANT=',(1.0D0-(3.0D0*XK/SIGRF)*XDT*TRAPI)
  558. print*,'HEHEF=',(TRA-XDT*TRAPI-TRADN)
  559. print*,'XEXEF=',TRAPI
  560. print*,'XEXEF=',TRADN
  561. print*,'FEFEF=',XDT
  562. print*,'TRA=',TRA
  563. print*,'DEPST',(DEPST(I),I=1,6)
  564. ENDIF
  565. C
  566. C-----CALCUL DE L ENDOMMAGEMENT EFFECTIF (GRANDEUR INDICATIVE)
  567. C-----PAS DE SENS THERMO
  568. C
  569. XFSIG = 1.0D0-(SIGMA1(1,1)+SIGMA1(2,2)+SIGMA1(3,3))/SIGRF
  570.  
  571. IF ((XFSIG.LE.1.0D0).AND.(XFSIG.GE.0.0D0)) THEN
  572. XDEFF = XDT*XFSIG
  573. ELSE
  574. XDEFF = XDT
  575. ENDIF
  576.  
  577. 444 CONTINUE
  578.  
  579. IF (ABS(SIGMA1(3,3)).LE.1.0D0) THEN
  580.  
  581. GOTO 678
  582.  
  583. ELSE
  584.  
  585. DEPST(3) = -SIGMA1(3,3)/XE
  586. IF (ABS(DEPST(3)).GE.1.0D10) THEN
  587. print*,'SIGMA1(3,3)=',SIGMA1(3,3)
  588. print*,'soucis d actualisation'
  589. print*,'XE=',XE
  590. print*,'DPEST3',DEPST(3)
  591.  
  592.  
  593. ENDIF
  594.  
  595. VAR0(18) = EPSI(3,3)
  596.  
  597. C Tenseur de glissement interne
  598. VAR0(2) = EPSPI(1,1)
  599. VAR0(3) = EPSPI(2,2)
  600. VAR0(4) = EPSPI(3,3)
  601. VAR0(5) = EPSPI(1,2)
  602. VAR0(6) = EPSPI(1,3)
  603. VAR0(7) = EPSPI(2,3)
  604.  
  605. C Tenseur de la contrainte de rappel
  606. VAR0(8) = X(1,1)
  607. VAR0(9) = X(2,2)
  608. VAR0(10) = X(3,3)
  609. VAR0(11) = X(1,2)
  610. VAR0(12) = X(1,3)
  611. VAR0(13) = X(2,3)
  612.  
  613. C Tenseur des deformations plastiques
  614. VAR0(35) = DEFN(1,1)
  615. VAR0(36) = DEFN(2,2)
  616. VAR0(37) = DEFN(3,3)
  617. VAR0(38) = DEFN(1,2)
  618. VAR0(39) = DEFN(1,3)
  619. VAR0(40) = DEFN(2,3)
  620.  
  621. C Endommagement
  622. VAR0(14) = XDT
  623.  
  624. C Force associee a l ecrouissage isotrope
  625. VAR0(15) = XFFT
  626.  
  627. C Contraintes tests
  628. VAR0(22) = SIGMA1(1,1)
  629. VAR0(23) = SIGMA1(2,2)
  630. VAR0(24) = SIGMA1(3,3)
  631. VAR0(25) = SIGMA1(1,2)
  632. VAR0(26) = SIGMA1(1,3)
  633. VAR0(27) = SIGMA1(2,3)
  634.  
  635. C Ecrouissage plastique
  636. VAR0(41) = XECO
  637.  
  638. C Deformation plastique cumulee
  639. VAR0(42) = XDPP
  640.  
  641. ENDIF
  642.  
  643. ENDDO
  644.  
  645. 678 CONTINUE
  646.  
  647. C
  648. C-----STOCKAGE DES CONTRAINTES AVEC SYMMETRISATION POUR ROBUSTESSE
  649. C-----ON MET SIG_33 BIEN QU ELLE SOIT CALCULEE A UNE TOL. PRES
  650. C
  651. SIGF(1)= SIGMA1(1,1)
  652. SIGF(2)= SIGMA1(2,2)
  653. SIGF(3)= SIGMA1(3,3)*0.0D0
  654. SIGF(4)= (SIGMA1(1,2)+SIGMA1(2,1))/2.0D0
  655.  
  656. C
  657. C-----STOCKAGE DES DEFORMATIONS DE GLISSEMENT
  658. C
  659. VARF(2)=EPSPI(1,1)
  660. VARF(3)=EPSPI(2,2)
  661. VARF(4)=EPSPI(3,3)
  662. VARF(5)=EPSPI(1,2)
  663. VARF(6)=EPSPI(1,3)
  664. VARF(7)=EPSPI(2,3)
  665. C
  666. C-----TENSEUR DES DEFORMATIONS PLASTIQUES
  667. C
  668. VARF(35) = DEFN(1,1)
  669. VARF(36) = DEFN(2,2)
  670. VARF(37) = DEFN(3,3)
  671. VARF(38) = DEFN(1,2)
  672. VARF(39) = DEFN(1,3)
  673. VARF(40) = DEFN(2,3)
  674. C
  675. C-----Ecrouissage plastique
  676. C
  677. VARF(41) = XECO
  678. C
  679. C-----Deformation plastique cumulee
  680. C
  681. VARF(42) = XDPP
  682. C
  683. C-----STOCKAGE DE L ECROUISSAGE CINEMATIQUE
  684. C
  685. VARF(8)=X(1,1)
  686. VARF(9)=X(2,2)
  687. VARF(10)=X(3,3)
  688. VARF(11)=X(1,2)
  689. VARF(12)=X(1,3)
  690. VARF(13)=X(2,3)
  691. C
  692. C-----STOCKAGE DES CONTRAINTES TESTS
  693. C
  694. VARF(22)=SIGMA1(1,1)
  695. VARF(23)=SIGMA1(2,2)
  696. VARF(24)=SIGMA1(3,3)
  697. VARF(25)=SIGMA1(1,2)
  698. VARF(26)=SIGMA1(1,3)
  699. VARF(27)=SIGMA1(2,3)
  700. C
  701. C-----STOCKAGE DES CONTRAINTES TESTS
  702. C
  703. VARF(29)=SIGPI(1,1)
  704. VARF(30)=SIGPI(2,2)
  705. VARF(31)=SIGPI(3,3)
  706. VARF(32)=SIGPI(1,2)
  707. VARF(33)=SIGPI(1,3)
  708. VARF(34)=SIGPI(2,3)
  709. C
  710. C-----ENDOMMAGEMENT EFFECTIF
  711. C
  712. VARF(28) = XDEFF
  713. C
  714. C-----STOCKAGE DE L ENDOMMAGEMENT
  715. C
  716. VARF(14) = XDT
  717. C
  718. C-----STOCKAGE DES FORCES THERMODYNAMIQUES ASSOCIEES A L ECROUISSAGE ISOTROPE
  719. C
  720. VARF(15) = XFFT
  721. C
  722. C-----STOCKAGE DES DEFORMATIONS TOTALES
  723. C
  724. VARF(16)=EPSI(1,1)
  725. VARF(17)=EPSI(2,2)
  726. VARF(18)=EPSI(3,3)
  727. VARF(19)=EPSI(1,2)
  728. VARF(20)=EPSI(1,3)
  729. VARF(21)=EPSI(2,3)
  730. C
  731. C-----BALISE DE SORTIE POUR LE NON LOCAL
  732. C
  733. 2000 CONTINUE
  734.  
  735. RETURN
  736. END
  737.  
  738.  
  739.  
  740.  
  741.  
  742.  
  743.  
  744.  
  745.  
  746.  
  747.  
  748.  
  749.  
  750.  
  751.  
  752.  
  753.  
  754.  
  755.  
  756.  
  757.  
  758.  
  759.  
  760.  
  761.  
  762.  

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