*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
*solid section,elset=base,material=stiff
*solid section,elset=slider,material=stiff
*material,name=stiff
*hyperelastic,neo
100,.01
*contact
*contact inclusion, all exterior
*mpc
beam,slidertop,top
*boundary
basebot,encastre
*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
*end step
*step,unsymm=yes,inc=1000
*static
1.,375.,,1.
*boundary
top,1,1,7.45
*output,field,variable=PRESELECT
*print,contact=yes
*end step