板坯轧制的刚粘塑性有限元探析.pdfVIP

第8卷第3期 塑性工程学报 V01．8N0．3 1111 1塑 兰!’ !!!! 1 』Q堕塞盥垒垦Q旦!至△§!!￡!!Y垦堕鱼!堕曼曼星!堕鱼 板坯轧制的刚粘塑性有限元分析。 张晓明姜正义 ‘东北大学轧制技术及连轧自动化国家重点实验室沈阳110006’ 刘相华 王国栋 摘要：作者将变位的不完全乔里斯基分解刚度矩阵和共轭梯度法相结合，即变位的不完全乔里斯基分解和共轭梯 度法(shifted1CCG法)，求解速度增量方程，进而对板坯轧制过程进行刚粘塑性有限元分析。井对该算法的迭代 次数、变位参数廿进行了检测，同时对工艺参数的模拟结果进行了分析。验证了变位的不完全乔里斯基分解和共 轭梯度法在板坯轧制中的有效性和稳定性。 关键词：三维刚牯塑性有限元；变位的不完全乔里斯基分解和共轭梯度法；板坯轧制 —一整体速度场 l概述 ，～节点力向量 联立方程组的维数依赖于实际问题复杂程度和 为了提高控制板坯轧制的模型精度，近儿十年 空间离散化后的单元和节点的数量，它也依赖于用 来刚粘塑性有限元分析方法得到了广泛应用。一些 于单元和节点的编号方式。 研究者已经采用有限元法对板坯进行了平辊轧 如果A是对称正定矩阵．那么它就可以用共轭 法“2一和立辊水平辊轧制模拟-…。不过，由于三维 梯度法解。为了提高该解法的效率，矩阵A可通过 模型和其他假设条件过于简化，因而在精度和效果 不完全乔里斯基分解为预条件矩阵，然后与共轭梯 t还存在一些问题。 度法联合求解。这种方法对于大的稀疏系统显得十 本研究中，在三维刚粘塑性有限元法的基础上， 分有效[6-。93。如果A为一个M矩阵，那么它就可以 通过引入变位的不完全乔里斯基分解和共轭梯度法 分解成如下形式 对板坯轧制过程进行了分析，该方法考虑了高度压 A=M—N=UJ一^|(2) 下量、初始板坯宽度、立辊直径、摩擦条件和变位 式中L、u～—分别为下和上三角形矩阵 参数．』，的影响。研究了立辊直径的变化对轧制力的 N——对称矩阵 影响，压下量对轧制力的影响以及参数m对轧制力 如果A是对称的，那么M也将是对称的，此时 及迭代次数的影响。为了进一步检验采用变位的不 式(2)变成 完全乔里斯基分解和共轭梯度法进行刚粘塑性有限 A=CC7N (3) 元计算的有效性和稳定性，本研究还将其模拟结果 式中．T表示转置，此分解可与常用的共轭梯 与用牛顿一拉普逊法得到的结果进行了比较。 度法联合求解具有M矩阵方程组“。 在金属成形过程分析中，由有限单元的离散化 2变位的不完全乔里斯基分解和 而得到的刚度矩阵A是一个大的稀疏矩阵”j，R并 der Vorst”o， 共轭梯度法运算法则 非一M矩阵。通过对Me