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.
A collection of procedures allowing the calculation of free energy change accompanying the transformation from austenite to ferrite of the same chemical composition, and its derivative with respect to the carbon concentration at the T0 boundary. The result incorporates a compensation for Zener ordering effects.
Language: | FORTRAN 90 |
Product form: | Source code. |
MODULE MAP_Steel_Free_Energy PRIVATE PUBLIC :: MAP_STEEL_FTO1,MAP_STEEL_G91 CONTAINS FUNCTION MAP_STEEL_FTO1(F,H,S,T,W,W1,X) USE MAP_Util_Global_Data IMPLICIT NONE REAL(KIND=real_8_30) :: MAP_STEEL_FTO1 REAL(KIND=real_8_30),INTENT(IN) :: F,H,S,T,W,W1,X END FUNCTION MAP_STEEL_FTO1 FUNCTION MAP_STEEL_G91(F,H,S,T,W,W1,XTO) USE MAP_Util_Global_Data IMPLICIT NONE REAL(KIND=real_8_30),INTENT(IN) :: F,H,S,T,W,W1,XTO REAL(KIND=real_8_30) MAP_STEEL_G91 END FUNCTION MAP_STEEL_G91 SUBROUTINE ZEN_CALC(T,X,ZENER,ZEN1,ZEN2) USE MAP_Util_Global_Data IMPLICIT NONE REAL(KIND=real_8_30),INTENT(IN) :: T,X REAL(KIND=real_8_30),INTENT(OUT) :: ZENER,ZEN1,ZEN2 END SUBROUTINE ZEN_CALC END MODULE MAP_Steel_Free_Energy
MAP_Steel_Free_Energy contains two public and one private procedures.
The PUBLIC procedures are :
The PRIVATE procedure, ZEN_CALC, calculates a Zener ordering term, which describes the ordering of carbon atoms, that leads the body-centered cubic lattice of ferrite becoming body-centered tetragonal; the degree of ordering increases with increasing carbon concentration and descreasing transformation temperature [1]. ZEN_CALC is used by MAP_STEEL_FTO1 and MAP_STEEL_G91.
MAP_STEEL_FTO1 and MAP_STEEL_G91:
ZEN_CALC:
For ZEN_CALC:
None.
No information supplied.
See the corresponding Fortran 77 procedures, MAP_STEEL_FTO1 and MAP_STEEL_G91.
A value for F could be obtained using, for example, MAP function MAP_STEEL_ENERGY.
PROGRAM MAP_FTO_TEST USE MAP_Util_Global_Data USE MAP_Steel_Free_Energy IMPLICIT NONE ! ! Variable declarations ! REAL(KIND=real_8_30) :: H,S,X,T,W,W1,F,& FTO, G9 READ(5,*) F,H,S,T,W,W1,X FTO = MAP_STEEL_FTO1(F,H,S,T,W,W1,X) G9 = MAP_STEEL_G91(F,H,S,T,W,W1,X) WRITE(6,1001) FTO WRITE(6,1002) G9 STOP 1001 FORMAT('Free energy change is ',F12.3) 1002 FORMAT('Differential is ',F12.3) END PROGRAM MAP_FTO_TEST
-2.8557824601061275E+03 3.8575000000000000E+04 13.4799995422363281 4.7300000000000000E+02 8.7324743196967374E+03 4.8570000000000000E+04 1.7560130583041827E-02
Free energy change is -1902.860 Differential is 66529.702
carbon, t0 boundary, ferrite, austenite, free energy change
MAP originated from a joint project of the National Physical Laboratory and the University of Cambridge.
MAP Website administration / map@msm.cam.ac.uk