Télécharger ouglof.eso

Retour à la liste

Numérotation des lignes :

  1. C OUGLOF SOURCE BR232186 16/09/12 12:43:41 9078
  2. SUBROUTINE OUGLOF(XMAT,DEPS,SIG0,VAR0,SIGF,VARF)
  3.  
  4. C
  5.  
  6. C====&===1=========2=========3=========4=========5=========6=========7==
  7.  
  8. C Commentaires : Subroutine permettant de mettre en oeuvre le
  9.  
  10. C modele OUGLOVA pour representer
  11.  
  12. C le comportement de l'acier corrodé
  13.  
  14. C
  15.  
  16. C
  17.  
  18. C Auteurs : R. Paredes
  19.  
  20. C : CEA-DEN/DANS/DM2S/SEMT/EMSI
  21.  
  22. C : Romili.Paredes@cea.fr
  23.  
  24. C
  25.  
  26. C Date : Aout 2016
  27.  
  28. C====&===1=========2=========3=========4=========5=========6=========7==
  29.  
  30. C
  31.  
  32. C-----DECLARATION GENERALE----------------------------------------------
  33.  
  34. C
  35.  
  36. IMPLICIT INTEGER(I-N)
  37. IMPLICIT REAL*8(A-H,O-Z)
  38.  
  39. C
  40.  
  41. C-----DECLARATION DES VARIABLES-----------------------------------------
  42.  
  43. C
  44.  
  45. DIMENSION XMAT(*),DEPS(*),SIG0(*),VAR0(*),SIGF(*),VARF(*)
  46.  
  47.  
  48.  
  49. C Parametre pour le nombre d iteration locales internes
  50.  
  51. IMAX = 1000
  52.  
  53. C
  54.  
  55. C-----DONNEES MATERIAUX-------------------------------------------------
  56.  
  57. C
  58.  
  59. XYG = XMAT(1)
  60.  
  61. XNU = XMAT(2)
  62.  
  63. XGC = XYG/(2.0D0*(XNU + 1.0D0))
  64.  
  65. XSIGY = XMAT(5)
  66.  
  67. XK = XMAT(6)
  68.  
  69. XM = XMAT(7)
  70.  
  71. XTC = XMAT(8)
  72.  
  73. XDC = XMAT(9)
  74.  
  75. C
  76.  
  77. C-----SEUIL DE CORROSION------------------------------------------------
  78.  
  79. C
  80.  
  81. C-----Deformation plastique a rupture
  82.  
  83. IF (XTC.LE.15.0D0) THEN
  84.  
  85. EPSR = -0.0111D0*XTC + 0.2345D0
  86.  
  87. ELSE
  88.  
  89. EPSR = -0.0006D0*XTC + 0.051D0
  90.  
  91. ENDIF
  92.  
  93.  
  94.  
  95. C-----Seuil de rupture et d'endommagement
  96.  
  97. XPR = EPSR
  98.  
  99. XPD = 0.8D0*XPR
  100.  
  101. C
  102.  
  103. C-----ETAT INITIAL------------------------------------------------------
  104.  
  105. C
  106.  
  107. XD = VAR0(1)
  108.  
  109. XR = VAR0(2)
  110.  
  111. XP = VAR0(3)
  112.  
  113. XZT = VAR0(4)
  114.  
  115. XNRUP = VAR0(5)
  116.  
  117. XEPSP = VAR0(6)
  118.  
  119. XEPSF = DEPS(1) + VAR0(10)
  120.  
  121. C
  122.  
  123. C-----CALCUL------------------------------------------------------------
  124.  
  125. C
  126.  
  127.  
  128.  
  129. C-----Prediction elastique
  130.  
  131. SIG = (1.0D0 - XD)*XYG*(XEPSF - XEPSP)
  132.  
  133. SIGEQ = SIG
  134.  
  135.  
  136.  
  137. C-----Boucle d'endommagement
  138.  
  139. DO I = 1,IMAX
  140.  
  141. C--------Evaluation de la fonction seuil
  142.  
  143. FP = (SIGEQ/(1.0D0 - XD)) - (XR + XSIGY)
  144.  
  145. FP0 = FP
  146.  
  147.  
  148.  
  149. C--------Boucle de plasticite
  150.  
  151. IF (FP.GT.1.0D0) THEN
  152.  
  153.  
  154.  
  155. DO J = 1,IMAX
  156.  
  157. C--------------Calcul des derives
  158.  
  159. DFDSG = SIG/(SIGEQ*(1.0D0 - XD))
  160.  
  161. DPHDSG = DFDSG
  162.  
  163.  
  164.  
  165. C--------------Calcul du multiplicateur plastique
  166.  
  167. TAU = XYG*DPHDSG
  168.  
  169. AUX1 = DFDSG*TAU
  170.  
  171. AUX2 = SIGEQ/(XK*(1.0D0 - XD))
  172.  
  173. AUX3 = XSIGY/XK
  174.  
  175. AUX4 = 1.0D0 - XM
  176.  
  177. AUX5 = XK/XM
  178.  
  179. XLAMBDA = FP/(AUX1 + AUX5*((AUX2 - AUX3)**AUX4))
  180.  
  181.  
  182.  
  183. C--------------Mise à jour des variables
  184.  
  185. SIG = SIG - XLAMBDA*TAU
  186.  
  187. SIGEQ = SIG
  188.  
  189. XEPSP = XEPSP + XLAMBDA*DPHDSG
  190.  
  191. XP = XEPSP
  192.  
  193. XR = XK*(XP**(1.0D0/XM))
  194.  
  195.  
  196.  
  197. C--------------Reevaluation de la fonction seuil
  198.  
  199. FP = (SIGEQ/(1.0D0 - XD)) - (XR + XSIGY)
  200.  
  201.  
  202.  
  203. CRIT1 = DABS(FP/FP0)
  204.  
  205. CRIT2 = DABS(XLAMBDA)
  206.  
  207.  
  208.  
  209. IF ((CRIT1.LE.1.0D-8).OR.(CRIT2.LE.1.0D-10)) THEN
  210.  
  211. GOTO 910
  212.  
  213. ENDIF
  214.  
  215.  
  216.  
  217. ENDDO
  218.  
  219.  
  220.  
  221. C-----Sortie de la boucle de plasticité
  222.  
  223. 910 CONTINUE
  224.  
  225.  
  226.  
  227. C-----------Evaluation du seuil d'endommagement
  228.  
  229. FD = XP - (XPD + XZT)
  230.  
  231. IF (I.EQ.1) THEN
  232.  
  233. FD0 = FD
  234.  
  235. ENDIF
  236.  
  237.  
  238.  
  239. CRIT = DABS(FD/FD0)
  240.  
  241.  
  242.  
  243. IF ((FD.GT.0.0D0).AND.(CRIT.GT.1.0D-5)) THEN
  244.  
  245.  
  246.  
  247. C--------------Calcul de l'endommagement
  248.  
  249. XD = (XDC/(XPR - XPD))*(XP - XPD)
  250.  
  251. XZT = XP - XPD
  252.  
  253.  
  254.  
  255. C--------------Endommagement critique
  256.  
  257. IF (XD.GT.XDC) THEN
  258. XD = XDC
  259.  
  260. ENDIF
  261.  
  262.  
  263.  
  264. ELSE
  265.  
  266. GOTO 900
  267.  
  268. ENDIF
  269.  
  270.  
  271.  
  272. IF (CRIT.LE.1.0D-5) THEN
  273.  
  274. GOTO 900
  275.  
  276. ENDIF
  277.  
  278.  
  279.  
  280. ELSE
  281.  
  282. GOTO 900
  283.  
  284. ENDIF
  285.  
  286.  
  287.  
  288. ENDDO
  289.  
  290.  
  291.  
  292. C-----Sortie de la boucle d'endommagement
  293.  
  294. 900 CONTINUE
  295.  
  296.  
  297.  
  298. C-----Rupture
  299.  
  300. IF (XP.GE.XPR) THEN
  301.  
  302. XNRUP = 1
  303.  
  304. XP = XPR
  305.  
  306. SIG = 1.0D0
  307.  
  308. ENDIF
  309.  
  310. C
  311.  
  312. C-----SOCKAGE EN SORTIE-------------------------------------------------
  313.  
  314. C
  315.  
  316. C-----Les variables internes
  317.  
  318. VARF(1) = XD
  319.  
  320. VARF(2) = XR
  321.  
  322. VARF(3) = XP
  323.  
  324. VARF(4) = XZT
  325.  
  326. VARF(5) = XNRUP
  327.  
  328. VARF(6) = XEPSP
  329.  
  330. VARF(10) = XEPSF
  331.  
  332.  
  333.  
  334. C-----Les contraintes
  335.  
  336. SIGF(1) = SIG
  337.  
  338. SIGF(2) = SIG0(2) + XGC*DEPS(2)
  339.  
  340. SIGF(3) = SIG0(3) + XGC*DEPS(3)
  341.  
  342.  
  343.  
  344. C-----Fin de l integration
  345.  
  346. RETURN
  347.  
  348. END
  349.  

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