History Variables for Certain Material Models

Below, a table of history variables of certain material models available in LS-DYNA 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 LS-DYNA distributor.

In the following description, the notation 'eqp' corresponds to the variable 'effective plastic strain' (eps) in LS-Prepost.

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 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 ) Remark: alpha is in-plane 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 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 non-strain-rate 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 non-strain-rate option.
*MAT_025			1 2 3 4 5 6 7	hardening variable K cap - J1-axis 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 6-8 9-11 28 29	strain in a-direction strain in b-direction strain in ab-direction 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 in-plane 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 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 a-direction local strain b-direction local shear strain (ab-plane) transverse shear damage Remark: alpha is in-plane angle between material coordinate system and element coordinate system.	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 a-direction local strain b-direction local shear strain (ab-plane)
*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 ) Remark: alpha is in-plane angle between material coordinate system and element coordinate system.
*MAT_058	1 2 3 6 7 8 10 11 12 15 16 17 18 19 20 21 22 23 *	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 a-direction local strain b-direction local shear strain (ab-plane) effective strain local strain c-direction local transverse shear strain (bc-plane) local transverse shear strain (ca-plane) damage variable for transverse shear behavior strain-rate in a-direction strain-rate in b-direction strain-rate in ab-plane eff. strain-rate Remarks: 1. alpha is in-plane angle between material coordinate system and element coordinate system. 2. strain-rates 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 20 21 22 23 24 25 27 29 30 31 32 33 34 35	w1 - damage in a-direction w2 - damage in b-direction w3 - damage in c-direction w4 - damage in ab-shear w5 - damage in bc-shear w6 - damage in ca-shear r1 - threshold a-dir r2 - threshold b-dir r3 - threshold c-dir r4 - threshold ab-dir r5 - threshold bc-dir r6 - threshold ca-dir 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 a-dir d2tot - local strain b-dir d3tot - local strain c-dir d4tot - local shear strain (ab-plane) d5tot - local shear strain (bc-plane) d6tot - local shear strain (ca-plane) efail - failure flag (1=not failed, 0=failed)) effstn - effective strain strain-rate in aa-direction strain-rate in bb-direction strain-rate in cc-direction strain-rate in ab-plane strain-rate in bc-plane strain-rate in ca-plane eff. strain-rate
*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: alpha is in-plane angle between material coordinate system and element coordinate system.	eqp 1 2 3 4 15 16 17 18	check (avg. of 8 failure flags) longitudinal tension failure flag (0=failed) transverse tension failure flag thru-thick. shear/long.tension failure flag thru-thick. tension (delam.) failure flag thru-thick. shear/transv. tension failure flag longitudinal compression failure flag transv. compression failure flag thru-thick. compression failure flag
*MAT_063 (Model=0)			eqp 1	volumetric strain ln(V) V-1
*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: alpha is in-plane angle between material coordinate system and element coordinate system.	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 cut-off 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 1-direction of crack normal 2-direction of crack normal 3-direction 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: alpha is in-plane angle between material coordinate system and element coordinate system.
*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 3 4	bulking pressure increment due to energy loss damage plastic strain increment effective strain rate
*MAT_111			1 2 3 4 5 10 18 19 22 23 24 25	damage material density measure (1=fully dense) max. volumetric strain modified volumetric strain (mu bar) volumetric strain increment equivalent stress volumetric strain (mu) normalized pressure volumetric plastic strain volumetric elastic strain plastic strain to fracture strain rate
*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 1-6 8-13 14-19 21-23	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: alpha is in-plane angle between material coordinate system and element coordinate system.
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) 11-component of backstress 22-component of backstress 12-component 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 8 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) Tsai-Wu or Tsai-Hill failure Remark: alpha is in-plane angle between material coordinate system and element coordinate system.	1 2 3 4 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) Tsai-Wu or Tsai-Hill failure
*MAT_158	1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 * **	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 a-direction local strain b-direction 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: When a threshold variable reaches the value of 1.0, the corresponding failure criterion is reached. Note: The threshold failure value rs is only active if FS=-1.0 (faceted failure surface) is used. Remark: alpha is in-plane angle between material coordinate system and element coordinate system.
*MAT_159			eqp 3 4	misc. (see manual) ductile damage brittle damage
*MAT_163			eqp 1 9	logarithmic relative volume ln(V) negative volumetric strain V-1 volumetric strain rate
*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 1-direction stiffness in 2-direction stiffness in 3-direction
*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 ab-plane shear damage in bc-plane shear damage in ac-plane
*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: alpha is in-plane angle between material coordinate system and element coordinate system.
*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 10 11 12 13 14 15 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: alpha is in-plane angle between material coordinate system and element coordinate system.
*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 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_255	1 2 5	Youngs's modulus Poisson's ratio temperature	1 2 5	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 damage element size/thickness-1 (le/te-1) 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_271			1	relative density
*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