*HEADING
 CANTILEVER BEAM MODELED USING SHELL ELEMENTS 
 ELASTIC/PLASTIC CASE WITH CONTACT ON
 ANALYTICAL RIGID SURFACE 
*NODE
  1, 0., 0., 0.
 21,.50, 0., 0.
301, 0.,.05, 0.
321,.50,.05, 0.
*NGEN,NSET=FREE
1,21,1
*NGEN,NSET=CENTER
301,321,1
*NFILL,NSET=ALLN
FREE,CENTER,3,100
*NSET,NSET=CLAMP,GEN
1,301,100
*ELEMENT,TYPE=S4R
1, 1,2,102,101
*ELGEN,ELSET=BEAM
1, 20,1,1, 3,100,20
*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
 999, .30, .035,-.07
*BOUNDARY
999,1,6
*SURFACE, NAME=CYLINDER, TYPE=CYLINDER
0.3,-0.01,-0.07, 0.4,-0.01,-0.07
0.3,0.06,-0.07
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
 21,321,100
***NSET,NSET=QA_TESTCLAMP
**CLAMP,
*ELSET,ELSET=QA_TESTEL,GEN
1,20,1
*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