*HEADING
40 MM NOTCHED BEAM 3 POINT BENDING, C3D8R
*** MEDIUM MESH
*** UNITS: NEWTON, METER
** THICKNESS IS CHANGED TO MAKE THE CRACKED ELEMENTS CUBIC
*NODE
*NGEN,NSET=LHS
*NGEN,NSET=L900
*NGEN,NSET=LNOTCH
*NGEN,NSET=CENT
*NFILL,NSET=BACK
LHS,L900,18,20
L900,LNOTCH,4,20
LNOTCH,CENT,2,20
*NCOPY,OLD=BACK,CHANGE=1000,SHIFT,NEW=FRONT
*NSET,NSET=NOUT
*NSET,NSET=CENSYM,GEN
*ELEMENT,TYPE=C3D8R
*ELGEN,ELSET=BEAM
*ELSET,GEN,ELSET=EOUT
*SOLID SECTION,ELSET=BEAM,MATE=CONC
*MATERIAL,NAME=CONC
*DENSITY
2400.,
*ELASTIC
30.E9,.2
*CONCRETE DAMAGED PLASTICITY
*CONCRETE COMPRESSION HARDENING
*CONCRETE TENSION DAMAGE, TYPE=DISPLACEMENT
*CONCRETE TENSION STIFFENING,TYPE=DISPLACEMENT
*BOUNDARY
CENSYM,1
1,2,3
1001,2
**
*STEP,INC=100,NLGEOM
*STATIC,RIKS
0.001,1.0,,0.02,,493,2,-1.5E-03
*BOUNDARY
493,2,2,-1.5E-03
1493,2,2,-1.5E-03
*MONITOR,DOF=2,NODE=493
*RESTART,WRITE
*OUTPUT, FIELD, VARIABLE=PRESELECT
*ELEMENT OUTPUT
PEEQ, PEEQT, DAMAGET, DAMAGEC, SDEG
*OUTPUT,HISTORY,VARIABLE=PRESELECTED
*ELEMENT OUTPUT,ELSET=EOUT
S,E
*NODE OUTPUT,NSET=NOUT
U,RF
*END STEP