1Gpa下1m/s的滑動摩擦
加載過程中，我的F控溫和F不控溫狀態下，溫度都能達到一個平衡值，但是一到滑動就會出現EF，CD層溫度下降的趨勢，而且也不能最終平衡到某一溫度附近，不知道是我的代碼原因，還是速度太小？
代碼需要實現的功能：1Gpa下的1m/s的滑動摩擦
各位幫我看看唄，哪兒有問題，私信或郵箱我，謝謝！！！！！！！！
#-----------------------EF速度恆定
velocity E_F set NULL ${input_velocity} NULL sum yes
#-----------------------B層恆溫300K
fix 1 B nvt temp 300.0 300.0 50.0
reset_timestep 0
timestep 0.5
run 1000
#-----------------------C_D層nve
fix 2 C_D nve
fix 3 C_D wall/reflect zhi ${reflect2} units box #儘量保持密度一定
#-----------------------E_F層nve
fix 4 E_F nve/limit 0.05
fix 5 E_F addforce 0.0 0.0 -${f_to_addforce}
fix 6 F temp/rescale 10 300.0 300.0 10.0 1.0
reset_timestep 0
timestep 0.3
run 140000
