*HEADING
 THICK CYLINDER - INTERNAL PRESSURE    (CAXA81)  
*NODE
1,10.
11,20.
101,10.,1.
111,20.,1.
201,10.,2.
211,20.,2.
*NGEN,NSET=SIDE1
1,11
201,211
*NGEN
101,111
*NSET,NSET=BORE1
1,101,201
*ELEMENT,TYPE=CAXA81,OFFSET=1000
1,1,3,203,201,2,103,202,101
*ELGEN,ELSET=EALL
1,5,2
*ELSET,ELSET=PRIN
3,
*NSET,NSET=SIDE,GENERATE
1,11
201,211
1001,1011
1201,1211
*NSET,NSET=BORE
1,101,201
1001,1101,1201
*SOLID SECTION,MATERIAL=A1,ELSET=EALL
*MATERIAL,NAME=A1
*ELASTIC
30.E6,.3
*PLASTIC
30000., 
** Axisymmetric constraints
*EQUATION
2
1001,1,1., 1,1,-1.
2
1101,1,1.,101,1,-1.
2
1201,1,1.,201,1,-1.
2
1011,1,1.,11,1,-1.
2
1111,1,1.,111,1,-1.
2
1211,1,1.,211,1,-1.
*BOUNDARY
SIDE,2
*STEP,NLGEOM,INC=50
*STATIC
.1,1.
*BOUNDARY
BORE,1,1,20.
*PRINT,RESIDUAL=NO
*NODE PRINT
U,RF
*EL PRINT,ELSET=PRIN,POSITION=CENTROIDAL,SUMMARY=NO
S,MISES,
E, 
*NODE FILE,NSET=BORE
U, 
RF, 
*EL FILE,ELSET=PRIN,POSITION=CENTROIDAL
S, 
** Selective output
*OUTPUT,HISTORY,FREQ=6
*NODE OUTPUT,NSET=SIDE
U1,U2,RF1,RF2
*OUTPUT,FIELD,FREQ=6
*ELEMENT OUTPUT,ELSET=EOUT
S,LE,SINV,PE
*ELSET,ELSET=EOUT
1,3,5
*END STEP