[MAP Logo]

Materials Algorithms Project
Program Library


Program MAP_STEEL_MUCG46_90 F90

  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

H.K.D.H. Bhadeshia,
Phase Transformations Group,
Department of Materials Science and Metallurgy,
University of Cambridge,
Cambridge, U.K.

Conversion to Fortran 90 by S. Cardie
Address as above.

Top | Next

Purpose

A powerful suite of software for modelling of the thermodynamics and kinetics of solid-state transformations in steels.

Top | Next | Prev

Specification

Complete program

Language: FORTRAN 90
Product form: Source code

Top | Next | Prev

Description

MAP_STEEL_MUCG46_90 is a very powerful program for the modelling of transformations in steels. Many people use it to calculate a TTT diagram, but there is in fact an enormous amount of phase diagram and kinetic and thermodynamic information generated in the file OUT.DAT which is created each time the program is run. This program is essential for any of the solid state phase transformations from austenite to martensite, bainite, WidmanstŠtten ferrite, allotriomorphic ferrite etc.

The TTT diagram consists of two C-curves for the initiation of transformation [1]. The higher temperature C-Curve is for reconstructive ("diffusional") reactions such as allotriomorphic ferrite and bainite. The lower C-curve is for displacive reactions such as WidmanstŠtten ferrite, bainite and acicular ferrite [2, pp. 175-185 & 124-134, Chapter 6].

The program asks for compositions (in wt.%) of seven constituents to be supplied. There are maximum and minimum limits imposed on each constituent as shown :-

Input
Order		 Element Max. (wt.%) Min. (wt.%)
1 Carbon 2.0 0.001
2 Silicon 2.5 0.0
3 Manganese 3.5 0.0
4 Nickel 3.5 0.0
5 Molybdenum 1.5 0.0
6 Chromium 3.5 0.0
7 Vanadium 1.5 0.0

The program will return an error message and ask for the data to be entered again if these limits are exceeded.

The maximum number of different alloys which can be analysed in a single run is set by the variable J1.

Top | Next | Prev

References

  1. H.K.D.H. Bhadeshia, A thermodynamic analysis of isothermal transformation diagrams, Metal Science, 16 (1982) 159-165.
  2. H.K.D.H. Bhadeshia, Bainite in Steels, Institute of Materials, London (1992) 1-450

Top | Next | Prev

Parameters

Note that the KIND type for all real variables/constants is determined by the parameter real_8_30, which is set within the module MAP_Util_Global_Data.

Input parameters

The used has the option to enter data either via the keyboard, or from an ASCII data file, called IN.DAT. In the latter case, the data should appear as a single column of values. See below for an example.

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 (calculated).

J6 - integer
J6 is the identifying number of each alloy to be examined.

Output parameters

All output is to two data files, OUT.DAT and PLOT.DAT

The results contained therein are:

OUT.DAT :

Name Description
FPRO Driving force for allotriomorphic ferrite formation by a paraequilibrium mechanism. All driving forces are in units of J/mol and all calculations are for the paraequilibrium state.
FPROA  
GMAX Driving force for nucleation.
CTEMP Temperature in ˇC
FSON  
XEQ Paraequilibrium carbon concentration of austenite in mole fraction.
XEQ2 As above, but allowing for 50 J/mol of strain energy in the ferrite in order to do calculations for Widmanstatten ferrite.
FTO Driving force for diffusionless transformation.
XTO T-zero carbon concentration in mole fraction.
X44  
XTO44 As XTO. but with 400 J/mol of stored energy in the ferrite, to allow calculations on bainite
SHEART Times in seconds of the lower C-curve of the TTT diagram.
DIFFT Times in seconds of the upper C-curve of the TTT diagram.

PLOT.DAT contains SHEART, DIFFT and CTEMP, as defined above.

Top | Next | Prev

Error Indicators

Allowable composition ranges limited - see Description.

Top | Next | Prev

Accuracy

No information supplied.

Top | Next | Prev

Further Comments

See also MAP_STEEL_MUCG46, a Fortran 77 version of this program.

Top | Next | Prev

Example

1. Program text

       Complete Program

2. Program data

Example of file IN.DAT

0.39
2.05
0.00
4.08
0.00
0.00
0.00

3. Program results

Contents of file OUT.DAT

-------------------------------------------------------------------------------------------------------------------------------

                *******NUMBER    1           *********
      C=  0.3900      SI=  2.0500      MN=  0.0000      NI=  4.0800      MO=  0.0000      CR=  0.0000      V=  0.0000
      C=  0.0176      SI=  0.0395      MN=  0.0000      NI=  0.0376      MO=  0.0000      CR=  0.0000      V=  0.0000

                 CARBON CONTENT=   0.01756     T10=  0.367363     T20= -0.379855     WGAMMA=  8732.
-------------------------------------------------------------------------------------------------------------------------------
   FPRO  FPROA   GMAX  CTEMP  X NUCLEUS  FSON      XEQ     XEQ50  FTO      XTO        X44     XTO400      SHEART     DIFFT
  -2303. -2527. -2788.  200.  0.25D-07  -2752.   0.1880   0.1867 -1903.  0.0669  4  0.0505  0.0542  6   0.31D+05   0.38D+24
  -2167. -2385. -2637.  220.  0.57D-07  -2610.   0.1821   0.1807 -1780.  0.0653  3  0.0478  0.0513  5   0.77D+04   0.22D+22
  -2032. -2244. -2494.  240.  0.12D-06  -2468.   0.1761   0.1746 -1672.  0.0621  6  0.0452  0.0484  5   0.22D+04   0.19D+20
  -1898. -2104. -2352.  260.  0.24D-06  -2327.   0.1699   0.1684 -1552.  0.0588  6  0.0425  0.0455  5   0.76D+03   0.26D+18
  -1764. -1964. -2210.  280.  0.46D-06  -2185.   0.1636   0.1620 -1433.  0.0557  5  0.0398  0.0426  5   0.30D+03   0.51D+16
  -1631. -1825. -2068.  300.  0.84D-06  -2044.   0.1572   0.1554 -1313.  0.0525  5  0.0370  0.0397  5   0.13D+03   0.14D+15
  -1500. -1687. -1926.  320.  0.15D-05  -1902.   0.1505   0.1487 -1193.  0.0492  5  0.0343  0.0367  5   0.68D+02   0.52D+13
  -1369. -1549. -1784.  340.  0.25D-05  -1761.   0.1436   0.1417 -1073.  0.0459  5  0.0315  0.0338  4   0.38D+02   0.25D+12
  -1240. -1413. -1642.  360.  0.41D-05  -1620.   0.1365   0.1344  -953.  0.0427  5  0.0288  0.0308  4   0.25D+02   0.16D+11
  -1112. -1277. -1500.  380.  0.64D-05  -1479.   0.1290   0.1269  -833.  0.0394  5  0.0259  0.0278  4   0.17D+02   0.13D+10
   -999. -1157. -1374.  400.  0.99D-05  -1353.   0.1220   0.1197  -729.  0.0365  5  0.0237  0.0251  4   0.13D+02   0.12D+09
   -899. -1051. -1264.  420.  0.15D-04  -1242.   0.1154   0.1130  -639.  0.0340  5  0.0221  0.0234  3   0.10D+02   0.14D+08
   -800.  -946. -1150.  440.  0.22D-04  -1131.   0.1086   0.1061  -548.  0.0315  5  0.0205  0.0213  4   0.88D+01   0.19D+07
   -703.  -842. -1038.  460.  0.31D-04  -1019.   0.1016   0.0989  -458.  0.0290  5  0.0189  0.0192  4   0.84D+01   0.30D+06
   -608.  -739.  -926.  480.  0.44D-04   -908.   0.0943   0.0914  -367.  0.0265  5  0.0173  0.0169  4   0.88D+01   0.59D+05
   -514.  -637.  -815.  500.  0.60D-04   -797.   0.0867   0.0836  -277.  0.0251  4  0.0156  0.0145  4   0.10D+02   0.14D+05
   -423.  -537.  -702.  520.  0.82D-04   -686.   0.0788   0.0755  -186.  0.0229  4  0.0140  0.0119  4   0.14D+02   0.39D+04
   -335.  -440.  -590.  540.  0.11D-03   -576.   0.0704   0.0669   -96.  0.0205  4  0.0123  0.0093  4   0.22D+02   0.14D+04
   -250.  -344.  -479.  560.  0.14D-03   -465.   0.0617   0.0579    -5.  0.0178  5  0.0107  0.0063  5   0.43D+02   0.60D+03
   -171.  -252.  -372.  580.  0.19D-03   -355.   0.0525   0.0484    85.  0.0152  5  0.0090  0.0032  5   0.11D+03   0.34D+03
    -99.  -165.  -255.  600.  0.24D-03   -245.   0.0426   0.0382     0.  0.0123  5  0.0073  0.0000  5   0.51D+03   0.35D+03
    -39.   -84.  -144.  620.  0.31D-03   -136.   0.0321   0.0272     0.  0.0093  5  0.0055  0.0000  1   0.67D+04   0.85D+03
    -37.   -80.  -144.  640.  0.38D-03   -131.   0.0314   0.0265     0.  0.0094  1  0.0096  0.0000  1   0.50D+04   0.22D+03
     -7.   -27.   -54.  660.  0.47D-03    -49.   0.0229   0.0177     0.  0.0070  5  0.0095  0.0000  1   0.52D+06   0.31D+04

    ***** FTO VERSUS TEMPERATURE ****   

 WIDMANSTATTEN START TEMPERATURES    460. &    460.
 NUCLEATION LIMITED BAINITE START TEMP=   460. C
     MARTENSITE START TEMPERATURE=   319. C
-------------------------------------------------------------------------------------------------------------------------------

Top | Next | Prev

Auxiliary Routines

Modules:
MAP_Steel_Acquire
MAP_Steel_Activity
MAP_Steel_Free_Energy
MAP_Util_Global_Data

Top | Next | Prev

Keywords

widmanstatten, bainite, martensite, start temperature

Top | Next | Prev

Download

Download source code

Download compiled executable for Windows

Top | Prev

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

Top | MAP Homepage