IMR OpenIR
板材轧制过程刚塑性有限元快速算法研究
其他题名Research on fast arithmetic of rigid plastic FEM in strip rolling
张光亮
学位类型博士
导师张士宏
2008-12-23
学位授予单位中国科学院金属研究所
学位授予地点金属研究所
学位专业材料加工工程
关键词板材轧制 刚塑性有限元法 初始速度场 线性搜索 并行算法 快速算法
摘要基于传统轧制理论的在线数学模型计算精度已经得到了最大限度的发挥,人们将进一步提高在线控制模型精度的期望寄予了人工智能和有限元等方法。有限元法已经被广泛应用于板材轧制过程离线分析,然而由于计算时间长不能在线应用。开展有限元快速算法的研究,实现有限元的在线应用具有很大挑战性以及重要的理论价值和实际意义。本文进行了板材轧制过程刚塑性有限元在线快速计算理论的研究,实现轧制力、轧制力矩、轧制功率和前滑值等轧制参数的高精度快速计算。 首先,建立了板材轧制过程刚塑性有限元快速计算基本理论。简化了有限元模型,只对轧制接触区下方划分了有限元网格。采用反正切摩擦模型解决轧制中性点问题。采用在速度不连续面附加剪切塑性变形能耗率项的方法解决第一类速度奇异点问题。采用建立统一塑性变形能耗率的方法处理刚性区问题。开发了快速计算刚塑性有限元程序FAR-2D ParL,分析了有限元计算时间的分布规律,确定把最耗费时间的Hessian矩阵的集成和一维线性搜索作为提高求解速度的主要研究方向。 针对刚塑性有限元Newton-Rapshon迭代求解过程,初始速度场越接近真实速度场迭代越容易收敛的特点,开展了加快收敛的初始速度场设定的研究。讨论了传统设定初始速度场的工程法和G函数法,并提出了神经网络法。为建立准确和高效的神经网络模型,对输入数据进行了相关性分析,去除了相关性不明显的输入参数,减少了变量的个数。神经网络模型用动量修正的反向传播算法对刚塑性有限元程序提供的大量训练样本进行学习获得。三种初始速度场设定方法进行了比较,结果表明神经网络法设定的初始速度场和初始总能耗率泛函都非常接近真实值,计算只需少量迭代步即可收敛,大量提高了有限元的求解速度,为实现有限元的快速计算提供了较好的初始速度场设定方法。 开展了加快阻尼因子线性搜索算法的研究。首先,引进了黄金—抛物法搜索包含极小值点的初始区间。其次,引进了缩小搜索区间求解阻尼因子的数学一维最优化算法:Fibonacci法和不求导数的Brent法。与传统的试探法搜索初始区间和黄金分割法缩小区间的线性搜索算法进行了比较,结果表明黄金—抛物法搜索初始区间和不求导数的Brent法缩小搜索区间可以明显减少能耗率泛函的计算次数,采用这两者组合得到的黄金—抛物加Brent法减少了线性搜索所消耗的计算时间,提高了有限元求解过程的计算速度。 进行了刚塑性有限元求解板材轧制过程并行算法的研究。提出了刚塑性有限元Newton-Raphson迭代过程Hessian矩阵和能耗率泛函求解的共享存储并行算法、板材连轧过程多道次分布式存储并行算法和多道次分布式共享存储混合并行算法。采用现场多道次连轧数据测试各种并行算法的性能,结果表明Newton-Raphson迭代过程并行处理中 Hessian矩阵的并行化可以获得很高的加速比,能耗率泛函的并行化提高有限元求解速度较弱。板材连轧多道次并行计算表现出很好的加速作用,加速比接近进程数。多道次分布式共享存储的混合并行计算时,即在道次内Newton-Raphson迭代采用共享存储并行和道次间采用分布式存储并行,可以获得大于进程数的加速比,是提高有限元计算速度效果最好的方式。 进行了提高刚塑性有限元求解速度的软硬件平台优化测试,提出了板材轧制过程有限元在线快速计算方案。为提供高效的有限元在线快速计算软硬件平台,对操作系统、编程语言和软件以及硬件等进行了优化测试。提出了提高轧制力计算数值稳定性的能量法。建立了4种板材轧制过程刚塑性有限元在线快速计算理论方案,并采用板材生产现场7道次连轧数据进行了测试。结果表明4种快速计算方案中最快的方案完成7道次的计算时间小于60ms,最慢的方案也小于190ms,均能满足板材连轧过程在线设定对计算时间的要求。预测轧制力和实测轧制力的比较体现了较高的预测精度,可以达到提高在线控制模型设定精度和适用性的目标。
其他摘要The precision of the mathematical models based on the traditional rolling theories has reached its limited point. In order to improve the precision of the online model, most researchers focus on the artificial intelligence and the finite element method (FEM). FEM has been widely used for the offline analysis in strip rolling, but it is unavailable for the online application due to the long computational time involved. The research on the fast algorithm of FEM for online application is a great challenge task with important theoretical value and practical significance. In this work, the fast online calculation theory for rigid plastic FEM (RPFEM) in strip rolling is researched, and then the fast online calculation for the various rolling parameters is carried out including rolling force, rolling torque, rolling power and forward slip etc.. Firstly, the basic fast calculation theory of RPFEM for strip rolling is established. The finite element model is simplified by only meshing the deformation zone of the workpiece under the contact area with roll. The arc tangent frictional shear stress model is used to deal with the neutral point problem in strip rolling. The shear deformation work rate on the discontinuous velocity sections is adopted in the total energy functional of RPFEM to resolve the first singular point. The fast calculation software of RPFEM in strip rolling FAR-2D ParL has been developed, and the distribution rule of computational time is investigated, as a result the major ways to increase the solution speed of RPFEM are fixed on the most time-consuming parts including the integration of the Hessian matrix and the one dimensional line search. During the Newton-Rapshon solution of RPFEM, the closer the initial guess approaches the real velocity field, the easier the convergence can be obtained. Thus the setup of initial guess to speed up the convergence rate is studied systematically. The traditional methods including the Engineering Method and G Functional to set up the initial guess are discussed, and the Neural Network method is newly proposed. In order to establish an accurate and efficient neural network model, a correlation analysis of the input data is performed to remove the insignificant input parameters, thus the number of the input variables for the Neural Network model is reduced. The neural network model has been trained using the momentum amended back propagation algorithm from a huge amount of samples calculated by the rigid plastic finite element code. Comparisions for the three methods to set up the initial guess are performed, and the results show that the initial guess and the initial energy functional are very close to the real values, thus FEM can get the convergence through several steps iteration. Therefore, the solution speed of FEM can be increased greatly by using the Neural Network model; this provides a better method to set up the initial guess for the fast calculation of RPFEM. A fast line searh algorithm for the relaxation factor is developed. Firstly, the Golden –Parabola method is introduced to bracket an initial interval containing the minimum point. Secondly, the one dimensional optimization algorithms including the Fibonacci’s method and the Brent’s method without derivatives are introduced to reduce the interval and to obtain the relaxation factor. Compared with the traditional trial method for bracketing the minimum point and the Golden Section line search, the new fast line search using the Golden–Parabola method and the Brent’s method presents great advantages of reducing the frequency for evaluating the energy functional. By adopting the new fast line search algorithm, the computational time of the line search can be reduced greatly, thus the solution speed of FEM can be increased remarkably. Parallel algorithms for the rigid plastic finite element solution of strip rolling are developed. The computation for the Hessian matrix and the energy functional during the Newton-Raphson solution of RPFEM is parallelized using the shared memory parallelism. During a multi-pass solution for continuous strip rolling, a multi-pass parallel algorithm with the distributed memory parallelism is developed, and a multi-pass mixed parallel algorithm by adopting the shared memory parallelization for the Newton-Raphson solution within a pass and the distributed memory parallelization across passes is also proposed. The parallel performance has been tested by the numerical examples of the practical rolling data for the multi-pass continuous strip rolling. The results show that the parallelization for the Hessian matrix can obtain high speedup, while the parallelization for the energy function presents a small ability to increase the solution speed of FEM. The multi-pass distributed memory parallelization presents a great effectiveness to speed up the solution of RPFEM; the speedup approaching to the number of processes can be obtained. The multi-pass mixed parallel algorithm is the best way to speed up the solution of RPFEM, the speedup exceeding the number of processes can be obtained. The opitimization tests for the computer software and hardware platform to increase the solution speed of RPFEM have been performed, and the theoretical schemes of the fast online finite element calculation are proposed. In order to obtain a high efficient computer software and hardware platform for the fast online finite element calculation, the optimization tests have been performed on the operation system, programming language and compiler software, hardware etc.. An energy method for calculating the rolling force with a high numerical stability is newly proposed. Four theoretical schemes for the fast online calcaulation by RPFEM are constructed, and the seven-pass continuous strip rolling data from the production line have been used to test the feasibility of the four theoretical schemes. In order to finish the computation for seven passes, the fastest scheme consumes less than 60ms, the the slowest scheme consumes less than 190ms, thus the computational time for the four schemes can meet the requirements of online control. The predicted rolling force is compared with the measured values and presents a high precision; hence, the online application of FEM can improve the precision and flexibility of the online control model.
页数148
语种中文
文献类型学位论文
条目标识符http://ir.imr.ac.cn/handle/321006/17225
专题中国科学院金属研究所
推荐引用方式
GB/T 7714
张光亮. 板材轧制过程刚塑性有限元快速算法研究[D]. 金属研究所. 中国科学院金属研究所,2008.
条目包含的文件
条目无相关文件。
个性服务
推荐该条目
保存到收藏夹
查看访问统计
导出为Endnote文件
谷歌学术
谷歌学术中相似的文章
[张光亮]的文章
百度学术
百度学术中相似的文章
[张光亮]的文章
必应学术
必应学术中相似的文章
[张光亮]的文章
相关权益政策
暂无数据
收藏/分享
所有评论 (0)
暂无评论
 

除非特别说明,本系统中所有内容都受版权保护,并保留所有权利。