Télécharger la plaquette Cast3M


fr en

Télécharger dil2sat.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * dilatométrie et essai satoh (contrainte déplacement nul) éprouvette axisymmétrique
  2. * calcul /essai
  3. graph = 'N' ;
  4. opti dime 2 mode axis elem qua4 ;
  5.  
  6. p_ori = 0. 0. ; e_r = 1. 0. ; e_z = 0. 1. ;
  7.  
  8. ray1 = 5.e-3 ; h1 = 1.e-2 ;
  9. px_1 = p_ori plus p_ori;
  10. pe_1 = px_1 plus (ray1*e_r) ; l_1 = d px_1 1 pe_1 ;
  11. s_1 = l_1 trans 4 (e_z * h1) ; l_3 = s_1 cote 3 ;
  12. px_2 = l_3 point final ;
  13.  
  14. ************************************************************************
  15. * chargement temperature
  16. ch_T = manu chpo s_1 T 1. ;
  17. ev_T = evol manu temp (prog 0. pas 10. 200.) T (prog 0. pas 141. 1410. pas -141. 0.) ;
  18. cha_1 = char 'T' ev_T ch_T ;
  19. ************************************************************************
  20. * acier inox
  21. * coef dilatation
  22. l_temp1 = prog
  23.       20. 50. 150. 250. 350. 450. 550. 650. 750. 850. 875. ;
  24. al_in = evol manu t l_temp1  alph (1.e-6*(prog 
  25.  13.75 15.415 15.835 16.565 17.0 17.43 17.735 17.775 17.775
  26.    17.775 17.775)) ;
  27. a2i304l = evol manu t l_temp1 alph (1.e-6*(prog
  28.   17.775 17.775 18.33 18.53 18.884 18.914 18.967 18.695 17.999 
  29.    17.775 17.775 )) ;
  30. * module d young
  31. l_temp2 = prog 20. 200. 400. 600. 900. ;
  32. evly_in =  evol manu t l_temp2 youn (1.e9*(prog
  33.   192.6 207. 197.9 197.9 81.9)) ;
  34. y2i304l = evol manu t (prog 20. 200. 400. 600. 800. 900.)
  35.      youn (1.e9*(prog 198.85 198.85 198.85 198.85 81.9 68.)) ;
  36. ltmpt8 = prog 0. 100. 326. 604. 
  37.               800. 1000. 1500. 2000. ;
  38. lyoun2 = (prog 195. 195. 170. 145. 
  39. *             105. 50. 1. 1.) * 1.e9 ;
  40.              105. 5. .001 .0001) * 1.e9 ;
  41. lyoun1 = (prog 215. 215. 215. 180.
  42. *             105. 50. 1. 1.) * 1.e9 ;
  43.              105. 50. 0.001 .0001) * 1.e9 ;
  44. evly1 = evol manu 'T' ltmpt8 'YOUN' lyoun1 ;                          
  45. evly2 = evol manu 'T' ltmpt8 'YOUN' lyoun2 ; 
  46. evly_in = evly1 et evly2 ;
  47.  
  48. * seuil elasticite a 0,2%
  49. ltmpt5 = prog 0. 28. 326. 604.      
  50.          800.  1000. 1100. 1200. 1400. 1450. ;
  51. lsigy2 = (prog 220. 220. 162. 140.
  52.           106.  47. 26. 10. 1. 1.)*1.e6 ;
  53. evsigy2 = evol manu 'T' ltmpt5 'SIGY' lsigy2 coul jaune ; 
  54. lsigy1 = (prog 300. 300. 300.  250.
  55.             80.  25. 5. 1. 1. 1. )*1.e6 ;                   
  56. evsigy1 = evol manu 'T' ltmpt5 'SIGY' lsigy1 coul jaune ;  
  57. evsy_in = evsigy1 et evsigy2 ;                     
  58.  
  59.  
  60. l_h2 = (prog (0.01*100.) (0.01*100.) (0.01*100.) (0.01*100.) (0.02*100.) 
  61. *                (0.02*50) (0.01*10.) (0.01*1.))*1.e9 ;
  62.                0. 0. 0. ) * 1.e9 ;
  63.  
  64. l_h1 = (prog (0.08*195.) (0.02*215.) (0.02*215.) (0.01*215.) 
  65. *           (0.01*50.) (0.01*50.) (0.01*10.) (0.01*1.))*1.e9  ;            
  66.               0. 0. 0. 0. ) * 1.e9 ;
  67. ltmpt8 = prog 0. 100. 326. 604. 
  68.               800. 1000. 1500. 2000. ;
  69. evlh1 = evol manu 'T' ltmpt8 'H' l_h1 ;
  70. evlh2 = evol manu 'T' ltmpt8 'H' l_h2 ;
  71. evlh_in = evlh1 et evlh2 ;
  72.  
  73. mo_in = mode s_1 mecanique elastique plastique cinematique fusion cons 'avec';
  74. ca_in = mate mo_in youn evly_in nu 0.3 rho 7.8e9 alph al_in sigy evsy_in h evlh_in talp 20. tref 20.
  75.    tfus 1380.  xtma 20.;
  76. ************************************************************************
  77. * la courbe de dilatométrie de l'inox est linéaire réversible
  78. ************************************************************************
  79. * SATOH
  80. t_sin = table 'PASAPAS'                                          ;
  81. t_sin . modele = mo_in                                  ;
  82. t_sin . caracteristiques = ca_in                        ;
  83. t_sin . blocages_mecaniques = (bloque (l_1 et l_3) uz) et (bloq (px_1 et px_2) ur) ;
  84. t_sin . chargement = cha_1                               ;
  85. t_sin . temps_calcules = prog 0. pas 10. 200.            ;
  86.  
  87. pasapas t_sin ;
  88.  
  89.   ind = 0 ; ord_1 = prog ; ord_2 = prog ; abs_T = prog ;
  90.   repeter b_z (dime (t_sin.temps)) ;
  91.     abs_T = abs_T et (prog (maxi (tire cha_1 t_sin . temps . ind))) ;
  92.    ord_2 = ord_2 et (prog (extr t_sin . contraintes . ind smzz 1 2 2)) ;
  93.     ind = ind + 1 ;
  94.   fin b_z ;
  95.  
  96.   ev_fm_in = evol vert manu abs_T (ord_2 / 1.e6) ;
  97.  
  98. * inox : satoh 1200
  99. ray_0 = 5.e-3 ; coe_0 = 1./(pi*ray_0*ray_0*1.e6) ;
  100. *
  101. m_8b_f = prog 253.18878 -10261.238 -16346.998 -20703.32  
  102.   -24870.158 -26649.551 -26923.436 -26758.072 -25262.9 -17800.826
  103.   -12700.364 -8132.1699 -5364.0371 -3725.8948 -2652.74 -1546.8728
  104.   -354.8703 1863.7703 2716.4312 5029.812 6321.7222 7448.2676 
  105.   8263.032 9336.1787 10262.91 10932.98 11894.161 12621.076 
  106.   13301.482 14541.717 15465.001 16600.158 17699.145 17842.115 ;
  107.  
  108. *49.01936 H 29.
  109. m_8b_t = prog 22.557192 73.672592 155.48824 243.35089 
  110.   327.8024 399.28125 483.42267 560.276 639.76 736.10 
  111.    807.17072 891.67395 969.76831 1031.5824 1086.8325 1130.3503
  112.   1221.8308 1151.0756 1047.6044 901.75232 841.85065 794.7149 
  113.   748.76788 675.01495 586.89386 516.03522 427.29395 362.12051 
  114.   294.77637 200.65997 151.92204 100.18644 60.234772 29. ;
  115. *
  116. * pour avoir un bon debut de courbe
  117. m_8b_f2 = prog -118.88135 -1161.0223 -19997.074 -24756.471 
  118.   -26473.852 -27464.314 -26923.436 -26758.072 -25262.9 -17800.826
  119.   -12700.364 -8132.1699 -5364.0371 -3725.8948 -2652.74 -1546.8728
  120.   -354.8703 1863.7703 2716.4312 5029.812 6321.7222 7448.2676 
  121.   8263.032 9336.1787 10262.91 10932.98 11894.161 12621.076 
  122.   13301.482 14541.717 15465.001 16600.158 17699.145 17842.115 ;
  123.  
  124. m_8b_t2 = prog 20. 40.  102.35716 186.86037 
  125.   263.81772 349.81973 483.42267 560.276 639.76 736.10 
  126.    807.17072 891.67395 969.76831 1031.5824 1086.8325 1130.3503
  127.   1221.8308 1151.0756 1047.6044 901.75232 841.85065 794.7149 
  128.   748.76788 675.01495 586.89386 516.03522 427.29395 362.12051 
  129.   294.77637 200.65997 151.92204 100.18644 60.234772 29.01936 ;
  130.  
  131. es_8b = evol manu 'T (C)' m_8b_t2 'sm (MPa)' (m_8b_f2 * coe_0)
  132.      coul rose ;
  133.  
  134. si (ega graph 'OUI') ;
  135.  titre ' calcul SATOH inox avec TFUS et XTMA (vert)/ essai (rose)' ;
  136.   dess (ev_fm_in et es_8b) ;
  137. finsi ;
  138.     esi4 = (extr m_8b_f2 4)*coe_0;
  139.   err5 =  abs ( (((extr ord_2 2)/1.e6) - esi4) / esi4) ; 
  140.  
  141.     esi33 = (extr m_8b_f2 33)*coe_0;
  142.   err6 =  abs ( (((extr ord_2 20)/1.e6) - esi33) / esi33) ; 
  143.  
  144.   si ((err5 > 5.e-2) ou (err6 > 5.e-2)) ; erre 5 ; sinon  ; erre 0 ; finsi ;
  145.  
  146. ************************************************************************
  147. * acier ferritique martensitique
  148. * coef de dilatation
  149. l_temp1 = prog
  150.       20. 50. 150. 250. 350. 450. 550. 650. 750. 850. 875. ;
  151. a_16mnd5 = evol manu t l_temp1 alph (1.e-6*(prog 
  152.  9.694 11.01 11.597 12.087 12.493 12.677 12.783 12.783 12.783
  153.    12.783 12.783 ));
  154. a216mnd5 = evol manu t l_temp1 alph (1.e-6*(prog 
  155.  12.88 12.88 13.32 13.48 13.62 13.63 13.47 13.42 13.42 13.42 13.42)) ;
  156. a2_fm = evol manu 
  157.   'T' (prog 20. 250. 320. 500. 600. 800. 1200.)
  158.   'ALPH' (prog 1. 1. 0.7 0. 0.35 0.7 1.1 )*1.e-5 ;  
  159. aj16mnd5 = evol manu 'T' (prog 20. 1200.) 
  160.                    'ALPH' (prog 0.8 0.78)*1.e-5 ;
  161. a1_fm = evol manu 
  162.   'T' (prog 20. 250. 320. 750. 830. 1200.)
  163.   'ALPH' (prog 1. 1.1 1.2 1.2 .8 1.1)*1.e-5 ;                    
  164.  
  165. al_fm = a1_fm et a2_fm ;
  166.  
  167. * module d young
  168. l_temp2 = prog 20. 200. 400. 600. 700. 900. ;
  169. e1yo_fm = evol manu t l_temp2 youn (1.e9*(prog
  170.   183.4 198.3 140.4 158.7  100. 69.)) ;
  171. e2yo_fm = evol manu t (prog 20. 200. 300. 400. 500. 600. 700.
  172.    800. 900.) youn (1.e9*(prog 210.85 210.85 150.35 100.45 
  173.    90.55 77.7 70.85 60.0 33.0)) ;
  174.  
  175. evly_fm = e1yo_fm et e2yo_fm ;
  176. *
  177. * seuil elasticite a 02%
  178. e1sy_fm = evol manu t l_temp2 sigy (1.e6*(prog
  179.   417. 418. 402. 313. 50. 22.)) ;
  180. e2sy_fm = evol manu t (prog 20. 200. 400. 600. 800. 900. 1000.)
  181.    sigy (1.e6*(prog 400. 100. 100. 100. 50. 30. 2.)) ;
  182.  
  183. evsy_fm = e1sy_fm et e2sy_fm ;
  184. *
  185.  
  186.  
  187. mo_fm = mode s_1 mecanique elastique plastique parfait fusion cons 'avec';
  188. ca_fm = mate mo_fm youn evly_fm nu 0.3 rho 7.8e9 alph al_fm sigy evsy_fm talp 20. tref 20.
  189.    tfus 1380.  xtma 20.;
  190. ************************************************************************
  191. * dilatométrie 
  192. t_dil = table 'PASAPAS'                                  ;
  193. t_dil . modele = mo_fm                                   ;
  194. t_dil . caracteristiques = ca_fm                         ;
  195. t_dil . blocages_mecaniques = (bloque l_1 uz) et (bloq (px_1 et px_2) ur) ;
  196. t_dil . chargement = cha_1                               ;
  197. t_dil . temps_calcules = prog 0. pas 5. 200.            ;
  198.  
  199. pasapas t_dil ;
  200.  
  201.   ind = 0 ; ord_1 = prog ; ord_2 = prog ; abs_T = prog ;
  202.   repeter b_z (dime (t_dil.temps)) ;
  203.  
  204.     abs_T = abs_T et (prog (maxi (tire cha_1 t_dil . temps . ind))) ;
  205.     ord_2 = ord_2 et (prog (maxi (exco t_dil. deplacements . ind uz))) ;
  206.     ind = ind + 1 ;
  207.   fin b_z ;
  208.  
  209.   ev_fm_di = evol vert manu abs_T (ord_2/h1*1.e2)  ;
  210.  
  211. di_c37_2 = evol rose manu
  212.   'T' (prog 20. 250. 710. 750. 830. 1200. 500. 320. 250. 20.) 
  213.  'eps(%)' (prog 0. 0.25 0.85 0.9 0.64 1.3 0. 0.27 0.25 0.) ;  
  214.  
  215.  
  216. si (ega graph 'OUI') ;
  217.  titre ' calcul dilato ferritique martensite avec TFUS et XTMA (vert)/ essai (rose)' ;
  218.   dess (ev_fm_di et di_c37_2 ) ;
  219. finsi ;
  220.  
  221.   ez3 = (extr (extr di_c37_2 ordo) 3) ;
  222.   err1 =  abs ( (((extr ord_2 11)/h1*1.e2) - ez3) / ez3) ; 
  223.  
  224.   ez7 = (extr (extr di_c37_2 ordo) 7) ; 
  225.   err2 =  abs (((extr ord_2 34)/h1*1.e2) - ez7); 
  226.  
  227.   si ((err1 > 5.e-2) ou (err2 > 5.e-2)) ; erre 5; sinon ; erre 0 ; finsi ;
  228.  
  229. ************************************************************************
  230. * SATOH
  231. t_sat = table 'PASAPAS'                                          ;
  232. t_sat . modele = mo_fm                                  ;
  233. t_sat . caracteristiques = ca_fm                        ;
  234. t_sat . blocages_mecaniques = (bloque (l_1 et l_3) uz) et (bloq (px_1 et px_2) ur) ;
  235. t_sat . chargement = cha_1                               ;
  236. t_sat . temps_calcules = prog 0. pas 5. 200.            ;
  237.  
  238. pasapas t_sat ;
  239.  
  240.   ind = 0 ; ord_1 = prog ; ord_2 = prog ; abs_T = prog ;
  241.   repeter b_z (dime (t_sat.temps)) ;
  242.     abs_T = abs_T et (prog (maxi (tire cha_1 t_sat . temps . ind))) ;
  243.    ord_2 = ord_2 et (prog (extr t_sat . contraintes . ind smzz 1 2 2)) ;
  244.     ind = ind + 1 ;
  245.   fin b_z ;
  246.  
  247.   ev_fm_sa = evol vert manu abs_T (ord_2 / 1.e6) ;
  248.  
  249. ab2_t = prog  21.6145   79.515  167.503  249.236  332.933 
  250.    418.078  587.385  666.74    732.447  810.25  907.544  
  251.        992.896  1073.65  1159.31  1045.01  933.651  840.235
  252.   770.754  652.782  567.585  508.547  487.197  470.395  
  253.         432.191  409.6  387.111  345.13  301.553  246.703  205.862
  254.   151.891   103.14  49.5825  21.7179 ;
  255.  
  256. ord_s = prog  -9.06295 -10671. -22849.1  -32222.2  -32384.1
  257.   -31486.4  -24048.3   -12564.5  -3458.5  183.912  -2969.17
  258.           -2076.66  -1325.43  35.7349  1855.22  2575.43  3374.9
  259.   4315.65  5945.61  6605.51  4284.64  1438.25  -5.61696  
  260.        -4183.88  -3825.5  -5097.06  -8029.6  -5720.79  -0.447977  5700.94
  261.   13282.1  20215.4  27377.9  33623.8  ;
  262. *
  263. ev_sa = evol rose manu 'T' ab2_t 'sig(MPa)' (ord_s/7.843e-5/1.e6) ;
  264. *dess ev_sa titre ' essai Satoh / 16MND5 / yc38_2 ' ;
  265. si (ega graph 'OUI') ; 
  266.  titre ' calcul SATOH ferritique martensite avec TFUS et XTMA (vert)/ essai (rose)' ;
  267.   dess (ev_fm_sa et ev_sa) ;
  268. finsi ;
  269.  
  270.     es5 = (extr ord_s 5)/7.843e-5/1.e6 ;
  271.   err3 =  abs ( (((extr ord_2 4)/1.e6) - es5) / es5) ; 
  272.  
  273.     es27 = (extr ord_s 27)/7.843e-5/1.e6 ;
  274.   err4 =  abs ( (((extr ord_2 37)/1.e6) - es27) / es27) ; 
  275.  
  276.   si ((err3 > 5.e-2) ou (err4 > 5.e-2)) ; erre 5 ; sinon  ; erre 0 ; finsi ;
  277.  
  278. ************************************************************************
  279. fin ;
  280.  
  281.  
  282.  
  283.  
  284.  
  285.  

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