*HEADING
Biaxial tension with Mohr-Coulomb tension cutoff
*NODE, NSET=QA_TEST_ALLNODES
1, 0.0,0.,0.
2, 1.0,0.,0.
3, 1.0,1.,0.
4, 0.0,1.,0.
5, 0.0,0.,1.
6, 1.0,0.,1.
7, 1.0,1.,1.
8, 0.0,1.,1.
**
*nset,nset=nx0
1,4,5,8
*NSET,NSET=nx1
2,3,6,7
*nset,nset=ny0
1,2,5,6
*NSET,NSET=nz0
1,2,3,4
*NSET,NSET=nz1
5,6,7,8
*ELEMENT,TYPE=C3D8R,ELSET=C3D_1
1, 1,2,3,4,5,6,7,8
*SOLID SECTION,ELSET=C3D_1,MATERIAL=M1
*MATERIAL,NAME=M1
*DENSITY
1000. 
*ELASTIC
2.E7,0.3
*mohr coulomb
30.,20.
*mohr coulomb HARDENING
1.E3,0.0
1.E3,1.0
*tension cutoff
1.e3,0.
1.e3,1.
**
*boundary
nx0,1,1
ny0,2,2
nz0,3,3
*STEP,nlgeom=yes,inc=100
*static
0.0001,1
*print,solve=no
*BOUNDARY
nx1,1,1,0.1
nz1,3,3,0.1
**
*OUTPUT, FIELD,freq=100
*ELEMENT OUTPUT
S,PEEQ,PEEQT
*NODE OUTPUT
U,
*OUTPUT,HISTORY
*element output,elset=c3d_1
peeq,peeqt,pener
*END STEP
*STEP,nlgeom=yes,inc=100
*static
0.1,1
*BOUNDARY
nx1,1,1,-0.1
nz1,3,3,-0.1
*END STEP