*HEADING
 CANTILEVER BEAM MODELED USING SHELL ELEMENTS
 ELASTIC/PLASTIC CASE WITH PENALTY CONTACT
 CYLINDRICAL REFINEMENT OF RIGID SURFACE BECAUSE IT IS
   PARTIALLY WEIGHTED AS A SLAVE SURFACE.
*Preprint, Contact=Yes, Model=Yes, History=Yes
*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
*NGEN,NSET=NR1,LINE=C
*NCOPY,CHANGE NUMBER=100,OLD SET=NR1,SHIFT
*NCOPY,CHANGE NUMBER=200,OLD SET=NR1,SHIFT
*NCOPY,CHANGE NUMBER=300,OLD SET=NR1,SHIFT
*NCOPY,CHANGE NUMBER=400,OLD SET=NR1,SHIFT
*NCOPY,CHANGE NUMBER=500,OLD SET=NR1,SHIFT
*NCOPY,CHANGE NUMBER=600,OLD SET=NR1,SHIFT
*NCOPY,CHANGE NUMBER=700,OLD SET=NR1,SHIFT
*ELEMENT,TYPE=R3D4,ELSET=CYLINDER
*ELGEN,ELSET=CYLINDER
*BOUNDARY
999,1,6
*RIGID BODY,ELSET=CYLINDER,REF NODE=999
*SURFACE,TYPE=ELEMENT,NAME=BOTTOM
BEAM,SNEG
*************************************************************
**General contact definition
*Contact
*Contact Inclusions, AllElementBased
*************************************************************
*STEP
*DYNAMIC,EXPLICIT
,1.0E-02
*DSLOAD
BOTTOM,P,-1.E5
** 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