80Q77.板材轧制中的改进循环中点求积牛顿算法-东大数学系王海平.pdfVIP

2007 年全国塑性加工理论与新技术学术研讨会 2007 年 5 月 沈阳 * 板材轧制中的改进循环中点求积牛顿算法 1 1 2 2 王海平 ,宋叔尼 ,李长生 ,刘相华 (1 ．东北大学数学系 沈阳；2. 东北大学轧制技术及连轧自动化国家重点实验室 沈阳 110004) 摘 要：本文针对常用的 Newton-Raphson(N-R)迭代法对初始值限制较多,提出了一种改善初值的算法——改进的循 环中点求积 Newton 法: 即在 Newton-Raphson(N-R)迭代法与同伦法结合的基础上修改了原算法的步长和迭代终止的 判定标准. 实验结果表明,这一算法是可行的。 关键词：刚塑性有限元；初速度场；中点求积法；Newton-Raphson 迭代法 Improvement for Rotative Mid-point Mensuration and Newton-Raphson in Strip Rolling 1 1 2 2 WANG Haiping ，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: Newton-Raphson method has a strict limit to the initial value. An improvement mid-point mensuration is proposed. In this method, both pedometer and the termination of iteration in the primary alogithm are modified on the basis of the combination of Newton-Raphson(N-R) method and homotopy method. The results of experiment show that the improved algorithm is usefull. Key words: Rigid plastic FEM; Muzzle velocity field; Mid-point mensuration; Newton-Raphson method. 1 概述 随着现代科学技术和经济的发展，用户对钢铁产品的需求越来越趋于多样化、个性化和优质 化[1]。为了提高控制板坯轧制的模型精度，近几十年刚塑性有限元分析方法得到了广泛应用。不过， 计算效率上还存在一些问题，初始速度场的设定就是其中之一。如果处理不当，可能导致迭代过程 不易收敛或计算结果精度降低等问题。 数学上的最优化理论已经证明[2]，当目标函数的梯度和 Hessian 矩阵的表达式容易得到时，使 用 Newton-Raphson 迭代法求解一旦顺利，其收敛速度是很快的，条件是能给出一个较好的初始点。 然而，Newton-Raphson 迭代法对初始点又要求比较严格。为此，对本文提出了一种改进的循环中 点求积Newton 法。 求积法是数值分析中常用的算法，如 Euler 法、Runge-Kutta 法等。通过对工作量的研究，文献 [3]中给出了循环中点求积公式，并将它与 Newton 法相结合形成了一种新的算法循环中点求积 Newton 法。