| 非线性动力学常微分方程组高精度数值积分方法 |
| |
| 引用本文： | 郑兆昌,沈松,苏志霄.非线性动力学常微分方程组高精度数值积分方法[J].力学学报,2003,35(3):284-295. |
| |
| 作者姓名： | 郑兆昌 沈松 苏志霄 |
| |
| 作者单位： | 清华大学工程力学系,北京,100084 |
| |
| 基金项目： | 国家重点基础研究发展规划资助项目(998020316)，国家自然科学基金资助项目(19972029) |
| |
| 摘 要： | 建立了一种求解非线性动力学常微分方程组初值问题的新方法．若非线性函数一阶导数存在，则给出解的积分方程表达式，计算得到按规定误差要求的高精度数值解．引入一般自治或非自治非线性系统的首次近似Jacobi矩阵，不作任何假设重构等价的非线性常微分方程组，简捷而有广泛的适应性，不改变方程的本质，但其主项构成线性化方程组，其它项则代表非线性函数高阶余项而不涉及Taylor级数展开计算，给出该方程组初值问题的Duhamel卷积分解析表达式，在时间步长内进行数值积分选代求解，在指定误差内快速收敛，逐步递推获得非线性常微分方程的瞬态响应和全时域高精度数值解．积分解连续满足微分方程组而不是在离散的步长端点上满足代数方程组，打破了传统用增量法在离散点上建立的代数方程组迭代求解，从而使传统Euler型逐步积分法的各种差分格式算法改变成真正的积分格式算法．数值计算中给出指数矩阵递增展开式，变矩阵乘法为乘积系数的加法，避免了大量矩阵自乘而大大提高计算效率．算法验证为无条件稳定，则保证对线性常微分方程而言，计算中舍入误差的传播不会扩散，不出现计算机字长有限而引起舍入误差导致计算不确定性问题．基于以上理论和数值方法，计算了线性非线性算例并进行了分析，验证了本方法简捷而有广泛的适应性，可以有足够的精确性．
|
| 关 键 词： | 非线性动力学 有限差分法 常微分方程组 初值问题 高精度数值积分方法 高阶余项 |
| 修稿时间： | 2001年8月31日 |
ACCURATE SOLUTION OF NON-LINEAR DYNAMIC SYSTEMS BY NUMERICAL INTEGRATION METHOD |
| |
| Zheng Zhaochang Shen Song Su Zhixiao.ACCURATE SOLUTION OF NON-LINEAR DYNAMIC SYSTEMS BY NUMERICAL INTEGRATION METHOD[J].chinese journal of theoretical and applied mechanics,2003,35(3):284-295. |
| |
| Authors: | Zheng Zhaochang Shen Song Su Zhixiao |
| |
| Abstract: | |
| |
| Keywords: | non-linear dynamics FEM reconstructed equivalent nonlinear differential equation high-order remainder duhamel convolution integration |
| 本文献已被 CNKI 维普 万方数据 等数据库收录！ |
