基于修正的Magnus方法的高振荡动力系统的数值积分方法.pdfVIP

应用数学和力学，第27卷第10期 MathematicsandMechanics Applied 2006年10月15日出版 v01．27。No．10．Oct．15．2006 文章编号：lOOO-0887(2006)lO一1211—08 1000—0887 ⑥应用数学和力学编委会，ISSN 基于修正的Magnus方法的高振荡 动力系统的数值积分方法。 李文成1，邓子辰2¨，黄永安2 (1．西北工业大学理学院应用数学系，西安710072； 2．西北工业大学力学与土建学院工程力学系，西安710072： 3．大连理工大学工业装备结构分析国家重点实验室，大连116023) (岳珠峰推荐) 摘要：基于建立于一般线性动力系统上的Magnus数值积分方法，针对随时闻而高频率振荡的二 阶动力系统，给出了有效的修正Magnus数值积分算法．首先，将二阶动力系统重新表示为一阶系 统的形式，通过引进新变量进行参考坐标变换，使动力系统的高振荡性质保留在新形式内；进而基 于局部线性化技术用修正的Magnus方法求解新形式下的系统方程；最后，通过一系列数值实验说 明了文中方法的有效性． 关键词：动力系统；高振荡；Magnus方法；Hamilton系统 中图分类号：0322；0241文献标识码：A 引 言 考虑具有初值问题的二阶微分方程 y(0)=Yo，Y(O)=Yo， (】) Y+g(t)Y=0， 并且其解的振荡性态的时间尺度远小于积分区间．本文讨论这类高振荡动力系统的数值计算 问题．高振荡动力系统广泛存在于许多领域，不仅在诸如分子动力系统，电路仿真系统，柔性 体系统等领域，甚至振荡现象渗透于太空中卫星天线的折叠情形，以及车辆悬轴导向系统的定 位过程中．我们更感兴趣的问题是如何对形如式(1)的高振荡问题的模型方程进行有效的数 值求解，该问题是解决高振荡问题的首要条件． 高振荡微分方程的数值积分方法，历来是微分方程的数值求解问题中，最富有挑战性的 难点问题之一．经典数值算法效果甚微，这也激励人们探讨新颖而富有创造性的数值方 为固定频率)的微分方程的三角函数积分方法；他的方法很快推广到形如：多+Ay=g(t，Y)， *收稿13期：2005—08．16；修订日期：2006-04-06 资金项目： 国家自然科学基金资助项目；教育部新世纪优秀人才计划资助项目(NCET-04- 0958)；大连理工大学工业装备结构分析国家重点实验室开放基金资助项目 作者简介： mail：wenchengli@nwpu．edu．on)． 12ll 万方数据 李文成 邓子辰 黄永安 Impulse 将研究推进到与时间相关的A(t)情形，发展了完全不同的方法来数值求解高振荡线性微分方 程．基于Iserles的思想，本文开展了以下工作：首先在第1节给出一般线性微分方程的基于 Magnus展式的数值逼近格式，然后给出高振荡微分方程的修正Magnus积分方法；第2节进一 步基于Cayl