*HEADING
 CORNERING ROLLING INERTIA. 
 S4R ELEMENTS; USES 2D DEPENDENCY FILE  pstsa2shhia.res
*SYMMETRIC MODEL GENERATION,REVOLVE,NODEOFF=1000,ELEMENT=1000
 0.0, 0.0, 0.0,    0.0, 0.0, 1.0
 0.0, 1.0, 0.0
 360.0, 24
*NSET,NSET=SECTR,GEN
 100, 399
*NSET,NSET=SECTS,GEN
 400, 699
*NSET,NSET=SECT
 SECTR, SECTS
*ELSET,ELSET=SECTR,GEN
 100, 299
*ELSET,ELSET=SECTS,GEN
 300, 499
*ELSET,ELSET=SECT
 SECTR, SECTS
*****************************
*STEP,INC=100,NLGEOM
 PLANE STRAIN BOUNDARY CONDITIONS
*STEADY STATE TRANSPORT,INERTIA=YES
*BOUNDARY
 NSIDE, 3
 NRIM,  1, 3
*DLOAD
 ESHELLR, CENT, 3.6,   0.0, 0.0, 20.5,  0.0, 1.0, 0.0
*MOTION,TYPE=VELOCITY,ROTATION
 NTIRES, 10.0,   0.0, 0.0, 22.0,  0.0, 1.0, 22.0
*NODE PRINT,FREQ=0
*EL PRINT,ELSET=SECTR,F=100
 S, 
 E,
*EL PRINT,ELSET=SECTS,F=100
 S, 
*EL FILE,ELSET=SECT,F=100
 S, E
*OUTPUT,FIELD,FREQ=100
*ELEMENT OUTPUT,ELSET=SECT
S,E
*END STEP
*****************************
*STEP,INC=100,NLGEOM
 FREE AXIAL DIRECTION
*STEADY STATE TRANSPORT,INERTIA=YES
 0.25, 1.0
*BOUNDARY,OP=NEW
 NRIM,  1, 3
*DLOAD
 ESHELLR, CENT, 3.6,   0.0, 0.0, 20.5,  0.0, 1.0, 0.0
*MOTION,TYPE=VELOCITY,ROTATION
 NTIRES, 10.0,   0.0, 0.0, 22.0,  0.0, 1.0, 22.0
*NODE PRINT,NSET=SECTR,F=100
 U, 
*NODE PRINT,NSET=SECTS,F=100
 U, 
*NODE FILE,NSET=SECT,F=100
 U, 
*OUTPUT,FIELD,FREQ=100
*NODE OUTPUT,NSET=SECT
U,
*END STEP