LAMMPS静水压力与冯·米塞斯应力计算脚本

LAMMPS 后阶段的应力脚本，挺适合拿来直接上手用的。它从文件里读应力张量数据，每行六个分量那种，像σxx, σyy, σzz, σxy, σyz, σzx，格式比较标准。脚本会算两个指标：静水压力和冯·米塞斯应力，也就是常说的体积平均应力和材料失效准则那一套。

静水压力的计算挺直接，就是主应力的平均值。冯·米塞斯应力那块就稍微绕点，用的是公式，但脚本都封装好了，你只管跑。输出结果格式也清爽，拿来画图或者写论文用都方便。

这个脚本最适合你在模拟完之后局部单元的应力场，比如研究材料变形、裂纹分布啥的。要注意的是，数据文件要准备好，单位统一，不然算出来会偏差挺大。

如果你经常在 LAMMPS 里做力学，这段代码真的可以收着。搭配一些论文，比如 非静水压力条件下巷道围岩破碎区应力分布，还能互相印证一下，研究会更深入。