*heading
spherical indentation problem - axisymmetric - finite sliding
********************
**
** blank definition
**
*preprint, contact=yes
*node
1,0.0,2.0
10,6.0,2.0
51,0.0,0.0
60,6.0,0.0
*ngen,nset=top
1,10,1
*ngen,nset=bottom
51,60,1
*nfill,nset=nplate
top,bottom,5,10
*nset,nset=left,gen
1,51,10
*element,type=cax4h
1,1,11,12,2
*elgen,elset=eplate
1,9,1,1,5,10,9
*elset,elset=fplate,gen
1,4,1
*solidsection,elset=eplate,material=rubber
*material,name=rubber
*hyperelastic,n=1
80.0,20.0
********************
**
** sphere definition
**
*node
100,0.0,2.05
105,2.0,4.05
106,0.0,3.05
111,1.0,4.05
200,0.0,4.05
*ngen,line=c,nset=outer
100,105,1,200
*ngen,line=c,nset=inner
106,111,1,200
*element,type=cax4
100,100,101,107,106
*elgen,elset=sphere
100,5,1,1
*elset,elset=fsphere,gen
100,103,1
*element,type=cax3,elset=sphere
105,106,107,200
106,107,108,200
107,108,109,200
108,109,110,200
109,110,111,200
*solidsection,elset=sphere,material=steel
*material,name=steel
*elastic
30.e6,.3
********************
**
** fill in contact definition here
**
*elset,elset=sset, gen
100, 104, 1
*elset,elset=eset, gen
1,5,1
*surface, name=ssphere, trim=yes
sset,
*surface, name=splate, trim=yes
eset,
*surface interaction, name=press
2.,
*nset,nset=nodeprnt
bottom, 105,111,200
*boundary
200,2,2
105, 2,2
111, 2,2
100,1,1
106,1,1
200,1,1
bottom,2,2
left,1,1
********************
**
** step definition
**
*CONTACT INITIALIZATION DATA, NAME=SHRINK_FIT, SEARCH ABOVE=0.1
*CONTACT
*CONTACT INCLUSIONS, ALL EXTERIOR
*CONTACT PROPERTY ASSIGNMENT
ssphere,splate,press
*CONTACT INITIALIZATION ASSIGNMENT
ssphere, splate, SHRINK_FIT
*contact formulation, type=master slave roles
ssphere, splate, slave
*step,nlgeom,unsymm=yes
*static
.2,1.0
*boundary
200, 2,2,-1.0
105, 2,2,-1.0
111, 2,2,-1.0
*print, freq=0
*node print, freq=0
*el print, freq=0
*output,field,freq=1,var=preselect
*output,history,freq=1,var=preselect
*end step