Télécharger la plaquette Cast3M


fr en

Télécharger biblio.procedur

Retour à la liste

Numérotation des lignes : 

  1. * BIBLIO    PROCEDUR  PASCAL    22/12/07    21:15:02     11514          
  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 : le 2e argument doit etre le mot-cle REFE' ;
  12.     quit biblio ;
  13.   fins ;
  14. sino ;
  15.   argu mot2/'ENTIER' ;
  16.   si (exis mot2) ;
  17.     erre '***** ERREUR : le 2e argument doit etre le mot-cle REFE' ;
  18.     quit biblio ;
  19.   fins ;
  20.   ent1     = 0 ;
  21. fins ;
  22.  
  23. * imot1  : indice du mot dans la base biblio, 0 si pas dans la base.
  24. imot1    = 0 ;
  25.  
  26. *----------------------------------------------------------------------*
  27. *                            MOT-CLE : 316L                            *
  28. *----------------------------------------------------------------------*
  29.  
  30. si (ega mot1 '316L') ;
  31.   imot1    = 1 ;
  32.  
  33. * Initialisation table de sortie :
  34.   TBIB1    = table ;
  35.  
  36. *------------------------------ 316L [1] ------------------------------*
  37.  
  38. * Proprietes thermomecaniques du 316L :
  39. * Lionel Depradeux.
  40. * "Simulation numerique du soudage - Acier 316L : validation sur cas
  41. * tests de complexite croissante".
  42. * These de doctorat.
  43. * Institut National des Sciences Appliquees de Lyon.
  44. * 2004.
  45.  
  46. * Reference par defaut :
  47.   si (ent1 ega 0) ; ent1 = 1 ; fins ;
  48.  
  49.   si (ent1 ega 1) ;
  50.  
  51. *---- THERMIQUE ---- :
  52.  
  53. *   K1     : Conductivite thermique (W/m/K) :
  54.     lt1      = prog 20. 100. pas 100. 1000. 1200. 1500. ;
  55.     lk1      = prog 14. 15.2 16.6 17.9 19. 20.6 21.8 23.1 24.3 26. 27.3 29.9 60. ;
  56.     K1       = evol manu 'T' lt1 'K' lk1 ;
  57.     TBIB1 . 'K' = K1 ;
  58.  
  59. *   Rho1   : Masse volumique (kg/m^3)
  60.     lt1      = prog 20. 100. pas 100. 1000. 1200. 1500. ;
  61.     lrho1    = prog 8000. 7970. 7940. 7890. 7850. 7800. 7750. 7700. 7660. 7610. 7570. 7450. 7200. ;
  62.     Rho1     = evol manu 'T' lt1 'RHO' lrho1 ;
  63.     TBIB1 . 'RHO' = Rho1 ;
  64.  
  65. *   Cp1    : Capacite calorifique massique (J/K/kg)
  66.     lt1      = prog 20. 100. pas 100. 1000. 1200. ;
  67.     lcp1     = prog 450. 490. 525. 545. 560. 570. 580. 595. 625. 650. 660. 677. ;
  68.     Cp1      = evol manu 'T' lt1 'C' lcp1 ;
  69.     TBIB1 . 'C' = Cp1 ;
  70.  
  71. *   Tfus1  : Temperaure de fusion (degC) :
  72.     Tfus1    = 1450. ;
  73.     TBIB1 . 'TFUS' = Tfus1 ;
  74.  
  75. *   Qlat1  : Chaleur latente massique de fusion (J/kg)
  76.     Qlat1    = 2.5e5 ;
  77.     TBIB1 . 'QLAT' = Qlat1 ;
  78.  
  79. *---- MECANIQUE ---- :
  80.  
  81. *   Nu1    : coefficient de Poisson
  82.     TBIB1 . 'NU' = 0.3 ;
  83.  
  84. *   Ym1    : module de Young (Pa) :
  85.     lt1      = prog 20. 100. pas 100. 1200. ;
  86.     lym1     = prog 197000. 191500. 184000. 176500. 168000. 160000. 151500. 142500. 130000. 108000. 81500. 32000. 7400. ;
  87.     lym1     = 1.e6 * lym1 ;
  88.     ym1      = evol manu T lt1 youn lym1 ;
  89.     TBIB1 . 'YOUN' = ym1 ;
  90.  
  91. *   Alph1  : dilatation thermique (/K)
  92.     lt1      = prog 20. 100. 300. 400. 600. 700. 800. 1000. 1400. ;
  93.     lalph1   = prog 15.5 16.0 17.1 17.5 18.4 18.7 19.0 19.4 19.6 ;
  94.     lalph1   = 1.e-6 * lalph1 ;
  95.     alph1    = evol manu T lt1 ALPH lalph1 ;
  96.     TBIB1 . 'ALPH' = alph1 ;
  97.  
  98. *   Sy1    : limite d'elasticite (Pa)
  99.     lt1      = prog 20. 100. pas 100. 1400. ;
  100.     lsy1     = prog 287. 237. 198. 172. 157. 151. 145. 136. 127. 115. 79. 38. 24. 18. 2. ;
  101.     lsy1     = 1.e6 * lsy1 ;
  102.     sy1      = evol manu T lt1 SIGY lsy1 ;
  103.     TBIB1 . 'SIGY' = sy1 ;
  104.  
  105. *   Trac1  : Courbes de traction (Pa)
  106. *   Ecro1  : Courbes d'ecrouissage (Pa)
  107.     Trac1    = vide nuage ;
  108.     Ecro1    = vide nuage ;
  109. *   A 20 degC :
  110.     T1       = 20. ;
  111.     lept1    = prog 0.3459 1.1599 5.2129 20.2847 ;
  112.     lept1    = 1.e-2 * lept1 ;
  113.     lsig1    = prog 287.33 314.94 419.42 560.97 ;
  114.     lsig1    = 1.e6 * lsig1 ;
  115.     trac20   = (evol manu 'EPST' lept1 'SIG' lsig1) coul 1 ;
  116.     Trac1    = Trac1 et (nuag comp 'T' T1 comp 'TRAC' trac20) ;
  117.     ymT1     = ipol ym1 T1 ;
  118.     lepp2    = lept1 - (lsig1 / ymT1) ;
  119.     Ecro20   = (evol manu 'EPSP' lepp2 'SIG' lsig1) coul 1 ;
  120.     Ecro20   = prol line Ecro20 0. ;
  121.     Ecro1    = Ecro1 et (nuag comp 'T' T1 comp 'ECRO' Ecro20) ;
  122. *   A 200 degC :
  123.     T1       = 200. ;
  124.     lept1    = prog 0.3077 1.1257 5.1816 20.2556 ;
  125.     lept1    = 1.e-2 * lept1 ;
  126.     lsig1    = prog 198.16 231.36 334.1 470.24 ;
  127.     lsig1    = 1.e6 * lsig1 ;
  128.     trac200  = (evol manu 'EPST' lept1 'SIG' lsig1) coul 2 ;
  129.     Trac1    = Trac1 et (nuag comp 'T' T1 comp 'TRAC' trac200) ;
  130.     ymT1     = ipol ym1 T1 ;
  131.     lepp2    = lept1 - (lsig1 / ymT1) ;
  132.     Ecro200  = (evol manu 'EPSP' lepp2 'SIG' lsig1) coul 2 ;
  133.     Ecro200  = prol line Ecro200 0. ;
  134.     Ecro1    = Ecro1 et (nuag comp 'T' T1 comp 'ECRO' Ecro200) ;
  135. *   A 300 degC :
  136.     T1       = 300. ;
  137.     lept1    = prog 0.2977 1.1137 5.1740 20.2558 ;
  138.     lept1    = 1.e-2 * lept1 ;
  139.     lsig1    = prog 172.47 200.62 307.12 451.42 ;
  140.     lsig1    = 1.e6 * lsig1 ;
  141.     trac300  = (evol manu 'EPST' lept1 'SIG' lsig1) coul 3 ;
  142.     Trac1    = Trac1 et (nuag comp 'T' T1 comp 'TRAC' trac300) ;
  143.     ymT1     = ipol ym1 T1 ;
  144.     lepp2    = lept1 - (lsig1 / ymT1) ;
  145.     Ecro300  = (evol manu 'EPSP' lepp2 'SIG' lsig1) coul 3 ;
  146.     Ecro300  = prol line Ecro300 0. ;
  147.     Ecro1    = Ecro1 et (nuag comp 'T' T1 comp 'ECRO' Ecro300) ;
  148. *   A 400 degC :
  149.     T1       = 400. ;
  150.     lept1    = prog 0.2939 1.1115 5.1755 20.2593 ;
  151.     lept1    = 1.e-2 * lept1 ;
  152.     lsig1    = prog 157.83 187.32 294.89 435.66 ;
  153.     lsig1    = 1.e6 * lsig1 ;
  154.     trac400  = (evol manu 'EPST' lept1 'SIG' lsig1) coul 4 ;
  155.     Trac1    = Trac1 et (nuag comp 'T' T1 comp 'TRAC' trac400) ;
  156.     ymT1     = ipol ym1 T1 ;
  157.     lepp2    = lept1 - (lsig1 / ymT1) ;
  158.     Ecro400  = (evol manu 'EPSP' lepp2 'SIG' lsig1) coul 4 ;
  159.     Ecro400  = prol line Ecro400 0. ;
  160.     Ecro1    = Ecro1 et (nuag comp 'T' T1 comp 'ECRO' Ecro400) ;
  161. *   A 500 degC :
  162.     T1       = 500. ;
  163.     lept1    = prog 0.2948 1.1120 5.1765 20.2620 ;
  164.     lept1    = 1.e-2 * lept1 ;
  165.     lsig1    = prog 151.65 179.12 282.34 419.14 ;
  166.     lsig1    = 1.e6 * lsig1 ;
  167.     trac500  = (evol manu 'EPST' lept1 'SIG' lsig1) coul 5 ;
  168.     Trac1    = Trac1 et (nuag comp 'T' T1 comp 'TRAC' trac500) ;
  169.     ymT1     = ipol ym1 T1 ;
  170.     lepp2    = lept1 - (lsig1 / ymT1) ;
  171.     Ecro500  = (evol manu 'EPSP' lepp2 'SIG' lsig1) coul 5 ;
  172.     Ecro500  = prol line Ecro500 0. ;
  173.     Ecro1    = Ecro1 et (nuag comp 'T' T1 comp 'ECRO' Ecro500) ;
  174. *   A 600 degC :
  175.     T1       = 600. ;
  176.     lept1    = prog 0.2960 1.1121 5.1730 20.2556 ;
  177.     lept1    = 1.e-2 * lept1 ;
  178.     lsig1    = prog 145.47 169.87 262.17 387.23 ;
  179.     lsig1    = 1.e6 * lsig1 ;
  180.     trac600  = (evol manu 'EPST' lept1 'SIG' lsig1) coul 6 ;
  181.     Trac1    = Trac1 et (nuag comp 'T' T1 comp 'TRAC' trac600) ;
  182.     ymT1     = ipol ym1 T1 ;
  183.     lepp2    = lept1 - (lsig1 / ymT1) ;
  184.     Ecro600  = (evol manu 'EPSP' lepp2 'SIG' lsig1) coul 6 ;
  185.     Ecro600  = prol line Ecro600 0. ;
  186.     Ecro1    = Ecro1 et (nuag comp 'T' T1 comp 'ECRO' Ecro600) ;
  187. *   A 700 degC :
  188.     T1       = 700. ;
  189.     lept1    = prog 0.2956 1.1128 5.1485 20.2389 ;
  190.     lept1    = 1.e-2 * lept1 ;
  191.     lsig1    = prog 136.2 160.7 211.65 340.49 ;
  192.     lsig1    = 1.e6 * lsig1 ;
  193.     trac700  = (evol manu 'EPST' lept1 'SIG' lsig1) coul 7 ;
  194.     Trac1    = Trac1 et (nuag comp 'T' T1 comp 'TRAC' trac700) ;
  195.     ymT1     = ipol ym1 T1 ;
  196.     lepp2    = lept1 - (lsig1 / ymT1) ;
  197.     Ecro700  = (evol manu 'EPSP' lepp2 'SIG' lsig1) coul 7 ;
  198.     Ecro700  = prol line Ecro700 0. ;
  199.     Ecro1    = Ecro1 et (nuag comp 'T' T1 comp 'ECRO' Ecro700) ;
  200. *   A 800 degC :
  201.     T1       = 800. ;
  202.     lept1    = prog 0.2980 1.1167 5.1529 20.2260 ;
  203.     lept1    = 1.e-2 * lept1 ;
  204.     lsig1    = prog 127. 151.75 198.75 293.75 ;
  205.     lsig1    = 1.e6 * lsig1 ;
  206.     trac800  = (evol manu 'EPST' lept1 'SIG' lsig1) coul 8 ;
  207.     Trac1    = Trac1 et (nuag comp 'T' T1 comp 'TRAC' trac800) ;
  208.     ymT1     = ipol ym1 T1 ;
  209.     lepp2    = lept1 - (lsig1 / ymT1) ;
  210.     Ecro800  = (evol manu 'EPSP' lepp2 'SIG' lsig1) coul 8 ;
  211.     Ecro800  = prol line Ecro800 0. ;
  212.     Ecro1    = Ecro1 et (nuag comp 'T' T1 comp 'ECRO' Ecro800) ;
  213. *   A 900 degC :
  214.     T1       = 900. ;
  215.     lept1    = prog 0.3068 1.1273 5.1493 20.1566 ;
  216.     lept1    = 1.e-2 * lept1 ;
  217.     lsig1    = prog 115.37 137.5 161.25 169.12 ;
  218.     lsig1    = 1.e6 * lsig1 ;
  219.     trac900  = (evol manu 'EPST' lept1 'SIG' lsig1) coul 11 ;
  220.     Trac1    = Trac1 et (nuag comp 'T' T1 comp 'TRAC' trac900) ;
  221.     ymT1     = ipol ym1 T1 ;
  222.     lepp2    = lept1 - (lsig1 / ymT1) ;
  223.     Ecro900  = (evol manu 'EPSP' lepp2 'SIG' lsig1) coul 11 ;
  224.     Ecro900  = prol line Ecro900 0. ;
  225.     Ecro1    = Ecro1 et (nuag comp 'T' T1 comp 'ECRO' Ecro900) ;
  226. *   A 1000 degC :
  227.     T1       = 1000. ;
  228.     lept1    = prog 0.2966 1.1009 5.1184 20.1227 ;
  229.     lept1    = 1.e-2 * lept1 ;
  230.     lsig1    = prog 78.75 82.25 96.5 100. ;
  231.     lsig1    = 1.e6 * lsig1 ;
  232.     trac1000 = (evol manu 'EPST' lept1 'SIG' lsig1) coul 10 ;
  233.     Trac1    = Trac1 et (nuag comp 'T' T1 comp 'TRAC' trac1000) ;
  234.     ymT1     = ipol ym1 T1 ;
  235.     lepp2    = lept1 - (lsig1 / ymT1) ;
  236.     Ecro1000 = (evol manu 'EPSP' lepp2 'SIG' lsig1) coul 10 ;
  237.     Ecro1000 = prol line Ecro1000 0. ;
  238.     Ecro1    = Ecro1 et (nuag comp 'T' T1 comp 'ECRO' Ecro1000) ;
  239.     TBIB1 . 'TRAC' = Trac1 ;
  240.     TBIB1 . 'ECRO' = Ecro1 ;
  241.  
  242. *trac0 = trac20 et trac200 et trac300 et trac400 et trac500 et trac600 et trac700 et trac800 et trac900 et trac1000 ;
  243. *ecro0 = ecro20 et ecro200 et ecro300 et ecro400 et ecro500 et ecro600 et ecro700 et ecro800 et ecro900 et ecro1000 ;
  244. *dess (trac0 et ecro0) ;
  245.     TBIB1 . 'CREATEUR' = mot 'BIBLIO' ;
  246.     TBIB1 . 'REFERENCE' = chai mot1 '_' ent1 ;
  247.  
  248.     resp TBIB1 ;
  249.     quit biblio ;
  250.  
  251.   fins ;
  252. * Fin [1], "316L" 
  253.  
  254. *------------------------------ 316L [2] ------------------------------*
  255.  
  256. * Proprietes thermomecaniques du 316L :
  257. * Camille Cambon, Issam Bendaoud, Sebastien Rouquette, Fabien Soulie.
  258. * "Influence of the first weld bead on strain and stress states in wire+arc additive manufacturing”.
  259. * The 12th International Seminar ”Numerical Analysis of Weldability",
  260. * Institute for Materials Science, Joining and Forming (IMAT),Sep 2018, Seggau, Austria.
  261. * hal-01954354 https://hal.archives-ouvertes.fr/hal-01954354
  262.  
  263.   si (ent1 ega 2) ;
  264.  
  265. *---- THERMIQUE ---- :
  266.  
  267. *   K1     : Conductivite thermique (W/m/K) :
  268.     lt1      = prog 20. 100. pas 100. 1000. 1200. 1400. ;
  269.     lk1      = prog 14. 15.2 16.6 17.9 19. 20.6 21.8 23.1 24.3 26. 27.3 29.9 60. ;
  270.     K1       = evol manu 'T' lt1 'K' lk1 ;
  271.     TBIB1 . 'K' = K1 ;
  272.  
  273. *   Enth1  : Enthalpie volumique (J/m^3)
  274.     lt1      = prog 0. pas 100. 1500. ;
  275.     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 ;
  276.     lenth1   = 1.e9 * lenth1 ;
  277.     Enth1    = evol manu 'T' lt1 'ENTH' lenth1 ;
  278.     TBIB1 . 'ENTH' = Enth1 ;
  279.  
  280. *   Tfus1  : Temperature de fusion (degC) :
  281.     Tfus1    = 1450. ;
  282.     TBIB1 . 'TFUS' = Tfus1 ;
  283.  
  284. *   Qlat1  : Chaleur latente volumique de fusion (J/m^3)
  285.     Qlat1    = 2.e9 ;
  286.     TBIB1 . 'QLAT' = Qlat1 ;
  287.  
  288. *---- MECANIQUE ---- :
  289.  
  290. *   Ym1    : module de Young (Pa) :
  291.     lt1      = prog 20. 100. pas 100. 1200. ;
  292.     lym1     = prog 197000. 191500. 184000. 176500. 168000. 160000. 151500. 142500. 130000. 108000. 81500. 32000. 7400. ;
  293.     lym1     = 1.e6 * lym1 ;
  294.     ym1      = evol manu T lt1 youn lym1 ;
  295.     TBIB1 . 'YOUN' = ym1 ;
  296.  
  297. *   Sy1    : limite d'elasticite (Pa)
  298.     lt1      = prog 20. 100. pas 100. 1400. ;
  299.     lsy1     = prog 287. 237. 198. 172. 157. 151. 145. 136. 127. 115. 79. 38. 24. 18. 2. ;
  300.     lsy1     = 1.e6 * lsy1 ;
  301.     sy1      = evol manu T lt1 SIGY lsy1 ;
  302.     TBIB1 . 'SIGY' = sy1 ;
  303.     ltsy1    = lt1 ;
  304.  
  305. *   Alph1  : dilatation thermique (/K)
  306.     lt1      = prog 20. 100. 300. 400. 600. 700. 800. 1000. 1400. ;
  307.     lalph1   = prog 15.5 16.0 17.1 17.5 18.4 18.7 19.0 19.4 19.6 ;
  308.     lalph1   = 1.e-6 * lalph1 ;
  309.     alph1    = evol manu T lt1 ALPH lalph1 ;
  310.     TBIB1 . 'ALPH' = alph1 ;
  311.  
  312. *   H1     : module tangent d'elasticite (Pa)
  313.     lt1      = prog 20. 700. pas 100. 1300. ;
  314.     lh1      = prog 2400. 2400. 2350. 1500. 800. 725. 150. 10. ;
  315.     lh1      = 1.e6 * lh1 ;
  316.     H1       = evol manu T lt1 H lh1 ;
  317.     TBIB1 . 'H' = H1 ;
  318.  
  319. *   Ecro1  : Courbes d'ecrouissage :
  320.     lt1      = ((lt1 et ltsy1) uniq) ordo ;
  321.     nb1      = dime lt1 ;
  322.     lep1     = prog 0. 1. ;
  323.     ecro1    = vide nuage ;
  324.     ic1      = 0 ;
  325.     repe b1 nb1 ;
  326.       ic1      = ic1 + 1 ;
  327.       si (ic1 ega 9) ; ic1 = ic1 + 1 ; fins ;
  328.       ti1      = extr lt1 &b1 ;
  329.       syi1     = ipol sy1 ti1 ;
  330.       hi1      = ipol H1  ti1 ;
  331.       lmsi1    = prog syi1 (syi1+hi1) ;
  332.       ecroi1   = evol manu epeq lep1 smeq lmsi1 ;
  333.       ecroi1   = ecroi1 coul ic1 ;
  334.       ecro1    = ecro1 et (nuag comp T ti1 comp ECRO ecroi1) ;
  335.     fin b1 ;
  336.     TBIB1 . 'ECRO' = ecro1 ;
  337.  
  338.     TBIB1 . 'CREATEUR' = mot 'BIBLIO' ;
  339.     TBIB1 . 'REFERENCE' = chai mot1 '_' ent1 ;
  340.  
  341.     resp TBIB1 ;
  342.     quit biblio ;
  343.  
  344.   fins ;
  345. * Fin [2], "316L"
  346.  
  347. *------------------------------ 316L [4] ------------------------------*
  348.  
  349.   si (ent1 ega 4) ;
  350.  
  351. *--- Conductivite thermique  [w/mm.K] :
  352.   LT1      = prog 20 100 200 300 400 500 600 700 800 900 1000 1100 1200 1300 1400 ;
  353.   LK1      = prog 0.0141 0.0152 0.0167 0.0181 0.0195 0.0209 0.0224 0.0238 0.0252 0.0267 0.0281 0.0295 0.0309 0.0324 0.0338 ;
  354.   EVK1     = EVOL MANU 'T' LT1 'K' LK1 ;
  355.   comm /!\  conversion [w/mm.K] -> [w/m.K] ;
  356.   K_316L   = 1.e3 * EVK1 ;
  357.   TBIB1 . 'K' = K_316L ;
  358.  
  359.  
  360.  
  361. *--- Capacite calorifique massique [J/Kg.K] :
  362.   LT1      = prog 20 100 200 300 400 500 600 700 800 900 1000 1100 ;
  363.   LCP1     = prog 492 502 514 526 538 550 562 575 587 599 611 623  ;
  364.   EVCP1    = evol manu 'T' LT1 'C' LCP1 ;
  365.   C_316L   = EVCP1 ;
  366.   TBIB1 . 'C' = C_316L ;
  367.  
  368.  
  369.  
  370. *--- Module d'Young (GPa) :
  371.   LT1      = prog 20 100 200 300 400 500 600 700 800 900 1000 1100 1200 1300 1400 4000 ;
  372.   LYM1     = prog 195.6 191.2 185.7 179.6 172.6 164.5 155.0 144.1 131.4 116.8 100.0 80.00 57.00 30.00 2.000 2.000 ;
  373.   YM_316L = evol manu 'T' LT1 'YOUN' LYM1 ;
  374.   comm /!\ conversion [MPa] -> [Pa] ;
  375.   TBIB1 . 'YOUN' = 1.e9 * YM_316L ;
  376.  
  377.  
  378.  
  379. *--- Coefficient de Poisson :
  380.   TBIB1 . 'NU'   = 0.29 ;
  381.  
  382.  
  383.  
  384. *--- Coeficient de dilatation thermique [1/K] :
  385.   LT1      = prog 20 100 200 300 400 500 600 700 800 900 1000 1100 1200 1300 1400 4000 ;
  386.   LALPHA1  = prog 1.46E-05 1.54E-05 1.62E-05 1.69E-05 1.74E-05 1.78E-05 1.81E-05 1.84E-05 1.87E-05 1.90E-05 1.93E-05 1.95E-05 1.98E-05 2.00E-05 2.02E-05 2.02E-05 ;
  387.   AL_316L = evol manu 'T' LT1 'ALPH' LALPHA1 ;
  388.   TBIB1 . 'ALPH' = AL_316L ;
  389.  
  390.  
  391.  
  392. *--- limite élastique [MPa]
  393.   LT1      = prog 20 275 550 750 900 1000 1100 1400 ;
  394.   LSY1     = prog 125.6 97.60 90.90 71.40 66.20 31.82 19.73 2.100 ;
  395.   EVSY1    = evol manu 'T' LT1 'SIGY' LSY1 ;
  396.   comm /!\ conversion [MPa] -> [Pa] ;
  397.   SY_316L = 1.e6 * EVSY1 ;
  398.   TBIB1 . 'SIGY' = SY_316L ;
  399.  
  400.  
  401.  
  402. *--- Coefficients ecrouissage cinematique modele EP Chaboche :
  403. *    LCiMur : coefficient Ci article Muransky
  404. *    LGiMur : coefficient Gammai article Muransky
  405. *    LCiC3M : coefficient Ci loi EP Chaboche Cast3M
  406. *    LAiC3M : coefficient Ai loi EP Chaboche Cast3M
  407.   LT1      = prog 20 275 550 750 900 1000 1100 1400 ;
  408.   LC1Mur   = prog 156435.0 100631.0 64341.00 56232.00 0.050000 0 0 0 ;
  409.   LG1Mur   = 1410.9 ;
  410.   LC2Mur   = prog 6134.0 5568.0 5227.0 4108.0 292.00 0 0 0 ;
  411.   LG2Mur   = 47.19 ;
  412.   comm /!\ conversion [MPa] -> [Pa] ;
  413.   LC1Mur   = 1.e6 * LC1Mur ;
  414.   LC2Mur   = 1.e6 * LC2Mur ;
  415.   LC1C3M   = LG1Mur ;
  416.   LC2C3M   = LG2Mur ;
  417.   LA1C3M   = LC1Mur / LG1Mur ;
  418.   LA2C3M   = LC2Mur / LG2Mur ;
  419.   EVA1     = evol manu 'T' LT1 'A1' LA1C3M ;
  420.   EVA2     = evol manu 'T' LT1 'A2' LA2C3M ;
  421.   R0_316L = evol manu 'T' (extr SY_316L absc) 'R0' (extr SY_316L ordo) ;
  422.   TBIB1 . 'R0' = R0_316L ;
  423.   TBIB1 . 'A1' = EVA1 ;
  424.   TBIB1 . 'C1' = LC1C3M ;
  425.   TBIB1 . 'A2' = EVA2 ;
  426.   TBIB1 . 'C2' = LC2C3M ;
  427.  
  428.  
  429.  
  430. *--- Coefficients ecrouissage isotrope modele EP Chaboche :
  431.   LT1      = prog 20 275 550 750 900 1000 1100 1400 ;
  432.   LRM1     = prog 153.4 154.7 150.6 57.90 0 0 0 0 ;
  433.   EVRM1    = evol manu 'T' LT1 'RM' LRM1 ;
  434.   comm /!\ conversion [MPa] -> [Pa] ;
  435.   EVRM1    = 1.e6 * EVRM1 ;
  436.   TBIB1 . 'RM' = EVRM1 + (evol manu 'T' (extr SY_316L absc) 'RM' (extr SY_316L ordo)) ;
  437.   TBIB1 . 'B'  = 6.9 ;
  438.  
  439.  
  440.  
  441.     TBIB1 . 'CREATEUR' = mot 'BIBLIO' ;
  442.     TBIB1 . 'REFERENCE' = chai mot1 '_' ent1 ;
  443.  
  444.     resp TBIB1 ;
  445.     quit biblio ;
  446.  
  447.   fins ;
  448. * Fin [4], "316L"
  449.  
  450. * Reference ent1 n'existe pas pour mot1
  451.   erre (chai '*****  La reference' ' ' ent1 ' ne concerne pas' ' ' mot1) ;
  452.   quit biblio ;
  453.  
  454. * Fin mot-cle "316L"
  455. fins ;
  456.  
  457. *----------------------------------------------------------------------*
  458. *                           MOT-CLE : 9Cr ODS                          *
  459. *----------------------------------------------------------------------*
  460.  
  461. si (ega mot1 '9Cr ODS') ;
  462.   imot1    = 2 ;
  463.  
  464. * Initialisation table de sortie :
  465.   TBIB1    = table ;
  466.  
  467. *----------------------------- 9Cr ODS [3] ----------------------------*
  468.  
  469. * Brendan Le Gloannec.
  470. * "Modifications microstructurales sous sollicitations thermomécaniques sévères :
  471. * application au soudage par résistance des gaines de combustibles en aciers ODS".
  472. * Thèse de doctorat.
  473. * Université de Bordeaux.
  474. * 2016.
  475.  
  476. * Numero reference bib. par defaut :
  477.   si (ent1 ega 0) ; ent1 = 3 ; fins ;
  478.  
  479.   si (ent1 ega 3) ;
  480.  
  481.  
  482. *--- Masse volumique [Kg/m3] :
  483.   RHO_9ODS =      7760. ;
  484.   TBIB1 . 'RHO' = RHO_9ODS ;
  485.  
  486.  
  487.  
  488. *--- Conductivite thermique  [w/mm.K] :
  489.   CONDTH_9ODS1  = EVOL MANU 'T' (PROG  2.5e+01  5e+01     1e+02    1.5e+02
  490.                                        2e+02    2.5e+02   3e+02    3.5e+02
  491.                                        4e+02    4.5e+02   5e+02    5.5e+02)
  492.                             'K' (PROG  2.3e-02  2.27e-02  2.23e-02
  493.                                        2.36e-02 2.49e-02  2.54e-02  2.63e-02
  494.                                        2.6e-02  2.65e-02  2.71e-02  2.74e-02
  495.                                        2.72e-02) ;
  496.   CONDTH_9ODS2  = EVOL MANU 'T' (PROG 6e+02     6.5e+02   7e+02    7.5e+02
  497.                                       8e+02     8.5e+02   9e+02    9.5e+02
  498.                                       1e+03     1.05e+03  1.1e+03  1.15e+03
  499.                                       1.18e+03  1.45e+03)
  500.                             'K' (PROG 2.73e-02  2.65e-02    2.57e-02
  501.                                         2.14e-02  2.53e-02  4.44e-02  2.9e-02
  502.                                         2.65e-02  2.65e-02  2.65e-02  2.73e-02
  503.                                         2.82e-02  2.88e-02  3.28e-02) ;
  504.   comm /!\  conversion [w/mm.K] -> [w/m.K] ;
  505.   CONDTH_9ODS   = 1.e3 * (CONDTH_9ODS1 + CONDTH_9ODS2) ;
  506.   TBIB1 . 'K' = CONDTH_9ODS ;
  507.  
  508.  
  509.  
  510. *--- Capacite calorifique massique [J/Kg.K] :
  511.   CAPAC_9ODS1 = EVOL MANU 'T' (PROG 2.5e+01  5.0e+01  1.0e+02  1.5e+02  2.0e+02
  512.                                     2.5e+02  3.0e+02  3.5e+02  4.0e+02
  513.                                     4.5e+02  5.0e+02  5.5e+02  7.0e+02
  514.                                     7.5e+02  8.0e+02  8.5e+02  9.0e+02)
  515.                           'C' (PROG 4.4e+02  4.4e+02  4.4e+02  4.75e+02 5.1e+02
  516.                                     5.3e+02  5.6e+02  5.7e+02  6.0e+02
  517.                                     6.4e+02  6.8e+02  7.2e+02  9.5e+02
  518.                                     9.6e+02  7.8e+02  1.14e+03 6.8e+02) ;
  519.   CAPAC_9ODS2 = EVOL MANU 'T' (PROG 9.5e+02   1.0e+03   1.05e+03  1.10e+03
  520.                                     1.15e+03  1.45e+03  1.46e+03  1.47e+03
  521.                                     1.48e+03  1.5e+03)
  522.                           'C' (PROG 6.2e+02   6.2e+02   6.2e+02   6.4e+02
  523.                                     6.6e+02   7.95e+02  9.94e+03  9.942e+03
  524.                                     8.38e+02  8.38e+02) ;
  525.   CAPAC_9ODS  = (CAPAC_9ODS1 + CAPAC_9ODS2) ;
  526.   TBIB1 . 'C' = CAPAC_9ODS ;
  527.  
  528.  
  529.  
  530. *--- Conductivite electrique 9Cr ODS [Ohm.mm]^-1 :
  531.   CONDEL_9ODS  = EVOL MANU 'T'  (PROG  2.2e+01   6.7e+02  1.5e+03   1.6e+03)
  532.                            'KD' (PROG  0.2e+04   1.0e+03  0.606e+03 0.606e+03) ;
  533.   comm /!\ conversion [Ohm.mm]^-1 -> [Ohm.m]^-1 ;
  534.   TBIB1 . 'KE' = 1.e3*CONDEL_9ODS ;
  535.  
  536.  
  537.  
  538. *--- Module d'Young (MPa) :
  539.   YONG_9CrODS = EVOL MANU 'T' (PROG 2.e+01    1.2e+02  2.2e+02  4.0e+02   4.7e+02   6.15e+02  6.6e+02  7.6e+02
  540.                                     8.3e+02   9.6e+02  1.e+03   1.5e+03)
  541.  
  542.                        'YOUN' (PROG  2.16e+5  2.1e+5   2.05e+5  1.9e+5    1.8e+5    1.6e+5    1.5e+5   1.2e+5
  543.                                      1.15e+5  1.15e+5  1.14e+5  1.     ) ;
  544.   comm /!\ conversion [MPa] -> [Pa] ;
  545.   TBIB1 . 'YOUN' = 1.e6 * YONG_9CrODS ;
  546.  
  547.  
  548.  
  549. *--- Coefficient de Poisson :
  550.   NU_9CrODS = EVOL MANU 'T' (PROG 2.e+01   1.45e+03  1.5e+03 )
  551.                        'NU' (PROG 2.9e-01  2.9e-01   0.49e-01) ;
  552.   comm /!\ : on prend une valeur constante car le choix de modelisation du comportement en phase liquide ;
  553. *       n'est pas du ressort des donnees :
  554.   TBIB1 . 'NU'   = 0.29 ;
  555.  
  556.  
  557.  
  558. *--- Coeficient de dilatation thermique [1/K] :
  559.   ALPH_9CrODS1 = EVOL MANU 'T'(PROG  2.e+01  1.e+02    1.5e+02  2.e+02  2.5e+02  3.e+02   3.5e+02  4.e+02  4.5e+02
  560.                                      5.e+02  5.5e+02 )
  561.                       'ALPH' (PROG 1.01e-05  1.07e-05  1.1e-05   1.13e-05  1.17e-05  1.21e-05  1.24e-05 1.28e-05
  562.                                    1.31e-05  1.36e-05  1.39e-05  );
  563.   ALPH_9CrODS2 = EVOL MANU 'T'(PROG  6.e+02   6.5e+02  7.e+02   7.5e+02  8.e+02   8.5e+02 9.5e+02
  564.                                      1.e+03  1.05e+03  1.1e+03  1.5e+03)
  565.                       'ALPH' (PROG 1.42e-05  1.46e-05  1.5e-05  1.53e-05 1.57e-05
  566.                                    1.6e-05   -3.e-06   2.4e-05  2.4e-05   2.4e-05   2.91e-05) ;
  567.   ALPH_9CrODS =   ALPH_9CrODS1+ALPH_9CrODS2;
  568.   TBIB1 . 'ALPH' = ALPH_9CrODS ;
  569.  
  570.  
  571.  
  572. *--- limite élastique [MPa]
  573.   KK_9CrODS = EVOL MANU 'T'  (PROG   2.e+01     4.e+02    5.e+02   6.e+02   6.5e+02   7.5e+02  9.5e+02  1.25e+03 1.45e+03 1.5e+03)
  574.                         'KK' (PROG   1.083e+03  8.87e+02  7.03e+02 4.32e+02 1.96e+02  1.65e+02  1.e+02  2.2e+01  1.        1.);
  575.   comm /!\ conversion [MPa] -> [Pa] ;
  576.   KK_9CrODS = 1.e6 * KK_9CrODS ;
  577.   TBIB1 . 'KK' = KK_9CrODS ;
  578.  
  579.  
  580.  
  581. *--- Coefficient d'ecrouissage isotrope a saturation Rmin & Rmax (Rmax = KK + Rmin) [MPa] lisse :
  582.   RMIN_9CrODS =  EVOL MANU 'T'    (PROG  2.0e+01  4.e+02   5.e+02   6.e+02   7.5e+02   9.5e+02  1.45e+03)
  583.                            'RMAX' (PROG  1.33e+02 6.6e+01  4.7e+01  4.3e+01  3.8e+01   1.       0.) ;
  584.   comm /!\ conversion [MPa] -> [Pa] ;
  585.   RMIN_9CrODS = 1.e6 * RMIN_9CrODS ;
  586.   TBIB1 . 'RMIN' = RMIN_9CrODS ;
  587.   KK_9CrODS      = EVOL MANU 'T' (EXTR KK_9CrODS ABSC) 'RMAX' (EXTR KK_9CrODS ORDO) ;
  588.   TBIB1 . 'RMAX' = KK_9CrODS + RMIN_9CrODS ;
  589.  
  590.  
  591.  
  592. *--- Coefficient d'ecrouissage isotrope BR [MPa]
  593.   BR_9CrODS =  EVOL MANU 'T'  (PROG  2.e+01    4.e+02   5.e+2    6.e+02    6.5e+02  7.5e+2    9.5e+02  1.50e+03)
  594.                          'BR' (PROG  1.24e+02  1.6e+02  2.53e+02  1.26e+02 6.7e+01  1.08e+02  0.       0.      ) ;
  595.   comm /!\ conversion [MPa] -> [Pa] ;
  596.   BR_9CrODS = 1.e6 * BR_9CrODS ;
  597.   TBIB1 . 'BR' = BR_9CrODS ;
  598.  
  599.  
  600.  
  601. *--- Exposant viscosite N :
  602.   N_9CrODS = EVOL MANU 'T' (PROG  2.e+01    5.e+02   6.5e+02   7.5e+02  9.5e+02   1.25e+03  1.45e+03  )
  603.                        'N' (PROG  6.43e+00  6.43e+00 4.32e+00  2.77e+00 1.81e+01  1.136e+01 1.1360e+01) ;
  604.   TBIB1 . 'N'    = N_9CrODS ;
  605.  
  606.  
  607.  
  608. *--- Parametre viscosite K0 [MPa.s]
  609.   K0_9CrODS = EVOL MANU 'T'  (PROG   2.e+01    5.e+02    6.5e+02   7.5e+02   9.5e+02    1.25e+03  1.50e+03)
  610.                         'K0' (PROG   3.11e+02  3.11e+02  3.03e+02  4.61e+02  1.412e+02  1.02e+02  1.02e+02 ) ;
  611.   comm /!\ conversion [MPa] -> [Pa] ;
  612.   K0_9CrODS = 1.e6 * K0_9CrODS ;
  613.   TBIB1 . 'K0' = K0_9CrODS ;
  614.  
  615.  
  616.  
  617.     TBIB1 . 'CREATEUR' = mot 'BIBLIO' ;
  618.     TBIB1 . 'REFERENCE' = chai mot1 '_' ent1 ;
  619.  
  620.     resp TBIB1 ;
  621.     quit biblio ;
  622.  
  623.   fins ;
  624. *--- Fin [3], "9Cr ODS"
  625.  
  626.  
  627.  
  628. * Numero reference bib. ent1 n'existe pas pour mot1
  629.   erre (chai '*****  La reference' ' ' ent1 ' ne concerne pas' ' ' mot1) ;
  630.   quit biblio ;
  631.  
  632. * Fin mot-cle "9Cr ODS"
  633. fins ;
  634.  
  635. *----------------------------------------------------------------------*
  636. *                           MOT-CLE : Inconel 82                       *
  637. *----------------------------------------------------------------------*
  638.  
  639. si (ega mot1 'Inconel 82') ;
  640.   imot1    = 3 ;
  641.  
  642. * Initialisation table de sortie :
  643.   TBIB1    = table ;
  644.  
  645. *---------------------------- Inconel 82 [4] --------------------------*
  646. * Reference : voir notice BIBLIO
  647.  
  648. * Numero reference bib. par defaut :
  649.   si (ent1 ega 0) ; ent1 = 4 ; fins ;
  650.  
  651.   si (ent1 ega 4) ;
  652.  
  653.  
  654. *--- Masse volumique [Kg/m3] :
  655. *  RHO1     =      'XX' ;
  656. *  TBIB1 . 'RHO' = RHO1 ;
  657.  
  658.  
  659.  
  660. *--- Conductivite thermique  [w/mm.K] :
  661.   LT1      = prog 27. 77. 127. 177. 227. PAS 100. 1327. ;
  662.   LK1      = prog 0.0142 0.0151 0.0160 0.0168 0.0177 0.0195 0.0213 0.0230 0.0248 0.0265 0.0283 0.0301 0.0318 0.0336 0.0353 0.0371 ;
  663.   EVK1     = EVOL MANU 'T' LT1 'K' LK1 ;
  664.   comm /!\  conversion [w/mm.K] -> [w/m.K] ;
  665.   K_INCO82   = 1.e3 * EVK1 ;
  666.   TBIB1 . 'K' = K_INCO82 ;
  667.  
  668.  
  669.  
  670. *--- Capacite calorifique massique [J/Kg.K] :
  671.   LT1      = prog 20. 100. PAS 100. 900. 1204. 1316. ;
  672.   LCP1     = prog 444 465 486 502 519 536 578 595 611 628 733 766 ;
  673.   EVCP1    = evol manu 'T' LT1 'C' LCP1 ;
  674.   C_INCO82 = EVCP1 ;
  675.   TBIB1 . 'C' = C_INCO82 ;
  676.  
  677.  
  678.  
  679. *--- Module d'Young (GPa) :
  680.   LT1      = prog 21 93 149 204 260 316 371 427 482 538 600 PAS 100. 1000. 1150 1260 1316 1371 ;
  681.   LYM1     = prog 213.7 208.2 206.2 203.4 200.0 197.9 194.4 190.3 186.2 182.0 180.0 172.0 164.0 154.0 143.0 115.0 78.30 39.30 2.000 ;
  682.   YM_INCO82 = evol manu 'T' LT1 'YOUN' LYM1 ;
  683.   comm /!\ conversion [MPa] -> [Pa] ;
  684.   TBIB1 . 'YOUN' = 1.e9 * YM_INCO82 ;
  685.  
  686.  
  687.  
  688. *--- Coefficient de Poisson :
  689.   TBIB1 . 'NU'   = 0.31 ;
  690.  
  691.  
  692.  
  693. *--- Coeficient de dilatation thermique [1/K] :
  694.   LT1      = prog 38 66 93 121 149 178 204 232 260 288 316 343 371 399 428 500 PAS 100. 900. 1093 1371 ;
  695.   LALPHA1  = prog 1.24E-05 1.27E-05 1.30E-05 1.32E-05 1.33E-05 1.35E-05 1.36E-05 1.38E-05 1.39E-05 1.40E-05 1.41E-05 1.42E-05 1.43E-05 1.44E-05 1.45E-05 1.49E-05 1.53E-05 1.58E-05 1.61E-05 1.64E-05 1.70E-05 1.78E-05 ;
  696.   AL_INCO82 = evol manu 'T' LT1 'ALPH' LALPHA1 ;
  697.   TBIB1 . 'ALPH' = AL_INCO82 ;
  698.  
  699.  
  700.  
  701. *--- limite élastique [MPa]
  702.   LT1      = prog 22 316 538 760 982 1200 1400 ;
  703.   LSY1     = prog 206.8 161.0 171.0 144.7 34.80 21.00 2.000 ;
  704.   EVSY1    = evol manu 'T' LT1 'SIGY' LSY1 ;
  705.   comm /!\ conversion [MPa] -> [Pa] ;
  706.   SY_INCO82 = 1.e6 * EVSY1 ;
  707.   TBIB1 . 'SIGY' = SY_INCO82 ;
  708.  
  709.  
  710.  
  711. *--- Coefficients ecrouissage cinematique modele EP Chaboche :
  712. *    LCiMur : coefficient Ci article Muransky
  713. *    LGiMur : coefficient Gammai article Muransky
  714. *    LCiC3M : coefficient Ci loi EP Chaboche Cast3M
  715. *    LAiC3M : coefficient Ai loi EP Chaboche Cast3M
  716.   LT1      = prog 22.000 316.00 538.00 760.00 982.00 1200.0 1400.0 ;
  717.   LC1Mur   = prog 2.54865E+05 1.92628E+05 1.04132E+05 1.31679E+05 1.79993E+05 0.0000 0.0000 ;
  718.   LG1Mur   = 2496.2 ;
  719.   LC2Mur   = prog 1731.0 2304.6 2571.2 2668.4 623.60 0.0000 0.0000 ;
  720.   LG2Mur   = 3.9560 ;
  721.   comm /!\ conversion [MPa] -> [Pa] ;
  722.   LC1Mur   = 1.e6 * LC1Mur ;
  723.   LC2Mur   = 1.e6 * LC2Mur ;
  724.   LC1C3M   = LG1Mur ;
  725.   LC2C3M   = LG2Mur ;
  726.   LA1C3M   = LC1Mur / LG1Mur ;
  727.   LA2C3M   = LC2Mur / LG2Mur ;
  728.   EVA1     = evol manu 'T' LT1 'A1' LA1C3M ;
  729.   EVA2     = evol manu 'T' LT1 'A2' LA2C3M ;
  730.   R0_INCO82 = evol manu 'T' (extr SY_INCO82 absc) 'R0' (extr SY_INCO82 ordo) ;
  731.   TBIB1 . 'R0' = R0_INCO82 ;
  732.   TBIB1 . 'A1' = EVA1 ;
  733.   TBIB1 . 'C1' = LC1C3M ;
  734.   TBIB1 . 'A2' = EVA2 ;
  735.   TBIB1 . 'C2' = LC2C3M ;
  736.  
  737.  
  738.  
  739. *--- Coefficients ecrouissage isotrope modele EP Chaboche :
  740.   LT1      = prog 22.000 316.00 538.00 760.00 982.00 1200.0 1400.0 ;
  741.   LRM1     = prog 83.5    75.3   83.2    0.     0.      0.     0. ;
  742.   EVRM1    = evol manu 'T' LT1 'RM' LRM1 ;
  743.   comm /!\ conversion [MPa] -> [Pa] ;
  744.   EVRM1    = 1.e6 * EVRM1 ;
  745.   TBIB1 . 'RM' = EVRM1 + (evol manu 'T' (extr SY_INCO82 absc) 'RM' (extr SY_INCO82 ordo)) ;
  746.   TBIB1 . 'B'  = 20. ;
  747.  
  748.  
  749.  
  750.     TBIB1 . 'CREATEUR' = mot 'BIBLIO' ;
  751.     TBIB1 . 'REFERENCE' = chai mot1 '_' ent1 ;
  752.  
  753.     resp TBIB1 ;
  754.     quit biblio ;
  755.  
  756.   fins ;
  757. *--- Fin [4], "Inconel 82"
  758.  
  759.  
  760.  
  761. * Numero reference bib. ent1 n'existe pas pour mot1
  762.   erre (chai '*****  La reference' ' ' ent1 ' ne concerne pas' ' ' mot1) ;
  763.   quit biblio ;
  764.  
  765. * Fin mot-cle "Inconel 82"
  766. fins ;
  767.  
  768. *----------------------------------------------------------------------*
  769. *                            MOT-CLE : A508                            *
  770. *----------------------------------------------------------------------*
  771.  
  772. si ((ega mot1 'A508') ou (ega mot1 '16MND5')) ;
  773.   imot1    = 4 ;
  774.  
  775. * Initialisation table de sortie :
  776.   TBIB1    = table ;
  777.  
  778. * Numero reference bib. par defaut :
  779.   si (ent1 ega 0) ; ent1 = 4 ; fins ;
  780.  
  781. *------------------------------ A508 [4] ------------------------------*
  782. * Reference : voir notice BIBLIO
  783.  
  784.   si (ent1 ega 4) ;
  785.  
  786.  
  787. *--- Masse volumique [Kg/m3] :
  788. *  RHO1     =      'XX' ;
  789. *  TBIB1 . 'RHO' = RHO1 ;
  790.  
  791.  
  792.  
  793. *--- Conductivite thermique  [w/mm.K] :
  794.   LT1      = prog 25.000 50.000 100.00 150.00 200.00 250.00 300.00 350.00 400.00 450.00 500.00 550.00 600.00 650.00 700.00 750.00 800.00 850.00 900.00 950.00 1000.0 1100.0 1200.0 1300.0 1400.0 ;
  795.   LK1      = prog 3.84000E-02 3.94000E-02 4.11000E-02 4.12000E-02 4.09000E-02 4.04000E-02 3.95000E-02 3.91000E-02 3.83000E-02 3.74000E-02 3.64000E-02 3.58000E-02 3.48000E-02 3.42000E-02 3.42000E-02 5.90000E-02 3.67000E-02 2.69000E-02 2.84000E-02 2.78000E-02 2.83000E-02 2.95000E-02 3.09000E-02 3.24000E-02 3.38000E-02 ;
  796.   EVK1     = EVOL MANU 'T' LT1 'K' LK1 ;
  797.   comm /!\  conversion [w/mm.K] -> [w/m.K] ;
  798.   K_A508   = 1.e3 * EVK1 ;
  799.   TBIB1 . 'K' = K_A508 ;
  800.  
  801.  
  802.  
  803. *--- Capacite calorifique massique [J/Kg.K] :
  804.   LT1      = prog 25.000 50.000 100.00 150.00 200.00 250.00 300.00 350.00 400.00 450.00 500.00 550.00 600.00 650.00 700.00 750.00 800.00 850.00 900.00 950.00 1000.0 1327.0 1527.0 ;
  805.   LCP1     = prog 430.00 450.00 487.00 508.00 526.00 542.00 560.00 581.00 605.00 630.00 659.00 696.00 746.00 823.00 946.00 1510.0 873.00 619.00 633.00 611.00 615.00 735.00 835.00 ;
  806.   EVCP1    = evol manu 'T' LT1 'C' LCP1 ;
  807.   C_A508   = EVCP1 ;
  808.   TBIB1 . 'C' = C_A508 ;
  809.  
  810.  
  811.  
  812. *--- Module d'Young (GPa) :
  813.   LT1      = prog 20.000 100.00 200.00 300.00 400.00 500.00 600.00 700.00 800.00 900.00 1000.0 1100.0 1200.0 1300.0 1400.0 5000.0 ;
  814.   LYM1     = prog 212.00 207.20 201.30 194.70 187.10 178.30 167.90 156.20 142.40 126.60 108.40 86.700 61.800 32.500 2.2000 2.2000 ;
  815.   YM_A508 = evol manu 'T' LT1 'YOUN' LYM1 ;
  816.   comm /!\ conversion [GPa] -> [Pa] ;
  817.   TBIB1 . 'YOUN' = 1.e9 * YM_A508 ;
  818.  
  819.  
  820.  
  821. *--- Coefficient de Poisson :
  822.   TBIB1 . 'NU'   = 0.294 ;
  823.  
  824.  
  825.  
  826. *--- Coeficient de dilatation thermique [1/K] :
  827.   LT1      = prog 25.000 50.000 100.00 150.00 200.00 250.00 300.00 350.00 400.00 450.00 500.00 550.00 600.00 650.00 700.00 750.00 800.00 850.00 900.00 950.00 1000.0 1100.0 1200.0 1300.0 1400.0 ;
  828.   LALPHA1  = prog 7.83000E-06 7.82000E-06 7.81000E-06 7.80000E-06 7.78000E-06 7.77000E-06 7.75000E-06 7.73000E-06 7.72000E-06 7.70000E-06 7.68000E-06 7.67000E-06 7.65000E-06 7.70000E-06 7.68000E-06 7.66000E-06 7.64000E-06 7.62000E-06 7.59000E-06 7.57000E-06 7.55000E-06 7.55000E-06 7.55000E-06 7.55000E-06 7.55000E-06 ;
  829.   AL_A508 = evol manu 'T' LT1 'ALPH' LALPHA1 ;
  830.   TBIB1 . 'ALPH' = AL_A508 ;
  831.  
  832.  
  833.  
  834. *--- limite élastique [MPa]
  835.   LT1      = prog 20.000 316.00 538.00 760.00 982.00 1200.0 1400.0 ;
  836.   LSY1     = prog 529.00 471.30 395.70 62.400 30.800 16.400 2.0000 ;
  837.   EVSY1    = evol manu 'T' LT1 'SIGY' LSY1 ;
  838.   comm /!\ conversion [MPa] -> [Pa] ;
  839.   SY_A508 = 1.e6 * EVSY1 ;
  840.   TBIB1 . 'SIGY' = SY_A508 ;
  841.  
  842.  
  843.  
  844. *--- Coefficients ecrouissage cinematique modele EP Chaboche :
  845. *    LCiMur : coefficient Ci article Muransky
  846. *    LGiMur : coefficient Gammai article Muransky
  847. *    LCiC3M : coefficient Ci loi EP Chaboche Cast3M
  848. *    LAiC3M : coefficient Ai loi EP Chaboche Cast3M
  849.   LT1      = prog 20.000 316.00 538.00 760.00 982.00 1200.0 1400.0 ;
  850.   LC1Mur   = prog 1.36102E+05 2.55002E+05 5.08928E+05 2.15895E+05 14497. 0.0000  0.0000 ;
  851.   LG1Mur   = 5130.9 ;
  852.   LC2Mur   = prog 7308.9 10672. 5589.7 2053.8 1034.3 0.0000 0.0000 ;
  853.   LG2Mur   = 53.220 ;
  854.   comm /!\ conversion [MPa] -> [Pa] ;
  855.   LC1Mur   = 1.e6 * LC1Mur ;
  856.   LC2Mur   = 1.e6 * LC2Mur ;
  857.   LC1C3M   = LG1Mur ;
  858.   LC2C3M   = LG2Mur ;
  859.   LA1C3M   = LC1Mur / LG1Mur ;
  860.   LA2C3M   = LC2Mur / LG2Mur ;
  861.   EVA1     = evol manu 'T' LT1 'A1' LA1C3M ;
  862.   EVA2     = evol manu 'T' LT1 'A2' LA2C3M ;
  863.   R0_A508  = evol manu 'T' (extr SY_A508 absc) 'R0' (extr SY_A508 ordo) ;
  864.   TBIB1 . 'R0' = R0_A508 ;
  865.   TBIB1 . 'A1' = EVA1 ;
  866.   TBIB1 . 'C1' = LC1C3M ;
  867.   TBIB1 . 'A2' = EVA2 ;
  868.   TBIB1 . 'C2' = LC2C3M ;
  869.  
  870.  
  871.  
  872. *--- Coefficients ecrouissage isotrope modele EP Chaboche :
  873.   TBIB1 . 'RM' = 0. + (evol manu 'T' (extr SY_A508 absc) 'RM' (extr SY_A508 ordo)) ;
  874.   TBIB1 . 'B'  = 20. ;
  875.  
  876.  
  877.  
  878.     TBIB1 . 'CREATEUR' = mot 'BIBLIO' ;
  879.     TBIB1 . 'REFERENCE' = chai mot1 '_' ent1 ;
  880.  
  881.     resp TBIB1 ;
  882.     quit biblio ;
  883.  
  884.   fins ;
  885. *--- Fin [4], "A508"
  886.  
  887. *------------------------------ A508 [5] ------------------------------*
  888. * Reference : voir notice BIBLIO
  889.  
  890.   si (ent1 ega 5) ;
  891.  
  892. *--- Module d'Young (GPa) :
  893.   ltc1     = prog 20. 50. PAS 50. 900. 1000. 1050. 1250. 1500. ;
  894.   ltYm1    = prog 204. 203. 201. 199. 196. 192. 187. 182. 176. 171. 166. 158. 142. 105. 82. 64. 52. 42. 33. 15. 9.3 2.5 1.5 ;
  895.   comm /!\ conversion [GPa] -> [Pa] ;
  896.   ltYm1    = 1.e9 * ltYm1 ;
  897.   YM_A508  = evol manu 'T' ltc1 'YOUN' ltYm1 ;
  898.   TBIB1 . 'YOUN' = YM_A508 ;
  899.  
  900.  
  901.  
  902. *--- Coefficient de Poisson :
  903.   TBIB1 . 'NU'   = 0.3 ;
  904.  
  905.     TBIB1 . 'CREATEUR' = mot 'BIBLIO' ;
  906.     TBIB1 . 'REFERENCE' = chai mot1 '_' ent1 ;
  907.  
  908.     resp TBIB1 ;
  909.     quit biblio ;
  910.  
  911.   fins ;
  912. *--- Fin [5], "A508"
  913.  
  914. * Numero reference bib. ent1 n'existe pas pour mot1
  915.   erre (chai '*****  La reference' ' ' ent1 ' ne concerne pas' ' ' mot1) ;
  916.   quit biblio ;
  917.  
  918. * Fin mot-cle "A508"
  919. fins ;
  920.  
  921. *----------------------------------------------------------------------*
  922. *                          MOT-CLE : Austenite                         *
  923. *----------------------------------------------------------------------*
  924.  
  925. si (ega mot1 'Austenite') ;
  926.   imot1    = 5 ;
  927.  
  928. * Initialisation table de sortie :
  929.   TBIB1    = table ;
  930.  
  931. * Numero reference bib. par defaut :
  932.   si (ent1 ega 0) ; ent1 = 5 ; fins ;
  933.  
  934. *---------------------------- Austenite [5] ---------------------------*
  935. * Reference : voir notice BIBLIO
  936.  
  937.   si (ent1 ega 5) ;
  938.  
  939. *   Appel a BIBLIO pour proprietes elastiques :
  940.     t16MND5  = biblio 16MND5 refe 5 ;
  941.  
  942. *--- Module d'Young :
  943.     TBIB1 . 'YOUN' = t16MND5 . 'YOUN' ; 
  944.  
  945. *--- Coefficient de Poisson :
  946.     TBIB1 . 'NU'   = t16MND5 . 'NU' ; 
  947.  
  948. *--- Dilatation thermique :
  949.     TBIB1 . 'TALP' = 25. ;
  950.     ltc1     = prog 25. 1500. ;
  951.     lAlph1   = (2.52e-9 * ltc1) + 22.6e-6 ;
  952.     evCaAus1 = evol manu 'T' ltc1 'ALPH' lAlph1 ;
  953.     TBIB1 . 'ALPH' = evCaAus1 ;
  954.  
  955. *--- Difference de compacite avec la ferrite :
  956.     TBIB1 . 'DCOM' = -0.01 ;
  957.  
  958. *--- Limite d'elasticite :
  959.     ltc1     = prog 20. 200. 400. 500. 600. PAS 50. 850. 1000. 1100. 1250. 1500. ;
  960.     lSyAus1  = prog 250. 200. 150. 120. 95. 85. 75. 67. 60. 52. 30. 18. 5. 5. ;
  961.     lSyAus1  = 1.e6 * lSyAus1 ;
  962.     evCSyAu1 = evol manu 'T' ltc1 'SIGY' lSyAus1 ;
  963.     TBIB1 . 'SIGY' = evCSyAu1 ;
  964.  
  965. *--- Courbes d'ecrouissage :
  966.     lep1     = prog 0. 0.05 0.1 0.15 0.2 ;
  967.     T1       = 20. ;
  968.     lsm1     = prog 0. 345. 490. 570. 630. ;
  969.     lsm1     = 1.e6 * lsm1 ;
  970.     sigyT1   = ipol evCSyAu1 T1 ;
  971.     lsm1     = lsm1 + sigyT1 ;
  972.     evecAus1 = evol manu 'EPSP' lep1 'SIG' lsm1 ;
  973.     evecAus1 = evecAus1 coul 1 ;
  974.     nuagAus1 = nuag comp 'T' T1 comp 'ECRO' evecAus1 ;
  975.     T1       = 200. ;
  976.     lsm1     = prog 0. 345. 490. 570. 630. ;
  977.     lsm1     = 1.e6 * lsm1 ;
  978.     sigyT1   = ipol evCSyAu1 T1 ;
  979.     lsm1     = lsm1 + sigyT1 ;
  980.     evecAus1 = evol manu 'EPSP' lep1 'SIG' lsm1 ;
  981.     evecAus1 = evecAus1 coul 2 ;
  982.     nuagAus1 = nuagAus1 et (nuag comp 'T' T1 comp 'ECRO' evecAus1) ;
  983.     T1       = 500. ;
  984.     lsm1     = prog 0. 145. 250. 345. 440. ;
  985.     lsm1     = 1.e6 * lsm1 ;
  986.     sigyT1   = ipol evCSyAu1 T1 ;
  987.     lsm1     = lsm1 + sigyT1 ;
  988.     evecAus1 = evol manu 'EPSP' lep1 'SIG' lsm1 ;
  989.     evecAus1 = evecAus1 coul 3 ;
  990.     nuagAus1 = nuagAus1 et (nuag comp 'T' T1 comp 'ECRO' evecAus1) ;
  991.     T1       = 700. ;
  992.     lsm1     = prog 0. 120. 215. 310. 405. ;
  993.     lsm1     = 1.e6 * lsm1 ;
  994.     sigyT1   = ipol evCSyAu1 T1 ;
  995.     lsm1     = lsm1 + sigyT1 ;
  996.     evecAus1 = evol manu 'EPSP' lep1 'SIG' lsm1 ;
  997.     evecAus1 = evecAus1 coul 4 ;
  998.     nuagAus1 = nuagAus1 et (nuag comp 'T' T1 comp 'ECRO' evecAus1) ;
  999.     T1       = 800. ;
  1000.     lsm1     = prog 0. 95. 165. 230. 295. ;
  1001.     lsm1     = 1.e6 * lsm1 ;
  1002.     sigyT1   = ipol evCSyAu1 T1 ;
  1003.     lsm1     = lsm1 + sigyT1 ;
  1004.     evecAus1 = evol manu 'EPSP' lep1 'SIG' lsm1 ;
  1005.     evecAus1 = evecAus1 coul 5 ;
  1006.     nuagAus1 = nuagAus1 et (nuag comp 'T' T1 comp 'ECRO' evecAus1) ;
  1007.     T1       = 900. ;
  1008.     lsm1     = prog 0. 36. 66. 96. 126. ;
  1009.     lsm1     = 1.e6 * lsm1 ;
  1010.     sigyT1   = ipol evCSyAu1 T1 ;
  1011.     lsm1     = lsm1 + sigyT1 ;
  1012.     evecAus1 = evol manu 'EPSP' lep1 'SIG' lsm1 ;
  1013.     evecAus1 = evecAus1 coul 6 ;
  1014.     nuagAus1 = nuagAus1 et (nuag comp 'T' T1 comp 'ECRO' evecAus1) ;
  1015.     T1       = 1000. ;
  1016.     lsm1     = prog 0. 20. 40. 60. 80. ;
  1017.     lsm1     = 1.e6 * lsm1 ;
  1018.     sigyT1   = ipol evCSyAu1 T1 ;
  1019.     lsm1     = lsm1 + sigyT1 ;
  1020.     evecAus1 = evol manu 'EPSP' lep1 'SIG' lsm1 ;
  1021.     evecAus1 = evecAus1 coul 7 ;
  1022.     nuagAus1 = nuagAus1 et (nuag comp 'T' T1 comp 'ECRO' evecAus1) ;
  1023.     T1       = 1100. ;
  1024.     lsm1     = prog 0. 5. 10. 15. 20. ;
  1025.     lsm1     = 1.e6 * lsm1 ;
  1026.     sigyT1   = ipol evCSyAu1 T1 ;
  1027.     lsm1     = lsm1 + sigyT1 ;
  1028.     evecAus1 = evol manu 'EPSP' lep1 'SIG' lsm1 ;
  1029.     evecAus1 = evecAus1 coul 8 ;
  1030.     nuagAus1 = nuagAus1 et (nuag comp 'T' T1 comp 'ECRO' evecAus1) ;
  1031.     T1       = 1250. ;
  1032.     lsm1     = prog 0. 2.5 5. 7.5 10. ;
  1033.     lsm1     = 1.e6 * lsm1 ;
  1034.     sigyT1   = ipol evCSyAu1 T1 ;
  1035.     lsm1     = lsm1 + sigyT1 ;
  1036.     evecAus1 = evol manu 'EPSP' lep1 'SIG' lsm1 ;
  1037.     evecAus1 = evecAus1 coul 10 ;
  1038.     nuagAus1 = nuagAus1 et (nuag comp 'T' T1 comp 'ECRO' evecAus1) ;
  1039.     TBIB1 . 'ECRO' = nuagAus1 ;
  1040.  
  1041.  
  1042.     TBIB1 . 'CREATEUR' = mot 'BIBLIO' ;
  1043.     TBIB1 . 'REFERENCE' = chai mot1 '_' ent1 ;
  1044.  
  1045.     resp TBIB1 ;
  1046.     quit biblio ;
  1047.  
  1048.   fins ;
  1049. *--- Fin [5], "Austenite"
  1050.  
  1051. * Numero reference bib. ent1 n'existe pas pour mot1
  1052.   erre (chai '*****  La reference' ' ' ent1 ' ne concerne pas' ' ' mot1) ;
  1053.   quit biblio ;
  1054.  
  1055. * Fin mot-cle "Austenite"
  1056. fins ;
  1057.  
  1058. *----------------------------------------------------------------------*
  1059. *                          MOT-CLE : Ferrite                           *
  1060. *----------------------------------------------------------------------*
  1061.  
  1062. si (ega mot1 'Ferrite') ;
  1063.   imot1    = 5 ;
  1064.  
  1065. * Initialisation table de sortie :
  1066.   TBIB1    = table ;
  1067.  
  1068. * Numero reference bib. par defaut :
  1069.   si (ent1 ega 0) ; ent1 = 5 ; fins ;
  1070.  
  1071. *----------------------------- Ferrite [5] ----------------------------*
  1072. * Reference : voir notice BIBLIO
  1073.  
  1074.   si (ent1 ega 5) ;
  1075.  
  1076. *   Appel a BIBLIO pour proprietes elastiques :
  1077.     t16MND5  = biblio 16MND5 refe 5 ;
  1078.  
  1079. *--- Module d'Young :
  1080.     TBIB1 . 'YOUN' = t16MND5 . 'YOUN' ; 
  1081.  
  1082. *--- Coefficient de Poisson :
  1083.     TBIB1 . 'NU'   = t16MND5 . 'NU' ; 
  1084.  
  1085. *--- Dilatation thermique :
  1086.     TBIB1 . 'TALP' = 25. ;
  1087.     ltc1     = prog 25. 350. ;
  1088.     lalpFer1 = 7.7e-9 * ltc1 + 12.35e-6 ;
  1089.     ltc1     = ltc1     et 1500.  ;
  1090.     lalpFer1 = lalpFer1 et 15.e-6 ;
  1091.     evCaFer1 = evol manu 'T' ltc1 'ALPH' lalpFer1 ;
  1092.     TBIB1 . 'ALPH' = evCaFer1 ;
  1093.  
  1094. *--- Limite d'elasticite :
  1095.     ltc1     = prog 20. 200. 400. 500. 600. PAS 50. 850. 1000. 1100. 1250. 1500. ;
  1096.     lSyMB1   = prog 600. 530. 453. 398. 294. 208. 150. 105. 66. 50. 30. 18. 5. 5. ;
  1097.     lSyMB1   = 1.e6 * lSyMB1 ;
  1098.     evCSyMB1 = evol manu 'T' ltc1 'SIGY' lSyMB1 ;
  1099.     TBIB1 . 'SIGY' = evCSyMB1 ;
  1100.  
  1101. *--- Courbes d'ecrouissage :
  1102.     lep1     = prog 0. PAS 0.005 0.055 0.1 0.15 0.2 ;
  1103.     T1       = 20. ;
  1104.     lsm1     = prog 0. 52. 71. 85. 100. 111. 125. 138. 148. 159. 169. 177. 243. 284. 322. ;
  1105.     lsm1     = 1.e6 * lsm1 ;
  1106.     sigyT1   = ipol evCSyMB1 T1 ;
  1107.     lsm1     = lsm1 + sigyT1 ;
  1108.     evecFer1 = evol manu 'EPSP' lep1 'SIG' lsm1 ;
  1109.     evecFer1 = evecFer1 coul 1 ;
  1110.     nuagFer1 = nuag comp 'T' T1 comp 'ECRO' evecFer1 ;
  1111.     T1       = 100. ;
  1112.     lsm1     = prog 0. 22. 48. 68. 84. 96. 107. 119. 126. 133. 140. 148. 200. 235. 260. ;
  1113.     lsm1     = 1.e6 * lsm1 ;
  1114.     sigyT1   = ipol evCSyMB1 T1 ;
  1115.     lsm1     = lsm1 + sigyT1 ;
  1116.     evecFer1 = evol manu 'EPSP' lep1 'SIG' lsm1 ;
  1117.     evecFer1 = evecFer1 coul 2 ;
  1118.     nuagFer1 = nuagFer1 et (nuag comp 'T' T1 comp 'ECRO' evecFer1) ;
  1119.     T1       = 200. ;
  1120.     lsm1     = prog 0. 22. 41. 55. 70. 83. 95. 109. 121. 129. 138. 145. 195. 230. 255. ;
  1121.     lsm1     = 1.e6 * lsm1 ;
  1122.     sigyT1   = ipol evCSyMB1 T1 ;
  1123.     lsm1     = lsm1 + sigyT1 ;
  1124.     evecFer1 = evol manu 'EPSP' lep1 'SIG' lsm1 ;
  1125.     evecFer1 = evecFer1 coul 3 ;
  1126.     nuagFer1 = nuagFer1 et (nuag comp 'T' T1 comp 'ECRO' evecFer1) ;
  1127.     T1       = 300. ;
  1128.     lsm1     = prog 0. 10. 30. 50. 75. 95. 111. 125. 138. 148. 160. 169. 210. 230. 235. ;
  1129.     lsm1     = 1.e6 * lsm1 ;
  1130.     sigyT1   = ipol evCSyMB1 T1 ;
  1131.     lsm1     = lsm1 + sigyT1 ;
  1132.     evecFer1 = evol manu 'EPSP' lep1 'SIG' lsm1 ;
  1133.     evecFer1 = evecFer1 coul 4 ;
  1134.     nuagFer1 = nuagFer1 et (nuag comp 'T' T1 comp 'ECRO' evecFer1) ;
  1135.     T1       = 400. ;
  1136.     lsm1     = prog 0. 30. 57. 73. 83. 92. 99. 105. 112. 117. 123. 127. 156. 175. 183. ;
  1137.     lsm1     = 1.e6 * lsm1 ;
  1138.     sigyT1   = ipol evCSyMB1 T1 ;
  1139.     lsm1     = lsm1 + sigyT1 ;
  1140.     evecFer1 = evol manu 'EPSP' lep1 'SIG' lsm1 ;
  1141.     evecFer1 = evecFer1 coul 5 ;
  1142.     nuagFer1 = nuagFer1 et (nuag comp 'T' T1 comp 'ECRO' evecFer1) ;
  1143.     T1       = 500. ;
  1144.     lsm1     = prog 0. 33. 49. 60. 66. 71. 75. 79. 82. 84. 86. 87. 100. 115. 130. ;
  1145.     lsm1     = 1.e6 * lsm1 ;
  1146.     sigyT1   = ipol evCSyMB1 T1 ;
  1147.     lsm1     = lsm1 + sigyT1 ;
  1148.     evecFer1 = evol manu 'EPSP' lep1 'SIG' lsm1 ;
  1149.     evecFer1 = evecFer1 coul 6 ;
  1150.     nuagFer1 = nuagFer1 et (nuag comp 'T' T1 comp 'ECRO' evecFer1) ;
  1151.     T1       = 600. ;
  1152.     lsm1     = prog 0. 17. 24. 27. 30. 32. 34. 35. 36. 36.4 36.8 37.2 39. 41. 43. ;
  1153.     lsm1     = 1.e6 * lsm1 ;
  1154.     sigyT1   = ipol evCSyMB1 T1 ;
  1155.     lsm1     = lsm1 + sigyT1 ;
  1156.     evecFer1 = evol manu 'EPSP' lep1 'SIG' lsm1 ;
  1157.     evecFer1 = evecFer1 coul 7 ;
  1158.     nuagFer1 = nuagFer1 et (nuag comp 'T' T1 comp 'ECRO' evecFer1) ;
  1159.     T1       = 650. ;
  1160.     lsm1     = prog 0. 0.25 0.5 0.75 1. 1.25 1.5 1.75 2. 2.25 2.5 2.7 4.4 5.6 6.5 ;
  1161.     lsm1     = 1.e6 * lsm1 ;
  1162.     sigyT1   = ipol evCSyMB1 T1 ;
  1163.     lsm1     = lsm1 + sigyT1 ;
  1164.     evecFer1 = evol manu 'EPSP' lep1 'SIG' lsm1 ;
  1165.     evecFer1 = evecFer1 coul 8 ;
  1166.     nuagFer1 = nuagFer1 et (nuag comp 'T' T1 comp 'ECRO' evecFer1) ;
  1167.     T1       = 700. ;
  1168.     lsm1     = prog 0. 0.25 0.5 0.75 1. 1.25 1.5 1.75 2. 2.25 2.5 2.78 4.2 5.3 6.2 ;
  1169.     lsm1     = 1.e6 * lsm1 ;
  1170.     sigyT1   = ipol evCSyMB1 T1 ;
  1171.     lsm1     = lsm1 + sigyT1 ;
  1172.     evecFer1 = evol manu 'EPSP' lep1 'SIG' lsm1 ;
  1173.     evecFer1 = evecFer1 coul 10 ;
  1174.     nuagFer1 = nuagFer1 et (nuag comp 'T' T1 comp 'ECRO' evecFer1) ;
  1175.     T1       = 750. ;
  1176.     lsm1     = prog 15 * 1.e-3 ;
  1177.     lsm1     = 1.e6 * lsm1 ;
  1178.     sigyT1   = ipol evCSyMB1 T1 ;
  1179.     lsm1     = lsm1 + sigyT1 ;
  1180.     evecFer1 = evol manu 'EPSP' lep1 'SIG' lsm1 ;
  1181.     evecFer1 = evecFer1 coul 11 ;
  1182.     nuagFer1 = nuagFer1 et (nuag comp 'T' T1 comp 'ECRO' evecFer1) ;
  1183.     TBIB1 . 'ECRO' = nuagFer1 ;
  1184.  
  1185.  
  1186.     TBIB1 . 'CREATEUR' = mot 'BIBLIO' ;
  1187.     TBIB1 . 'REFERENCE' = chai mot1 '_' ent1 ;
  1188.  
  1189.     resp TBIB1 ;
  1190.     quit biblio ;
  1191.  
  1192.   fins ;
  1193. *--- Fin [5], "Ferrite"
  1194.  
  1195. * Numero reference bib. ent1 n'existe pas pour mot1
  1196.   erre (chai '*****  La reference' ' ' ent1 ' ne concerne pas' ' ' mot1) ;
  1197.   quit biblio ;
  1198.  
  1199. * Fin mot-cle "Ferrite"
  1200. fins ;
  1201.  
  1202. *------------------------- Mot-cle non trouve -------------------------*
  1203.  
  1204. si (imot1 ega 0) ;
  1205.   erreur (chai '*****  Pas de donnee disponible pour' ' ' mot1) ;
  1206.   quit biblio ;
  1207. fins ;
  1208.  
  1209. FINP ;
  1210.  
  1211.  
  1212.  
  1213.  

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