*HEADING
 CANTILEVER BEAM MODELED USING SHELL ELEMENTS --- ELASTIC/PLASTIC CASE WITH
 CONTACT
*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,ELSET=BEAM,MATERIAL=STEEL,POISSON=0.3
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
*BOUNDARY
999,1,6
*CONTACT PAIR,INTERACTION=CYL
BOTTOM,CYLINDER
*SURFACE INTERACTION,NAME=CYL
*SURFACE BEHAVIOR,AUGMENTED LAGRANGE
*RESTART,WRITE,FREQ=25
*NSET,NSET=NOUT
*RIGID BODY,ANALYTICAL SURFACE=CYLINDER,REFNODE=999
*SURFACE,TYPE=REVOLUTION,NAME=CYLINDER
START,0.04,0.07
LINE,0.04,0.0
*SURFACE,TYPE=ELEMENT,NAME=BOTTOM
BEAM,SNEG
*STEP,INC=2000,NLGEOM
*DYNAMIC,HAFTOL=1.0e3
2.5E-6,1.0E-02
*DLOAD
BEAM,P,-1.E5
*NODE PRINT,FREQ=0
*EL PRINT,FREQ=0
*OUTPUT, FIELD, FREQUENCY=10
*NODE OUTPUT
 U,V,RF
*OUTPUT, HISTORY, FREQUENCY=10
*ENERGY OUTPUT
*OUTPUT, FIELD, FREQUENCY=9999
*ELEMENT OUTPUT
 STH,SF,SE
 PE,LE
*END STEP