*heading
Shrink fit of matched rings
*preprint,model=yes
*node, nset=base1
1,10,0,0
301,12.01,0,0
401,10,0,5
701,12.01,0,5
50,-10,0,0
350,-12.01,0,0
450,-10,0,5
750,-12.01,0,5
99,10,0,0
399,12.01,0,0
499,10,0,5
799,12.01,0,5
*ngen, line=c, nset=in_in
1,50,1,,0,0,0,0,0,1
50,99,1,,0,0,0,0,0,1
401,450,1,,0,0,5,0,0,1
450,499,1,,0,0,5,0,0,1
*ngen, line=c, nset=in_out
301,350,1,,0,0,0,0,0,1
350,399,1,,0,0,0,0,0,1
701,750,1,,0,0,5,0,0,1
750,799,1,,0,0,5,0,0,1
*nfill, nset=in_nodes
in_in,in_out,3,100
*element, type=c3d8
1,1,101,102,2,401,501,502,402
98,98,198,101,1,498,598,501,401
*elgen, elset=in_elems
1,97,1,1,3,100,100
*elgen, elset=in_elems
98,3,100,100
*solid section, elset=in_elems, material=elas
*material, name=elas
*elastic
1.e6,.3
*elset,elset=in_out,generate
201,298
*surface, name=in_out
in_out,s4
*node, nset=base2
1001,12,0,0
1301,14,0,0
1401,12,0,5
1701,14,0,5
1019,-12,0,0
1319,-14,0,0
1419,-12,0,5
1719,-14,0,5
1037,12,0,0
1337,14,0,0
1437,12,0,5
1737,14,0,5
*ngen, line=c, nset=out_in
1001,1019,1,,0,0,0,0,0,1
1019,1037,1,,0,0,0,0,0,1
1401,1419,1,,0,0,5,0,0,1
1419,1437,1,,0,0,5,0,0,1
*ngen, line=c, nset=out_out
1301,1319,1,,0,0,0,0,0,1
1319,1337,1,,0,0,0,0,0,1
1701,1719,1,,0,0,5,0,0,1
1719,1737,1,,0,0,5,0,0,1
*nfill, nset=out_nodes
out_in,out_out,3,100
*element, type=c3d8
1001,1001,1101,1102,1002,1401,1501,1502,1402
1036,1036,1136,1101,1001,1436,1536,1501,1401
*elgen, elset=out_elems
1001,35,1,1,3,100,100
*elgen, elset=out_elems
1036,3,100,100
*solid section, elset=out_elems, material=elas
*elset,elset=out_in,generate
1001,1036
*surface, name=out_in
out_in,s6
*surface interaction, name=smooth
*nset,nset=all_nodes
in_nodes,out_nodes
*equation
2
1,1,1.,50,1,1.
*boundary
all_nodes,3,3
1,2,2
50,2,2
1001,2,2
**
*CONTACT INITIALIZATION DATA, NAME=shrink_fit, INTERFERENCE FIT
**
*CONTACT
*CONTACT INCLUSIONS, ALL EXTERIOR
*CONTACT PROPERTY ASSIGNMENT
in_out, out_in, smooth
*SURFACE PROPERTY ASSIGNMENT, PROPERTY=GEOMETRIC CORRECTION
in_out,CIRCUM,0,0,0,0,0,1
out_in,CIRCUM,0,0,0,0,0,1
*CONTACT INITIALIZATION ASSIGNMENT
in_out,out_in,shrink_fit 
*step,nlgeom,unsymm=yes
*static,direct
0.1,1.0
*controls,parameter=field
1.e-15
*output,field, freq=1
*node output
u,
*element output
s,
*contact output
cstress,
*print,contact=yes
*end step