反思02

时间步长与运行步数

之前的我的运行时间过短，基本都是0.3 x 150000的组合。对EF加载（我目前改为15MPa，作为练习的话建议先1Gpa运行看看与论文对比）。运行时间如果过短，可能会出现EF层“压”不住CD的现象，这是正常的，辣鸡UP在此花了很长的时间，现在还在跑程序。。。。。。一有结果就会总结在博客，目前的运行状态是1.0 x 500000，祝我好运吧！

定压滑动

目前的定压滑动采用的是move linear+addforce，由于目前还没有什么好结果，但是这俩命令在我一个测试小系统中都是可以运行的！！速度和力都是起效果的，这里有个细节，就是处理施加速度时温度的处理，相关可以看我命令更新部分关于compute temp/partial的介绍。

fix			13 E_F move linear NULL ${input_velocity} NULL
fix			14 E_F addforce 0.0 0.0 -${f_to_addforce}

之前的效果不好主要还是其他的一些细节问题。。。UP正在疯狂调试中!!!！

摩擦力，摩擦系数 ，chunk温度 etc.

这部分正确性没法儿保证，但是大家可以先熟悉相关命令，在摩擦这方面应该可以起到很大的作用！！！

#-----------------------摩擦系数计算
compute			1 D reduce ave fx fy fz
compute			2 E reduce ave fx fy fz
compute			3 C_D reduce ave fx fy fz
compute			4 E_F reduce ave fx fy fz

fix 			1 D ave/time 100 3 1000 c_1[*] file cof1.txt mode vector
fix 			2 E ave/time 100 3 1000 c_2[*] file cof2.txt mode vector
fix 			3 C_D ave/time 100 3 1000 c_3[*] file cof3.txt mode vector
fix 			4 E_F ave/time 100 3 1000 c_4[*] file cof4.txt mode vector

#-----------------------滑动块温度计算
compute 		5 E_F chunk/atom type
compute    		mytemp1 E_F temp/chunk 5 temp com yes
fix			5 E_F ave/time 100 3 1000 c_mytemp1[*] file temp1.txt mode vector

compute 		6 E chunk/atom type
compute    		mytemp2 E temp/chunk 6 temp com yes
fix			6 E ave/time 100 3 1000 c_mytemp2[*] file temp2.txt mode vector

#-----------------------滑动位移计算
compute 		7 E_F msd com yes average yes
fix			7 E_F ave/time 100 3 1000 c_7[*] file dis1.txt mode vector
compute 		8 E msd com yes average yes
fix			8 E ave/time 100 3 1000 c_8[*] file dis2.txt mode vector

#-----------------------表层速度计算
compute			9 C_D reduce ave vx vy vz
fix			9 C_D ave/time 100 3 1000 c_9[*] file v1.txt mode vector

nve/limit

该命令不能用在数据提取阶段！！！！！！！！！！！！！
啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊啊！！！！！！！！！！！！
详情见基础命令！！