-----------------------------------------------------------------------------
Inputs parameter
-----------------------------------------------------------------------------
** MATERIALS
** 
*parameter
rho0 = 2.4e-09
G = 16400
A = 0.30
N = 0.79
B = 1.73
M = 0.61
C = 0.005
edot0 = 1.0
T = 2.4
sigIMax = 9.5
sigFMax = 1.1
HEL = 80
PHEL = 45
beta = 1.0
D1 = 0.04
D2 = 1.0
efMax = 1.0
efMin = 0.001
K1 = 116000
K2 = -243000
K3 = 506000
FS = 0.2
lDamage = 0 
-----------------------------------------------------------------------------
Defining materail parameter
-----------------------------------------------------------------------------
*Material, name=abq_jh2_concrete
*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