实心坯二辊斜轧过程三维热—力耦合分析.doc

实心坯二辊斜轧过程三维热—力耦合分析 第22卷第1期北京科技大学学报 2000年2月JournalofUniversj廿ofScienceandTechno|ogyaei;ing Vb1.22NO.1 Feh.2DDD 实心坯二辊斜轧过程三维热一力耦合分析 .zI一李胜祗”孙中建拿连诗t/巧/7 }华东冶金学院清金1系,j鞍h2430022)北京科技大学材料乖l学与_ll程学院,北京100083 摘要用MARC/Autoforge程亭埘一辊斜轧穿4Ltf圃管坯在穿孔准备区的轧制过程进行 _I维弹单性热一耦台分析得到了,应变场温度分布,证实丁物理模拟实验得到的变形强 度汾管坯断血的w}分.彤态ji详细分忻了w形分的成因;澄清了孔腔形成理论I.的一 分歧毁现了帮的分形态:管坯中心存在着较大的正值静水压力. …”熟一力 关于:=辊斜轧实心圆,时变强度分布, 卢于述等--进行过研究,仇L述研究都难反 映斜轧变形的复杂过程例如,轧什螺旋前进的 运动学特点,斜轧变形区的几何特征,导『n工具 的作削效粜,更无法反映轧辊倾角t包括送进角 和辗轧角)的影响,故其结论尚卟能代表圃坯斜 轧的真实情况.将圆坯板缩,坯的横锻和 横轧实验结果应用到斜轧穿孔过程显然是不恰 当的,也影响r人们对中心孔腔形成本质的认 识. 本文采用MARC/Autoforge软件,在完全实 际的参数条件下模拟吲坯二辊斜轧过程,以揭 示穿孔顶头前金属的真实应力与变形状态,深 化对内撕裂形成机理的认识. 1数值模拟的模型与边界条件 1.1孔型构成 实心坯斜轧孔型由轧辊和导盘构成,如图 1所示. 变形工具的主要尺寸,转速及调整参数(取 自宝钢集公司140连轧管机组狄舍尔穿-~L,F)L J实际参数j如卜: 轧辊L桶形)直释Dr-!ln.入口锥角F2.5., 【{JU锥角-2.75.,转速HFI20r/mln;导盘直径 D18Ill,转速n~-26r/min,送进角a-12.,辗轧角 0.;轧辊距离154.6mm,导盘距离 L~-173I1]in 9994)4.09收稿李胜祗.44.副教授,陴』 暑垤 1.2圆坯尺寸,材质及单元划分 圆坯直径=l78mm,村质为以C22,相当 十中国的20#钢.圆坯长度的取定,一方面麻保 O.5L▲ .- 图1斜轧变形区构成示意图 FiglSketchofdeformationzoneconfigurationfor rotaryrolling 证轧件充满变形区并稍有赢余,另一方面要考 虑缩短计算时间,这里取Lb=225mill,用1800个 六面体单元离散轧件,横断面h单元数为120 个,长度方向单元数为15个.材料的热物理性 能(导热系数,热容,线膨胀系数)是随温度变化 的,变形抗力是热力学参数(变形程度,变形速 度和变形温度)的函数,这些材料数据均由 MARC/Autoforge的材料库中读取. 1.3初始条件和边界条件 轧辊和导盘均视作恒温刚性体,网坯开轧 温度为1250℃,边界条件如: (1)位移边界条件.圆坯中心所有节点的位 移”=酢=0. (2)摩擦条件.变形工具与轧件接触面间遵 循库仑摩擦定律,轧辊与轧件间摩擦因数 Vol~2No-1乍肚祗等:实心坯辊斜轧过一维热耦台分析 = 04,导盘与轧件问摩擦数=0.3. (3)传热边界条件.求解这类问题的有关理 论描述详见文献[4】,求解本问题的边界条件见 文献[5].轧制时轧件的白Fh表面1竽在热传导,对 流和热辐射类热边界条什,对热辐射和对流 换热情况,边界条件可统一表小为: :T—)(1) 式中:q一热流密度:等效导热系数,=+甩: 一 对流换热系数,这里取h=0.02kW/(m:℃);h 一 辐射换热系数,由下式计算: h:EoIT~咒)(+)(2) 黑度,这里取E=06;波尔然曼常数;, 一 分别为轧件表面和环境热力学温度. 接触热传导可描述为: = h(T-)(3) 式中:～接触导热系数,这里取 20kW/(m”‘℃);分别为轧件表面和与之接 触的变形_『_具表面热力学温度. 接触面卜的摩擦边界条件为: qr-F~-AV『41 式中:一接触面上摩擦力:△接触面卜相对 滑动速度. 模拟中忽略了圆坯端部与环境的传热.对 于变形热效应的影响,Werthehner认为对大多 数金属而言,变形功的90%转化为热量. 2模拟结果的讨论 在Pentium200双CPU微机进行模拟,以 F就模拟得到的结果进行分析讨论. 2.1应力和应变分布 国2为数值模拟得到的存不同横截面(Z) 处,各正应力分量沿直径(与轴方向一致)的分 /.r}f /:0!\ ,3,,0 ,\”, —— 3 , ,I\| ]/rl f\\ ill, f-ll 布,然在靠近轧什与轧辊接触面处为横向 J卡席力,而在圆坯中心区域为横