*HEADING
 ESC8SXR3: Isoparametric rebar under bending  [SC8R]
*PREPRINT,ECHO=NO,HISTORY=NO,MODEL=YES
*NODE,NSET=ALL
  1, 0.0, 0.0,-5.0
  2,10.0, 0.0,-5.0
  3,10.0,10.0,-5.0
  4, 0.0,10.0,-5.0
 91, 0.0, 0.0,+5.0
 92,10.0, 0.0,+5.0
 93,10.0,10.0,+5.0
 94, 0.0,10.0,+5.0
*NSET,NSET=LEFT
 1,4
 91,94
*NSET,NSET=BOT
 1,2
 91,92
*ELEMENT,TYPE=SC8R,ELSET=SHELL
 1, 1,2,3,4, 91,92,93,94
*SHELL SECTION,MATERIAL=ELMAT, ELSET=SHELL, POISSON=ELASTIC
 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,3
*STEP,PERTURBATION
*STATIC
*EL FILE
 S,
 E,
*OUTPUT,FIELD
*ELEMENT OUTPUT
 S,
 E,
*ELEMENT OUTPUT,REBAR
 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,1,+10000.
  3,1,+10000.
 92,1,-10000.
 93,1,-10000.
*END STEP
*STEP,NLGEOM
 LARGE DISPLACEMENT ANALYSIS
*STATIC
 1.,1.
*EL FILE
 S,
 E,
*OUTPUT,FIELD
*ELEMENT OUTPUT
 S,
 E,
*ELEMENT OUTPUT,REBAR
 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,1,+10000.
  3,1,+10000.
 92,1,-10000.
 93,1,-10000.
*END STEP