一种发动机振动不平衡相位的实时测量方法
【技术领域】
[0001] 本发明属于航空发动机振动监测技术领域，具体涉及一种发动机振动不平衡相位 的实时测量方法。
【背景技术】
[0002] 航空发动机的工作条件十分恶劣，发动机健康状况关系到飞机及乘坐人员的安 全。对发动机的状态和健康状况进行监控可以保障飞机及人员安全的有效手段。振动采集 监控是发动机状态监控的重要功能，作为典型的旋转机械设备，发动机振动幅值、频率、相 位反映了发动机状态的重要信息。
[0003] 进行发动机振动信号采集监测有三方面的应用：振动超限监测、发动机故障诊断 及发动机转子配平。振动超限监测实时监测发动机振动大小，当发动机健康状况出现变化， 反映到发动机振动水平上时，可以提前采取相应检查及维护措施，避免造成故障的恶化及 重大事故;发动机部件出现不同的损伤或故障，其振动信号有相应的振动特征频谱。利用发 动机不同频次的振动信号可以实现基于发动机振动信号的发动机故障诊断;发动机运行在 较大振动水平下会对发动机的健康状况带来损伤，当发动机振动水平较大且并未出现故 障，可以对发动机进行配平，降低发动机振动水平。
[0004] 对发动机配平，常用方法是三圆法，三圆法不需要求取发动机振动不平衡相位，但 三圆法需要发动机开车次数较多，耗时耗力也严重影响发动机的运营效率。求取发动机转 子不平衡量的幅值和相位可以实现发动机振动的快速动平衡，直接计算出平衡重量的大小 和安装角度。转子振动不平衡相位通过转子转速相位基准信号和振动信号来计算。
[0005] 对振动信号直接进行FFT运算实现傅里叶变换，求取发动机振动信号的幅值和相 位，由于非整周期采样会存在频谱泄露，使得求取的振动信号幅值和相位值偏离实际值，频 谱泄露对相位值的影响最大可达90°。理论上可以实时调整振动信号的采样频率或进行插 值运算保证参与FFT运算的信号满足整周期条件，进而避免频谱泄露。但因为发动机转速是 连续变化的，根据转速信号实时调整采样频率，存在信号延迟及使得系统复杂的问题。使用 数学插值也因为发动机转速连续变化，存在插值计算精度有限的问题。
[0006] 另一方面，发动机振动信号中除高低压转子基频、倍频信号外，还包含大量的噪声 信号。受发动机转子基频振动外的无关信号影响，使得使用FFT算法求取振动信号相位精度 较差，难以满足发动机动平衡的需求。
【发明内容】
[0007] 本发明的目的在于克服上述振动不平衡相位测量方法的不足，提出了一种基于互 功率谱估计的航空发动机振动不平衡相位测量方法，以避免由于频谱泄露、噪声等因素的 影响，提高测量精度。
[0008] 为实现上述目的，本发明的计算方法如下：
[0009] 该发动机振动不平衡相位的实时测量方法通过检测相位基准信号中的脉冲，生成 相位基准参考信号，求取相位基准参考信号与振动信号的互功率谱;然后根据互功率谱寻 找发动机振动不平衡相位频域计算点，利用互功率谱实部、虚部估计振动不平衡相位。
[0010] 具体包括以下步骤：
[0011] 1)检测相位基准信号中的脉冲，以脉冲为起始，提取长度为N的相位基准信号序列
和振动信号序列，计算相位基准信号序列中相邻两个脉冲之间Tk， 正弦信号，拼接成长度为N的正弦信号序列，N=MAX(2n)&2n<f S，fs为所述脉冲信号的采样 频率；
[0012] 2)对步骤1)产生的正弦信号序列进行加汉宁窗处理，然后进行点数为N的快速傅 里叶变换得到x(k);对步骤1)提取的振动信号序列进行加汉宁窗处理，然后进行点数为N的 快速傅里叶变换得到Y(k);
[0013] 3)按下式计算求取互功率谱密度Gxy;
[0014] Gxy = X*(k) · Y(k)
[0015] 这里的*的含义是"求实部"。
[0016] 4)对步骤3)求取的互功率谱，按下式找出幅频曲线中发动机转子转速频率附近的 幅值最大点，并记录该频率点k:
[0018] An,P为转子基频处的振动加速度峰值，fni为转子基频，N为FFT变换的点数，floor (X)表示小于X的最大整数；
[0019] 5)根据下式利用步骤4)记录的频率点k计算所得相位即为发动机转子振动不平衡 相位
[0020] Pk = atan2(imag(GXy(k))，real(GXy(k)))〇 [0021 ]本发明与现有技术相比，有以下优点：
[0022] 1.通过求取相位基准信号与振动信号的互功率谱，将振动信号中与发动机转子同 频的信号加强，噪声和其它频率成分都得到了抑制，提高了振动不平衡相位的估计精度； [0023] 2.对振动信号进行加窗FFT求与参考信号的互功率谱，振动信号的FFT结果也可用 于转子振动超限检测，避免了重复计算，节约了机载设备的计算资源；
[0024] 3.本发明关键运算FFT，在DSP等嵌入式系统中均有成熟工程化的模块可调用，算 法实现简单，适合实时实现。
[0025] 4.运算量小、稳定性好、易于工程实现。
【附图说明】
[0026] 图1是本发明方法的应用环境框图；
[0027] 图2是振动不平衡相位求取流程图；
[0028]图3是由相位基准脉冲信号生成同频正弦参考信号的示意图；
[0029]图4是汉宁窗时域特性图；
[0030]图5是正弦参考信号与振动信号互功率谱幅频图；
[0031]图6是正弦参考信号与振动信号互功率谱相频图；
[0032] 图7是使用本发明方法求取某次发动机试车振动不平衡相位的结果。
【具体实施方式】
[0033] 下面通过具体的实施方式对本发明作进一步的详细说明。
[0034]本发明包括的处理步骤及相关参数如下：
[0035]步骤1，记录相位基准脉冲信号，采样频率fs为10K，以脉冲信号为起始，提取长度 为8192点的相位基准信号和转子振动信号，计算相位基准信号中相邻两个脉冲之间的时间
为频率生成正弦信号，拼接成长度为8192点的正弦信号序列。
[0036]步骤2,按照下式生成8192点的汉宁窗函数序列，分别与正弦参考序列和振动信号 序列进行对应项相乘，汉宁窗幅值修正系数C = 0.5。图4为汉宁窗的时域特性图。
[0038] 步骤3,对加窗后的正弦参考序列和振动信号序列分别进行8192点的FFT，对FFT结 果除以汉宁窗幅值修正系数C，得到正弦参考序列傅里叶变换X(k)和振动信号序列傅里叶 变换Y(k);
[0039] 步骤4,按照下式根据正弦参考序列和振动信号序列FFT结果求取正弦参考信号与 振动信号的互功率谱；
[0040] Gxy = X*(k) · Y(k)
[0041] 步骤5,对步骤4求取的互功率谱，按如下公式找出幅频曲线中发动机转子转速频 率附近的幅值最大点，并记录该频率点k:
[0043] An,P为转子基频处的振动加速度峰值，fni为转子基频，N为FFT变换的点数，符号 floor (X)表示小于X的最大整数。
[0044] 步骤6,根据下式利用步骤5记录的频率点k计算所得相位就为发动机转子振动不 平衡相位。
[0045] Pk = atan2(imag(Gxy(k)) ,real(Gxy(k)))
[0046] 循环执行以上步骤，实现对发动机振动不平衡相位的实时测量。
[0047] 图7为在某次发动机振动试车过程中，试用本发明计算求取的发动机转子振动不 平衡相位，并且使用该相位结果，制定发动机动平衡方案，对发动机进行了振动配平，发动 机振动量级得到有效改善，进而验证了本发明方法的正确性和工程适用性。
[0048]综上，本发明综合考虑发动机振动超限监测和振动配平的需求，计算转速相位基 准参考信号和振动信号的互功率谱求取振动不平衡相位，使得实现快速动平衡成为可能， 可大大缩短振动配平时间，减少工作量，提高了飞机运营效率。
【主权项】
1. 一种发动机振动不平衡相位的实时测量方法，其特征在于:检测相位基准信号中的 脉冲，生成相位基准参考信号，求取相位基准参考信号与振动信号的互功率谱;然后根据互 功率谱寻找发动机振动不平衡相位频域计算点，利用互功率谱实部、虚部估计振动不平衡 相位。2. 根据权利要求1所述的发动机振动不平衡相位的实时测量方法，其特征在于，具体包 括W下步骤： 1) 检测相位基准信号中的脉冲，W脉冲为起始，提取长度为N的相位基准信号序列和振 动信号序列，计算相位基准信号序列中相邻两个脉冲之间化，然后WlAk生成同频的正弦信 号，拼接成长度为N的正弦信号序列，N=MAX(2D化2n<fs，fs为所述脉冲信号的采样频率； 2) 对步骤1)产生的正弦信号序列进行加汉宁窗处理，然后进行点数为N的快速傅里叶 变换得到X化）；对步骤1)提取的振动信号序列进行加汉宁窗处理，然后进行点数为N的快速 傅里叶变换得到Y化）； 3 )按下式计算求取互功率谱密度Gxy ; Gxy=料化）· Y化） 4) 对步骤3)求取的互功率谱，按下式找出幅频曲线中发动机转子转速频率附近的幅值 最大点，并记录该频率点k:Αη,ρ为转子基频处的振动加速度峰值，fni为转子基频，賦JFFT变换的点数，floor(x)表 示小于X的最大整数； 5) 根据下式利用步骤4)记录的频率点k计算所得相位即为发动机转子振动不平衡相位 Pk = a1:an2(imag(Gxy(；k)) 'real (Gxy(；k)))。
【专利摘要】本发明提出一种基于互功率谱估计的发动机振动不平衡相位的实时测量方法，主要解决航空发动机振动超限配平振动不平衡相位求取问题。本发明通过求取相位基准信号与振动信号的互功率谱，将振动信号中与发动机转子同频的信号加强，噪声和其它频率成分都得到了抑制，提高了振动不平衡相位的估计精度；对振动信号进行加窗FFT求与参考信号的互功率谱，振动信号的FFT结果也可用于转子振动超限检测，避免了重复计算，节约了机载设备的计算资源；其中的关键运算FFT，在DSP等嵌入式系统中均有成熟工程化的模块可调用，算法实现简单，适合实时实现。
【IPC分类】G01M15/12
【公开号】CN105571874
【申请号】CN201510925852
【发明人】冯健朋, 马振华, 马萌, 樊芊
【申请人】中国航空工业集团公司西安航空计算技术研究所
【公开日】2016年5月11日
【申请日】2015年12月11日