简单断面型钢热轧过程数值模拟.docVIP

J 1999年 第11卷 第3期 Vol.11 No.3 1999 简单断面型钢热轧过程的数值模拟 张　鹏　　鹿守理　　高永生 摘　要：以在椭圆孔型中轧制方坯为例，采用全三维弹塑性大变形热力耦合有限元法来模拟简单断面型钢的热轧过程，分析并建立了热传导、接触和摩擦等边界条件模型，计算得到轧件形状、温度变化及变形等结果。而且数值模拟结果与实验结果吻合较好。关键词：热轧，型钢，变形，有限元法中图分类号：TG335 Numerical Simulation for Hot Shape Rolling Processof Shape Steel with Simple SectionZhang Peng　　Lu Shouli　　Gao YongshengUniversity of Science and Technology Beijing Beijing 100083 ABSTRACT：A full 3-D elastoplastic thermo-mechanical coupled FEM for large strain is applied to simulate the hot shape rolling process from square to oval. Boundary conditions, such as heat transfer, contact and friction, are analysed and the models about them are established. The billet shape, temperature change and strain are obtained by calculation. The calculated results are in agreement with those from the experiments.KEY WORDS：hot shape rolling,shape steel,strain,FEM　　由于金属塑性变形过程的复杂性、影响参数的多样性和交互性，需要借助计算机模拟才能更好地揭示其本质和影响参数之间的关系，从而优化工艺过程，降低生产成本，缩短产品开发研制周期。有限元法是实现上述目标的最有利工具。目前，由于三维变形接触情况复杂，实体造型和网格划分技术要求较高，计算及存储量大，采用有限元法研究型钢热轧问题是比较有限的［1～5］。为了减少计算量，研究者对变形区内速度场和温度场都作了不同程度的假设。但是，在热轧过程，特别是在型钢的热轧过程中轧件内的变形、变形速率、温度等都存在显著的不均匀分布，从而导致组织性能的不均匀分布。因此，为提高模拟的精度，作者借助商业有限元软件包MARC/Autoforge，在分析建立热传导、摩擦以及材料模型基础上，采用弹塑性大变形热力耦合有限元法对方坯在辊道上散热和在椭圆孔型中轧制的过程进行解析。计算结果与实验结果相吻合，证明了该模型的正确性以及应用这种方法求解轧制问题的可靠性。1　有限元模型及边界条件　　由于轧件具有对称性，取1/4进行网格划分。总单元数为8×8×100=6 400 个。材料的热物理性能(热传导系数、热容和热膨胀系数)取为温度的函数，变形抗力取为变形量、变形速率和温度的函数。以上材料数据由MARC/Autoforge的材料库提供。1.1　热边界条件　　轧制时轧件的表面存在热传导、对流和热辐射三类热边界条件，对于热辐射和热对流情况，边界条件可统一写为： 　　　　　q=H(t-t∞) (1) 式中　q——热流密度，kW/m2；　　　H——等效热传导系数(H=h+hr)，kW/(m2.℃)；　　　h——对流换热系数［h=3×10-3(t-t∞)0.25］［6］,kW/(m2.℃)；　　　hr——辐射换热系数［hr=Eσ(t+t∞)(t2+t2∞)］,kW/(m2.℃)；　　　E——黑度(采用试错法［6］确定E=0.76)；　　　σ——波尔兹曼常数；　　　t，t∞——轧件表面和环境温度，℃。　　轧件与轧辊之间的热传导边界条件可描述为： 　　　　　　qc=hc(t-td) (2) 式中　hc——接触热传导系数，kW/(m2.℃)；　　　td——轧件与轧辊接触处轧辊的表面温度，℃。　　Chen W C［7］通过热轧板材试验，获得3个品种的钢在14种变形条件下的接触热传导系数值。本文作者对这些数据进行处理得到接触热传导系数回归模型，即：　　　　　　hc=0.696p-34.4 (3) 式中　p——接触表面压力，MPa。　　在轧件的上下、左右对称面，热流密度q满足如下条件：　　　　　　　　　q=0 (4) 　　与轧辊接触面的摩擦热为：