一种榫接结构激光冲击强化的多重精度优化方法技术

本发明专利技术涉及一种榫接结构激光冲击强化的多重精度优化方法，(1)采用循环对称技术建立榫接结构模型；(2)对榫接结构模型进行二维模型的抽取和简化以提高计算效率；(3)对榫接结构二维简化模型进行激光冲击强化处理和载荷的施加，并通过有限元分析结果确定危险区域；(4)分析时间对计算精度的影响，确定多重精度优化的三种精度模型；(5)选取脉冲宽度、搭接率、光斑尺寸和峰值压力为设计参数，最大残余压应力和残余拉应力覆盖率为优化变量，结合三种精度模型，对榫接结构激光冲击强化的工艺参数进行多重精度优化，获得针对榫接结构的激光冲击强化最佳工艺参数。

【技术实现步骤摘要】
一种榫接结构激光冲击强化的多重精度优化方法
本专利技术是一种面向榫接结构激光冲击强化的多重精度优化方法，能够有效提高榫接构激光冲击强化的优化效率，获得优化后的工艺参数，属于航空航天发动机

技术介绍
航空发动机是一种极限产品，在高温、高压、高转速等的复杂载荷/环境下工作。航空发动机的涡轮盘和涡轮叶片通过榫接结构进行连接和固定，榫接结构不仅要承受涡轮本身的离心载荷和热载荷，还要承受叶片传来的离心载荷、气动载荷、热应力等。在如此苛刻的工作环境下，涡轮榫接部位极易产生疲劳裂纹，造成重大事故。大量研究表明金属构件的疲劳破坏问题与其表面完整性关系较大，因为在交变载荷的作用下构件表面产生裂纹并逐渐扩展，会导致结构件整体的破坏。对榫接部位进行激光冲击强化处理可以诱导产生残余压应力场，在一定程度上抑制疲劳裂纹萌生和扩展，从而延长构件的疲劳寿命。为了选取合适的激光冲击强化工艺参数值，对其进行优化是有必要的。激光冲击强化是高度非线性的处理过程，利用ABAQUS/Explicit对激光冲击强化过程进行显式动力学分析会有较大的时间消耗，尤其是结构的尺寸较大或者结构几何较为复杂时。即便是通过试验设计的方式来构造响应面方程，由于参数较多，构造全局响应面方程所需的试验设计点较多，计算的时间成本同样难以接受。多重精度优化可以充分地发挥各种精度分析模型的特点，在确保得到合理优化结果的前提下实现计算成本的最小化。现有文献“蔡兰,张永康.激光冲击抗金属疲劳断裂的激光参数优化试验研究[J].中国激光,1996(12):1117-1120”通过试验的方法，以合金的疲劳寿命为优化目标，脉冲宽度和能量为研究对象，对激光冲击强化技术进行参数优化，并得出了当脉冲宽度为30ns、能量为16J时，激光冲击强化对材料具有非常好的强化效果。但由于试验次数的限制和试验结果分散性的影响，试验的方法很难获得比较好的优化结果。现有文献“CorreaC,LaraLRD,DíazM,etal.Effectofadvancingdirectiononfatiguelifeof316Lstainlesssteelspecimenstreatedbydouble-sidedlasershockpeening[J].InternationalJournalofFatigue,2015,79:1-9”通过有限元模拟对激光冲击强化的工艺参数进行优化，并结合试验发现Ti-6Al-2Sn-4Zr-2MoCT试件的弯曲疲劳寿命得到了显著改善，验证了激光冲击强化的有效性。但研究对象为简单的平板试样，仍缺乏对复杂构件的激光冲击强化优化工作，且采用的三维建模计算时耗时巨大。
技术实现思路
本专利技术技术解决方案：克服现有技术的不足，提供一种面向榫接结构激光冲击强化的多重精度优化方法，能够给出针对榫接结构的激光冲击强化最佳工艺参数。工艺参数优化后榫接结构的激光冲击强化不仅使材料表面的应力大幅度降低，而且使得应力分布更加均匀。本专利技术技术解决方案：一种面向榫接结构激光冲击强化的多重精度优化方法，建立了榫接结构模型，对涡轮榫接结构复杂构型进行二维模型抽取和简化并引入激光冲击强化工艺，利用多重精度优化的方法实现针对榫接结构的激光冲击强化最佳工艺参数的获取。实现步骤如下：(1)采用循环对称技术建立榫接结构模型：模拟的某涡轮盘总共有68个榫槽，考虑到涡轮盘的榫槽是沿周向均匀分布的，且几何体、载荷均循环对称，为了简化模型，采取循环对称的技术建立模型，且不考虑振动应力、热应力、腐蚀及蠕变等因素的影响。(2)对涡轮榫接结构进行二维模型抽取和简化：抽取涡轮榫接结构的一个扇区，再平面化成为榫接结构二维简化模型，只保留榫头榫齿榫槽等关键部位，其他部分以载荷或边界条件的形式施加在榫接结构二维简化模型的边界上。榫头结构设定网格尺寸为0.5-1mm，自由分网；为避免榫接结构二维简化模型的网格尺寸因突变造成较大误差，在榫齿表面，沿着接触面和过渡区截取1-2mm深度的体，作为重点考察的部位并均匀分网，网格尺寸约为0.08-0.10mm。(3)有限元计算分析；对步骤(2)建立的榫接结构二维简化模型进行激光冲击强化处理和载荷的施加。在榫接结构底部施加径向位移约束，在三个侧面(由榫槽顶部沿涡轮盘径向方向的剖面、由榫槽底部沿涡轮盘径向方向的剖面、榫槽顶部弧面)施加轴向位移约束，在需要强化的榫齿表面施加冲击载荷，利用隐式回弹实现榫接结构二维简化模型内部的应力平衡。整个榫接结构的而为简化模型共有四对榫齿结构接触，并定义榫头齿面为主面(Master)，榫槽齿面为从面(Slave)，设置榫头榫齿间摩擦系数为0.3。定义榫槽和榫齿间的接触及离心载荷，得到榫接结构激光冲击强化的数值模型，最终通过有限元分析结果确定危险区域；(4)确定多重精度优化的精度模型：针对有限元分析软件的显式动态分析时间对计算结果精度的影响，将冲击载荷的动能耗散较为彻底(耗散为冲击载荷初始动能的1％)的榫接结构二维简化模型作为高等精度模型，将能量耗散率略差(耗散为冲击载荷初始动能的2％)的榫接结构二维简化模型作为中等精度模型，响应面模型为低精度模型，充分利用三种精度分析方法实现榫接结构激光冲击强化处理过程的参数优化。根据试验设计方法预先定义中等精度数值响应数据库；分别应用高等精度模型和中等精度模型对优化循环初始点x0处的目标函数和约束值进行分析计算；求解标度函数并修正中等精度数值响应数据库；使用修正后的中等精度数值响应数据库分析构造或更新响应面近似方程；使用响应面近似方程进行优化搜索；对优化结果应用高等精度模型进行评估；满足收敛条件则计算终止，否则更新标度函数，将计算点更新为并开始重新计算。(5)获取针对榫接结构的激光冲击强化最佳工艺参数：选取脉冲宽度、搭接率、光斑尺寸和峰值压力为设计参数，以最大残余压应力和残余拉应力覆盖率为优化变量，要求残余拉应力覆盖率不能超过榫接结构过渡齿面的一半，利用三种精度模型，对榫接结构激光冲击强化的工艺参数进行多重精度优化，经过多次重复步骤(4)的迭代过程，获得步骤(3)中确定的危险区域的最大残余压应力随迭代次数的变化图，得到针对榫接结构的最佳工艺参数。优化后，结构的最大应力从310.6MPa降低到了273MPa。榫接结构经过激光冲击强化参数优化后不仅材料表面的应力大幅度降低，应力分布也更加均匀。所述步骤(1)中采用循环对称技术建立榫接结构模型的具体为：(11)涡轮盘共有多个榫接结构，每个榫接结构具有多个榫槽(本专利技术为68个榫槽)，榫槽是沿涡轮盘周向均匀分布的，几何体、载荷均循环对称，每个榫槽中各装有一个叶片；(12)选取涡轮盘的某个榫接结构，截取所述榫接结构的一半建立模型。所述步骤(2)中，榫头结构需要预先设定网络，设定网格尺寸为0.5-1mm，自由分网；为避免榫接结构二维简化模型的网格尺寸因突变造成较大误差，在榫齿表面沿着与另一榫齿表面的接触面和过渡区截取1-2mm深度的体均匀分网，网格尺寸约为0.08-0.10mm。所述步骤(3)中，有限元计算分析时，榫接结构二维简化模型的激光冲击强化处理和载荷的施加需要分两个阶段进行，第一阶段是显式动力学分析，第二个阶段是隐式回弹和加载，即先导入显式动力学分析结果至隐式分析求解器中进行隐式回弹，然本文档来自技高网...

【技术保护点】
1.一种榫接结构激光冲击强化的多重精度优化方法，其特征在于，实现步骤如下：(1)采用循环对称技术建立榫接结构模型，基于涡轮盘的榫槽沿周向均匀分布，几何体及载荷均对称的因素，模拟涡轮盘的三维涡轮榫接结构；(2)对涡轮盘的三维涡轮榫接结构进行二维模型抽取和简化，抽取涡轮盘的三维涡轮榫接结构中的一个扇区，再将所述扇区平面化为榫接结构的二维简化模型，只保留榫接结构中的榫头、榫齿和榫槽的关键部位，其他部分以载荷或边界条件的形式施加在榫接结构二维简化模型的边界上；(3)有限元计算分析；对榫接结构二维简化模型进行激光冲击强化处理和载荷的施加，榫接结构的底部施加径向位移约束，榫接结构的三个侧面，即由榫槽顶部沿涡轮盘径向方向的剖面、由榫槽底部沿涡轮盘径向方向的剖面、由榫槽顶部弧面，施加轴向位移约束，在需要强化的榫齿表面施加冲击载荷，利用隐式回弹实现榫接结构二维简化模型内部的应力平衡；整个榫接结构二维简化模型共有四对榫齿结构接触，以榫头齿面为主面(Master)，榫槽齿面为从面(Slave)，设置榫头榫槽间的摩擦系数，定义榫槽和榫齿间的接触及离心载荷，得到榫接结构激光冲击强化的数值模型，最终通过有限元分析结果确定危险区域；(4)确定多重精度优化的精度模型，所述多重精度优化是指基于同一个模型，通过不同的处理方式得到不同精度的模型应用到优化过程；针对有限元分析软件的显式动态分析时间对计算结果精度的影响，将冲击载荷的动能耗散较为彻底的榫接结构二维简化计算模型作为高等精度模型，将冲击载荷的动能耗散率略差的榫接结构二维简化计算模型作为中等精度模型，将响应面模型作为低等精度模型，利用三种精度分析方法实现榫接结构激光冲击强化处理参数多重精度优化；(5)获取针对榫接结构的激光冲击强化最佳工艺参数：选取脉冲宽度、搭接率、光斑尺寸和峰值压力为设计参数，以最大残余压应力和残余拉应力覆盖率为优化变量，利用三种精度模型，采用非线性规划模型对榫接结构激光冲击强化的工艺参数进行多重精度优化，经过多次重复步骤(4)的迭代过程，获得步骤(3)中确定的危险区域的最大残余压应力随迭代次数的变化图，得到针对榫接结构的激光冲击强化最佳工艺参数。...

【技术特征摘要】
1.一种榫接结构激光冲击强化的多重精度优化方法，其特征在于，实现步骤如下：(1)采用循环对称技术建立榫接结构模型，基于涡轮盘的榫槽沿周向均匀分布，几何体及载荷均对称的因素，模拟涡轮盘的三维涡轮榫接结构；(2)对涡轮盘的三维涡轮榫接结构进行二维模型抽取和简化，抽取涡轮盘的三维涡轮榫接结构中的一个扇区，再将所述扇区平面化为榫接结构的二维简化模型，只保留榫接结构中的榫头、榫齿和榫槽的关键部位，其他部分以载荷或边界条件的形式施加在榫接结构二维简化模型的边界上；(3)有限元计算分析；对榫接结构二维简化模型进行激光冲击强化处理和载荷的施加，榫接结构的底部施加径向位移约束，榫接结构的三个侧面，即由榫槽顶部沿涡轮盘径向方向的剖面、由榫槽底部沿涡轮盘径向方向的剖面、由榫槽顶部弧面，施加轴向位移约束，在需要强化的榫齿表面施加冲击载荷，利用隐式回弹实现榫接结构二维简化模型内部的应力平衡；整个榫接结构二维简化模型共有四对榫齿结构接触，以榫头齿面为主面(Master)，榫槽齿面为从面(Slave)，设置榫头榫槽间的摩擦系数，定义榫槽和榫齿间的接触及离心载荷，得到榫接结构激光冲击强化的数值模型，最终通过有限元分析结果确定危险区域；(4)确定多重精度优化的精度模型，所述多重精度优化是指基于同一个模型，通过不同的处理方式得到不同精度的模型应用到优化过程；针对有限元分析软件的显式动态分析时间对计算结果精度的影响，将冲击载荷的动能耗散较为彻底的榫接结构二维简化计算模型作为高等精度模型，将冲击载荷的动能耗散率略差的榫接结构二维简化计算模型作为中等精度模型，将响应面模型作为低等精度模型，利用三种精度分析方法实现榫接结构激光冲击强化处理参数多重精度优化；(5)获取针对榫接结构的激光冲击强化最佳工艺参数：选取脉冲宽度、搭接率、光斑尺寸和峰值压力为设计参数，以最大残余压应力和残余拉应力覆盖率为优化变量，利用三种精度模型，采用非线性规划模型对榫接结构激光冲击强化的工艺参数进行多重精度优化，经过多次重复步骤(4)的迭代过程，获得步骤(3)中确定的危险区域的最大残余压应力随迭代次数的变化图，得到针对榫接结构的激光冲击强化最佳工艺参数。2.根据权利要求1所述的榫接结构激光冲击强化的多重精度优化方法，其特征在于：所述步骤(1)中采用循环对称技术建立榫接结构模型的具体为：(11)涡轮盘共有多个榫接结构，每个榫接结构具有多个榫槽，榫槽是沿涡轮盘周向均匀分布的，几何体、载荷均循环对称，每个榫槽中各装有一个叶片；(12)选取涡轮盘的某个...

【专利技术属性】
技术研发人员：胡殿印，王荣桥，李金俊，高晔，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：北京,11