[MAP Logo]

Materials Algorithms Project
Program Library



Subroutine MAP_STEEL_TOPOLOGY

  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,
Materials Science and Metallurgy, University of Cambridge, U. K.
Graduate Institute of Ferrous Technology (GIFT), POSTECH, Korea

E-mail: hkdb@postech.ac.kr

The original separate programs for each deformation have been combined into one by:

JaeYong Chae
Computational Metallurgy Laboratory,
Graduate Institute of Ferrous Technology (GIFT), POSTECH, Korea

E-mail: exfeel@postech.ac.kr

Added to MAP: July, 2006.

Top | Next

Purpose

Calculation of the change in grain surface area per unit volume and
grain edge length per unit volume as a function of a variety of deformations.

Top | Next | Prev

Specification

Language: FORTRAN

Top | Next | Prev

Description

The effect of plastic deformation on the grain boundary surface area per unit volume and edge length per unit volume is examined using two methods. First, by applying homogeneous deformations to tetrakaidecahedra in a variety of orientations, and then by using the principles of stereology. This program contains nine subroutines. Each one of subroutines calculates the change in grain surface area per unit volume and grain edge length per unit volume using one of the nine deformation methods.
The nine deformation types that this program can support are as follows :

  1. Plane strain deformation of plates
  2. Plane strain deformation of oriented plates
  3. Axisymmetric compression of tetrakaidecahedra
  4. Plane strain deformation of tetrakaidecahedra
  5. Cross rolling of tetrakaidecahedra
  6. Shear of tetrakaidecahedra
  7. Axisymmetric tension of tetrakaidecahedra
  8. Plane strain deformation combined with shear of tetrakaidecahedra
  9. Plane strain deformation of needle

Users can select one of the nine deformation types and then the program runs a corresponding routine. Each routine initially defines vectors of a plate or tetrakaidecahedra. Result vectors are calculated by producting a deformation matrix(vectors) to that vectors. The deformation matrix is determined by each one of the deformation methods. After that, the routine calculates area and length from result vectors and calculates their ratio. Each routine is run in many times(number of loop is depended on cases).

The name of source code is as follows :

topology.for

This is a source code of this program. For the execution, it needs to be compiled.
compile example :
g77 topology.for -o name.out

Top | Next | Prev

References

Top | Next | Prev

Parameters

A# - real array of dimension 3, # represents numbers.
A is a vector to represent initial plate or tetrakaidecahedron. plate has three vectors(A1~A3) and tetrakaidecahedron has six vectors.

S# - deformation matrix (3 by 3)
S is a 3 by 3 matrix. It is used to product with plate or tetrakaidecahedron vectors.

B# - real array of dimension 3
B is a vector that records deforamtion result.

C#, D# - real array of dimension 3
C#, D# is temporary vectors needed in various calculating.

AREA# - reals, # : 1~15(tetra) or 1~3(plate)
AREA# used to record area due to squares, hexagons, or rectangles.

L# - reals, # : 1~6
L# is used to record length.

AREATOT, AREARATIO - real
total area, area ratio.

LENTOT, LRATIO - real
total length, length ratio.

R - real array of demension 9
Rotation matrix

SEL - integer
SEL is used to select a menu.

Top | Next | Prev

Error Indicators

None.

Top | Next | Prev

Accuracy

No information.

Top | Next | Prev

Further Comments

None.

Top | Next | Prev

Example

1. Program Compile

 g77 [sourcecode name] -o [executable file name]

2. Program menu

 Select deformation type:
    1 : Plane strain deformation of plates
    2 : Plane strain deformation of oriented plates
    3 : Axisymmetric compression of tetrakaidecahedra
    4 : Plane strain deformation of tetrakaidecahedra
    5 : Cross rolling of tetrakaidecahedra
    6 : Shear of tetrakaidecahedra
    7 : Axisymmetric tension of tetrakaidecahedra
    8 : Plane strain deformation combined with Shear
        of tetrakaidecahedra
    9 : Plane strain deformation of Needle
 0 : QUIT
 

3. Program data

 Auto generated.

4. Program results

(This is the case of menu number 7)
****************************************************
 Result for
  Axisymmetric tension                           
  of randomly oriented tetrakaidecahedra
 The first tetrakaidecahedron is oriented as in
 the following reference:
  S. B. Singh and H. K. D. H. Bhadeshia
  Materials Science and Technology 14 (1998) 832-834.
 ****************************************************
    L/L_0     A/A_0     S_11      S_22      S_33        Eq
      1.00      1.00      0.00      0.00      0.00      0.00
      1.23      1.17      0.69     -0.35     -0.35      0.69
      1.65      1.39      1.10     -0.55     -0.55      1.10
      2.11      1.59      1.39     -0.69     -0.69      1.39
      2.59      1.78      1.61     -0.80     -0.80      1.61
      3.08      1.94      1.79     -0.90     -0.90      1.79
      3.57      2.10      1.95     -0.97     -0.97      1.95
      4.07      2.24      2.08     -1.04     -1.04      2.08
      4.56      2.38      2.20     -1.10     -1.10      2.20
      5.06      2.51      2.30     -1.15     -1.15      2.30
      5.55      2.63      2.40     -1.20     -1.20      2.40
      6.05      2.75      2.48     -1.24     -1.24      2.48
      6.55      2.86      2.56     -1.28     -1.28      2.56
      7.05      2.96      2.64     -1.32     -1.32      2.64
      7.55      3.07      2.71     -1.35     -1.35      2.71
...

Top | Next | Prev

Auxiliary Routines

No auxiliary routines

Top | Next | Prev

Keywords

steel, topology, metallography, deformation

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.

Top | MAP Homepage