*HEADING
: ISOPARAMETRIC REBAR UNDER BENDING  [S4]
*PREPRINT,ECHO=YES,HISTORY=NO,MODEL=NO
*NODE,NSET=ALL
1, 
2,10.
3,10.,10.
4,0.,10.
*NSET,NSET=LEFT
1,4
*NSET,NSET=BOT
1,2
*ELEMENT,TYPE=S4,ELSET=SHELL
1, 1,2,3,4
*SHELL SECTION,MATERIAL=ELMAT, ELSET=SHELL
10.,3
*REBAR LAYER
REBAR1,0.1,2.5,-2.5,RBMAT1,0,1
*MATERIAL,NAME=ELMAT
*ELASTIC
3.E6, 
*MATERIAL,NAME=RBMAT1
*ELASTIC
30.E6, 
*NSET,NSET=RIGHT
2,3
*BOUNDARY
LEFT,1,6
RIGHT,4
RIGHT,6
*STEP,PERTURBATION
*STATIC
*EL FILE
 S,
 E,
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
E,
*EL FILE,REBAR=REBAR1
 S,E
*NODE FILE
 U,RF,COORD
*OUTPUT,FIELD
*NODE OUTPUT
U,RF,COORD
*EL PRINT
 S,
 E,
*EL PRINT,REBAR=REBAR1
 S,E
*NODE PRINT
 U,
 RF,
 COORD,
*CLOAD
2,5,-100000.
3,5,-100000.
*END STEP
*STEP,NLGEOM
LARGE DISPLACEMENT ANALYSIS
*STATIC
1.,1.
*EL FILE
 S,
 E,
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
E,
*EL FILE,REBAR=REBAR1
 S,E
*NODE FILE
 U,RF,COORD
*OUTPUT,FIELD
*NODE OUTPUT
U,RF,COORD
*EL PRINT
 S,
 E,
*EL PRINT,REBAR=REBAR1
 S,E
*NODE PRINT
 U,
 RF,
 COORD,
*CLOAD ,OP=NEW
2,5,-100000.
3,5,-100000.
*END STEP