干货分子动力学 干货详情

LAMMPS剪切模拟三把斧:变形、速度与力量

会魔法的猫68
LAMMPS模拟分子动力学剪切模拟非平衡分子动力学NVT/SLLOD系综全局形变剪切速度边界剪切外力加载剪切流变行为粘度计算应力应变曲线周期性边界条件摩擦与界面滑移受限体系塑性变形蠕变临界剪切应力

同样是做剪切模拟,面前摆着好几种方法,到底该选哪一个?它们之间到底有什么区别?命令又该怎么写?今天咱们就来把LAMMPS里最核心的三种剪切方法——全局形变剪切速度边界剪切外力加载剪切——彻底捋清楚.

核心思路:两种“折腾”体系的方式

在深入三种方法之前,先记住一个根本点:所有剪切模拟,本质上都是通过两种思路来“折腾”你的原子体系:

主动控“形变”:强制让体系产生一个固定的剪切速率或应变。这是最常用的思路,好比用模具硬压,看材料怎么“反抗”。

主动控“力”或“速度”:通过给边界原子施加力或速度,让它们“拖着”整个体系动起来。这更贴近真实实验里怎么“拉”或“推”材料。接下来,我们看看这三种方法分别是怎么实现这两条思路的。

方法一:全局形变剪切(fix deform)—— “均匀拉伸”的王者

这是LAMMPS官方最推荐、也是最标准的方法。它的核心操作是直接改变模拟盒子的形状。盒子一变形,里面的所有原子就会乖乖地跟着产生均匀的剪切变形。

就像这样:想象你有一块橡皮泥,你用手从两边均匀地推它,让它整体倾斜。fix deform干的就是这个活儿。

最佳用武之地

研究晶体、非晶体、聚合物、流体的均匀剪切行为。

想要计算粘度、剪切模量、绘制应力-应变曲线这些宏观力学参数。

使用周期性边界条件时,这个方法没有边界效应,精度最高。

一句话优点剪切均匀、参数可控、代码简洁、通用性强

关键提醒:做流体剪切和粘度计算时,一定要用 nvt/sllod 系综,它能修正剪切流动带来的速度分布,避免温度计算失真。

方法二:速度边界剪切—— “摩擦拖拽”的高手

这种方法不改变盒子形状。它的做法是:固定住上下(或左右)边界的原子,然后给它们施加一个恒定的切向速度,让边界原子像“传送带”一样,通过层与层之间的摩擦,带动中间的材料发生剪切。

就像这样:你用手掌压住一本书的封面,另一只手压住封底,然后两只手反向搓动。书页之间就会产生相对滑动。速度边界剪切就是这个原理。

最佳用武之地

研究摩擦、界面滑移、壁面附近的流变行为

模拟薄膜、多层结构、受限体系(比如夹在两块板之间的液体)。

特别想观察边界本身对材料剪切行为的影响。

一句话特点:剪切从边界到中心是不均匀的(梯度分布),更贴近真实实验的加载方式,但均匀性不如方法一。

核心逻辑:分组固定边界原子 → 给边界原子一个恒定速度 → 让中间原子自由运动。

方法三:外力加载剪切—— “恒定拉力”的模拟器

前两种方法是控制“变形”或“速度”,而这一种方法直接对边界原子施加一个恒定的剪切力或剪切应力。相当于你用恒定的力去推材料的一侧,观察它怎么慢慢变形。

就像这样:你用一根绳子以恒定的力量去拉一个重物,重物在绳子的牵引下开始运动。外力剪切就是用恒定的“拉力”去驱动原子运动。

最佳用武之地

研究材料的剪切屈服、塑性变形、蠕变(随时间缓慢变形)行为。

测试材料的临界剪切应力(多大劲儿能破坏它)和失稳阈值。

模拟“恒力加载,观察变形响应”的真实实验场景。

一句话特点应力恒定,应变随材料的“脾气”自发变化。这与前两种“应变可控、应力自发变化”的思路完全互补。

核心逻辑:固定一侧边界 → 对另一侧边界原子施加恒定切向力 → 让体系自由变形。

关键补充addforce 命令是在原子上叠加一个恒定的外力。这个力的大小需要根据你的体系尺寸和目标剪切应力来调整,适合做材料的“极限测试”。

总结:怎么选?看这张表就够了


评论

0/1000发布评论
全部评论

会魔法的猫

关注
TA的主页

干货推荐