*heading
 First order element test with elements falling off
*node
1,-6,0
25,6,0
1001,-6,5
1025,6,5
2001,-2,5
2009,2,5
2801,-2,9
2809,2,9
*node,nset=top
5000,0,9
*ngen,nset=basebot
1,25
*ngen,nset=basetop
1001,1025
*ngen,nset=sliderbot
2001,2009
*ngen,nset=slidertop
2801,2809
*nfill,nset=base
basebot,basetop,10,100
*nfill,nset=slider
sliderbot,slidertop,8,100
*element,type=cpe8h
1,1,3,203,201,2,103,202,101
1001,2001,2003,2203,2201,2002,2103,2202,2101
*elgen,elset=base
1,12,2,1,5,200,100
*elgen,elset=slider
1001,4,2,1,4,200,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
linear,2002,2001,2003
linear,2004,2003,2005
linear,2006,2005,2007
linear,2008,2007,2009
*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,freq=10
*element output 
s,e
*node output 
u,
*end step
*step,unsymm=yes
*static
1.,72.,,1.
*boundary
top,1,1,7.3
*end step
*step,inc=200
*static
1.,144,,1.
*boundary
top,1,1,-7.3
*end step
*step
*static
1.,72.,,1.
*boundary
top,1,1,0.
*end step