[MAP Logo]

Materials Algorithms Project
Program Library

Subroutine MAP_STEEL_RATE2

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

To calculate the parabolic rate constant.

Language:	FORTRAN
Product form:	Source code

SUBROUTINE MAP_STEEL_RATE2(WDIFF,XBAR,TC,XQ,TEMPC,ALP,AGSIZE)

DOUBLE PRECISION WDIFF,XBAR,TC(10),XQ(10),TEMPC(10),ALP(10),AGSIZE

None supplied.

Gilmour et al., Metall. Trans., 3, (1972), 3213.
H.K.D.H. Bhadeshia et al., Physical Metallurgy of Steel Welds, Conference proceedings, Denver, (1987).

WDIFF - real: WDIFF is the carbon-carbon interaction energy in austenite, not allowing for the presence of molybdenum, chromium, or vanadium. WDIFF may be calculated using the routine MAP_STEEL_OMEGA.
XBAR - real: XBAR is the mole fraction of carbon in the alloy.
TC - real array of dimension 10: TC contains the set of temperatures at which calculations are to be made (in centigrade).
XQ - real array of dimension 10: XQ contains the equilibrium mole fractions of carbon in austenite for the temperatures in TC.
AGSIZE - real: AGSIZE is the austenite grain size (2a) in microns.

ANS - real: ANS is the diffusion coefficient is the parent phase - and the integral divided by XEQ - XBAR, where XEQ is the equilibrium mole fraction of carbon in austenite (in cm²s^-1).
ERROR - real: ERROR is the error (if using NAG routine D01GAF)
SOFTT - real: SOFTT is the time for soft impingement (in seconds).
DIFF() - real: DIFF() is the diffusivity of carbon in austenite (in cm²s^-1).
XALPHA - real: XALPHA is the equilibrium mole fraction of carbon in ferrite.