*HEADING
: PURE BENDING OF A THICK WALLED CYLINDER  [CAXA41]
** 6-IN LONG CYLINDER OF 2-IN & 6-IN INNER AND OUTER RADII  (8X12 MESH) 
** RIGID BODY MODE IN THE 1 DIRECTION IS ELIMINATED BY THE FOLLOWING EQUATION
*EQUATION
2, 
1,1,1.,1001,1,-1.
*NODE
1,2.,0.
9,6.,0.
121,2.,6.
129,6.,6.
*NGEN,NSET=BOT
1,9
*NGEN,NSET=TOP
121,129
*NFILL,NSET=P000
BOT,TOP,12,10
*ELEMENT,TYPE=CAXA41,OFFSET=1000,ELSET=ALL
1,1,2,12,11
*ELGEN,ELSET=ALL
1,8,1,1,12,10,10
*SOLID SECTION,ELSET=ALL,MATERIAL=STEEL
*MATERIAL,NAME=STEEL
*ELASTIC
30.E6,.33
*EXPANSION
0.0001, 
*NSET,NSET=BOTTOM,GEN
1,9
1001,1009
*NSET,NSET=L1,GEN
1,121,10
1001,1121,10
*NSET,NSET=L2,GEN
5,125,10
1005,1125,10
*NSET,NSET=L3,GEN
9,129,10
1009,1129,10
*ELSET,ELSET=END,GEN
111,118            
*BOUNDARY
BOTTOM,2,2
*NSET,NSET=NFILE
1,9,121,129, 1001,1009,1121,1129
*ELSET,ELSET=EFILE
1,8,111,118
*STEP
*STATIC
*ENERGY PRINT
*ENERGY FILE
*DLOAD
END,P3NU
*EL PRINT,ELSET=EFILE,POSITION=AVERAGED AT NODES
S, 
E, 
*NODE PRINT,NSET=NFILE,TOTAL=NO
COORD,U1,U2,RF1,RF2
** *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,POSITION=AVERAGED AT NODES
S, 
*NODE FILE,NSET=NFILE
U, 
*END STEP