*heading
 CONCRETE DAMAGED PLASTICITY: SHEAR CYCLIC TEST (C3D8R)
*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=nlside
 1,4,5,8
*nset, nset=nbase
 1,2,5,6
*nset, nset=nfront
 5,6,7,8
*nset, nset=nrside
 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
*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
*density
 2400.0
*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.nbase,1,2
 assembly-1.part-1-1.ntop,2,2
**
**
*step, nlgeom=no
*dynamic,explicit
 ,1.0
*boundary,type=velocity
 assembly-1.part-1-1.ntop,1,1,0.002
*output, field
*element output, elset=assembly-1.part-1-1.eall
 S, LE, PE, PEEQ, PEEQT, DAMAGET, DAMAGEC, SDEG
*node output, nset=assembly-1.part-1-1.nall
 U, V, RF
*output, history, variable=preselect
*element output, elset=assembly-1.part-1-1.eall
 S, PE, E, PEEQ, PEEQT, DAMAGET, DAMAGEC, SDEG,
 ENER, ELEN
*energy output
 ALLDMD
*end step
**
**
*step, nlgeom=no
*dynamic,explicit
 ,1.0
*boundary,type=velocity
 assembly-1.part-1-1.ntop,1,1,-0.004
*end step
**
**
*step, nlgeom=no
*dynamic,explicit
 ,1.0
*boundary,type=velocity
 assembly-1.part-1-1.ntop,1,1,0.004
*end step