===== Appendix 12: *.MET file ===== File number = 35 ; Generally called *IN.MET\\ This *.met file is read in the Prepro by METLAW.F\\ \\ This file contains all the data necessary to use effectively the laws [[laws:meta|META]], [[laws:metamec|METAMEC]], [[laws:elamet|ELAMET]], [[laws:arbthmet|ARBTHMET]] and [[laws:thmet|THMET]]. It must always exist to perform a metallurgical thermal analysis. Sections 1 to 8 are repeated with increasing ILAWN if more than one steel is described. ==== 1. Title ==== ^Title (A70)^^ |Any comment that will be reproduced on the output listing. Try to characterise your steel (60NCD11, ARBED, 42CD4, ...)|| ==== 2. General data ==== ^General data (10I5, G10.0, 2I5)^^ |ILAWN|Number of the steel described. This number is entered under the reference number IMETA by the law [[laws:meta|META]]| |IMPER|= 0 No impression\\ = 1 Impression on file number 36 generally called *IN.OUM| |NTPCA|Number of parameters in section [[appendices:a12# Characteristic temperatures (TTT + equilibrium F-C diagrams) |3.]] \\ = 8 for [[laws:meta|META]]; \\ = 20 for [[laws:metamec|METAMEC]]| |NPA|= 5; Number of parameters described by polynomials (section [[appendices:a12#2. Parameters described by polynomials of temperature |4.]])| |NDPO|Maximum degree of polynomials (maximum value = 7)| |NVM|= 0 No mechanical parameters described| |NT1| Maximum number of proeutectoid| |NT2| temperatures in pearlite| |NT3| the data tables related to bainite| |NTEMP|No mechanical parameters depending on the temperature| |DT|Temperature used during the simulation = temperature given in the .MET file + DT (a non-null value can be used if the temperature values in the .DAT file are expressed sing an unity that is different from the temperatures in the .MET file, for instance celcius in one file and Kelvin in the other)| |IPOLY| 1 Thermo-physical parameters $\lambda, \rho, C, H_v$ and L are given as polynomials function of the temperature. (see [[appendices:a12#Parameters described by polynomials of temperature |Parameters described by polynomials of temperature]]) \\ 0 Thermo-physical parameters $\lambda, \rho, C, H_v$ and L are given as data tables, functions of the temperature (see [[appendices:a12#Parameters described by polynomials of temperature |Parameters described by polynomials of temperature]]) \\ The preprocessor displays explicit information on-screen about this parameter.| |IET|= 1: Definition of the tangent modulus according to the strain level for each phase and temperature| ^If IET = 1 (I5)^^ |NEPS| Number of strain levels| __Remark:__ Some values are already defined in [[laws:meta# Integer parameters|Integer parameters]], pay attention to give the same value. Values as NTPCA, NPA, NVM, NTEMP are defined a first time automatically in LAWPRE. If you want to change these values, change also in the FORTRAN source file called LMETA.F. ==== 3. Characteristic temperatures (TTT + equilibrium F-C diagrams)==== ^Title (A5)^^ |Title|TPCAR written in columns 1 to 5| ^Parameters (7G10.0/7G10.0/6G10.0) - Only NTPCA parameters are read^^ |$A_3$ or $A_{cm}$|$A_3$: equilibrium temperature for the beginning of the ferrite transformation\\ $A_{cm}$:equilibrium temperature for the beginning of the cementite transformation| |$A_1$|equilibrium temperature for the eutectoïd transformation| |TH|Under the temperature TH, the pearlitic transformation is not preceded by the proeutectoïd transformation| |$B_s$|Temperature of the possible beginning of the bainitic transformation| |$B_f$|Under this temperature the bainitic transformation is complete| |$M_s$|Beginning temperature for the martensite transformation| |AM|Coefficient of the Marburger law for the martensite transformation| |FINCU|If no transformation has occurred when the temperature $B_s$ is reached, the SCHEIL's sum is multiplied by FINCU (generally FINCU = 0.0)| |CP_e| Values defining the shift in the diagram TTT : $D = C\ \sigma_{equivalent}$ for the ferrite and the pearlite | |CB_a| Values defining the shift in the diagram TTT : $D = C\ \sigma_{equivalent}$ for the bainite | |A| Values that gives the variation of M_s | |B| $\Delta M_s = A \sigma_{moi} + B \sigma_{equivalent}$ | |EXPR| $\gamma \rightarrow $ Pr | |EXPE| $\gamma \rightarrow $ Pe : Dilatation due to the austenite transformation | |EXBA| $\gamma \rightarrow $ Ba (the reference volume is the austenite at 0E C) | |EXMA| $\gamma \rightarrow $ MA | |K4=K3| Coefficient in the plasticity transformation formulae : ferrite, cementite, pearlite | |K5| Coefficient in the plasticity transformation formulae : bainite | |K6| Coefficient in the plasticity transformation formulae : martensite | |TLIQUID| Temperature where the steel is considered to be fully liquid. Beyond this temperature, the preprocessor will automatically set the thermal dilatation coefficient to null values.\\ \textbf{Important : put an initial value even if you don't model liquid state}| All the characteristic temperatures are defined on the figure below: {{ :laws:meta.png?500 |}} __Remark__ : Some additional parameters can occur depending on the steel and its plasticity transformation formula or the modification of the formula of the shift ($D=C\sigma$). If you want to change, you must adopt: - NTPCA (section 2); - Subroutine METLAW that read and write with comments the parameters of section 3; - Subroutine ARMEA that read the great vector PAMET where are stored the parameters of section 3 and where the formulae of $D$ and $\varepsilon^{pt}$ are implemented. ==== 3bis. If IET = 1 ==== ^Title (A5)^^ |TITLE|STLVL in the first 5 columns| ^Repeated NEPS times (G10.0)^^ |EPS|Values of the NEPS strain levels (variable tangent modulus)| ====4. Parameters described by polynomials of temperature==== ^Title (A5)^^ |Title|POCOE written from columns 1 to 5| **If IPOLY =1** ^Parameter definition (A5)^^ |A1∴A2|See explanation below| ^Polynomial coefficients (NDPO G10.0)^^ |A(I) \\ I=1,NDPO+1|See explanation below| ^End of section^^ |The end of this section is detected by writing "FI" followed by a blank card|| The conductivity $\lambda$, the mass density $\rho$, the heat capacity C and the hardness$ H_v$ have to be defined for each phase. The latent heat L is defined for each transformation. You can choose the order in which you want to define these parameters: - Firstly you give five letters (A5): 'A1^A2' where:\\ - A1 defines the thermal parameter\\ * $\lambda \rightarrow$ LA * C$\rightarrow$ CA * $\rho \rightarrow$ RO * $H_v \rightarrow$ HV * L $\rightarrow$ TR\\ - $A_2$ defines the phase; For $\lambda, \rho, C, H_v$: * Austenite $\rightarrow$ AU * Bainite$\rightarrow$ BA * Proeutectoïd $\rightarrow$ PR * Martensite$\rightarrow$MA * Pearlite$\rightarrow$PE \\ For L (latent heat) A2 = PR, PE, BA or MA to define in which phase the austenite is transformed. - Secondly you give the polynomial coefficient (A(I), I=1,NDPO+1) which defines the following polynomial: $$A(1)+A(2)T+A(3)T^{2}+ … +A(NDPO+1)T^{NDPO}$$ - Finally the non-defined parameters will be initialized to zero. So if you know that only certain phases will be present you do not need to define the other phase parameters. ** If IPOLY = 0 ** \\ One must write FI followed by a blank card.\\ :!: If input parameters (temperature dependent) are given as a table, all the tables must have the same length. Otherwise, it does not work properly. ==== 5. Mechanical parameters (Only for laws coupled with mechanics) ==== - Firstly you give five letters : 'A1^A2' and NTER where : * A1 defines the mechanical parameter: * YOUNG modulus: YO * POISSON ratio: NU * Thermal dilatation: AC or AP (AC for the $\alpha$ coefficient of classical type and AP for the $\alpha$ coefficient of partial type) * Yield stress $\sigma_y$: SY * Plastic slope: ET * If IPOLY=0: * Thermal conductivity: LA * Mass density: RO * Heat capacity: CA * Vickers hardness: HV * Latent heat of transformation: TR * '^' is a space; * A2 defines the phase concerned; * Austenite: AU (Except for A1=TR) * Proeutectoïd: PR * Pearlite: PE * Bainite: BA * Martensite: MA * NTER as the number of temperature used to describe the evolution of the parameter.\\ - Secondly, **ONLY IF** AC is chosen: \\ (A5,G10.0) 3 spaces and 'TO' or 'T0' \\ VALUE: Value of $T_0$ (usually $T_0$ is the room temperature so 20°C or 293K, be careful there is no default value for this parameter, so you should enter a value, otherwise the preprocessor will crash) - Thirdly, you repeat NTER times: * TEMPE: Temperature * VALUE: Value of the parameter __Remarks :__ - No defined tables are initialized to zero - A1 = FI followed by a blank card indicates the end of section 5. - For the table describing ALPHA, the first temperature __must be__ zero otherwise the integration of $\int_0^{T_{\alpha}} dT$ will not be correct. \\ __Remarks about the thermal coefficient $\alpha$:__ - The $\alpha_C$ coefficient of classical type is defined by : \[\alpha_C(T) = \frac{1}{L_0}\frac{L_{(T)}-L_0}{T-T_0}\] where $L_0=L(T_0)$.\\ The $\alpha_P$ coefficient of partial type is defined by : \[\alpha_P(T)=\frac{1}{L_{(T)}}\frac{dL}{dT}\] The user could give the classical type or the partial type $\alpha$ coefficient (with AC or AP).\\ - If $\alpha$ is of classical type ($\alpha_C$) : * The unity, chosen by the user, of $T_0$ has to be the same as the unity of the temperatures at which the $\alpha_C$ coefficient is given. * The temperatures, at which the $\alpha_C$ coefficient is given, have to be given in increasing order. * It is necessary to give at least 2 values of $\alpha_C$ at two temperatures. * Each $\alpha_C$ coefficient has to verify the following relation : \[1+\alpha_C(T).(T-T_0)>0\;\text{ and not equal to 0}\] - If the user gives the classical type $\alpha_C$ coefficient (AC), the pre-processor will calculate the partial type $\alpha_P$ coefficient. If the user gives the partial one, the pre-processor keeps these values. In the case of the calculus of the partial type $\alpha_P$ coefficient, the following relations are used (see the file ALPHAPARTIEL.F) for given couples ($T_i,\alpha_{Ci}$) with $i=1,...,n$ : \[\alpha_{P1} = \frac{\alpha_{C1}+(T_1-T_0)\frac{\alpha_{C2}-\alpha_{C1}}{T_2-T_1}}{1+\alpha_{C1}(T_1-T_0)}\]\[\alpha_{Pn} = \frac{\alpha_{Cn}+(T_n-T_0)\frac{\alpha_{Cn}-\alpha_{Cn-1}}{T_n-T_{n-1}}}{1+\alpha_{Cn}(T_n-T_0)}\] and for $i$ such as $1