*heading
 CONCRETE DAMAGED PLASTICITY: CYCLIC TEST (T3D2)
*preprint, model=yes
*part, name=part-1
*node, nset=nall
 1, 0.0, 0.0, 0.0
 2, 1.0, 0.0, 0.0
*element, type=t3d2, elset=eall
 1, 1, 2
*solid section, elset=eall, material=alle
 1.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
*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.1,pinned
 assembly-1.part-1-1.2,2,3
*amplitude,name=v1
 0.0,1.0,1.0,1.0
*amplitude,name=v2
 0.0,-1.0,5.0,-1.0
*amplitude,name=v3
 0.0,1.0,4.0,1.0
**
**
*step
*dynamic,explicit
 ,1.0
*boundary,amplitude=v1,type=velocity
 assembly-1.part-1-1.2,1,1,0.001
*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,LE,PEEQ,PEEQT, DAMAGET, DAMAGEC, SDEG,
 ENER, ELEN
*energy output
 ALLDMD
*end step
**
**
*step
*dynamic,explicit
 ,5.0
*boundary,amplitude=v2,type=velocity
 assembly-1.part-1-1.2,1,1,0.001
*end step
**
**
*step
*dynamic,explicit
 ,4.0
*boundary,amplitude=v3,type=velocity
 assembly-1.part-1-1.2,1,1,0.001
*end step