*HEADING
: PURE BENDING OF A THICK WALLED CYLINDER  [CAXA4R3]
** 6-IN LONG CYLINDER OF 2-IN & 6-IN INNER AND OUTER RADII  (16X24 MESH) 
** RIGID BODY MODE IN THE 1 DIRECTION IS ELIMINATED BY THE FOLLOWING EQUATION
*EQUATION
2, 
1,1,1.,3001,1,-1.
*NODE
1,2.,0.
17,6.,0.
481,2.,6.
497,6.,6.
*NGEN,NSET=BOT
1,17
*NGEN,NSET=TOP
481,497
*NFILL,NSET=P000
BOT,TOP,24,20
*ELEMENT,TYPE=CAXA4R3,OFFSET=1000,ELSET=ALL
1,1,2,22,21
*ELGEN,ELSET=ALL
1,16,1,1,24,20,20
*SOLID SECTION,ELSET=ALL,MATERIAL=STEEL
*MATERIAL,NAME=STEEL
*ELASTIC
30.E6,.33
*EXPANSION
0.0001, 
*NSET,NSET=BOTTOM,GEN
1,17
1001,1017
2001,2017
3001,3017
*NSET,NSET=L1,GEN
1,481,20
1001,1481,20
2001,2481,20
3001,3481,20
*NSET,NSET=L2,GEN
9,489,20
1009,1489,20
2009,2489,20
3009,3489,20
*NSET,NSET=L3,GEN
17,497,20
1017,1497,20
2017,2497,20
3017,3497,20
*ELSET,ELSET=END,GEN
461,476        
*BOUNDARY
BOTTOM,2,2
*NSET,NSET=NFILE
1,17,481,497, 
1001,1017,1481,1497
2001,2017,2481,2497
3001,3017,3481,3497
*ELSET,ELSET=EFILE
1,16,461,476
*STEP
*STATIC
*ENERGY PRINT
*DLOAD
END,P3NU
*EL PRINT,ELSET=EFILE
S, 
E, 
*NODE PRINT,NSET=NFILE,TOTAL=NO
COORD,U1,U2,RF1,RF2
***EL PRINT,POSITION=AVERAGED AT NODES
**S
**E
***NODE PRINT,NSET=L1,TOTAL=NO
**COORD,U1,U2,RF1,RF2
***NODE PRINT,NSET=L2,TOTAL=NO
**COORD,U1,U2,RF1,RF2
***NODE PRINT,NSET=L3,TOTAL=NO
**COORD,U1,U2,RF1,RF2
*EL FILE,ELSET=EFILE
S, 
*NODE FILE,NSET=NFILE
U,RF
*END STEP