解非线性方程的自动调节阻尼法.pdfVIP

( ) 　应用数学和力学 ,第 19 卷 第 2 期 1998 年 2 月 应用数学和力学编委会编 　　　　　　　　　　　 　 　　Applied Mathematics and Mechanics 重 庆 出 版 社 出 版 解非线性方程的自动调节阻尼法 ① ① 常海萍 　　黄太平 ( ) 钟万勰推荐 ,1995 年 6 月 20 日收到,1996 年 6 月 10 日收到修改稿 摘 　　要 　　解非线性方程组的一般方法是将其线性化 ,形成各种形式的迭代程序进行数值近似计算 ·　对 于复杂强非线性问题 ,在迭代过程中往往不易收敛 ,甚至数值失稳而发散 ·　不能满足工程要求·　 常规的牛顿法及改进的牛顿法均未彻底解决这一问题 ,因而使得复杂强非线性问题的数值模拟计 算受到了限制 ·　本文提出一种新的方法 ———自动调节阻尼法,是对带阻尼因子的牛顿法的进一步 改进 ·　引进阻尼因子向量,在迭代过程中 ,通过判断与调整 ,不断地自动调节阻尼因子向量 , 引用 有效收敛系数与加速系数 ,改善对赋初值的要求 ,加速求解的迭代过程 ,保证了复杂强非线性方程 求解的稳定性 ·　采用这一新的方法, 已成功地数值模拟了飞机中的一些复杂的传热问题 ,可进一 步推广用于非线性流动、传热、结构动力响应等各种复杂强非线性的工程问题的数值模拟计算 ·　 　　关键词 　非线性方程 　稳定性 　牛顿法 　自动调节阻尼法 　阻尼因子 §1 引　　言 工程上大量问题都属于非线性问题 ,如具有非线性辐射 —对流边界的 ,热容量有限的传热 问题,非线性流体力学问题 ,非线性结构动力响应问题等 ·　除了极为特殊的非线性方程组外, 直接解法几乎是不可能的 ,一般都是采用线性化的方法去构成各种形式的迭代程序 ·　然而在 迭代过程中往往会出现解的振荡 ,不易收敛 ,甚至发散的现象 ·　特别是一些较复杂的问题,如 计算燃烧室火焰筒壁温 , 由于筒壁很薄,壁面上有气流对流换热和高温燃气的辐射换热 ,当有 气膜冷却时 ,壁面上还存在热流方向相反的辐射换热和射流冷却 ·　这种情况下迭代过程极易 出现解的发散·　又如航空发动机的加力燃烧室筒体,飞机发动机短舱后部框架等的传热问题 , ( 都会遇到类似的情况 ·　常规的迭代法如牛顿迭代法及牛顿法的几种改进方案 带松驰因子的 ) 牛顿法和带阻尼因子的牛顿法 都因为对初始近似值的要求仍太严 ,收敛速度不快而告失败 ·　 长期以来 , 由于没有找到一种妥善解决的方法以致于在如预测具有非线性辐射 —对流边界的 传热等这类问题上 ,数值模似计算受到极大的限制 ·　本文提出了一种新的方法 ———自动调节 阻尼法 ,可以有效地用于强非线性问题的数值计算 ,处理比较复杂的问题 ,对赋初值的要求可 以放宽,迭代过程收敛速度加快[1 ] [2 ] [3 ] ·　 为了改善非线性问题数值计算中对赋初值的严格要求 ,为了防止由于受非线性参数影响 而在迭代过程中矩阵出现病态 ,引进阻尼因子向量 ,而且在迭代过程中 , 阻尼因子可以不断地 ① 南京航空航天大学 205 教研室 ,南京 210016 151 © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 152 常 　　海 　　萍 　　　黄 　　太 　　平