Home

Run "cubic_to_mono.m" if want twin and habit elements. That is, if want:

nI, aI, shearI, nII, aII, shearII,...
bplusI, mplusI, bminusI, mminusI, shapestrainI,...
bplusII, mplusII, bminusII, mminusII, shapestrainII, Ui, Uj, lambdaI, lambdaII, Q

for a given set of lattice parameters and for a given correspondence variant (CV) pair (e.g., for a twin composed of CVs 1 and 2). This code reports the elements in terms of the B2 cubic phase. If want things reported in B19' phase, go into code and uncomment lines 41, 47 and comment lines 42, 48. Note - this will report twin elements in the B19' phase but the habit elements won't work (haven't written "habit_elements_B19p.m").

Run "BhatConversions.m" if want the B19' monoclinic orientations of all possible austenite-compatible martensite twins (i.e., twins that can form compatible habit planes). Will rotate orientations according to an input austenite grain orientation. Will also report transformation strain in y-direction and work based on a uniaxial load in y-direction. Go into code and update variable "Stress" if know actual stress state.

Run "BhatConversionsTwinRelOnly.m" if want the B19' monoclinic orientations of all possible twins (i.e., not just the ones that can form compatible habit planes). This is for, e.g., self-accommodated martensite. Will also report the twinning strain in the y-direction.

Note: "BhatConversions.m" and "BhatConversionsTwinRelOnly.m" uses a special code called "CV2Hexrd.m" which converts the lattices outputted by the Crystallographic Theory of Martensite (CTM) calculations to appropriate B19' orientations in terms of the real B9' lattice vector system (a < b < c, beta ~= pi/2). This is done assuming the coordinate system that CHESS and HEXRD uses (y is loading direction, beam propagates in -z direction).