*HEADING
 CANTILEVER BEAM MODELED USING SHELL ELEMENTS 
 ELASTIC/PLASTIC CASE WITH CONTACT
*NODE
  1, 0., 0., 0.00125
 21,.50, 0., 0.00125
301, 0.,.05, 0.00125
321,.50,.05, 0.00125
*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,OFFSET=SPOS
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,-.10
1010, .33,-.01,-.07
1019, .30,-.01,-.04
1028, .27,-.01,-.07
1037, .30,-.01,-.10
1100, .30, .06,-.07
1101, .30, .06,-.10
1110, .33, .06,-.07
1119, .30, .06,-.04
1128, .27, .06,-.07
1137, .30, .06,-.10
*NGEN,NSET=RIGID,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
*ELEMENT,TYPE=R3D4,ELSET=CYLINDER
1001, 1001,1101,1102,1002
1036, 1036,1136,1101,1001
*ELGEN,ELSET=CYLINDER
1001, 35,1,1
*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
0.02, 
*STEP
*DYNAMIC,EXPLICIT
,1.0E-02
*DSLOAD
BOTTOM,P,-1.E5
*CONTACT PAIR
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
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