一种基于惯性回转中心调控最优化及智能学习的大型高速回转装备多级零部件选配方法与流程

本发明涉及一种基于惯性回转中心调控最优化及智能学习的大型高速回转装备多级零部件选配方法，属于机械装配技术领域。

背景技术：

航空发动机是飞机的“心脏”，其性能很大程度上决定了飞机的性能。我国要建设成为航空强国，研发和制造高性能的航空发动机是关键之一。多级转静子是航空发动机的主要部件，提升多级转静子的制造和装配精度能够有效的提升整机的性能，减小航空发动机自身振动。多级转静子的优化是目前航空发动机优化设计的研究热点，单级转子或静子的加工制造质量以及多级转静子的装配质量等多方面都对航空发动机性能有较大影响，目前急需一种多级转静子多目标优化方法。目前现有的优化方法中存在的问题为：方法或模型仅实现了多级转静子装配几何量优化测量，建立的同轴度模型未考虑绕x和y轴的旋转误差，未建立多级转静子不平衡量装配模型，不能综合优化多级转静子装配几何量、相对间隙量以及不平衡量多种参数，最终导致实用性有限。

技术实现要素：

本发明为了解决现有技术中不能综合优化多级转静子装配几何量、相对间隙量以及不平衡量多种参数，最终导致实用性有限的问题，提出了一种基于惯性回转中心调控最优化及智能学习的大型高速回转装备多级零部件选配方法，以解决航空发动机转静子圆轮廓测量精度低，装配同轴度低、装配不平衡量大的问题，改善航空发动机的性能。所采取的技术方案如下：

一种基于惯性回转中心调控最优化及智能学习的大型高速回转装备多级零部件选配方法，其特征在于，所述方法包括：

第一步、建立多级转子装配同轴度的预测模型；

第二步、结合所述预测模型，通过多级转静子装配后的偏心误差，获得转子与静子的偏心误差，并根据偏心误差求得所述转子和静子的相对跳动量，进而实现多级转静子装配后相对间隙的获取；

第三步、利用定位、定向误差引起的偏心误差传递矩阵，获得装配后各级转静子定位、定向误差引起的第n级转静子不平衡量；

第四步、将单级转静子自身不平衡量和装配后各级转静子定位、定向误差引起的第n级转静子不平衡量进行矢量相加，获得各转静子单级不平衡量；然后将各转静子单级不平衡量进行矢量叠加获得多转静子初始不平衡量，实现多级转静子初始不平衡量的计算；

第五步、建立bp神经网络，实现多级转静子初始不平衡量预测后的修正。

第六步、依据所述预测模型中的多级转静子装配后同轴度、多转静子初始不平衡量、多级转静子装配后相对间隙与角向安装位置之间的关系，建立基于各级转静子角向安装位置的多级转静子装配几何量、不平衡量和相对间隙的多目标优化模型；利用遗传算法优化各批次转静子角向安装位置，实现多参数多级转静子选配。

进一步地，第一步所述预测模型的建立过程包括：

步骤一、根据误差在装配过程中传递累积原理，确定n级转静子装配后由各级转静子定位和定向误差引起的偏心误差传递矩阵；

步骤二、确定装配后第n级转静子偏心量；

步骤三、利用同轴度的iso标准定义，确定n级转静子装配后同轴度；

结合步骤一至步骤三获得的偏心误差传递矩阵、第n级转静子偏心量以及n级转静子装配后同轴度即完成多级转子装配同轴度的预测模型的建立。

进一步地，步骤一所述偏心误差传递矩阵t0-n为：

其中，pi为第i级转静子径向测量面圆心的理想位置向量；dpi为第i级转静子径向测量面圆心位置的加工误差向量；sri为第i级转静子绕z轴旋转θri角度的旋转矩阵；sr1为单位矩阵；sxi为第i级转静子基准平面绕x轴旋转θxi角度的旋转矩阵；syi为第i级转静子基准平面绕y轴旋转θyi角度的旋转矩阵；sxj-1为第j-1级转静子基准平面绕x轴旋转θxj-1角度的旋转矩阵；syj-1为第j-1级转静子基准平面绕y轴旋转θyj-1角度的旋转矩阵；srj-1为第j-1级转静子绕z轴旋转θrj-1角度的旋转矩阵。

进一步地，步骤二所述偏心量为：

其中，dx0-n为装配后第n级转静子测量面圆心在x轴方向上的累积偏移量，dy0-n为装配后第n级转静子测量面圆心在y轴方向上的累积偏移量；pi为第i级转静子径向测量面圆心的理想位置向量；dpi为第i级转静子径向测量面圆心位置的加工误差向量；sri为第i级转静子绕z轴旋转θri角度的旋转矩阵；sr1为单位矩阵；sxj-1为第j-1级转静子基准平面绕x轴旋转θxj-1角度的旋转矩阵；syj-1为第j-1级转静子基准平面绕y轴旋转θyj-1角度的旋转矩阵；srj-1为第j-1级转静子绕z轴旋转θrj-1角度的旋转矩阵。

进一步地，步骤三所述n级转静子装配后同轴度的表达式为：

其中，dx0-i为n级转静子装配后第i级转静子测量面圆心在x轴方向上的累积偏心误差；dy0-i为n级转静子装配后第i级转静子测量面圆心在y轴方向上的累积偏心误差。

进一步地，第三步所述第n级转静子不平衡量为：

其中，ux0-n为装配后第n级转静子测量面在x轴方向上的不平衡量；uy0-n为装配后第n级转静子测量面在y轴方向上的不平衡量；m0-n为装配后第n级转静子的质量；pi为第i级转静子径向测量面圆心的理想位置向量；dpi为第i级转静子径向测量面圆心位置的加工误差向量；sri为第i级转静子绕z轴旋转θri角度的旋转矩阵；sr1为单位矩阵；sxj-1为第j-1级转静子基准平面绕x轴旋转θxj-1角度的旋转矩阵；syj-1为第j-1级转静子基准平面绕y轴旋转θyj-1角度的旋转矩阵；srj-1为第j-1级转静子绕z轴旋转θrj-1角度的旋转矩阵。

进一步地，第五步中bp神经网络建立过程包括网络参数选择及模型训练过程。

网络参数选择包括学习因子选择、隐含层层数选择及隐含层网络层数选择。学习因子采用从0.001开始，保证收敛情况下，每次3倍的速率进行递增。利用交叉验证代价函数与测试样本代价函数曲线，保证网络不处于过拟合的情况下，逐渐增加网络层数及隐含层节点数。

模型训练过程包括参数初始化、数据导入、网络正向计算、代价函数计算及权值、阈值更新五个步骤。

本发明有益效果：

本发明提出了一种基于惯性回转中心调控最优化及智能学习的大型高速回转装备多级零部件选配方法，所述大型高速回转装备多级零部件选配方法分析了航空发动机转静子的定位误差及定向误差在装配中的传递过程，确定n级转静子装配后的偏心误差的传递关系，得到多级转静子装配后同轴度预测模型；利用转静子的偏心误差可以进而得到转静子的相对跳动量，实现相对间隙的预测；同时，依据n级转静子装配后的偏心误差的传递关系，得到装配后各级转静子定位误差和定向误差引入的不平衡量，再复合单级转静子自身不平衡量，依据动平衡公式合成不平衡量，得到多级转静子装配后不平衡量预测模型；其次，考虑现场装配中影响转子初始不平衡量的各因素，基于bp神经网络建立不平衡量预测修正模型，实现转子装配初始不平衡量精确预测；最后，建立了基于各级转静子角向安装位置的多级转静子同轴度、相对间隙以及不平衡量的多目标优化模型，利用遗传算法对各批次转静子角向安装位置优化，实现了多参数多级转静子的选配。

附图说明

图1为所述大型高速回转装备多级零部件选配方法的流程图。

具体实施方式

下面结合具体实施例对本发明做进一步说明，但本发明不受实施例的限制。

实施例1：

一种基于惯性回转中心调控最优化及智能学习的大型高速回转装备多级零部件选配方法，其特征在于，所述方法包括：

第一步、建立多级转子装配同轴度的预测模型；

第二步、结合所述预测模型，通过多级转静子装配后的偏心误差，获得转子与静子的偏心误差，并根据偏心误差求得所述转子和静子的相对跳动量，进而实现多级转静子装配后相对间隙的获取；

第三步、利用定位、定向误差引起的偏心误差传递矩阵，获得装配后各级转静子定位、定向误差引起的第n级转静子不平衡量；所述第n级转静子不平衡量为：

其中，ux0-n为装配后第n级转静子测量面在x轴方向上的不平衡量；uy0-n为装配后第n级转静子测量面在y轴方向上的不平衡量；m0-n为装配后第n级转静子的质量；pi为第i级转静子径向测量面圆心的理想位置向量；dpi为第i级转静子径向测量面圆心位置的加工误差向量；sri为第i级转静子绕z轴旋转θri角度的旋转矩阵；sr1为单位矩阵；sxj-1为第j-1级转静子基准平面绕x轴旋转θxj-1角度的旋转矩阵；syj-1为第j-1级转静子基准平面绕y轴旋转θyj-1角度的旋转矩阵；srj-1为第j-1级转静子绕z轴旋转θrj-1角度的旋转矩阵。

第四步、将单级转静子自身不平衡量和装配后各级转静子定位、定向误差引起的第n级转静子不平衡量进行矢量相加，获得各转静子单级不平衡量；然后将各转静子单级不平衡量进行矢量叠加获得多转静子初始不平衡量，实现多级转静子初始不平衡量的计算；

第五步、建立bp神经网络，实现多级转静子初始不平衡量预测后的修正。

其中，bp神经网络建立过程包括网络参数选择及模型训练过程：

网络参数选择包括学习因子选择、隐含层层数选择及隐含层网络层数选择。学习因子采用从0.001开始，保证收敛情况下，每次3倍的速率进行递增。利用交叉验证代价函数与测试样本代价函数曲线，保证网络不处于过拟合的情况下，逐渐增加网络层数及隐含层节点数。

模型训练过程包括参数初始化、数据导入、网络正向计算、代价函数计算及权值、阈值更新五个步骤。

第六步、依据所述预测模型中的多级转静子装配后同轴度、多转静子初始不平衡量、多级转静子装配后相对间隙与角向安装位置之间的关系，建立基于各级转静子角向安装位置的多级转静子装配几何量、不平衡量和相对间隙的多目标优化模型；利用遗传算法优化各批次转静子角向安装位置，实现多参数多级转静子选配。

其中，第一步所述预测模型的建立过程包括：

步骤一、根据误差在装配过程中传递累积原理，确定n级转静子装配后由各级转静子定位和定向误差引起的偏心误差传递矩阵；所述偏心误差传递矩阵t0-n为：

其中，pi为第i级转静子径向测量面圆心的理想位置向量；dpi为第i级转静子径向测量面圆心位置的加工误差向量；sri为第i级转静子绕z轴旋转θri角度的旋转矩阵；sr1为单位矩阵；sxi为第i级转静子基准平面绕x轴旋转θxi角度的旋转矩阵；syi为第i级转静子基准平面绕y轴旋转θyi角度的旋转矩阵；sxj-1为第j-1级转静子基准平面绕x轴旋转θxj-1角度的旋转矩阵；syj-1为第j-1级转静子基准平面绕y轴旋转θyj-1角度的旋转矩阵；srj-1为第j-1级转静子绕z轴旋转θrj-1角度的旋转矩阵。

步骤二、确定装配后第n级转静子偏心量；所述偏心量为：

其中，dx0-n为装配后第n级转静子测量面圆心在x轴方向上的累积偏移量，dy0-n为装配后第n级转静子测量面圆心在y轴方向上的累积偏移量；pi为第i级转静子径向测量面圆心的理想位置向量；dpi为第i级转静子径向测量面圆心位置的加工误差向量；sri为第i级转静子绕z轴旋转θri角度的旋转矩阵；sr1为单位矩阵；sxj-1为第j-1级转静子基准平面绕x轴旋转θxj-1角度的旋转矩阵；syj-1为第j-1级转静子基准平面绕y轴旋转θyj-1角度的旋转矩阵；srj-1为第j-1级转静子绕z轴旋转θrj-1角度的旋转矩阵。

步骤三、利用同轴度的iso标准定义，确定n级转静子装配后同轴度；所述n级转静子装配后同轴度的表达式为：

其中，dx0-i为n级转静子装配后第i级转静子测量面圆心在x轴方向上的累积偏心误差；dy0-i为n级转静子装配后第i级转静子测量面圆心在y轴方向上的累积偏心误差。

结合步骤一至步骤三获得的偏心误差传递矩阵、第n级转静子偏心量以及n级转静子装配后同轴度即完成多级转子装配同轴度的预测模型的建立。

本发明提出了一种基于惯性回转中心调控最优化及智能学习的大型高速回转装备多级零部件选配方法，所述大型高速回转装备多级零部件选配方法分析了航空发动机转静子的定位误差及定向误差在装配中的传递过程，确定n级转静子装配后的偏心误差的传递关系，得到多级转静子装配后同轴度预测模型；利用转静子的偏心误差可以进而得到转静子的相对跳动量，实现相对间隙的预测；同时，依据n级转静子装配后的偏心误差的传递关系，得到装配后各级转静子定位误差和定向误差引入的不平衡量，再复合单级转静子自身不平衡量，依据动平衡公式合成不平衡量，得到多级转静子装配后不平衡量预测模型；其次，考虑现场装配中影响转子初始不平衡量的各因素，基于bp神经网络建立不平衡量预测修正模型，实现转子装配初始不平衡量精确预测；最后，建立了基于各级转静子角向安装位置的多级转静子同轴度、相对间隙以及不平衡量的多目标优化模型，利用遗传算法对各批次转静子角向安装位置优化，实现了多参数多级转静子的选配。

虽然本发明已以较佳的实施例公开如上，但其并非用以限定本发明，任何熟悉此技术的人，在不脱离本发明的精神和范围内，都可以做各种改动和修饰，因此本发明的保护范围应该以权利要求书所界定的为准。