A complete model of a glued connection exposed to tensile loading.
*UNIT_SYSTEM
SI
*PARAMETER
tol = 1.0e-4, "tolerance"
delta = 2.0e-3, "discretization"
w = 4.0e-3, "width"
h = 1.0e-4, "film thickness"
rho = 1000.0, "density"
E = 1.0e6, "Young's modulus"
pr = 0.45, "Poisson's ratio"
sig_f = 1.0e7, "failure stress in tension"
tau_f = 5.0e6, "failure stress in shear"
G_I = 1.0e4, "delamination energy - modus I"
G_II = 1.0e4, "delamination energy - modus II"
Lp = 0.1, "flange size"
Rp = 0.02, "pipe radius"
hp = 0.002, "sheet thickness"
disp = 0.01, "pipe displacement"
tend = 0.01, "termination time"
*TIME
[%tend]
#
# MESH
*COMPONENT_BOX
"flange L"
1, 1, 1, 10, 10
[-%hp], [-%Lp/2], [-%Lp/2], 0, [%Lp/2], [%Lp/2]
*COMPONENT_BOX
"flange R"
2, 2, 1, 10, 10
0, [-%Lp/2], [-%Lp/2], [%hp], [%Lp/2], [%Lp/2]
*COMPONENT_PIPE
"pipe L"
3, 3, 10, 12, 1
[-%Lp-%hp], 0, 0, [-%hp], 0, 0, [%Rp], [%Rp+%hp]
*COMPONENT_PIPE
"pipe R"
4, 4, 10, 16, 1
[%hp], 0, 0, [%Lp+%hp], 0, 0, [%Rp], [%Rp+%hp]
*CHANGE_P-ORDER
ALL, 0, 3
*SMOOTH_MESH
ALL, 0, 45.0
#
# MATERIAL
*MAT_METAL
1, 2700.0, 70.0e9, 0.3
1
*FUNCTION
1
150.0e6 + 100.0e6*(1 - exp(-5*epsp))
#
# PARTS
*PART
"flange L"
1, 1
"flange R"
2, 1
"pipe L"
3, 1
"pipe R"
4, 1
#
# CONNECTIONS
*MERGE
"pipe L to flange L"
P, 3, P, 1
*MERGE
"pipe R to flange R"
P, 4, P, 2
*CONNECTOR_GLUE_LINE
"flange"
1
P, 1, P, 2, 22, [%tol], [%delta], [%w]
[%h], [%rho], [%E], [%pr], [%sig_f], [%tau_f], [%G_I], [%G_II]
*PATH
"glue line"
22
0, [-0.4*%Lp], [-0.4*%Lp]
0, [0.4*%Lp], [-0.4*%Lp]
0, [0.4*%Lp], [0.4*%Lp]
0, [-0.4*%Lp], [0.4*%Lp]
0, [-0.4*%Lp], [-0.4*%Lp]
#
# BOUNDARY CONDITIONS
*BC_MOTION
"left"
1
G, 123
V, X, 34, -1
*BC_MOTION
"right"
2
G, 234
V, X, 34
*FUNCTION
34
smooth_v(%disp, 0, %tend)
#
# CONTACT
*CONTACT
"general"
1
ALL, 0, ALL, 0
#
# PART SET
*SET_PART
12
1, 2
#
# GEOMETRIES
*GEOMETRY_BOX
123
[-%Lp-1.1*%hp], 0, 0, [-%Lp-0.9*%hp], 0, 0
*GEOMETRY_BOX
234
[%Lp+0.9*%hp], 0, 0, [%Lp+1.1*%hp], 0, 0
*END