*HEADING
 ANNULAR RING UNDER APPLIED TWIST ON INNER AND OUTER SURFACES
 ELEMENT : CGAX6M
*RESTART, WRITE
*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
*NFILL, NSET=ALL
INNER, OUTER, 20, 1
*NSET,NSET=BOTTOM,GENERATE
1, 21
*ELEMENT, TYPE=CGAX6M, ELSET=EALL
1,1,3,203,2,103,102
2,1,203,201,102,202,101
3, 3, 5, 203, 4, 104, 103
4, 5, 205, 203, 105, 204, 104
*ELGEN, ELSET=EALL
1, 5, 4, 4, 2, 202, 20
2, 5, 4, 4, 2, 202, 20
3, 5, 4, 4, 2, 198, 20
4, 5, 4, 4, 2, 198, 20
*ELCOPY,SHIFT NODES=400,OLD SET=EALL,ELEMENT SHIFT=40,NEW SET=E2ALL
*ELSET,ELSET=ROW1,GENERATE
1,20,1
*ELCOPY,SHIFT NODES=800,OLD SET=ROW1,ELEMENT SHIFT=80,NEW SET=ROW5
*ELSET,ELSET=EALL
E2ALL, ROW5
*SOLID SECTION, ELSET=EALL, MATERIAL=MALL
*MATERIAL, NAME=MALL
*ELASTIC
1.E6, 0.3
*BOUNDARY
BOTTOM,1,2
**
** STEP I -- SIMPLE LINEAR PERTURBATION
**
*STEP, PERTURBATION
*STATIC
*BOUNDARY
INNER, 5,5, 0.01
OUTER, 5,5, -0.01
*END STEP
**
** STEP II -- SIMPLE LINEAR STEP
**
*STEP
*STATIC
0.5, 1.0
*BOUNDARY
INNER, 5,5, 0.01
OUTER, 5,5, -0.01
*EL PRINT, ELSET=EALL
S, MISES
E,
COORD,
*EL PRINT, ELSET=EALL, POSITION=AVERAGED AT NODES
E13, S13
*NODE PRINT
U,
*NODE PRINT, NSET=OUTER,TOTALS=YES
RF, 
*NODE PRINT, NSET=INNER,TOTALS=YES
RF, 
*EL FILE
S, 
E,
COORD,
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
E,
COORD,
*NODE FILE
U,RF
*OUTPUT,FIELD
*NODE OUTPUT
U,RF
*END STEP