lammps应力计算

Question

在lammps的real单位里，如何进行总体体系的应力计算，real单位下应力的单位是什么

Answer 1

上午好☀️☀️☀️️
本答案参考ChatGPT-3.5

在 LAMMPS 中，用于计算整个体系的应力（stress）以及压力（pressure）的命令是“compute stress/atom”，计算完毕后，可以用“thermo_style”命令来输出实时结果。

步骤如下：

1. 在 LAMMPS 中，设置模拟系统的单位为 real，命令为“units real”。

2. 在输入文件中添加 compute stress/atom 命令以计算应力。根据定义，应力单位为能量或力除以面积，real 单位下应力的单位是：

（能量单位 / 长度的平方）=（kJ / mol / nm²）

或

（力单位 / 长度）=（bar，可以转换为 MPa）。

下面是 compute stress/atom 的示例代码：

compute myStress all stress/atom NULL
thermo_style custom step press pxx pyy pzz pe etotal

在此示例中，“myStress”是计算的 ID 号，"all" 表示对所有原子进行计算，“NULL”表示不对组群进行分析。

在运行时，LAMMPS 会计算完整个体系的应力，并输出压力和晶格应变（strain）等附加数据。

1. 添加 thermo_style 命令以输出力、应力和压力的结果。

在此示例中，“thermo_style” 命令的“custom”选项指定了输出多个变量。输出格式为：

Step Pxx Pyy Pzz PE E_total

这意味着，在每个时间步长上，LAMMPS将输出 step、Pxx、Pyy、Pzz、势能（PE）和总能量（E_total）的值。

希望这些信息对你有帮助！