*HEADING
 Mohr-Coulomb, plane strain tension/compression, simple shear
 Material 1: hadening on MC and softening on TC
 Elements tested CPE4R
*NODE,NSET=ALLN
  1, 2., 0.
  2, 3., 0.
  3, 3., 1.
  4, 2., 1.
*node,nset=n_ss
 21, 0., 0.
 22, 1., 0.
 23, 1., 1.
 24, 0., 1.
*nset,nset=nx0
1,4
*nset,nset=ny0
1,2
*nset,nset=ny1
3,4
*nset,nset=ny0_ss
21,22
*nset,nset=ny1_ss
23,24
*ELEMENT,TYPE=Cpe4r,ELSET=e11
 1, 1,2,3,4
*ELEMENT,TYPE=Cpe4r,ELSET=e31
21, 21,22,23,24
*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
*ELASTIC
 300.E3,.3
*MOHR COULOMB
 30.,20.
*MOHR COULOMB HARDENING
<sigc>, 0.
<sigh>, 1.0
*TENSION CUTOFF
<sigt>, 0.
<sigs>, 1.0
**
*BOUNDARY
n_ss,2,2
ny0_ss,1,1
nx0,1,1
ny0,2,2
*STEP,nlgeom=yes,unsymm=yes
** Load in compression
** Elements should yield on MC
*STATIC
 1.,20.,,1.
*BOUNDARY
ny1,2,2,-0.1
ny1_ss,1,1,0.1
*print,solve=no
*output,field,op=add
*element output,elset=alle
peeqt
*output,history,op=add
*element output,elset=alle
peeq,peeqt,pener
*END STEP
*STEP
** Load in the opposite direction
*STATIC
 1.,40.,,1.
*BOUNDARY
ny1,2,2,0.1
ny1_ss,1,1,-0.1
*END STEP