基于有限元分析的薄壁件铣削加工变形仿真.docxVIP

基于有限元分析的薄壁件铣削加工变形仿真 在现代航空工业中，广泛使用薄边框架、梁、壁等框架作为主要支撑结构。这些航空薄边框架的结构复杂，相对坚固，尺寸精度要求高。在数值加工中，由于剩余力、切削力、切屑热和安装力的影响，大多数弯曲变形会导致加工变形，甚至被遗弃。目前，国内外对薄边框架加工变形的研究很多。主要基于有限模拟技术在切割过程中对薄边框架侧壁、腹部螺钉和整体变形进行工艺模拟，并预测其加工变形规律。许多研究人员认为，由于切削力的弹性变形是导致薄边框架侧壁和腹部加工变形的主要因素，而薄边框架的整体加工变形主要受剩余力的影响。此外，黄志刚等人还综合考虑了剩余力、锁削力和锁瑞龙热的因素，并研究了锁加工序列对框架壁切削结构整体加工变形的影响。在这项工作中，我们建立了一个有限模型，以研究不同毛种初量的分布、剩余力和切削力的结合对加工过程中薄边框架整体变形的影响规律。 1 薄膜研磨与加工的有限模型 1.1 面体单元和六面体单元的数值计算 选择悬臂板结构的典型薄壁件作为加工模拟对象.图1(a)所示为薄壁件毛坯的几何模型,该方形毛坯经粗加工切除多余材料Ⅰ和Ⅱ,再经精加工切除材料Ⅲ后,最终得到如图1(b)所示的成形薄壁件.工件的初始毛坯网格模型由六面体单元和五面体单元混合组成(如图2).为了比较真实地建模精加工阶段刀具-工件的接触面,用五面体楔形单元剖分多余材料Ⅲ,而毛坯的其他部分采用六面体单元.五面体单元和六面体单元在刀具进给方向(+Y)和刀轴方向(+Z)的长度相等.为简化计算,网格单元沿进给方向的长度尺寸等于瞬时加工表面宽度dw,其值由下列公式计算得到: dw=R2?(R?ae)2????????????√=2Rae?(ae)2??????????√.(1) 式中:R为螺旋齿圆柱立铣刀半径,ae为精加工阶段刀具径向切削宽度.网格单元沿刀具轴向的尺寸与dw相等或相近. 考虑薄壁件实际加工过程中广泛采用真空吸盘吸附夹紧的方式,工件底面与工作台台面充分接触.因而,在有限元建模时约束工件底面所有网格节点的平动自由度.另外,在建立工件材料模型时,充分考虑加工过程中可能产生的塑性变形, 假设在整个加工过程中材料处于弹塑性变形状态. 1.2 残余应力的分层施加 残余应力是指在没有外力作用的情况下,物体内部保持平衡而存留的应力.加工前,工件内的残余应力处于自平衡状态,即工件内任一截面内的合力与合力矩满足 ∫σrdA=0, (2) ∫dM=0. (3) 式中:σr为残余应力,M为力矩,A为截面面积.在零件铣削加工过程中,随着多余材料的不断切除,初始残余应力场的平衡状态不断被打破,工件只有通过变形才能达到新的平衡状态,因而在研究薄壁件加工变形时,应考虑残余应力因素的影响,使所建立的分析模型更准确. 要研究残余应力引起加工变形的规律,必须首先测试出零件加工前毛坯初始残余应力的分布情况,进而作为加工仿真的初始条件施加到有限元模型中.目前,残余应力的测试方法很多,但均存在局限性,特别是在测量材料内部三维残余应力场分布的时候.如果毛坯几何形体复杂,残余应力分布没有规律,即便能够准确测量毛坯的初始残余应力场,也很难将其准确施加到有限元模型中.文献视铝合金预拉伸板7050T7451内部残余应力沿厚度方向均匀分布,因而在建模时,将测量得到的残余应力分层施加到有限元模型中.分层施加法适用于残余应力分布规则、几何形体简单的模型,且施加精度与毛坯厚度方向网格划分疏密相关.本文采用有限元法模拟零件加工前材料的淬火热处理及预拉伸工艺,计算得到零件内的初始残余应力场分布. 对于有限元数值模拟淬火残余应力的形成过程, 可分为热和力两部分进行准耦合分析:先进行热传导分析,采用瞬态热传导方程求解淬火过程的温度场,并根据已求得的瞬态温度场来修正该瞬时的材料弹塑性性质,忽略塑性变形对瞬态温度场的影响;再读入前面分析得到的温度随时间变化的温度场,进行瞬态热应力场分析,从而实现准耦合模拟.淬火残余应力场就是工件温度冷却至室温时热应力场的最终状态.预拉伸消除残余应力模拟则是在淬火模拟结果的基础上,约束工件垂直于拉伸方向的一个端面在该方向上的自由度,并给予另一个端面沿拉伸方向一个位移值(预拉伸量),同时保持淬火模拟过程中的三点约束.经过上述模拟,可以得到工件毛坯淬火以及淬火后对应不同预拉伸量的残余应力场分布. 1.3 不同旋转角度下的切削力试验 采用刚性力学模型进行切削力预报,尽管柔性力学模型预测精度更好,但不同切削力模型对同一切削过程的切削力和表面误差预测结果具有很好的一致性,因而刚性力学模型亦能很好地预测薄壁件加工变形规律.如图3所示,建立刀具坐标系oxyz,z轴正向为刀轴矢量方向,x轴正向为刀具进给方向.从微观角度出发,将螺旋切削刃划分为轴向等距的微段切削刃(刀齿微元