*HEADING
PINCHED CYL. WITH DIAPH. AT ENDS, FINE IRREG. MESH - STRI65 - TYPE 1    
*PREPRINT,ECHO=YES,MODEL=NO,HISTORY=NO
*NODE ,NSET=P1
1,300.,0.,0.,1.,0.,0.
9,300.,63.2616,0.,1.,0.,.0
25,300.,90.,0.,1.,0.,0.
1601,300.,0.,300.,1.,0.,0.
1609,300., 63.2616,140.,1.,0.,0.
1625,300.,90.,140.,1.,0.,0.
2409,300.,0.,300.,1.,0.,0.
2425,300.,90.,300.,1.,0.,0.
*NGEN,NSET=BND1A
1,9
*NGEN,NSET=BND1B
9,25
*NSET,NSET=BND1
BND1A,BND1B
*NGEN,NSET=BND2
1,1601,100
*NGEN,NSET=BND3
2409,2425
*NGEN,NSET=BND4A
25,1625,100
*NGEN,NSET=BND4B
1625,2425,100
*NSET,NSET=BND4
BND4A,BND4B
*NGEN,NSET=NODES1
9,1609,100
*NGEN,NSET=NODES2
1609,1625
*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=STRI65
    1,     1,     3,   203,     2,   103,   102
50001,     1,   203,   201,   102,   202,   101
*ELGEN,ELSET=EALL
    1,    12,     2,     1,     8,   200,   100
50001,    12,     2,     1,     8,   200,   100
*ELEMENT,TYPE=STRI65
 1001,  1609,  1611,  1811,  1610,  1711,  1710
51001,  1609,  1811,  1809,  1710,  1810,  1709
*ELGEN,ELSET=EALL
 1001,     8,     2,     1,     4,   200,   100
51001,     8,     2,     1,     4,   200,   100
*ELSET,ELSET=PRINT
 1, 
*TRANSFORM,TYPE=C,NSET=ALL
.0,.0,.0,300,.0,.0
*SHELL SECTION,ELSET=EALL,MATERIAL=MAT
3.0, 
*MATERIAL,NAME=MAT
*ELASTIC
3.0E6,0.3
*BOUNDARY
BND2,YSYMM
BND1,ZSYMM
BND4,YSYMM
BND3,2
BND3,1
*MPC
 TIE,1709,1608
 TIE,1809,1607
 TIE,1909,1606
 TIE,2009,1605
 TIE,2109,1604
 TIE,2209,1603
 TIE,2309,1602
 TIE,2409,1601
*STEP
*STATIC
*CLOAD
25,1,-0.25
*EL PRINT,ELSET=PRINT
*EL FILE,ELSET=PRINT
SF,
*OUTPUT,FIELD
*ELEMENT OUTPUT,ELSET=PRINT
SF,
*NODE FILE,NSET=P1
U,RF
*OUTPUT,FIELD
*NODE OUTPUT,NSET=P1
U,RF
*END STEP