*HEADING
 Test of REBAR, MGAX2 in substructure   
 SUBSTRUCTURE GENERATION FILE psupreb8_gen1.inp SHOULD BE RUN FIRST.           
**
*NODE
 1, 2., 0.
 6, 2., 5.
*ELEMENT,TYPE=Z1,ELSET=SUPER, FILE=psupreb8_gen1
1,1,6
*SUBSTRUCTURE PROPERTY,ELSET=SUPER
*BOUNDARY
 1, 1 
 1, 2
 1, 5 
*STEP
 STEP 1: Axial pulling
*STATIC
*SLOAD
SUPER,AXIAL,1.0
*OUTPUT,FIELD,VAR=PRESELECT
*SUBSTRUCTURE PATH,ENTER ELEMENT=1
*OUTPUT,FIELD,VAR=PRESELECT
*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
*SUBSTRUCTURE PATH,LEAVE
*MONITOR,NODE=6,DOF=5,FRE=1
**
*END STEP
**
**
*STEP
 STEP 2: Twisting
*STATIC
*SLOAD
SUPER,TWIST,1.0
*SUBSTRUCTURE PATH,ENTER ELEMENT=1
*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
*SUBSTRUCTURE PATH,LEAVE
*MONITOR,NODE=6,DOF=5,FRE=1
**
*END STEP