CN110457730A - 一种用于涡轴发动机外置机匣传感器布局方法 - Google Patents
一种用于涡轴发动机外置机匣传感器布局方法 Download PDFInfo
- Publication number
- CN110457730A CN110457730A CN201910444930.9A CN201910444930A CN110457730A CN 110457730 A CN110457730 A CN 110457730A CN 201910444930 A CN201910444930 A CN 201910444930A CN 110457730 A CN110457730 A CN 110457730A
- Authority
- CN
- China
- Prior art keywords
- formula
- equation
- stress
- rank
- follows
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000005452 bending Methods 0.000 claims description 16
- 230000001788 irregular Effects 0.000 claims description 8
- 238000009795 derivation Methods 0.000 claims 3
- 238000004364 calculation method Methods 0.000 claims 1
- 238000009434 installation Methods 0.000 abstract description 5
- 238000005259 measurement Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 1
- 238000005299 abrasion Methods 0.000 description 1
- 230000000712 assembly Effects 0.000 description 1
- 238000000429 assembly Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L1/00—Measuring force or stress, in general
- G01L1/16—Measuring force or stress, in general using properties of piezoelectric devices
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种用于涡轴发动机外置机匣传感器布局方法。依据欧拉—贝努力方程完成涡轴发动机外置机匣的运动方程的建立,采用数值求解法计算方程的前4阶特征根,获取各阶的固有频率,并得到与固有频率对应的振型函数,结合压电传感器测量到的信号为应力信息,得到应力与挠度的对应关系公式,然后,引入熵信息理论,在频域内考虑响应信号的各阶频率的信息量,计算各阶频率信息量沿横轴的分布,最终根据分布量获得传感器的最佳测点安装位置。
Description
技术领域
本发明属于系统测试方法应用领域,具体是一种用于涡轴发动机外置机匣传感器布局的新方法。
背景技术
作为直升机上重要的动力输出系统——涡轴发动机,是直升机上的典型机械旋转系统,是由大量相互联系但工作过程又彼此不同的部件组成,存在着强烈的非线性、非平稳性、不确定性等复杂系统特征,因此,易发生故障,且一旦发生故障,将导致灾难性事故的发生。涡轴发动机主要由转子系统、轴承系统、传动齿轮系统等部件/系统组成,其中各种组件内部靠齿轮及轴进行动力传递。如转子系统包括燃气发生器转子、自由涡轮转子等组成,其工作环境恶劣,且工作强度大,导致转子叶片易发断裂、磨损等多种故障,且状态变化规律难以掌握,性能难以准确预测,所以对涡轴发动机性能状态的监测一直是提高航空武器装备安全运行的一个无法逾越的技术瓶颈,目前最有效的技术手段是涡轴发动机上安装各种不同型号、不同规格的传感器,然而传感器的安装不是随便找个测点进行安装那么简单,且涡轴发动机内部不能安装传感器,仅有外置机匣的相应部位可以安装,再者,根据振动信号的传递规律,信号传递到外置机匣处,其衰减特别厉害,如果安装位置不精确,则有可能导致测量的信息不准确甚至测量不到相应的有用信号,无法实现涡轴发动机的准确监测目标。因此,需要大力开展涡轴发动机外置机匣传感器布局技术,提高涡轴发动机信号采集有效性及可靠性。
发明内容
本发明的目的在于提供一种用于涡轴发动机外置机匣传感器布局安装准确性的新方法,提升监测涡轴发动机实时信号的能力。
实现本发明目的的技术解决方案为:一种用于涡轴发动机外置机匣传感器布局方法,包括以下步骤:
步骤1,将涡轴发动机外置机匣视为一段弯曲不规则的横向梁,运用欧拉—贝努力方程建立弯曲横向梁的运动方程;
步骤2,运用数值求解法求取该运动方程的前4阶特征根;
步骤3,获取对应各阶的固有频率,并得到与各阶固有频率对应的振型函数;
步骤4,在发动机外部结构处任意位置贴压电传感器,并通过外部的采集设备采集压电信号,并将压电传感器测量到的信号作为为应力信息,获取应力与挠度的关系;
步骤5,运用熵信息理论,在频域内结合响应信号的各阶频率信息量,计算各阶频率信息量沿弯曲梁的横向的分布,从而确定传感器的位置布局。
本发明与现有技术相比,其显著优点:(1)本发明的方法可以快速实现在涡轴发动机上寻找最佳的安装点,从而根据最佳位置,可以实现涡轴发动机的准确监测。(2)该方法简单实用,不仅能够在发动机上,在其他的设备上,也能实现对传感器测点的快速测量与确定。
附图说明
图1是本发明的涡轴发动机外部机匣模拟钢板横向弯曲度的坐标系。
图2为涡轴发动机外置机匣传感器测点模拟试验验证原理示意图。
具体实施方式
本发明的主要实现步骤如下:
步骤1.将涡轴发动机外置机匣看做一段弯曲不规则的横向梁,运用欧拉—贝努力方程,建立弯曲横向梁振动的运动方程;
步骤2.根据数值求解法,求取运动方程的前4阶特征根;
步骤3.获取对应各阶的固有频率,并得到与各阶固有频率对应的振型函数;
步骤4.由压电传感器测量到的信号,得到应力信息与挠度的关系;
步骤5.将熵信息理论引入,在频域内考虑响应信号的各阶频率信息量,计算各阶频率信息量沿轴向的分布,根据应力信息熵大小,熵越大,说明该位置采集到的信号富含的信息越多,最终根据各阶频率信息量可以确定最佳传感器安装位置。
下面结合附图对本发明作进一步说明。
参见图1和2,一种用于涡轴发动机外置机匣传感器布局方法,由于涡轴发动机的外置机匣是不规则的弯曲梁,将其看做一段弯曲不规则的横向梁,依据欧拉—贝努力方程完成涡轴发动机外置机匣的运动方程的建立,采用数值求解法计算方程的前4阶特征根,获取各阶的固有频率,并得到与固有频率对应的振型函数,结合压电传感器测量到的信号为应力信息,得到应力与挠度的对应关系公式,然后,引入熵信息理论,在频域内考虑响应信号的各阶频率的信息量,计算各阶频率信息量沿横轴的分布,最终根据分布量,可以获得传感器的最佳测点安装位置。
步骤1,将涡轴发动机外置机匣看做一段弯曲不规则的横向梁,根据欧拉—贝努力方程,建立弯曲梁振动的运动方程。各参数及方程建立如下:
假设涡轴发动机的外置机匣为一段弯曲不规则的横向梁,梁长为l,单位长度的质量为ρ,抗弯刚度EI,建立如图1所示的坐标系,坐标x处微元dx在任一时刻的横向力为p(x,t),弯曲梁振动的运动方程为
当p(x,t)≡0时:
y(x,t)=Y(x)Y(t) (2)
振型函数为:
Y(x)=c1chβx+c2shβx+c3cosβx+c4sinβx (3)
欧拉—贝努力方程完成弯曲横向梁的运动方程,其主要是为了实现将弯曲横向梁的应力分布通过数学公式表达出来,从而根据应力分布最终找出传感器的最佳布局。
步骤2,主要采用数值法求得前4阶特征根:
β1l=1.875’ β2l=4.694
β3l=7.855’ β4l=10.966
运用数值求解法完成应力方程的前4阶特征根的计算,其主要是为了得到各阶的固有频率。
步骤3,根据步骤2得到的特征根,获取对应各阶的固有频率,并得到与各阶固有频率对应的振型函数。
对应的各阶固有频率为:
与ωi对应的振型函数为:
Yi=(x)=chβix-cosβix+ki(shβix-sinβix) (5)
式中:
运用熵信息,主要是为了计算各阶频率信息量沿弯曲梁的横向分布,根据分布量可以确定传感器的最佳布局。
步骤4,由压电传感器测量得到的信号为应力信息,应力与挠度的关系为:
式中:tb为梁的厚度。
步骤5,引用信息熵的概念,在频域内考虑响应信号的各阶频率信息量,计算各阶频率信息量沿x轴的分布。根据式(2),上式可以简化为:
式中:k为常量。
根据式(13),应力信息的概率密度为:
根据式(6)应力信息熵为:
指的是n阶横向力,此时计算出的E(x)即为应力信息熵,其为一组在弯曲梁长度内的变化曲线,找出其峰值点,即为所需要找出的传感器最佳布点。
首先取一块长度为1米的弯曲钢板,由于知道梁长为l,单位长度的质量为ρ,抗弯刚度EI,,则根据上述公式可以计算出应力信息熵。然后,搭建一套以压电传感器为测试信号的应力测试系统,验证系统如图2所示,在计算好的响应位置点安装压电传感器,然后根据压电传感器,完成振动台随机激励对应的信息熵,与计算出来的进行对比,以确定传感器布局方法的准确性。
Claims (5)
1.一种用于涡轴发动机外置机匣传感器布局方法,其特征在于包括以下步骤:
步骤1,将涡轴发动机外置机匣视为一段弯曲不规则的横向梁,运用欧拉—贝努力方程建立弯曲横向梁的运动方程;
步骤2,运用数值求解法求取该运动方程的前4阶特征根;
步骤3,获取对应各阶的固有频率,并得到与各阶固有频率对应的振型函数;
步骤4,在发动机外部结构处任意位置安装压电传感器,通过外部的采集设备采集压电信号,并将压电传感器测量到的信号作为应力信息,获取应力与挠度的关系;
步骤5,运用熵信息理论,在频域内结合响应信号的各阶频率信息量,计算各阶频率信息量沿弯曲梁的横向的分布,从而确定传感器的位置布局。
2.根据权利要求1所述的传感器布局方法,其特征在于,步骤1运用欧拉—贝努力方程建立弯曲横向梁的运动方程的具体方法如下:
设涡轴发动机的外置机匣为一段弯曲不规则的横向梁,梁长为l,单位长度的质量为ρ,抗弯刚度EI,建立单位长度内横向力坐标系,坐标x处微元dx在任一时刻的横向力为p(x,t),弯曲梁振动的运动方程为
式中,x和y代表的是单位长和对应的横向力;
当p(x,t)≡0时:
y(x,t)=Y(x)Y(t) (2)
所构建的振型函数为:
Y(x)=c1chβx+c2shβx+c3cosβx+c4sinβx (3)
式中,c代表建立振型函数的系数,β是特征根。
3.根据权利要求1所述的传感器布局方法,其特征在于:所述步骤3中获取对应的各阶固有频率为:
式中,EI为抗弯刚度均值;βi为四阶特征根,此时的i取值为1~4;
得到与ωi对应的振型函数为:
Yi(x)=chβix-cosβix+ki(shβix-sinβix) (5)
式中:
4.根据权利要求1所述的传感器布局方法,其特征在于:步骤4中,应力与挠度的关系为:
式中:tb为梁的厚度,y″(x,t)为应力的二阶求导。
5.根据权利要求1、2、3或4中任意一项所述的传感器布局方法,其特征在于:所述步骤5中,引用信息熵的概念,在频域内结合响应信号的各阶频率信息量,计算各阶频率信息量沿x轴的分布;
根据式(2),式(6)简化为:
σ(x)=k*Y″(x) (7)
式中:k为常量;
根据应力信息的概率密度计算公式及式(7),应力信息的概率密度为:
式中,Yi″指的是横向二阶求导横向力二阶求导;
根据式(6)应力信息熵为:
式中,指的是n阶横向力。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910444930.9A CN110457730A (zh) | 2019-05-27 | 2019-05-27 | 一种用于涡轴发动机外置机匣传感器布局方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910444930.9A CN110457730A (zh) | 2019-05-27 | 2019-05-27 | 一种用于涡轴发动机外置机匣传感器布局方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110457730A true CN110457730A (zh) | 2019-11-15 |
Family
ID=68481049
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910444930.9A Pending CN110457730A (zh) | 2019-05-27 | 2019-05-27 | 一种用于涡轴发动机外置机匣传感器布局方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110457730A (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070087756A1 (en) * | 2005-10-04 | 2007-04-19 | Hoffberg Steven M | Multifactorial optimization system and method |
CN104374468A (zh) * | 2014-10-24 | 2015-02-25 | 哈尔滨工业大学深圳研究生院 | 一种细长物体固有频率测量方法 |
EP3236384A1 (en) * | 2008-11-21 | 2017-10-25 | Exxonmobil Upstream Research Company | Methods and systems for modeling, designing, and conducting drilling operations that consider vibrations |
-
2019
- 2019-05-27 CN CN201910444930.9A patent/CN110457730A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070087756A1 (en) * | 2005-10-04 | 2007-04-19 | Hoffberg Steven M | Multifactorial optimization system and method |
EP3236384A1 (en) * | 2008-11-21 | 2017-10-25 | Exxonmobil Upstream Research Company | Methods and systems for modeling, designing, and conducting drilling operations that consider vibrations |
CN104374468A (zh) * | 2014-10-24 | 2015-02-25 | 哈尔滨工业大学深圳研究生院 | 一种细长物体固有频率测量方法 |
Non-Patent Citations (1)
Title |
---|
魏民祥: "基于信息熵原理的结构振动传感器的布局", 《仪器仪表学报》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8951011B2 (en) | Wind turbine and a method for monitoring a wind turbine | |
CN101016879B (zh) | 用于机舱风速修正的方法和系统 | |
US8332164B2 (en) | Method for determining fatigue damage in a power train of a wind turbine | |
EP2868921B1 (en) | Wind turbine and method for evaluating health state of blade thereof | |
US8511182B2 (en) | Determining an accumulated load of a wind turbine in angular sectors | |
CN102622458B (zh) | 一种风力发电机组振动与载荷综合评估系统及评估方法 | |
CN103105296A (zh) | 轴承动态特性参数测试装置 | |
CA1329847C (en) | Turbine blade fatigue monitor | |
CN106052743B (zh) | 一种评估传感器质量对频响函数影响大小的方法 | |
EP2369292A1 (en) | A Strain Gauge Assembly and Method | |
WO2013129989A1 (en) | Condition monitoring of a rotating system based on a time stamped signal | |
CN108225785B (zh) | 具有传感器布置的装置模块 | |
US20180363633A1 (en) | Abnormality monitoring apparatus and abnormality monitoring method for wind farm | |
JP2016188612A (ja) | 疲労評価システム及びこれを備えた風力発電装置、並びに、風力発電装置の疲労評価方法 | |
Rebelo et al. | Structural monitoring of a wind turbine steel tower—Part I: System description and calibration | |
CN104462862A (zh) | 一种基于三次b样条尺度函数的机械结构动载荷识别方法 | |
CN108363832A (zh) | 一种计算输电线路覆冰厚度的方法 | |
CN103245913A (zh) | 大型发电机组次同步振荡信号检测与分析的方法及系统 | |
CN107014486B (zh) | 一种输电线路微风振动监测装置的核查单元及方法 | |
JP2018031610A (ja) | 累積疲労損傷度推定システムまたは累積疲労損傷度推定方法 | |
CN110457730A (zh) | 一种用于涡轴发动机外置机匣传感器布局方法 | |
Antunes | Monitoring sea level change at Cascais tide gauge | |
CN108204342A (zh) | 风力发电机的叶片结冰识别方法和装置 | |
CN113495169B (zh) | 风力机风轮前风速超声测量系统 | |
CN110220688A (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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20191115 |