*HEADING
40 MM NOTCHED BEAM 3 POINT BENDING, CPS4R
*** MEDIUM MESH
*** UNITS: NEWTON, METER
*NODE
1,
361,.9,
441,.98,
481,1.,
5,,.1
365,.9,.1
445,.98,.1
485,1.,.1
13,,.2
373,.9,.2
453,.98,.2
493,1.,.2
*NGEN,NSET=LHS
1,5,1
5,13,1
*NGEN,NSET=L900
361,365,1
365,373,1
*NGEN,NSET=LNOTCH
441,445,1
445,453,1
*NGEN,NSET=CENT
481,485,1
485,493,1
*NFILL,NSET=BACK
LHS,L900,18,20
L900,LNOTCH,4,20
LNOTCH,CENT,2,20
*NSET,NSET=NOUT
1,485,493
*NSET,NSET=CENSYM,GEN
485,493,1
*ELEMENT,TYPE=CPS4R
1,1,21,22,2
300,445,465,466,446
*ELGEN,ELSET=BEAM
1,12,1,1,22,20,12
300,8,1,1,2,20,8
*ELSET,GEN,ELSET=EOUT
300,316,1
*SOLID SECTION,ELSET=BEAM,MATE=CONC
0.05,
*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.75, 7.4474474E-5
*CONCRETE TENSION STIFFENING,TYPE=DISPLACEMENT
3.33E6,0.0
3.33E4,7.4474474E-5
*BOUNDARY
CENSYM,1
1,2
*RESTART,WRITE
**
*STEP,INC=100,NLGEOM
*STATIC,RIKS
0.001,1.0,,0.02,,493,2,-1.5E-03
*BOUNDARY
493,2,2,-1.5E-03
*MONITOR,DOF=2,NODE=493
*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