*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
*Ncopy,old set=c3d41,change number=10,shift,multiple=1,new set=c3d42
*Ncopy,old set=c3d42,change number=10,shift,multiple=1,new set=c3d43
*Nset,nset=c3d4
C3D41,c3d42,c3d43
*Element,elset=solid,type=c3d10m
*Solid Section, elset=solid, material=hyper
1.,
*Nset,Nset=top_out
*Nset,Nset=middle
** ----------------------------------------------------------------
**
** PART INSTANCE: RIGID DIE
**
**
** PART INSTANCE: Part-2-1
**
*System
*Node
*Surface, type=REVOLUTION, name=master
START,           2.,         -0.5
CIRCL,           0.,     -0.85235,           0.,           5.
*RigidBody, anal=master, ref=ref
*Nset, nset=ref
*System
*Nset, Nset=base
**
** 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
** Interaction: Int-1
*Contact Pair, interaction=smooth, smallsliding, TYPE=SURFACE TO SURFACE
slave, master
*Nset, Nset=topbotm
*Elset,Elset=Resel
** ----------------------------------------------------------------
**
** 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