baoluyao 发表于 2025-9-12 21:51
你用的公式和你的命令不对应。你的公式是用体系的总体压强进行计算，你的in文件中的命令使用压强分布进行计 ...

谢谢老师，我是计算了每个区块的应力张量然后积分计算的，具体的计算过程和公式确实有些差别，主要的问题是在lammps生成的原始数据处理上，因为compute stress/atom 计算的是单个原子的应力，需要变化到整个区块上的应力分量，就是这一步的计算，我没弄清楚。

你用的公式和你的命令不对应。你的公式是用体系的总体压强进行计算，你的in文件中的命令使用压强分布进行计算。你得到的结果需要用压强分布的积分进行计算，你取找一下公式就懂了。