*HEADING
40 MM NOTCHED BEAM 3 POINT BENDING, C3D8R
*** FINE MESH
*** UNITS: NEWTON, METER
** THICKNESS IS CHANGED TO MAKE THE CRACKED ELEMENTS CUBIC
*NODE
1,
1081,.9,
1321,.98,
1441,1.,
9,,.1
1089,.9,.1
1329,.98,.1
1449,1.,.1
25,,.2
1105,.9,.2
1345,.98,.2
1465,1.,.2
*NGEN,NSET=LHS
1,9,1
9,25,1
*NGEN,NSET=L900
1081,1089,1
1089,1105,1
*NGEN,NSET=LNOTCH
1321,1329,1
1329,1345,1
*NGEN,NSET=CENT
1441,1449,1
1449,1465,1
*NFILL,NSET=BACK
LHS,L900,36,30
L900,LNOTCH,8,30
LNOTCH,CENT,4,30
*NCOPY,OLD=BACK,CHANGE=2000,SHIFT,NEW=FRONT
0.,0.,.005

*NSET,NSET=NOUT
1,1449,1465
*NSET,NSET=CENSYM,GEN
1449,1465,1
3449,3465,1
*ELEMENT,TYPE=C3D8R
1,1,31,32,2,2001,2031,2032,2002
1100,1329,1359,1360,1330,3329,3359,3360,3330
*ELGEN,ELSET=BEAM
1,24,1,1,44,30,24
1100,16,1,1,4,30,16
*ELSET,GEN,ELSET=EOUT
1100,1163,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
2001,2
**
*STEP,INC=100,NLGEOM
*STATIC,RIKS
0.01,1.0,,0.02,,1465,2,-1.5E-03
*BOUNDARY
1465,2,2,-1.5E-03
3465,2,2,-1.5E-03
*MONITOR,DOF=2,NODE=1465
*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