GEOMETRY_PART

Geometries
*GEOMETRY_PART
"Optional title"
gid, csysid
pid
Parameter definition
VariableDescription
gid Geometry identification number
csysid Local coordinate system ID ($\gamma$SPH only)
pid Part ID
Description

This command is used to define a geometry with the shape of a part.

Only $\gamma$SPH:

A local coordinate system can be defined, to move $\gamma$SPH parts generated with the axisymmetric filling routine.

Example
Geometry defined by *GEOMETRY_PART and filled with water particles

The commands below defines a part, and a geometry with the shape of this part. The geometry is then filled with water particles.


*UNIT_SYSTEM
SI
*PARAMETER
%lx = 1.0
%ly = 0.1
%s = 0.2
%ang = 10
%lz = 1.2*tan(%ang)
# Dummy part nr 200 consisting of one element defined by 8 nodes
*COMPONENT_BOX_IRREGULAR
1, 1, 1, 1, 1
0.0, 0.0, 0.0, 0.0, [%ly/2], 0.0
[%lx], [%ly/2], 0.0, [%lx], 0.0, 0.0
0.0, 0.0, [%s], 0.0, [%ly/2], [%s]
[%lx], [%ly/2], [%s+%lx*tan(%ang)], [%lx], [0], [%s+%lx*tan(%ang)]
*MAT_RIGID
1, 7800
*MAT_FLUID
"Water"
2, 1000.0, 2.25e9
*PART
1, 1
*GEOMETRY_PART
"Geometry to be filled with water particles"
2
1
# Define domain for which the SPH-particles live
*PARTICLE_DOMAIN
, , 1e5
0.0, 0, 0.0, [%lx], [%ly/2], [%s+1.2*%lx*tan(%ang)]
1, 1, 1, 1, 1, 0
*PARTICLE_SPH
1
2, 2
*END
Water bottle

In the example below GEOMETRY_PART is used to fill the inside of a water bottle with SPH particles with material properties of water.


*UNIT_SYSTEM
SI
*PARAMETER
R1 = 0.035, "bottle radius"
R2 = 0.007, "bottle neck radius"
R3 = 0.006, "impactor radius"
t = 0.003, "bottle wall thickness"
h1 = 0.1, "bottle height 1"
h2 = 0.23, "bottle height 2"
h3 = 0.25, "bottle height 3"
N = 4, "mesh density parameter"
Np = 100000, "number of SPH particles"
*TIME
0.001, 0, 2.0e-8
#
# --- MESH ---
#
*COMPONENT_CYLINDER
"bottle 1"
1, 1, 1, [%N]
0, 0, 0, 0, 0, [%t], [%R1]
*COMPONENT_PIPE
"bottle 2"
2, 1, 1, [4*%N], 1
0, 0, 0, 0, 0, [%t], [%R1], [%R1+%t]
*COMPONENT_PIPE
"bottle 3"
3, 1, 70, [4*%N], 1
0, 0, [%t], 0, 0, [%h3], [%R1], [%R1+%t]
*MERGE_DUPLICATED_NODES
P, 1, P, 1, [%t/10]
#
*COMPONENT_SPHERE
"impactor"
4, 2, 4
0, 0, [%h3+10*%R3], [%R3]
*COMPONENT_BOX
"table"
5, 3, 10, 10, 1
[-2*%R1], [-2*%R1], [-%t], [2*%R1], [2*%R1], 0
*COMPONENT_PIPE
"dummy water"
6, 4, 60, [8*%N], 1
0, 0, [%t], 0, 0, [%h2], 0, [%R1]
#
*TRANSFORM_MESH_CYLINDRICAL
"bottle"
1, P, 1, 1, 111, 111
*TRANSFORM_MESH_CYLINDRICAL
"dummy water"
2, P, 4, 1, 0, 111
*COORDINATE_SYSTEM_CYLINDRICAL
1, 0, 0, 0
0, 0, 1, 1, 0, 0
*FUNCTION
111
0.5 * H(z-%h1) * (cos(180*min(1,(z-%h1)/(%h2-%h1))) - 1) * (%R1 - %R2)
#
*CHANGE_P-ORDER
ALL, 0, 3
*SMOOTH_MESH
ALL, 0, 45.0
#
# --- MATERIAL ---
#
*MAT_METAL
1, 7800.0, 210.0e9, 0.3
1
*FUNCTION
1
5.0e8 + 5.0e8*epsp^0.3
*MAT_RIGID
2, 7800.0
*MAT_FLUID
"water"
3, 1000.0, 2.0e9
#
# --- PART ---
#
*PART
"rubber bottle"
1, 1
"bullet"
2, 2
"table"
3, 2
"water dummy"
4
#
# --- BC ---
#
*BC_MOTION
"table"
3
P, 3, XYZ, XYZ
*INITIAL_VELOCITY
P, 2, 0, 0, -1000.0
#
# --- CONTACT ---
#
*CONTACT
"general contact"
1
PS, 123, PS, 123
*SET_PART
123
1, 2, 3
#
# --- PARTICLE ---
#
*PARTICLE_DOMAIN
PS, 123, [%Np]
-0.1, -0.1, -0.01, 0.1, 0.1, 0.3

0, 1
*PARTICLE_SPH
"water"
1
3, 1
*GEOMETRY_PART
"water"
1
4
*END
$\gamma$SPH example
*UNIT_SYSTEM
SI
*PARAMETER
R = 50.0e-3
dx = 5.0e-3
Nb = %R/%dx
angle = 30
*PART
"Reference mesh"
1,1
"SPH Axi -1"
2,2
*COMPONENT_PIPE
"Reference Mesh"
1,1,1,200,4
0,1.5*%R,-0.5*%dx,0,1.5*%R,0.5*%dx,%R,%R+%dx
*GEOMETRY_PART
"Geometry from mesh"
2,2
1
*COORDINATE_SYSTEM_FIXED
"local csys"
2,0.0,0.0,3.0*%R
cos(%angle),0,-sin(%angle),0,1,0
*MAT_RIGID
1, 7800.0
*MAT_JC
2, 2700, 70e9, 0.3
350e6
*SPH
0
2, 2, %dx, -1.0, -1.0
*END