三辊行星轧制运动和管坯变形规律仿真模拟.docxVIP

三辊行星轧制运动和管坯变形规律的仿真模拟作者：李冰杨志刘化民张士宏摘要：本文应用弹塑性大变形有限元法，对连铸铜管坯三辊行星轧制过程进行了三维有限元模拟研究。通过坐标变换的方法建立了三辊行星轧制有限元模型，该模型考虑了轧辊的倾斜角β和偏转角α。模拟分析得出三辊行星轧制过程中管坯断面要经受一个由圆形到三角形再归到圆形的变形过程。获知在三辊行星轧制过程中起到主要作用的是三向压应力状态，有利于坯料的大变形和提高轧制管材性能。同时模拟分析获知三辊行星轧制过程中坯料质点的运动轨迹和相应规律。该模拟计算对于探索复杂的三辊行星轧制规律具有重要意义。关键词：三辊行星轧制；铜管；有限元分析；变形规律；应力状态1 引言三辊行星轧机（Planet Schr?ge Walzwerk , 简称PSW）最初是由德国施罗曼－西马克公司于1974年研制成功，用于轧制钢棒材，后应用到轧制钢管材。90年代初，芬兰奥托昆普公司将该型轧机应用到铜管材生产的铸轧法中，成功地轧制出完全再结晶的细晶铜管材，使“水平连铸－行星轧制”工艺成为现今典型的铜管铸轧生产线。然而，目前三辊行星轧制在模具和工艺设计等方面主要基于有限的经验和试错基础上，产品的组织性能和工艺缺陷不能有效控制，生产的稳定性和成品率还有待提高。90年代初期三辊行星轧制的理论研究集中于三辊行星轧管机（PSW）的结构和运动特点的分析，建立了轧制过程中轧辊的速度方程[1]；通过平面断面假设对三辊行星轧制力能参数进行理论推导计算[2]；在理论分析三辊行星轧制力和轧制功率的基础上，进行轧制实验研究轧出预期管材[3]等等。董顺德等对于轧机芯棒运动的控制及管坯宏观组织的变化进行了分析[4]。台湾的Chi-Kang Shih等利用LS-DYNA3D建立起三辊行星轧制钢棒材数值模拟模型，对钢棒材变形过程中的应力应变分布规律进行了分析[5]。Chi-Kang Shih等还利用塑性泥作为轧制实验原料对三辊行星轧制棒材进行实验模拟并就轧件变形尺寸、螺旋纹的间距和偏转角度等与模拟结果进行了对比分析[6]。本文作者对于铜管坯三辊行星轧制的有限元模型建立和咬入阶段变形情况已进行了一定的模拟分析[7]。由于三辊行星轧制的运动和坯料的变形流动非常复杂，无法直接观测，通过应用有限元模拟的方法来研究和再现三辊行星轧制这一复杂的变形过程有很大的生产实际和理论研究意义。本文结合生产实际情况，利用商业化CAE软件MSC.MARC建立起三辊行星轧制仿真模拟模型，对行星轧制稳态阶段进行了有限元模拟分析。2 三维有限元模型的建立2.1 三辊行星轧机模型的建立在有限元模拟过程中，空间轧辊辊形的建立是一个非常复杂的问题，需要建立起合理的空间坐标体系。由于实际生产过程中轧辊为3Cr2W8V合金钢，其相比于紫铜坯料强度大、硬度高，在轧制成形过程中只会产生相对很小的弹性变形，因而在模拟计算分析过程中可将轧辊定义为刚性。在建立轧辊模型时可用刚性的外表面来表示生产中应用的轧辊。建模过程中必须利用坐标变换的方法，建立起既有倾斜角β又有偏转角α的空间轧辊模型。坐标变换推导方程如下。对于三辊行星轧机来说，设备的结构参数即行星轮轴线至轧制线的中心距A是一定值，与轧制线距离为A的O1X1轴线与轧辊轴线O1XR的交点O1，为轧辊局部坐标系的中心点。当整体坐标系沿Y轴移动距离A，绕O1X1轴偏转角度α，再绕O1Z2轴倾斜角度β，最后所得到的局部坐标系O1XRYRZR，即为轧辊轴线所确立的局部坐标系统，如图1所示。图1 三辊行星轧制坐标变换图当已知轧辊局部坐标系上一点RR的坐标：通过坐标变换即可推导出整体坐标系与局部坐标系间的转换方程：其中D0R为转换矩阵，其中为α偏转角，β为倾斜角，通过转换矩阵D0R即可得到RR对应全局坐标系R0的坐标为：其中：根据如上方程作者建立起围绕轧制轴线互相成120°角的三个轧辊模型。三辊行星轧制有限元模拟模型如图2所示。在模拟过程中由于轧辊的公转带来了轧辊自转轴线的旋转，使模拟计算的接触问题更加复杂，并需要大量模拟计算时间，因此在模拟过程中使轧辊的自转轴线不变，这样轧辊的公转速度就转化为轧出铜管的旋转速度。坯料共划分成13320个单元，16418个节点。图2 三辊行星轧制有限元模型a)侧视图；b)正视图2.2 轧制工艺参数的确定计算中采用更新的Lagrange算法、Prandtl-Reuss流动方程以及von Mises屈服准则处理行星轧制过程三维大变形成形问题。选择8节点6面体等参单元划分有限元网格。轧辊与轧件之间采用剪切摩擦模型,铜管坯的泊松比0.35，密度8.94g/cm3,弹性模量110GPa, 初始屈服强度100MPa。其它相关轧制工艺参数如表1所示。表1 轧制工艺参数倾斜角β(°)偏转角α(°)轧辊转速rad/sec小车速度mm/sec轧辊与坯料间