[MAP Logo]

Materials Algorithms Project
Program Library


Program MAP_STEEL_MUCG46B

  1. Provenance of code.
  2. Purpose of code.
  3. Specification.
  4. Description of subroutine's operation.
  5. References.
  6. Parameter descriptions.
  7. Error indicators.
  8. Accuracy estimate.
  9. Any additional information.
  10. Example of code
  11. Auxiliary subroutines required.
  12. Keywords.
  13. Download source code.
  14. Links.

Provenance of Source Code

Hong-Seok Yang, H. K. D. H. Bhadeshia
Computational Metallurgy Laboratory,
Graduate Institute of Ferrous Technology,
Pohang University of Science and Technology,
Pohang, South Korea

Top| Next

Purpose

MUCG46 has been enhanced in order to calculate martensite-start temperature including the effect of austenite grain size, stress and plastic strain in the austenite prior to its transformation.

Top| Next| Prev

Specification

Language: FORTRAN
Product form: Source code

PROGRAM MAP_STEEL_MUCG46B

DOUBLE PRECISION BDIF(40), BSH(40), BTEM(40), C(8), DFPRO(40), DDFTO(40),
& DT4(40), DXQ(40)

DOUBLE PRECISION A, A1, A44, AEQ, AFE, AFE44, AFEQ, AJ, AJ1, BS, CONST,
& CORR,CTEMP, DA44, DAFE44, DF441, DIFFT, ETEQ, ETEQ2, F, F44, FPRO,
& FPROA, FSON, FTO, GMAX, H, H1, MS, R, S, S1, SHEART, SLOPE, STRAIN, T,
& T10, T20, TEQ, TEQ2, V14, W, W1, WS, WS1, X, XA, X1, X44, XBAR, XEQ,
& XEQ2, XMS, XTO, XTO400

DOUBLE PRECISION MAP_STEEL_AFEG, MAP_STEEL_CG, MAP_STEEL_DAFEG,
& MAP_STEEL_DCG, MAP_STEEL_ENERGY, MAP_STEEL_FTO1, MAP_STEEL_WSFUN

INTEGER I, IK, IMAX, J1, J2, J5, J6, J8, J98, J99

Top| Next| Prev

Description

The program is same with MUCG46 but calculates martensite-start temperature as a function of chemical composition, prior austenite grain size, external stress and pre-strain

Top| Next| Prev

References

  1. Hong-Seok Yang and H.K.D.H. Bhadeshia, Austenite Grain Size and the Martensite-Start Temperature, Scripta Materialia, Vol. 60, 2009, 493-495.
  2. H.K.D.H. Bhadeshia, Theory of Transformations in Steels, 2021, CRC Press, pp.1-557.
  3. H.K.D.H. Bhadeshia, A thermodynamic analysis of isothermal transformation diagrams, Metal Science, 16(1982) 159-165.
  4. H.K.D.H. Bhadeshia, Bainite in Steels , Institute of Materials, London (2015), 1-450.
  5. H.K.D.H. Bhadeshia, Rationalisation of Shear Transformations in Steels, Acta Metallurgica, Vol. 29, 1981, pp. 1117-1130.
  6. H.K.D.H. Bhadeshia, Critical Assessment of Growth of Plates in Steels, Materials Science and Technology, Vol. 1, 1985, pp. 497-504.
  7. Hong-Seok Yang and H.K.D.H. Bhadeshia, 2011

Top| Next| Prev

Parameters

Input parameters

C- real array of dimension 8
C holds the composition of each alloy:-
C(1) is the Carbon concentration (in weight percent).
C(2) is the Silicon concentration (in weight percent).
C(3) is the Manganese concentration (in weight percent).
C(4) is the Nickel concentration (in weight percent).
C(5) is the Molybdenum concentration (in weight percent).
C(6) is the Chromium concentration (in weight percent).
C(7) is the Vanadium concentration (in weight percent).
C(8) is the Ferrite concentration (not entered).

Gs- Prior austenite grain size in um
Stress- External stress in MPa
Prestr- Pre-Strain


Output parameters

None
All output is to STDOUT.

Top| Next| Prev

Error Indicators

Allowable composition ranges limited - see Description.

Top| Next| Prev

Accuracy

No information supplied.

Top| Next| Prev

Further Comments

None.

Top| Next| Prev

Example

1. Program text

       N/A - complete program

2. Program data

1
0.1  0  0.5  0  0  0  0  50  2  100  1

3. Program results


-------------------------------------------------------------------------------------------------------------------------------
    
 *******NUMBER      1  *********
    C=  1.0000   SI=  0.2000   MN=  0.5000   NI=  0.0000   MO=  0.0000   CR=  0.5000    V=  0.0000
    C=  0.0447   SI=  0.0038   MN=  0.0049   NI=  0.0000   MO=  0.0000   CR=  0.0052    V=  0.0000

    Austenite grain size (micrometres) =  1.
    Applied stress (MPa) = 1000.
    Plastic strain in austenite =0.00
    
  CARBON CONTENT=   0.04475  T10= -0.007293  T20= -0.311364   WGAMMA=  8345.
-------------------------------------------------------------------------------------------------------------------------------
   FPRO  FPROA   GMAX  CTEMP  X NUCLEUS  FSON      XEQ    XEQ50   FTO      XTO      X44   XTO400     SHEART     DIFFT   
  -1858. -2395. -2882.  200.  0.92D-07  -2864.   0.1959   0.1946  -973.  0.0794   0.0600  0.0620    0.27D+05   0.34D+24
  -1732. -2253. -2740.  220.  0.21D-06  -2716.   0.1903   0.1890  -873.  0.0766   0.0574  0.0602    0.64D+04   0.19D+22
  -1609. -2112. -2596.  240.  0.44D-06  -2569.   0.1846   0.1833  -773.  0.0730   0.0549  0.0572    0.18D+04   0.17D+20
  -1487. -1972. -2442.  260.  0.88D-06  -2421.   0.1789   0.1774  -673.  0.0691   0.0524  0.0538    0.63D+03   0.23D+18
  -1367. -1833. -2294.  280.  0.17D-05  -2274.   0.1729   0.1715  -574.  0.0664   0.0498  0.0506    0.25D+03   0.44D+16
  -1249. -1697. -2146.  300.  0.30D-05  -2127.   0.1669   0.1653  -476.  0.0630   0.0473  0.0473    0.11D+03   0.12D+15
  -1133. -1561. -1998.  320.  0.53D-05  -1979.   0.1607   0.1590  -378.  0.0589   0.0448  0.0443    0.56D+02   0.45D+13
  -1020. -1427. -1850.  340.  0.90D-05  -1832.   0.1543   0.1525  -280.  0.0551   0.0422  0.0413    0.32D+02   0.22D+12
   -909. -1295. -1702.  360.  0.15D-04  -1685.   0.1476   0.1458  -183.  0.0513   0.0396  0.0383    0.21D+02   0.14D+11
   -801. -1165. -1554.  380.  0.23D-04  -1538.   0.1408   0.1388   -86.  0.0479   0.0369  0.0352    0.15D+02   0.11D+10
   -699. -1042. -1413.  400.  0.35D-04  -1397.   0.1339   0.1319     5.  0.0449   0.0345  0.0322    0.11D+02   0.11D+09
   -616.  -941. -1295.  420.  0.53D-04  -1280.   0.1279   0.1257     0.  0.0425   0.0331  0.0300    0.90D+01   0.12D+08
   -536.  -841. -1177.  440.  0.77D-04  -1163.   0.1216   0.1193     0.  0.0400   0.0317  0.0278    0.79D+01   0.17D+07
   -458.  -743. -1059.  460.  0.11D-03  -1046.   0.1152   0.1128     0.  0.0374   0.0302  0.0255    0.76D+01   0.28D+06
   -384.  -647.  -942.  480.  0.15D-03   -929.   0.1085   0.1060     0.  0.0348   0.0289  0.0234    0.81D+01   0.55D+05
   -313.  -554.  -824.  500.  0.21D-03   -813.   0.1017   0.0990     0.  0.0323   0.0275  0.0212    0.97D+01   0.13D+05
   -246.  -463.  -714.  520.  0.29D-03   -696.   0.0946   0.0918     0.  0.0299   0.0262  0.0189    0.13D+02   0.37D+04
   -185.  -375.  -590.  540.  0.38D-03   -580.   0.0872   0.0842     0.  0.0276   0.0249  0.0164    0.22D+02   0.14D+04
   -128.  -290.  -474.  560.  0.50D-03   -463.   0.0795   0.0763     0.  0.0252   0.0236  0.0136    0.45D+02   0.62D+03
    -79.  -209.  -362.  580.  0.65D-03   -347.   0.0715   0.0681     0.  0.0227   0.0224  0.0108    0.12D+03   0.38D+03
    -39.  -133.  -241.  600.  0.84D-03   -231.   0.0630   0.0594     0.  0.0201   0.0214  0.0079    0.68D+03   0.44D+03
    -11.   -63.  -125.  620.  0.11D-02   -116.   0.0542   0.0503     0.  0.0174   0.0189  0.0048    0.13D+05   0.15D+04
     -8.   -53.  -110.  640.  0.13D-02    -99.   0.0527   0.0488     0.  0.0175   0.0245  0.0046    0.19D+05   0.66D+03
     -0.    -3.   -18.  660.  0.16D-02     -6.   0.0452   0.0411     0.  0.0153   0.0275  0.0019    0.12D+09   0.24D+06
    
    ***** FTO VERSUS TEMPERATURE ****   
    
 WIDMANSTATTEN START TEMPERATURES    300. &    300.
 NUCLEATION LIMITED BAINITE START TEMP=   300. C
 MARTENSITE START TEMPERATURE=   194. C
-------------------------------------------------------------------------------------------------------------------------------

Top| Next| Prev

Auxiliary Routines

Steel Subroutines

MAP_STEEL_AXTO
MAP_STEEL_GMAAX
MAP_STEEL_OMEGA
MAP_STEEL_TTTT
MAP_STEEL_WSTINE

Steel Functions

MAP_STEEL_AFEG
MAP_STEEL_CG
MAP_STEEL_DAFEG
MAP_STEEL_DCG
MAP_STEEL_ENERGY
MAP_STEEL_FTO1
MAP_STEEL_WSFUN

Utility Subroutines

MAP_UTIL_ANALY
MAP_UTIL_REED
MAP_UTIL_REEDI

Top| Next| Prev

Keywords

Widmanstatten, bainite, martensite, start temperature

Top| Next| Prev

Download

Download source code where MS is function of composition, stress, plastic strain, austenite grain shape.

Download source code where change in MS due to austenite grain shape affected by strain, neglected (S11=S33=1).

Top| Prev

MAP originated from a joint project of the National Physical Laboratory and the University of Cambridge.

Top | MAP Homepage