*Heading
 Self contact with c3d10h
**
** PARTS
**
*Part, name=Part-1
*End Part
*Part, name=Part-2
*End Part
*Part, name=Part-3
*End Part
**
** ASSEMBLY
**
*Assembly, name=Assembly
**
*Instance, name=Part-1-1, part=Part-1
*Node
*Element, type=C3D10H
** Region: (Section-1:Picked)
*Elset, elset=_PickedSet2, internal, generate
** Section: Section-1
*Solid Section, elset=_PickedSet2, material=Material-1
1.,
*End Instance
**
*Instance, name=Part-3-1, part=Part-3
*Node
*Nset, nset=Part-3-1-RefPt_, internal
*Surface, type=CYLINDER, name=r-base
START,          -15.,          -4.
 LINE,           4.,          -4.
*Rigid Body, ref node=Part-3-1-RefPt_, analytical surface=r-base
*End Instance
**
*Instance, name=Part-2-1, part=Part-2
*Node
*Nset, nset=Part-2-1-RefPt_, internal
*Surface, type=REVOLUTION, name=r-top
START,           1.,          1.6
 LINE,           1.,         -1.6
*Rigid Body, ref node=Part-2-1-RefPt_, analytical surface=r-top
*End Instance
*Nset, nset=_PickedSet12, internal, instance=Part-2-1
*Nset, nset=_PickedSet13, internal, instance=Part-3-1
*Nset, nset=_PickedSet14, internal, instance=Part-2-1
*Nset, nset=_PickedSet15, internal, instance=Part-3-1
*Nset, nset=boun-1, instance=Part-1-1
*Elset, elset=boun-1, instance=Part-1-1
*Nset, nset=boun-3, instance=Part-1-1
*Elset, elset=boun-3, instance=Part-1-1
*Elset, elset=_inside_S3, internal, instance=Part-1-1
*Elset, elset=_inside_S2, internal, instance=Part-1-1
*Elset, elset=_inside_S1, internal, instance=Part-1-1
*Elset, elset=_inside_S4, internal, instance=Part-1-1
*Surface, type=ELEMENT, name=inside
_inside_S3, S3
_inside_S2, S2
_inside_S4, S4
_inside_S1, S1
*Elset, elset=_outside_S3, internal, instance=Part-1-1
*Elset, elset=_outside_S4, internal, instance=Part-1-1
*Elset, elset=_outside_S1, internal, instance=Part-1-1
*Elset, elset=_outside_S2, internal, instance=Part-1-1
*Surface, type=ELEMENT, name=outside
_outside_S3, S3
_outside_S4, S4
_outside_S2, S2
_outside_S1, S1
*End Assembly
**
** MATERIALS
**
*Material, name=Material-1
*Hyperelastic
1000.,0.,0.01
**
** INTERACTION PROPERTIES
**
*SURFACE INTERACTION, NAME=BASE
*Friction
0.5,
*SURFACE INTERACTION, NAME=TOP
*Friction
0.5,
*Surface Interaction, name=inside
**
** INTERACTIONS
**
** Interaction: Int-1
*Contact Pair, interaction=top, type=surface
outside, Part-2-1.r-top
** Interaction: Int-2
*Contact Pair, interaction=base, type=surface
outside, Part-3-1.r-base
** Interaction: Int-3
*Contact Pair, interaction=inside, type=surface
inside,
** ----------------------------------------------------------------
**
**
** STEP: Step-1
**
*Step, name=Step-1, nlgeom=YES, unsymm=yes
*Static
0.01, 1.
**
** BOUNDARY CONDITIONS
**
*Boundary
_PickedSet12, 1, 1
_PickedSet12, 2, 2, -6.2
_PickedSet12, 3, 6
_PickedSet13, 1, 6
boun-1, 1, 1
boun-3, 3, 3
**
** OUTPUT REQUESTS
**
**
** FIELD OUTPUT: F-Output-1
**
*Output, field, variable=PRESELECT, freq=10
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT, freq=10
*End Step
** ----------------------------------------------------------------
**
** STEP: Step-2
**
*Step, name=Step-2, nlgeom=YES, unsymm=yes, inc=1000
*Static
0.05, 1.
**
** BOUNDARY CONDITIONS
**
*Boundary,op=new
**
** comment out the next two lines
** and uncomment out the succeeding two lines
** to run the full analysis
**
_PickedSet14, 1, 1, -1.25
_PickedSet14, 6, 6, -1.0
**_PickedSet14, 1, 1, -10.
**_PickedSet14, 6, 6, -8.
_PickedSet14, 2, 2, -6.2
_PickedSet14, 3, 5
_PickedSet15, 1, 6
boun-3, 3, 3
*End Step