*HEADING
 PINCHED CYL. WITH DIAPH. AT ENDS - COARSE IRREGULAR MESH - SC6R,TYPE 1  
*PREPRINT,ECHO=YES,MODEL=NO,HISTORY=NO
*NODE
    1,298.5,0.,0., 1.,0.,0.
    9,298.5,63.2616,0., 1.,0.,0.
   25,298.5,90.,0., 1.,0.,0.
 1601,298.5,0.,300., 1.,0.,0.
 1609,298.5,63.2616,140., 1.,0.,0.
 1625,298.5,90.,140., 1.,0.,0.
 2409,298.5,0.,300., 1.,0.,0.
 2425,298.5,90.,300., 1.,0.,0.
10001,301.5,0.,0., 1.,0.,0.
10009,301.5,63.2616,0., 1.,0.,0.
10025,301.5,90.,0., 1.,0.,0.
11601,301.5,0.,300., 1.,0.,0.
11609,301.5,63.2616,140., 1.,0.,0.
11625,301.5,90.,140., 1.,0.,0.
12409,301.5,0.,300., 1.,0.,0.
12425,301.5,90.,300., 1.,0.,0.
*NGEN,NSET=BND1A
 1,9
 10001,10009
*NGEN,NSET=BND1B
 9,25
 10009,10025
*NSET,NSET=BND1
 BND1A,BND1B
*NGEN,NSET=BND2
 1,1601,100
 10001,11601,100
*NGEN,NSET=BND3
 2409,2425
 12409, 12425
*NGEN,NSET=BND4A
 25,1625,100
 10025,11625,100
*NGEN,NSET=BND4B
 1625,2425,100
 11625, 12425,100
*NSET,NSET=BND4
 BND4A,BND4B
*NGEN,NSET=NODES1
 9,1609,100
 10009,11609,100
*NGEN,NSET=NODES2
 1609,1625
 11609, 11625
*NFILL,NSET=ALL
 BND2,NODES1,8
 NODES2,BND3,8,100
 NODES1,BND4A,16
*NMAP,TYPE=CYLINDRICAL,NSET=ALL
 0.,0.,0., 300.,0.,0.
 0.,300.,0.

*ELEMENT,TYPE=SC6R
 1,1,3,203,10001,10003,10203
 2,1,203,201,10001,10203,10201
*ELGEN,ELSET=EALL
 1,12,2,2,8,200,100
 2,12,2,2,8,200,100
*ELEMENT,TYPE=SC6R
 1001,1609,1611,1811,11609,11611,11811
 1002,1609,1811,1809,11609,11811,11809
*ELGEN,ELSET=EALL
 1001,8,2,2,4,200,100
 1002,8,2,2,4,200,100
*ELSET,ELSET=PRINT
 1, 
*TRANSFORM,TYPE=C,NSET=ALL
.0,.0,.0,300,.0,.0
*SHELL SECTION,ELSET=EALL,MATERIAL=MAT, POISSON=ELASTIC
 3.0, 
*MATERIAL,NAME=MAT
*ELASTIC
 3.0E6,0.3
*BOUNDARY
 BND2,YSYMM
 BND1,ZSYMM
 BND4,YSYMM
 BND3,2
 BND3,1
*MPC
 TIE,1809,1607
 TIE,2009,1605
 TIE,2209,1603
 TIE,2409,1601
 TIE,11809,11607
 TIE,12009,11605
 TIE,12209,11603
 TIE,12409,11601
*STEP
*STATIC
*CLOAD
 25,1,-0.125
 10025,1,-0.125
*EL PRINT,ELSET=PRINT
*ENERGY PRINT
*EL FILE,ELSET=PRINT
 S,E
*OUTPUT,FIELD
*ELEMENT OUTPUT,ELSET=PRINT
 S,E
*NODE FILE
 U,RF
*OUTPUT,FIELD
*NODE OUTPUT
 U,RF
*END STEP