含高阶余项的非线性动力系统数值计算法  

A high-precision numerical method for non-linear dynamic systems

作　　者：李永强[1] 王薇[1] 刘杰[2] 刘宇[2] 金志强[1] 

机构地区：[1]东北大学理学院力学系,辽宁沈阳110004 [2]东北大学机械与自动化学院,辽宁沈阳110004

出　　处：《计算力学学报》2007年第5期555-559,共5页Chinese Journal of Computational Mechanics

基　　金：国家自然科学基金(598335050)资助项目

摘　　要：建立了一种求解非线性动力系统高精度数值计算的新方法,重构了等价的非线性动力系统方程,该方程考虑了非线性函数的任意高阶项,并给出了该方程的Duhamel积分表达式,在时间步长内用Newton-Raphson法进行数值迭代求解,该方法能连续满足微分方程而不只是在离散的步长端点满足方程,从而打破了传统的Euler型有限差分法。计算实例表明,该方法计算精度高于传统的Runge-Kutta,Newmark-β和Wilson-θ等方法。　A new high-precision numerical arithmetic for solving the non-linear dynamic system is proposed.The ordinary nonlinear dynamic equation is reconstructed,and the new equivalent equation taining arbitrary high order remainder.The arithmetic presents Duhamel integration expression,using Newton-Raphson iterative arithmetic to seek the numerical solution,satisfying the differential equation continuously rather than at discrete spots,therefore,the arithmetic exceeds the traditional Euler finite differential method.Compared with the traditional method,such as Runge-Kutta method,Newmark-β method and Wilson-θ et al.,the calculation precision of this method is much higher.

关 键 词：非线性动力系统 高阶余项 Duhamel积分 Newton-Raphson法 

分 类 号：O322[理学—一般力学与力学基础]