Télécharger t4m.eso

Retour à la liste

Numérotation des lignes :

t4m
  1. C T4M SOURCE JK148537 23/03/20 21:15:07 11638
  2. C t4m SOURCE G-M ROUX CHAUFFAGE+REFROIDISSEMENT V6
  3. SUBROUTINE t4m(iwrk52,iwrk53,iele,igau)
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8(A-H,O-Z)
  6. C
  7.  
  8. -INC PPARAM
  9. -INC CCOPTIO
  10. -INC CCHAMP
  11. -INC DECHE
  12. REAL*8 TINC0,TINC1
  13.  
  14.  
  15. TINC0 = 0.D0
  16. TINC1 = 0.D0
  17. wrk52 = iwrk52
  18. wrk53 = iwrk53
  19. C>>>>>>Initialisations des champs
  20. T0 = ture0(1)
  21. T1 = turef(1)
  22. Tm=(T0 + T1) / 2.
  23. TPS0=TEMP0
  24. TPS1=TEMPF
  25. cJK148537 28/02/23
  26. c PHA10 = valma0(1)
  27. c PHA20 = valma0(2)
  28. c PHA30 = valma0(3)
  29. PHA10 = rhas0(1)
  30. PHA20 = rhas0(2)
  31. PHA30 = rhas0(3)
  32. TINC0 = valma0(1)
  33. PHA40 = rhas0(4)
  34. rhasF(1)=PHA10
  35. rhasF(2)=PHA20
  36. rhasF(3)=PHA30
  37. XMATF(1)=TINC0
  38. XMATF(2)=0.
  39. rhasF(4)=PHA40
  40. XMATF(12)=valma0(12)
  41. C>>>>>>Initialisation des constantes
  42. PARR_Ar=valma0(9)
  43. PARR_Br=valma0(10)
  44. PARR_Cr=valma0(11)
  45. AC1=valma0(2)
  46. PARJMA_k=valma0(3)
  47. PARJMA_m=valma0(4)
  48. PARKM_ms=valma0(7)
  49. PARKM_ks=valma0(8)
  50. PARHB_k=valma0(5)
  51. PARHB_e=valma0(6)
  52. PARHB_n=valma0(13)
  53. TF1=valma0(14)
  54. TF2=valma0(15)
  55. C
  56. C>>>>>germination de la phase austénitique
  57. IF (Tm.GT.AC1.AND.TINC0.LT.1.) THEN
  58. t_inc1= Tm - AC1
  59. t_inc2 = (PARR_Cr / t_inc1)
  60. t_inc3 = EXP(t_inc2)
  61. t_inc4 = PARR_Br - Tm
  62. IF (t_inc4.LE.0.) THEN
  63. XMATF(1)=1.
  64. ELSE
  65. t_inc=t_inc3 * t_inc4 * PARR_Ar
  66. TINC1 = TINC0 + ((TPS1 - TPS0) / t_inc)
  67. XMATF(1)=TINC1
  68. ENDIF
  69. IF (TINC1.GT.1) THEN
  70. XMATF(1) = 1.
  71. ENDIF
  72. ENDIF
  73. C ->si la germination est terminé, tinc=1
  74. IF (Tm.GE.AC1.AND.TINC0.EQ.1.OR.PHA20.GT.0) THEN
  75. XMATF(1) = 1.
  76. ENDIF
  77. C ->si T<Ms et yg=0
  78. IF (Tm.LT.PARKM_ms.AND.PHA20.EQ.0) THEN
  79. XMATF(1) = 0.
  80. ENDIF
  81. C ->si T>Assat
  82. IF (Tm.GT.PARR_Br) THEN
  83. XMATF(1) = 1.
  84. ENDIF
  85. C
  86. C>>>>>>croissance de la phase austénitique (méthode de Heun)
  87. IF (T1.GT.AC1.AND.XMATF(1).EQ.1.AND.PHA20.LT.1) THEN
  88. IF (T1.LE.1000) THEN
  89. C ->on passes les caractéristiques thermiques en Kelvin
  90. AC1=AC1 + 273.
  91. T1b=T1 + 273.
  92. T0b=T0 + 273.
  93. C ->calcul de Teq(t=ti)
  94. Teq_11 = AC1 * PARJMA_k
  95. Teq_12 = 1. - PHA20
  96. Teq_13 = LOG(Teq_12)
  97. Teq_14 = Teq_11 + ((-1. * Teq_13) ** (1. / PARJMA_m))
  98. Teq_1 = (Teq_14 / PARJMA_k)
  99. C ->calcul de dyA/dt(t=ti)
  100. dPHA2_1_1 = (-1. * PARHB_e) / (8.31 * T0b)
  101. dPHA2_1_2 = PARHB_k * EXP(dPHA2_1_1)
  102. dPHA2_1_3=ABS(T0b - Teq_1)
  103. dPHA2_1 = dPHA2_1_2 * (dPHA2_1_3 ** PARHB_n) * (1. - PHA20)
  104. C ->calcul de yA(t=ti)
  105. PHA2_1=PHA20 + dPHA2_1 * (TPS1 - TPS0)
  106. IF (PHA2_1.GT.9.999d-01) THEN
  107. PHA21=1.
  108. ELSE
  109. C ->calcul de Teq(t=ti+1)
  110. Teq_21 = AC1 * PARJMA_k
  111. Teq_22 = 1. - PHA2_1
  112. Teq_23 = LOG(Teq_22)
  113. Teq_24 = Teq_21 + ((-1. * Teq_23) ** (1. / PARJMA_m))
  114. Teq_2 = (Teq_24 / PARJMA_K)
  115. C ->calcul de dyA/dt(t=ti+1)
  116. dPHA2_2_1 = (-1. * PARHB_e) / (8.31 * T1b)
  117. dPHA2_2_2 = PARHB_k * EXP(dPHA2_2_1)
  118. dPHA2_2_3=ABS(T1b - Teq_2)
  119. dPHA2_2 = dPHA2_2_2 * (dPHA2_2_3 ** PARHB_n) * (1. - PHA2_1)
  120. C ->calcul de yA(t)(ti+1)
  121. PHA21=PHA20 + ( (dPHA2_1 + dPHA2_2) / 2.)*(TPS1 - TPS0)
  122. IF (PHA21.GT.9.9d-01) THEN
  123. PHA21=1.
  124. ENDIF
  125. ENDIF
  126. IF (PHA30.GT.0.AND.PHA10.GT.0.) THEN
  127. C ->calcul des nouvelles phases si martensite+martensite revenue
  128. PHA11=PHA10-(PHA21-PHA20)/2.
  129. PHA31=PHA30-(PHA21-PHA20)/2.
  130. IF (PHA11.GT.0.) THEN
  131. PHA11=0.
  132. PHA31=1 - PHA21
  133. ENDIF
  134. IF (PHA31.GT.0.) THEN
  135. PHA31=0.
  136. PHA11=1 - PHA21
  137. ENDIF
  138. rhasF(2)=PHA21
  139. rhasF(1)=PHA11
  140. rhasF(3)=PHA31
  141. XMATF(12)=PHA21
  142. ELSE
  143. IF (PHA10.GT.0.) THEN
  144. C ->calcul des nouvelles phases si martensite revenue
  145. rhasF(2)=PHA21
  146. rhasF(1)=1. - PHA21
  147. rhasF(3)=0.
  148. XMATF(12)=PHA21
  149. ELSE
  150. C ->calcul des nouvelles phases si martensite
  151. rhasF(2)=PHA21
  152. rhasF(3)=1. - PHA21
  153. rhasF(1)=0.
  154. XMATF(12)=PHA21
  155. ENDIF
  156. ENDIF
  157. ENDIF
  158. ENDIF
  159. IF (T1.GT.1000.AND.T1.LT.TF1) THEN
  160. rhasF(2)=1.
  161. rhasF(3)=0.
  162. rhasF(1)=0.
  163. XMATF(12)=1.
  164. ENDIF
  165. C>>>>>transformation liquide-solide(loi linéaire)
  166. IF (T1.GE.TF1.AND.T1.LE.TF2) THEN
  167. rhasF(4)=(T1 - TF1) * (1.d0 / (TF2 - TF1))
  168. rhasF(2)=1.d0 - rhasF(4)
  169. rhasF(1)=0.
  170. rhasF(3)=0.
  171. XMATF(12)=1.
  172. ENDIF
  173. IF (T1.GE.TF2) THEN
  174. rhasF(4)=1.
  175. rhasF(2)=0.
  176. rhasF(1)=0.
  177. XMATF(12)=1.
  178. XMATF(1)=1.
  179. rhasF(3)=0.
  180. ENDIF
  181. IF (T1.LT.TF1.AND.PHA40.GT.0.) THEN
  182. rhasF(4)=0.
  183. rhasF(2)=1.
  184. rhasF(1)=0.
  185. XMATF(12)=1.
  186. rhasF(3)=0.
  187. ENDIF
  188. C
  189. C>>>>>transformation martensitique
  190. IF (T1.LT.T0.AND.PHA20.GT.0.AND.T1.LE.PARKM_ms) THEN
  191. PHA31_1 = -1. * PARKM_ks *(PARKM_ms - T1)
  192. PHA31_2 = EXP(PHA31_1)
  193. PHA31_3=XMATF(12) * (1-PHA31_2)
  194. PHA11=PHA10
  195. PHA21=XMATF(12) - PHA31_3
  196. PHA31=1 - (PHA11 + PHA21)
  197. IF (PHA21.LT.1.d-03) THEN
  198. PHA21=0.
  199. PHA11=PHA10
  200. PHA31=1 - PHA11
  201. rhasF(4)=0.
  202. ENDIF
  203. C ->sorties
  204. rhasF(1)=PHA11
  205. rhasF(2)=PHA21
  206. rhasF(3)=PHA31
  207. ENDIF
  208. C
  209. C>>>>>transformation martensite revenue
  210. rhasF(5)=0.
  211. C
  212. c if(iele.eq.1.and.igau.eq.3) write(6,*) (rhasF(jj),jj= 1,5)
  213. do 80 jj = 2,15
  214. if (jj.eq.12) goto 80
  215. xmatf(jj) = valma0(jj)
  216. 80 continue
  217. C
  218. RETURN
  219. END
  220.  
  221.  
  222.  
  223.  
  224.  
  225.  
  226.  
  227.  
  228.  
  229.  
  230.  
  231.  
  232.  
  233.  
  234.  

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