*HEADING
: PURE BENDING OF A THICK WALLED CYLINDER  [CAXA83]
** 6-IN LONG CYLINDER OF 2-IN & 6-IN INNER AND OUTER RADII  (1X1 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.
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=CAXA83,OFFSET=1000,ELSET=ALL
1,1,9,129,121,5,69,125,61
*SOLID SECTION,ELSET=ALL,MATERIAL=STEEL
*MATERIAL,NAME=STEEL
*ELASTIC
30.E6,.33
*EXPANSION
0.0001, 
*NSET,NSET=BOTTOM,GEN
1,9
1001,1009
2001,2009
3001,3009
*NSET,NSET=L1,GEN
1,121,10
1001,1121,10
2001,2121,10
3001,3121,10
*NSET,NSET=L2,GEN
5,125,10
1005,1125,10
2005,2125,10
3005,3125,10
*NSET,NSET=L3,GEN
9,129,10
1009,1129,10
2009,2129,10
3009,3129,10
*ELSET,ELSET=END
1    ,      ,    
*BOUNDARY
BOTTOM,2,2
*NSET,NSET=NFILE
1,9,121,129, 3001,3009,3121,3129
*ELSET,ELSET=EFILE
1,8,111,118
*STEP
*STATIC
*ENERGY PRINT
*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