Télécharger simtrc.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : simtrc.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. SAUT PAGE ;
  5. complet = faux ;
  6. GRAPH = 'N' ;
  7. *********************************************************************
  8. * simtrc.dgibi *
  9. * Simulation du trc avec l'operateur comp *
  10. * Prise en compte de la concentration en carbone (cte) *
  11. * version initiale : Martinez le 30/07/98 *
  12. *********************************************************************
  13. * mettre complet = vrai pour generer le trc du 16MND5
  14. **** Parametres du calcul :
  15. **** Vitesse de refroidissement
  16. *
  17. lis_ref1 = prog -0.05 -0.1 -0.15 -0.2 -0.3 -0.5 -1. -2. -3. -4. -5.
  18. -6. -7. -8. -9. -10. -12. -15. -18. -20. -25. -30.
  19. -40. -50. -60. -80. -100. -400. ;
  20. lis_ref1 = prog -0.05 -0.1 -0.3 -0.4 -0.5 -1. -2. -3. -4. -5. -7. -10.
  21. -15. -25. -35. -40. -60. -80. -150. ;
  22. lis_ref2 = prog -0.07 -0.11 -0.18 -0.23 -0.4 -0.57 -1.3 -2.8 -3.1
  23. -4.2 -5.5 -6.4 -7.3 -8.7 -9.9 -10.01 -12.3 -15.5
  24. -19. -23. -26. -34. -45. -57. -61. -88. -300. -390. ;
  25. lis_ref = lis_ref1 ;
  26. *
  27. **** Delta Temperature
  28. d_tempe = 2. ;
  29. *
  30. ******** INITIALISATION DU DEPOUILLEMENT ***************************
  31. *
  32. dfer = 100. ;
  33. dbai = 100. ;
  34. dmar = 100. ;
  35. ffer = 100. ;
  36. fbai = 100. ;
  37. fmar = 100. ;
  38. tdf = 1. ; tff = 1. ;
  39. tdb = 1. ; tfb = 1. ;
  40. tdm = 1. ; tfm = 1. ;
  41. lis_dfer = prog ;
  42. lis_dbai = prog ;
  43. lis_dmar = prog ;
  44. lis_tdf = prog ;
  45. lis_tdb = prog ;
  46. lis_tdm = prog ;
  47. lis_ffer = prog ;
  48. lis_fbai = prog ;
  49. lis_tff = prog ;
  50. lis_tfb = prog ;
  51. lis_refr = prog ;
  52. lis_pfer = prog ;
  53. lis_pbai = prog ;
  54. lis_pmar = prog ;
  55. *
  56. **** LECTURE DES DONNEES DU DIAGRAMME TRC **************************
  57. *
  58. DTRC_1 = nuage 'TP'*flottant 'ZFF'*flottant 'ZFB'*flottant
  59. 'TDF'*flottant 'TFF'*flottant 'TDB'*flottant 'TFB'*flottant
  60. -0.05 0.67 0.33 769. 644. 548. 395.
  61. -0.1 0.45 0.55 752. 644.0 566.0 395.0
  62. -0.3 0.16 0.84 714.0 639.0 588.0 395.0
  63. -0.4 0.06 0.94 700.0 659.0 591.0 395.0
  64. -0.5 0. 1. 700.0 690.0 593.0 395.0
  65. -1.0 0. 0.96 700.0 690. 600.0 395.0
  66. -2.0 0. 0.91 700.0 690.0 588.0 395.0 ;
  67. DTRC_1 = DTRC_1 et (nuage 'TP'*flottant 'ZFF'*flottant 'ZFB'*flottant
  68. 'TDF'*flottant 'TFF'*flottant 'TDB'*flottant 'TFB'*flottant
  69. -3.0 0. 0.87 700.0 390.0 581.0 396.0
  70. -4.0 0.0 0.69 700.0 690.0 573.0 396.0
  71. -5.0 0.0 0.55 700.0 690.0 568.0 398.0
  72. -7.0 0.0 0.46 700.0 690.0 559.0 400.0
  73. -10.0 0.0 0.34 700.0 690.0 543.0 402.0
  74. -15.0 0.0 0.22 700.0 690.0 525.0 405.0
  75. -25.0 0.0 0.05 700.0 690.0 498.0 415.0) ;
  76. DTRC_1 = DTRC_1 et (nuage 'TP'*flottant 'ZFF'*flottant 'ZFB'*flottant
  77. 'TDF'*flottant 'TFF'*flottant 'TDB'*flottant 'TFB'*flottant
  78. -35.0 0.0 0. 710.0 700.0 430.0 425.0
  79. -40.0 0.0 0. 710.0 700.0 430.0 425.0
  80. -60.0 0.0 0. 710.0 700.0 430.0 425.0
  81. -80.0 0.0 0.0 710.0 700.0 430.0 425.0
  82. -150.0 0.0 0.0 710.0 700.0 430.0 425.0 );
  83. *
  84. * Chauffage (modele de leblond)
  85. *
  86. DTCH_1 = nuage
  87. 'COMP' 'TE' 730.0 750.0 770.0 790.0 810.0 830.0 840.0 860.0
  88. 880.0 900.0
  89. 'COMP' 'CK' 0.0 0.22 0.53 1.05 2.02 4.55 5.6 7.37 10.77 20.0
  90. 'COMP' 'CL' 1.0 0.97 0.94 0.87 0.76 0.45 0.0 0.0 0.0 0.0 ;
  91.  
  92. **** LECTURE DU MAILLAGE *******************************************
  93. *
  94. opti dime 2 elem qua8 mode axis ;
  95. elmat = table ;
  96. piece = table ;
  97. a00 = 0. 0. ;
  98. a01 = 0.01 0. ;
  99. a02 = 0.01 0.01 ;
  100. a03 = 0. 0.01 ;
  101. d1 = droi 1 a00 a01 ;
  102. d2 = droi 1 a01 a02 ;
  103. d3 = droi 1 a02 a03 ;
  104. d4 = droi 1 a03 a00 ;
  105. l_ext = (d1 et d2 et d3) ;
  106. tot = dalle d1 d2 d3 d4 plan ;
  107. *
  108. ******** DEFINITION D'UN MODELE *************************************
  109. *
  110. p_mode1 = mode tot 'THERMIQUE' 'ISOTROPE' 'CONS' me1 ;
  111. si faux ;
  112. p_mode2 = mode tot 'TAILGRAI' 'ISOTROPE' 'CONS' me2 ;
  113. p_mode3 = mode tot 'DIFFUSION' 'CONS' carb1 ;
  114. finsi ;
  115.  
  116. momic_1 = mode tot melange cerem ;
  117. mamic_1 = mate momic_1
  118. * changement de type de certaines donnees
  119. 'AC1' 730. 'AR1' 760. 'MS0' 422.
  120. 'BETA' -0.0307 'AC' 9146. 'AA' 37.5 'ZS' 0. 'TPLM' -0.5 'CARB' 0.0074
  121. 'ACAR' 540. 'DG0' 0.00001 'AGRA' 11200.
  122. 'TIHT' 1000. 'TFHT' 200. 'DTHT' 2. 'NHTR' DTRC_1 'NLEB' DTCH_1
  123. 'AUST' 1. 'FERR' 0. 'BAIN' 0. 'MART' 0. 'MS' 400. ;
  124. *********************************************************************
  125. ******** BOUCLE SUR LES HISTOIRES THERMIQUES ************************
  126. *********************************************************************
  127. *
  128. n_bou1 = dime lis_ref ;
  129. repeter bou_1 n_bou1 ;
  130. *
  131. ******** CHARGEMENT THERMIQUE **************************************
  132. *
  133. t0tempe = 1000.0 ;
  134. t1tempe = 100.0 ;
  135. *
  136. t_chtem = manu CHPO tot 1 'T' 1. ;
  137. t0chtem = t0tempe * t_chtem ;
  138. *
  139. app_fer = 0 ;
  140. app_bai = 0 ;
  141. app_mar = 0 ;
  142. *
  143. ******** INITIALISATION *********************************
  144. *
  145. si faux ;
  146. chcarb = manu 'CHML' p_mode3 'ACC' 0.0074 'TYPE' 'ACTIVITE' ;
  147. chdg = manu 'CHML' p_mode2 'DG' 10.e-6 'TYPE' 'GRAIN' ;
  148. finsi ;
  149. p_table = table ;
  150. p_table.'TEMPERATURES' = table ;
  151. p_table.'PHASES' = table ;
  152. p_table.'TRC' = table ;
  153. *
  154. p_table.TEMPERATURES.0 = t0chtem ;
  155. p_table.PHASES.0 = mamic_1 ;
  156. *
  157. lis_tps=prog ;
  158. lis_te =prog ;
  159. lis_pa =prog ;
  160. lis_pf =prog ;
  161. lis_pb =prog ;
  162. lis_pm =prog ;
  163. *
  164. *
  165. *********************************************************************
  166. ******** BOUCLE SUR LE TEMPS ****************************************
  167. *********************************************************************
  168. *
  169. t_point = extr lis_ref &bou_1 ;
  170. *
  171. t0temps = 0. ;
  172. t1temps = (t1tempe - t0tempe) / t_point ;
  173. atpnt = ABS t_point ;
  174. deltt = d_tempe / atpnt ;
  175. n_bou2 = enti((t1temps - t0temps) / deltt) + 1 ;
  176. *
  177. t_ps = t0temps - deltt ;
  178. *
  179. mess '=============================================================' ;
  180. mess '=============================================================' ;
  181. mess ' Vitesse de refroidissement : ' t_point ;
  182. *
  183. repeter bou_2 n_bou2 ;
  184. *
  185. chdtps1 = manu chml momic_1 'DTPS' deltt ;
  186. chtemp0 = manu chml tot temp t_ps ;
  187. chtemp1 = manu chml tot temp (t_ps + deltt) ;
  188. t_ps = t_ps + deltt ;
  189. p_table.temperatures.&bou_2 = p_table.temperatures.(&bou_2-1) +
  190. (t_chtem * t_point * deltt) ;
  191. n1 = (dime p_table.temperatures) - 1 ;
  192. si (n1 > 0) ;
  193. cot1 = p_table.temperatures.&bou_2 ;
  194. cot0 = p_table.temperatures.(&bou_2-1) ;
  195. dt1 = cot1 - cot0 ;
  196. dtp1 = dt1 / deltt ;
  197. max_dt = maxi dt1 ;
  198. cht1 = chan 'CHAM' cot1 tot ;
  199. cht0 = chan 'CHAM' cot0 tot ;
  200. chp0 = p_table.'PHASES'.(&bou_2-1) ;
  201. chtp1 = chan cham dtp1 tot ;
  202. htab1 = p_table.'TRC' ;
  203.  
  204. cho1 = comp momic_1 (cht0 et chtemp0 et chp0 et chdtps1)
  205. (cht1 et chtemp1) ;
  206.  
  207. lma_1 = extr mamic_1 comp ;
  208. chp1 = exco cho1 lma_1 ;
  209.  
  210. p_table.'PHASES'.&bou_2 = chp1 ;
  211. finsi ;
  212. *=====================================================================
  213. **** verification des temperatures de debut et de fin de chgt de phase
  214. *=====================================================================
  215. *
  216. p_fe0 = maxi (exco 'FERR' chp0) ;
  217. p_ba0 = maxi (exco 'BAIN' chp0) ;
  218. p_au = maxi (exco 'AUST' chp1) ;
  219. p_fe = maxi (exco 'FERR' chp1) ;
  220. p_ba = maxi (exco 'BAIN' chp1) ;
  221. p_ma = maxi (exco 'MART' chp1) ;
  222. t_app = maxi (exco T cot1) ;
  223. tr_ma = maxi (exco chp1 MS) ;
  224. *
  225. lis_tps= lis_tps et (prog t_ps) ;
  226. lis_te = lis_te et (prog t_app) ;
  227. lis_pa = lis_pa et (prog p_au) ;
  228. lis_pf = lis_pf et (prog p_fe) ;
  229. lis_pb = lis_pb et (prog p_ba) ;
  230. lis_pm = lis_pm et (prog p_ma) ;
  231. *
  232. **** Temperatures d'apparition :ferrite, bainite et martensite
  233. *
  234. si (p_fe > 0.001) ;
  235. si (app_fer EGA 0) ;
  236. app_fer = 1 ;
  237. t_app = maxi (exco T cot1) ;
  238. dfer = t_app ;
  239. tdf = t_ps ;
  240. finsi ;
  241. finsi ;
  242. si (p_ba > 0.001) ;
  243. si (app_bai EGA 0) ;
  244. app_bai = 1 ;
  245. t_app = maxi (exco T cot1) ;
  246. dbai = t_app ;
  247. tdb = t_ps ;
  248. finsi ;
  249. finsi ;
  250. si (p_ma > 0.001) ;
  251. si (app_mar EGA 0) ;
  252. app_mar = 1 ;
  253. t_app = maxi (exco T cot1) ;
  254. dmar = t_app ;
  255. tdm = t_ps ;
  256. finsi ;
  257. finsi ;
  258. *
  259. **** Temperatures de fin de transformation de ferrite et bainite
  260. *
  261. si (app_fer EGA 1) ;
  262. d_fe = p_fe - p_fe0 ;
  263. si (d_fe < 0.000001) ;
  264. app_fer = 2 ;
  265. t_app = maxi (exco T cot1) ;
  266. ffer = t_app ;
  267. tff = t_ps ;
  268. finsi ;
  269. finsi ;
  270. si (app_bai EGA 1) ;
  271. d_ba = p_ba - p_ba0 ;
  272. si (d_ba < 0.000001) ;
  273. app_bai = 2 ;
  274. t_app = maxi (exco T cot1) ;
  275. fbai = t_app ;
  276. tfb = t_ps ;
  277. finsi ;
  278. finsi ;
  279. *
  280. si ((non complet) et (&bou_2 > 150)) ;
  281. quitter bou_2 ;
  282. finsi ;
  283. fin bou_2 ;
  284. *
  285. lis_dfer = lis_dfer et (prog dfer) ;
  286. lis_dbai = lis_dbai et (prog dbai) ;
  287. lis_dmar = lis_dmar et (prog dmar) ;
  288. lis_tdf = lis_tdf et (prog tdf) ;
  289. lis_tdb = lis_tdb et (prog tdb) ;
  290. lis_tdm = lis_tdm et (prog tdm) ;
  291. *
  292. lis_ffer = lis_ffer et (prog ffer) ;
  293. lis_fbai = lis_fbai et (prog fbai) ;
  294. lis_tff = lis_tff et (prog tff) ;
  295. lis_tfb = lis_tfb et (prog tfb) ;
  296. *
  297. **** presentation des resultats ***************************************
  298. *
  299. mess 'ferrite' dfer ffer ;
  300. mess 'bainite' dbai fbai ;
  301. mess '==========================================================' ;
  302. *
  303. *
  304. n_f = (dime p_table.PHASES) - 1 ;
  305. ch4 = p_table.PHASES.n_f ;
  306. **** austenite
  307. p_au = maxi (exco 'AUST' ch4) ;
  308. **** ferrite
  309. p_fe = maxi (exco 'FERR' ch4) ;
  310. **** bainite
  311. p_ba = maxi (exco 'BAIN' ch4) ;
  312. **** martensite
  313. p_ma = maxi (exco 'MART' ch4) ;
  314. **** Verification des proportions *****************************
  315. p_tot = p_au + p_fe + p_ba + p_ma ;
  316. *
  317. mess ' INTG des fractions volumiques : ' p_tot ;
  318. mess '===============================================================';
  319. mess 'Tpoint %austenite %ferrite %bainite %martensite';
  320. mess t_point p_au p_fe p_ba p_ma ;
  321. *
  322. si (non complet) ;
  323. quitter bou_1 ;
  324. finsi ;
  325. fin bou_1 ;
  326. *
  327. ************************************************************************
  328. ******** TRACE DU TRC **************************************************
  329. ************************************************************************
  330. *
  331. ev_dfer = evol bleu manu TPS lis_tdf T lis_dfer ;
  332. ev_dbai = evol vert manu TPS lis_tdb T lis_dbai ;
  333. ev_dmar = evol manu TPS lis_tdm T lis_dmar ;
  334. ev_ffer = evol bleu manu TPS lis_tff T lis_ffer ;
  335. ev_fbai = evol vert manu TPS lis_tfb T lis_fbai ;
  336. *
  337. tabdess = table;
  338. tabdess . 1 = 'MARQ TRIA NOLI' ;
  339. tabdess . 2 = 'MARQ LOSA NOLI' ;
  340. tabdess . 3 = 'MARQ CARR NOLI' ;
  341. tabdess . 4 = 'MARQ TRIA NOLI' ;
  342. tabdess . 5 = 'MARQ LOSA NOLI' ;
  343. tabdess . 'TITRE' = table ;
  344. tabdess . 'TITRE' . 1 = MOT 'ferrite' ;
  345. tabdess . 'TITRE' . 2 = MOT 'bainite' ;
  346. tabdess . 'TITRE' . 3 = MOT 'martensite' ;
  347. tabdess . 'TITRE' . 4 = MOT 'ferrite' ;
  348. tabdess . 'TITRE' . 5 = MOT 'bainite' ;
  349. *
  350. 'SI' ('EGA' GRAPH 'O') ;
  351. opti trac psc ;
  352. dess (ev_dfer et ev_dbai et ev_dmar et ev_ffer et ev_fbai) 'LOGX'
  353. LEGE TITX 'temps' TITY 'temperature' tabdess ;
  354. 'FINSI' ;
  355. *
  356. *
  357. si (non complet) ;
  358. p_au0 = .65196 ;
  359.  
  360. RESI=ABS((p_au-p_au0)/p_au0);
  361. SI(RESI &lt;EG 1E-2);
  362. ERRE 0;
  363. SINO;
  364. ERRE 5;
  365. FINSI;
  366. finsi ;
  367.  
  368. si complet ;
  369. opti sauv 'simtrc.k2000' ;
  370. finsi ;
  371. fin;
  372.  
  373.  
  374.  

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