*HEADING
 Mohr-Coulomb, Uniaxial compression followed by Uniaxial tension
 Material 1: hadening on MC and softening on TC
 Elements tested C3D8R and CAX4R
*NODE,NSET=ALLN
**
*nset,nset=nx0
*nset,nset=ny0
*nset,nset=nz0
*nset,nset=ny1
*ELEMENT,TYPE=C3D8r,ELSET=e11
*ELEMENT,TYPE=Cax4r,ELSET=e31
*elset,elset=e1
e11,e31
*else,elset=alle
e1
*SOLID SECTION,ELSET=e1,MATERIAL=m1
**
*parameter
 sigt = 1000
 sigs = 0.1*sigt
 sigc = sqrt(3)*sigt/2
 sigh = 2.0*sigc
**
** M1: hardening MC, softening TC
**
*MATERIAL,NAME=m1
*DENSITY
 1.0E-3
*ELASTIC
 300.E3,.3
*MOHR COULOMB
*MOHR COULOMB HARDENING
<sigc>, 0.
<sigh>, 1.0
*TENSION CUTOFF
<sigt>, 0.
<sigs>, 1.0
**
*BOUNDARY
nx0,1,1
ny0,2,2
nz0,3,3
*STEP,nlgeom=yes
** Load in compression
** Elements should yield on MC
*DYNAMIC,EXPLICIT
 ,1.
*BOUNDARY,TYPE=VELOCITY
ny1,2,2,-0.1
*output,field
*element output,elset=alle
peeqt
*node output
u
*output,history,variable=preselect
*element output,elset=alle
peeq,peeqt,pener
*END STEP
*STEP
** Load in tension
** Elements should yield on TC
*DYNAMIC,EXPLICIT
 ,2.
*BOUNDARY,TYPE=VELOCITY
ny1,2,2,0.1
*END STEP