求解微分方程初值问题一种弧长法.pdfVIP

( ) 　应用数学和力学 ,第 20 卷 第 8 期 1999 年 8 月 应用数学和力学编委会编 　　　　　　　　　　　 　 　　Applied Mathematics and Mechanics 重 庆 出 版 社 出 版 ( ) 文章编号 :1000-0887 1999 08-0875-06 求解微分方程初值问题的一种弧长法 1 2 1 武际可 , 　许为厚 , 　丁红丽 ( 1 北京大学 力学与工程科学系 ,北京 100871 ; 2 香港科技大学 数学系 ,香港) (黄永念推荐) 摘要 : 　对于连续介质力学问题中导出的微分方程初值问题 ,常常具有解奇异性 ,如不连续、Stiff 性质或激波间断 ·　本文通过在相应空间, 引入一个或数个弧长参数变量 ,克服解的奇异性 ·　对于 常微分方程组引入弧长参数变量后 ,奇异性得以消除和削弱 ,应用一般的解常微分方程组的方法 ( ) 如 Runge- Kutta 法 求解 ·　对于偏微分方程引入弧长参数变量后,在相应的空间离散成常微分方 程组 ,用解奇异性常微分方程组相同的方法即可求解 ·　本文给出了两个算例·　 关 　键 　词 : 　微分方程 ; 　数值方法 ; 　弧长法 ; 　Stiff 方程 ; 　Burgers’方程 中图分类号: 　O322 ;O2418 　　　文献标识码 : 　A 引　　言 ( ) ( ) 对于连续介质力学问题导出的偏微分方程 组 或常微分方程 组 ,常常具有解的不连续、 Stiff 性质或激波间断等奇异性 ·　在国际上,连续介质力学问题中解的奇异性问题 ,一直是一 个尚未很好解决的问题 , 目前常用的方法如奇异摄动法和子波法,为克服奇异性而使计算量成 数量阶地增长 ·　 有限单元法是结构计算最有效的方法之一 ,在实际应用中 ,人们发现载荷增量法因其刚度 矩阵的奇异性而无法穿过载荷的极值点 ,Riks[1 ] 首次应用弧长法求解了带有参数的非线性方 程中的极值问题 ,20 年来已应用弧长法解决了大量的非线性问题[2 ,3 ] ·　伪弧长法的基本思想 是在解曲线上引入弧长参数 ,并增加一个约束方程 ·　本文把伪弧长法应用于求解带有奇异性 的偏微分方程或常微分方程组的初值问题·　对于常微分方程组引入弧长参数变量后,奇异性 ( ) 得以消除和削弱 ,应用一般的解常微分方程组的方法 如四阶 Runge- Kutta 法 求解 ·　对于偏 微分方程引入弧长参数变量后 ,在相应的空间离散成常微分方程组 ,用解奇异性常微分方程组 相同的方法即可求解 ·　这一方法因为仅在相应的空间引入弧长参量,来克服解的不连续、Stiff 性质和激波间断 ,与现有的方法相比 ,思路简单 ,程序实现容易 ,并且精度高 ·　本文给出了两个 算例 ,此方法也可推广应用于高维空间的微分方程的初值问题中去 ·　 1 　弧 长 方 法 考虑常微分方程初值问题 收稿日期 : 　1998-04-20 ; 修订日期 : 　1999-02-14 ( ) 基金项目 : 　国家自然科学基金资助项目 ( )