收 稿日 期: 2003-10-14
基 金项 目: 国 家自然 科学 基金资 助项 目(59995440)
·
作 者简介: 时 旭(1975 - ),女,辽宁沈 阳人,东 北大学 博士 研究生; 刘 相华(1953 - ),男,黑龙 江 双鸭 山 人,东 北 大学 教 授,博 士生
导 师; 王 国栋(1942 - ),男,辽宁 大连人,东 北大 学教授,博士生 导师
·
第 25 卷第10 期
2 004 年 10 月
东 北 大 学 学 报 ( 自 然 科 学 版 )
Journal of Northeastern University(Natural Science)
Vol.25,No .10
Oct . 2 0 0 4
文章编号: 1005-3026(2004)10-0957-04
四 辊 冷 轧 机 轧 辊 弯 曲 和
压 扁 变 形 的 有 限 元 分 析
时 旭
1
, 刘相华
1
, 王国栋
1
, 李山青
2
(1 .东北大 学 轧制 技术及 连轧 自动化 国家重 点实 验室, 辽宁 沈 阳 110004; 2 .宝 山钢铁 股份 有限公 司, 上 海 201900)
摘 要: 借助 Marc 有限元软件,采用三维弹塑性有限元法对四辊冷轧机冷轧过程进行了模
拟,同时对轧辊变形进行了分析
·
计算模型中将辊系变形与带钢变形统一 考虑,并 解决了轧件 与辊
系之间的耦合问题,避免了采用假定或迭 代方法 确定轧 制力分 布时产 生的误 差
·
采用逐 步收敛 的
求解过程使计算结果精确、可靠
·
在不同 的轧制 条件下,得出 了带钢 宽度、弯辊力 等参数 对辊系 弯
曲、工作辊接触弧上的压扁变形、板宽方向的压扁变形和有载辊缝的影响,为板形 分析与控制 提供
了一种新的计算方法和参考数据
·
关 键 词: 冷轧;四辊轧机;弹塑性有限元;辊系变形
中图分类号: TG 335 文献标识码: A
冷轧产品向更薄、精度更高趋势的发展,促进
了板带材的高精度轧制和板形控制原理研究的深
入
·
板带材的高精度主要指板带材纵向和横向的
厚度精度,而轧辊的弹性变形直接影响到最终产
品的断面形状
[1]
,这就对轧机辊系弹性变形的计
算精度提出了更高的要求
·
采用二维有限元法计
算辊系变形
[2 ,3]
, 单元厚度的确定使计算模型与
实际的物理模型存在差异
·
三维有限元网格为所
研究的物理模型提供了最确切的表示
[4]
,被广泛
应用在轧制方面
[5 ~8]
·
其中刚塑性有限元法具有
计算效率高的特点, 广泛应用于热轧中
[9]
·
弹塑
性有限元可以更加真实地反映出轧制中的各种状
态,应用在冷轧中可以全面地考虑轧辊与带钢的
变形,但计算成本高
·
采用三维有限元方法研究辊系变形,普遍的
做法是将轧辊的弹性变形和轧件的塑性变形作为
两个独立的问题分别计算
·
但由于轧件与轧辊的
变形是其相互作用的结果,它们之间的耦合存在
着复杂的摩擦和接触,使得辊系与轧件之间的力
的关系很难获得,通常采取假设的方法
[1,4,10]
,降
低了结果的精确度
·
本文在大型非线性有限元分
析软件 Marc 提供的平台上进行轧辊弹性变形的
模拟
·
计算模型将辊系变形与带钢变形统一考虑,
解决了轧件与辊系之间的耦合问题
·
1 模型的建立
1 .1 几何建模及单元网格划分
按照轧制的实际尺寸建立几何模型
·
由于轧
制条件及几何参数的对称性,取四分之一作为模
拟对象,模型的几何参数及物理性能见表 1
·
在划
分单元网格时,为了减少单元数目,变形体的可能
接触区划分得较密,如图 1 所示
·
表 1 模型参数
Table 1 Parameters of model
类别 参 数 数 值
支
撑
辊
辊身直径/ mm 1 200
辊身长度/ mm 1 350
辊径直径/ mm 700
工
作
辊
辊身直径/ mm 470
辊身长度/ mm 1 510
辊径直径/ mm 285
带
钢
轧前厚度/ mm 0 A.8
宽度/ mm 700 ,960,1 230
轧辊和带钢均选用八节点六面体单元,采用
修正的 Lagrange 法描述的大变形弹塑性有限元
模型, Von Mises 屈服准则
·
摩擦条件选用 Shear
for rolling
·
用于轧制的“friction for rolling”针对轧