转子动平衡教程PPT课件.ppt

2020 4 9 1 第三章转子动平衡 低速动平衡 刚性转子动平衡 工艺平衡装配平衡一步平衡多步平衡本机平衡整机 台架 平衡国际标准 ISO1940高速动平衡 柔性转子动平衡 模态平衡法影响系数法混合法参考标准 DIS5406 柔性转子动平衡 标准草案DIS5343 评价柔性转子平衡的准则 参考 2020 4 9 2 第一节概述 一 刚性转子动平衡1 静平衡目标 平衡静力方法 随遇平衡法2 动平衡目标 平衡力与力矩方法 动平衡机 低速平衡工艺平衡装配平衡一步平衡多步平衡本机平衡整机平衡国际标准 ISO19403 动平衡基本要求1 至少两个平衡面 平衡面距离要远 并尽能靠近支点 平衡配重半径位置要尽可能大 以便达到最大效果 2020 4 9 3 2 一步平衡 多为短寿命军用发动机采用3 多步平衡 多为长寿命民用发动机采用4 平衡方法 寻找重点寻找轻点 频闪法 影响系数法极坐标矢量图法三元平衡法5 原理 不平衡力Pj产生支反力FP1与FP2不平衡力矩RL产生支反力FR1与FR2则在支点有合力动平衡 动平衡精度1 me G0 g cm 工程实际应用2 e G0 mm s 国际标准 ISO1940将平衡品质分为11个等级 按比值为2 5的等比级数递增排列 2020 4 9 4 第二节柔性转子动平衡 一 柔性转子平衡特点1 柔性转子 n ncr1 转轴产生弯曲变形2 高速动平衡 多平面 多转速平衡过程目的 1 将不平衡力与不平衡力偶降到许可范围2 将n阶固有振型不平衡量降到许可范围3 标准 1 国际标准草案DIS5406 柔性转子动平衡 2 参考标准5343 柔性转子动平衡 4 方法 1 振型 模态 平衡法2 影响系数法3 混合法等 2020 4 9 5 5 平衡特点1 刚性转子 低速平衡后 在工作转速以下运行平稳 2 柔性转子 低速平衡后 仅平衡了低速下支承动反力 高速下轴产生弯曲变形 弯矩将随转速发生变化 支承动反力也将发生变化 3 柔性转子动平衡目的 在工作转速下 尽可能消除支承动反力 并使转子沿轴长的弯矩最小如图3 1所示 刚性转子有对柔性转子有 2020 4 9 6 F为转子变形产生的离心力 4 影响因素多 a 不同转速下挠度影响b 各阶振型对平衡的影响5 实际发动机只有少数几个平面可用于平衡 只能在有限个转速上得到平衡 6 问题 如何利用少数几个平面来获得一定转速范围内转子的良好平衡 7 假设条件 a 在一定平衡条件下 轴承振幅与转子不平衡量成正比 b 轴承振幅与不平衡力之间的相位不变 c 转子中非线性因素 如油膜 等影响 不影响上述假设条件 2020 4 9 7 二 转子在不平衡力作用下的运动方程设一转子为等截面轴 简支在各向同性的支承上轴的面积为A 单位长度质量为 A 截面质心为G z 截面偏心距为 z 质心连线为一空间曲线 如图所示 根据牛顿运动定律 得到yoz平面内的运动方程 其中则有由材料力学可知代入运动方程得到 2020 4 9 8 同理可得到xoz平面内的运动方程为引入复数表达式 令则有式中 为质心空间曲线1 设 z 0 即无质量偏心的情况 运动方程为设解为代入运动方程中 2020 4 9 9 并令得到特征方程为则所以代入边界条件 z 0 s 0 0 z l s l 0 解得 c2 c3 c4 0 c1sin kl 0要求非零解 则c1 0 所以sin kl 0因此有 kl n 得到固有频率为 2020 4 9 10 各阶主振型为 前三阶振型为2 设 z 0 即有质量偏心的情况 且质心按第n阶主振型函数 平面 分布 运动方程为设解为代入运动方程得根据假设 z 常数 则有 2020 4 9 11 式中 An为系数 sn z 为第n阶主振型由运动微分方程 得到设特解为Dn为待定系数代入运动方程得方程的齐次通解为sn z 且有故有特解方程为得到系数故转轴的振型为由此得到如下结论 2020 4 9 12 1 若质心按第n阶振型分布 只激起第n阶主振动2 转轴振型为一平面曲线 振幅为倍3 当 n时 振幅 产生第n阶主振型共振3 z 0 且质心为任意空间分布曲线 设为按主振型分解得即有质心分布示意图见图3 4所示 2020 4 9 13 式中代入运动方程有设转轴振型为代入运动方程得式中Sn z 为第n阶振型函数 也是对应齐次方程解所以有特解为 2020 4 9 14 利用固有振型的正交性 得解得系数转子振动为或 2020 4 9 15 三 柔性转子运动特点1 柔度曲线s z 随转速 而变变化1 0 6 c1 系数将增大 转子振型s z 是各阶主振型合成的空间曲线 3 cn时 第n阶主振型幅值系数明显增大 其它各阶则小很多 若 c1 此时振型近似有4 随着转速增加 各阶主振型依次突现出来 一般转子 主要是前三阶主振型的影响 比较挠度曲线与不平衡量的关系 它们展开项相同 幅值相差一个倍率 考虑阻尼有 2020 4 9 16 式中 cr 为无阻尼时系统的固有频率 r为挠度曲线各阶分量与该阶不平衡分量的相位差 由于阻尼影响 即使在临界转速下 转子振型也不是一根平面曲线 但实际进行动平衡时 仍以无阻尼的主振型平面加以考虑 3 转子主振型的正交性不平衡分布力在x y方向的分量为 2020 4 9 17 转子挠曲线在x y轴上的投影为各阶不平衡力在yoz平面和xoz平面上对k阶振型做功之和为由主振型正交性 2020 4 9 18 可知 1 各阶主振动之间不发生能量传递 2 n阶不平衡分量只能激起n阶主振型 不会激起其它各阶振型 3 利用主振型的正交性 可对转子进行逐阶平衡 完成柔性转子动平衡 2020 4 9 19 第三节模态平衡法 振型平衡法 一 模态平衡法及平衡条件根据主振型的正交性 可采用逐阶平衡的办法进行柔性转子动平衡 对于一般转子 主要是前三阶振型 以等截面轴为例进行分析 见图3 5设距起始端z1处有一集中重量w1位于半径R1上 集中重量均匀分布在2b的范围内 以U z 表示其分布 则式中 2020 4 9 20 取单位长度质量为m A 则有上式代表集中重量矩折合成单位长轴段质心偏移 按各阶主振型展开成式中 Cn1 n阶主振型系数 第二个下标表示所加平衡重量编号 sn z 各阶主振型函数 假设为已知 利用正交性 对折合轴段质心偏移展开式两边乘以sn z 并沿轴长积分 等式左边为 等式左边为 由此可得 2020 4 9 21 若在不同位置z1 z2 zk截面上 分别在半径R1 R2 Rk处加平衡配重W1 W2 Wk k个平衡重量引起转子质心的偏移为式中 为了平衡转子第n阶主振型分量 要求平衡重量形成的第n阶振型质心偏移和转子自身第n阶主振型质心偏移在同一平面上 大小相等 方向相反 即满足即若有一组k个最小的不平衡重量Uj 与n阶不平衡量相当 即 2020 4 9 22 式中 U z 转子不平衡量分布函数 其中 值应为最小 称这组量Uj j 1 k 为第n阶振型不平衡当量Une 即柔性转子的平衡不考虑阻尼情况下应满足下列三个力学平衡方程 2020 4 9 23 方程组中 第一 第二两式为刚性平衡条件 第三式为柔性平衡条件 二 配重面的选择及矢量平衡原理1 柔性转子平衡为多平面多转速平衡 2 平衡面选取 有N平面及N 2平面法两种 N平面法 平衡N阶振型 选用N个平衡面 N 2平面法 平衡N阶振型 选用N 2个平衡面 一般N平面法不能完全平衡支承动反力 但两种方法都有使用 平衡面选择很重要 选择不当将使平衡配重增大 原因 平衡面选择主要依据转子振型 实际发动机平衡面选择受到限制 图3 6为N 2平面法的平衡面选取 I II平衡面消除III IV V平衡面对低速动平衡的影响 2020 4 9 24 通常选择在紧靠支承的位置 以免影响高速时III IV V三个平面对振型不平衡量的校正 但由于在临界转速时 支承位移较大 I II平面的校正量对III IV V平面仍有一定干扰 图3 6 a 为平衡一阶振型时的三个平面的校正量 平面III的校正量对二阶振型不起作用 图3 6 b c 为平衡二阶及三阶振型的校正量组 测量柔性转子振型比较困难 可以轴承处的振动代替测量转子挠度 即矢量平衡法 图3 7为矢量平衡三角形 矢量为转子测点相对某一角向参考坐标测得的振动 矢量为转子上某点加试配重后同转速下测点与参考坐标下测得的振动 则矢量 为试重P的响应 为消除原始振动 加试配重平面上所需校正量为 2020 4 9 25 式中 称为影响系数矢量 用于影响系数法 称为反应系数矢量 用于模态平衡法 试重在原方位反时针旋转 角 其重量按OM对MN之比放大 即为校正量 平衡步骤 1 在第一阶临界转速附近测得两轴承处振动矢量 分解为对称矢量 该分量由一阶振型分量引起 2 加试配重后 在同一转速下测得振动 则矢量为试重引起的对称振动矢量 2020 4 9 26 3 平衡一阶振型分量的校正重量为 4 平衡二阶振型分量时 在二阶临界转速nc2附近测得两轴承振动及 其反对称分量为 它由二阶不平衡量引起 加反对称试重后 测得两轴承处的振动矢量为及 则矢量即为引起的反对称振动分量 故应加校正量为 2020 4 9 27 三 柔性转子平衡时的支承动反力柔性转子动平衡目的 1 消除支承动反力 2 消除转子挠度与弯矩 难于同时满足 则以最少的配重使转子在轴向 水平及垂直三方向振动在整个转速范围内最小 柔性转子挠曲振型为 设各阶振型函数为 简支梁情况 则转子振型为转子原始不平衡n阶分量可写成 则转子变形为 2020 4 9 28 支承动反力刚性部分由力矩平衡关系得设m z m 常数 等截面轴 上式积分整理得柔性部分支承动反力为积分整理得因此 一个轴承上所受到的总动反力为 2020 4 9 29 将代入得或由材料力学 通过振型函数求导得1 平衡一阶振型分量后的支承动反力设 简支梁 一阶振型分量为C1sin n l 其中则一阶挠曲振型为 2020 4 9 30 设采用位于中部的一个集中质量校正 即z1 l 2 校正量为W1R1 由 3 31 式得由于n 1 z l 2 故有若所选校正量满足C1 C11 即或此时 转子中部的一个校正量W1R1可以使一阶不平衡分量获得平衡 消除了柔性部分的动反力 转子一阶振型不平衡分量引起刚性部分的动反力为校正量W1R1加在中部后 一个支承上的动反力为 2020 4 9 31 比较 3 47 和 3 48 得可见 1 转子中部的一个集中平衡配重可使转子挠曲得以平衡 但不能全部消除转子的动力 2 支承处刚性部分动反力只能平衡掉78 5 3 为消除支承处一阶振型全部动反力 应在支承处同侧平面上加平衡配重W2及W3 若加重半径相同 则配重之比为 此时配重可消除转子一阶不平衡分量引起的转子挠曲 弯矩及支承动反力 但不能消除高阶振型分量影响 若工作转速远离高阶工作转速 则影响不大 若在转子两侧加校正量W1R1 W2R2 见图3 9 则有 2020 4 9 32 式中 z 1 当z1 b z z1 b及l z1 b z l z1 b z 1 z在其余各处 由 3 30 得由 3 31 得即设z1 l 4 取n 1得 为使两校正量能平衡一阶振型不平衡分量 应使C11 C1 即 2020 4 9 33 或两个校正量引起的支承处动反力为 由 3 47 可得同样为完全平衡刚性部分动反力 应在支承处两平面上相反方向再加两个校正重量W3及W4 若加重半径相同 则比值为 可见 平衡一阶振型不平衡量 转子中部加配重所需重量最少 平衡面越靠近两支承 需加配重越大 2 平衡二阶振型时的支承动反力设二阶不平衡分量为C2sin 2 z l 其中 由此引起转子二阶挠曲振型为 2020 4 9 34 引起两支承处的动反力大小相等 方向相反 为平衡二阶不平衡量应在两侧加反对称校正量W1R1 W2R2 如图3 10所示 则有式中 同理 利用正交性可得即 2020 4 9 35 取z1 l 4 n 2 则为平衡二阶振型 应使C22 C2 此时有校正量W2R2引起的刚性部分动反力为 由 3 42 得两者之比为 表明 加反对称校正量只能平衡支承动反力的78 5 为平衡全部刚性部分动反力 有力矩平衡方程得为平衡二阶振型分量 应使 2020 4 9 36 有以上两式解得 四 模态平衡的N平面法与N 2平面法1 转子偏心沿轴向及周向分布是随机的 2 1957年克劳尔 菲特恩 KlausFedern德 提出N 2平面法W 凯琳贝尔格 Kellenberge瑞士 加以充实 3 1959年R E D皮肖帕 Bishop 提出N平面法4 1970年皮肖帕和凯琳贝尔格对N平面法和N 2平面法进行了评论 基本获得共识 1 基本原理等截面运动方程为式中 运动方程也可用矢量形式表示为 2020 4 9 37 式中 m z 单位轴段长质量 作用在转子上不平衡力及校正力之和 设 0 可得齐次微分方程的通解 即特征值与特征向量 由主振型正交条件得将按sn z 展开后求运动微分方程的特解 由两部分组成 一是分布的不平衡力 另一部分是集中作用的不平衡力及校正力 其中为集中不平衡量及校正量 如图3 11 将按主振型展开 2020 4 9 38 且利用正交特性 可求得运动方程特解为式中第一项为式 3 22 第二项为集中不平衡量Uk引起 由此可知 当 接近某 cn时 转子的变形将是该阶主振型 其它各阶可不考虑 分析支承处动反力 由力矩平衡可得 2020 4 9 39 将特解代入得当 c1 或 0 上式与刚性转子相应式子完全相同 表明柔性转子在低速时 其特性与一般刚性转子相同 根据动平衡要求 应使支承动反力为零 2020 4 9 40 即由 3 62 得将上式矢量在y z及y x平面分解 仅讨论其中一个分量 2020 4 9 41 此时 矢量便成为标量 若保留高阶影响 仅消除前N阶振型 则有 可以得到满足上述条件的N个有限方程 2020 4 9 42 式中的集中质量分为两部分 M个原始存在的不平衡量 K个待确定的校正量 等式右边第一项是集中的不平衡量 第二项是转子连续分布的相应阶不平衡量 由于平衡面是有限的N个 满足式 3 66 式 3 64 就不能得到满足 即支承动反力不为零 因此 N平面法能减少转子变形 但不能平衡支承的动反力 影响了转子刚性平衡 因此 应补充式 3 64 可写为 因此 为了将式 3 66 与 3 67 一起求解 必需要有N 2个方程 即N 2平面法 2020 4 9 43 因此 N 2平面法由更高平衡精度 一般认为 N 3时 可采用N平面法 N 3时 宜采用N 2平面法 2 应用举例模型见图3 12首先在低于第一阶临界转速nc130 范围内进行刚性转子平衡 然后升速进行柔性转子平衡 取N 3 1 N平面法平衡面 I III V式 3 66 右端以 1 2 3表示 为未知数 由计算或测量得前三阶主振型 2020 4 9 44 s1 zk s2 zk s3 zk zk为三个加重面轴向位置 k 1 2 3 设振型已规格化 式 3 66 可写成将转速升到nc1附近 近似由 2 3 0 可得一阶振型平衡量组式中 A1 A2 A3 求解上述方程组的具体数值 相对值为平衡重量大小可用前述试加重量法求得 重量组配置见图3 13 a 平衡一阶振型后 然后升速至nc2附近 此时 同样近似有 1 3 0 同样可得第二组平衡量组 2020 4 9 45 二阶振型平衡量组相对值为若 U23 U21 则取 U23 为分母 平衡重量配置见图3 13 b 同理 升速至第三阶临界转速nc3附近 有 1 2 0 可得三阶振型平衡量组由此 可求得平衡量组相对值以上各平衡量组位于不同轴向平面上 须按矢量运算进行叠加 显然 所加的一组不平衡量并不满足 因此 N平面法部分地破坏了刚性转子平衡 2 N 2平面法 2020 4 9 46 平衡面数为N 2 5个 即图3 13中的I II III IV V共五个平面 其余步骤类似N平面法 联合式 3 66 及 3 67 得求解上述方程组可得U1 U2 U3 U4 U5 可分别令 2 3 0 1 3 0 1 2 0解上述方程组得三组U1 U2 U3 U4 U5 平衡量组配置见图3 13 显然 这是满足 因此 不破坏原先完成的 刚性转子 平衡 注意点 2020 4 9 47 1 加重面数不少于要平衡的主振型阶数 2 加重数值分配比例应满足正交条件 不产生附加的不平衡分量 即要求计算或测量得到准确的振型曲线 否则就不满足正交条件 2020 4 9 48 第四节影响系数法 对于线性系统 在一定转速下 不平衡与响应存在关系 式中 i处响应 j处不平衡量 影响系数 j点处不平衡量与i点处振动响应关系 为一矢量 且影响系数是转速的函数 与转速有关 测取影响函数 应使转速稳定 影响系数确定 首先在一定转速下 测取i点原始不平衡处响应 振幅及相位 然后 在j平面上加一已知不平衡量 再测取i点处的振动 由引起的振动为 即 2020 4 9 49 故影响系数为 上式按复数运算 设转子有M个平衡面 共有N个测点 则有式中 为转子原有的相当于再平衡面1 2 M上的不平衡量 测试转速为 1 再平面1上加不平衡量后 将转子启动 仍在 1下测取各点响应 记为 2020 4 9 50 两式相减得由此可得同理 在2平面上加试重可以求得 加M次试重 可求得所有M N个 幅值与相位 影响系数 改变转速为 2 重新测试 可得M N个影响系数 2020 4 9 51 二 动平衡方程组 矢量联立方程组 1 刚性转子支承在两支点的刚性转子 只需两个平衡面 分别在其上加试重可得到两支承 平衡面 处的影响系数若在一 二两个平衡面上应加校正量为 使得两支承处的振动为零 即满足方程式中 为第一 第二测点的原始振动值 由克莱姆法则 可得 2020 4 9 52 2 柔性转子设加重面为j 1 2 3 M 选定测点 包含所有转速下的测点 为i 1 2 3 N 通过加试重可得个测点影响系数 设各测点原始振动为 i 1 2 N 则应在平衡面加校正量后 各测点振动为零 由此得矢量联立方程组 用矩阵形式表示为 简写为 2020 4 9 53 柔性平衡转速通常不止一个 为使方程有精确解 必需使M N 而N H n 其中n为一个转速下的测点数 H为平衡转速数 通常平衡转速应包含临近临界转速及工作转速 由于加重面数M受到限制 一般N M 即方程数多于未知数 称为矛盾方程组 故不存在精确解 只能得到近似解 若N M 则为不定方程 无唯一解 平衡面选择 1 影响系数大 灵敏度高 2 测试精度高 重复性好 3 排除相关平衡面 4 关键 测量的影响系数必需准确 三 解复数矛盾方程组的最小二乘法1 在N个测量方程中 减去N M个方程 使得方程组成为正规方程组求解 缺点 难于满足大量舍去其它测点结果要求 甚至导致平衡状态恶化 2020 4 9 54 2 六十年代 P T Goodman提出 计算平衡校正量的最小二乘法 设想是将个测点剩余 残余 振动普遍下降 即加上一组平衡配重后 各测点振动平方和最小 特点 测点数不受限制 满足了柔性转子多平衡面多转速下进行平衡的要求 设原始振动为 各平衡转速下的影响系数为 j 1 2 M 为平衡面编号 i 1 2 N 为测点编号 待求校正量为 令 i 1 2 N 为各测点剩余振动 在线性条件下 满足方程组 2020 4 9 55 简写为 令s为剩余振动的平方和 则式中 为的共轭复数 若一组校正量满足上式 则为最有近似解 由上式可知 s是自变量 j 1 2 M 的二次函数 因此 上式为求解二次函数s的极值问题 由于s是的连续函数 且故必存在一组 或使得s值为最小 即满足 2020 4 9 56 由于 则有将实部与虚部展开得代入极值条件中可得 2020 4 9 57 由上两式得将式 3 82 代入 并整理得此即为由最小二乘法导出的对应于矢量矛盾方程组的复数正规方程组 写成矩阵形式为 由此可知 将复数矩阵的共轭转置矩阵左乘原矛盾方程组 可直接求得对应复数正规方程组 即 2020 4 9 58 设为矛盾方程组近似最优解 可用矩阵表示为 最后 用所求结果代入 3 78 式 校验剩余振动值是否在允许范围内 四 加权迭代在选定配重方案下总的剩余振动为若存在剩余振动偏大 即 max R 可进行加权迭代 使得最大剩余振动下降 但其它测点剩余振动有可能增加 但剩余振动平方和以不加权迭代的数值最小 设N个测振点的加权因子为 i 分别乘以 3 79 两边得 2020 4 9 59 由最小二乘法得其中加权因子根据经验得出 一般推荐 如果剩余振动过大 可进行第二次迭代 一般迭代1 2次即可 因此 影响系数法平衡柔性转子的步骤为 1 确定平衡面数 2 给定转速下 测得N个测点的原始不平衡振动值 3 依次在M个平衡面上加已知试重 在同样给定转速下测出N个测点的振动值 4 根据测量结果 编程 计算影响系数 5 利用最小二乘法或加权最小二乘法求解矛盾方程组 得到校正量 6 校核剩余振动 满足即结束 若不满足 进入第5 步骤 2020 4 9 60 模态平衡法与影响系数法的比较 如果平衡转速选择在N个临界转速附近 并选M N个校正面 一个测振点 校正平面zj试重qj引起的振动为当平衡转速 接近临界转速 ck时 上式中以k阶响应项为主 其它项可忽略 此时影响系数为 此时影响系数矩阵的每一行与sk z1 sk z2 sk zM 成比例 此时 两种方法等效 若平衡转速不选在临界转速附近 则结论不成立 2020 4 9 61 第五节柔性转子其它平衡方法 一 混合平衡法 振型影响系数法 1 影响系数法优点 可同时平衡几阶振型 平衡方便 可利用计算机辅助平衡 便于实现数据处理自动化 缺点 高速下平衡启动次数多 高阶振型敏感性降低 可能存在相关平面 得到不正确的校正量 2 模态平衡法优点 高速平衡启动次数少 敏感性高 低阶振型不受影响 缺点 振型不易侧准 系统阻尼大时不够有效 对于轴系平衡时 在临界转速附近不易获得单一振型 3 混合平衡法结合两者优点 在影响系数法基础上 充分利用模态平衡 2020 4 9 62 法中的振型分离特点选择各项系数 效果较好 以图3 14所示转子为例设在转速范围内出现两阶临界转速 按模态平衡法中N 2平面法 选取四个校正面 即N 4 校正面按模态平衡法选取 如图中1 2 3 4四个平面 由振型平衡方程可得 2020 4 9 63 式中 简写为上式校正量有唯一解 但是随机分布的 无法直接求解出校正量 利用影响系数 可列出四个线性方程若选取合适 上述方程组有唯一解 且满足式 3 91 由式 3 92 得式中 2020 4 9 64 若值满足式 3 95 则式 3 93 的解必满足式 3 91 取两支点为测点 按振型分离步骤升速 先进行低速动平衡 转速为 0 由式 3 95 得式中上角标 0 表示原始状态 低速平衡 振型不平衡分量作用很小 上式可简化为 可求得两支承振动校正量为 2020 4 9 65 由此得到校正量 然后升速至 1 进行一次高速动平衡 测得两支承振动为 式中 2020 4 9 66 因 1接近一阶临界转速 主要为一阶振型不平衡分量的作用 上式可近似简化为方程组线性相关 故只能取一个振动值 若 一般取较大的值 按N 2平面法 校正重由下式获得式中 为相应于加一次配重后在转速 0时二支承的振动 该方程相容 满足平衡方程式 3 91 中前三个方程 即转子可以平稳地通过一阶临界转速 又不破坏刚性转子平衡 同理 将转速升到二阶临界转速附近 进行高速动平衡 平衡转速为 2 二支点振动可近似简化为 2020 4 9 67 同样 方程组线性相关 只能取振动较大值 此时校正量为四个 可由下式求得式中 相应为转子加上二次配重后二支承振动 该方程满足平衡方程式 3 91 中的四个方程 最后得到总校正量为该校正量可使转子在全速范围内达到平衡 2020 4 9 68 二 振型圆平衡法 模态响应圆平衡法 模态平衡法要求在各阶临界转速下进行平衡 但会引起过大振动 且不稳定 振型圆平衡法 利用模态分析技术 通过振型圆 在临界转速时分离出该阶模态 尤其适用于多跨轴系 1 基本原理振型圆 模态圆 已知Jeffcott转子运动方程为 式中 z x iy 令n c 2m 其稳态解为 2020 4 9 69 由式 3 100 有由此得故有上式为复平面上的振型圆方程 圆心在虚轴一侧并与实轴相切 见图3 15 a 经过圆心和临界转速点的直径OA称为共振直径 实际振型圆不封闭 见图3 15 b 特点 2020 4 9 70 1 具有一般模态圆特征 2 在临界转速附近 相位变化剧烈 当 c时 ds d 03 c时 2 表明不平衡激振力垂直共振直径 并超前 因而确定了不平衡激振力方向 4 单自由度系统 矢端圆近似为一个圆 多自由度系统 图形比较复杂 但在各阶临界转速附近 矢端图仍接近一个圆 图3 16为一双支点对称转子轴承测得的振型圆图 不平衡位于平面1 轴承A 一侧 靠近不平衡量一侧的两阶临界转速时的振型圆是同相的 另一侧两个振型圆则是反相的 2 平衡方法影响因素 其它振型 初始弯曲 热变形等 1 A法 根据预先计算结果 在临界转速附近测得转子或轴承的振幅 相位 绘制振型图 见图3 17 2020 4 9 71 2 B法 借试重求出影响矢量 特点 1 测量精度高 数据可靠 2 测量速度快 不需稳速测量 3 平衡精度高 可有效分离振型 4 需具备多点连续自动检测和分析装置 三 限制最大平衡配重法高速微型发动机转子 由于结构限制 平衡配重受到一定限制 由此 对转子配重的最大值做出限制 以便在规定平衡配重的条件下达到平衡精度要求的目的 1 首先 选择可用于平衡的面 进行转子动力学分析 剔除相关平面 确定平衡面数N 2 对最小二乘法进行改进 限定最大配重量 通过优化配置 如遗传算法等 获得在限定最大配重条件下的最优 2020 4 9 72 平衡配重组合 3 校核平衡结果 四 谐波小波在柔性转子平衡中的应用1 动平衡信号的特点不平衡引起的振动信号主要是与转速同频率的正弦波 是振动信号的主要成分 实际测量信号中可能还有2 3 4等倍频 l 2的亚倍频 随机振动等成分 因此实际信号可以用下面的模型来表示 3 4 1 式中e t 是测量实际信号 E0是测量仪输出的直流信号 是不平衡量引起的基频信号 是各谐波分量 2 谐波小波滤波原理谐波小波为复小波 在频域紧支 有明确的函数表达式 2020 4 9 73 3 4 2 图3 4 1谐波小波时域图图3 4 2谐波小波频域图根据小波理论对谐波小波进行伸缩 平移就生成谐波小波函数族及其伸缩平移函数族构成信号的正交基 因而 将谐波小波作为基函数系可以将信号既不重叠 又无遗漏地分解到相互独立的频域空间 重构转速频段的信号即实现了滤波功能 图3 4 3二进谐波小波对频域空间的划分 2020 4 9 74 二进谐波小波在高频段存在频段带宽太大的缺点 需要对其进行改进 写为更一般的形式 3 4 3 其中m n可以不是整数 取一定步长k n m 上式变为 3 4 4 定义一对复小波系数 3 4 5 则离散信号f r 在频段2m 2n 的重构为 3 4 6 2020 4 9 75 不平衡量信号是与转速同频的 因此可以根据转速信息来调整m和n的值 从而达到控制通带的带宽和频率中心的目的 完全符合自适应滤波的要求 m和n的值由下式确定 3 4 7 其中N为采样点数 fs为采样频率 fmin和fmax分别为带通滤波器下截止频率和带通滤波器上截止频率 3 仿真试验为了验证上述方法的有效性 根据式 3 4 1 建立含噪不平衡量信号的模型 3 4 8 不平衡量信号频率f1 480Hz 倍频f2 960Hz 噪声n x 是均值为0 方差为1的高斯白噪声信号 采样频率12 2kHz 采样长度1024 滤波带宽32Hz 下截止频率464Hz 上截止频率496Hz 滤波试验结果见图3 4 4 2020 4 9 76 图3 4 4滤波试验结果图 2020 4 9 77 由图3 4 4 a 可见不平衡量信号完全淹没在噪声信号中 无法获得信号幅值和相位信息 图3 4 4 b 和图3 4 4 c 分别为无噪声干扰的不平衡量信号和谐波小波滤波后的信号 两图比较可以看出 推广的谐波小波滤波能够有效地消除噪声干扰 较好地还原出了不平衡量信号 幅值存在微小误差 相位无滞后 2020 4 9 78 第六节柔性转子平衡精度 1 刚性转子 国际标准ISO 1940 采用偏心速度评定 偏心速度Ve可表示为 Ve e 1000 毫米 秒 式中 e 不平衡偏心 微米 其值为其中 URA 转子允许剩余不平衡量 克 毫米 M 转子质量 公斤 乘积e 代表了转子转速为 时的质心速度 国际标准中按不同的e 值将平衡精度分为11个等级 各等级之间公比为2 5 对于燃气轮机转子 标准规定平衡精度应达到G2 5级 即偏心速度Ve 1 2 5毫米 秒 转速较高的转子 Ve相应取较小值 确定偏心速度Ve后 可按下式求得转子允许的剩余不平衡量URA 即 2020 4 9 79 其中 取最大工作转速 另外 还可以从作用在轴承上的不平衡力F与转子重量W之比值