*heading First order element test with elements falling off *preprint, model=yes,history=yes,contact=yes *node 1,-6,0 13,6,0 501,-6,5 513,6,5 1001,-2,5 1005,2,5 1401,-2,9 1405,2,9 *node,nset=top 2000,0,9 *ngen,nset=basebot 1,13 *ngen,nset=basetop 501,513 *ngen,nset=sliderbot 1001,1005 *ngen,nset=slidertop 1401,1405 *nfill,nset=base basebot,basetop,5,100 *nfill,nset=slider sliderbot,slidertop,4,100 *element,type=cpe4h 1,1,2,102,101 1001,1001,1002,1102,1101 *elgen,elset=base 1,12,1,1,5,100,100 *elgen,elset=slider 1001,4,1,1,4,100,100 *elset,elset=basetop,generate 401,412 *elset,elset=sliderbot,generate 1001,1004 *nset,nset=sliderbot,generate 1001,1005 *nset,nset=basetop,generate 501,513 *solid section,elset=base,material=stiff *solid section,elset=slider,material=stiff *material,name=stiff *hyperelastic,neo 100,.01 *surface,name=basetop,property=WEAR basetop,s3 *surface,name=sliderbot,property=WEAR sliderbot,s1 *contact pair,type=surface,interaction=contact sliderbot,basetop *surface interaction,name=contact *mpc beam,slidertop,top *boundary basebot,encastre *surfaceproperty, name=WEAR *wearsurfaceproperties, friccoefdependent=NO, unitlesswearcoef=NO 2.5e-10 *AMPLITUDE,NAME=AMP1,DEFINITION=TABULAR,TIME=STEPTIME 0.0,0.0 50.0,1.0 100.0,0.0 *AMPLITUDE,NAME=AMP2,DEFINITION=TABULAR,TIME=STEPTIME 0.0,0.0 50.0,-1.0 100.0,0.0 *step, nlgeom *static 1,1 *boundary top,1,1 top,6,6 *cload top,2,-100. *node print,nset=top u,rf *output,field,variable=PRESELECT *contact output cwear *end step ** ** *STEP CYCLING,START 10,1 10000,25000 ** *step,unsymm=yes,inc=1000 *static 2.,100.,,2. *boundary,AMPLITUDE=AMP1 top,1,1,3.5 *output,field,variable=PRESELECT *contact output cwear *print,contact=yes *end step ** *step,unsymm=yes,inc=1000 *static 2.,100.,,2. *STEPCYCLINGCONTROL,CRITERION=WEAR-BASED,ACTION=END CYCLING,NAME=TEST sliderbot,ABSMAX,0.012,0.006 basetop,ABSMAX,0.012,0.006 *boundary,AMPLITUDE=AMP2 top,1,1,3.5 *output,field,variable=PRESELECT *contact output cwear *print,contact=yes *end step *STEP CYCLING,END