LAMMPS剪切模拟三把斧:变形、速度与力量
同样是做剪切模拟,面前摆着好几种方法,到底该选哪一个?它们之间到底有什么区别?命令又该怎么写?今天咱们就来把LAMMPS里最核心的三种剪切方法——全局形变剪切、速度边界剪切和外力加载剪切——彻底捋清楚.
核心思路:两种“折腾”体系的方式
在深入三种方法之前,先记住一个根本点:所有剪切模拟,本质上都是通过两种思路来“折腾”你的原子体系:
主动控“形变”:强制让体系产生一个固定的剪切速率或应变。这是最常用的思路,好比用模具硬压,看材料怎么“反抗”。
主动控“力”或“速度”:通过给边界原子施加力或速度,让它们“拖着”整个体系动起来。这更贴近真实实验里怎么“拉”或“推”材料。接下来,我们看看这三种方法分别是怎么实现这两条思路的。
方法一:全局形变剪切(fix deform)—— “均匀拉伸”的王者
这是LAMMPS官方最推荐、也是最标准的方法。它的核心操作是直接改变模拟盒子的形状。盒子一变形,里面的所有原子就会乖乖地跟着产生均匀的剪切变形。
就像这样:想象你有一块橡皮泥,你用手从两边均匀地推它,让它整体倾斜。fix deform干的就是这个活儿。
最佳用武之地:
研究晶体、非晶体、聚合物、流体的均匀剪切行为。
想要计算粘度、剪切模量、绘制应力-应变曲线这些宏观力学参数。
使用周期性边界条件时,这个方法没有边界效应,精度最高。
一句话优点:剪切均匀、参数可控、代码简洁、通用性强。
关键提醒:做流体剪切和粘度计算时,一定要用 nvt/sllod 系综,它能修正剪切流动带来的速度分布,避免温度计算失真。
方法二:速度边界剪切—— “摩擦拖拽”的高手
这种方法不改变盒子形状。它的做法是:固定住上下(或左右)边界的原子,然后给它们施加一个恒定的切向速度,让边界原子像“传送带”一样,通过层与层之间的摩擦,带动中间的材料发生剪切。
就像这样:你用手掌压住一本书的封面,另一只手压住封底,然后两只手反向搓动。书页之间就会产生相对滑动。速度边界剪切就是这个原理。
最佳用武之地:
研究摩擦、界面滑移、壁面附近的流变行为。
模拟薄膜、多层结构、受限体系(比如夹在两块板之间的液体)。
特别想观察边界本身对材料剪切行为的影响。
一句话特点:剪切从边界到中心是不均匀的(梯度分布),更贴近真实实验的加载方式,但均匀性不如方法一。
核心逻辑:分组固定边界原子 → 给边界原子一个恒定速度 → 让中间原子自由运动。
方法三:外力加载剪切—— “恒定拉力”的模拟器
前两种方法是控制“变形”或“速度”,而这一种方法直接对边界原子施加一个恒定的剪切力或剪切应力。相当于你用恒定的力去推材料的一侧,观察它怎么慢慢变形。
就像这样:你用一根绳子以恒定的力量去拉一个重物,重物在绳子的牵引下开始运动。外力剪切就是用恒定的“拉力”去驱动原子运动。
最佳用武之地:
研究材料的剪切屈服、塑性变形、蠕变(随时间缓慢变形)行为。
测试材料的临界剪切应力(多大劲儿能破坏它)和失稳阈值。
模拟“恒力加载,观察变形响应”的真实实验场景。
一句话特点:应力恒定,应变随材料的“脾气”自发变化。这与前两种“应变可控、应力自发变化”的思路完全互补。
核心逻辑:固定一侧边界 → 对另一侧边界原子施加恒定切向力 → 让体系自由变形。
关键补充:addforce 命令是在原子上叠加一个恒定的外力。这个力的大小需要根据你的体系尺寸和目标剪切应力来调整,适合做材料的“极限测试”。
总结:怎么选?看这张表就够了
评论
会魔法的猫

干货推荐

小胶质细胞
Mr弘🔬
63 8

有髓神经纤维超微结构
Mr弘🔬
57 5

【代码分享】费米面嵌套(Fermi surface nesting)函数计算
张阳小号1
651 2

lammps中文教程5
春风得意马蹄疾
48 2

计算机材料设计Materials-Studio教程11
活着
35 1

计算机材料设计Materials-Studio教程5
活着
33 2
