*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.
*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
1000, .30,-.01,-.07
1001, .30,-.01,-.11
1010, .34,-.01,-.07
1019, .30,-.01,-.03
1028, .26,-.01,-.07
1037, .30,-.01,-.11
1100, .30, .06,-.07
1101, .30, .06,-.11
1110, .34, .06,-.07
1119, .30, .06,-.03
1128, .26, .06,-.07
1137, .30, .06,-.11
*NGEN,NSET=NR1,LINE=C
1001,1010,1,1000
1010,1019,1,1000
1019,1028,1,1000
1028,1037,1,1000
1101,1110,1,1100
1110,1119,1,1100
1119,1128,1,1100
1128,1137,1,1100
*NCOPY,CHANGE NUMBER=100,OLD SET=NR1,SHIFT
0.,0.01,0.

*NCOPY,CHANGE NUMBER=200,OLD SET=NR1,SHIFT
0.,0.02,0.

*NCOPY,CHANGE NUMBER=300,OLD SET=NR1,SHIFT
0.,0.03,0.

*NCOPY,CHANGE NUMBER=400,OLD SET=NR1,SHIFT
0.,0.04,0.

*NCOPY,CHANGE NUMBER=500,OLD SET=NR1,SHIFT
0.,0.05,0.

*NCOPY,CHANGE NUMBER=600,OLD SET=NR1,SHIFT
0.,0.06,0.

*NCOPY,CHANGE NUMBER=700,OLD SET=NR1,SHIFT
0.,0.07,0.

*ELEMENT,TYPE=R3D4,ELSET=CYLINDER
1001, 1001,1101,1102,1002
1036, 1036,1136,1101,1001
*ELGEN,ELSET=CYLINDER
1001, 35,1,1, 7,100,100
1036,  7,100,100
*BOUNDARY
999,1,6
*SURFACE,TYPE=ELEMENT,NAME=BOTTOM
BEAM,SNEG
*SURFACE,TYPE=ELEMENT,NAME=OUTSIDE
CYLINDER,SPOS
*RIGID BODY,ELSET=CYLINDER,REF NODE=999
*STEP
*DYNAMIC,EXPLICIT
,1.0E-02
*DSLOAD
BOTTOM,P,-1.E5
*CONTACT PAIR,MECHANICAL CONSTRAINT=PENALTY
OUTSIDE,BOTTOM
**
** 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
21,40,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_TESTNH
U,
***
*NODE OUTPUT
RF, 
*OUTPUT,HISTORY,TIME INTERVAL=.0025
*ENERGY OUTPUT,VAR=PRESELECT
*END STEP