*heading
 test contact controls with free block
*preprint,model=yes,history=yes
*node
1,0,0,.01
2,1,0,.01
3,1,1,.01
4,0,1,.01
5,0,0,1
6,1,0,1
7,1,1,1
8,0,1,1
11,0,0,-.99
12,1,0,-.99
13,1,1,-.99
14,0,1,-.99
15,0,0,-.01
16,1,0,-.01
17,1,1,-.01
18,0,1,-.01
21,0,0,-2
22,1,0,-2
23,1,1,-2
24,0,1,-2
25,0,0,-1.01
26,1,0,-1.01
27,1,1,-1.01
28,0,1,-1.01
*nset,nset=ntop1
5,6,7,8
*nset,nset=nbot1
11,12,13,14
*nset,nset=ntop2
15,16,17,18
*nset,nset=nbot2
21,22,23,24
*element, type=c3d8, elset=top
1,1,2,3,4,5,6,7,8
*element, type=c3d8, elset=mid
11,11,12,13,14,15,16,17,18
*element, type=c3d8, elset=bot
21,21,22,23,24,25,26,27,28
*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,s2
*surface,name=stop2
mid,s1
*surface,name=sbot2
bot,s2
*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
ntop1,3,3,-.05
*contact controls, stabilize
,1.
*contact controls, stabilize=10., slave=sbot2, master=stop2
,.1
*print, contact=yes
*energy print
*contact print
cpress,cshear1,cshear2,cdpress,cdshear1,cdshear2
*node print,summ=yes
u,rf
*output, field
*element output
s,e
*node output
u,rf
*contact output
cstress,cdstress,cdisp
*contact file
cstress,cdstress,cdisp
*end step
*step,unsymm=yes
slide
*static
.1,1.,.1,.1
*boundary
ntop1,1,1,.05
ntop1,2,2
ntop1,3,3,-.05
*contact controls, stabilize
*end step
*step
open
*static
.1,1.,.1,.1
*boundary
ntop1,1,1,.05
ntop1,2,2
ntop1,3,3,.04
*contact controls, stabilize
 ,0.,0.072
*end step