材料模拟软件：LIGGGHTS二次开发_LIGGGHTS计算模型与方法.docx

PAGE1

PAGE1

LIGGGHTS计算模型与方法

1.计算模型概述

LIGGGHTS（LaminarINcompressibleGeneral-GeometryGranularHeatTransferandTransportSimulator）是一款基于离散元方法（DEM）的材料模拟软件，主要用于模拟颗粒系统的动力学行为。DEM是一种数值方法，通过跟踪每个颗粒的运动和相互作用来模拟颗粒系统的宏观行为。LIGGGHTS的核心计算模型基于牛顿第二定律，即每个颗粒的运动方程。这些方程通过数值积分方法求解，从而得到颗粒的位移、速度和加速度。

1.1牛顿第二定律

牛顿第二定律是LIGGGHTS计算模型的基础，描述了每个颗粒的运动状态。对于一个颗粒(i)，其运动方程可以表示为：

[m_i={i}^{ext}+{ji}_{ij}]

其中：-(m_i)是颗粒(i)的质量。-(i)是颗粒(i)的位置向量。-({i}^{ext})是作用在颗粒(i)上的外部力，如重力、流体阻力等。-(_{ij})是颗粒(i)和颗粒(j)之间的接触力。

1.2接触力模型

接触力模型是DEM的核心部分，用于描述颗粒之间的相互作用。LIGGGHTS支持多种接触力模型，常见的包括Hertz-Mindlin模型和线性弹簧模型。这些模型通过定义颗粒接触的力和力矩来模拟颗粒之间的碰撞和摩擦。

1.2.1Hertz-Mindlin模型

Hertz-Mindlin模型是一种非线性接触模型，适用于模拟硬颗粒之间的接触。该模型考虑了颗粒的弹性形变和摩擦力。接触力({ij})可以分解为法向力({ij}^n)和切向力(_{ij}^t)：

[_{ij}^n=k_n_n]

[_{ij}^t=k_t_t]

其中：-(k_n)和(k_t)分别是法向和切向的弹簧刚度。-(_n)和(_t)分别是法向和切向的形变。-()和()分别是法向和切向的单位向量。

1.3数值积分方法

数值积分方法用于求解颗粒的运动方程。LIGGGHTS常用的数值积分方法包括Verlet算法和Predictor-Corrector算法。这些方法通过时间步长将连续的运动方程离散化，从而逐步计算颗粒的位置和速度。

1.3.1Verlet算法

Verlet算法是一种简单而有效的数值积分方法，适用于模拟颗粒系统的运动。该算法的核心思想是通过前一个时间步的位移和速度来预测当前时间步的位移，然后通过当前时间步的力来校正速度。具体步骤如下：

预测位置：[_i^{n+1}=_i^n+_i^nt+_i^nt^2]

计算加速度：[_i^{n+1}=]

校正速度：[_i^{n+1}=_i^n+(_i^n+_i^{n+1})t]

其中：-(_i^n)和(_i^n)分别是第(n)个时间步的位移和速度。-(_i^n)是第(n)个时间步的加速度。-(t)是时间步长。

1.4多体动力学模型

多体动力学模型用于描述多个颗粒之间的相互作用。LIGGGHTS通过定义颗粒之间的接触力和力矩来模拟多体动力学系统。此外，还可以考虑颗粒与壁面之间的相互作用。

1.4.1颗粒-颗粒接触

颗粒-颗粒接触的力和力矩可以通过Hertz-Mindlin模型来计算。例如，对于两个颗粒(i)和(j)，其接触力可以表示为：

[{ij}={ij}^n+_{ij}^t]

其中，法向力({ij}^n)和切向力({ij}^t)通过形变和弹簧刚度来计算。

1.4.2颗粒-壁面接触

颗粒与壁面之间的接触力可以通过类似的模型来计算。例如，对于一个颗粒(i)与壁面的接触，其法向力(_{i}^{wall,n})可以表示为：

[{i}^{wall,n}=k_n{i}^{wall}]

其中：-(k_n)是法向弹簧刚度。-(_{i}^{wall})是颗粒与壁面的法向形变。-()是法向单位向量。

1.5热传导模型

LIGGGHTS还支持热传导模型，用于模拟颗粒系统中的热传递过程。热传导模型通过定义每个颗粒的温度和热导率来计算颗粒之间的热交换。

1.5.1粒子温度更新

每个颗粒的温度可以通过以下方程更新：

[=({ji}{ij}-_i^{ext})]

其中：-(T_i)是颗粒(i)的温度。-(C_i)是颗粒(i)的热容。-(_{ij})是