*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 smoothing, name=smooth
in_out,out_in,circum,0,0,0,0,0,1
*contact pair, interaction=smooth, type=surface, geometric=smooth, small sliding
in_out,out_in
*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
*step,nlgeom,unsymm=yes
*static
0.1,1.0
*contact interference, shrink
in_out,out_in
*controls,parameter=field
1.e-15
*output,field
*node output
u,
*element output
s,
*contact output, master=out_in, slave=in_out
CSTRESS,
*end step