Mesh commands

"Optional title"
coid, entype, enid, $m_{add}$, distribution

Parameter definition

Command ID
Entity type
options: N, G, P, PS
Entity ID
Added mass or added mass per unit area (constant, CURVE or FUNCTION)
options: constant, fcn
Mass distribution type
0 $\rightarrow$ node mass weighted distribution ($m_{add}$ defines total mass)
1 $\rightarrow$ area weighted distribution ($m_{add}$ defines total mass)
2 $\rightarrow$ area weighted distribution ($m_{add}$ defines mass per unit area)
3 $\rightarrow$ directional area weighted distribution ($m_{add}$ defines mass per unit area)


This command is used to add mass to a node or geometrical region.

Directional distribution

The option distribution=3 can be used to model added mass effects for bodies submerged in a liquid. The added node mass is then a vector $\mathbf{v}_m = \left[ v_{m,x} \; v_{m,y} \; v_{m,z} \right]^T$:

$\displaystyle{ \mathbf{v}_m = \left\{ \begin{array}{c} v_{m,x} \\ v_{m,y} \\ v_{m,z} \end{array} \right\} = m_{add} \cdot A \cdot \hat{\mathbf{n}} = m_{add} \cdot A \cdot \left\{ \begin{array}{c} n_x \\ n_y \\ n_z \end{array} \right\}}$

where $A$ is the surface area represented by the node and $\hat{\mathbf{n}}$ is the surface normal vector at the node. That is, the added mass effect is assumed to be proportional to the area (of the body) projected in the acceleration direction. The node acceleration $\mathbf{a}$ at a given node force $\mathbf{F} = \left[ F_x \; F_y \; F_z \right]^T$ becomes:

$\displaystyle{ \mathbf{a} = \left\{ \begin{array}{c} a_x \\ a_y \\ a_z \end{array} \right\} = \left\{ \begin{array}{c} F_x / (m + v_{m,x}) \\ F_y / (m + v_{m,y}) \\ F_z / (m + v_{m,z}) \end{array} \right\} = \left\{ \begin{array}{c} F_x / (m + m_{add} A n_x) \\ F_y / (m + m_{add} A n_y) \\ F_z / (m + m_{add} A n_z) \end{array} \right\} }$

where $m$ is the physical node mass. Note that it remains for the user to define $m_{add}$. As a rule of thumb it should be proportional to the smallest projected side length of the body (in the acceleration direction) multiplied by the liquid density.

One can take a rigid square plate with side length $L$ submerged in a liquid with density $\rho$ as an example. According to potential flow theory, $m_{add}$ should be (Patton K.T., Tables of Hydrodynamic Mass Factors for Translational Motion, ASME, 1965):

$\displaystyle{ m_{add} = 0.3758 \rho L }$


Node mass

The following command adds a point mass of $0.13 \mathrm{kg}$ to a node with ID=1001.

*UNIT_SYSTEM SI *ADD_MASS "point mass" 1, N, 1001, 0.13
Time dependent mass

A mass of $2.5 \mathrm{kg}$ is distributed to a surface and then removed after $0.01 \mathrm{s}$.

*UNIT_SYSTEM SI *ADD_MASS "distributed mass" 1, G, 55, fcn(345), 1 *GEOMETRY_SEED_COORDINATE 55 0.023, 0.45, 0.01 *FUNCTION 345 H(0.01-t) * 2.5