CN113408357A - 一种颤振紊流响应信号到脉冲响应信号的生成方法 - Google Patents

一种颤振紊流响应信号到脉冲响应信号的生成方法 Download PDF

Info

Publication number
CN113408357A
CN113408357A CN202110558839.7A CN202110558839A CN113408357A CN 113408357 A CN113408357 A CN 113408357A CN 202110558839 A CN202110558839 A CN 202110558839A CN 113408357 A CN113408357 A CN 113408357A
Authority
CN
China
Prior art keywords
impulse response
response signal
turbulence
generating
flutter
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
Application number
CN202110558839.7A
Other languages
English (en)
Other versions
CN113408357B (zh
Inventor
段世强
郑华
吴政龙
周江涛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN202110558839.7A priority Critical patent/CN113408357B/zh
Publication of CN113408357A publication Critical patent/CN113408357A/zh
Application granted granted Critical
Publication of CN113408357B publication Critical patent/CN113408357B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B64—AIRCRAFT; AVIATION; COSMONAUTICS
    • B64F—GROUND OR AIRCRAFT-CARRIER-DECK INSTALLATIONS SPECIALLY ADAPTED FOR USE IN CONNECTION WITH AIRCRAFT; DESIGNING, MANUFACTURING, ASSEMBLING, CLEANING, MAINTAINING OR REPAIRING AIRCRAFT, NOT OTHERWISE PROVIDED FOR; HANDLING, TRANSPORTING, TESTING OR INSPECTING AIRCRAFT COMPONENTS, NOT OTHERWISE PROVIDED FOR
    • B64F5/00—Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for
    • B64F5/60—Testing or inspecting aircraft components or systems
    • 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
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M7/00—Vibration-testing of structures; Shock-testing of structures
    • G01M7/02—Vibration-testing by means of a shake table
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/04—Architecture, e.g. interconnection topology
    • G06N3/044—Recurrent networks, e.g. Hopfield networks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/04—Architecture, e.g. interconnection topology
    • G06N3/045—Combinations of networks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/08—Learning methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/02—Preprocessing
    • G06F2218/04—Denoising
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/08—Feature extraction
    • 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

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Artificial Intelligence (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Biomedical Technology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Manufacturing & Machinery (AREA)
  • Transportation (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种颤振紊流响应信号到脉冲响应信号的生成方法，首先构建基于颤振试验的脉冲响应生成系统，该系统为深度学习模型，包括encoder、middle和decoder三部分，其中encoder和decoder以一维卷积、ReLU激活函数和BatchNormal为基本结构组成，middle由LSTM网络结构组成，用于脉冲响应信号的频域特征提取；然后将新紊流响应信号和对应的脉冲响应信号作为脉冲响应生成系统的训练数据，采用Pytorch深度学习框架，得到最终的脉冲响应生成系统；最后将最终的脉冲响应生成系统部署在实际工程中。本发明在保证算法实时性的前提下可以获得较好的模态参数估计结果。

Description

一种颤振紊流响应信号到脉冲响应信号的生成方法
技术领域
本发明属于信号处理技术领域，具体涉及一种脉冲响应信号的生成方法。
背景技术
由于飞行器或者气动弹性结构的应用环境要求，大气紊流激励总是作用于结构系统中，因此，在气动弹性系统的试验试飞中，通过大气紊流激励的响应信号进行结构的模态分析是一种比较常见的结构模态分析方法，然而，由于紊流响应信噪比低，模态参数估计需要复杂的计算；最为理想的模态参数估计是通过脉冲响应实现的，但是，在实际的颤振试验试飞中，脉冲响应的获取难度大、风险高以及数据有效性差等问题，使得基于脉冲响应信号的模态参数辨识在实际的工程试验中难以作为唯一的模态参数辨识方法。进一步地，大气紊流激励通常表现为随机激励，在实际工程中该激励信号无法测量，因此，通过激励和响应信号对结构系统进行分析的方法在紊流激励条件下无法实现，通常需要在仅知响应的情况下进行结构系统的模态参数估计。
现阶段常见的自然激励信号模态参数识别方法主要以随机减量技术、随机子空间以及自回归建模及其扩展方法等为主，以上方法主要存在以下问题：
(1)随机减量技术对单一模态的紊流响应信号而言，可以有效计算对应信号的脉冲响应信号，但是，实际的颤振试飞试验中往往存在密集模态问题，随机减量技术难以获得有效的脉冲响应信号；
(2)随机子空间方法需要遍历多个阶数的模态参数，然后以各阶数的模态参数为特征判定稳定的模态参数(频率、阻尼)，但是，随着遍历阶数的增加，计算模态参数的算法耗时非常高，难以做到实时处理，算法的时间复杂度和空间复杂度对实时的模态参数分析该来的诸多困难，是该方法无法在线处理颤振试飞紊流响应数据的主要缺点。另一方面，由于需要通过多个阶数的模态参数进行稳定模态的分析，现阶段对于稳定模态的问题主要以聚类方法为主，但是针对稳定点的判定，基于响应信号定阶是目前该类方法的主要难点，需要通过人为的判定实现阶数的确定。
(3)自回归建模等相关的方法的主要问题是阶数如何确定，自回归建模方法首先需要确定建模的阶数，然后基于阶数建立系数方程进行参数优化，该过程的第一步便是确定要建模的阶数，目前该步骤往往是通过相关数据处理人员人为实现的。
发明内容
为了克服现有技术的不足，本发明提供了一种颤振紊流响应信号到脉冲响应信号的生成方法，首先构建基于颤振试验的脉冲响应生成系统，该系统为深度学习模型，包括encoder、middle和decoder三部分，其中encoder和decoder以一维卷积、ReLU激活函数和BatchNormal为基本结构组成，middle由LSTM网络结构组成，用于脉冲响应信号的频域特征提取；然后将新紊流响应信号和对应的脉冲响应信号作为脉冲响应生成系统的训练数据，采用Pytorch深度学习框架，得到最终的脉冲响应生成系统；最后将最终的脉冲响应生成系统部署在实际工程中。本发明在保证算法实时性的前提下可以获得较好的模态参数估计结果。
本发明解决其技术问题所采用的技术方案包括如下步骤：
步骤1：将由颤振飞行试验采集的紊流响应信号通过频率分析区间带通滤波之后生成带通滤波后的紊流响应信号；
步骤2：构建基于颤振的脉冲响应生成系统；
步骤2-1：通过仿真生成紊流响应信号以及对应的脉冲响应信号，将紊流响应信号和对应的脉冲响应信号作为脉冲响应生成系统的训练数据；
步骤2-2：将仿真产生的紊流响应信号作为输入，依次经过encoder、middle和decoder，最终的输出与该仿真产生的紊流响应对应的脉冲响应进行误差计算；
其中encoder和decoder均以一维卷积、ReLU激活函数和BatchNormal为基本结构组成，middle由LSTM网络结构组成，用于脉冲响应信号的频域特征提取；
步骤2-3：采用Pytorch深度学习框架，以步骤2-1生成的系统训练数据对脉冲响应生成系统进行网络结构优化，得到最终的脉冲响应生成系统组成和参数信息；
步骤3：将最终的脉冲响应生成系统用于实际工程中颤振试验试飞的紊流响应数据的脉冲响应信号生成。
本发明的有益效果如下：
本发明通过数据驱动的方法，将紊流响应信号转化为对应的脉冲响应信号，通过数据驱动模型获取自然激励响应的脉冲响应信号，然后通过脉冲响应信号获取系统的模态参数，保证算法实时性的前提下可以获得较好的模态参数估计结果。
附图说明
图1为本发明方法脉冲响应信号生成系统流程图。
图2为本发明方法的脉冲响应生成系统框图。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
本发明提出一种脉冲响应信号生成方法，由于结构系统的紊流响应信号与脉冲响应信号之间在数学原理上被认为是卷积过程，且脉冲响应信号很难测量，因此，通过可以测量的紊流响应信号基于生成模型生成脉冲响应信号，对于一个训练完成的生成模型，仅通过模型的加载和计算可实现对实际试验试飞的紊流响应信号生成脉冲响应信号。
如图1所示，一种颤振紊流响应信号到脉冲响应信号的生成方法，包括如下步骤：
步骤1：将由颤振飞行试验采集的紊流响应信号通过分析频率区间带通滤波之后获得带通滤波后的紊流响应信号，该步骤作为实际工程中信号处理的常见手段主要目的在于滤除无关频带的信息，获得分析频率相关的有用信息；
步骤2：构建基于颤振的脉冲响应生成系统，计算紊流响应信号对应的脉冲响应信号，其中脉冲响应信号生成系统作为深度学习模型，需要进行模型预训练，便于实际的工程应用，具体步骤如下：
步骤2-1：通过仿真生成紊流响应信号以及对应的脉冲响应信号，将仿真的紊流响应信号和对应的脉冲响应信号作为脉冲响应生成系统的训练数据；
步骤2-2：如图2所示，将仿真的紊流响应信号作为输入，依次经过encoder、middle和decoder，最终的输出与该仿真产生的紊流响应对应的脉冲响应进行误差计算；
其中encoder和decoder均以一维卷积、ReLU激活函数和BatchNormal为基本结构组成，middle由LSTM网络结构组成，用于脉冲响应信号的频域特征提取；
步骤2-3：采用Pytorch深度学习框架，以步骤2-1生成的仿真训练数据对脉冲响应生成系统进行网络结构优化，得到最终的脉冲响应生成系统组成和参数信息；
步骤3：将最终的脉冲响应生成系统用于实际工程中颤振试验试飞的紊流响应数据的脉冲响应信号生成，包括对部署脉冲响应生成系统和工程试验数据的处理。
脉冲响应信号生成系统的流程如图1所示，图中基于仿真信号训练的脉冲响应信号生成系统根据步骤2中所述的内容建立，通过仿真数据进行脉冲响应信号生成系统的模型参数确定是由于实践中工程试验数据往往不能获得理想的脉冲激励响应信号，无法作为深度学习模型的数据集进行模型的优化和训练，因此，基于仿真数据进行模型的优化，然后部署该模型用于颤振试验试飞的紊流响应信号对应的脉冲响应信号的生成。
具体的实施方式分两部分，包括脉冲响应信号生成系统的网络参数优化和模型的部署应用。
首先，按照深度学习算法的训练与部署原则，针对设计的脉冲响应信号生成系统网络模型建立数据集，对模型进行优化与参数优化，其中，数据集通过仿真建立脉冲响应与紊流响应构建数据集，用于脉冲响应生成系统的网络参数优化，网络优化的过程基于Pytorch深度学习框架实现。
其次，脉冲响应信号生成系统的应用则主要是对预训练的网络进行应该用部署，该步骤主要是在于对实际颤振试验试飞的紊流响应信号进行脉冲响应的生成，用于颤振试验试飞中的模态参数估计。

Claims (1)

1.一种颤振紊流响应信号到脉冲响应信号的生成方法，其特征在于，包括以下步骤：
步骤1：将由颤振飞行试验采集的紊流响应信号通过频率分析区间带通滤波之后生成带通滤波后的紊流响应信号；
步骤2：构建基于颤振的脉冲响应生成系统；
步骤2-1：通过仿真生成紊流响应信号以及对应的脉冲响应信号，将紊流响应信号和对应的脉冲响应信号作为脉冲响应生成系统的训练数据；
步骤2-2：将仿真产生的紊流响应信号作为输入，依次经过encoder、middle和decoder，最终的输出与该仿真产生的紊流响应对应的脉冲响应进行误差计算；
其中encoder和decoder均以一维卷积、ReLU激活函数和BatchNormal为基本结构组成，middle由LSTM网络结构组成，用于脉冲响应信号的频域特征提取；
步骤2-3：采用Pytorch深度学习框架，以步骤2-1生成的系统训练数据对脉冲响应生成系统进行网络结构优化，得到最终的脉冲响应生成系统组成和参数信息；
步骤3：将最终的脉冲响应生成系统用于实际工程中颤振试验试飞的紊流响应数据的脉冲响应信号生成。
CN202110558839.7A 2021-05-21 2021-05-21 一种颤振紊流响应信号到脉冲响应信号的生成方法 Active CN113408357B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110558839.7A CN113408357B (zh) 2021-05-21 2021-05-21 一种颤振紊流响应信号到脉冲响应信号的生成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110558839.7A CN113408357B (zh) 2021-05-21 2021-05-21 一种颤振紊流响应信号到脉冲响应信号的生成方法

Publications (2)

Publication Number Publication Date
CN113408357A true CN113408357A (zh) 2021-09-17
CN113408357B CN113408357B (zh) 2024-02-13

Family

ID=77679233

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110558839.7A Active CN113408357B (zh) 2021-05-21 2021-05-21 一种颤振紊流响应信号到脉冲响应信号的生成方法

Country Status (1)

Country Link
CN (1) CN113408357B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101650883A (zh) * 2009-02-13 2010-02-17 中国人民解放军空军航空大学 大气紊流在飞行模拟器上的一种仿真方法
CN105629237A (zh) * 2015-12-25 2016-06-01 中国航天空气动力技术研究院 利用机载雷达对风场进行危险性评估的方法
CN106874535A (zh) * 2016-12-28 2017-06-20 南京航空航天大学 一种高阶大气紊流成型滤波器的设计方法及仿真方法
CN107860548A (zh) * 2017-09-12 2018-03-30 南京航空航天大学 一种近似在线的颤振边界预测方法
WO2019013992A2 (en) * 2017-06-29 2019-01-17 The Regents Of The University Of Colorado, A Body Corporate PHONONIC REGULATORS FOR CONTROLLING A TURBULENT FLOW
CN109885854A (zh) * 2018-11-23 2019-06-14 南京航空航天大学 基于arma模型的颤振边界实时预测系统及预测方法
CN111898327A (zh) * 2020-06-30 2020-11-06 西北工业大学 一种气动弹性系统颤振信号异常数据扩展方法
CN112580588A (zh) * 2020-12-29 2021-03-30 西北工业大学 一种基于经验模态分解的颤振信号智能识别方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101650883A (zh) * 2009-02-13 2010-02-17 中国人民解放军空军航空大学 大气紊流在飞行模拟器上的一种仿真方法
CN105629237A (zh) * 2015-12-25 2016-06-01 中国航天空气动力技术研究院 利用机载雷达对风场进行危险性评估的方法
CN106874535A (zh) * 2016-12-28 2017-06-20 南京航空航天大学 一种高阶大气紊流成型滤波器的设计方法及仿真方法
WO2019013992A2 (en) * 2017-06-29 2019-01-17 The Regents Of The University Of Colorado, A Body Corporate PHONONIC REGULATORS FOR CONTROLLING A TURBULENT FLOW
CN107860548A (zh) * 2017-09-12 2018-03-30 南京航空航天大学 一种近似在线的颤振边界预测方法
CN109885854A (zh) * 2018-11-23 2019-06-14 南京航空航天大学 基于arma模型的颤振边界实时预测系统及预测方法
CN111898327A (zh) * 2020-06-30 2020-11-06 西北工业大学 一种气动弹性系统颤振信号异常数据扩展方法
CN112580588A (zh) * 2020-12-29 2021-03-30 西北工业大学 一种基于经验模态分解的颤振信号智能识别方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
唐炜;史忠科;: "扫频激励下的飞机颤振模态参数小波辨识研究", 振动与冲击, no. 02, 25 February 2009 (2009-02-25), pages 172 - 177 *
李扬;周丽;杨秉才;: "飞机紊流激励响应的模态参数识别", 振动工程学报, no. 06, 15 December 2016 (2016-12-15), pages 963 - 969 *

Also Published As

Publication number Publication date
CN113408357B (zh) 2024-02-13

Similar Documents

Publication Publication Date Title
CN113707176B (zh) 一种基于声信号及深度学习技术的变压器故障检测方法
CN113609955B (zh) 一种基于深度学习和数字孪生的三相逆变器参数辨识方法
CN106451498B (zh) 一种基于改进广义形态滤波的低频振荡模态辨识方法
CN108562811A (zh) 基于双向长短期记忆的复杂电能质量扰动分析方法
CN110224394A (zh) 适用于非平稳功率振荡信号特征提取的傅里叶分解算法
CN116861544B (zh) 一种基于边云协同的建筑异常振动源定位方法及相关设备
CN114218778A (zh) 一种用于声爆试验数据的分析方法及装置
Li et al. An interpretable aerodynamic identification model for hypersonic wind tunnels
Alguri et al. Transfer learning of ultrasonic guided waves using autoencoders: A preliminary study
CN114662414A (zh) 一种基于图小波神经网络模型的油藏生产预测方法
CN110782041B (zh) 一种基于机器学习的结构模态参数识别方法
CN112200238A (zh) 基于声响特征的硬岩拉剪破裂识别方法与装置
CN111365624A (zh) 一种输卤管道泄漏检测的智能终端与方法
CN111541255A (zh) 基于动力学系统的低频振荡模态识别方法及系统
CN115270239A (zh) 基于动力特性和智能算法响应面法的桥梁可靠性预测方法
Dan et al. An automatic real-time cable modal frequency identification and tracking algorithm by combining recursive band-pass filter and recursive Hilbert transform
CN113408357B (zh) 一种颤振紊流响应信号到脉冲响应信号的生成方法
US6507798B1 (en) Time-frequency dependent damping via Hilbert damping spectrum
CN109632071A (zh) 一种基于时频特征产生水声环境噪声数据的方法及装置
CN111898327B (zh) 一种气动弹性系统颤振信号异常数据扩展方法
CN114838924A (zh) 一种基于风致振动非平稳响应的结构阻尼比识别方法
Duan et al. A novel adaptive fault diagnosis method for wind power gearbox
Chauhan Parameter estimation and signal processing techniques for operational modal analysis
CN113886992A (zh) 一种基于多保真度数据的数字孪生建模方法
Govers et al. Using stochastic experimental modal data for identifying stochastic finite element parameters of the AIRMOD benchmark structure

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