Commands list

×Contents and search
EXPAND ALL SORT A-Z

PARTICLE_SPH

Particle

*PARTICLE_SPH
"Optional title"
sid, $N_p^{sid}$
mid, gid, correction, $s$, $dsf$, E2S, slide, $t_{end}$
$\Delta$

Parameter definition

Variable
Description
sid
Subdomain ID
$N_p^{sid}$
Number of particles
default: not used
mid
Material ID
gid
Geometry ID
correction
Dynamic SPH node coordinate correction
options:
0 $\rightarrow$ weak correction
1 $\rightarrow$ strong correction
2 $\rightarrow$ no correction
$s$
Smoothing length scale factor
default: 2.5
$dsf$
Particle density scale factor (particles per unit volume)
default: 1
E2S
SPH particles are activated as elements are eroded (element to SPH mapping)
options:
0 $\rightarrow$ no
1 $\rightarrow$ yes
slide
Sliding contact flag
options:
0 $\rightarrow$ full merge
1 $\rightarrow$ sliding only
$t_{end}$
Particle deactivation time or FUNCTION (fcn)
options: time or fcn
default: not used
$\Delta$
Target particle spacing (only used if map=1)
default: automatic

Description

SPH subdomain definition.

The subdomain ID (sid) determines in which order particles are filled into the global domain. In case subdomains are overlapping, the domain with the largest ID will overwrite (remove) particles belonging to domains with lower domain ID's.

Material models presented with an SPH badge in the command manual are supported by SPH.

The smoothing length (radius) is defined as:

$\displaystyle{ R = s \cdot h }$

where $h$ is the initial SPH node spacing and the default value of $s=2.5$.

There is, by default, a dynamic repositioning of SPH coordinates that helps supressing tensile instability modes. This feature can be turned off by setting correction=2.

The SPH particle density scale factor $dsf$ is optional and it is used to control the relative particle density (number of particles per unit volume) in the different subdomains. Increasing $dsf$ for a subdomain means that a larger fraction of the total number of particles will be allocated to that specific subdomain.

$N_p^{sid}$ is optional. If not defined, then this sub-id will use a fraction of the the particles $N_p$ defined in PARTICLE_DOMAIN.

By default, SPH particles using different material ID's will merge when coming into contact. However, this interaction can be changed to friction free sliding by setting slide=1.

If an SPH particle becomes isolated from its neighbours, it is automatically deactivated.

With E2S=1 SPH particles are generated dynamically to replace eroded elements. This option requires that gid refers to a GEOMETRY_PART.