[MAP Logo]

Materials Algorithms Project
Program Library


Function MAP_STEEL_YB

  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

G.I. Rees,
Phase Transformations Group,
Department of Materials Science and Metallurgy,
University of Cambridge,
Cambridge, U.K.

Top | Next

Purpose

Called by the function MAP_STEEL_WKINETIC in the calculation of the volume fraction of Widmanstätten ferrite or bainite formed, when the nucleation rate and growth rate are known.

Top | Next | Prev

Specification

Language:FORTRAN
Product form:Source code

DOUBLE PRECISION FUNCTION MAP_STEEL_YB(S, T, TC, X)

DOUBLE PRECISION S, T, TC, X, MAP_STEEL_YB

Top | Next | Prev

Description

MAP_STEEL_YB is used to calculate an element in the numerical integration for the case t > tc, where:

Equation

where Vmax is the maximum extent to which transformation to bainite or Widmanstätten ferrite is possible, G is the growth rate, rho is the plate tip radius, I is the nucleation rate, and Sv is the surface area of austenite grain boundary per unit volume.

Equation

Top | Next | Prev

References

  1. M. Umemoto, A. Hiramatsu, A. Moriya, T. Watanabe, S. Nanba, N. Nakajima, G. Anan, and Y. Higo, I.S.I.J. Int., 32, (1992), 306-315.
  2. G.I. Rees, H.K.D.H. Bhadeshia, and T. Maurickx, Progress Report for SOLLAC, August 1995.

Top | Next | Prev

Parameters

Input parameters

S - real
S = 2.sqrt 2I(G*rho)0.5Sv (see above).

T - real
T is the time of reaction (in seconds).

TC - real
TC is the time to grow across a grain (in seconds).

X - real
X is a measure of the position in the integration, e.g. n/N, where N is the total number of steps, and n the the number of the given step.

Output parameters

MAP_STEEL_YB- real
MAP_STEEL_YB = {1 - exp {-S ({2/3}(TC(1 - X))1.5} + (TC - X*TC)1.5(T - TC)) }}

Top | Next | Prev

Error Indicators

None.

Top | Next | Prev

Accuracy

If the exponent in the expression for Yb is less than -173 then Yb is taken to be 1.

Top | Next | Prev

Further Comments

None.

Top | Next | Prev

Example

1. Program text

       DOUBLE PRECISION S, T, TC, X, YBX, MAP_STEEL_YB
       READ (5,*) S, T, TC, X
       YBX = MAP_STEEL_YB(S, T, TC, X)
       WRITE (6,10) YBX
10     FORMAT ('YB = ', D13.5)
       STOP
       END

2. Program data

0.45754D-03   0.27125D+02   0.95721D+01   0.39216D-01

3. Program results

YB = 0.32328D-01

Top | Next | Prev

Auxiliary Routines

None.

Top | Next | Prev

Keywords

trapezium rule

Top | Next | Prev

Download

Download source code

Top | Prev

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

MAP Website administration / map@msm.cam.ac.uk

Top | Index | MAP Homepage Valid HTML 3.2!