CN109117460A - 一种基于端跳测量计算转子装配轴线偏斜的方法 - Google Patents

Abstract

一种基于端跳测量计算转子装配轴线偏斜的方法，包括三接触点的计算、三角形判定准则和偏斜矩阵齐次坐标变换算法三部分；三点接触计算主要包括初始测量数据的表征、接触点的确定、偏转方向的确定及表征方法；三角形判定准则用于判定所得三点是否符合实际，详细描述了判定原则、计算方法与舍弃点的确定过程；偏斜矩阵齐次坐标变换算法包括倾角方向及大小的矩阵表征，最终装配结果整体轴线偏斜的计算方法与矩阵表示方法等。本方法基于生产实际中的实测端跳数据，实现了在装配前对轴线偏斜的预测，提高了装配后的转子同轴度，同时大大提高了一次装配通过率，对航空发动机转子件装配过程中的轴线预测、装配相位调整与优化有着重要的现实指导意义。

Description

一种基于端跳测量计算转子装配轴线偏斜的方法

技术领域

本发明属于转子装配轴线偏斜计算方法，可应用于如航空发动机高压压气机转子、高压涡轮盘、低压转子组件等重要零组件装配过程中的轴线预测、装配相位优化、装配指导过程之中。

背景技术

在装备制造过程中，装配作为一个十分重要的环节，对产品的性能及可靠性有着直接影响。而对于航空发动机转子组件而言，保证其装配后的同轴度满足要求非常重要，为了降低由于“试错”装配产生的时间和人力成本，提高装配后的同轴度，装配相位的预测及优化至关重要，而要提高一次装配通过率、降低反复拆装中产生的零件磨损及时间成本，就需要结合工厂实际装配流程，寻求各级盘在装配前端面的形貌与装配后所产生的轴线偏斜之间的关系。

转子装配轴线偏斜的测定是应用于装配相位预测优化中的重要步骤。拟合平面法是一种基于所测转子端面跳动数据进行最小二乘拟合，通过拟合平面方程计算装配堆叠偏斜的方法，这种处理算法程序简单，能够从整体上把握两结合面的形貌特征，但对于具有“双高点”或不规则形貌的转子结合面，其拟合平面与源数据间误差较大，所预测轴线偏斜的精度也难以保证；实际回转中心线法是Axiam公司所采用的轴线偏斜预测方法，其本质是通过对各级盘端面形貌跳动数据的测量，探求其装配后产生的实际回转轴线，并基于此轴线进行装配相位的调整与优化，此方法主要应用于国外，关键技术尚未引进国内，具有较高的预测准确度，但受我国高端零部件制造水平限制，零件超差现象普遍存在，本方法应用过程中没有圆度仪基本的调心调倾过程，较难适应于我国航空发动机转子的生产装配过程中。

本发明是基于对转子装配前的法兰盘结合面端面跳动的测量，不考虑法向弹性变形，提出了一种装配后轴线偏斜的计算方法，该方法可以实现对某一装配相位下各级盘的轴线偏斜预测，而后可以进一步实现装配相位的优化调整，具有重要的现实意义。

发明内容

为满足航空发动机转子件装配过程中的同轴度要求，本发明结合工程实践，基于转子件各单级盘实测跳动数据，提出了一种计算各级间装配偏斜的计算方法。

本方法的计算原理如下：

在刚性假设下，考虑两级盘装配中某一确定的安装相位，计算端面跳动与偏斜量的关系。由于不共线的三点确定一个平面，本方法的核心在于找到刚性条件下两级盘的三个初始接触点，通过这三点确定上一级盘相对下一级盘的偏斜量。需要指出的是，由于机械加工过程使得零件的表面形貌具有“宏观连续性”的特征——即虽然相邻测点测得的跳动值是离散的，但在小范围内的各测点波动较小，在较大尺度上来看是连续的特征——因而所找出的三个点可以表征刚性接触下的实际初始接触点的位置。

本方法需首先采用圆度仪对各级盘法兰结合面的端面跳动进行测量，而后以每相邻两级盘的端跳数据为输入参数，通过本算法的计算，得出相邻两级盘的相对偏斜矩阵；将每相邻两级盘的偏斜矩阵叠乘，即可得到反映装配体同轴度的总偏斜矩阵；需要指出的是调整不同的装配相位会得到不同的偏斜矩阵，但相邻两级盘的装配相位受螺栓孔、定位孔等的限制是有限的个数。以运用本方法计算出的总偏斜矩阵作为评定参数，可在装配前寻求最优装配相位，大大降低反复拆装次数，提高一次装配合格率。

附图说明

图1是端面跳动对装配偏斜影响示意图。

图1中：转子A的上、下端面分别为A_O2、A_O1；转子B上下端面分别为B_O2、B_O1；装配过程中由于两接触面A_O2、B_O1表面不平所引起的转子B相对转子A的轴线偏斜如图所示；

图2是齐次坐标变换矩阵示意图。

具体实施方式

以下结合附图和技术方案，进一步说明本方法的具体计算方式。

实施例

一种基于端跳测量计算转子装配轴线偏斜的方法，步骤如下：

步骤A：相邻两级盘，转子A、转子B的上下面中心分别为A_O2、A_O1，转子B上下面中心分别为：B_O2、B_O1，两转子配合面为转子A上端面A₂及转子B 下端面B₁。两接触面A₂、B₁用矩阵表征，数据的形式是一个圆环，分别为A(α,z)及 B(α,z)，用极坐标表示法表征α处某点的跳动值z；圆心O在全局坐标系中的位置为O(0,0),已知转子半径R；以转子A上端面A₂为基面，找到转子B下端面B₁在与A₂三点接触时的三点，此三点可确定接触后的平面；

步骤B：计算第一接触点：下端转子A固定不动，上端转子B逐渐向下平移逼近，即转子B的下端面B₁平移逼近转子A的上端面A₂，假设平移距离为d后产生第一个接触点c₁，则c₁即为两端面实际距离最近的两点。此处需指出c₁为最先接触的两点，未必是最终的三个接触点之一；

此阶段为平移接触过程，输入为上端面A₂及下端面B₁的端跳数据，A(α,z)及 B(α,z′)。找出相距最近的两点，即第一个接触点c₁，在此对两组端跳数据求和，获得求和后的数据矩阵SUM(α,z‘’)，找出数据矩阵SUM(α,z‘’)端跳和最大值对应的点即为第一个接触点，即z‘’最大值z‘’_max对应的α，即为第一个接触点c₁(α)；

计算式：z‘’＝z+z′； (1)

步骤C：计算第二接触点：使转子B的下端面B₁绕点c₁旋转继续与转子A的上端面A₂靠近，旋转方向为c₁与下端面B₁圆心O连线方向，旋转一定角度θ₁后，产生第二个接触点c₂；

此阶段为单点旋转接触过程，输入为：最大距离z‘’_max、第一个接触点c₁(α)、上端面A₂及下端面B₁的端跳数据，A(α,z)及B(α,z′)。根据第一阶段得出的接触点c₁(α)，可确定旋转方向为接触点c₁(α)与下端面B₁圆心O′的连线方向，旋转朝向圆心O′。根据第一阶段得出的最大距离z‘’_max及上端面A₂及下端面B₁的端跳数据，A(α,z)及B(α,z′)可以计算出每一组点的接触剩余距离在旋转方向上投影所对应的角度θ_i，找出θ_i最小值所对应的点，此点即为第二个接触点c₂(α₂)；

计算式：d＝z‘’_max-z-z′； (2)

步骤D：计算第三接触点：下端面B₁绕点c₁与点c₂连线的垂直平分线方向，朝向圆心O方向旋转，继续向上端面A₂靠近，旋转一定角度θ₂后，产生第三个接触点c₃；

此阶段为连线旋转接触过程，输入为：θ_i最小值第二个接触点c₂(α₂)、最大距离z‘’_max、第一个接触点c₁(α)、上端面A₂及下端面B₁端跳数据 A(α,z)、B(α,z′)。根据第二阶段得出的接触点c₁(α)、c₂(α₂)，可确定旋转方向为接触点c₁(α)与c₂(α₂)的连线垂直平分线方向，旋转朝向圆心O"。根据第一阶段得出的最大距离z‘’_max、上端面A₂及下端面B₁的端跳数据，A(α,z)及B(α,z′)、第二阶段得出的θ_i最小值第二个接触点c₂(α₂)，可以计算出每一组点的接触剩余距离在旋转方向上投影所对应的角度θ_i2，找出θ_i2最小值所对应的点，此点即为第三个接触点c₃(α₃)；

计算式：d₂＝z‘’_max-z-z′-d×l_i/l′； (4)

步骤E：三角形判定：根据获得的三点c₁、c₂、c₃判定接触三点是否符合三角判定准则。若符合，三角判定程序结束，获得三接触点。若不符合，则进入下一步骤；

三角形判定准则说明：可根据局部坐标系中的三点所构成的三角形是否包含圆心来判定三点的有效性，若圆心在三角形内，则此三角形必为锐角三角形；三角形判定准则的意义在于避免由于测点在小区域上的集中偏高(或偏低)所导致的三个接触点非常靠近的情况，这种情况不符合刚性假设下的接触稳定性，在实际接触中会使得三点所构劣弧中间的点(同时也是钝角所对应的点)脱离接触，不符合实际情况，同时也说明所选的装配相位不合适；

计算方法：可根据c₁、c₂、c₃三点连线三角形内角来判断：若三个内角均为锐角，则O"在三角形内，三点满足实际情况，三个接触点坐标即可确定；否则， O"不在三角形内，需要舍弃三点所构成的劣弧中间的点，重新寻找接触点，执行步骤F；

步骤F：此阶段为重新寻找第三接触点阶段，此阶段输入量为：第一个接触点c₁(α)、第二个接触点c₂(α₂)、第三个接触点c₃(α₃)。基于实际情况出发，当圆心O"不在三点c₁、c₂、c₃所构三角形内时，上端件B无法稳定，会继续倾斜寻找另一个点，同时三个接触点中会有一个点脱离接触，即三点中间的点脱离接触，可确定旋转方向为另外两点的连线垂直平分线方向，旋转朝向圆心O"，重新执行步骤D，以未脱离接触两点作为新的c₁(α)、c₂(α₂)接触点进行计算，而后执行步骤E进行三角形判定，直至满足三角形判定准则得出最终接触点。

步骤G：三个接触点确定后，轴线偏斜即为步骤C及步骤D计算第二、第三接触点时转子B的两次轴线倾斜的矢量叠乘；另外若在执行步骤E三角判定时反复执行了步骤F，则在计算整体轴线偏斜时需将每次执行步骤E中的偏斜矩阵都代入叠乘；如图2所示，测量面的齐次坐标变换矩阵可表示为：

式中，u、v为平移量，z为截面高度，β、α相当于平面法向量(-A,-B,1)中的 A、B，从数值上看，-A分量为绕y轴转角，-B分量为绕x轴转角的负值；按式 (6)所示计算方法将步骤C、D、F中的每次轴线偏斜用坐标变换矩阵H_i表示，如若整个过程产生了n次坐标变换，则将这些坐标变换矩阵按执行步骤对应的变换顺序叠乘即可得到整体轴线偏斜变换矩阵：

整体轴线偏斜变换矩阵H可按式(6)转化为轴线偏移方向及偏斜大小两个参数，即为最终结果。

Claims (1)

1.一种基于端跳测量计算转子装配轴线偏斜的方法，其特征在于，步骤如下：

步骤A：相邻两级盘，转子A的上下面中心分别为A_O2、A_O1，转子B上下面中心分别为：B_O2、B_O1，两转子配合面为转子A上端面A₂及转子B下端面B₁；两接触面A₂、B₁用矩阵表征，数据的形式是一个圆环，分别为A(α,z)及B(α,z)，用极坐标表示法表征α处某点的跳动值z；圆心O在全局坐标系中的位置为O(0,0)，已知转子半径R；以转子A上端面A₂为基面，找到转子B下端面B₁在与A₂三点接触时的三点，此三点确定接触后的平面；

步骤B：计算第一接触点：下端转子A固定不动，上端转子B逐渐向下平移逼近，即转子B的下端面B₁平移逼近转子A的上端面A₂，假设平移距离为d后产生第一个接触点c₁，则c₁即为两端面实际距离最近的两点；c₁为最先接触的两点，未必是最终的三个接触点之一；

此阶段为平移接触过程，输入为上端面A₂及下端面B₁的端跳数据，A(α,z)及B(α,z′)；找出相距最近的两点，即第一个接触点c₁，在此对两组端跳数据求和，获得求和后的数据矩阵SUM(α,z‘’)，找出数据矩阵SUM(α,z‘’)端跳和最大值对应的点即为第一个接触点，即z‘’最大值z‘’_max对应的α，即为第一个接触点c₁(α)；

计算式：z‘’＝z+z′； (1)

步骤C：计算第二接触点：使转子B的下端面B₁绕点c₁旋转继续与转子A的上端面A₂靠近，旋转方向为c₁与下端面B₁圆心O连线方向，旋转一定角度θ₁后，产生第二个接触点c₂；

此阶段为单点旋转接触过程，输入为：最大距离z‘’_max、第一个接触点c₁(α)、上端面A₂及下端面B₁的端跳数据，A(α,z)及B(α,z′)；根据第一阶段得出的接触点c₁(α)，确定旋转方向为接触点c₁(α)与下端面B₁圆心O′的连线方向，旋转朝向圆心O′；根据第一阶段得出的最大距离z‘’_max及上端面A₂及下端面B₁的端跳数据，A(α,z)及B(α,z′)计算出每一组点的接触剩余距离在旋转方向上投影所对应的角度θ_i，找出θ_i最小值所对应的点，此点即为第二个接触点c₂(α₂)；

计算式：d＝z‘’_max-z-z′；(2)

步骤D：计算第三接触点：下端面B₁绕点c₁与点c₂连线的垂直平分线方向，朝向圆心O方向旋转，继续向上端面A₂靠近，旋转一定角度θ₂后，产生第三个接触点c₃；

此阶段为连线旋转接触过程，输入为：θ_i最小值第二个接触点c₂(α₂)、最大距离z‘’_max、第一个接触点c₁(α)、上端面A₂及下端面B₁端跳数据A(α,z)、B(α,z′)；根据第二阶段得出的接触点c₁(α)、c₂(α₂)，确定旋转方向为接触点c₁(α)与c₂(α₂)的连线垂直平分线方向，旋转朝向圆心O"；根据第一阶段得出的最大距离z‘’_max、上端面A₂及下端面B₁的端跳数据，A(α,z)及B(α,z′)、第二阶段得出的θ_i最小值第二个接触点c₂(α₂)，计算出每一组点的接触剩余距离在旋转方向上投影所对应的角度θ_i2，找出θ_i2最小值所对应的点，此点即为第三个接触点c₃(α₃)；

计算式：d₂＝z‘’_max-z-z′-d×l_i/l′； (4)

步骤E：三角形判定：根据获得的三点c₁、c₂、c₃判定接触三点是否符合三角判定准则；若符合，三角判定程序结束，获得三接触点；若不符合，则进入下一步骤；

三角形判定准则说明：根据局部坐标系中的三点所构成的三角形是否包含圆心来判定三点的有效性，若圆心在三角形内，则此三角形必为锐角三角形；

计算方法：根据c₁、c₂、c₃三点连线三角形内角来判断：若三个内角均为锐角，则O"在三角形内，三点满足实际情况，三个接触点坐标即确定；否则，O"不在三角形内，需要舍弃三点所构成的劣弧中间的点，重新寻找接触点，执行步骤F；

步骤F：此阶段为重新寻找第三接触点阶段，此阶段输入量为：第一个接触点c₁(α)、第二个接触点c₂(α₂)、第三个接触点c₃(α₃)；基于实际情况出发，当圆心O"不在三点c₁、c₂、c₃所构三角形内时，上端件B无法稳定，继续倾斜寻找另一个点，同时三个接触点中有一个点脱离接触，即三点中间的点脱离接触，确定旋转方向为另外两点的连线垂直平分线方向，旋转朝向圆心O"，重新执行步骤D，以未脱离接触两点作为新的c₁(α)、c₂(α₂)接触点进行计算，而后执行步骤E进行三角形判定，直至满足三角形判定准则得出最终接触点；

步骤G：三个接触点确定后，轴线偏斜即为步骤C及步骤D计算第二、第三接触点时转子B的两次轴线倾斜的矢量叠乘；另外若在执行步骤E三角判定时反复执行步骤F，则在计算整体轴线偏斜时需将每次执行步骤E中的偏斜矩阵都代入叠乘；测量面的齐次坐标变换矩阵表示为：

式中，u、v为平移量，z为截面高度，β、α相当于平面法向量(-A,-B,1)中的A、B，从数值上看，-A分量为绕y轴转角，-B分量为绕x轴转角的负值；按式(6)所示计算方法将步骤C、D、F中的每次轴线偏斜用坐标变换矩阵H_i表示，如若整个过程产生n次坐标变换，则将这些坐标变换矩阵按执行步骤对应的变换顺序叠乘即可得到整体轴线偏斜变换矩阵：

整体轴线偏斜变换矩阵H按式(6)转化为轴线偏移方向及偏斜大小两个参数，即为最终结果。