*HEADING
 ANNULAR RING UNDER APPLIED TWIST ON INNER 
 AND OUTER SURFACES SUBSTRUCTURE USAGE
 SUBSTRUCTURE GENERATION FILE torsholcyl_cgax8_substruct_gen1.inp SHOULD 
 BE RUN FIRST.
*NODE
1, 1.0, 0.0
21, 2.0, 0.0
1021, 2.0, 1.0
1001, 1.0, 1.0
*NGEN, NSET=INNER
1, 1001, 100
*NGEN, NSET=OUTER
21,1021, 100
*ELEMENT,TYPE=Z1,ELSET=SUPER, FILE=torsholcyl_cgax8_substruct_gen1
1,1,21,101,121,201,221,301,321,401,421,501,521,601,621,701,721,
801,821,901,921,1001,1021
*SUBSTRUCTURE PROPERTY,ELSET=SUPER
*RESTART, WRITE
** STEP I -- SIMPLE LINEAR PERTURBATION
**
*STEP, PERTURBATION
*STATIC
*BOUNDARY
INNER, 5,5, 0.01
OUTER, 5,5, -0.01
*OUTPUT,VAR=PRESELECT,FIELD
*END STEP
**
** STEP II -- SIMPLE LINEAR STEP
**
*STEP
*STATIC
0.5, 1.0
*BOUNDARY
INNER, 5,5, 0.01
OUTER, 5,5, -0.01
*NODE PRINT, NSET=OUTER,TOTALS=YES
RF, 
*NODE PRINT, NSET=INNER,TOTALS=YES
RF, 
*OUTPUT,VAR=PRESELECT,FIELD
*SUBSTRUCTURE PATH,ENTER ELEMENT=1
*EL PRINT, ELSET=EALL
S, MISES
E,
COORD,E13,S13
*EL PRINT, ELSET=EALL, POSITION=AVERAGED AT NODES
E13, S13
*NODE PRINT
U,
*EL FILE
S, 
E,
COORD,
*NODE FILE
U,
*OUTPUT,VAR=PRESELECT,FIELD
*SUBSTRUCTURE PATH,LEAVE
*END STEP