*heading
 test contact controls with free block
*preprint,model=yes,history=yes
*node
1,0,.01
2,1,.01
3,1,1
4,0,1
11,0,-.99
12,1,-.99
13,1,-.01
14,0,-.01
21,0,-2
22,1,-2
23,1,-1.01
24,0,-1.01
*nset,nset=ntop1
3,4
*nset,nset=nbot2
21,22
*element, type=cpe4, elset=top
1,1,2,3,4
*element, type=cpe4, elset=mid
11,11,12,13,14
*element, type=cpe4, elset=bot
21,21,22,23,24
*elset, elset=eall
top,mid,bot
*solid section, elset=eall, material=mat
*material, name=mat
*elastic
1000.,
*surface,name=stop1
top,s1
*surface,name=sbot1
mid,s3
*surface,name=stop2
mid,s1
*surface,name=sbot2
bot,s3
*contact pair, interaction=cp1, small sliding
sbot1,stop1
*contact pair, interaction=cp2, small sliding
sbot2,stop2
*surface interaction, name=cp1
*surface interaction, name=cp2
*surface behavior, no separation
*friction
0.0002
*boundary
nbot2,encastre
*step
close the gaps
*static
.1,1.,.1,.1
*boundary
ntop1,1,1
ntop1,2,2,-.05
*contact controls, stabilize
,1.
*contact controls, stabilize=10., slave=sbot2, master=stop2
,.1
*print, contact=yes
*energy print
*contact print
cstress, cdstress, cdisp
*node print,summ=yes
u,rf
*output, field
*element output
s,e
*node output
u,rf
*contact output
cstress,cdstress,cdisp
*end step
*step,unsymm=yes
slide
*static
.1,1.,.1,.1
*boundary
ntop1,1,1,.05
ntop1,2,2,-.05
*contact controls, stabilize
*end step
*step
open
*static
.1,1.,.1,.1
*boundary
ntop1,1,1,.05
ntop1,2,2,.04
*contact controls, stabilize
0.,0.,0.072
*end step