非线性动力学常微分方程组高精度数值积分方法

被引:22
作者
郑兆昌
沈松
苏志霄
机构
[1] 清华大学工程力学系,清华大学工程力学系,清华大学工程力学系北京,北京,北京
关键词
非线性动力学; 有限差分法; 重构等价非线性微分方程; 高阶余项; Duhamel卷积分;
D O I
暂无
中图分类号
O241.4 [数值积分法、数值微分法];
学科分类号
070102 [计算数学];
摘要
建立了一种求解非线性动力学常微分方程组初值问题的新方法.若非线性函数一阶导数存在,则给出解的积分方程表达式,计算得到按规定误差要求的高精度数值解.引入一般自治或非自治非线性系统的首次近似Jacobi矩阵,不作任何假设重构等价的非线性常微分方程组,简捷而有广泛的适应性,不改变方程的本质,但其主项构成线性化方程组,其它项则代表非线性函数高阶余项而不涉及Taylor级数展开计算,给出该方程组初值问题的Duhamel卷积分解析表达式,在时间步长内进行数值积分迭代求解,在指定误差内快速收敛,逐步递推获得非线性常微分方程的瞬态响应和全时域高精度数值解.积分解连续满足微分方程组而不是在离散的步长端点上满足代数方程组,打破了传统用增量法在离散点上建立的代数方程组迭代求解,从而使传统:Euler型逐步积分法的各种差分格式算法改变成真正的积分格式算法.数值计算中给出指数矩阵递增展开式,变矩阵乘法为乘积系数的加法,避免了大量矩阵自乘而大大提高计算效率.算法验证为无条件稳定,则保证对线性常微分方程而言,计算中舍入误差的传播不会扩散,不出现计算机字长有限而引起舍入误差导致计算不确定性问题.基于以上理论和数值方法,计算了线性非线性算例并进行了分析,验证了本方法简捷而有广泛的适应性,可以
引用
收藏
页码:284 / 295
页数:12
相关论文
共 6 条
[1]
非线性动力系统线性模型数值计算的Taylor变换法 [J].
苏志霄 ;
郑兆昌 ;
高永毅 .
力学学报, 2002, (04) :586-593
[2]
NUMERICAL METHOD IN DYNAMIC RESPONSE OF NONLINEAR SYSTEMS [J].
郑兆昌 ;
谭明一 .
Applied Mathematics and Mechanics(English Edition), 1985, (01) :93-101
[3]
非线性系统动力分析的模态综合技术 [J].
郑兆昌 .
应用数学和力学, 1983, (04) :563-572
[4]
Nineteen Dubious Ways to Compute the Exponential of a Matrix.[J].Cleve Moler;Charles Van Loan.SIAM Review.1978, 4
[5]
RANDOMLY TRANSITIONAL PHENOMENA IN THE SYSTEM GOVERNED BY DUFFINGS EQUATION [J].
UEDA, Y .
JOURNAL OF STATISTICAL PHYSICS, 1979, 20 (02) :181-196
[6]
非线性振动分析.[M].褚亦清;李翠英编著;.理工大学出版社.1996,