[发明专利]一种基于端跳测量计算转子装配轴线偏斜的方法

权利要求书

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

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

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

此阶段为连线旋转接触过程，输入为：θ_i最小值第二个接触点c₂(α₂)、z‘’最大值z‘’_max、第一个接触点c₁(α)、上端面A₂及下端面B₁端跳数据A(α,z)、B(α，z′)；根据第一阶段、第二阶段得出的接触点c₁(α)、c₂(α₂)，确定旋转方向为接触点c₁(α)与c₂(α₂)的连线垂直平分线方向，旋转朝向圆心O；根据第一阶段得出的z‘’最大值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为平移量，m为截面高度，β、δ相当于平面法向量(-A,-B,1)中的A、B，从数值上看，-A分量为绕y轴转角，-B分量为绕x轴转角的负值；按式(6)所示计算方法将步骤C、D、F中的每次轴线偏斜用坐标变换矩阵H_i表示，如若整个过程产生n次坐标变换，则将这些坐标变换矩阵按执行步骤对应的变换顺序叠乘能得到整体轴线偏斜变换矩阵：

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