CN114812794A - 一种扫频速率引起的叶尖定时参数辨识误差修正方法 - Google Patents
一种扫频速率引起的叶尖定时参数辨识误差修正方法 Download PDFInfo
- Publication number
- CN114812794A CN114812794A CN202210029094.XA CN202210029094A CN114812794A CN 114812794 A CN114812794 A CN 114812794A CN 202210029094 A CN202210029094 A CN 202210029094A CN 114812794 A CN114812794 A CN 114812794A
- Authority
- CN
- China
- Prior art keywords
- sweep
- frequency
- blade
- parameter
- eta
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 54
- 230000004044 response Effects 0.000 claims abstract description 33
- 238000013016 damping Methods 0.000 claims abstract description 12
- 238000012937 correction Methods 0.000 claims abstract description 10
- 230000001133 acceleration Effects 0.000 claims abstract description 8
- 210000003918 fraction a Anatomy 0.000 claims abstract description 4
- 230000001360 synchronised effect Effects 0.000 claims abstract description 4
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000010408 sweeping Methods 0.000 abstract description 2
- 238000005259 measurement Methods 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000036039 immunity Effects 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005316 response function Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/72—Wind turbines with rotation axis in wind direction
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Turbine Rotor Nozzle Sealing (AREA)
Abstract
本发明涉及一种扫频速率引起的叶尖定时参数辨识误差修正方法,针对使用单参数法辨识叶片共振参数的参数辨识方法进行误差修改,执行如下步骤:获取变速扫频时转子旋转加速度;利用单参数法辨识叶片同步振动参数,包括扫频响应峰值振幅Asweep、扫频响应峰值频率fsweep、扫频响应品质因数Qsweep;使用迭代法,求取无量纲扫频参数η值;计算峰值振幅分数A(η)、归一化频率误差f(η)和阻尼比分数ζ(η)的值;计算误差修正后的叶片峰值振幅、峰值频率和品质因数。
Description
技术领域
本发明涉及基于一种扫频速率引起的叶尖定时参数辨识误差修正方法。
背景技术
叶片振动监测对保障旋转机械的正常运行非常重要,叶尖定时是一种有效地旋转叶片振动测量技术。 利用叶尖定时,能实现叶片振动事件辨识和故障诊断。与传统的应变片测量相比,叶尖定时具有非接触、 低介入,能实现全叶片测量的优点。但是,叶尖定时是一种欠采样测量方法,如何通过欠采样信号重构叶 片振动信息,辨识叶片振动参数是叶尖定时的研究热点和难点。利用叶尖定时数据进行叶片振动参数辨识 的主要方法有周向傅里叶拟合、单参数法、双参数法、自回归法、空间类算法及近年提出的稀疏重构和压 缩感知的方法。其中,单参数法和基于单参数法的多传感器振动倍频辨识方法具有良好的鲁棒性和抗噪性, 在工程上得到广泛应用。
准确辨识叶片振动参数对叶片状态监测及故障诊断尤为重要。以裂纹故障为例,其引起的叶片共振频 率变化为2~3Hz左右,这就要求叶片振动参数辨识精度要优于此值。在叶尖定时测振技术中,通过变速 扫频获取叶片振动参数是一种通用且有效的方法。变速扫频是改变转子转速,使叶片激振力频率发生变化, 从而叶片振动位移响应发生相应的变化。单支叶尖定时传感器可获取不同转速下叶尖振动位移值,使用单 参数法进行曲线拟合,可获取叶片共振幅值、共振转速、相位、常偏量等参数。再利用多支叶尖定时传感 器(一般3支以上)测得的叶片振动相位差,实现叶片振动倍频数和振动频率的辨识。
推导单参数法曲线拟合的目标函数时,使用了单自由度系统的稳态响应结果近似表示了叶片全响应。 当单自由度系统阻尼较小(即品质因数较大),或扫频速率(转速加速度)较大时,单自由度系统瞬态响 应则不能忽略。由于稳态和瞬态响应的叠加,会使得观测到的系统响应峰值振幅减小,峰值频率偏移(偏 移方向与扫频方向相同)。严重时,还会出现“振铃”现象,形成多个共振峰。过大的扫频速率会造成实际振 动信号包络与幅频响应曲线偏差较大,使用单参数法辨识叶片振动参数辨识时,将产生较大误差。如何减 小扫频速率引起的叶尖定时参数辨识误差,实现叶片振动参数的高精度辨识,是现有技术中待解决的问题。
参考文献;
[1]Peter Russhard.Development of a Blade Tip Timing Based EngineHealth Monitoring System[D].University of Manchester(United Kingdom),2010.
[2]欧阳涛,郭文力,段发阶,李孟麟.基于叶尖定时的旋转叶片同步振动辨识新方法[J].振动与冲 击,2011,30(08):249-252+257.
发明内容
本发明的目的是针对扫频速率引起的叶尖定时参数辨识误差问题,引入了无量纲扫频参数η和扫频响 应特性函数,给出了扫频响应特性函数与无量纲扫频参数η的关系曲线和拟合公式,进而利用迭代法实现 了无量纲扫频参数η的求解。基于叶尖定时单参数法的振动参数辨识结果,提出一种扫频速率引起的叶尖 定时参数辨识误差修正方法,可用于叶尖定时单参数法参数辨识误差修正,实现叶片共振参数的高精度辨 识。技术方案如下:
一种扫频速率引起的叶尖定时参数辨识误差修正方法,针对使用单参数法辨识叶片共振参数的参数辨 识方法进行误差修改,执行如下步骤:
(1)获取变速扫频时转子旋转加速度a,单位为rpm/s;
(2)利用单参数法辨识叶片同步振动参数,包括扫频响应峰值振幅Asweep、扫频响应峰值频率fsweep、 扫频响应品质因数Qsweep;
(3)使用迭代法,求取式(1)中无量纲扫频参数η值:
(4)将无量纲扫频参数η值计算结果带入式(2)和式(3),计算峰值振幅分数A(η)、归一化频率误 差f(η)和阻尼比分数ζ(η)的值:
(5)将(2)和(3)得到的结果带入式(4),计算误差修正后的参数值;
其中,A0、f0、Q0即为修正后的叶片峰值振幅、峰值频率和品质因数。
本发明由于采取以上技术方案,其具有以下优点:
(1)本发明引入了无量纲扫频参数η和扫频响应特性函数,给出了扫频响应特性函数的拟合公式, 利用迭代法实现了无量纲扫频参数η的求解。
(2)本发明基于叶尖定时单参数法的振动参数辨识结果,提出一种扫频速率引起的叶尖定时参数辨 识误差修正方法,实现了叶片振动参数的高精度辨识。
附图说明
图1为扫频响应特性函数峰值振幅分数曲线图。
图2为扫频响应特性函数归一化频率误差曲线图。
图3为扫频响应特性函数阻尼比分数曲线图。
图4为叶尖定时参数辨识误差修正方法流程图。
图5为变速扫频叶片振动响应曲线。
图6为变速扫频传感器采样曲线。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
首先对一些参数和背景知识进行说明。
扫频速率,即旋转机械转速加速度;叶尖定时,利用若干支安装在机匣上的叶尖定时(BTT)传感器 和1支安装在转轴上的转速同步(OPR)传感器,测量叶片共振参数;参数辨识,使用单参数法[1-2]辨识叶 片共振参数;误差修正,利用修正公式修正单参数法的参数辨识误差。
为实现扫频速率引起的叶尖定时参数辨识误差修正,引入无量纲扫频参数η和扫频响应特性函数A(η), f(η),ζ(η),定义为:
其中,a为扫频速率(即转子加速度,单位为rpm/s)。A(η),f(η),ζ(η)分别称为峰值振幅分数、归一化 频率误差、阻尼比分数,其曲线图分别如图2、图3和图4所示。Asweep,fsweep,ζsweep,Qsweep分别为变速扫频 时利用单参数法辨识得到的扫频响应峰值振幅、扫频响应峰值频率、扫频响应阻尼比和扫频响应品质因数; A0,f0,ζ0,Q0为叶片真实的峰值振幅、峰值频率、阻尼比和品质因数。
对A(η),f(η)进行曲线拟合,拟合公式为:
由单自由度系统振动理论可知,叶片共振峰值振幅与阻尼比成反比,故阻尼比分数与峰值振幅分数亦 成反比,有:
通过传统单参数法辨识得到一组振动参数[1,2],设辨识得到的扫频响应峰值振幅、扫频响应峰值频率、 扫频响应品质因数分别为Asweep,fsweep,Qsweep.。由式(5)和式(6)可得,无量纲扫频参数η可表示为:
上式简化为:
将式(7)和式(8)带入式(10),无量纲扫频参数η求解表达式可表示为:
基于单参数法的辨识结果,利用迭代法,根据式(11)可求解无量纲扫频参数η。修正后的振动参数值计 算公式如式(12)所列。
扫频速率引起的叶尖定时参数辨识误差修正方法流程如图4所示,具体实施过程如下:
(1)模拟某转子实验台叶片质量m=1kg,阻尼c=8Ns/m,刚度k=64×104N/m。设变速扫频范围 为6600~9000rpm(110~150Hz),扫频加速度a=600rpm/s,经过了叶片一阶共振频率。叶片振动响应 和叶尖定时传感器采样曲线图5和图6所示。使用单参数法通过变速扫频辨识得到振动参数为:扫频响应 峰值振幅Asweep=0.7885mm,扫频响应峰值频率fsweep=128.9069Hz,扫频响应品质因数Qsweep=70.1431。
(2)基于(1)辨识得到的叶片振动参数和转子加速度,使用迭代法,求取式(11)中无量纲扫频参数η 值为5.6710。
(3)将无量纲扫频参数η值计算结果带入式(7)。计算得到峰值振幅分数A(η)=0.7332和归一化频率 误差f(η)=1.4075。
表1
(4)将(1)和(3)得到的结果带入式(12),得到误差修正后的峰值振幅A0=1.0799mm,峰值频率 f0=127.0456Hz,品质因数Q0=96.0684。计算修正前和修正后各参数的相对误差,如表1所列。相比修 正前,修正后的参数相对误差明显减小。
Claims (1)
1.一种扫频速率引起的叶尖定时参数辨识误差修正方法,针对使用单参数法辨识叶片共振参数的参数辨识方法进行误差修改,执行如下步骤:
(1)获取变速扫频时转子旋转加速度a,单位为rpm/s;
(2)利用单参数法辨识叶片同步振动参数,包括扫频响应峰值振幅Asweep、扫频响应峰值频率fsweep、扫频响应品质因数Qsweep;
(3)使用迭代法,求取式(1)中无量纲扫频参数η值:
(4)将无量纲扫频参数η值计算结果带入式(2)和式(3),计算峰值振幅分数A(η)、归一化频率误差f(η)和阻尼比分数ζ(η)的值:
(5)将(2)和(3)得到的结果带入式(4),计算误差修正后的参数值;
其中,A0、f0、Q0即为修正后的叶片峰值振幅、峰值频率和品质因数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210029094.XA CN114812794B (zh) | 2022-01-12 | 2022-01-12 | 一种扫频速率引起的叶尖定时参数辨识误差修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210029094.XA CN114812794B (zh) | 2022-01-12 | 2022-01-12 | 一种扫频速率引起的叶尖定时参数辨识误差修正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114812794A true CN114812794A (zh) | 2022-07-29 |
CN114812794B CN114812794B (zh) | 2024-03-22 |
Family
ID=82527694
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210029094.XA Active CN114812794B (zh) | 2022-01-12 | 2022-01-12 | 一种扫频速率引起的叶尖定时参数辨识误差修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114812794B (zh) |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103917851A (zh) * | 2011-09-08 | 2014-07-09 | 三菱重工业株式会社 | 误差频率分量获取装置、旋转角度获取装置、电机控制装置以及旋转角度获取方法 |
US20140365166A1 (en) * | 2013-05-30 | 2014-12-11 | Rolls-Royce Plc | Blade tip timing |
CN104697623A (zh) * | 2014-10-24 | 2015-06-10 | 南京航空航天大学 | 变速激励下叶片异步振动参数识别方法 |
CN105547619A (zh) * | 2015-12-04 | 2016-05-04 | 东北大学 | 一种薄壁构件高阶模态频率及阻尼测试方法及系统 |
WO2017187333A1 (en) * | 2016-04-26 | 2017-11-02 | University Of Pretoria | A method and system for monitoring rotor blades of a turbomachine using blade tip timing (btt) |
EP3290653A1 (en) * | 2016-09-02 | 2018-03-07 | United Technologies Corporation | Dynamical system parameter identification for turbomachine |
CN108871543A (zh) * | 2018-08-30 | 2018-11-23 | 天津大学 | 恒速下叶片异步振动频率非调和傅里叶分析方法 |
CN109790757A (zh) * | 2016-06-27 | 2019-05-21 | 比勒陀利亚大学 | 使用叶片尖端定时(btt)测量转子叶片尖端偏转的方法和系统 |
CN111579060A (zh) * | 2020-06-28 | 2020-08-25 | 大连理工大学 | 一种基于叶尖定时技术高精度的旋转叶片振动测量方法 |
CN112179480A (zh) * | 2020-08-26 | 2021-01-05 | 西安交通大学 | 一种谐波叠加拟合叶尖定时振动参数识别的方法、系统及电子设备 |
CN112903271A (zh) * | 2021-01-22 | 2021-06-04 | 中国航发沈阳发动机研究所 | 一种转子叶片非接触式异步振动参数辨识方法 |
CN113586177A (zh) * | 2021-08-10 | 2021-11-02 | 西安交通大学 | 基于单叶端定时传感器的叶片固有频率识别方法 |
-
2022
- 2022-01-12 CN CN202210029094.XA patent/CN114812794B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103917851A (zh) * | 2011-09-08 | 2014-07-09 | 三菱重工业株式会社 | 误差频率分量获取装置、旋转角度获取装置、电机控制装置以及旋转角度获取方法 |
US20140365166A1 (en) * | 2013-05-30 | 2014-12-11 | Rolls-Royce Plc | Blade tip timing |
CN104697623A (zh) * | 2014-10-24 | 2015-06-10 | 南京航空航天大学 | 变速激励下叶片异步振动参数识别方法 |
CN105547619A (zh) * | 2015-12-04 | 2016-05-04 | 东北大学 | 一种薄壁构件高阶模态频率及阻尼测试方法及系统 |
WO2017187333A1 (en) * | 2016-04-26 | 2017-11-02 | University Of Pretoria | A method and system for monitoring rotor blades of a turbomachine using blade tip timing (btt) |
CN109790757A (zh) * | 2016-06-27 | 2019-05-21 | 比勒陀利亚大学 | 使用叶片尖端定时(btt)测量转子叶片尖端偏转的方法和系统 |
EP3290653A1 (en) * | 2016-09-02 | 2018-03-07 | United Technologies Corporation | Dynamical system parameter identification for turbomachine |
CN108871543A (zh) * | 2018-08-30 | 2018-11-23 | 天津大学 | 恒速下叶片异步振动频率非调和傅里叶分析方法 |
CN111579060A (zh) * | 2020-06-28 | 2020-08-25 | 大连理工大学 | 一种基于叶尖定时技术高精度的旋转叶片振动测量方法 |
CN112179480A (zh) * | 2020-08-26 | 2021-01-05 | 西安交通大学 | 一种谐波叠加拟合叶尖定时振动参数识别的方法、系统及电子设备 |
CN112903271A (zh) * | 2021-01-22 | 2021-06-04 | 中国航发沈阳发动机研究所 | 一种转子叶片非接触式异步振动参数辨识方法 |
CN113586177A (zh) * | 2021-08-10 | 2021-11-02 | 西安交通大学 | 基于单叶端定时传感器的叶片固有频率识别方法 |
Non-Patent Citations (4)
Title |
---|
SPADA R P等: "Applying compressed sensing to blade tip timing data: A parametric analysis", PROCEEDINGS OF THE 10TH INTERNATIONAL CONFERENCE ON ROTOR DYNAMICS–IFTOMM, vol. 3, no. 10, 31 December 2019 (2019-12-31) * |
ZHANG J等: "High-precision extraction method for blade tip-timing signal with eddy current sensor", INTERNATIONAL JOURNAL OF ROTATING MACHINERY, 31 December 2020 (2020-12-31) * |
郭浩天;段发阶;汪猛;: "基于叶尖定时单参数法的叶片共振倍频数辨识", 天津大学学报(自然科学与工程技术版), no. 09, 30 September 2016 (2016-09-30) * |
高强等: "正弦扫频速率对结构响应的影响分析", 航天器环境工程, vol. 37, no. 3, 30 June 2020 (2020-06-30) * |
Also Published As
Publication number | Publication date |
---|---|
CN114812794B (zh) | 2024-03-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
RU2449252C2 (ru) | Способ обнаружения повреждения опорного подшипника качения двигателя | |
US9016132B2 (en) | Rotating blade analysis | |
US9341511B2 (en) | Timing analysis | |
US20170067860A1 (en) | Control Area Network Machine Diagostic | |
CN105675113B (zh) | 基于微型传感器的旋转机械角域振动信号获取装置和方法 | |
CN111397877B (zh) | 一种旋转机械拍振故障检测与诊断方法 | |
CN112105907B (zh) | 用于监测齿轮系统的方法和设备 | |
CN102967414B (zh) | 基于频谱校正的微速差双转子系统不平衡分量的提取方法 | |
CN110686764A (zh) | 基于全相位差分原理的恒速叶片异步振动频率测量方法 | |
Lin et al. | A review and strategy for the diagnosis of speed-varying machinery | |
CN105865793A (zh) | 一种提高多转子航空发动机振动监测精度的方法 | |
CN103116032A (zh) | 用于获取风力发电机组的转速的方法和设备 | |
EP1285274A1 (en) | Methof of determining speed of rotation of a motor and a computer software product to carry out the method | |
CN101368870B (zh) | 一种用于机械转子单截面轴振分析的幅值频谱的制作方法 | |
Li et al. | Application of a Method of Identifiying Instantaneous Shaft Speed from Spectrum in Aeroengine Vibration Analysis | |
CN117686232A (zh) | 一种燃气轮机振动基频实时提取方法、装置及存储介质 | |
CN114812794A (zh) | 一种扫频速率引起的叶尖定时参数辨识误差修正方法 | |
CN112465068A (zh) | 一种基于多传感器数据融合的旋转设备故障特征提取方法 | |
CN112597969A (zh) | 一种滚动轴承故障诊断方法、系统及介质 | |
CN117419923A (zh) | 适用于发动机的流水型硬件相位解算方法以及解算系统 | |
CN116317781A (zh) | 基于微速差的同轴双转子系统不平衡振动控制方法和装置 | |
CN110779723B (zh) | 一种基于霍尔信号的变速工况电机轴承精确故障诊断方法 | |
CN113358210A (zh) | 一种基于压力脉动的增压器涡轮叶片振动监测方法 | |
CN112903271A (zh) | 一种转子叶片非接触式异步振动参数辨识方法 | |
CN114674412B (zh) | 一种无转速同步信号的叶尖定时测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |