*HEADING
 esc8sxr4: ISOPARAMETRIC AND SKEW REBAR UNDER BENDING  [SC8R]
*PREPRINT,ECHO=YES,HISTORY=NO,MODEL=NO
*NODE,NSET=ALL
  1,  0.0, 0.0, -1.0
  2, 10.0, 0.0, -1.0
  3, 10.0,10.0, -1.0
  4,  0.0,10.0, -1.0
 91,  0.0, 0.0, +1.0
 92, 10.0, 0.0, +1.0
 93, 10.0,10.0, +1.0
 94,  0.0,10.0, +1.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
 2.,3
*REBAR LAYER
 REBAR1,1.0,2.5,0,RBMAT1,0,1
 REBAR2,1.0,2.5,0,RBMAT2,90,1
 REBAR3,1.0,3.5355,0,RBMAT3,45,1
 REBAR4,1.0,3.5355,0,RBMAT4,135,1
*MATERIAL,NAME=ELMAT
*ELASTIC
 1., 
*MATERIAL,NAME=RBMAT1
*ELASTIC
 30.E6, 
*MATERIAL,NAME=RBMAT2
*ELASTIC
 30.E6, 
*MATERIAL,NAME=RBMAT3
*ELASTIC
 30.E6, 
*MATERIAL,NAME=RBMAT4
*ELASTIC
 30.E6, 
*BOUNDARY
 ALL,3,6
 LEFT,1
 BOT,2
*STEP,PERTURBATION
*STATIC
*EL FILE,REBAR=REBAR1
 S,E,
*OUTPUT,FIELD
*ELEMENT OUTPUT
 S,E,
*ELEMENT OUTPUT,REBAR
 S,E
*EL FILE,REBAR=REBAR2
 S,E,
*EL FILE,REBAR=REBAR3
 S,E,
*EL FILE,REBAR=REBAR4
 S,E,
*NODE FILE
 U,RF,COORD
*OUTPUT,FIELD
*NODE OUTPUT
 U,RF,COORD
*EL PRINT,REBAR=REBAR1
 S,E,
*EL PRINT,REBAR=REBAR2
 S,E,
*EL PRINT,REBAR=REBAR3
 S,E,
*EL PRINT,REBAR=REBAR4
 S,E,
*NODE PRINT
 U,
 RF,
 COORD,
*CLOAD
 2,1,5000.
 3,1,5000.
 3,2,5000.
 4,2,5000.
 92,1,5000.
 93,1,5000.
 93,2,5000.
 94,2,5000.
*ENERGY PRINT
*END STEP
*STEP,NLGEOM
 LARGE DISPLACEMENT ANALYSIS
*STATIC
 1.,1.
*EL FILE,REBAR=REBAR1
 S,E,
*EL FILE,REBAR=REBAR2
 S,E,
*EL FILE,REBAR=REBAR3
 S,E,
*EL FILE,REBAR=REBAR4
 S,E,
*NODE FILE
 U,RF,COORD
*OUTPUT,FIELD
*ELEMENT OUTPUT,REBAR
 S
*NODE OUTPUT
 U,RF,COORD
*EL PRINT,REBAR=REBAR1
 S,E,
*EL PRINT,REBAR=REBAR2
 S,E,
*EL PRINT,REBAR=REBAR3
 S,E,
*EL PRINT,REBAR=REBAR4
 S,E,
*NODE PRINT
 U,
 RF,
 COORD,
*CLOAD ,OP=NEW
  2,1,5000.
  3,1,5000.
  3,2,5000.
  4,2,5000.
 92,1,5000.
 93,1,5000.
 93,2,5000.
 94,2,5000.
*END STEP