History Variables for Certain Material Models
Below, a table of history variables of certain material models available in LSDYNA is given. One must set NEIPH and/or NEIPS in *DATABASE_EXTENT_BINARY in order to write these history variables, excepting "eqp".
There is no guarantee for completeness and correctness in connection with the release you use. For binding information please contact your local LSDYNA distributor.
In the following description, the notation 'eqp' corresponds to the variable 'effective plastic strain' (eps) in LSPrepost.
material 
no.  shells (and thick shells type=1, 2, 6) 
no. 
solids (and thick shells type=3, 5, 7) 

*MAT_001  eqp  volumetric strain  eqp  volumetric strain 
*MAT_002  eqp 1 2 3 4 5 6 7 8 9 10 11 12 
stiffness component C11 stiffness component C12 stiffness component C13 stiffness component C14 stiffness component C22 stiffness component C23 stiffness component C24 stiffness component C33 stiffness component C34 stiffness component C44 stiffness component C55 stiffness component C56 stiffness component C66 
1..9  deformation gradient 
*MAT_003  1 2 3 4 5 
back stress component xx back stress component yy back stress component xy back stress component yz back stress component zx 
1 2 3 4 5 
back stress component xx back stress component yy back stress component xy back stress component yz back stress component zx 
*MAT_004  1 3 
reference temperature current temperature 
1 3 
reference temperature current temperature 
*MAT_005  eps  plastic volumetric strain  
*MAT_007  1..9  deformation gradient  
*MAT_015  1 3 4 5 6 
failure value current pressure cutoff JC damage parameter temperature change JC failure strain 
5 6 
temperature change JC damage parameter 
*MAT_017  4 11 12 
fracture flag fracture stress element size 

*MAT_019 
1 
effective (plastic) strain rate failure ratio = hisvar#3/failure value current value of component upon which failure is based failure value of component (avail. for VP=0 only) failure flag (1=not failed; 0=failed) 
same as shells  
*MAT_021  15  temperature  10 11 12 13 14 15 20 
first component of first vector to define material coordinate system (based on global system) second component of first vector to define material coordinate system (based on global system) third component of first vector to define material coordinate system (based on global system) first component of second vector to define material coordinate system (based on global system) second component of second vector to define material coordinate system (based on global system) third component of second vector to define material coordinate system (based on global system) current temperature 
*MAT_022 (LAMSHT=1) 
eqp 1 2 3 4 * 
flag for longitudinal tensile failure mode flag for transverse tensile failure mode flag for transverse compressive failure mode cosine(alpha) (see remark *) sine(alpha) (see remark *)

1 
flag for longitudinal tensile failure mode flag for transverse tensile failure mode flag for transverse compressive failure mode flag for delamination failure mode 
*MAT_024  1 4 5 
effective strain rate (VP=0); effective plastic strain rate (VP=1) current hardening slope current yield stress Remark: 4 and 5 apply only to nonstrainrate option. 
1 4 5 
effective strain rate (VP=0); effective plastic strain rate (VP=1) current hardening slope current yield stress Remark: 4 and 5 apply only to nonstrainrate option. 
*MAT_025  1 2 3 4 5 6 7 
hardening variable K cap  J1axis interception, X(K) volumetric plastic strain first stress invariant second stress invariant response mode number number of iterations 

*MAT_026  eqp  relative volume  
*MAT_030  7 8 
martensite fraction xi_s loading function F 
18 19 
martensite fraction xi_s loading function F 
*MAT_034 (FORM =+14/14/24) 
eqp 1 2 3 4 5 68 911 28 29 
strain in adirection strain in bdirection strain in abdirection local stress fiber a local stress fiber b local shear stress internal strains some damping stresses effective strain (von Mises) max. effective strain (von Mises) 

*MAT_036  1 3 4 6 7 8 9 10 11 * 
reverse of old yield stress cosine(alpha) (see remark *) sine(alpha) (see remark *) effective strain rate current yield stress current hardening slope back stress component alpha_11 back stress component alpha_22 back stress component alpha_12


*MAT_037  2 3 
current hardening slope current yield stress 

*MAT_054 
1 2 3 4 5 6 8 9 10 11 12 16 * 
flag for longitudinal tensile failure mode flag for longitudinal compressive failure mode flag for transverse tensile failure mode flag for transverse compressive failure mode total failure damage parameter (SOFT) cosine(alpha) (see remark *) sine(alpha) (see remark *) local strain adirection local strain bdirection local shear strain (abplane) transverse shear damage

1 2 3 4 5 6 15 16 17 
flag for longitudinal tensile failure mode flag for longitudinal compressive failure mode flag for transverse tensile failure mode flag for transverse compressive failure mode total failure damage parameter (SOFT) local strain adirection local strain bdirection local shear strain (abplane) 
*MAT_055 
1
2 3 4 5 6 8 9 * 
flag for longitudinal tensile failure mode flag for longitudinal compressive failure mode flag for transverse tensile failure mode flag for transverse compressive failure mode total failure damage parameter (SOFT) cosine(alpha) (see remark *) sine(alpha) (see remark *)


*MAT_058 
1

damage in longitudinal direction
2. strainrates are smoothed values based on DT. 
1 
w1  damage in adirection 
*MAT_059  4 5 6 7 8 
cosine(alpha) (see remark) sine(alpha) (see remark) failure value (1=intact, 0=failed) damage (1=intact, 10^8=crashfront,1=failed) failure parameter f for ellipsoidal failure surface (applies only when SC<0)

eqp 
check (avg. of 8 failure flags) 
*MAT_063 (Model=0)  eqp 1 
volumetric strain ln(V) V1 

*MAT_063 (Model=1)  eqp 1 
plastic strain effective strain rate 

*MAT_065  6  temperature change  6  temperature change 
*MAT_072R3  1 2 3 4 5 6 
internal energy pressure from bulk viscosity volume from previous time step plastic volumetric strain slope of eta vs. lambda curve lambda 

*MAT_077 
eqp  strain energy density  
*MAT_081  1 4 5 
damage current hardening slope current yield stress 
1 4 5 
damage current hardening slope current yield stress 
*MAT_082  1 4 5 6 7 8 10 11 12 13 * 
damage 1 current hardening slope current yield stress damage 2 plastic strain in 1 plastic strain in 2 cosine(alpha) (see remark *) sine(alpha) (see remark *) Wilkins damage D (only _RCDC) Fracture fraction F (only _RCDC) Remark: 
1 2 3 6 7 8 9 10 23 24 
damage 1 damage 2 damage 3 current hardening slope current yield stress plastic strain in 1 plastic strain in 2 plastic strain in 3 Wilkins damage D (only _RCDC, 11/1/12) Fracture fraction F (only _RCDC, 11/1/12) 
*MAT_083  eqp 14 15 16 17 18 19 20 
strain energy density tension cutoff stress TC 1st principal strain rate 2nd principal strain rate 3rd principal strain rate 1st principal nominal strain 2nd principal nominal strain 3rd principal nominal strain 

*MAT_089 
eqp
1 4 5 7 19 21 
high tide value of total eff strain
effective strain rate filtered strain rate maximum strain rate failure strain total effective strain damage 
eqp
1 4 5 7 19 21 
high tide value of total eff strain
effective strain rate filtered strain rate maximum strain rate failure strain total effective strain damage 
*MAT_096  eqp 1 2 3 4 
normalized internal variable alpha internal variable alpha 1direction of crack normal 2direction of crack normal 3direction of crack normal 

*MAT_098  1 1 
max(0,ln(normalized eff strain rate)) (VP=0) failure flag (VP=1) 

*MAT_100  8 10 
damage effective strain rate for OPT=1,2,7,8,12,22 

*MAT_100_DA 
1 
density 

*MAT_101  5 7 9 
strain rate failure flag current unloading modulus 

*MAT_103  13 14 * 
cosine(alpha) (see remark *) sine(alpha) (see remark *) Remark: 

*MAT_104  3 4 
damage strain localization indicator 
3  damage 
*MAT_105  1  damage  1  damage 
*MAT_106  1 2 3 4 5 6 9 
temperature Young's modulus Poisson's ratio yield stress isotropic scaling factor kinematic scaling factor effective plastic strain rate 
1 2 3 4 5 6 9 
temperature Young's modulus Poisson's ratio yield stress isotropic scaling factor kinematic scaling factor effective plastic strain rate 
*MAT_110 
1 
bulking pressure increment due to energy loss 

*MAT_111  1  damage  
*MAT_120  1 4 5 6 7 11 13 14 
void volume fraction triaxiality effective strain rate growth of voids nucleation of voids material damage deviatoric microscopic pl. strain volumetric macroscopic pl. strain 
1 2 3 4 5 11 13 14 
void volume fraction triaxiality effective strain rate growth of voids nucleation of voids material damage deviatoric microscopic pl. strain volumetric macroscopic pl. strain 
*MAT_120_JC  1 4 5 6 7 8 9 10 11 13 14 
void volume fraction triaxiality effective strain rate growth of voids nucleation of voids JC failure strain JC damage domain variable material damage deviatoric microscopic pl. strain volumetric macroscopic pl. strain 
1 2 3 4 5 6 7 8 11 13 14 
void volume fraction triaxiality effective strain rate growth of voids nucleation of voids JC failure strain JC damage domain variable material damage deviatoric microscopic pl. strain volumetric macroscopic pl. strain 
*MAT_120_RC..  1 2 3 6 7 8 12 
void volume fraction damage critical damage fracture fraction trial stress weighting term asymmetric strain weighting term damage slope 
1 2 3 4 5 6 10 
void volume fraction damage critical damage fracture fraction trial stress weighting term asymmetric strain weighting term damage slope 
*MAT_123  1 6 7 
effective strain rate thinning strain major in plane strain 
1 7 
effective strain rate major principal strain 
*MAT_126 
eqp 
relative volume 

*MAT_133 
1 
back stress component alpha_11 

*MAT_136 *MAT_136_STD *MAT_136_2017 
1 2 3 4 5 6 7 15 
equivalent stress theta (angle between principal stress system and rolling direction) cosine(alpha) sine(alpha) 11component of backstress 22component of backstress 12component of backstress strain rate (von Mises equivalent strain) 

*MAT_138 
eqp 
damage 

*MAT_143  7 8 9 10 
damage, parallel tension damage, parallel compression damage, perpendicular tension damage, perpendicular compression 

*MAT_147  2 3 4 5 6 
current friction angle
effective strain damage damage criterion threshold current damage criterion 

*MAT_155  4 5 6 7 
tensile pressure cutoff (set to 0 if PCUTT or PCUTC are reached and SCALEP=0). pressure cutoff flag (initially 1 but set to 0 if PCUTT or PCUTC are reached). failure mode (0=no failure; 1=compressive failure; 2=tensile failure; 3=failure by plastic strain). current flow stress. 

*MAT_157 
1 
cosine(alpha) (see remark *) 
1 
effective total strain rate  as running average (VP=0) 
*MAT_158 
1 ** 
damage in longitudinal direction


*MAT_159  eqp 3 4 
misc. (see manual) ductile damage brittle damage 

*MAT_169 
1 
damage 

*MAT_181  13 eqp 
strain rate strain energy density 
13 eqp 
strain rate strain energy density 
*MAT_183  28  strain rate  13  strain rate 
*MAT_186 (TES=0/1) 
3 4 
separation in tangential direction separation in normal direction 

*MAT_186 (TES=2) 
1 3 4 5 
effective separation parameter stiffness in 1direction stiffness in 2direction stiffness in 3direction 

*MAT_187 
1 
accumulated plastic multiplier 
1 
accumulated plastic multiplier 
*MAT_188 
1 
effective creep strain 
1 
effective creep strain 
*MAT_221  eqp 1 2 3 4 5 6 7 8 
damage in traction along a damage in traction along b damage in traction along c damage in compression along a damage in compression along b damage in compression along c shear damage in abplane shear damage in bcplane shear damage in acplane 

*MAT_224 
1 16 
plastic strain rate 
5 
plastic strain rate 
*MAT_225  1 4 5 6 7 8 
plastic strain rate xx back stress yy back stress xy back stress yz back stress zx back stress 
1 4 5 6 7 8 
plastic strain rate xx back stress yy back stress xy back stress yz back stress zx back stress 
*MAT_234  1 2 3 18 19 * 
fill yarn angle (in degrees) warp yarn angle (in degrees) degradation coefficient cosine(alpha) (see remark *) sine(alpha) (see remark *) Remark: 

*MAT_235  1 2 3 4 9 10 13 15 16 
braid angle (fill) braid angle (warp) undulation angle (fill) undulation angle (warp) strain in fill direction strain in warp direction discount factor stress in fill direction stress in warp direction 

*MAT_240 
eqp 
damage 

*MAT_243  1 3 4 6 7 8 * 
reverse of old yield stress cosine(alpha) (see remark *) sine(alpha) (see remark *) effective strain rate current yield stress current hardening slope Remark: 

*MAT_244  1 2 3 4 5 6 7 8 9 10 11 12 13 14 17 19 20 25 
amount austenite amount ferrite amount pearlite amount bainite amount martensite Vickers hardness yield stress ASTM grain size number Youngs's modulus plastic strain austenite plastic strain ferrite plastic strain pearlite plastic strain bainite plastic strain martensite temperature rate current temperature reference temperature plastic strain rate 
1 2 3 4 5 6 7 8 9 10 11 12 13 14 17 19 20 25 
amount austenite amount ferrite amount pearlite amount bainite amount martensite Vickers hardness yield stress ASTM grain size number Youngs's modulus plastic strain austenite plastic strain ferrite plastic strain pearlite plastic strain bainite plastic strain martensite temperature rate current temperature reference temperature plastic strain rate 
*MAT_249 (version R9) (Starting with R10, variable POSTV affects history variable output) 
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 
ID of 1st fiber ID of 2nd fiber ID of 3rd fiber shear angle (rad) btw 1st and 2nd fibers shear angle (rad) btw 2nd and 3rd fibers current direction of 1st fiber in global coordinates current direction of 1st fiber in global coordinates current direction of 1st fiber in global coordinates current direction of 2nd fiber in global coordinates current direction of 2nd fiber in global coordinates current direction of 2nd fiber in global coordinates current direction of 3rd fiber in global coordinates current direction of 3rd fiber in global coordinates current direction of 3rd fiber in global coordinates number of fibers length change of 1st fibers length change of 2nd fibers length change of 3rd fibers fiber angle fiber angle 

*MAT_252 
1 
damage plastic arc length effective strain rate temperature yield stress damaged yield stress triaxiality gammac gammaf 

*MAT_255  1 2 5 
Youngs's modulus Poisson's ratio temperature 
1 
Youngs's modulus Poisson's ratio temperature 
*MAT_258  1 2 3 4 5 8 9 11 12 13 14 15 16 
plastic strain rate 

*MAT_272  2 3 4 
internal energy per unit volume alpha (porosity value) damage value 

*MAT_273 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27

kappa equivalent strain plastic strain xx plastic strain yy plastic strain zz plastic strain xy plastic strain zy plastic strain xz kappa tension kdt kappa tension kdt1 kappa tension kdt2 kappa compression kdc kappa compression kdc1 kappa compression kdc2 damage variable tension omegaT damage variable compression omegaC strain rate factor compression factor, alphac equivalent strain  tension equivalent strain  compression total xx strain total yy strain total zz strain total xy strain total yz strain total xz strain equivalent strain w/out rate factor 

*MAT_277 
1 2 3 4 5 
temperature degree of cure chemical strain thermal strain initial temperature 
1 2 3 4 5 
temperature degree of cure chemical strain thermal strain initial temperature 
.
EOS 
material  no.  

EOS_JWL  *MAT_008  epspl 1 2 3 4 5 
Lighting time Internal Energy (E * ro_ini / mass * vol_element) Bulk Viscosity Previous Volume Burn Fraction (=1 when fully burned) Applies to AFTERBURN option only; see Manual 
EOS_IGNITION_AND_GROWTH_OF_REACTION_IN_HE  *MAT_009  4 7 9 10 
temperature burn fraction 1 / Ve 1 / Vp 
EOS_IGNITION_AND_GROWTH_OF_REACTION_IN_HE  *MAT_010  5 8 10 11 
temperature burn fraction 1 / Ve 1 / Vp 