基于残余应力重分布的导轨梁加工变形仿真与变形控制工艺方案研究.pdf

研究论文RESEARCH引文格式：刘雪丽，王宇航，李超然，等.基于残余应力重分布的导轨梁加工变形仿真与变形控制工艺方案研究.航空制造技术，2023,66(13):102-112.LIU Xueli,WANG Yuhang,LI Chaoran,et al.Research on machining deformation simulation and deformation control technologyscheme of guideway beam based on residual stress redistributionJj.Aeronautical Manufacturing Technology,2023,66(13):102-112.基于残余应力重分布的导轨梁加工变形仿真与变形控制工艺方案研究刘雪丽，王宇航”,李超然，李彩红，张晟玮，孙剑飞2.3（1.中国航发航空科技股份有限公司，成都6 10 50 3；2.北京航空航天大学机械工程及自动化学院，北京10 0 19 1；3.北京航空航天大学江西研究院景德镇分院，景德镇3330 0 0）【摘要为解决大型薄壁铸件导轨梁在材料去除过程中因残余应力的释放与重分布导致变形超差的工艺难题，对导轨梁零件加工工艺进行分析，实现零件模型的简化与子结构分割；同时开展零件毛坏表面残余应力测量，成功建立零件毛坏初始应力模型。在此基础上结合实际加工工艺开展零件加工有限元仿真，模拟加工过程中由于材料去除引起的残余应力释放，预测了加工过程中残余应力重分布规律和加工变形情况。总结了零件加工变形的有限元仿真结果，提出抑制零件加工变形的工艺方案。经验证，改进后的工艺顺序使零件最大变形量由0.48 5mm降至0.0 8 1mm，降低8 3.3%，避免了零件在加工过程中的尺寸超差。同时该平面作为后续加工的基准，保证了后续加工的精度，为生产工艺优化提供了有效的理论依据。关键词：导轨梁；残余应力；有限元仿真；加工变形预测；工艺优化Research on Machining Deformation Simulation and Deformation Control Technology Scheme ofGuideway Beam Based on Residual Stress RedistributionLIU Xuel,WANG Yuhang,LI Chaoran,LI Caihong,ZHANG Shengwei,SUN Jjianfei?3(1.AECC Aero Science and Technology Co.,Ltd.,Chengdu 610503,China;2.School of Mechanical Engineering and Automation,Beihang University,Beijing 100191,China;3.Jingdezhen Branch of Jiangxi Research Institute,Beihang University,Jingdezhen 333000,China)ABSTRACT In order to solve the process problem of over-deformation caused by the release and redistribution ofresidual stress in the material removal process of large thin-walled casting guideway beam,the processing technology ofguide beam parts was analyzed to realize the simplification of part model and substructure segmentation.At the same time,the surface residual stress measurement of the part blank was carried out,and the initial stress model of the part blank wassuccessfully established.On this basis,the finite element simulation of part machining is carried out in combination withthe actual machining process.The residual stress release caused by material removal during machining is simulated,andthe residual stress redistribution law and machining deformation during machining are predicted.Summarize the finiteelement simulation results of the machining deformation of the parts,and propose a process plan to suppress the machiningdeformation of the parts.It is verified that the improved process sequence reduces the maximum deformation of the partsfrom 0.485 mm to 0.081 mm,which is reduced by 83.3%,and avoids the oversize of the parts during processing.At thesame time,this plane is used as the benchmark for subsequent processing,which ensures the accuracy of subsequentprocessing and provides an effective theoretical basis for the optimization of production process.*基金项目：大飞机专项（JPPT-F2019-12-1）。102航空制造技术2023年第6 6 卷第13期RESEARCH研究论文Keywords:Guideway beam;Residual stress;Finite element simulation;Machining deformation prediction;Process optimizationD0I:10.16080/j.issn1671-833x.2023.13.102导轨梁是发动机外涵道反推装置的重要组成部分，其主要作用是推动滑出外涵道，将燃气流或风扇出口的空气流转折成一定角度，斜向前方45喷出，产生与飞机飞行方向相反的推力，以达到缩短飞机着陆滑行距离的目的2-3。导轨梁模型如图1所示。在飞机降落过程中，导轨梁上的吊耳等部位将会承受反推装置重力、高速飞行时的风阻力、推出后与飞机形成的振动，以及反复推出、收回的摩擦力等。导轨梁属于典型的大尺寸薄壁零件，具有刚性低、材料去除率大等不利于铣削加工的特点，有必要对其加工变形量进行预测和控制4。为了获取毛坏内部残余应力分布模型，预测材料去除产生的加工变形规律，本文将对ZL116铸造铝合金毛坏关键表面进行残余应力测量，对其材料去除阶段进行加工变形仿真，通过仿真变形结果提出合理的加工工艺改进方案，为企业实际生产提供合理有据的帮助。1零件加工工艺分析1.1研究对象分析图2 所示为导轨梁零件的A基准面（A面）B面、三角框（C面）内外侧3个主要加工表面，该零件具备几个尺寸特征：整体尺寸大，长度跨度为1450 mm；B面为垂直于A面的平面，且内部需要加工主、副导轨槽；对位图1导轨梁模型Fig.1Guideway beam model置精度要求较高；壁厚较小、加工去除余量大，精加工后零件最薄处仅为3mm，主要加工表面材料去除比率大于50%；空腔结构不对称，B面的背面布满加强筋等，铸造应力大且分布不均，在材料去除过程中易因应力释放而产生较大的加工变形5-6 。O1.2零件工艺分析本文主要研究内容为导轨梁毛坏阶段3个区域的主要加工工艺流程，如图3所示，即加工A面（误差0.1mm）加工B面（误差 0.5mm）加工C面内外侧(误差 0.5mm)。1.3加工难点分析材料的去除过程中，毛坏内部残余应力逐步释放和重分布产生加工变形，A面作为后续加工过程的基准，B面图2 导轨梁主要加工表面Fig.2Main processing surface of guideway beam粗加工A面（材料去除量1.5mm）半精加工A面（材料去除量1mm）精加工A面（材料去除量0.5mm）加工C面外侧加工C面内侧粗加工B面（材料去除量1.5mm）半精加工B面（材料去除量1mm）精加工B面（材料去除量0.5mm）图3零件加工工艺流程Fig.3Process flow of parts processing2023年第6 6 卷第13期航空制造技术10 3三角框内外侧A基准面研究论文RESEARCH其产生的加工变形必然会影响到下一加工部位的加工精度,因此其在整个导轨梁的加工过程中最为关键；在加工A面的同时去除了部分C面外侧的材料，随后进行C面外、内侧的材料去除，经现场调研，C面外侧加工后超差率较高，需要同时关注这两个部位的加工变形情况；B面作为导轨梁中最大的平面，跨度极大,其平面度与A面的垂直度难以保证。1.4零件简化与子结构分割在不影响残余应力测量和材料去除仿真结果的前提下，为了方便构建仿真模型，避免仿真过程中的不收敛情况，在原始模型中删除后续加工步骤中生成的特征,如孔、圆角、倒角、密封槽等，同时也删除一些对全局分析无明显影响的无规则特征。零件模型简化前后的对比如图4 所示。由于零件毛坏尺寸跨度较大，且结构复杂，进一步简化研究对象,图5（a）中圈注的部位不属于上述研究对象，本身不涉及加工，其内部残余应力并不会引起加工B面时产生的加工变形,因此这一部分不作为后续研究对象，在三维模型中将此部位拆分出去，为了保证试验中的毛坏与简化模型一致，用线切割将此部分切下，简化后的模型如图5（b）所示。由于第1步加工的A面和后续加工的C面位于B面的一端,虽然在工艺基准和尺寸要求上相互关联，但从加工变形的角度分析,加工A面产生的变形不会影响B面产生变形，只会对后续加工的基准产生影响，同时，加工B面产生的变形对A面产生的影响也微乎其微。而A面和C面存在明显的关联性,其结构相连,加工过程产生的变形也必然相互影响。因此,为了同时关注两个部位的加工变形情况，便于进行加工表面的残余应力测量和后续的仿真计算，在此引人子结构法，将导轨梁的C面和A面连接部分从导轨梁的整体中独立成1个子结构进行研究，三维模型中的分割位置如图5（c）所示，左侧为后续研究的子结构，为了保证试验中的毛坏与简化模型一致，利用线切割将实际毛坏沿此模型中的分割区域切开。2零零件毛坏初始残余应力测量2.1试验设备和仪器试验使用的仪器设备主要有Prism残余应力测量仪及数控加工中心VL-866。P r is m 残余应力测量仪如图6 所示。钻孔是最常用的通过应力释放来测量残余(a）无关区域（b）无关区域切除（c）子结构划分图5 子结构分割示意图Fig.5Schematic diagram of substructure division气压箱转速控制器一激光发生仪空气吹屑嘴高速气动钻30(a）简化前夹具一45(b）简化后图4 零件模型的简化Fig.4Simplified part model104航空制造技术2023年第6 6 卷第13 期激光发生仪摄像头图6 Prism残余应力测量系统Fig.6 Prism residual stress measurement system图7 残余应力测量现场RESEARCH研究论文应力的方法 7。通过在指定材料表面钻一个小盲孔，小角度并对焦。孔周围会自发地建立一个新的应力平衡，使孔附近表面（4）对刀。设定待测零件测量表面处于零点位置，发生位移与形变，通常采用贴应变片来测量表面发生的并进行基础设置和光学设置。变化。本文采用Prism电子散斑干涉钻孔残余应力测（5）测量。电脑自动控制进给与干涉图形的拍摄量仪（芬兰Stresstechgroup公司）来测量毛坏表层的残工作，根据进给设置参数，钻头每钻到设定深度，都会进余应力，该设备基于钻孔法，并结合数字成像和电子散行拍照，并留下当前深度下孔周围表面的散斑图，测量斑干涉技术（Electronic speckle pattern interferometry，现场如图7 所示，残余应力测量参数如表1所示。ESPI),相比其他残余应力测量设备，Prism采用最新的（6）测量结束分析计算数据，通过对比各个深度拍ESPI技术，单层测量时间小于5 min，测量绝对误差小照得到的散斑图，求解不同深度X、Y 方向和切向的残于7 MPa，相比于传统贴应变片的钻孔法，省去了贴应余应力分布。变片的时间和设备调试时间，而基于电子散斑的非接触为了较好地表征待测零件的残余应力，以特征为参式测量原理，排除了传统测量的诸多干扰因素，钻孔与考将测量表面分成若干个测量区域,用测点处所测得的测量全程由电脑控制，极大提高了测量效率和测量精深度方向残余应力代表该区域同一深度的平均残余应度。同时,Prism的钻孔系统可用于多种晶体、非晶体及力。由于A、C 面紧邻且有相同平面,故在下文的试验单相和多相材料。分析过程中将A、C 面放在一起讨论。2.2测量原理与区域划分2.2.2A、C 面区域划分在试验过程中，由于Prism残余应力测量系统单理论上，同一测量面上选取的测量点越多,区域划次测量深度最大为2 mm，为获得导轨梁毛坏各个加工分越密集，其平均残余应力代表整体残余应力就越精面上大于2 mm深度的残余应力，需采用剥层法，剥层确。考虑测量效率和时间,将A、C 面分为12 个测量区所使用设备为数控加工中心VL-866。由于每一次剥域,测量区域划分及测量点的选取如图8（a）所示。该层后，零件内部的残余应力都会因材料去除而重新分测量法是对某局部范围内的残余应力平均值进行测量，布，因此只有第1层测量的残余应力才是初始残余应而在不同深度范围内的残余应力为逐层测量。力。在测量余下部分的残余应力时，此时零件内部的2.2.3B面区域划分残余应力不再是初始残余应力，所测得残余应力需要将B面分成8 个测量区域,区域的划分和测量点的进行修正。该测量方法是对某局部范围内的残余应力选取如图8（b）所示。平均值进行测量，而且在不同深度的残余应力是逐层2.3测量结果与分析测量，因此在测量点的残余应力可以看作是按层分布。测量完成后，对测量数据进行求解计算，每个测点首先利用Prism激光干涉仪测量前2 mm深度方向的残余应力，然后采用局部剥层法，将小孔周围的材料去除，剥层厚度为2 mm，再利用Prism激光干涉仪从剥层后的2 mm深度处开始测量，此时测量2 4 mm深度的残余应力，如此循环操作完成0 6 mm厚的零件毛坏残余应力测量。由于采用局部剥层法，剥层区域远远小于残余应力划分区域，因此对剥层后获得的残余应力不再予以修正。2.2.1测量过程在本试验中，被测材料为ZL116铝合金，其弹性模量为7 0 0 0 0 MPa，泊松比为0.3 9 18 。Prism残余应力测量系统单次测量具体过程如下。（1）首先按照指定顺序打开各个设备，将储气罐冲压，激光发生仪预热10 min，使激光发生器的功率稳定。（2）安装待测零件，调试零件位置，将零件待测表面垂直对准钻头，并手动进给主轴，检查钻头行程是否能够完成整个盲孔深度的测量。（3）打开辅助照明设备，调整摄像头、激光发生仪Fig.7Scene of residual stress measurement表1残余应力测量参数Table 1 Parameters of residual stress measurement钻削条件测量材料钻头材料钻头直径/mm钻头转速/（rmin）进给速度/（mms）2023年第6 6 卷第13 期航空制造技术10 5参数ZL116硬质合金3.17590000.005研究论文RESEARCH的单次测量深度为2 mm，在2 mm的深度范围内,将残余应力测量层分为4 层，每0.5 mm得到一组X向和Y向的残余应力测量数据。导轨梁A面和B面的加工余量为3 mm，C 面的加工余量为9 mm，A 面和B面的毛坏状态壁厚为9 mm，C 面毛坏状态壁厚为12mm,考虑到C面的结构可以近似为内外对称结构,可以将12 mm壁厚的中性层作为对称面,其剩余6 mm未测得的残余应力层等效为已测6 mm残余应力的对称应力，这样就可以获得C面毛坏整个区域的全部应力状态 9 。对于A面的残余应力同样做对称等效。由于A面毛坏壁厚为9 mm,将其分为6 层,分别为2 mm、2mm、0.5 m m、0.5 m m、2 m m、2 m m,残余应力测量深度为4.5 mm。这样的等效处理可以大幅减少测量所花费的时间提高试验效率,且要处理的数据也相应减少。B面毛坏壁厚9 mm，如果在已测量6 mm深度的前提下进行4 6 mm深度的剥层,然后测量6 8 mm深度的残余应力,其剩余壁厚仅剩1mm，P r i s m 测量系统是基于盲孔法测量，其测量结果将不准确；同样剥层次数过多，剥层区域残余应力释放导致下一层残余应力重分布产生的累计误差将不能被忽略。综上两点原因，对于剩余3mm残余应力层未予以测量。且其为薄板结构,外侧是加工表面，内侧连接筋板，越靠近内侧越不容易产生加工变形,因此,将剩余3 mm等效为无应力状态。2.3.1A、C 面测量结果图9 所示为A、C 面各区域的残余应力测量值。其中A面各区域平均残余应力基本在-2 0 0 2 0 0 MPa浮动，其中区域5 和区域9 深度方向平均残余应力波动较大。区域5 位于C面与导轨梁主体的过渡部位,且位于A面的中央区域，而区域9 位于A面的边缘，在后续10（a）A、C面区域划分（b）B面区域划分图8 测量区域划分及测点位置选取Fig.8Division of measuring area and selection of measuring pointlocation106航空制造技术2023年第6 6 卷第13 期的加工变形仿真过程要注意与C面连接部位和边缘部位的加工变形情况；C面整体残余应力幅值较小，最大残余应力在5 0 MPa左右，且沿深度方向呈对称分布，X方向和Y方向残余应力虽然数值不同，但分布规律类似，由于C面材料去除比率达7 5%，且内外材料均需要去除，其中内外交替的去除顺序会影响残余应力的演化过程，因此在后续仿真中需要关注去除顺序不同产生的变形。2.3.2B 面测量结果图10 为B面各区域的残余应力测量值，这3 个区域深度方向应力波动较大，区域1是B面上最远离A面的部分，区域4 和5 分别位于副导轨槽的两侧，位于B面的中间区域,从A面各区域和B面的残余应力测量结果可以得出，导轨梁两个平面的中间区域应力波动较大，且应力值较大。其原因可能是较大平面上的材料逐渐去除,其中间区域和边缘区域的材料冷却速度不同，造成其应力分布不均或应力波动较大。3零件加工过程仿真及变形预测本文采用ANSYS作为加工变形仿真的分析工具。为了使仿真模型顺利收敛,使用非线性分析功能对模型的A面、B面、C面毛坏施加上述测得的初始残余应力，模拟材料去除过程，预测各个加工阶段完成后零件的加工变形情况。3.1仿真流程加工变形仿真的主要流程为：将上文中简化后并分成两部分的零件毛坏三维模型导入有限元软件；划分网格，并将上文所测得的残余应力值代人有限元模型模拟零件的残余应力并施加约束；根据不同工序创建多个载荷步，每个工序对应1个载荷步；使用生死单元法模拟零件在加工过程中的材料去除,其原理就是将每个加工11步骤中要去除的材料单元刚度矩阵乘以一个很小的因子，使其各物理参量都趋近于0，不会对计算结果产生2影响；最后得出加工变形云图。具体加工变形仿真流程9如图 11 所示。3.2仿真结果3.2.1A、C 面子结构仿真结果子结构包含导轨梁的A面和C面部分。网格划分与坐标建立如图12 所示。由于残余应力层的厚度决定了划分单元的大小，所以在划分网格时,将施加应力层部分的网格单元大小划分为0.5 mm。按照加工工艺中各加工阶段的去除余量对加工表面进行分层去除，加工顺序分为5 个阶段：（1）去除A面1.5mm余量；（2）去除A面1mm余量；（3）去除A面0.5mm余量；（4）去除C面外侧余量；（5）去除C面内侧RESEARCH研究论文150100500-50100-150-2000150+X100Y500-50-1000100500-50-100040X200-2040-600100-XY12深度/mm112深度/mm112深度/mm123456深度/mmFig.9Residual stress values in each area of plane A and C100505000-50XY-100-1503434XY34X-50-10050200-X150-Y100500-5050200150100500-5050500-50区-10 00图9 A、C 面各区域残余应力测量值11-XY12深度/mm(a)A面-XY123456深度/mm(b)C面2深度/mm23深度/mm334445550500-50-1000200150X-Y100500-50-100040X20Y0-20-4001X-Y112深度/mm112深度/mm12深度/mm2深度/mm3334563444555余量。从工艺中可以得出，A面作为后续加工的基准面，其Z方向的变形会严重影响后续加工精度，因此提取每一加工阶段后剩余A面在Z方向的变形量，如图13所示。从仿真结果可知，A面材料去除1.5 mm时，Z方向的变形主要集中在C面连接A面一侧中心区域和没有支撑连接的A面弧形边缘区域，当精加工A面完成后，加工C面时，A面基本不产生加工变形。图14 为C面材料去除后的加工变形。由于导轨梁的A面是后续加工基准面，因此对Z方向加工精度要求较高，图15 描述了A面弧形边缘部位和C面连接A面一侧随着材料去除最大变形量的变化情况。如图15 所示，A面区域8 和9 在加工阶段均比较稳定,最大变形量维持在0.1mm左右,区域6 虽然有一定波动,但总体也稳定在-0.3 -0.1 mm之间。因为C面在未加工之前,其毛坏壁厚为12 mm,本身结构较为稳定，因此C面连接A面一侧的区域10 在前4 个加工阶段时变形量较为稳定，在-0.0 5 mm左右。在第5 个加工阶段，C面连接A面一侧变形量从-0.0 4 2 mm变化到-0.4 8 5 mm,其原因可能是第4 阶段加工完成后，C面整体应力水平较高，C面的最终的保留部分和内侧待去除余量部分产生平衡，而零件整体壁厚仍较厚，最薄2023年第6 6 卷第13 期航空制造技术10 7研究论文RESEARCH300250200150100500501502002503000123456深度/mm(a)区域1300100500501001502002503000123456深度/mm(d)区域425020025050050-100-150-200250-3000123456深度/mm(b)区域2X250Y50-100150200250-3000123456深度/mm(e）区域5250200250X20015010050050100150200250300025020015010050050-100-15020025013000123456深度/mm(f)区域6XX123456深度/mm(c)区域3X150050-100-150-200250-3000123456深度/mm(g)区域7图10B面各区域残余应力测量值Fig.10Residual stress values in each area of plane B导人简化后二次建模、的三维模型编制命令流Fig.l1 Process of finite element machining deformation simulation处为9 mm，第5 阶段完成后，随着C面内侧余量去除，C面最薄壁厚仅为3 mm,结构强度迅速下降，由于材料的去除是C面连接A面一侧的背面,随着C面内部应力平衡被打破，零件朝着相反的方向产生变形，最大变形量为-0.4 8 5 mm。3.2.2B面仿真结果B面加工变形仿真流程与3.2.1节相同，图16 所示为B面网格划分和坐标建立。按照加工工艺要求，与实际加工中以A面作为基准相同,在仿真中对A面施加全约束，分3 个加工阶段对B面余量分层去除：（1）去除B面1.5 mm余量；（2）去除B面1mm余量；（3）去除B面0.5 mm余量。从工艺中可以看出，B面是导轨梁零件中最大的平面,需要在此面加工出主导轨槽和副导轨槽108航空制造技术2023年第6 6 卷第13 期050-100-X150Y20025013000123456深度/mm(h)区域8划分网格、对子结构体施加残余应力底面施加约束图11有限元加工变形仿真流程Fig.12Substructure meshing of plane A and C等特征,其对平面度和与A面的垂直度都有很高的要求，因此应主要关注各个加工阶段导轨梁在Y方向和Z方向的变形。而仿真结果显示，导轨梁在X方向几乎不产生生死单元法模拟材料去除图12 A、C面子结构网格划分多载荷步输出变形求解云图RESEARCH研究论文NODAL SOLUTIONNODAL SOLUTIONSTEP=6STEP=3SUB-1SUB=1TIME=6TIME=3UZ(AVG)UZ(AVG)RSYS=11RSYS-11DMX-0.286E-03DMX=0.116E-03SMN=-0.285E03SMN=-0.111E-03SMX=0.135E-03SMX=0.772E-04变形量/m-0.111E030.691E-04-0.900E-040.482E-04（a）粗加工后零件变形变形量/m0.273E04-0.145E-040.638E05-0.354E040.563E040.772E040.285E03-0.238E03-0.192E-03-0.145E-03(a)C面外侧去除-0.984E-040.517E04-0.508E-050.416E04-0.882E04-0.135E-03NODAL SOLUTIONNODALSOLUTIONSTEP=4SUB=1TIME=4UZ(AVG)RSYS-11DMX=0.144E-03SMN=-0.136E-03SMX-0.137E-03-0.136E030.754E04-0.106E03-0.451E04（b）半精加工后零件变形NODAL SOLUTIONSTEP=5SUB=1TIME=5UZ(AVG)RSYS=11DMX=0.287E-03SMN=-0.285E-03SMX=0.127E-03-0.285E03-0.193E-03-0.239E03-0.148E03（c）精加工后零件变形图13 A面逐步加工Z向变形Fig.13Gradually processed deformation in Z-direction of plane A变形,因此在后续分析中不予赘述。各加工阶段Y方向和Z方向的整体变形情况如图17 所示。从图17 可以看出，在加工导轨梁B面时，导轨梁远离A面一端一直存在向Y+方向偏移的变形,粗加工B面完成后，变形量为0.2 7 8 mm；半精加工完成后，变形量为0.3 9 1mm，精加工完成后最终稳定在了0.3 8 8 mm，说明导轨梁的Y方向加工变形主要产生于粗加工阶段，其余量去除比率超过5 0%以后，其远离A面一端Y方向的变形相对稳定。而对于导轨梁Z方向上的变形,各个阶段的变形量均小于0.1mm,整体变形较小,其原因可能是相对于导轨梁在Z方向的尺寸B面去除的余量STEP=7SUB=1TIME-7UZ(AVG)RSYS-11DMX=0.485E-03SMN=0.485E-03SMX=0.144E-03变形量/m-0.485E03-0.147E-040.460E04-0.156E04-0.763E04变形量/m-0.102E03-0.101E-04-0.559E-04-0.357E04变形量/m0.345E-03-0.205E030.107E030.415E030.137E03Fig.14Deformation of plane C after removal of internal and0.20.10-0.1-0.2-0.815E04-0.30.127E03-0.4-0.5Fig.15Z-direction deformation of key area of plane A and CY图16 B面网格划分与坐标建立Fig.16Mesh generation and coordinate establishment of plane B尺寸较小，B面的下方是筋板等结构，能够增加导轨梁在Z方向的强度来抵抗变形,虽然在加工过程中产生了Z方向的微小变形，但能通过加工余量来予以修正。2023年第6 6 卷第13 期航空制造技术10 90.653E040.275E03-0.135E03（b）C面内侧去除图14 C面内外侧材料去除后的变形externalmaterials区域2+区域6区域8、9+区域10图15 A、C 面重点区域Z向变形量-0.745E04-0.416E-05-0.144E03研究论文RESEARCHNODAL SOLUTIONSTEP=3SUB=1DMX=0.327E-03NODAL SOLUTIONTIME=3UY(AVG)SMN=-0.475E-05SMX=0.278E-03RSYS=11STEP=3SUB-1DMX=0.327E-03TIME=3SMN=-0.172E-03UZ(AVG)SMX=0.139E-04RSYS=11变形量/m-0.475E-05-0.582E-04-0.267E04-0.896E-04（a）粗加工Y方向变形NODAL SOLUTIONSTEP=4SUB=1DMX=0.393E-03变形量/m-0.121E03-0.184E-03-0.153E03-0.215E03TIME=4UY(AVG)SMN=-0.255E-05SMX=0.391E-03-0.247E030.278E03RSYS=11-0.172E-030.152E03NODALSOLUTIONSTEP=4SUB=1DMX-0.393E-03-0.131E03-0.110E03（b）粗加工Z方向变形TIME=4UZ(AVG)SMN=-0.163E-04-0.896E04-0.689E-04-0.482E-04-0.275E04RSYS=11SMX=0.405E-04-0.676E05-0.139E-04变形量/m-0.255E-05-0.849E04-0.412E04-0.129E03（c）半精加工Y方向变形NODAL SOLUTIONSTEP-5SUB=1DMX=0.388E-03变形量/m0.172E03-0.260E-03-0.216E03-0.304E-03TIME=5UY(AVG)SMN=-0.301E-05SMX=0.388E-030.347E030.391E-03RSYS=110.163E04-0.994E-05NODALSOLUTIONSTEP-5SUB-1DMX-0.388E-03-0.363E05-0.268E-05（d）半精加工Z方向变形TIME-5UZ(AVG)SMN=-0.320E-04-0.899E050.153E-040.216E04-0.279E04RSYS=11SMX=0.645E-050.342E-040.405E-04变形量/m-0.301E-05-0.838E04-0.404E-04-0.127E03（e）精加工Y方向变形Fig.17 Gradually processed deformation of plane B变形量/m-0.171E03-0.257E030.214E030.301E03图17 B面逐步加工变形-0.344E03-0.388E03-0.320E04-0.277E-04-0.234E040.192E04（f）精加工Z方向变形-0.149E04-0.106E-04-0.636E050.209E05-0.218E05-0.645E054零件加工变形控制工艺方案零件加工变形控制工艺方案主要针对区域加工顺序、切削参数、工艺路线等3 个方面 10 。由于导轨梁的3个加工区域A、B、C 面都属于较为平整的平面,故工艺路线较为简单,故延续粗加工半精加工精加工的工艺路线，主要针对前文划分的残余应力测量区域进行加工顺序的调整。4.1A、C面工艺优化方案A面以及与C面相连面为后续加工装夹的基准,如图18 所示。因此对该面Z方向的加工变形要求较高。由3.2 节的仿真结果可知，原始工艺下，C面的区域10110航空制造技术2023年第6 6 卷第13 期图18 导轨梁装夹方案Fig.18Clamping scheme of guideway beam研究论文RESEARCH在最后一道工序完成后产生较大变形,且变形方向朝向基于这一规律制定了导轨梁子结构加工改进工艺方案，已加工表面，若在导轨梁的后续加工中继续以A面作为该方案的主要思路是在A面精加工之前,提前完成C面工序基准,则会严重影响导轨梁的加工精度,因此需要对的内外侧材料去除，使区域10 的变形得以释放，而由于现有工艺进行优化。根据毛坏初始残余应力获取试验得变形方向正好朝向A面一侧，这样就可以在精加工A面到的应力数据，基于去应力的加工思想，让毛坏初始应力时，修正区域10 产生的翘曲变形，使其变形量达到最小，尽早释放,以便在加工过程中对其进行逐阶段修正,达到改进后的工艺顺序为：（1）去除A面1.5 mm余量；（2）使实际加工变形量趋向于最小的目的。由于区域10 会去除A面1mm余量；（3）去除C面外侧余量；（4）去除在所有工序完成时产生乙负方向的变形，而理论上，在无C面内侧余量；（5）去除A面0.5 mm余量。A、C面Z向约束条件下让零件自由释放，无论零件各部分材料的去加工变形仿真结果如图19 所示。除顺序如何,最终零件产生的加工变形情况是近似的,计算可得优化后的工艺在C面内侧材料去除及精NODAL SOLUTIONSTEP=3SUB=1TIME=3UZ(AVG)RSYS=11DMX-0.112E-03SMN=-0.106E-03SMX-0.769E-04-0.106E-03-0.657E-04-0.861E04-0.453E-04(a)A面粗加工NODAL SOLUTIONSTEP=4SUB=1TIME=4UZ(AVG)RSYS=-11DMX-0.131E-03SMN=-0.109E-03SMX=0.124E-03变形量/m-0.249E-040.158E-04-0.456E050.362E-04变形量/m0.566E-04-0.109E-030.769E-04-0.827E04-0.569E-04-0.311E04(b）A 面半精加工-0.530E050.205E-040.463E040.721E-040.980E-040.124E-03NODAL SOLUTIONSTEP=5SUB=1TIME=5UZ(AVG)RSYS=11DMX=0.162E-03SMN=-0.161E-03SMX=0.130E-030.612