*HEADING
 PIPE CRUSHING BENCHMARK [SC8R]
*PREPRINT, MODEL=YES               
*NODE
  101, 1.901
  117, 0.   ,1.901
  501, 2.25
  517, 0.   ,2.25
 1101, 1.901,0    ,1
 1117, 0.   ,1.901,1
 1501, 2.25 ,0    ,1
 1517, 0.   ,2.25 ,1
*NGEN,NSET=YSYMM
  101, 501,100
 1101,1501,100
*NGEN,NSET=XSYMM
 117, 517,100
1117,1517,100
*NGEN,LINE=C, NSET=NALL
  101, 117,1,,0.,0.,0.
 1101,1117,1,,0.,0.,1.
  201, 217,1,,0.,0.,0.
 1201,1217,1,,0.,0.,1.
  301, 317,1,,0.,0.,0.
 1301,1317,1,,0.,0.,1.
  401, 417,1,,0.,0.,0.
 1401,1417,1,,0.,0.,1.
  501, 517,1,,0.,0.,0.
 1501,1517,1,,0.,0.,1.
*ELEMENT,TYPE=SC8R
 101,101,201,202,102,1101,1201,1202,1102
*ELGEN,ELSET=TUBE
 101,4,100,1,16,1,1000
*SHELL SECTION,ELSET=TUBE,MATERIAL=A1, POISSON=ELASTIC, STACK=1
 1.0,
*MATERIAL,NAME=A1
*ELASTIC
 27.E6,.3
*PLASTIC
 35000.,0.
 40000.,.0035
 43750.,.0083
 46250.,.0133
 50000.,.0182
 52500.,.0281
 55000.,.0380
 215000.,1.000
*NODE,NSET=RSNODE
 10000,5.,2.25
*ELSET,ELSET=ETOP,GENERATE
 104,15104,1000
** 
** Define a rigid analytical surface
** *rigid body,analytical surface=BSURF,REFNODE=10000
** *surface,TYPE=CYLINDER,NAME=BSURF
** 0.0,2.25,0.0,5.0,2.25,0.0
** 0.0,2.25,1.0
** START,0.0,0.0
** LINE,5.0,0.0
**
*NODE
  2000,-1.0,2.25,-1.0
  3000, 5.0,2.25,-1.0
  3010, 5.0,2.25, 1.0
  2010,-1.0,2.25, 1.0
 12000,-1.0,2.30,-1.0
 13000, 5.0,2.30,-1.0
 13010, 5.0,2.30, 1.0
 12010,-1.0,2.30, 1.0
*NGEN, NSET=PLATESIDE1
  2000, 3000,100
 12000,13000,100
*NGEN, NSET=PLATESIDE2
  2010, 3010,100
 12010,13010,100
*NFILL, NSET=PLATE
 PLATESIDE1,PLATESIDE2,10,1
*ELEMENT, TYPE=SC8R
 20000,2000,2001,2101,2100,12000,12001,12101,12100
*ELGEN, ELSET=PLATE
 20000,10,100,1,10,1,10
*SHELL SECTION,ELSET=PLATE,MATERIAL=A1
 0.05
*SURFACE,NAME=BSURF,TRIM=YES
 PLATE, S1
*RIGID BODY,ELSET=PLATE,REFNODE=10000
*SURFACE,NAME=ASURF,TRIM=YES
 ETOP
*CONTACT PAIR,INTERACTION=SMOOTH
 ASURF,BSURF
*SURFACE INTERACTION,NAME=SMOOTH
*BOUNDARY
 XSYMM,1,1
 YSYMM,2,2
 NALL, 3,3
*RESTART,WRITE,FREQUENCY=25
*STEP,NLGEOM,INC=200
*STATIC
 0.002 , 1.0   ,0.0001 , 0.01
*BOUNDARY
 RSNODE,1,1,0
 RSNODE,2,2,-1.3
 RSNODE,3,6,0
*MONITOR,NODE=10000,DOF=2,FREQUENCY=2
*OUTPUT,FIELD,FREQ=25
*CONTACT OUTPUT,VARIABLE=PRESELECT,SLAVE=ASURF
*OUTPUT,FIELD,FREQ=5
*NODE OUTPUT
 U,
*ELEMENT OUTPUT
 PEEQ,
 SF,
 STH,
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=RSNODE
 U, 
 RF, 
*END STEP