lammps案例：CNT拉伸案例代码解析

大家好，我是小马老师。
本文和大家共同欣赏一个CNT拉伸的lammps案例，案例代码由Simon Gravelle编写，原文：

https://github.com/simongravelle/LAMMPS/tree/main/CNTdeformationAIREBO

关于拉伸的案例，公众号也推送了多个案例，但这个案例所用方法与平时所用方法有点区别。

lammps拉伸模拟最常见的方式是deform和velocity拉伸，使用compute stress/atoms或者pxx等方式统计拉伸应力。
案例作者并没有使用以上两种方式，而是循环调用displace_atoms命令逐渐向两侧移动固定端实现拉伸。
在计算拉伸力时，调用fix setforce命令的结果，方法比较独特。
对于固定两端的拉伸，在计算应变时应该只计算中间mobile部分的伸长量，本文给出了比较具体的计算方法，值得大家参考。
下面是案例in文件全文ÿ