*HEADING
40 MM NOTCHED BEAM 3 POINT BENDING, C3D8R
*** COARSE MESH
*** UNITS: NEWTON, METER
** THICKNESS IS CHANGED TO MAKE THE CRACKED ELEMENTS CUBIC
*NODE
1,
91,.9,
111,.98,
121,1.,
3,,.1
93,.9,.1
113,.98,.1
123,1.,.1
7,,.2
97,.9,.2
117,.98,.2
127,1.,.2
*NGEN,NSET=LHS
1,3,1
3,7,1
*NGEN,NSET=L900
91,93,1
93,97,1
*NGEN,NSET=LNOTCH
111,113,1
113,117,1
*NGEN,NSET=CENT
121,123,1
123,127,1
*NFILL,NSET=BACK
LHS,L900,9,10
L900,LNOTCH,2,10
LNOTCH,CENT,1,10
*NCOPY,OLD=BACK,CHANGE=1000,SHIFT,NEW=FRONT
0.,0.,.02

*NSET,NSET=NOUT
1,123,127
*NSET,NSET=CENSYM,GEN
123,127,1
1123,1127,1
*ELEMENT,TYPE=C3D8R
1,1,11,12,2,1001,1011,1012,1002
113,113,123,124,114,1113,1123,1124,1114
*ELGEN,ELSET=BEAM
1,6,1,1,11,10,10
113,4,1,1
*ELSET,GEN,ELSET=EOUT
113,116,1
*SOLID SECTION,ELSET=BEAM,MATE=CONC
*MATERIAL,NAME=CONC
*DENSITY
2400., 
*ELASTIC
30.E9,.2
*CONCRETE DAMAGED PLASTICITY
38.0
*CONCRETE COMPRESSION HARDENING
2.0E7,0.0
3.0E7,0.0015
*CONCRETE TENSION DAMAGE, TYPE=DISPLACEMENT
0.0, 0.0
0.9, 7.4474474E-05
*CONCRETE TENSION STIFFENING,TYPE=DISPLACEMENT
3.33E6,0.0
3.33E4,7.4474474E-5
*BOUNDARY
CENSYM,1
1,2,3
1001,2
**
*STEP,INC=100,NLGEOM
*STATIC,RIKS
0.001,1.0,,0.02,,127,2,-1.5E-03
*BOUNDARY
127,2,2,-1.5E-03
1127,2,2,-1.5E-03
*MONITOR,DOF=2,NODE=127
*RESTART,WRITE
*OUTPUT, FIELD, VARIABLE=PRESELECT
*ELEMENT OUTPUT
PEEQ, PEEQT, DAMAGET, DAMAGEC, SDEG
*OUTPUT,HISTORY,VARIABLE=PRESELECT
*ELEMENT OUTPUT,ELSET=EOUT
S,E
*NODE OUTPUT,NSET=NOUT
U,RF
*END STEP