Télécharger metallurgie_08.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : metallurgie_08.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5. *__________________________________________________________________*
  6. * TEST METALLURGIE_08 *
  7. * *
  8. * CALCUL DES PROPORTIONS DE PHASE METALLURGIQUE *
  9. * *
  10. * *
  11. * Le modele metallurgique est code en ESOPE (Cast3m) *
  12. * *
  13. *__________________________________________________________________*
  14. * *
  15. * Vitesse de refroidissement : 100 C/s *
  16. * - Transformation (type Koistinen-Marburger) de l'autenite *
  17. * en martensite seulement *
  18. * - La valeur calculee de la proportion d'AUSTENITE est *
  19. * comparée avec la valeur analytique *
  20. *__________________________________________________________________*
  21. *
  22. 'OPTI' 'TRAC' PSC 'EPTR' 10 ;
  23. IDESS = VRAI ;
  24.  
  25. OPTI DIME 2 ELEM 'TRI3';
  26.  
  27. P1 = 0. 0. ;
  28. P2 = 1. 0. ;
  29. P3 = 0. 1. ;
  30.  
  31. MAILT = MANU 'TRI3' P1 P2 P3;
  32.  
  33. LISTPHA = MOTS 'AUST' 'MART' ;
  34. LISTREAC = MOTS 'AUST' ;
  35. LISTPROD = MOTS 'MART' ;
  36. LISTTYPE = MOTS 'KOIS' ;
  37. NOMCONS = '16MND5';
  38.  
  39. MODCP1 = MODE MAILT 'METALLURGIE' 'PHASES' LISTPHA
  40. 'REACTIFS' LISTREAC
  41. 'PRODUITS' LISTPROD
  42. 'TYPE' LISTTYPE
  43. 'CONS' NOMCONS ;
  44.  
  45. TINI = 900. ;
  46. TFIN = 20. ;
  47.  
  48. MS1 = 380. ;
  49. KM1 = 0.0247 ;
  50.  
  51. MATCP1 = 'MATE' MODCP1 'MS1' MS1 'KM1' KM1 ;
  52.  
  53. *Declaration de la vitesse de refroidissement : en Degres/s
  54. DT_vit = -100. ;
  55.  
  56. *Calcul de la duree du refroidissement :
  57. DUREE = -1. * (TINI - TFIN) / DT_vit ;
  58.  
  59. *Declaration du Temps initial :
  60. TEMPINI = 1. ;
  61. *Calcul du Temps final :
  62. TEMPFIN = TEMPINI + DUREE ;
  63. *Calcul du pas de temps en fonction de la vitesse de refroidissement :
  64. DTPS = -10 / DT_vit ;
  65.  
  66. LTCAL ='PROG' TEMPINI 'PAS' DTPS TEMPFIN ;
  67. T = TINI + (DT_vit * (LTCAL - TEMPINI)) ;
  68.  
  69. ALPHA = 1. ;
  70. *Initialisation de CHE1 : proportions de phases initiales :
  71. PROP = 'PROG' ALPHA ;
  72. CTYP = 'CHAI' 'VARIABLES_INTERNES' ;
  73. CHE1 ='MANU' 'CHML' MODCP1 'AUST' ('EXTR' PROP 1)
  74. 'MART' 0.
  75. 'TYPE' CTYP 'NOEUD' ;
  76.  
  77. nbpas = (DIME LTCAL) - 1;
  78. REPE CTEMP nbpas;
  79.  
  80. * On recupere le temps debut de pas et le temps fin de pas :
  81. temp0 = EXTR LTCAL &CTEMP ;
  82. temp1 = EXTR LTCAL (&CTEMP + 1);
  83. CHETPS0 = 'MANU' 'CHML' MAILT 'TEMP' temp0 ;
  84. CHETPS1 = 'MANU' 'CHML' MAILT 'TEMP' temp1 ;
  85. * On calcule la temperature debut de pas et la temperature fin de pas :
  86. T0 = EXTR T &CTEMP ;
  87. T1 = EXTR T (&CTEMP + 1) ;
  88. CHT0 = 'MANU' 'CHML' MAILT 'T' T0 ;
  89. CHT1 = 'MANU' 'CHML' MAILT 'T' T1 ;
  90.  
  91. CHE1 = COMP MODCP1 (CHE1 'ET' CHETPS0 'ET' CHT0 'ET' MATCP1 ) ( CHETPS1 'ET' CHT1 'ET' MATCP1 ) ;
  92.  
  93. * Extraction pour tracer :
  94. Pi ='EXCO' CHE1 AUST 'AUST';
  95. PiEvol = 'EXTR' Pi AUST 1 1 1;
  96. PROP = PROP 'ET' PiEvol ;
  97.  
  98. FIN CTEMP;
  99.  
  100. *Evolution calculee :
  101. EVO = 'EVOL' 'NOIR' 'MANU' 'TEMP' LTCAL 'AUST' PROP;
  102.  
  103. *Evolution analytique :
  104. TBORN = BORN T 'MAXI' MS1;
  105. PANALY = EXP (-1. * KM1 * (MS1 - TBORN)) ;
  106. PANALY = BORN PANALY 'MAXI' ALPHA ;
  107. EVOANA= 'EVOL' 'VERT' 'MANU' 'TEMP' LTCAL 'AUSTENITE' PANALY ;
  108.  
  109. *Erreur absolue entre analytique et calculee :
  110. ERRABS = 'ABS' (PANALY '-' PROP) ;
  111. EVOERR = 'EVOL' 'ROUG' 'MANU' 'TEMP' LTCAL 'ERREUR' ERRABS ;
  112.  
  113. 'SI' IDESS ;
  114.  
  115. TL = TABL ;
  116. TL . 1 = 'TIRR' ;
  117. TL . 2 = 'MARQ ROND NOLI' ;
  118. TL . 'TITRE' = TABL ;
  119. TL . 'TITRE' . 1 = 'Evolution analytique' ;
  120. TL . 'TITRE' . 2 = 'Evolution calculee' ;
  121. DESS ( EVOANA ET EVO ) 'GRIL' 'TIRC' 'GRIS'
  122. 'LEGE' TL
  123. 'TITR' 'KOISTINEN-MARBURGER - Refroidissement acier 16MND5 -250 deg/s' ;
  124.  
  125. TL = TABL ;
  126. TL . 1 = 'TIRR' ;
  127. TL . 'TITRE' = TABL ;
  128. TL . 'TITRE' . 1 = 'Erreur';
  129. DESS EVOERR 'GRIL' 'TIRC' 'GRIS'
  130. 'LEGE' TL
  131. 'TITR' 'KOISTINEN-MARBURGER - Refroidissement acier 16MND5 -250 deg/s' ;
  132.  
  133. 'FINSI' ;
  134.  
  135. *Emission d'une erreur GIBIANE si l'erreur analytique - cast3m est trop grande :
  136. Diff1 = 'MAXI' ERRABS;
  137. Tol1 = 2.D-3; COMM 'Valeur fixee en 2018 par Thibault LINDECKER';
  138.  
  139. MESS (CHAI ' Erreur : ' Diff1 ' | Tolerence : ' Tol1);
  140. SI (Diff1 >EG Tol1);
  141. MESS 'Erreur dans le calcul de la proportion d AUSTENITE';
  142. ERRE 5;
  143. FINS;
  144.  
  145. MESS Diff1 ;
  146.  
  147. fin;
  148.  
  149.  
  150.  

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