基于非传统哈密顿变分原理的高阶辛算法.PDFVIP

第 32 卷 第 3 期 应 用 力 学 学 报 Vol.32 No.3 2015 年 6 月 CHINESE JOURNAL OF APPLIED MECHANICS Jun. 2015 文章编号：1000- 4939(2015) 03-0410-07 基于非传统哈密顿变分原理的高阶辛算法 富明慧 陆克浪 李纬华 黄策 张荧荧 ( 中山大学 工学院 510275 广州) 摘要：给出了非传统哈密顿变分原理的一种简化形式，并在此基础上利用拉格朗日多项式近似位 移和动量，采用高斯积分法对时间积分，建立了针对动力学初值问题的一类高阶辛算法。在建立 高阶辛算法的过程中，本文方法与基于传统哈密顿变分原理的辛算法不同，无需由端值问题向初 值问题转换，因此更加简捷有效。此外，给出了线性动力问题中本文算法保辛性的证明。当位移、 动量的插值次数和高斯积分点个数均为 m 时，本文算法是具有 2m 阶精度的辛算法，且是线性无 条件稳定的。通过数值算例结果表明，本文算法与辛算法性质吻合，并且计算效率比同阶辛龙格 库塔法提高了约 50%。 关键词：哈密顿系统；动力学初值问题；非线性；辛算法；非传统哈密顿变分原理 中图分类号：O322 文献标识码：A DOI ：10.11776/cjam.32.03.D021 的计算效率。基于离散的变分原理，文献[15-19]构 1 引 言 造了一系列能保辛和保动量或保辛且保能量和动量 的算法；文献[20]通过修正作用量给出了动力学初 辛算法[1-2]是求解哈密顿系统的一种有效的数 值问题的一系列保辛算法。 可以看出，变分原理是构造保辛方法的一个有 值积分方法，具有保辛结构的特点。大量实验已经 效途径，然而采用的变分原理不同，构造的保辛方 证明辛算法在哈密顿系统计算中有明显的优越性， [3] 法也是有所区别的。上述方法采用的都是哈密顿变 特别是在长时间仿真计算中。辛算法在经典力学 、 [4] [5] [6-7] 分原理或修正的哈密顿变分原理，属于端值型变分 多体动力学 、分子动力学 、量子力学 、天体 力学[8-9]、电磁场问题[10-11]、热传导问题[12-13]、大气 原理，须通过一定的技巧才能求解初值问题，因此 科学[14]等领域得到了广泛应用。 应用不够简便。为解决这一矛盾，本文首先给出对 应于动力学初值问题的非传统哈密顿变分原理。然 早期辛算法的构造主要基于哈密顿对偶方程， 后，用拉格朗日多项式近似位移和动量，采用高斯 这方面的研究已经取得相应的进展。然而基于微分 积分法对时间积分，并基于该变分原理推导出相应 方程构造高阶辛算法难度较大，因而人们逐步把关 保辛算法的计算公式。最后，通过数值算例讨论该 注的焦点转向了哈密顿变分原理。通过变分原理构 法的线性稳定性、保辛性、精度和数值性能。 造保辛算法，简便易行且具有很好的通用性和较高 基金项目：国家自然科学基金 ；