*HEADING
 PIPE CRUSHING BENCHMARK [SC8R]
*PREPRINT, MODEL=YES
*NODE
*NGEN,NSET=YSYMM
*NGEN,NSET=XSYMM
*NGEN,LINE=C, NSET=NALL
*ELEMENT,TYPE=SC8R
*ELGEN,ELSET=TUBE
*SHELL SECTION,ELSET=TUBE,MATERIAL=A1, POISSON=ELASTIC
 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
*ELSET,ELSET=ETOP,GENERATE
**
** 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
*NGEN, NSET=PLATESIDE1
*NGEN, NSET=PLATESIDE2
*NFILL, NSET=PLATE
 PLATESIDE1,PLATESIDE2,10,1
*ELEMENT, TYPE=SC8R
*ELGEN, ELSET=PLATE
*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