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
2 3
4
5 
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
* 
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 *) Remark: alpha is inplane angle between material coordinate system and element coordinate system. 
1
2
3
7 
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

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 Remark: alpha is inplane angle between material coordinate system and element coordinate system. 

*MAT_037  2 3 
current hardening slope current yield stress 

*MAT_054  1
2
3
4
5
* 
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 Remark: alpha is inplane angle between material coordinate system and element coordinate system. 
1
2
3
4
5 
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
* 
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 *)
Remark: 

*MAT_058  1 2 3 6 7 8
10
18
19
20 
damage in longitudinal direction damage in transverse direction damage in shear cosine(alpha) (see remark *) sine(alpha) (see remark *) total failure flag (1=not failed, 0=failed) local strain adirection local strain bdirection local shear strain (abplane) effective strain local strain cdirection local transverse shear strain (bcplane) local transverse shear strain (caplane) damage variable for transverse shear behavior strainrate in adirection strainrate in bdirection strainrate in abplane eff. strainrate.
Remarks:
2. strainrates are smoothed values based on DT. 
1 2 3 4 5 6 7 8 9 10 11 12 13
14
15
16
17
18
19
23
24
25
27 29 
w1 – damage in adirection w2 – damage in bdirection w3 – damage in cdirection w4 – damage in abshear w5 – damage in bcshear w6 – damage in cashear r1 – threshold adir r2 – threshold bdir r3 – threshold cdir r4 – threshold abdir r5 – threshold bcdir r6 – threshold cadir q11 first vector to define material coordinate system (based on elemental system) q12 first vector to define material coordinate system (based on elemental system) q13 first vector to define material coordinate system (based on elemental system) q31 second vector to define material coordinate system (base on elemental system) q32 second vector to define material coordinate system (base on elemental system) q33 second vector to define material coordinate system (base on elemental system) d1tot – local strain adir d2tot – local strain bdir d3tot – local strain cdir d4tot – local shear strain (abplane) d5tot – local shear strain (bcplane) d6tot – local shear strain (caplane) efail – failure flag (1=not failed, 0=failed)) effstn – effective strain strainrate in aadirection strainrate in bbdirection strainrate in ccdirection strainrate in abplane strainrate in bcplane strainrate in caplane eff. strainrate 
*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)
Remark: 
eqp 1
2
4
15
16
17 18 
check (avg. of 8 failure flags) longitudinal tension failure flag (0=failed) transverse tension failure flag thruthick. shear/long.tension failure flag thruthick. tension (delam.) failure flag thruthick. shear/transv. tension failure flag longitudinal compression failure flag transv. compression failure flag thruthick. compression failure flag 
*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 
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 8 10 11 12 13 18 19 20 21 
density damage effective strain rate tensile strengths shear strengths bending strengths failure function value failure function normal term failure function normal bending failure function normal shear 

*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
2 
bulking pressure increment due to energy loss damage plastic strain increment effective strain rate 

*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 16 813 1419 2123 
relative volume material direction cosines viscous stress components total strain components maximum shear strain components for LCA.LT.0 

*MAT_133  1 2 3 4 5 7
* 
back stress component alpha_11 back stress component alpha_22 cosine(alpha) (see remark *) sine(alpha) (see remark *) back stress component alpha_12 strain rate
Remark: 

*MAT_136
*MAT_136_STD *MAT_136_2017 
1 2
3 
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 1 2 
damage separation in tangential direction separation in normal direction 

*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 2 3 6
7
10
* 
cosine(alpha) (see remark *) sine(alpha) (see remark *) equivalent stress effective total strain rate – as running average (VP=0) effective plastic strain rate (VP=1) scale factor for anisotropic plastic strain rate (VP=1) TsaiWu or TsaiHill failure.
Remark: 
1
2
5
6
7
8
9
10 
effective total strain rate – as running average (VP=0) effective plastic strain rate (VP=1) current yield stress q11 first vector to define material coordinate system (based on elemental system) q12 first vector to define material coordinate system (based on elemental system) q13 first vector to define material coordinate system (based on elemental system) q31 second vector to define material coordinate system (base on elemental system) q32 second vector to define material coordinate system (base on elemental system) q33 second vector to define material coordinate system (base on elemental system) TsaiWu or TsaiHill failure 
*MAT_158  1 2 3 4
5
6
10
* 
damage in longitudinal direction damage in transverse direction damage in shear threshold failure value r1 (see remark *) threshold failure value r2 (see remark *) cosine(alpha) (see remark **) sine(alpha) (see remark **) total failure threshold failure value rs (see remark *) local strain adirection local strain bdirection local shear strain viscoelastic overstress elemental xx viscoelastic overstress elemental yy viscoelastic overstress elemental zz viscoelastic overstress elemental xy viscoelastic overstress elemental yz viscoelastic overstress elemental zx
Remark:
Remark: 

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

*MAT_169  1 2 3 4 5 6 
damage edge damage normal stress shear stress element size plastic strain rate (starting with release R6) 

*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 2 3 4 5 6 16 17 18 28 
accumulated plastic multiplier eq. plastic strain tension/compr. eq. plastic strain shear eq. plastic strain biaxial damage volumetric plastic strain plastic strain rate tension/compr. plastic strain rate shear plastic strain rate biaxial instabililty criterion 
1 2 3 4 5 6 16 17 18 28 
accumulated plastic multiplier eq. plastic strain tension/compr. eq. plastic strain shear eq. plastic strain biaxial damage volumetric plastic strain plastic strain rate tension/compr. plastic strain rate shear plastic strain rate biaxial instabililty criterion 
*MAT_188  1 2 
effective creep strain effective thermal strain 
1 2 
effective creep strain effective thermal 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 7 8 9 10 11 12 14 15
16 
plastic strain rate plastic work F element size temperature plastic failure strain triaxiality F2 (for FAILOPT=1 option) number of cycles F2>= 1.0 (for FAILOPT=1 option) fraction plas. work to heat 
5 8 9 10 11 12 13 14 16 18 19 
plastic strain rate plastic failure strain triaxiality Lode parameter plastic work F element size temperature F2 (for FAILOPT=1 option) fraction plas. work to heat number of cycles F2>= 1.0 (for FAILOPT=1 option) 
*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 1 2 3 4 5 6 7 8
9 11
12 13
16 17 
damage failure flag reference strain rate thickness GIIC strength S (Mode II) GIC strength T (Mode I) mixed mode separation at plateau start plastic separation in t1 direction plastic separation in t2 direction plastic separation in peel direction separation ut1 for increment separation ut2 for increment separation un for increment state variable (0…1 elastic, 1…2 plastic, 2…3 softening/damage) curve shape factor FG1 curve shape factor FG2 

*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 
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 2 3 4 5 6 7 8 9 
damage plastic arc length effective strain rate temperature yield stress damaged yield stress triaxiality gammac gammaf 

*MAT_258  1 2 3 4 5 8 9 11 12 13 14 15 16 
plastic strain rate damage element size/thickness1 (le/te1) bending indicator Omega CL fractur parameter Wc plastic thickness strain increment in thickness strain hardening variable R_1 hardening variable R_2 hardening variable R_3 element size/thickness (le/te) yield stress yield stress + isotropic hardening (sigma0+R1+R2+R3) 

*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 
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 
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 