*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