大家好,我是小马老师。
lammps常见的拉伸方法有两种:velocity拉伸和deform拉伸。
两种拉伸的原理不同,deform拉伸是通过调整box尺寸来重新设置原子位置,velocity拉伸通过固定一端不动、移动另一端的方式进行拉伸,从拉伸原理上看,velocity拉伸和拉伸实验更接近。
在实际模拟中,两种拉伸结果也不同,velocity也更容易出现“缩颈”,容易被拉断。
但velocity拉伸也有不足之处,拉伸变形不是均速的,特别是当拉伸速度比较小时,可能会有“拉不动”的现象。
曾经遇到过几例类似的案例,如SiC裂纹拉伸拉不动,设置拉伸速度后,裂纹仅仅扩大一点点就会闭合。
在聚乙烯/石墨烯拔出的案例中,拔出速度较小时,石墨烯无法从聚乙烯中拔出。
(1)原因分析
产生这种现象的主要原因是velocity命令仅设置原子的初始速度,原子速度被设定之后,velocity命令就不再起作用。
在模拟过程中,这部分原子的速度在系综的作用下会发生改变,原子速度较大时,这种改变可以忽略不计,原子初始速度较小时,这种速度的变化就显得比重较大。
如高速行驶的汽车轻踩刹车,速度会有所降低,但是总体速度还是挺快。如果是低速行驶的汽车,按同样力度踩刹车,可能车就停住了。
(2)解决方法
大部分情况下,velocity可以满足要求,当出现“拉不动”的现象时,可换一种速度设置方法,如fix move或fix addforce等命令。
设置fix move的原子不能再设其他系综,原子速度不会受到系综影响,能够保持匀速拉伸效果,基本能解决以上问题。
公众号:lammps加油站