*HEADING
 CANTILEVER BEAM MODELED USING SHELL ELEMENTS
 ELASTIC/PLASTIC CASE WITH CONTACT ON
 ANALYTICAL RIGID SURFACE
*NODE
*NGEN,NSET=FREE
*NGEN,NSET=CENTER
*NFILL,NSET=ALLN
FREE,CENTER,3,100
*NSET,NSET=CLAMP,GEN
*ELEMENT,TYPE=S4R
*ELGEN,ELSET=BEAM
*SHELL SECTION, SECTION INT=SIMP,ELSET=BEAM,MATERIAL=STEEL
0.0025,
*MATERIAL,NAME=STEEL
*DENSITY
7800.,
*ELASTIC
200.E9,.3
*PLASTIC
250.E6,
*BOUNDARY
CLAMP,1,1
CLAMP,3,6
301,2,2
CENTER,YSYMM
*NODE
*BOUNDARY
999,1,6
*SURFACE, NAME=CYLINDER, TYPE=CYLINDER
START,0.,0.04
CIRCL,0.04,0.,0.,0.
CIRCL,0.0,-0.04
CIRCL,-0.04,0.
CIRCL,0.,0.04
*SURFACE,TYPE=ELEMENT,NAME=BOTTOM
BEAM,SNEG
*RIGID BODY, REF NODE=999, ANALYTICAL SURFACE =CYLINDER
*STEP
*DYNAMIC,EXPLICIT
,1.0E-02
*DSLOAD
BOTTOM,P,-1.E5
*CONTACT PAIR
CYLINDER,BOTTOM
*PRINT,ALLKE=YES,CRITICAL ELEMENT=YES
**
** Output requests for postprocessing
**
*OUTPUT,FIELD,VAR=PRESELECT
**
** Output requests for qa testing
**
*NSET,NSET=QA_TESTNH,GEN
*NSET,NSET=QA_TESTCLAMP
CLAMP,
*ELSET,ELSET=QA_TESTEL,GEN
*NSET,NSET=QA_TESTN,ELSET=QA_TESTEL
*OUTPUT,FIELD,NUM=1
*ELEMENT OUTPUT,ELSET=QA_TESTEL
PEEQ,
*NODE OUTPUT,NSET=QA_TESTN
U,
*OUTPUT,FIELD,NUM=4,TIMEMARKS=YES
*NODE OUTPUT,NSET=QA_TESTCLAMP
RF,
*NODE OUTPUT,NSET=QA_TESTNH
U,
*OUTPUT,HISTORY,TIME INTERVAL=.0025
*ENERGY OUTPUT,VAR=PRESELECT
*END STEP