*PARAMETER
term       = 0.03                   # termination time
cx         = 0                      # center point x
cy         = 0                      # center point y
r_thick    = 0.0046                 # rubber thickness
r_radius   = 0.0268                 # rubber radius
r_dist     = 0.066384               # rubber distance from origo
s_thick    = 0.0005                 # structure thickness
s_radiusA1 = %r_radius-0.0001       # structure radiusA based on rubber radius
s_radiusB1 = 0.0230                 # structure radiusB
s_radiusA2 = %s_radiusA1-%s_thick   # structure radiusA based on rubber radius
s_radiusB2 = %s_radiusB1-%s_thick   # structure radiusB
s_height   = 0.0166                 # structure height
s_start_A1 = %r_dist+%r_thick       # help parameter
s_start_B1 = -%s_start_A1           # help parameter
s_start_A2 = %s_start_A1+%s_height  # help parameter
s_start_B2 = %s_start_B1-%s_height  # help parameter
b_radius  = 0.01                    # sphere radius
# Simulation controls ---------------------------------------------------------
*UNIT_SYSTEM
SI
*TIME
[%term]
*OUTPUT
[%term/100], [%term/1000], 1
# Structure definition --------------------------------------------------------
*COMPONENT_CYLINDER
"Top rubber"
1, 1, 1, 6
[%cx], [%cy], [%r_dist], [%cx], [%cy], [%s_start_A1], [%r_radius]
*COMPONENT_CYLINDER
"Bottom rubber"
2, 1, 1, 6
[%cx], [%cy], [-%r_dist], [%cx], [%cy], [%s_start_B1], [%r_radius]
*COMPONENT_SPHERE
"Sphere"
3, 2, 3
[%cx], [%cy], [-%r_dist+%b_radius], [%b_radius]
*COMPONENT_PIPE
"Top structure"
4, 3, 5, 50, 1
[%cx], [%cy], [%s_start_A1], [%cx], [%cy], [%s_start_A2], [%s_radiusA1], [%s_radiusA2]
[%s_radiusB1], [%s_radiusB2]
*COMPONENT_PIPE
"Bottom structure"
5, 3, 5, 50, 1
[%cx], [%cy], [%s_start_B1], [%cx], [%cy], [%s_start_B2], [%s_radiusA1], [%s_radiusA2]
[%s_radiusB1], [%s_radiusB2]
*SMOOTH_MESH
P, 2, 45
*CHANGE_P-ORDER
PS, 12, 3
# Part definition -------------------------------------------------------------
*PART
"Rubber"
1, 1, 0, 0, 45
*PART
"Sphere"
2, 2, 0, 0.001
*PART
"Structure"
3, 2, 0, 0.001
*SET_PART
12
1, 2
# Contact definition ----------------------------------------------------------
*CONTACT
"Sphere to Rubber"
2
P, 2, P, 1, 0.0, -100e9
*CONTACT
"Structure to Rubber"
1
P, 1, P, 3, 0.0, -100e9
# Load definition -------------------------------------------------------------
*LOAD_DAMPING
P, 1, 12
*FUNCTION
12
100.0
# Material definition ---------------------------------------------------------
*MAT_USER_3
"IMPETUS example user material"
1, 1000.0, 1.0e8
1.0e6, 1.0e6
*MAT_RIGID
2, 1.0e4, 100.0e9, 0.0
# Boundary condition ----------------------------------------------------------
*BC_MOTION
"Initial velocity for sphere"
P, 2, 0, 0, 0, 0, 0, 0.007
V, Z, 1,-1.0
*BC_MOTION
"Constraint sphere translation (XY) and rotation (XYZ)"
P, 2, XY, XYZ
*BC_MOTION
"Constraint structure translation (XYZ) and rotation (XYZ)"
P, 3, XYZ, XYZ
*BC_MOTION
"Constraint top rubber translation with nodes from geometry 1"
G, 1, XYZ
*BC_MOTION
"Constraint top rubber translation with nodes from geometry 2"
G, 2, XYZ
*CURVE
"Sphere velocity"
1
0.0,  0.0
0.01, 20.0
2.0,  20.0
*GEOMETRY_SEED_COORDINATE
"Clamp A"
1
7.70444e-9, -0.0268, 0.0694507
*GEOMETRY_SEED_COORDINATE
"Clamp B"
2
-7.70444e-9, -0.0268, -0.0679173
*END
