*HEADING
 Test of REBAR, MGAX2 in substructure   
 SUBSTRUCTURE GENERATION FILE psupreb13_gen1.inp SHOULD BE RUN FIRST.           
**
*preprint,model=yes
*NODE
 1, 2., 0.
 6, 2., 5.
*ELEMENT,TYPE=Z1,ELSET=SUPER, FILE=psupreb13_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,FREQ=1
*ELEMENT OUTPUT,REBAR=RB40_2
 S,E
*EL PRINT,REBAR=RB40_2,FREQ=1
 S,E,RBFOR,RBANG,RBROT
*EL FILE,REBAR=RB40_2,FREQ=1
 S,E,RBFOR,RBANG,RBROT
*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
*OUTPUT,FIELD,FREQ=1
*ELEMENT OUTPUT,REBAR
 S,E
*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
*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