*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=displacement,mode mix ratio=energy, mixed mode behavior=tabular, softening=linear
 0.8, 0., 0.
0.75, 1., 0.
 0.8, 0., 1.
 0.7, 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