带材轧制过程应力及变形的计算机仿真-机械工程学报.PDF

第 40 卷第 9 期 机 械 工 程 学 报 Vol.40 No.9 2004 年 9 月 CHINESE JOURNAL OF MECHANICAL ENGINEERING Sep. 2 0 0 4 带材轧制过程应力及变形的计算机仿真* 彭 艳 刘宏民 （燕山大学机械工程学院 秦皇岛 066004 ） 摘要：用三次样条流线条元法分析带材的三维塑性变形，用分割模型影响系数法分析辊系的弹性变形，用 BP 神经 网络方法对轧后带材板形进行模式识别，将三者联立，对900 mm HC 轧机冷轧带材轧制过程进行了计算机仿真， 得到了较为详细的应力及变形的仿真结果。仿真实例表明，提出的方法和模型符合实际，对三维轧制理论、板形 理论及控制技术的发展有重要意义。 关键词：条元法 轧制 仿真 HC 轧机 中图分类号：TG333 变形。根据板带轧制的特点，按照金属的流动迹线， * 0 前言 将变形区划分为许多流线条元；用康托洛维奇法构 造条元横向位移函数：沿纵向(轧制方向)用已知的 研究带材轧制过程中轧件和轧辊的应力及变形 函数表示，沿横向(板宽方向)用精度高、拟合能力 对了解轧机板形控制特性，发展轧制理论、板形理 强的三次样条插值函数表示；根据体积不变条件及 论及控制技术都有重要意义。在目前研究冷轧带材 几何关系导出流动速度和应变速度；根据 三维变形问题的方法中[1~5]，条元法[5~7]是一种很实 Levy-Mises 塑性流动方程和Von-Mises 塑性条件导 用的工程方法，能解决大宽厚比(大于 1 000)的带材 出三向应力模型，应用能量变分原理和优化方法确 轧制问题。这里采用三次样条流线条元法[6,7]分析仿 定条元节线上的待定位移参数[6,7] 。流线条元法迭 真带材三维塑性变形。在分析过程中，采用条元变 代计算的计算精度和收敛速度，取决于迭代计算 分法[8]确定三次样条流线条元法迭代计算的初值， 的初值选取是否合理及条元法数学模型本身的准 大大提高了计算收敛速度。 确性。 为了能使轧制过程中的应力及变形计算仿真结 1.1.1 条元变分法确定出口横向位移的初值 果与实际相符，必须将轧件三维塑性变形模型(轧件 将整个带宽划分为n 个纵向条元，如图 1 所示。 模型)与辊系弹性变形模型(辊系模型)耦合计算。轧 节线横向坐标用y (i 0,1,2,L,n ) 表示，节线上出口 i 件模型为辊系模型提供单位宽度轧制压力的横向分 横向位移用 u u( y ) 表示。单元宽度 b y − i i i i 布，为板形模式识别模型提供前张应力的横向分布； y i−1 ( i 1,2,L,n ) 。 辊系模型采用分割模型影响系数法建立[6,9]，为轧件 模型提供轧后带材厚度的横向分布；板形模式识别 模型采用 BP 神经网络方法，根据前张应力分布， 对轧后板形进行模式分类和提取。将这三个模型联 立，对 900 mm HC 轧机冷轧带材轧制过程进行了计 算机仿真，得到了轧件出口横向位移、前后张应力、 出口厚度、轧制压力和辊间压力的横向分布，以及 轧件内部应力应变在变形区分布规律等仿真结果。