[MAP Logo]

Materials Algorithms Project
Program Library


Subroutine MAP_STEEL_ALUM

  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.

Top | Next

Purpose

MAP_STEEL_ALUM estimates the amount of aluminium that remains in solid solution in ferrite in a low-alloy steel arc weld deposit.

Top | Next | Prev

Specification

Language:FORTRAN
Product form:Source code

SUBROUTINE MAP_STEEL_ALUM(OXYWT, ALTOT, ALSOL, MAXI, IU, IFAIL, CONST)

DOUBLE PRECISION OXYWT, ALTOT, ALSOL, CONST

INTEGER MAXI, IU, IFAIL

Top | Next | Prev

Description

The atomic weight of aluminium is 26.98 and that of oxygen is 16. At the beginning of reaction we have w(Al) of aluminium and w(O) of oxygen. A quantity x of aluminium reacts with the oxygen. For the reaction:-

2Al + 3O = Al2O3

We therefore have [w(Al) - x] of aluminium left unreacted, and [w(O) - {{16*3}/{26.98*2}}x] of oxygen left unreacted. The equilbrium constant is therefore:-

K = [1/{(w(Al) - x)2(w(O) - 0.89x)3}]

In the subroutine, w(O) = OXYWT, w(Al)=ALTOT and K0.5 = CONST. The value of CONST is derived by matching with data from [2-4].

Top | Next | Prev

References

1. H.K.D.H. Bhadeshia, Unpublished research, (1990).

2. G. Thewlis, Joining and Materials, 2, (1989), 25-31.

3. G. Thewlis, Joining and Materials, 2, (1989), 125-129.

4. F.C. Liao and S. Liu, American Welding Journal, 71, (1992), 94s-104s.

Top | Next | Prev

Parameters

Input parameters

ALTOT - real
ALTOT is the total weight percent of aluminium in weld.

OXYWT - real
OXYWT is the total weight percent of oxygen in weld.

MAXI - integer
MAXI is the maximum number of iterations to be performed.

CONST - real
CONST = K0.5. References [2-4] yield a value of 0.616D-04 for CONST.

Output parameters

ALSOL - real
ALSOL is the weight percent of aluminium in solid solution.

IU - integer
IU is the number of iterations performed.

IFAIL - integer
IFAIL is set to a value of 1 when the input data are outside the valid range.

Top | Next | Prev

Error Indicators

IFAIL set to 1 indicates that the calculations are outside the validated range of experimental data.

Top | Next | Prev

Accuracy

Predictiblity of soluble aluminium is ±50 parts per million by weight.

Top | Next | Prev

Further Comments

None.

Top | Next | Prev

Example

To calculate the soluble aluminium concentration when the total aluminium concentration is 0.02 by percentage weight, and the total oxygen concentration is 0.03 by percentage weight.

1. Program text

       DOUBLE PRECISION ALTOT, ALSOL, OXYWT, CONST
       INTEGER MAXI, IU, IFAIL
       READ (5,*) ALTOT, OXYWT, MAXI, CONST
       IFAIL=0
       CALL MAP_STEEL_ALUM(OXYWT, ALTOT, ALSOL, MAXI, IU, IFAIL, CONST)
       WRITE(6,1) ALTOT, ALSOL, OXYWT, IU
1      FORMAT('Total aluminium =',F8.5,' by weight percent'
&      'Dissolved aluminium =',F8.5,' by weight percent'
&      'Total oxygen =',F8.5,' by weight percent'
&      I4,' iterations were performed')
       STOP
       END

2. Program data

0.02   0.03   20   0.616D-04

3. Program results

Total aluminium = 0.0200 by weight percent
Dissolved aluminium = 0.00495 by weight percent
Total oxygen = 0.03000 by weight percent
6 iterations were performed

Top | Next | Prev

Auxiliary Routines

None.

Top | Next | Prev

Keywords

aluminium, solid solution, ferrite

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!

>