*Heading
** Job name: shell Model name: Model-1
**
** PARTS
**
*Part, name=Part-1
*End Part
*Part, name=Part-2
*End Part
**
** ASSEMBLY
**
*Parameter
thick=0.1
mu=0.2
Or=45.0
*Design parameter
thick,mu,Or
*Dsa Controls
Or, CD, 0.01
*Assembly, name=Assembly
**  
*Instance, name=Part-1-1, part=Part-1
0., 0., -1.
*Node
      1,           2.,           0.,           0.
      2,          1.5,           0.,           0.
      3,           1.,           0.,           0.
      4,          0.5,           0.,           0.
      5,           0.,           0.,           0.
      6,           2.,           0.,          0.5
      7,          1.5,           0.,          0.5
      8,           1.,           0.,          0.5
      9,          0.5,           0.,          0.5
     10,           0.,           0.,          0.5
     11,           2.,           0.,           1.
     12,          1.5,           0.,           1.
     13,           1.,           0.,           1.
     14,          0.5,           0.,           1.
     15,           0.,           0.,           1.
     16,           2.,           0.,          1.5
     17,          1.5,           0.,          1.5
     18,           1.,           0.,          1.5
     19,          0.5,           0.,          1.5
     20,           0.,           0.,          1.5
     21,           2.,           0.,           2.
     22,          1.5,           0.,           2.
     23,           1.,           0.,           2.
     24,          0.5,           0.,           2.
     25,           0.,           0.,           2.
*Element, type=S4R, elset=shells
 1,  1,  2,  7,  6
 2,  2,  3,  8,  7
 3,  3,  4,  9,  8
 4,  4,  5, 10,  9
 5,  6,  7, 12, 11
 6,  7,  8, 13, 12
 7,  8,  9, 14, 13
 8,  9, 10, 15, 14
 9, 11, 12, 17, 16
10, 12, 13, 18, 17
11, 13, 14, 19, 18
12, 14, 15, 20, 19
13, 16, 17, 22, 21
14, 17, 18, 23, 22
15, 18, 19, 24, 23
16, 19, 20, 25, 24
** Section: Section-1
*Shell Section, elset=shells, material=Material-1
<thick>, 5
*REBAR LAYER,ORIENT=ORI
RB,1.0,0.2,0,ELREB,0,1 
*Orientation,name=ori
1.0, 0.0, 0.0, 1.0, 0.0, -2.0
3, <Or>
*End Instance
**  
*Instance, name=Part-2-1, part=Part-2
*Node
      1,          1.5,          0.5
*Nset, nset=ref, internal
 1,
*Surface, type=CYLINDER, name=master, internal
START,          1.5,          0.5
PARAB,           1.,           0.,          0.5,          0.5
*Rigid Body, ref node=ref, analytical surface=master
*End Instance
*Nset, nset=edges, internal, instance=Part-1-1
  1,  5,  6, 10, 11, 15, 16, 20, 21, 25
*Nset, nset=middle, internal, instance=Part-1-1
  3,  8, 13, 18, 23
*Nset, nset=ref, internal, instance=Part-2-1
 ref,
*Elset, elset=__slave_SPOS, internal, instance=Part-1-1, generate
  1,  16,   1
*Surface, type=ELEMENT, name=_slave, internal
__slave_SPOS, SPOS
*Nset, nset=top, internal, instance=Part-1-1
22, 23, 24
*Nset, nset=bot, internal, instance=Part-1-1
2, 3, 4
*End Assembly
** 
** MATERIALS
** 
*Material, name=Material-1
*Elastic
 2e+09, 0.3
*Material, name=elreb
*Elastic
10000.0, 0.2
** 
** INTERACTION PROPERTIES
** 
*Surface Interaction, name=IntProp-1
1.,
*Friction
<mu>
** 
** INTERACTIONS
** 
** Interaction: Int-1
*Contact Pair, interaction=IntProp-1
_slave, Part-2-1.master
** ----------------------------------------------------------------
** 
** STEP: Step-1
** 
*Step, nlgeom, dsa, unsymm=yes
Step-1
*Static
0.1, 1., 1e-05, 1.
** 
** BOUNDARY CONDITIONS
** 
** Name: BC-1 Type: Symmetry/Antisymmetry/Encastre
*Boundary
edges, ENCASTRE
middle, 1,1
** Name: BC-2 Type: Displacement/Rotation
*Boundary
ref, 1, 1
ref, 3, 6
** 
** LOADS
** 
** Name: Load-1   Type: Concentrated force
*Cload
ref, 2, -5.E6
top, 3, 10000.0
bot, 3,-10000.0
** 
** OUTPUT REQUESTS
** 
*Output, field
*Element output
S
*Node output
U
*Contact output
Cstress
*Design response,freq=6
*Element response
S,
*Node response
U
*Contact response
Cstress
*controls, parameters=field, field=displacement
1e-05,1e-04
*controls, parameters=field, field=rotation
1e-05,1e-04
*End Step