材料模拟软件：LIGGGHTS二次开发_LIGGGHTS中的力场与势能函数.docx

PAGE1

PAGE1

LIGGGHTS中的力场与势能函数

在LIGGGHTS（Large-scaleAtomic/MolecularMassivelyParallelSimulatorforGranularandGranular-FluidSystems）中，力场和势能函数是模拟颗粒系统行为的核心部分。力场定义了颗粒之间的相互作用力，而势能函数则描述了这些力如何随着颗粒之间的相对位置变化。通过合理选择和定制力场与势能函数，可以准确模拟各种材料的行为，从简单的球形颗粒到复杂的多面体颗粒，从干燥的颗粒系统到含有流体的颗粒系统。

力场的基本概念

在LIGGGHTS中，力场（ForceField）是用于计算颗粒之间相互作用力的一组参数和公式。力场的选择直接影响模拟结果的准确性和可靠性。LIGGGHTS支持多种力场，包括经典的刚性球模型、软球模型、粘性模型等。每种力场都有其特定的应用场景和优势。

刚性球模型

刚性球模型是最简单的力场模型之一，假设颗粒是不可变形的刚性球。颗粒之间的接触力通过Hertz接触理论来计算。Hertz接触理论适用于弹性接触，可以描述颗粒之间的法向和切向力。

原理

Hertz接触理论基于弹性力学，假设颗粒在接触点处发生弹性变形。法向力(F_n)和接触变形()之间的关系可以通过以下公式表示：

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

其中：-(E^*)是颗粒材料的有效杨氏模量，定义为()，(E)为杨氏模量，()为泊松比。-(R)是颗粒的等效半径，对于两个半径分别为(R_1)和(R_2)的颗粒，(R)可以表示为()。

切向力(F_t)通过Coulomb摩擦定律来计算：

[F_t=(F_n,F_{})]

其中：-()是摩擦系数。-(F_{})是切向力的最大值，通常由切向变形(_t)和法向力(F_n)共同决定。

代码示例

以下是一个使用刚性球模型的LIGGGHTS输入文件示例，模拟两个球形颗粒的弹性接触。

#LIGGGHTSinputscriptforasimplerigidspheremodel

#Definethesimulationbox

regionboxblock010010010

create_box2box

#Definetheparticletypes

mass11.0

mass21.0

#Definetheinitialpositionsoftheparticles

create_atoms1random110012345box

create_atoms2random110012346box

#Definetheforcefield

pair_stylegran/hertz/history

pair_coeff111.00.30.10.001

pair_coeff121.00.30.10.001

pair_coeff221.00.30.10.001

#Definetheboundaryconditions

boundaryppp

#Definethetimeintegrationmethod

timestep0.001

run1000

代码解释

regionboxblock010010010：定义一个10x10x10的模拟盒子。

create_box2box：创建一个包含两种类型的颗粒的模拟盒子。

mass11.0和mass21.0：定义两种颗粒的质量。

create_atoms1random110012345box和create_atoms2random110012346box：在模拟盒子中随机创建各1个颗粒，类型分别为1和2，12345和12346是随机数种子。

pair_stylegran/hertz/history：选择刚性球模型的Hertz接触力场。

pair_coeff111.00.30.10.001、pair_coeff121.00.30.10.001和pair_coeff221.00.30.10.001：设置颗粒之间的接触参数，分别为杨氏模量、泊松比、摩擦系数和粘性系数。

boundaryppp：定义周期性边界条件。

timestep0.001：设置时间步长。

run1000：运行1000个时间步。

软球模型

软球模型（