CN111595433B - 航空发动机整机振动传感器的位置确定方法及系统 - Google Patents
航空发动机整机振动传感器的位置确定方法及系统 Download PDFInfo
- Publication number
- CN111595433B CN111595433B CN201910124612.4A CN201910124612A CN111595433B CN 111595433 B CN111595433 B CN 111595433B CN 201910124612 A CN201910124612 A CN 201910124612A CN 111595433 B CN111595433 B CN 111595433B
- Authority
- CN
- China
- Prior art keywords
- nodes
- node
- node set
- vibration
- vibration sensor
- 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.)
- Active
Links
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
-
- 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
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本公开涉及一种航空发动机整机振动传感器的位置确定方法及系统。方法包括:建立航空发动机的整机静子机匣的有限元模型,并设定整机振动传感器的备选节点集;对所述有限元模型进行模态分析,并对所述有限元模型进行谐响应分析,获得所述备选节点集中各个节点的参数响应曲线;在所述备选节点集中去除属于局部模态的节点;对所述备选节点集中各个节点的参数响应曲线进行统计,并根据得到的统计数据建立目标优化函数;计算所述备选节点集中各个节点对应的目标优化函数值,并按照所述目标优化函数值的数值大小对节点进行排序;根据排序结果确定所述整机振动传感器的设置位置。本公开实施例能够简化整机振动传感器的设置位置的确定过程。
Description
技术领域
本公开涉及航空发动机测试领域,尤其涉及一种航空发动机整机振动传感器的位置确定方法及系统。
背景技术
航空发动机整机振动(engine vibration)是航空发动机整机层面的振动,主要是指由发动机转子引起的振动。航空发动机整机振动测量是状态检测和故障诊断的重要手段,是航空发动机研制、试车、大修、服役等各个环节中不可或缺的重要内容。
振动传感器是振动测量的关键组件,将振动量(例如振动位移、振动速度、振动加速度)的变化转换成电量(例如电压、电流、电荷)或电参量(例如电阻、电感、电容)的变化。整机振动测量中的一个任务是选择振动传感器的测振点位置。转子不平衡力是整机振动最主要的激振源,振动传感器的测点要能监测到转子的不平衡振动,振动传感器测量的振动量应和转子不平衡力有明显和稳定的关系,能正确地反映发动机转子的振动。原则上说转子各支点轴承座位置最能反映转子的振动水平,但是轴承座位置内部空间狭小且环境恶劣,限制了传感器的安装。而整机振动测点位置选取合理与否,会影响测试结果的准确性,影响状态检测和故障诊断的有效应和正确性。
目前,整机振动传感器的具体设置位置的确定或基于经验或基于定性分析或基于试验逐步逼近,精度和可靠性不高,时间和经济成本高。
发明内容
有鉴于此,本公开实施例提供一种航空发动机整机振动传感器的位置确定方法及系统,能够简化整机振动传感器的设置位置的确定过程。
在本公开的一个方面,提供一种航空发动机整机振动传感器的位置确定方法,包括:
建立航空发动机的整机静子机匣的有限元模型,并设定整机振动传感器的备选节点集;
对所述有限元模型进行模态分析,并对所述有限元模型进行谐响应分析,获得所述备选节点集中各个节点的参数响应曲线;
在所述备选节点集中去除属于局部模态的节点;
对所述备选节点集中各个节点的参数响应曲线进行统计,并根据得到的统计数据建立目标优化函数;
计算所述备选节点集中各个节点对应的目标优化函数值,并按照所述目标优化函数值的数值大小对节点进行排序;
根据排序结果确定所述整机振动传感器的设置位置。
在一些实施例中,对所述备选节点集中各个节点的参数响应曲线进行统计的操作包括:
计算所述备选节点集中各个节点的振动传递率的均方根;
计算所述备选节点集中各个节点的振动传递率的标准差;和
计算所述备选节点集中各个节点与轴承座中心节点的振动传递率的相关系数。
在一些实施例中,建立目标优化函数的操作包括:
对所述均方根和所述标准差分别进行归一化,得到归一化均方根和归一化标准差;
设定所述归一化均方根、所述归一化标准差和所述相关系数的权值,并通过加权构造单目标优化函数。
在一些实施例中,建立目标优化函数的操作包括:
根据所述均方根、所述标准差和所述相关系数构造多目标优化函数。
在一些实施例中,所述参数响应曲线为振动传递率响应曲线。
在一些实施例中,在获得参数响应曲线后,还包括:
在所述备选节点集中去除整个频率段存在振动传递率小于第一阈值的节点;和/或
在所述备选节点集中去除整个频率段存在振动传递率大于第二阈值的节点,
其中,所述第一阈值小于所述第二阈值。
在一些实施例中,根据排序结果确定所述整机振动传感器的设置位置的操作包括:
根据所述整机振动传感器的需求数量,在所述排序结果中依序选择待选节点,并剔除所述待选节点中安装受限的待选节点;
将剔除后剩余的待选节点作为所述整机振动传感器安装的优选节点;
输出所述优选节点对应的节点编号和坐标。
在一些实施例中,根据排序结果确定所述整机振动传感器的设置位置的操作包括:
根据所述整机振动传感器的需求数量,在所述排序结果中依序选择待选节点,输出所述待选节点对应的节点编号和坐标;
剔除所述待选节点中安装受限的待选节点;
将剔除后剩余的待选节点作为所述整机振动传感器安装的优选节点。
在一些实施例中,设定整机振动传感器的备选节点集的操作包括:
设定所述整机振动传感器的备选节点集包括转子各支点轴承座的外表面节点和所述整机静子的承力机匣的安装边表面节点。
在本公开的一个方面,提供一种航空发动机整机振动传感器的位置确定系统,包括:
模型建立模块,用于建立航空发动机的整机静子机匣的有限元模型,并设定整机振动传感器的备选节点集;
模态分析模块,用于对所述有限元模型进行模态分析;
曲线获得模块,用于对所述有限元模型进行谐响应分析,获得所述备选节点集中各个节点的参数响应曲线;
节点去除模块,用于在所述备选节点集中去除属于局部模态的节点;
函数建立模块,用于对所述备选节点集中各个节点的参数响应曲线进行统计,并根据得到的统计数据建立目标优化函数;
函数计算模块,用于计算所述备选节点集中各个节点对应的目标优化函数值;
位置确定模块,用于按照所述目标优化函数值的数值大小对节点进行排序,并根据排序结果确定所述整机振动传感器的设置位置。
在一些实施例中,所述函数建立模块包括:
统计单元,用于计算所述备选节点集中各个节点的振动传递率的均方根和标准差,并计算所述备选节点集中各个节点与轴承座中心节点的振动传递率的相关系数;
归一化单元,用于对所述均方根和所述标准差分别进行归一化,得到归一化均方根和归一化标准差;
第一函数构造单元,用于设定所述归一化均方根、所述归一化标准差和所述相关系数的权值,并通过加权构造单目标优化函数。
在一些实施例中,所述函数建立模块包括:
统计单元,用于计算所述备选节点集中各个节点的振动传递率的均方根和标准差,并计算所述备选节点集中各个节点与轴承座中心节点的振动传递率的相关系数;
第二函数构造单元,用于根据所述均方根、所述标准差和所述相关系数构造多目标优化函数。
在一些实施例中,所述系统还包括:
节点调整模块,用于在所述曲线获得模块获得参数响应曲线后,调整所述备选节点集内的节点。
在一些实施例中,所述参数响应曲线为振动传递率响应曲线,所述节点调整模块用于在所述曲线获得模块获得参数响应曲线后,在所述备选节点集中去除整个频率段存在振动传递率小于第一阈值的节点,和/或整个频率段存在振动传递率大于第二阈值的节点,所述第一阈值小于所述第二阈值。
因此,根据本公开实施例,通过有限元模型对航空发动机的整机静子机匣进行模态分析和谐响应分析,并获得整机振动传感器的备选节点集中各节点的参数响应曲线,通过删除备选节点集内属于局部模态的节点来调整备选节点集内的节点,再通过统计参数响应曲线来建立目标优化函数,再通过函数值的排序来确定整机振动传感器的设置位置。这样可使确定的整机振动传感器的位置可量化,并且精度和可靠性较高。
附图说明
构成说明书的一部分的附图描述了本公开的实施例,并且连同说明书一起用于解释本公开的原理。
参照附图,根据下面的详细描述,可以更加清楚地理解本公开,其中:
图1是根据本公开航空发动机整机振动传感器的位置确定方法的一些实施例的流程示意图;
图2是根据本公开航空发动机整机振动传感器的位置确定方法的一些实施例中构造单目标优化函数的流程示意图;
图3是根据本公开航空发动机整机振动传感器的位置确定方法的一些实施例中构造多目标优化函数的流程示意图;
图4是根据本公开航空发动机整机振动传感器的位置确定系统的一些实施例的方框示意图。
应当明白,附图中所示出的各个部分的尺寸并不是按照实际的比例关系绘制的。此外,相同或类似的参考标号表示相同或类似的构件。
具体实施方式
现在将参照附图来详细描述本公开的各种示例性实施例。对示例性实施例的描述仅仅是说明性的,决不作为对本公开及其应用或使用的任何限制。本公开可以以许多不同的形式实现,不限于这里所述的实施例。提供这些实施例是为了使本公开透彻且完整,并且向本领域技术人员充分表达本公开的范围。应注意到:除非另外具体说明,否则在这些实施例中阐述的部件和步骤的相对布置、材料的组分、数字表达式和数值应被解释为仅仅是示例性的,而不是作为限制。
本公开中使用的“第一”、“第二”以及类似的词语并不表示任何顺序、数量或者重要性,而只是用来区分不同的部分。“包括”或者“包含”等类似的词语意指在该词前的要素涵盖在该词后列举的要素,并不排除也涵盖其他要素的可能。“上”、“下”、“左”、“右”等仅用于表示相对位置关系,当被描述对象的绝对位置改变后,则该相对位置关系也可能相应地改变。
在本公开中,当描述到特定器件位于第一器件和第二器件之间时,在该特定器件与第一器件或第二器件之间可以存在居间器件,也可以不存在居间器件。当描述到特定器件连接其它器件时,该特定器件可以与所述其它器件直接连接而不具有居间器件,也可以不与所述其它器件直接连接而具有居间器件。
本公开使用的所有术语(包括技术术语或者科学术语)与本公开所属领域的普通技术人员理解的含义相同,除非另外特别定义。还应当理解,在诸如通用字典中定义的术语应当被解释为具有与它们在相关技术的上下文中的含义相一致的含义,而不应用理想化或极度形式化的意义来解释,除非这里明确地这样定义。
对于相关领域普通技术人员已知的技术、方法和设备可能不作详细讨论,但在适当情况下,所述技术、方法和设备应当被视为说明书的一部分。
如图1所示,是根据本公开航空发动机整机振动传感器的位置确定方法的一些实施例的流程示意图。参考图1,在一些实施例中,航空发动机整机振动传感器的位置确定方法,包括步骤100-步骤600。在步骤100中,建立航空发动机的整机静子机匣的有限元模型,并设定整机振动传感器的备选节点集。由于转子支点轴承座附近和靠近轴承座截面的承力机匣安装边刚性较好,比较适宜安装,因此在设定备选节点集时,可使备选节点集包括转子各支点轴承座的外表面节点和所述整机静子的承力机匣的安装边表面节点。相比于更大范围的设置节点集,根据经验或者历史数据可选择适宜安装的上述部位的表面节点作为节点集可减少确定传感器位置时的计算量。当然,在另一些实施例中,也可以选择其他部位的节点或者选择更大范围的节点形成备选节点集。
在步骤200中,对所述有限元模型进行模态分析,并对所述有限元模型进行谐响应分析,获得所述备选节点集中各个节点的参数响应曲线。在模态分析时,可采用例如模态有效质量因子或模态参与因子比值等确定整体模态和局部模态,确定方法可采用将模态质量因子或模态参与因子比值发生突变的值设定为阈值,再通过该阈值区分整体模态和局部模态,例如将超过该阈值的模态定为整体模态,未超过的定为局部模态。除此之外,还可以采用其他方法设定上述阈值。
在完成模态分析后,可在整机静子机匣各个轴承座施加简谐载荷,以便实现谐响应分析,以获得备选节点集中各个节点的参数响应曲线。谐响应分析可采用模态叠加法或完全法。在已进行模态分析的基础上,优选通过模态叠加法对所述有限元模型进行谐响应分析。相比于其他的谐响应分析方法,模态叠加法计算效率更高。
通过谐响应分析获得的参数响应曲线的横轴为频率,纵轴为参数的幅值。例如,选择振动传递率(即测点振动幅值与激振力的比值)作为优化的基础物理量,对应得到的参数响应曲线为振动传递率响应曲线。该参数仅与静子机匣的动力学特性相关,而与激励大小无关,从而使通过该参数进行的整机振动传感器的位置优化更加准确。
在步骤300中,在所述备选节点集中去除属于局部模态的节点,以消除整体振动传感器对应的节点受到局部模态的影响。在另一些实施例中,还可以在步骤200之后,在步骤300之前、之后或同步地调整备选节点集内的节点,以免这些节电对后续的优化函数构造形成不利影响。例如在备选节点集中去除整个频率段存在振动传递率小于第一阈值的节点,以避免节点的某些频段振动信号过小而影响目标优化函数,和/或去除整个频率段存在振动传递率大于第二阈值的节点,该第二阈值大于第一阈值,以避免节点的某些频段振动过高而影响目标优化函数。
相应的,步骤400-步骤600中对备选节点集中各个节点的操作即为对去除某些节点后的备选节点集中所有节点的操作。
在步骤400中,对所述备选节点集中各个节点的参数响应曲线进行统计,并根据得到的统计数据建立目标优化函数。在获得参数响应曲线之后,可根据整机振动测点优化原则选择统计各类数据,并确定适合的目标优化函数。在步骤500中,计算所述备选节点集中各个节点对应的目标优化函数值,并按照所述目标优化函数值的数值大小对节点进行排序。再在步骤600中,根据排序结果确定所述整机振动传感器的设置位置。根据整机振动传感器布置的需求数量,可选择排序靠前的多个节点确定各个整机震动传感器的设置位置。
考虑到有些排序靠前的节点受到自身或周围环境的影响而无法正常安装,因此可在依序选择待选节点时,剔除待选节点中安装受限的待选节点,即选择节点时跳过安装受限的节点。然后,将剔除后剩余的待选节点作为所述整机振动传感器安装的优选节点,再输出所述优选节点对应的节点编号和坐标。
在另一些实施例中,也可先根据所述整机振动传感器的需求数量,在所述排序结果中依序选择待选节点,输出所述待选节点对应的节点编号和坐标。此时可以输出多于需求数量的节点,以方便人工挑选,在这些节点中可剔除所述待选节点中安装受限的待选节点,再将剔除后剩余的待选节点作为所述整机振动传感器安装的优选节点。
通过上述实施例的说明可知,相比于根据经验、或者定性分析、或者通过多次试验逐步逼近的传统方式,本实施例通过建立有限元模型来形成目标优化函数的方式,将整机静子机匣上可以设置整机振动传感器的各个位置作为备选节点进行函数值的计算,并根据计算结果排序来挑选适合的传感器安装位置,一方面简化了传感器安装位置的确定过程,节省了试验成本和时间,降低对经验的依赖,另一方面也使得位置的确定相比于定性分析更加精准和可靠。
在上述实施例中,统计量的选择对整机振动传感器的位置确定比较重要。在一些实施例中,选择备选节点集中各个节点的振动传递率的均方根(Root Mean Square,简称RMS)、标准差(Standard Deviation,简称STD)以及节点与轴承座中心节点的相关系数(correlation coefficient,简称CORRCOEF)来参与构造目标优化函数。
在这其中,节点的振动传递率的均方根的计算方式为将该节点中各个频率下的振动传递率平方后加起来再开方,可用于衡量振动传递率的总量。均方根越大说明振动传递率越高,该节点对输入的振动信号的衰减越小。节点的振动传递率的标准差的计算方式为将该节点中各个频率下的振动传递率与振动传递率的平均值的差平方后加起来再开方,可用于衡量振动传递率的离散性。标准差越小说明各频率段的振动传递率比较一致,振动传递率比较稳定。节点与轴承座中心节点的振动传递率的相关系数为该节点和轴承座中心节点的协方差与该节点和轴承座中心节点的方差开方乘积的比值,可用于衡量节点与轴承座中心节点的振动传递率之间的相关性。相关系数越高,该节点的振动特性与轴承座中心节点的振动特性一致性越高。
如图2所示,为根据本公开航空发动机整机振动传感器的位置确定方法的一些实施例中构造单目标优化函数的流程示意图。相比于前述各实施例,本实施例的步骤400具体包括步骤410-步骤430。在步骤410中,计算备选节点集中各个节点的振动传递率的均方根和标准差,并计算备选节点集中各个节点与轴承座中心节点的的振动传递率的相关系数。在步骤420中,对均方根和标准差分别进行归一化,得到归一化均方根(NRMS)和归一化标准差(NSTD)。再在步骤430中,设定归一化均方根、归一化标准差和相关系数的权值,并通过加权构造单目标优化函数。在归一化处理振动传递率的均方根时,可采用节点振动传递率的均方根除以整个节点集中振动传递率均方根的最大值的方式。在归一化处理振动传递率的标准差时,可采用节点振动传递率的标准差除以整个节点集中振动传递率标准差的最大值的方式。通过设定权值,可以获得对应的加权公式,例如:
F(x)=a*NRMS+b*NSTD+c*CORRCOEF,
其中,a、b、c分别为归一化均方根NRMS、归一化标准差NSTD和相关系数CORRCOEF的权值。权值的设定可根据历史数据或者多次试验确定。
如图3所示,是根据本公开航空发动机整机振动传感器的位置确定方法的一些实施例中构造多目标优化函数的流程示意图。与上一实施例相比,本实施例在步骤410之后包括步骤440,即根据均方根、标准差和相关系数构造多目标优化函数。例如构造均方根、标准差和相关系数每个统计量各自的目标优化函数。而对于多目标优化函数,可采用例如进化算法(Evolutionary Algorithm,简称EA)、粒子群算法(Particle Swarm Optimization,简称PSO)等进行求解。
如图4所示,是根据本公开航空发动机整机振动传感器的位置确定系统的一些实施例的方框示意图。参考图4,在一些实施例中,航空发动机整机振动传感器的位置确定系统可包括:模型建立模块10、模态分析模块20、曲线获得模块30、节点去除模块40、函数建立模块50、函数计算模块60和位置确定模块70。模型建立模块10用于建立航空发动机的整机静子机匣的有限元模型,并设定整机振动传感器的备选节点集。模态分析模块20用于对所述有限元模型进行模态分析。曲线获得模块30用于对所述有限元模型进行谐响应分析,获得所述备选节点集中各个节点的参数响应曲线。节点去除模块40用于在所述备选节点集中去除属于局部模态的节点
函数建立模块50用于对所述备选节点集中各个节点的参数响应曲线进行统计,并根据得到的统计数据建立目标优化函数。在一些实施例中,函数建立模块50可包括:统计单元、归一化单元和第一函数构造单元。统计单元用于计算所述备选节点集中各个节点的振动传递率的均方根和标准差,并计算所述备选节点集中各个节点与轴承座中心节点的振动传递率的相关系数。归一化单元用于对所述均方根和所述标准差分别进行归一化,得到归一化均方根和归一化标准差。第一函数构造单元用于设定所述归一化均方根、所述归一化标准差和所述相关系数的权值,并通过加权构造单目标优化函数。
在另一些实施例中,函数建立模块50可包括:统计单元和第二函数构造单元。统计单元,用于计算所述备选节点集中各个节点的振动传递率的均方根和标准差,并计算所述备选节点集中各个节点与轴承座中心节点的振动传递率的相关系数。第二函数构造单元,用于根据所述均方根、所述标准差和所述相关系数构造多目标优化函数。
函数计算模块60用于计算所述备选节点集中各个节点对应的目标优化函数值。位置确定模块70用于按照所述目标优化函数值的数值大小对节点进行排序,并根据排序结果确定所述整机振动传感器的设置位置。
在一些实施例中,上述系统实施例可进一步包括节点调整模块,该模块用于在所述曲线获得模块30获得参数响应曲线后,调整所述备选节点集内的节点。对于参数响应曲线为振动传递率响应曲线的情况,所述节点调整模块用于在所述曲线获得模块30获得参数响应曲线后,在所述备选节点集中去除整个频率段存在振动传递率小于第一阈值的节点,和/或整个频率段存在振动传递率大于第二阈值的节点。
本说明书中多个实施例采用递进的方式描述,各实施例的重点有所不同,而各个实施例之间相同或相似的部分相互参见即可。对于系统实施例而言,由于其整体以及涉及的模块、单元等与方法实施例中的内容存在对应关系,因此描述的比较简单,相关之处参见方法实施例的部分说明即可。
至此,已经详细描述了本公开的各实施例。为了避免遮蔽本公开的构思,没有描述本领域所公知的一些细节。本领域技术人员根据上面的描述,完全可以明白如何实施这里公开的技术方案。
虽然已经通过示例对本公开的一些特定实施例进行了详细说明,但是本领域的技术人员应该理解,以上示例仅是为了进行说明,而不是为了限制本公开的范围。本领域的技术人员应该理解,可在不脱离本公开的范围和精神的情况下,对以上实施例进行修改或者对部分技术特征进行等同替换。本公开的范围由所附权利要求来限定。
Claims (13)
1.一种航空发动机整机振动传感器的位置确定方法,包括:
建立航空发动机的整机静子机匣的有限元模型,并设定整机振动传感器的备选节点集;
对所述有限元模型进行模态分析,并对所述有限元模型进行谐响应分析,获得所述备选节点集中各个节点的参数响应曲线,所述参数响应曲线为振动传递率响应曲线;
在所述备选节点集中去除属于局部模态的节点;
对所述备选节点集中各个节点的参数响应曲线进行统计,并根据得到的统计数据建立目标优化函数;
计算所述备选节点集中各个节点对应的目标优化函数值,并按照所述目标优化函数值的数值大小对节点进行排序;
根据排序结果确定所述整机振动传感器的设置位置;
其中,对所述备选节点集中各个节点的参数响应曲线进行统计的操作包括:
计算所述备选节点集中各个节点的振动传递率的均方根;
计算所述备选节点集中各个节点的振动传递率的标准差;和
计算所述备选节点集中各个节点与轴承座中心节点的振动传递率的相关系数。
2.根据权利要求1所述的方法,其中,建立目标优化函数的操作包括:
对所述均方根和所述标准差分别进行归一化,得到归一化均方根和归一化标准差;
设定所述归一化均方根、所述归一化标准差和所述相关系数的权值,并通过加权构造单目标优化函数。
3.根据权利要求1所述的方法,其中,建立目标优化函数的操作包括:
根据所述均方根、所述标准差和所述相关系数构造多目标优化函数。
4.根据权利要求1所述的方法,其中,所述参数响应曲线为振动传递率响应曲线。
5.根据权利要求4所述的方法,其中,在获得参数响应曲线后,还包括:
在所述备选节点集中去除整个频率段存在振动传递率小于第一阈值的节点;和/或
在所述备选节点集中去除整个频率段存在振动传递率大于第二阈值的节点,
其中,所述第一阈值小于所述第二阈值。
6.根据权利要求1所述的方法,其中,根据排序结果确定所述整机振动传感器的设置位置的操作包括:
根据所述整机振动传感器的需求数量,在所述排序结果中依序选择待选节点,并剔除所述待选节点中安装受限的待选节点;
将剔除后剩余的待选节点作为所述整机振动传感器安装的优选节点;
输出所述优选节点对应的节点编号和坐标。
7.根据权利要求1所述的方法,其中,根据排序结果确定所述整机振动传感器的设置位置的操作包括:
根据所述整机振动传感器的需求数量,在所述排序结果中依序选择待选节点,输出所述待选节点对应的节点编号和坐标;
剔除所述待选节点中安装受限的待选节点;
将剔除后剩余的待选节点作为所述整机振动传感器安装的优选节点。
8.根据权利要求1所述的方法,其中,设定整机振动传感器的备选节点集的操作包括:
设定所述整机振动传感器的备选节点集包括转子各支点轴承座的外表面节点和所述整机静子的承力机匣的安装边表面节点。
9.一种航空发动机整机振动传感器的位置确定系统,包括:
模型建立模块,用于建立航空发动机的整机静子机匣的有限元模型,并设定整机振动传感器的备选节点集;
模态分析模块,用于对所述有限元模型进行模态分析;
曲线获得模块,用于对所述有限元模型进行谐响应分析,获得所述备选节点集中各个节点的参数响应曲线,所述参数响应曲线为振动传递率响应曲线;
节点去除模块,用于在所述备选节点集中去除属于局部模态的节点;
函数建立模块,用于对所述备选节点集中各个节点的参数响应曲线进行统计,并根据得到的统计数据建立目标优化函数;
函数计算模块,用于计算所述备选节点集中各个节点对应的目标优化函数值;
位置确定模块,用于按照所述目标优化函数值的数值大小对节点进行排序,并根据排序结果确定所述整机振动传感器的设置位置;
其中,所述函数建立模块包括:
统计单元,用于计算所述备选节点集中各个节点的振动传递率的均方根和标准差,并计算所述备选节点集中各个节点与轴承座中心节点的振动传递率的相关系数。
10.根据权利要求9所述的系统,其中,所述函数建立模块还包括:
归一化单元,用于对所述均方根和所述标准差分别进行归一化,得到归一化均方根和归一化标准差;
第一函数构造单元,用于设定所述归一化均方根、所述归一化标准差和所述相关系数的权值,并通过加权构造单目标优化函数。
11.根据权利要求9所述的系统,其中,所述函数建立模块还包括:
第二函数构造单元,用于根据所述均方根、所述标准差和所述相关系数构造多目标优化函数。
12.根据权利要求9所述的系统,还包括:
节点调整模块,用于在所述曲线获得模块获得参数响应曲线后,调整所述备选节点集内的节点。
13.根据权利要求12所述的系统,其中,所述参数响应曲线为振动传递率响应曲线,所述节点调整模块用于在所述曲线获得模块获得参数响应曲线后,在所述备选节点集中去除整个频率段存在振动传递率小于第一阈值的节点,和/或整个频率段存在振动传递率大于第二阈值的节点,所述第一阈值小于所述第二阈值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910124612.4A CN111595433B (zh) | 2019-02-20 | 2019-02-20 | 航空发动机整机振动传感器的位置确定方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910124612.4A CN111595433B (zh) | 2019-02-20 | 2019-02-20 | 航空发动机整机振动传感器的位置确定方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111595433A CN111595433A (zh) | 2020-08-28 |
CN111595433B true CN111595433B (zh) | 2022-07-08 |
Family
ID=72188818
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910124612.4A Active CN111595433B (zh) | 2019-02-20 | 2019-02-20 | 航空发动机整机振动传感器的位置确定方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111595433B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116132882B (zh) * | 2022-12-22 | 2024-03-19 | 苏州上声电子股份有限公司 | 一种扬声器的安装位置确定方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101799844A (zh) * | 2010-04-08 | 2010-08-11 | 北京航空航天大学 | 基于故障频率敏感的直升机齿轮箱振动传感器布置方法 |
CN104991982A (zh) * | 2015-04-16 | 2015-10-21 | 北京理工大学 | 一种飞行器气动弹性惯性传感器布局方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10360326B2 (en) * | 2016-05-18 | 2019-07-23 | Solar Turbines Incorporated | Method for determining vibratory contact stress at a blade attachment |
-
2019
- 2019-02-20 CN CN201910124612.4A patent/CN111595433B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101799844A (zh) * | 2010-04-08 | 2010-08-11 | 北京航空航天大学 | 基于故障频率敏感的直升机齿轮箱振动传感器布置方法 |
CN104991982A (zh) * | 2015-04-16 | 2015-10-21 | 北京理工大学 | 一种飞行器气动弹性惯性传感器布局方法 |
Non-Patent Citations (3)
Title |
---|
基于有限元仿真的航空发动机振动传感器布局优化方法研究;欧阳运芳等;《汽轮机技术》;20181025(第05期);第359-362页 * |
基于有限元分析的变压器振动测点优化;李阳海等;《广东电力》;20171025(第10期);第121-125页 * |
齿轮箱传感器优化布置研究;魏秀业等;《兵工学报》;20101115(第11期);第1508-1513页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111595433A (zh) | 2020-08-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109949290B (zh) | 路面裂缝检测方法、装置、设备及存储介质 | |
CN109933852B (zh) | 预测车辆尺寸偏差的方法、装置、存储介质及电子设备 | |
CN101504333B (zh) | 基于应变监测的索结构中索系统的递进式健康监测方法 | |
CN115166618B (zh) | 一种非平稳输出的电流互感器误差评估方法 | |
CN107885928A (zh) | 考虑测量误差的步进应力加速性能退化可靠性分析方法 | |
CN112199793B (zh) | 一种非高斯载荷的损伤等效试验谱优化方法 | |
CN105488352B (zh) | 基于长期挠度监测数据的混凝土桥梁刚度可靠度评估方法 | |
CN113313145B (zh) | 一种基于混合核相关向量机的快速道路交通事件检测方法 | |
CN110333078B (zh) | 一种滚动轴承退化状态阶段确定方法 | |
KR20200003664A (ko) | 공기질 데이터의 보정 로직을 생성하는 장치 및 방법 | |
CN109657720B (zh) | 一种电力变压器匝间短路故障的在线诊断方法 | |
CN103489034A (zh) | 预测与诊断在线海流监测数据的方法和装置 | |
CN108846200B (zh) | 一种基于迭代法的准静态桥梁影响线识别方法 | |
CN111595433B (zh) | 航空发动机整机振动传感器的位置确定方法及系统 | |
CN116242443A (zh) | 一种超声波计量仪表的动态计量方法 | |
CN115508770A (zh) | 一种基于kl-nb算法的电能表运行状态在线评估方法 | |
CN109982365B (zh) | 基于仿真和mro数据的天馈问题核查方法及装置 | |
US9771880B2 (en) | Real-time vehicle data acquisition and analysis | |
CN111062538B (zh) | 空间负荷预测中确定元胞负荷合理最大值的ceemd方法 | |
CN114546841A (zh) | 基于云计算的软件质量评估方法 | |
CN108053093A (zh) | 一种基于平均影响值数据变换的k-近邻故障诊断方法 | |
CN114357778A (zh) | 长期监测场景下的索力自适应识别方法 | |
JP3155954B2 (ja) | 電波伝搬損失特性推定方法 | |
CN116879121B (zh) | 一种基于光纤传感技术的空气颗粒物浓度实时监测系统 | |
CN116087622B (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 |