*HEADING
PINCHED CYL. WITH DIAPH. AT ENDS - COARSE IRREGULAR MESH - S4R,TYPE 2
*PREPRINT,ECHO=YES,MODEL=NO,HISTORY=NO
*NODE
*NGEN,NSET=BND1A
*NGEN,NSET=BND1B
*NSET,NSET=BND1
BND1A,BND1B
*NGEN,NSET=BND2
*NGEN,NSET=BND3
*NGEN,NSET=BND4A
*NGEN,NSET=BND4B
*NSET,NSET=BND4
BND4A,BND4B
*NGEN,NSET=NODES1
*NGEN,NSET=NODES2
*NFILL,NSET=ALL
BND2,NODES1,8
NODES2,BND3,8,100
NODES1,BND4A,8,2
*NMAP,TYPE=CYLINDRICAL,NSET=ALL
*ELEMENT,TYPE=S4R
*ELGEN,ELSET=EALL
*ELEMENT,TYPE=S4R
*ELGEN,ELSET=EALL
*ELSET,ELSET=PRINT
*TRANSFORM,TYPE=C,NSET=ALL
*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,1809,1607
 TIE,2009,1605
 TIE,2209,1603
 TIE,2409,1601
*STEP
*STATIC
*CLOAD
25,1,-0.25
*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