材料模拟软件：LIGGGHTS二次开发_LIGGGHTS中的颗粒动力学模拟.docx

PAGE1

PAGE1

LIGGGHTS中的颗粒动力学模拟

在上一节中，我们讨论了LIGGGHTS的基本功能和使用方法。本节将深入探讨LIGGGHTS中的颗粒动力学模拟，包括颗粒的运动方程、相互作用力的计算、边界条件的设置以及如何进行模拟设置和分析。

1.颗粒运动方程

在LIGGGHTS中，颗粒的运动遵循牛顿第二定律。每个颗粒的位置和速度通过以下方程进行更新：

[=m]

[=]

[(t+t)=(t)+t]

[(t+t)=(t)+(t+t)t]

其中：-()是作用在颗粒上的总力-(m)是颗粒的质量-()是颗粒的加速度-((t))和((t+t))分别是颗粒在时间(t)和(t+t)时的速度-((t))和((t+t))分别是颗粒在时间(t)和(t+t)时的位置-(t)是时间步长

LIGGGHTS使用Verlet算法来积分这些方程，这是一种常用的数值积分方法，具有较高的精度和稳定性。Verlet算法的基本形式如下：

[(t+t)=(t)+(t)t+(t)t^2]

[(t+t)=(t)+((t)+(t+t))t]

通过这些方程，LIGGGHTS能够准确地模拟颗粒在不同力作用下的运动。

2.颗粒相互作用力

颗粒间的相互作用力是颗粒动力学模拟的核心。LIGGGHTS支持多种颗粒间力的模型，包括但不限于：

线性弹簧-阻尼模型：颗粒间通过线性弹簧和阻尼器相互作用，适用于模拟弹性碰撞和粘性行为。

Hertz-Mindlin模型：基于赫兹接触理论，适用于模拟刚性颗粒的接触行为。

Coulomb摩擦模型：颗粒间通过库仑摩擦力相互作用，适用于模拟摩擦效应。

粘性力模型：颗粒间通过粘性力相互作用，适用于模拟流体环境中的颗粒行为。

线性弹簧-阻尼模型

线性弹簧-阻尼模型是最常用的颗粒间力模型之一。其基本方程为：

[=k(_1-2-d{})+(_1-_2)]

其中：-(k)是弹簧常数-()是阻尼系数-(_1)和(_2)是两个颗粒的位置-(_1)和(2)是两个颗粒的速度-(d{})是颗粒间的平衡距离

在LIGGGHTS中，可以通过以下命令设置线性弹簧-阻尼模型：

pair_stylegran/hooke

pair_coeff**1.0e51.0e31.0e-31.0e-3

pair_stylegran/hooke：选择线性弹簧-阻尼模型

pair_coeff**1.0e51.0e31.0e-31.0e-3：设置模型参数，依次为弹簧常数、阻尼系数、静摩擦系数和动摩擦系数

Hertz-Mindlin模型

Hertz-Mindlin模型基于赫兹接触理论，适用于模拟刚性颗粒的接触行为。其基本方程为：

[_n=E^*^{3/2}]

[_t=(_n,_t)]

其中：-(_n)是法向力-(_t)是切向力-(E^*)是有效弹性模量-(R)是有效半径-()是接触深度-()是摩擦系数-()是切向阻尼系数-(_t)是切向速度

在LIGGGHTS中，可以通过以下命令设置Hertz-Mindlin模型：

pair_stylegran/hertz

pair_coeff**1.0e71.0e-31.0e-31.0e-3

pair_stylegran/hertz：选择Hertz-Mindlin模型

pair_coeff**1.0e71.0e-31.0e-31.0e-3：设置模型参数，依次为弹性模量、法向阻尼系数、静摩擦系数和动摩擦系数

3.边界条件

边界条件的设置对于颗粒动力学模拟至关重要，它们决定了颗粒在模拟区域内的行为。LIGGGHTS支持多种边界条件，包括：

周期性边界条件：模拟区域在各个方向上是周期性的，适用于模拟无限大的颗粒系统。

固定边界条件：模拟区域的边界是固定的，适用于模拟有限大小的颗粒系统。

移动边界条件：边界可以在模拟过程中移动，适用于模拟流体-颗粒相互作用等场景。

周期性边界条件

周期性边界条件可以通过以下命令设置：

boundaryppp

boundaryppp：在x、y、z三个方向上设置周期性边界条件

固定边界条件

固定边界条件可以通过以下命令设置：

boundarysss

boundarysss：在x、y、z三个方向上设置固定边界条件

移动边界条件

移动边界条件可以通过以下命令设置：

fix1allwall/granzplane0.01