第四讲 运动方程的积分算法.ppt

第四讲 运动方程的积分算法;积分运动方程的注意点;积分运动方程的注意点;积分运动方程的注意点;标准Verlet算法;标准Verlet算法;解决方法一是，在第二个时间步中，把该步力看成常量，且使用普通运动学方程:;流程;标准Verlet算法;Leap-frog算法;Leap-frog算法;速度Verlet算法;速度Verlet算法;Gear预测校正法;Gear预测校正法;预测－校正法允许的时间步长比其它算法长两倍以上 每个积分步内要计算两次体系势能，以得到原子间相互作用力。 该算法的稳定差，能量波动较大，较verlet算法占用更多的内存。 ;时间步长;???间步长;RMS (Reservoir Modeling System保守系统的模拟) Energy deviation;时间步长;宏观物理量的计算;系统 的 物 理性质是系统中粒子坐标和速度的函数，对于任意一个时刻宏观概念物理量A，定义为;势能 部 分 可以按势函数计算，对于不同的原子势函数表达式有不同的计算公式。 动 能计 算 公式为;应变 是 表 示物体变形大小的测度，可以有不同的定义方式。 应变是一个相对量，反映了物体相对初始构型的变形程度。 应变又是纯几何量，与空间尺度变化无关。 因此宏观、纳观尺度下应变可以采用相同的定义方式;原子应力;原子应力;原子应力;原子应力;宏观应力实际上是原子应力在一定程度上的统计平均。 如果要计算体系中某个区域（由region定义，可以是整个模拟盒）所围成的“块”的应力，在Lammps只需将该区域里的所有原子的单原子应力值加起来，再除以这个区域的体积即可，无须进行单个原子体积的计算。 ;谢谢