码迷,mamicode.com
首页 > 其他好文 > 详细

Verlet Integration

时间:2016-03-18 23:29:00      阅读:321      评论:0      收藏:0      [点我收藏+]

标签:

 
 

Verlet Integration


Verlet 积分法是一种用于求解牛顿运动方程的数值方法,被广泛运用于动力学模拟以及视频游戏中。尔莱算法的优点在于:数值稳定性比简单的欧拉方法高很多,并保持了物理系统中的时间可逆性与相空间体积元体积守恒的性质。

基本韦尔莱算法

根据牛顿运动方程有

技术分享
代入到粒子的位移关于时间步的泰勒展开式中有:
技术分享
得到
技术分享
同理
技术分享
两式相加得
技术分享

技术分享
新位置的计算误差为四阶,技术分享 为时间步。因而韦尔莱算法中不涉及速度,如果希望得到速度,可以从前面的两式相减得出
技术分享

 

速度表示的韦尔莱算法

一般地,速度表示的韦尔莱算法更为常用,它可以给出同一时间变量下的速度和位置。它实际上与基本韦尔莱算法等价,精度相同。

首先对位置进行泰勒展开

技术分享
技术分享

 

对两式相减可得

技术分享
将最初的Verlet公式中的 技术分享 换成 技术分享
技术分享
代入前式,可得
技术分享
此式即为速度表示的韦尔莱算法。实际常用的计算步骤为

 

  1. 首先通过泰勒展开式 技术分享 计算得到位置 技术分享
  2. 技术分享 和系统的相互作用势条件(如果相互作用仅依赖位置 技术分享 )可以求得力场 技术分享
  3. 由速度表示的韦尔莱公式求出新的速度 技术分享
 
 

Verlet Integration

标签:

原文地址:http://www.cnblogs.com/crackpotisback/p/5293991.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!