求助计算多层二硫化钼弯曲时的层间滑移能

彭正

我打算通过多层二硫化钼弯曲时的层间能量变化来判断是那一层率先滑移。在实际向下中，中性层上半部分出现脱离。不知道应该怎么解决这个问题，是速度太快了，或者有什么其他的方法。
我的势函数是sw/mod和ilp/tmd,弯曲方式是将中性层设定为纯弯曲。希望大佬们可以给点建议。

variable        dx atom "((sin(v_theta)/v_kappa_safe) - v_x0) * v_ramp"
variable        dz atom "(-(1.0 - cos(v_theta))/v_kappa_safe) * v_ramp"
fix             4 L1 move variable v_dx NULL v_dz NULL NULL NULL