*Heading
** Incremental DSA for C3D10M with
** Small Sliding Contact
** Analytical Rigid Surface
** Shape and Size parameter used
*Parameter
shape=0.0
dp=1.0
c1=80.0*dp
*Design Parameter
shape,dp
*Parameter shape variation, parameter=shape
top_out,0.0,0.10,0.0
middle,0.0,0.05,0.0
*node,nset=c3d41
1, 0.0,  1.0, 0.0
2, 0.0, -1.0, 0.0
3, 2.0, -1.0, 0.0
4, 2.0,  1.0, 0.0
5, 0.0,  0.0, 0.0
6, 1.0, -1.0, 0.0
7, 2.0,  0.0, 0.0
8, 1.0,  1.0, 0.0
9, 1.0,  0.0, 0.0
*Ncopy,old set=c3d41,change number=10,shift,multiple=1,new set=c3d42
0.0, 0.0, 1.0
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*Ncopy,old set=c3d42,change number=10,shift,multiple=1,new set=c3d43
0.0, 0.0, 1.0
0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0
*Nset,nset=c3d4
C3D41,c3d42,c3d43
*Element,elset=solid,type=c3d10m
1, 21, 2, 23, 22, 15, 16, 29, 25, 12, 26
2, 21,23, 2, 3, 29, 16, 15, 19, 13, 6
3, 21, 3, 2, 1, 19, 6, 15, 11, 9, 5
4, 21, 23, 4, 24, 29, 17, 18, 28, 27, 14
5, 21, 4,23, 3, 18, 17, 29, 19, 7, 13
6, 21, 4, 3, 1, 18, 7, 19, 11, 8, 9
*Solid Section, elset=solid, material=hyper
1.,
*Nset,Nset=top_out
1,11,21,8,18,28,4,14,24
*Nset,Nset=middle
7,9,13,17,19,29
** ----------------------------------------------------------------
**
** PART INSTANCE: RIGID DIE
**
**
** PART INSTANCE: Part-2-1
**
*System
1.,     1.95235,          1.,          2.,     1.95235,          1.
*Node
36, -2.55814e-07,     -0.85235,  4.63245e-16
37,   -0.0476214,         -0.5,     -1.99943
38,           2.,         -0.5, -2.44929e-16
39,     -1.99891,         -0.5,    0.0660631
40,    0.0382983,         -0.5,      1.99963
41,       1.3475,    -0.695107, -8.51132e-17
42,     0.704092,    -0.809841,  7.24764e-17
43,   -0.0242672,    -0.782267,    -0.902664
44,      1.03453,         -0.5,     -1.71165
45,      1.75216,         -0.5,    -0.964329
46,     -1.79038,         -0.5,    -0.891363
47,     -1.11292,         -0.5,     -1.66175
48,    -0.733927,     -0.80609,     0.025908
49,     -1.34202,    -0.696218,    0.0457642
50,      1.76587,         -0.5,     0.938988
51,       1.0836,         -0.5,      1.68102
52,    0.0162298,    -0.780446,     0.914428
53,     -1.01898,         -0.5,      1.72096
54,     -1.73013,         -0.5,      1.00332
55,      1.29386,     -0.67679,    -0.591602
56,     0.760655,    -0.746399,     -0.80641
57,     -1.32045,    -0.675736,    -0.540757
58,    -0.806848,    -0.744498,    -0.774431
59,      1.30728,    -0.674875,     0.580349
60,     0.786546,       -0.744,     0.798633
61,     -1.28515,    -0.674945,     0.627204
62,    -0.763727,    -0.743058,     0.827043
63,           0.,     -0.85235,  2.44929e-16
*Surface, type=REVOLUTION, name=master
1., 1.95235, 1., 1., 101.952, 1.
START,           2.,         -0.5
CIRCL,           0.,     -0.85235,           0.,           5.
*RigidBody, anal=master, ref=ref
*Nset, nset=ref
63,
*System
*Nset, Nset=base
2,12,22,6,16,26,3,13,23
**
** MATERIALS
**
*Material, name=hyper
*Hyperelastic
<c1>,   20., 0.001
**
** INTERACTION PROPERTIES
**
*Surface Interaction, name=smooth
*Surface,type=element,name=slave
4,S4
6,S2
*Elset, elset=master, generate
13,  32,
** Interaction: Int-1
*Contact Pair, interaction=smooth, smallsliding
slave, master
*Nset, Nset=topbotm
1,4,24,21,14,8,11,28,23,13,26,12,
2,6
*Elset,Elset=Resel
3,4,5
** ----------------------------------------------------------------
**
** STEP: Step-1
**
*Step, nlgeom, dsa, inc=100
Step-1
*Static
0.1,1.0
** BOUNDARY CONDITIONS
**
** BC: BC-1 Type: Displacement
*Boundary
base, 1, 3, 0.0
ref, 1, 1, 0.0
ref, 3, 6, 0.0
ref, 2, 2, -0.3
**
** OUTPUT REQUESTS
**
*Output, Field, freq=6
*Node output, Nset=topbotm
U,RF
*Element output, Elset=Resel
S
*Contact output
Cstress,
*Design Response, freq=6
*Node Response,Nset=topbotm
U,RF
*Element Response, Elset=Resel
S
*Contact Response
Cstress
*End Step