转炉扭力杆缓冲装置三维非线性有限元分析.doc

1、中国冶金装备网中国冶金人的网转炉扭力杆缓冲装置三维非线性有限元分析童靳于 1 潘紫微 1 包家汉 1 潘天成 1 曹天明 2摘要 针对某厂扭力杆缓冲装置运行时出现的问题，运用有限元软件建立其整体装配非线性有限元模型，充分考虑螺栓预紧力、零件间过盈配合及接触非线性，分析扭力杆缓冲装置在最大倾动力矩与事故状态下的应力分布及强度，提出螺栓联接修改建议。研究结果表明：螺栓根部出现严重的应力集中现象，这与螺栓过早损坏是一致的；扭力杆在最大倾动力矩作用下安全裕度较大；在事故状态下最大应力出现在扭力杆键槽附近，其值接近材料疲劳强度极限。研究结果为螺栓联结的改进提供了依据，也为扭力杆的设计制造与检修提供理论依

2、据。关键词 扭力杆缓冲装置；整体装配模型；非线性有限元；预紧力；接触应力The Analysis of Converters Torsion System Based On 3-D Nonlinear FEMTong Jin-yu1, Pan Zi-wei1, Bao Jia-han1, Pan Tian-cheng1,Cao Tian-ming2(1.Mechanical Engineering College, Anhui University of Technology, Maanshan, Anhui 2303222，China；2.Maanshan Iron & Steel Comp

3、any Limited, Maanshan, Anhui 2303222，China)ABSTRACT Aimed at some problems about torsion system of a factory, the 3-D nonlinear FEM whole assembly model is established. The influence of bolt pretightening load 、the interference fit and the nonlinear contact is full considered. The stress distributio

4、n and strength is analyzed under the maximum tilting moment and the emergency condition and the proposal of bolt connection is made。It is found that stress concentration phenomenon exists in the root of bolt and it is in coincidence with the pre-mature destruction of it. The safety margin under the

5、maximum tilting moment condition is large, and under the emergency condition is approach the ultimate fatigue strength. Research Result provides a theoretical basis for improving the bolt connected and for designing and repairing torsion bar.KEYWORDS Torsion system Whole assembly model Nonlinear FEM

6、 Pretightening load Contact stress引言扭力杆缓冲装置是目前国内大型转炉倾动机构普遍采用的平衡缓冲装置，用来平衡转炉倾动时引起的悬挂减速机壳体的旋转力矩，其性能直接影响到转炉设备的正常工作，扭力杆缓冲装置的分析涉及到螺栓、螺柱的预紧，扭力杆与联接板间过盈配合，销与扭力杆和联接板间过盈配合，各零部件间的高度非线性接触。因而对扭力杆缓冲装置各接触界面上接触应力的精确研究,在理论上还存在较大的困难。但随着计算机技术、数值计算方法以及大型有限元软件的发展,利用非线性有限元软件对扭力杆缓冲装置使用中的高度非线性问题进行研究成为可能。纵览国内关于扭力杆的限元分析,虽然各自的

7、研究都具有一定的现实意义,但所作的模型距离实际的工况仍有相当的差距。主要表现在以下几个方面: 为分析方便,假设的力学模型仅为一个阶梯轴,且承受的载荷简化过于单一；没有考扭力杆缓冲装置各零部件间的非线性装配状态、接触运动状态及摩擦力、预紧力、转炉倾动力矩对其造成的影响，而实际上这些影响对扭力杆的应力是不可忽略的。某厂 300T 转炉扭力杆缓冲装置投入使用期间，其上下部联接座连接螺栓出现断裂，为了对现役扭力杆缓冲装置进行评估，并为同类产品的设计制造与检修提供理论依据，需要全面了解扭力杆缓冲装置的应力分布、变形，确定扭力杆的强度。本文采用大型通用有限元分析软件 ANSYS 对扭力杆缓冲装置整体装配模

8、型进行三维非线性有限元分析。1 扭力杆缓冲装置有限元分析1.1 几何模型的建立扭力杆缓冲装置主要由扭力杆、联接板、上下部联接座、轴承座、球铰、销、键、销轴、事故支座组成。扭力杆一方面通过球铰连接其两端的轴承座从而固定在基础上，另一方面通过键、销、销轴、球铰连接两端的联接板和上下部联接座，从而借助两对上下部联接座与二次减速机壳体相连1。根据图纸建立扭力杆缓冲装置的三维几何模型，为了能准确反映实际情况且达到研究目的,在合理的条件下,对某些次要因素和对分析结果影响微小的因素予以忽略,如某些小圆角、倒角用直线代替。考虑到扭力杆缓冲装置的对称性，只建立一半模型，如图 1 所示。a) 装配图 b) 爆炸图

9、图 1 扭力杆缓冲装置模型（一半）Fig.1 model of torsion system此外，根据该厂螺栓实际情况，螺栓根部未加倒角。扭力杆中部半径为 R1=190mm，锥形端最大半径为 R2=215mm，联接板作用面半径 R3=260mm，轴端半径为 R4=150mm。1.2 有限元模型的建立尽管扭力杆缓冲装置结构是对称的，但出钢侧和加料侧所受载荷不对称，且扭力杆缓冲装置中涉及到的零件数目多且各零件间均为非线性接触，所以为了完全模拟结构的真实状态，采用 4 节点三维实体单元 solid45 建立扭力杆缓冲装置整体的非线性有限元模型，如图 2 所示。为了使计算结果更加可靠，运用弹性力学和接

10、触力学知识，对接触部位网格进行人工干预，关键结构采用六面体网格划分。扭力杆缓冲装置有限元模型共有 141016 个单元，其中扭力杆部分单元数为 61868，占 43.87%，共有节点数 66517 个。扭力杆材料为 34CrNi3Mo，联接板、销、螺栓和螺柱材料均为 42CrMo，轴承座和轴承盖材料为 35CrMo，均为合金钢，所以该模型中所有零部件选取的弹性模量均为 E=2.1e5MPa，泊松比 。3.0图 2 扭力杆缓冲装置有限元模型Fig.2 The infinite element model of torsion system1.3 边界条件对于扭力杆及其相关零部件结构，采用扩展的拉

11、格朗日算法和库仑摩擦模型计算摩擦接触问题。扩展的拉格朗日算法为了找到精确的拉格朗日乘子，需要对罚函数修正项进行反复迭代，与罚函数方法相比，拉格朗日方法不易引起病态条件，对接触刚度变化的灵敏度较小。在基本的库仑摩擦模型中，两个接触面在开始相互滑动之前，在接触的界面上会产生达到某一量值以内的剪切力，这种状态叫做粘着状态2。库仑摩擦力定义了一个等效剪切力，一旦剪切力超过此值，两个表面之间将开始相互滑动，变为滑动状态。库仑摩擦的切向力为： 滑 动 状 态粘 着 状 态nttfKf式中 Kt 为切向刚度， 为接触节点相对于目标表面的弹性变形， 为滑动摩擦系数。本次计算取滑动摩擦系数为 0.15。在两端轴

12、承座底部施加垂直于底部即 z 方向的约束，限制其 z 方向的自由度；在地脚螺栓孔处，建立局部柱坐标系，对孔表面节点施加径向约束以限制其径向位移，这样整个轴承座就被固定了，对于扭力杆缓冲装置其余部分而言，可以通过接触对的设置来约束零件间的相对运动。扭力杆缓冲装置的接触类型为柔-柔接触，接触面单元类型选取四边形单元 CONTA174，目标面单元类型选取无中间节点的三节点单元 TARGE170，共建立了 58 个接触对。扭力杆与联接板为 H7/u6 的过盈配合，销与扭力杆配合为 H7/p6 的过盈配合，对于这两处过盈配合，如果仅用两个接触部件几何形状的偏差作为相应的过盈量来模拟过盈配合是不够准确的，

13、因为划分网格后几何上的过盈量并不等于实际过盈量，二者只是近似相等。所以这里采用了设置接触对选项的方法来模拟，在第一个载荷步中设置渐变的初始穿透，设置接触单元关键选项 KEYOPT(9)=4 来渐变施加接触面穿透，忽略由于几何模型造成的穿透，这样设置控制面-面接触的实常数 CNOF的值即为所需施加的过盈量。1.4 载荷的施加在螺栓、螺柱上施加预紧力会对变形和应力带来很大影响，为了真实地反映扭力杆的应力，采用预紧力单元 PRETS179 和 PSMESH 预拉伸网格命令对扭力杆缓冲装置的螺栓、螺柱施加预紧力F=1374690N；由于上下部联接座一个受拉，一个受压，故在两端销轴 1 与二次减速器相连

14、的位置施加大小相等，方向相反的载荷 P=M/L，其中 M 为转炉的倾动力矩，L 为扭力杆出钢侧与加料侧两个联接板之间的距离 5200mm。由根据设计图纸知，在最大倾动力矩为 6500KNm，则P=1.25103N，在事故状态下 P=3.75103N。2 有限元分析结果米塞斯等效应力是根据第四强度理论确定的（即畸变能密度理论） ，其值为，该理论认为形状改变比能是引起材料塑性变形破坏的主要原因。本文224131()()s通过有限元的后处理功能，提取扭力杆缓冲装置关键零件的 Misses 应力。2.1 螺栓应力分布由于出钢侧和加料侧螺栓的应力分布大致相同，仅给出出钢侧螺栓的应力分布如图 3 所示。

15、图 3 螺栓 Misses 应力(出钢侧)Fig.3 The Misses stress of bolt（Tapping Side）螺栓的材料为 42CrMo，在变载荷作用下受轴向载荷紧螺栓连接螺栓的许用拉应力：lp=s/n（其中 n 为螺纹预紧联接的安全系数，控制预紧力时 n=1.21.5；不控制预紧力时 n7.5） 。所以对于某厂该处螺栓联接，其许用拉应力 lp 应满足 lp124MPa。而由图 3 可以看出，联接板连接螺栓的最大应力 =376MPa，出现在螺栓根部，且该处存在严重的应力集中现象。这与该处螺栓过早损坏是一致的。建议更换为双头螺柱以后，没有再出现过类似断裂问题。2.2 扭力杆

16、应力分布在最大倾动力矩作用下，扭力杆的应力分布如图 4 所示。a）扭力杆缓冲装置整体应力分布 b）扭力杆整体应力分布c）扭力杆轴头应力 d)扭力杆联接板作用段e)锥形段和光轴段 f)光轴段图 4 最大倾动力矩作用下扭力杆 Misses 应力Fig.4 The Misses stress of torsion bar under the maximum tilting moment condition扭力杆米塞斯等效应力最大值出现在扭力杆联接板作用段，因为联接板作用段上除承受较大的弯曲应力外，由于联接板上开有键槽和销孔，且其与销与、联接板间均为过盈配合，故该段还承受了较大的扭转剪应力，尤其是在销

17、孔的表面、键槽的侧面剪应力较大。应力较大值出现在扭力杆轴头处，因为该处存在严重应力集中且承受较大拉-压应力及一定的弯曲应力。锥形段及光轴段应力较小且整个光轴段的应力分布为：由外表面向芯部逐渐减小，芯部应力值接近零，所以扭力杆光轴段可近似地视为等直圆杆纯扭转。在事故状态下扭力杆应力分布与在最大倾动力矩作用下基本相同，但整体应力水平有所提高，两种工况下扭力杆米塞斯应力值见表 1 所示。表 1 两种工况下扭力杆米塞斯应力值最大倾动力矩作用下应力（MPa）事故状态下应力（MPa）联接板作用段 221.1 681.8轴头 152.4 425.7锥形段 99.1 298.8光轴段 95.308 299.7

18、由于扭力杆作用在循环交变应力下，所以要保证其正常工作必须满足其强度小于材料的疲劳强度极限，扭力杆的材料 34CrNi3Mo 对应的疲劳强度极限为 s=686MPa。由以上分析可知扭力杆在最大倾动力矩作用下的安全系数为 3.1；在事故状态下的安全系数为 1。扭力杆属于重要零件，由3知其许用的安全系数为：1.22.5。故扭力杆在最大倾动力矩作用下安全裕度较大；而在事故状态下安全裕度较小，此时如果再增加转炉倾动力矩，二次减速器箱底就会与事故支座相碰，扭力杆所承受的扭矩就不再增加。结论通过对扭力杆缓冲装置在最大倾动力矩作用下及事故状态下的有限元分析，可以得出以下结论：（1）上下部联接座连接螺栓根部存在

19、严重的应力集中现象，可以将此处的连接螺栓更换为双头螺柱。（2）扭力杆所受应力为脉动循环变形力，其损坏的原因为多种因素作用下的疲劳破坏。（3）扭力杆主要承受剪切力，同时也承受一定的弯矩。扭力杆光轴段可近似地视为等直圆杆纯扭转整个光轴段的应力分布为：由外表面向芯部逐渐减小，芯部应力值接近零。（4）扭力杆的危险截面为联接板作用段，最易发生疲劳破坏的部位为键槽处。（5）在最大倾动力矩下扭力杆的安全裕度还很大，其安全系数完全符合结构安全要求；而在事故状态下的安全裕度较小。（6）本文编制了可以用于分析扭力杆缓冲装置应力及位移变化的参数化程序，能大大降低试验费用，节约设计成本,为扭力杆的设计提供了一个可行的

20、方法。参考文献1谭牧田氧气转炉炼钢设备M机械工业出版社，19832小飒工作室最新经典 ANSYS 及 Workbench 教程M 电子工业出版社， 20043成大先.机械设计手册（第四版）第二卷M.化学工业出版社，2001 4宋荣生.螺纹联接轴向预紧力试验分析J.天津理工学院学报，2001:24-275张玉忠.转炉倾动机构扭力杆有限元分析 J.机械管理开发，2006:33-356Frantiek Trebua，Failure analysis of mechanical elements in steelworks equipment by methods of experimental 应工 力况扭力杆mechanicsJEngineering Failure Analysis，20097 John D. Reid, Nicholas R. Hiser, Detailed modeling of bolted joints with slippageJFinite Elements in Analysis and Design，2005(文章来源：中国冶金装备网 )