*HEADING
 Test of REBAR in element MGAX1
**
**
*NODE
 1, 2., 0.
 6, 2., 5.
*NGEN,NSET=ALL
 1,6
**
*ELEMENT,TYPE=MGAX1
 1,1,2
*ELGEN,ELSET=MEMB
 1,1,,,5,1,1
*MEMBRANE SECTION,ELSET=MEMB,MATERIAL=MAT1
 0.05,
*REBAR LAYER
RB40,0.005,0.31416,0,MAT2,50,1  
*MATERIAL,NAME=MAT1
*ELASTIC
 1.E5,0.495
**
*MATERIAL,NAME=MAT2
*ELASTIC
 1.E8,0.495
**
*BOUNDARY
 1, 1
 1, 2
 1, 5
**
**
*STEP,NLGEOM,INC=100
 STEP 1: Axial pulling
*STATIC
*CLOAD
 6,2,200.
*EL PRINT,REBAR,FREQ=1
 S,E,RBFOR,RBANG,RBROT
*EL FILE,REBAR,FREQ=1
 S,E,RBFOR,RBANG,RBROT
*EL PRINT,ELSET=MEMB,FREQ=1
 S,E
 STH, 
*EL FILE,ELSET=MEMB,FREQ=1
 S,E
 STH, 
*NODE PRINT,FREQ=1
 U,CF,RF
*NODE FILE,FREQ=1
 U,CF,RF
*MONITOR,NODE=6,DOF=5,FRE=1
**
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=ALL
RF1,RF2
*OUTPUT,FIELD
*NODE OUTPUT
U,CF
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
*END STEP
**
**
*STEP,NLGEOM,INC=100
 STEP 2: Twisting
*STATIC
*CLOAD
 6,5,200.
*EL PRINT,REBAR,FREQ=1
 S,E,RBFOR,RBANG,RBROT
*EL PRINT,REBAR,POSITION=CENTROID,FREQ=1
 S,E,RBFOR,RBANG,RBROT
*EL FILE,REBAR,FREQ=1
 S,E,RBFOR,RBANG,RBROT
*EL FILE,REBAR,POSITION=CENTROID,FREQ=1
 S,E,RBFOR,RBANG,RBROT
*EL PRINT,ELSET=MEMB,FREQ=1
 S,E
 STH, 
*EL FILE,ELSET=MEMB,FREQ=1
 S,E
 STH, 
*NODE PRINT,FREQ=1
 U,CF,RF
*NODE FILE,FREQ=1
 U,CF,RF
*MONITOR,NODE=6,DOF=5,FRE=1
**
*OUTPUT,FIELD
*NODE OUTPUT
U,CF
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
*END STEP