*heading
 CONCRETE DAMAGED PLASTICITY: UNIAXIAL CYCLIC TEST (C3D8R)
 with ORIENTATION and PERTURBATION step.
*preprint, model=yes
*part, name=part-1
*node, nset=nall
 1, 0.0, 0.0, 0.0
 2, 1.0, 0.0, 0.0
 3, 1.0, 1.0, 0.0
 4, 0.0, 1.0, 0.0
 5, 0.0, 0.0, 1.0
 6, 1.0, 0.0, 1.0
 7, 1.0, 1.0, 1.0
 8, 0.0, 1.0, 1.0
*nset, nset=nback
 1,2,3,4
*nset, nset=nside
 1,4,5,8
*nset, nset=nbase
 1,2,5,6
*nset, nset=ndisp
 5,6,7,8
*nset, nset=ndisp2
 2,3,6,7
*nset, nset=ntop
 3,4,7,8
*nset, nset=nout
 5,
*element, type=c3d8r, elset=eall
 1, 1, 2, 3, 4, 5, 6, 7, 8
*solid section, elset=eall, material=alle, orientation=ori45
*orientation, name=ori45
 1.0, 1.0, 0.0, -1.0, 1.0, 0.0
 3, 0.0
*end part
*assembly, name=assembly-1
*instance, name=part-1-1, part=part-1
*end instance
*end assembly
*material,name=alle
*elastic
 2.648e+10, 0.167
*concrete damaged plasticity
 15.0,
*concrete compression hardening
 24.019e+6, 0.0000
 29.208e+6, 0.0004
 31.709e+6, 0.0008
 32.358e+6, 0.0012
 31.768e+6, 0.0016
 30.379e+6, 0.0020
 28.507e+6, 0.0024
 21.907e+6, 0.0036
 14.897e+6, 0.0050
  2.953e+6, 0.0100
*concrete compression damage
 0.0000, 0.0000
 0.1299, 0.0004
 0.2429, 0.0008
 0.3412, 0.0012
 0.4267, 0.0016
 0.5012, 0.0020
 0.5660, 0.0024
 0.7140, 0.0036
 0.8243, 0.0050
 0.9691, 0.0100
*concrete tension stiffening
 1.780e+6, 0.0000
 1.457e+6, 0.0001
 1.113e+6, 0.0003
 0.960e+6, 0.0004
 0.800e+6, 0.0005
 0.536e+6, 0.0008
 0.359e+6, 0.0010
 0.161e+6, 0.0020
 0.073e+6, 0.0030
 0.040e+6, 0.0050
*concrete tension damage
 0.0000, 0.0000
 0.3000, 0.0001
 0.5500, 0.0003
 0.7000, 0.0004
 0.8000, 0.0005
 0.9000, 0.0008
 0.9300, 0.0010
 0.9500, 0.0020
 0.9700, 0.0030
 0.9900, 0.0050
*boundary
 assembly-1.part-1-1.nback,3,3
 assembly-1.part-1-1.nside,1,1
 assembly-1.part-1-1.nbase,2,2
*step, nlgeom, unsymm=yes
*static
 0.05,1.0,,0.05
*boundary
 assembly-1.part-1-1.ntop,2,2,0.001
*el print, elset=assembly-1.part-1-1.eall
 PEEQ, PEEQT
 DAMAGET, DAMAGEC, SDEG
*output, field
*element output, elset=assembly-1.part-1-1.eall
 PEEQ, PEEQT, DAMAGET, DAMAGEC
*node output, nset=assembly-1.part-1-1.nall
u
*output, history
*element output, elset=assembly-1.part-1-1.eall
 PE, LE, PEEQ, PEEQT, DAMAGET, DAMAGEC
*node output, nset=assembly-1.part-1-1.nall
u
*end step
*step, nlgeom, unsymm=yes
*static
 0.05,1.0,,0.05
*boundary, op=new
 assembly-1.part-1-1.nback,3,3
 assembly-1.part-1-1.nside,1,1
 assembly-1.part-1-1.nbase,2,2
*end step
*step, nlgeom, unsymm=yes
*static
 0.05,1.0,,0.05
*boundary
 assembly-1.part-1-1.ntop,2,2,-0.005
*end step
*step, nlgeom, unsymm=yes
*static
 0.05,1.0,,0.05
*boundary, op=new
 assembly-1.part-1-1.nback,3,3
 assembly-1.part-1-1.nside,1,1
 assembly-1.part-1-1.nbase,2,2
*end step
*step, perturbation
*static
*boundary, op=new
 assembly-1.part-1-1.nback,3,3
 assembly-1.part-1-1.nside,1,1
 assembly-1.part-1-1.nbase,2,2
 assembly-1.part-1-1.ntop,2,2,0.001
*end step