*HEADING
 ANNULAR RING UNDER APPLIED TWIST ON INNER AND OUTER SURFACES 
 ELEMENT : CGAX8RT
 MPC: QUADRATIC, T LINEAR, AND TIE
*RESTART, WRITE
*NODE,NSET=NALL
1, 1.0, 0.0
19, 1.9, 0.0
20, 1.9, 0.0
22, 2.0, 0.0
2019, 1.9, 1.0
2020, 1.9, 1.0
2022, 2.0, 1.0
2001, 1.0, 1.0
*NGEN,NSET=BOTTOM
1,19
20,22
*NGEN,NSET=TOP
2001,2019
2020,2022
*NSET,NSET=INNER,GENERATE
1, 2001, 100
*NSET,NSET=OUTER,GENERATE
22,2022, 100
*NFILL, NSET=ALL
BOTTOM, TOP, 20, 100
*NSET,NSET=P1MPC,GENERATE
107,1707,400
*NSET,NSET=P2MPC,GENERATE
307,1907,400
*NSET,NSET=AMPC,GENERATE
7,1607,400
*NSET,NSET=BMPC,GENERATE
207,1807,400
*NSET,NSET=CMPC,GENERATE
407,2007,400
*NSET,NSET=QUADMPC
P1MPC,P2MPC,AMPC,BMPC,CMPC
*NSET,NSET=ATIE,GENERATE
19,2019,200
*NSET,NSET=BTIE,GENERATE
20,2020,200
*ELEMENT, TYPE=CGAX8RT, ELSET=EALL
1,1,3,203,201,2,103,202,101
4,7,9,409,407,8,209,408,207
10,20,22,422,420,21,222,421,220
*ELGEN, ELSET=EALL
1, 3, 2, 1, 10, 200, 10,
4, 6, 2, 1,  5, 400, 10,
10,1, 2, 1,  5, 400, 10,
*SOLID SECTION, ELSET=EALL, MATERIAL=MALL
*MATERIAL, NAME=MALL
*ELASTIC
1.E6, 0.3
*EXPANSION
0.0, 
*CONDUCTIVITY
1.0, 
*BOUNDARY
BOTTOM,1,2
*INITIAL CONDITIONS, TYPE=TEMPERATURE
NALL,20.0
*MPC
QUADRATIC,P1MPC,AMPC,BMPC,CMPC
QUADRATIC,P2MPC,AMPC,BMPC,CMPC
TLINEAR,BMPC,AMPC,CMPC
TIE,BTIE,ATIE
**
** STEP I
**
*STEP
*COUPLED TEMPERATURE-DISPLACEMENT, STEADY STATE
*BOUNDARY
INNER, 5,5, 0.01
OUTER, 5,5, -0.01
INNER, 11,11, 20.0
OUTER, 11,11, 40.0
*EL PRINT, ELSET=EALL
S, MISES
E,
COORD,E13,S13
*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, 
*NODE FILE
U,RF
*OUTPUT,FIELD
*NODE OUTPUT
U,RF
*EL FILE
S, 
E, 
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,
E,
*END STEP