Materials Algorithms Project
Program Library
Subroutine MAP_STEEL_TEXTURE
 Provenance of code.
 Purpose of code.
 Specification.
 Description of subroutine's operation.
 References.
 Parameter descriptions.
 Error indicators.
 Accuracy estimate.
 Any additional information.
 Example of code
 Auxiliary subroutines required.
 Keywords.
 Download source code.
 Links.
Saurabh Kundu and H. K. D. H. Bhadeshia
Materials Science and Metallurgy, University of Cambridge, U. K.
Email: sk459@cam.ac.uk
Added to MAP: September, 2006.
Top  Next
Calculation of the crystallographic texture obtained when austenite transforms into martensite. The austenite itself can be textured.
Top  Next  Prev
Top  Next  Prev
This program has the inputs consisting of the habit planes, displacement direction and the crystallographic orientation (gammaJalpha).
 It first calculates the angle theta (between applied stress and habit plane normal and beta (between the maximum shear stress and the shear stress component resolved on the habit plane).
 Then it calculates the mechanical free energy available for each martensite variant.
 It also measures the angle between trace of habit plane on the plane of observations with the applied stress direction. It must be noted that by default the sample XY plane is taken to be the plane of observations.
 The plot is made with respect to the principle sample axis which are taken as (100),(010) and (001)
 This program plots the martensite pole figure and also the pole figure for the habit planes active.
 User can plot as many variants as desired.
Inputs when running the program
IF AUSTENITE PUT 0 IF FERRITE PUT 1:
1
WHAT IS THE STARTING TEXTURE OF AUSTENITE:
0=RANDOM
1=COPPER
2=BRASS
3=GOSS
4=CUBE
5={110}<111>
6=MIXTURE
3
HOW MANY VARIANTS WORKING?
8
FIRST ELEMENT OF THE STRESS TENSOR:
200.0
FOR (100) POLE FIG PUT 100 :
FOR (111) POLE FIG PUT 111 :
FOR (110) POLE FIG PUT 110 :
100
HOW MANY GRAINS?
500
Name
crystal_habit_poly.for
This is a source code of this program. For the execution, it needs to be compiled.
compile example :
g77 crystal_habit_poly.for o name.out
Top  Next  Prev
Top  Next  Prev
 HABIT#  double precision array of dimension 24X3, # represents numbers.
 Matrix that stores 24 habit planes.
 DIRECTION#  Double precision array of dimension 24X3
 Matrix that stores 24 displacement directions for 24 variants of martensite.
 GM_J_AL#  Double precision array of dimension 72X3.
 It contains 24 matrix representing orientation relationship between martensite and austenite, following Bowles and MacKenzie notation can be written as (gamma J alpha).
 AL_J_GM#  Double precision array of dimension 72X3.
 It contains 24 matrix representing orientation relationship between martensite and austenite, following Bowles and MacKenzie notation can be written as (alpha J gamma).
 AL#  Double precision array of dimension 30X3.
 It stores all the plane indices for plotting pole figure.
 DEF_MAT#  Double precision array of dimension 3X3.
 DEF_MAT# This matrix stores stress tensor of the applied stress in the crystal axis.
 STRESS_TENSOR  Double precision array of dimension 3X3.
 Matrix represents the stress tensor in trhe sample axis.
 TRANS_MATRIX  Double precision array of dimension 3X3.
 Rotational matrix represnting the orientation relationship between any austenite grain with the sample axis.
 X_ST, Y_ST  DOUBLE PRECISION
 X and Y corordinate of any point in the pole figure.

Top  Next  Prev
None.
Top  Next  Prev
No information.
Top  Next  Prev
None.
Top  Next  Prev
1. Program Compile
g77 [sourcecode name] o [executable file name]
2. Program menu
3. Program data
Auto generated.
4. Program results
Top  Next  Prev
No auxiliary routines
Top  Next  Prev
Transformation texture, martensite, stainless steels, crysttalography
Top  Next  Prev
Download source code
Top  Prev
MAP originated from a joint project of the National Physical Laboratory and the University of Cambridge.
Top  MAP Homepage