*HEADING
 Mohr-Coulomb, plane strain tension/compression, simple shear
 Material 1: hadening on MC and softening on TC
 Elements tested CPE4R
*NODE,NSET=ALLN
*node,nset=n_ss
*nset,nset=nx0
*nset,nset=ny0
*nset,nset=ny1
*nset,nset=ny0_ss
*nset,nset=ny1_ss
*ELEMENT,TYPE=Cpe4r,ELSET=e11
*ELEMENT,TYPE=Cpe4r,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, DEPENDENCIES=1
*MOHR COULOMB HARDENING
<sigc>, 0.
<sigh>, 1.0
*TENSION CUTOFF
<sigt>, 0.
<sigs>, 1.0
**
*INITIAL CONDITIONS, TYPE=FIELD, VARIABLE=1
ALLN, 1.0
n_ss, 1.0
*BOUNDARY
n_ss,2,2
ny0_ss,1,1
nx0,1,1
ny0,2,2
*STEP,nlgeom=yes
** Load in compression
** Elements should yield on MC
*DYNAMIC,EXPLICIT
  ,1.
*BOUNDARY,TYPE=VELOCITY
ny1,2,2,-0.1
ny1_ss,1,1,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 the opposite direction
*DYNAMIC,EXPLICIT
 ,2.0
*BOUNDARY,TYPE=VELOCITY
ny1,2,2,0.1
ny1_ss,1,1,-0.1
*END STEP