Commands list

×Contents and search
EXPAND ALL SORT A-Z

GEOMETRY_PART

Geometries

*GEOMETRY_PART
"Optional title"
gid, csysid
pid

Parameter definition

Variable
Description
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