*HEADING
PINCHED CYL. WITH DIAPH. AT ENDS, COARSE IRREG. MESH - SC6R - TYPE 2   
*PREPRINT,ECHO=YES,MODEL=NO,HISTORY=NO
*NODE
1,298.5,0.,0.,1.,0.,0.
9,298.5,78.54,0.,1.,0.,.0
25,298.5,90.,0.,1.,0.,0.
1601,298.5,0.,300.,1.,0.,0.
1609,298.5,78.54,60.,1.,0.,0.
1625,298.5,90.,60.,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,78.54,0.,1.,0.,.0
10025,301.5,90.,0.,1.,0.,0.
11601,301.5,0.,300.,1.,0.,0.
11609,301.5,78.54,60.,1.,0.,0.
11625,301.5,90.,60.,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
*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