82Q79.板材轧制中的信赖域和牛顿混合迭代法-东大数学系尹丽萍.pdfVIP

2007 年全国塑性加工理论与新技术学术研讨会 2007 年 5 月 沈阳 * 板材轧制中的信赖域和牛顿混合迭代法 1 1 2 2 尹丽萍 ,宋叔尼 ,李长生 ,刘相华 (1 东北大学数学系，沈阳；2 东北大学轧制技术及连轧自动化国家重点实验室，沈阳) 摘 要：本文提出了一种加快刚塑性有限元的迭代收敛的新算法——信赖域和牛顿混合迭代法. 模拟计算结果表明, 该算法可满足计算要求, 能达到和 Newton-Raphson 法相同的收敛指标且有效提高计算速度. 关键词：刚塑性有限元; Newton-Raphson 法; 信赖域方法 Combining Trust-region and Newton-Raphson Algorithm in Strip Rolling 1 1 2 2 YIN Liping ，SONG Shuni ，LI Changsheng ，LIU Xianghua (1 Department of Mathematics, Northeastern University, Shenyang, 110004, 2 State Key Laboratory of Rolling and Automation, Northeastern University, Shenyang 110004) Abstract: A combining trust region and Newton-Raphson method for speeding up the convergence of iteration in rigid-plastic finite element analysis was discussed in the paper. The simulation results show that the calculation requirement can be assured. As well, the new method can reach the same convergence point as Newton-Raphson method and well speed up the calculation time. Key words: rigid-plastic finite element method; trust region method; Newton-Raphson method. 1 前言 迄今有很多非线性迭代算法运用于刚塑性有限元分析，其中包括简单迭代法和牛顿迭代法等。 众所周知，牛顿迭代法收敛很快，但缺点是它只具有局部收敛性。在各种改进牛顿法中， Newton-Raphson 法放松了对初始值的要求，改善了收敛域，但是收敛速度减慢了；并且每次迭代过 程中阻尼因子α是通过多次函数值比较而确定的。因而当方程组的维数较大时，仅确定α就需要相 当大的工作量。 信赖域方法[1]是非线性优化的一类具有全局收敛性的数值计算方法。自 Powell 首次提出后迅速 发展,已经和传统的线性搜索方法并列成为非线性规划的两类主要数值方法。组合信赖域与线性搜索 算法的思想在文献[2~4]中提出。本文介绍一种新的改进的信赖域方法——信赖域和牛顿混合迭代 法。先利用信赖域的全局搜索能力，减弱程序对初始点的要求，而靠近极小点的时候，信赖域法和 牛顿法奇偶交替进行求解，用局部寻优能力强的牛顿法对信赖域进行修正。 2 新算法的提出 2.1 信赖域的基本思想 在当前迭代点xk 的附近用一个二次函数逼近目标函数f . *基金项目