Ref: https://link.springer.com/article/10.1007%2Fs12541-017-0213-2
-----------------------------------------------------------------------------
Inputs parameter
-----------------------------------------------------------------------------
** MATERIALS
** 
*parameter
rho0 = 2.53e-09
G = 30400
A = 0.93
N = 0.77
B = 0.2
M = 1
C = 0.005
edot0 = 1.0
T = 2.4
sigIMax = 8
sigFMax = 1
HEL = 5950
PHEL = 2920
beta = 1.0
D1 = 0.043
D2 = 0.85
efMax = 1.0
efMin = 0.001
K1 = 45400
K2 = -138000
K3 = 290000
FS = 0.2
lDamage = 0 
-----------------------------------------------------------------------------
Defining materail parameter
-----------------------------------------------------------------------------
*Material, name=abq_jh2_glass
*density
 <rho0>
*user material, constants=32
 <rho0>, <G>, <A>, <N>, <B>, <M>, <C>, <edot0>
 <T>, <sigIMax>, <sigFMax>, <HEL>, <PHEL>, <beta>
 <D1>, <D2>, <efMax>, <efMin>, <FS>, <lDamage>
 <K1>, <K2>, <K3>
*DEPVAR, DELETE=8
8
1, PEEQ , "Equivalent plastic strain"
2, PEEQ_RATE, "Equivalent plastic strain rate"
3, DUCTCRT , "Ductile damage initiation criterion"
4, DAMAGE , "Damage variable"
5, DELTAP , "Pressure increment due to bulking"
6, YIELD , "Yield strength"
7, EVOL    , "Volumetric Strain, Mu"
8, MPSTATUS , "Material point status"
-----------------------------------------------------------------------------
Defining General Contact 
-----------------------------------------------------------------------------
**Interaction:Int-1
*Contact,op=NEW
*Contact Inclusions,ALL EXTERIOR
*contact controls assignment,nodal erosion=Yes