*heading
 CONCRETE DAMAGED PLASTICITY: BIAXIAL TEST (C3D8R)
*preprint, model=yes
*node, nset=nall
*nset, nset=nback
*nset, nset=nlside
*nset, nset=nbase
*nset, nset=nfront
*nset, nset=nrside
*nset, nset=ntop
*nset, nset=nout
*element, type=c3d8r, elset=eall
*solid section, elset=eall, material=alle
*material,name=alle
*density
 2400.0
*elastic
 2.648e+10, 0.167
*concrete damaged plasticity
*concrete compression hardening
*concrete compression damage
*concrete tension stiffening
*concrete tension damage
*boundary
 nback,3,3
 nlside,1,1
 nbase,2,2
*step, nlgeom, unsymm=yes
*static
 0.05,1.0,,0.05
*boundary
 ntop,2,2,0.001
 nrside,1,1,0.001
*el print, elset=eall
 PEEQ, PEEQT
 DAMAGET, DAMAGEC, SDEG
*output, field, variable=preselect
*element output, elset=eall
 PEEQ, PEEQT, DAMAGET, DAMAGEC, SDEG
*output, history, variable=preselect
*element output, elset=eall
 S, PE, LE, PEEQ, PEEQT, DAMAGET, DAMAGEC, SDEG
*restart, write
*end step
*step, nlgeom, unsymm=yes
*static
 0.05,1.0,,0.05
*boundary, op=new
 nback,3,3
 nlside,1,1
 nbase,2,2
*el print, elset=eall
 PEEQ, PEEQT
 DAMAGET, DAMAGEC, SDEG
*output, field, variable=preselect
*element output, elset=eall
 PEEQ, PEEQT, DAMAGET, DAMAGEC, SDEG
*output, history, variable=preselect
*element output, elset=eall
 S, PE, LE, PEEQ, PEEQT, DAMAGET, DAMAGEC, SDEG
*restart, write
*end step