*heading
 CONCRETE DAMAGED PLASTICITY (C3D8R)
 Import: Standard to Standard 
** First Standard analysis
*preprint, model=yes
*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
*material,name=alle
*density
 2400.0
*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
 nback,3,3
 nside,1,1
 nbase,2,2
*restart,write
*step, nlgeom, unsymm=yes
*static
 0.05,1.0,,0.05
*boundary
 ntop,2,2,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
*end step
*step, nlgeom, unsymm=yes
*static
 0.05,1.0,,0.05
*boundary
 ntop,2,2,-0.005
*end step