*node
1, 0., 0., 0.
2, 1., 0., 0.
3, 1., 1., 0.
4, 0., 1., 0.
5, 0., 0., 1.
6, 1., 0., 1.
7, 1., 1., 1.
8, 0., 1., 1.
*nset, nset=upper
5,6,7,8
*nset, nset=lower
1,2,3,4
*element, type=coh3d8, elset=QA_TEST_EALL
1, 1,2,3,4,5,6,7,8
*boundary
lower, 1,6
*cohesive section, elset=QA_TEST_EALL, material=mat, response=traction separation, thickness=geometry
*material, name=mat
*density
1.
*elastic, type=traction
2.e6, 1.e6, 1.e6
*damage initiation, criterion=maxs
1.1e5, 1.e5, 1.e5
*damage evolution,type=energy,mode mix ratio=energy, mixed mode behavior=tabular, softening=linear
100734.4, 0., 0.
  30000., 1., 0.
100734.4, 0., 1.
  20000., 1., 1.
*amplitude, name=disp
0., 0., 1., 1.
*step, nlgeom=yes
*dynamic, explicit
,1.
*boundary, type=displacement, amplitude=disp
upper,1,3,1.
*OUTPUT, HISTORY, TIME INTERVAL=0.01
*ELEMENT OUTPUT, ELSET=QA_TEST_EALL
S, SDEG, LE, NE, DMICRT, STATUS
MMIXDMI,MMIXDME
*end step