Télécharger biblio.procedur

Retour à la liste

Numérotation des lignes :

  1. * BIBLIO PROCEDUR PASCAL 21/03/09 21:15:00 10917
  2. DEBP BIBLIO mot1*'MOT' mot2/'MOT' ;
  3.  
  4. *------------------------- Lecture arguments -------------------------*
  5.  
  6. * Lecture indice reference biblio associee a mot1 :
  7. si (exis mot2) ;
  8. si (ega mot2 'REFE') ;
  9. argu ent1*'ENTIER' ;
  10. sino ;
  11. erre '***** ERREUR : on attend le mot-cle REFE' ;
  12. quit biblio ;
  13. fins ;
  14. sino ;
  15. ent1 = 1 ;
  16. fins ;
  17.  
  18. * imot1 : indice du mot dans la base biblio, 0 si pas dans la base.
  19. imot1 = 0 ;
  20.  
  21. *----------------------------------------------------------------------*
  22. * MOT-CLE : 316L *
  23. *----------------------------------------------------------------------*
  24.  
  25. si (ega mot1 '316L') ;
  26. imot1 = 1 ;
  27.  
  28. * Initialisation table de sortie :
  29. TBIB1 = table ;
  30.  
  31. *------------------------- 316L, Reference 1 --------------------------*
  32.  
  33. * Proprietes thermomecaniques du 316L :
  34. * Lionel Depradeux.
  35. * "Simulation numerique du soudage - Acier 316L : validation sur cas
  36. * tests de complexite croissante".
  37. * These de doctorat.
  38. * Institut National des Sciences Appliquees de Lyon.
  39. * 2004.
  40.  
  41. si (ent1 ega 1) ;
  42.  
  43.  
  44. *---- THERMIQUE ---- :
  45.  
  46. * K1 : Conductivite thermique (W/m/K) :
  47. lt1 = prog 20. 100. pas 100. 1000. 1200. 1500. ;
  48. lk1 = prog 14. 15.2 16.6 17.9 19. 20.6 21.8 23.1 24.3 26. 27.3 29.9 60. ;
  49. K1 = evol manu 'T' lt1 'K' lk1 ;
  50. TBIB1 . 'K' = K1 ;
  51.  
  52. * Rho1 : Masse volumique (kg/m^3)
  53. lt1 = prog 20. 100. pas 100. 1000. 1200. 1500. ;
  54. lrho1 = prog 8000. 7970. 7940. 7890. 7850. 7800. 7750. 7700. 7660. 7610. 7570. 7450. 7200. ;
  55. Rho1 = evol manu 'T' lt1 'RHO' lrho1 ;
  56. TBIB1 . 'RHO' = Rho1 ;
  57.  
  58. * Cp1 : Capacite calorifique massique (J/K/kg)
  59. lt1 = prog 20. 100. pas 100. 1000. 1200. ;
  60. lcp1 = prog 450. 490. 525. 545. 560. 570. 580. 595. 625. 650. 660. 677. ;
  61. Cp1 = evol manu 'T' lt1 'C' lcp1 ;
  62. TBIB1 . 'C' = Cp1 ;
  63.  
  64. * Tfus1 : Temperaure de fusion (degC) :
  65. Tfus1 = 1450. ;
  66. TBIB1 . 'TFUS' = Tfus1 ;
  67.  
  68. * Qlat1 : Chaleur latente massique de fusion (J/kg)
  69. Qlat1 = 2.5e5 ;
  70. TBIB1 . 'QLAT' = Qlat1 ;
  71.  
  72. *---- MECANIQUE ---- :
  73.  
  74. * Nu1 : coefficient de Poisson
  75. TBIB1 . 'NU' = 0.3 ;
  76.  
  77. * Ym1 : module de Young (Pa) :
  78. lt1 = prog 20. 100. pas 100. 1200. ;
  79. lym1 = prog 197000. 191500. 184000. 176500. 168000. 160000. 151500. 142500. 130000. 108000. 81500. 32000. 7400. ;
  80. lym1 = 1.e6 * lym1 ;
  81. ym1 = evol manu T lt1 youn lym1 ;
  82. TBIB1 . 'YOUN' = ym1 ;
  83.  
  84. * Alph1 : dilatation thermique (/K)
  85. lt1 = prog 20. 100. 300. 400. 600. 700. 800. 1000. 1400. ;
  86. lalph1 = prog 15.5 16.0 17.1 17.5 18.4 18.7 19.0 19.4 19.6 ;
  87. lalph1 = 1.e-6 * lalph1 ;
  88. alph1 = evol manu T lt1 ALPH lalph1 ;
  89. TBIB1 . 'ALPH' = alph1 ;
  90.  
  91. * Sy1 : limite d'elasticite (Pa)
  92. lt1 = prog 20. 100. pas 100. 1400. ;
  93. lsy1 = prog 287. 237. 198. 172. 157. 151. 145. 136. 127. 115. 79. 38. 24. 18. 2. ;
  94. lsy1 = 1.e6 * lsy1 ;
  95. sy1 = evol manu T lt1 SIGY lsy1 ;
  96. TBIB1 . 'SIGY' = sy1 ;
  97.  
  98. * Trac1 : Courbes de traction (Pa)
  99. * Ecro1 : Courbes d'ecrouissage (Pa)
  100. Trac1 = vide nuage ;
  101. Ecro1 = vide nuage ;
  102. * A 20 degC :
  103. T1 = 20. ;
  104. lept1 = prog 0.3459 1.1599 5.2129 20.2847 ;
  105. lept1 = 1.e-2 * lept1 ;
  106. lsig1 = prog 287.33 314.94 419.42 560.97 ;
  107. lsig1 = 1.e6 * lsig1 ;
  108. trac20 = (evol manu 'EPST' lept1 'SIG' lsig1) coul 1 ;
  109. Trac1 = Trac1 et (nuag comp 'T' T1 comp 'TRAC' trac20) ;
  110. ymT1 = ipol ym1 T1 ;
  111. lepp2 = lept1 - (lsig1 / ymT1) ;
  112. Ecro20 = (evol manu 'EPSP' lepp2 'SIG' lsig1) coul 1 ;
  113. Ecro20 = prol line Ecro20 0. ;
  114. Ecro1 = Ecro1 et (nuag comp 'T' T1 comp 'ECRO' Ecro20) ;
  115. * A 200 degC :
  116. T1 = 200. ;
  117. lept1 = prog 0.3077 1.1257 5.1816 20.2556 ;
  118. lept1 = 1.e-2 * lept1 ;
  119. lsig1 = prog 198.16 231.36 334.1 470.24 ;
  120. lsig1 = 1.e6 * lsig1 ;
  121. trac200 = (evol manu 'EPST' lept1 'SIG' lsig1) coul 2 ;
  122. Trac1 = Trac1 et (nuag comp 'T' T1 comp 'TRAC' trac200) ;
  123. ymT1 = ipol ym1 T1 ;
  124. lepp2 = lept1 - (lsig1 / ymT1) ;
  125. Ecro200 = (evol manu 'EPSP' lepp2 'SIG' lsig1) coul 2 ;
  126. Ecro200 = prol line Ecro200 0. ;
  127. Ecro1 = Ecro1 et (nuag comp 'T' T1 comp 'ECRO' Ecro200) ;
  128. * A 300 degC :
  129. T1 = 300. ;
  130. lept1 = prog 0.2977 1.1137 5.1740 20.2558 ;
  131. lept1 = 1.e-2 * lept1 ;
  132. lsig1 = prog 172.47 200.62 307.12 451.42 ;
  133. lsig1 = 1.e6 * lsig1 ;
  134. trac300 = (evol manu 'EPST' lept1 'SIG' lsig1) coul 3 ;
  135. Trac1 = Trac1 et (nuag comp 'T' T1 comp 'TRAC' trac300) ;
  136. ymT1 = ipol ym1 T1 ;
  137. lepp2 = lept1 - (lsig1 / ymT1) ;
  138. Ecro300 = (evol manu 'EPSP' lepp2 'SIG' lsig1) coul 3 ;
  139. Ecro300 = prol line Ecro300 0. ;
  140. Ecro1 = Ecro1 et (nuag comp 'T' T1 comp 'ECRO' Ecro300) ;
  141. * A 400 degC :
  142. T1 = 400. ;
  143. lept1 = prog 0.2939 1.1115 5.1755 20.2593 ;
  144. lept1 = 1.e-2 * lept1 ;
  145. lsig1 = prog 157.83 187.32 294.89 435.66 ;
  146. lsig1 = 1.e6 * lsig1 ;
  147. trac400 = (evol manu 'EPST' lept1 'SIG' lsig1) coul 4 ;
  148. Trac1 = Trac1 et (nuag comp 'T' T1 comp 'TRAC' trac400) ;
  149. ymT1 = ipol ym1 T1 ;
  150. lepp2 = lept1 - (lsig1 / ymT1) ;
  151. Ecro400 = (evol manu 'EPSP' lepp2 'SIG' lsig1) coul 4 ;
  152. Ecro400 = prol line Ecro400 0. ;
  153. Ecro1 = Ecro1 et (nuag comp 'T' T1 comp 'ECRO' Ecro400) ;
  154. * A 500 degC :
  155. T1 = 500. ;
  156. lept1 = prog 0.2948 1.1120 5.1765 20.2620 ;
  157. lept1 = 1.e-2 * lept1 ;
  158. lsig1 = prog 151.65 179.12 282.34 419.14 ;
  159. lsig1 = 1.e6 * lsig1 ;
  160. trac500 = (evol manu 'EPST' lept1 'SIG' lsig1) coul 5 ;
  161. Trac1 = Trac1 et (nuag comp 'T' T1 comp 'TRAC' trac500) ;
  162. ymT1 = ipol ym1 T1 ;
  163. lepp2 = lept1 - (lsig1 / ymT1) ;
  164. Ecro500 = (evol manu 'EPSP' lepp2 'SIG' lsig1) coul 5 ;
  165. Ecro500 = prol line Ecro500 0. ;
  166. Ecro1 = Ecro1 et (nuag comp 'T' T1 comp 'ECRO' Ecro500) ;
  167. * A 600 degC :
  168. T1 = 600. ;
  169. lept1 = prog 0.2960 1.1121 5.1730 20.2556 ;
  170. lept1 = 1.e-2 * lept1 ;
  171. lsig1 = prog 145.47 169.87 262.17 387.23 ;
  172. lsig1 = 1.e6 * lsig1 ;
  173. trac600 = (evol manu 'EPST' lept1 'SIG' lsig1) coul 6 ;
  174. Trac1 = Trac1 et (nuag comp 'T' T1 comp 'TRAC' trac600) ;
  175. ymT1 = ipol ym1 T1 ;
  176. lepp2 = lept1 - (lsig1 / ymT1) ;
  177. Ecro600 = (evol manu 'EPSP' lepp2 'SIG' lsig1) coul 6 ;
  178. Ecro600 = prol line Ecro600 0. ;
  179. Ecro1 = Ecro1 et (nuag comp 'T' T1 comp 'ECRO' Ecro600) ;
  180. * A 700 degC :
  181. T1 = 700. ;
  182. lept1 = prog 0.2956 1.1128 5.1485 20.2389 ;
  183. lept1 = 1.e-2 * lept1 ;
  184. lsig1 = prog 136.2 160.7 211.65 340.49 ;
  185. lsig1 = 1.e6 * lsig1 ;
  186. trac700 = (evol manu 'EPST' lept1 'SIG' lsig1) coul 7 ;
  187. Trac1 = Trac1 et (nuag comp 'T' T1 comp 'TRAC' trac700) ;
  188. ymT1 = ipol ym1 T1 ;
  189. lepp2 = lept1 - (lsig1 / ymT1) ;
  190. Ecro700 = (evol manu 'EPSP' lepp2 'SIG' lsig1) coul 7 ;
  191. Ecro700 = prol line Ecro700 0. ;
  192. Ecro1 = Ecro1 et (nuag comp 'T' T1 comp 'ECRO' Ecro700) ;
  193. * A 800 degC :
  194. T1 = 800. ;
  195. lept1 = prog 0.2980 1.1167 5.1529 20.2260 ;
  196. lept1 = 1.e-2 * lept1 ;
  197. lsig1 = prog 127. 151.75 198.75 293.75 ;
  198. lsig1 = 1.e6 * lsig1 ;
  199. trac800 = (evol manu 'EPST' lept1 'SIG' lsig1) coul 8 ;
  200. Trac1 = Trac1 et (nuag comp 'T' T1 comp 'TRAC' trac800) ;
  201. ymT1 = ipol ym1 T1 ;
  202. lepp2 = lept1 - (lsig1 / ymT1) ;
  203. Ecro800 = (evol manu 'EPSP' lepp2 'SIG' lsig1) coul 8 ;
  204. Ecro800 = prol line Ecro800 0. ;
  205. Ecro1 = Ecro1 et (nuag comp 'T' T1 comp 'ECRO' Ecro800) ;
  206. * A 900 degC :
  207. T1 = 900. ;
  208. lept1 = prog 0.3068 1.1273 5.1493 20.1566 ;
  209. lept1 = 1.e-2 * lept1 ;
  210. lsig1 = prog 115.37 137.5 161.25 169.12 ;
  211. lsig1 = 1.e6 * lsig1 ;
  212. trac900 = (evol manu 'EPST' lept1 'SIG' lsig1) coul 11 ;
  213. Trac1 = Trac1 et (nuag comp 'T' T1 comp 'TRAC' trac900) ;
  214. ymT1 = ipol ym1 T1 ;
  215. lepp2 = lept1 - (lsig1 / ymT1) ;
  216. Ecro900 = (evol manu 'EPSP' lepp2 'SIG' lsig1) coul 11 ;
  217. Ecro900 = prol line Ecro900 0. ;
  218. Ecro1 = Ecro1 et (nuag comp 'T' T1 comp 'ECRO' Ecro900) ;
  219. * A 1000 degC :
  220. T1 = 1000. ;
  221. lept1 = prog 0.2966 1.1009 5.1184 20.1227 ;
  222. lept1 = 1.e-2 * lept1 ;
  223. lsig1 = prog 78.75 82.25 96.5 100. ;
  224. lsig1 = 1.e6 * lsig1 ;
  225. trac1000 = (evol manu 'EPST' lept1 'SIG' lsig1) coul 10 ;
  226. Trac1 = Trac1 et (nuag comp 'T' T1 comp 'TRAC' trac1000) ;
  227. ymT1 = ipol ym1 T1 ;
  228. lepp2 = lept1 - (lsig1 / ymT1) ;
  229. Ecro1000 = (evol manu 'EPSP' lepp2 'SIG' lsig1) coul 10 ;
  230. Ecro1000 = prol line Ecro1000 0. ;
  231. Ecro1 = Ecro1 et (nuag comp 'T' T1 comp 'ECRO' Ecro1000) ;
  232. TBIB1 . 'TRAC' = Trac1 ;
  233. TBIB1 . 'ECRO' = Ecro1 ;
  234.  
  235. *trac0 = trac20 et trac200 et trac300 et trac400 et trac500 et trac600 et trac700 et trac800 et trac900 et trac1000 ;
  236. *ecro0 = ecro20 et ecro200 et ecro300 et ecro400 et ecro500 et ecro600 et ecro700 et ecro800 et ecro900 et ecro1000 ;
  237. *dess (trac0 et ecro0) ;
  238.  
  239. resp TBIB1 ;
  240. quit biblio ;
  241.  
  242. fins ;
  243. * Fin reference 1, 316L
  244.  
  245. *------------------------- 316L, Reference 2 --------------------------*
  246.  
  247. * Proprietes thermomecaniques du 316L :
  248. * Camille Cambon, Issam Bendaoud, Sebastien Rouquette, Fabien Soulie.
  249. * "Influence of the first weld bead on strain and stress states in wire+arc additive manufacturing”.
  250. * The 12th International Seminar ”Numerical Analysis of Weldability",
  251. * Institute for Materials Science, Joining and Forming (IMAT),Sep 2018, Seggau, Austria.
  252. * hal-01954354 https://hal.archives-ouvertes.fr/hal-01954354
  253.  
  254. si (ent1 ega 2) ;
  255.  
  256. *---- THERMIQUE ---- :
  257.  
  258. * K1 : Conductivite thermique (W/m/K) :
  259. lt1 = prog 20. 100. pas 100. 1000. 1200. 1400. ;
  260. lk1 = prog 14. 15.2 16.6 17.9 19. 20.6 21.8 23.1 24.3 26. 27.3 29.9 60. ;
  261. K1 = evol manu 'T' lt1 'K' lk1 ;
  262. TBIB1 . 'K' = K1 ;
  263.  
  264. * Enth1 : Enthalpie volumique (J/m^3)
  265. lt1 = prog 0. pas 100. 1500. ;
  266. lenth1 = prog 0.983 1.36 1.76 2.19 2.62 3.06 3.51 3.96 4.43 4.92 5.42 5.92 6.44 6.97 9.51 10.1 ;
  267. lenth1 = 1.e9 * lenth1 ;
  268. Enth1 = evol manu 'T' lt1 'ENTH' lenth1 ;
  269. TBIB1 . 'ENTH' = Enth1 ;
  270.  
  271. * Tfus1 : Temperature de fusion (degC) :
  272. Tfus1 = 1450. ;
  273. TBIB1 . 'TFUS' = Tfus1 ;
  274.  
  275. * Qlat1 : Chaleur latente volumique de fusion (J/m^3)
  276. Qlat1 = 2.e9 ;
  277. TBIB1 . 'QLAT' = Qlat1 ;
  278.  
  279. *---- MECANIQUE ---- :
  280.  
  281. * Ym1 : module de Young (Pa) :
  282. lt1 = prog 20. 100. pas 100. 1200. ;
  283. lym1 = prog 197000. 191500. 184000. 176500. 168000. 160000. 151500. 142500. 130000. 108000. 81500. 32000. 7400. ;
  284. lym1 = 1.e6 * lym1 ;
  285. ym1 = evol manu T lt1 youn lym1 ;
  286. TBIB1 . 'YOUN' = ym1 ;
  287.  
  288. * Sy1 : limite d'elasticite (Pa)
  289. lt1 = prog 20. 100. pas 100. 1400. ;
  290. lsy1 = prog 287. 237. 198. 172. 157. 151. 145. 136. 127. 115. 79. 38. 24. 18. 2. ;
  291. lsy1 = 1.e6 * lsy1 ;
  292. sy1 = evol manu T lt1 SIGY lsy1 ;
  293. TBIB1 . 'SIGY' = sy1 ;
  294.  
  295. * Alph1 : dilatation thermique (/K)
  296. lt1 = prog 20. 100. 300. 400. 600. 700. 800. 1000. 1400. ;
  297. lalph1 = prog 15.5 16.0 17.1 17.5 18.4 18.7 19.0 19.4 19.6 ;
  298. lalph1 = 1.e-6 * lalph1 ;
  299. alph1 = evol manu T lt1 ALPH lalph1 ;
  300. TBIB1 . 'ALPH' = alph1 ;
  301.  
  302. * H1 : module tangent d'elasticite (Pa)
  303. lt1 = prog 20. 700. pas 100. 1300. ;
  304. lh1 = prog 2400. 2400. 2350. 1500. 800. 725. 150. 10. ;
  305. lh1 = 1.e6 * lh1 ;
  306. H1 = evol manu T lt1 H lh1 ;
  307. TBIB1 . 'H' = H1 ;
  308.  
  309. * Ecro1 : Courbes d'ecrouissage :
  310. nb1 = dime lt1 ;
  311. lep1 = prog 0. 1. ;
  312. ecro1 = vide nuage ;
  313. ic1 = 0 ;
  314. repe b1 nb1 ;
  315. ic1 = ic1 + 1 ;
  316. si (ic1 ega 9) ; ic1 = ic1 + 1 ; fins ;
  317. ti1 = extr lt1 &b1 ;
  318. syi1 = ipol sy1 ti1 ;
  319. hi1 = ipol H1 ti1 ;
  320. lmsi1 = prog syi1 (syi1+hi1) ;
  321. ecroi1 = evol manu epeq lep1 smeq lmsi1 ;
  322. ecroi1 = ecroi1 coul ic1 ;
  323. ecro1 = ecro1 et (nuag comp T ti1 comp ECRO ecroi1) ;
  324. fin b1 ;
  325. TBIB1 . 'ECRO' = ecro1 ;
  326.  
  327. resp TBIB1 ;
  328. quit biblio ;
  329.  
  330. * Fin reference 2, 316L
  331.  
  332. * Message si REFE trop grand :
  333. sino ;
  334. mess '***** Pas de reference ' ' ' ent1 ' pour' ' ' mot1 ;
  335. saut 1 lign ;
  336. resp TBIB1 ;
  337. quit biblio ;
  338.  
  339. fins ;
  340.  
  341. * Fin mot-cle BIBLIO
  342. fins ;
  343.  
  344. *------------------------- Mot-cle non trouve -------------------------*
  345.  
  346. si (imot1 ega 0) ;
  347. mess '***** Pas de donnee disponible pour' ' ' mot1 ;
  348. saut 1 lign ;
  349. resp TBIB1 ;
  350. quit biblio ;
  351. fins ;
  352.  
  353. FINP ;
  354.  
  355.  
  356.  
  357.  

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