8. 基础知识
此部分描述了如何使用 LAMMPS 为用户和开发人员执行各种任务。术语表页面还列出了 MD 术语,以及相应 LAMMPS 手册页的链接。 LAMMPS 源代码分发的 examples
目录中包含的示例输入脚本以及示例脚本页面上突出显示的示例输入脚本还展示了如何设置和运行各种模拟。
8.1.通用基础知识
8.2. 设置入门
8.3. 分析入门
8.4. 力场入门
8.5. 软件包入门
8.5.1.有限尺寸球形和非球形粒子
8.5.2. 粒度模型
颗粒系统由具有直径的球形颗粒组成,而不是点颗粒。这意味着它们具有角速度,并且可以向它们施加扭矩以使它们旋转。
要运行粒度模型的模拟,您将需要使用以下命令:
-
atom_style sphere
-
fix nve/sphere 修复 nve/球体
-
fix gravity 固定重力
这个计算
-
compute erotate/sphere 计算旋转/球体
计算可与热力学信息一起输出的旋转动能。计算
-
compute fabric 计算结构
计算颗粒和非颗粒pair styles的结构张量的各种版本。
使用这 4 对势能之一,计算相互作用的粒子对之间的力和扭矩:
-
pair_style gran/history
-
pair_style gran/no_history
-
pair_style gran/hertzian
-
pair_style granular
这些命令实现特定于粒度系统的修复选项:
-
fix freeze
-
fix pour
-
fix viscous
-
fix wall/gran
-
fix wall/gran/region
fix style freeze将冻结原子的力和扭矩归零,并且应用于粒状系统而不是固定样式设置力。
为了模拟热传导,必须添加温度和热流原子变量:
- fix property/atom
温度积分fix
-
fix heat/flow
以及两者中定义的热传导选项
-
pair_style granular
-
fix wall/gran
为了提高计算效率,您可以使用以下命令消除冻结原子之间不必要的成对计算:
-
neigh_modify exclude
注:默认情况下,对于 2d 系统,颗粒粒子仍被建模为 3d 球体,而不是 2d 圆盘(圆),这意味着它们的转动惯量将与 3d 中的相同。如果您希望将 2d 颗粒粒子建模为 2d 圆盘,请参阅Howto 2d 文档页面上有关此主题的注释,其中讨论了 2d 模拟。
要添加自定义颗粒接触模型,请参阅modifying granular sub-models page.