Télécharger metallurgie_07.dgibi

Retour à la liste

Numérotation des lignes :

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

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