CN107685878B - 一种基于频响分析的飞机动力学监控方法 - Google Patents

一种基于频响分析的飞机动力学监控方法 Download PDF

Info

Publication number
CN107685878B
CN107685878B CN201710756853.1A CN201710756853A CN107685878B CN 107685878 B CN107685878 B CN 107685878B CN 201710756853 A CN201710756853 A CN 201710756853A CN 107685878 B CN107685878 B CN 107685878B
Authority
CN
China
Prior art keywords
frequency
airplane
order
weight
mode
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
Application number
CN201710756853.1A
Other languages
English (en)
Other versions
CN107685878A (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.)
Shenyang Aircraft Design and Research Institute Aviation Industry of China AVIC
Original Assignee
Shenyang Aircraft Design and Research Institute Aviation Industry of China AVIC
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 Shenyang Aircraft Design and Research Institute Aviation Industry of China AVIC filed Critical Shenyang Aircraft Design and Research Institute Aviation Industry of China AVIC
Priority to CN201710756853.1A priority Critical patent/CN107685878B/zh
Publication of CN107685878A publication Critical patent/CN107685878A/zh
Application granted granted Critical
Publication of CN107685878B publication Critical patent/CN107685878B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64FGROUND 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/00Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for
    • B64F5/60Testing or inspecting aircraft components or systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Transportation (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明涉及一种基于频响分析的飞机动力学监控方法,其属于气动弹性试验技术领域,其包括步骤一:获取飞机结构在第i阶模态的基准频率fri及检查频率fti;步骤二:根据权重误差公式得到飞机颤振特性的频率权重误差。本发明的基于频响分析的外场动力学监控方法通过应用频响函数分析及模态权重误差分析,可以在外场通过简单的频率检查进行动力学监控,从而获得飞机在使用过程中的动力学特性并应用于飞机颤振特性评估,最终给出飞机的使用包线,以保证飞机的飞行安全。

Description

一种基于频响分析的飞机动力学监控方法
技术领域
本发明属于气动弹性试验技术领域,尤其涉及一种基于频响分析的飞机动力学监控方法。
背景技术
随着客户需求的快速发展,对飞机平台的使用和维护成本提出了更高的要求,需要飞机平台在满足飞行使用的情况下更多降低使用成本。因此,在外场能够快速监测飞机动力学特性、保障飞机的使用效率同时降低使用成本,是飞机设计的发展方向之一。
同时,针对无人机平台系统,这些无人机对于机体多采用平时拆解、贮存,战时运至战区组装、执行飞行任务的模式。该种作战使用方式降低了平时的贮存保养成本,提高了战时执行任务的效率。但由于需要对飞机进行多次拆解、组装,可能会对飞机的动力学特性产生影响,同时在外场又缺乏充分的试验室条件进行测试,因此需要对外场条件下的动力学监控方法进行研究。
然而,受需求与经费等限制,目前国内缺少应用于飞机的外场监控手段,本发明专利能够提供一种在外场条件下解决飞机动力学监控的可行方法。
发明内容
本发明的目的是提供一种基于频响分析的飞机动力学监控方法,用于解决现有技术中,外场环境下无法测得测得的动力学特性不准的问题。
为达到上述目的,本发明采用的技术方案是:一种基于频响分析的飞机动力学监控方法,其包括
步骤一:获取飞机结构在第i阶模态的基准频率fri及检查频率fti
步骤二:根据权重误差公式得到飞机颤振特性的频率权重误差
Figure BDA0001392374890000021
上式中,ε为频率权重误差,ni为第i阶模态的误差权重,fri为第i阶模态的基准频率,fti为第i阶模态的检查频率。
在本发明优选实施例中,步骤一中,模态的阶数取值范围3阶至10阶,模态的阶数会导致最终计算的精度,高于10阶最会使得计算精度提高,但是增大里计算量,低于3阶则计算精度不够。
在本发明优选实施例中,模态的取值为5阶。
在本发明优选实施例中,5阶模态包括四阶弯曲及一阶扭转,由于弯曲和扭转时两种不同的状态,因此在模态中需每种都包括,即至少包括弯曲和扭转。
在本发明优选实施例中,5阶模态下的误差权重系数分别为:
一阶弯曲时,模态权重为40%~60%;
二阶弯曲时,模态权重为20%~40%;
三阶弯曲时,模态权重为90%~100%;
四阶弯曲时,模态权重为10~30%;
一阶扭转时,模态权重为90~100%。
本发明的基于频响分析的外场动力学监控方法通过应用频响函数分析及模态权重误差分析,可以在外场通过简单的频率检查进行动力学监控,从而获得飞机在使用过程中的动力学特性并应用于飞机颤振特性评估,最终给出飞机的使用包线,以保证飞机的飞行安全。本发明的方法理论基础可靠,监控方法简单,外场实施方便,可对飞机机体结构进行健康监控。依据本发明的方法可以预估飞机大修时间,降低飞机外场再次飞行准备时间,提高飞机利用效率,从而降低全寿命周期使用成本。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分,示出了符合本发明的实施例,并与说明书一起用于解释本发明的原理。
图1为本发明一实施例的实验室环境下得到的基准频响函数。
图2为本发明一实施例的外场环境下得到的检查频响函数。
具体实施方式
为使本发明实施的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行更加详细的描述。
对于飞机全机的动力学特性,通常是通过全机地面振动试验得到的。试验结果包括各阶的固有振动频率、模态、阻尼系数及广义质量,其试验系统都应包括被试飞机、支持系统、激励系统、测量系统、数据采集及处理系统,其试验结果可以应用全机的颤振分析,得到飞机的颤振特性,并最终确定飞机的颤振边界,以保证飞机在使用包线内的安全飞行。
通常,可以通过相位共振法得到较为精确的飞机固有振动频率和模态等,但在外场一般缺乏完善的试验系统,因此需要研究一种相对简便易行的方法对动力学特性进行评估以保障飞机的飞行使用。本发明通过简化的频率检查测出多自由度系统频响函数并分析,应用模态权重方法对外场条件下飞机的动力学特性进行评估。
对于一个多自由度的无阻尼系统,其运动方程为:
Figure BDA0001392374890000031
式中,M为系统的质量矩阵,K为系统的刚度矩阵,X为系统的响应,F为系统的激励力向量。
对式运动方程进行傅里叶变换,得到:
(K-ω2M)X(ω)=F(ω)
其阻抗矩阵为:
Z(ω)=K-ω2M
系统的频响函数H(ω)为阻抗矩阵的逆矩阵:
H(ω)=Z-1(ω)=(K-ω2M)-1
对于N阶自由度系统,其测量点l与激励点p之间的频响函数可表示为:
Figure BDA0001392374890000041
上式表示由激励点p单独引起的l点的响应。当
Figure BDA0001392374890000043
时,
Figure BDA0001392374890000044
系统即处于共振状态,相应的频率即为共振频率。
由于在外场的条件下不易通过模态共振的方法获得精确的飞机动力学特性,而通过频响函数来得到简化的动力学特性是相对可行的。但是,通常的动力学分析中频率的误差主要针对单一模态,而通过频响函数得到的模态频率准确性相对不高,因此在本发明中,进一步引入模态权重并通过整体频率误差来进行动力学分析。
根据颤振理论分析、风洞试验结果分析各阶模态对全机或某一部件颤振的影响,确定权重误差公式如下式所示:
以某大展弦比无人机的机翼颤振为例,在试验室条件下测试得到飞机的目标频响函数,以机翼的频响函数为例如图1所示。在飞机完成拆解、运输及外场组装后,测得的机翼的频响函数如图2所示。两次测试通过频响函数得到的模态频率点如表1所示,图1和图2中所标注的编号与表1中一一对应。
表1 机翼固有特性对比结果
Figure BDA0001392374890000042
然后,根据颤振计算结果分析各阶模态对机翼弯扭颤振的影响,由此确定各阶模态的权重,如表2所示。
表2 某大展弦比机翼模态权重示例
Figure BDA0001392374890000051
频率权重误差计算公式为:
Figure BDA0001392374890000052
通过应用本发明中权重误差公式计算含模态权重的频率误差仅为1.24%,从结果可以看出,实验室环境下测得的基准频率与外场拆装后又重新组装的检查频率的误差小于5%,这一误差水平在工程应用中能够满足要求。
本发明的基于频响分析的外场动力学监控方法通过应用频响函数分析及模态权重误差分析,可以在外场通过简单的频率检查进行动力学监控,从而获得飞机在使用过程中的动力学特性并应用于飞机颤振特性评估,最终给出飞机的使用包线,以保证飞机的飞行安全。本发明的方法理论基础可靠,监控方法简单,外场实施方便,可对飞机机体结构进行健康监控。依据本发明的方法可以预估飞机大修时间,降低飞机外场再次飞行准备时间,提高飞机利用效率,从而降低全寿命周期使用成本。
以上所述,仅为本发明的最优具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。

Claims (1)

1.一种基于频响分析的飞机动力学监控方法,其特征在于,包括
步骤一:获取飞机结构在第i阶模态的基准频率fri及检查频率fti
步骤二:根据权重误差公式得到飞机颤振特性的频率权重误差
Figure FDA0002470615300000011
上式中,ε为频率权重误差,ni为第i阶模态的误差权重,fri为第i阶模态的基准频率,fti为第i阶模态的检查频率;
其中,检查频率的模态阶数为5阶,模态包括四阶弯曲及一阶扭转,5阶模态下误差权重系数分别为:
一阶弯曲时,模态权重为40%~60%;
二阶弯曲时,模态权重为20%~40%;
三阶弯曲时,模态权重为90%~100%;
四阶弯曲时,模态权重为10%~30%;
一阶扭转时,模态权重为90%~100%。
CN201710756853.1A 2017-08-29 2017-08-29 一种基于频响分析的飞机动力学监控方法 Active CN107685878B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710756853.1A CN107685878B (zh) 2017-08-29 2017-08-29 一种基于频响分析的飞机动力学监控方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710756853.1A CN107685878B (zh) 2017-08-29 2017-08-29 一种基于频响分析的飞机动力学监控方法

Publications (2)

Publication Number Publication Date
CN107685878A CN107685878A (zh) 2018-02-13
CN107685878B true CN107685878B (zh) 2020-06-30

Family

ID=61155624

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710756853.1A Active CN107685878B (zh) 2017-08-29 2017-08-29 一种基于频响分析的飞机动力学监控方法

Country Status (1)

Country Link
CN (1) CN107685878B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110949689B (zh) * 2019-12-06 2022-07-22 江西洪都航空工业集团有限责任公司 一种飞机全动平尾旋转模态测量方法
CN112859592B (zh) * 2020-12-29 2022-08-09 中国航空工业集团公司西安飞机设计研究所 一种用于涡桨飞机结构模态频率控制的方法
CN114506471B (zh) * 2022-03-17 2023-06-20 北京机电工程研究所 一种飞机悬挂物的一阶弯曲频率确定方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103995918A (zh) * 2014-04-17 2014-08-20 中国航空工业集团公司沈阳飞机设计研究所 一种机翼变形和振动对飞机传递对准影响的分析方法
CN105912772A (zh) * 2016-04-08 2016-08-31 北京宇航系统工程研究所 一种集频率、振型和频响函数为一体的模型修正方法
CN106055733A (zh) * 2016-05-10 2016-10-26 中国人民解放军国防科学技术大学 多功能结构的动力学参数确定方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6947858B2 (en) * 2003-06-27 2005-09-20 The Boeing Company Methods and apparatus for analyzing flutter test data using damped sine curve fitting
US10634580B2 (en) * 2015-06-04 2020-04-28 The Boeing Company Systems and methods for analyzing flutter test data using damped sine curve fitting with the closed form shape fit

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103995918A (zh) * 2014-04-17 2014-08-20 中国航空工业集团公司沈阳飞机设计研究所 一种机翼变形和振动对飞机传递对准影响的分析方法
CN105912772A (zh) * 2016-04-08 2016-08-31 北京宇航系统工程研究所 一种集频率、振型和频响函数为一体的模型修正方法
CN106055733A (zh) * 2016-05-10 2016-10-26 中国人民解放军国防科学技术大学 多功能结构的动力学参数确定方法

Also Published As

Publication number Publication date
CN107685878A (zh) 2018-02-13

Similar Documents

Publication Publication Date Title
US9240083B2 (en) Rotor system health monitoring using shaft load measurements and virtual monitoring of loads
CN107685878B (zh) 一种基于频响分析的飞机动力学监控方法
Pak Wing shape sensing from measured strain
US20160109324A1 (en) Method and apparatus of multi-axis resonance fatigue test
CN106644458B (zh) 一种直升机尾减速器机匣疲劳试验方法
CN106802240B (zh) 一种旋翼桨毂连接件疲劳试验载荷调试方法
CN104699976A (zh) 一种包含平均应力影响的金属材料多轴高周疲劳失效预测方法
Peeters et al. Modern solutions for ground vibration testing of large aircraft
Chwalowski et al. Preliminary computational analysis of the (hirenasd) configuration in preparation for the aeroelastic prediction workshop
RU2599108C1 (ru) Способ мониторинга нагрузок и накопленной усталостной повреждаемости в условиях эксплуатации самолета
CN102288374A (zh) 同时识别多点随机载荷的试验平台及试验方法
Yeo et al. Rotor structural loads analysis using coupled computational fluid dynamics/computational structural dynamics
Griffith et al. Experimental modal analysis of 9-meter research-sized wind turbine blades
CN102175449A (zh) 基于风力机应变能响应的叶片故障诊断方法
CN104296924A (zh) 动力总成惯性参数测试方法
Peeters et al. Modern solutions for ground vibration testing of small, medium and large aircraft
Corle et al. Transient CFD/CSD tiltrotor stability analysis
Gonzalez et al. Components of a wind tunnel balance: Design and calibration
CN115169162B (zh) 飞机振动环境预计方法、装置及计算机可读存储介质
CN107577902A (zh) 一种基于ukf的飞机疲劳结构剩余寿命预测方法
Yao et al. Low-cost measurement of industrial shock signals via deep learning calibration
Sitaraman et al. Coupled CFD/CSD prediction of rotor aerodynamic and structural dynamic loads for three critical flight conditions
CN104849030B (zh) 用于往复式活塞发动机的连杆动态载荷测试装置
CN106644782B (zh) 一种基于功率谱密度的多轴疲劳裂纹萌生方向预测方法
Gasparetto et al. Experimental modal analysis of an aircraft wing prototype for SAE Aerodesign Competition

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