CN115169162B - 飞机振动环境预计方法、装置及计算机可读存储介质 - Google Patents
飞机振动环境预计方法、装置及计算机可读存储介质 Download PDFInfo
- Publication number
- CN115169162B CN115169162B CN202211083308.8A CN202211083308A CN115169162B CN 115169162 B CN115169162 B CN 115169162B CN 202211083308 A CN202211083308 A CN 202211083308A CN 115169162 B CN115169162 B CN 115169162B
- Authority
- CN
- China
- Prior art keywords
- vibration environment
- predicting
- aircraft
- strain
- target point
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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
Abstract
本发明提供飞机振动环境预计方法、装置及计算机可读存储介质,其中的方法将设备支架看做被动端,将支撑设备支架的飞机结构整体看做主动端,被动端目标点为设备安装位置,即为振动环境预计位置;获取计算模型,在该计算模型中目标点响应量通过各个路径输入点的相关系数、频响函数以及应变的乘积来计算;获取新工况下各个路径输入点的各方向的应变,并代入所述计算模型,用所得结果进行目标点处的振动环境预计。
Description
技术领域
本发明涉及飞机特定位置的振动环境预计方法、装置及计算机可读存储介质,尤其涉及机载设备安装位置的振动环境预计方法。
背景技术
飞机飞行过程中,其结构和机载产品会经受多种不同的使用环境。其中振动环境,由于其作用的严重性和持久性,是比较严酷的一种使用环境。实践表明,振动对飞机的寿命和技术效能有严重的影响,甚至危及飞机安全,许多机载产品的使用故障或者事故也均与振动问题有关。
对新研和改型飞机而言,飞行振动环境的预计与分析非常重要:一方面为飞机振动环境适应性和可靠性设计提供基础数据;另一方面为机载设备、系统研制试验及环境和可靠性鉴定试验提供数据。
业内对飞机振动环境预计发展出了三种预计方法:1)计算法;2)统计法;3)经验公式法。这三种预计方法或预计准确度较低,或数据来自以往型号,如采用这些方法来预计不同型号飞机飞行振动环境,可能导致后期振动环境试验过试验或欠试验。近年来发展出的一种基于神经网络的新方法,由于其预计速度快、精度高,有了一定的应用,然而预计过程中需要大量的样本训练,实际情况可能并不允许。本申请发明人在公开号为“CN114186349A”的中国发明专利申请中公开了名称为“一种特定工况下飞机结构振动情况预计方法”的发明,其基于试验工况下飞机结构上特定工况下易于直接测量振动响应点的时域信息、特定工况下不易直接测量振动响应点的时域信息,得到飞机结构试验工况下的振动传递率矩阵,其后以飞机结构试验工况下的振动传递率矩阵、特定工况下飞机结构上特定工况下易于直接测量振动响应点的时域信息,计算飞机结构上特定工况下不易直接测量振动响应点的时域信息,以此,实现对特定工况下飞机结构振动情况的预计。如果存在多种预计方法,则有利于通过不同预计方法的彼此验证来提高预计的准确性。
发明内容
本发明的目的在于提供一种飞机振动环境预计方法、装置及计算机可读存储介质,用于提高预计结果的准确性。
根据本发明一方面的飞机振动环境预计方法,其包括如下步骤:
将设备支架看做被动端,将支撑设备支架的飞机结构整体看做主动端,被动端目标点为设备安装位置,即为振动环境预计位置;
获取计算模型:
,其中为目标点响应量,i为整数变量,u为传递路径数量,、、分别为与第i个载荷识别位置的第一方向、第二方向、第三方向弹性模量以及应变、截面形状相关的系数,、、分别为第i个载荷识别位置第一方向、第二方向、第三方向的载荷到目标点的频响函数,、、分别为第i个载荷识别位置第一方向、第二方向、第三方向的应变,为频率;
测试新工况下各个载荷识别位置的第一方向、第二方向、第三方向的应变,并代入所述计算模型,用所得结果进行目标点处的振动环境预计。
在一个或多个实施例中,获取所述计算模型的步骤包括:测试从主动端到所述目标点的频响函数;测试从主动端到被动端的显示点的频响函数;测试目标点、显示点的响应;测试载荷识别位置的应变;依据各项测试获得数据进行参数识别,得到、、,进而得到所述计算模型。
在一个或多个实施例中,获取所述计算模型的步骤中,采用正弦扫频激励对所述主动端进行激励,以进行各项所述测试。
在一个或多个实施例中,采用激振器对所述主动端进行激励,以对所述计算模型进行验证。
在一个或多个实施例中,所述主动端包括支架梁,所述传递路径包括所述支架梁向所述被动端的各个载荷传递路径。
在一个或多个实施例中,所述主动端和所述被动端通过紧固件连接。
在一个或多个实施例中,所述主动端和所述被动端铆接。
在一个或多个实施例中,所述第一方向为飞机的航向,所述第二方向为飞机的垂向,所述第三方向垂直于所述第一方向和所述第二方向。
根据本发明另一方面的计算机可读存储介质,其上存储有计算机程序,所述程序被处理器执行时实现任一所述的飞机振动环境预计方法。
根据本发明又一方面的飞机振动环境预计装置,包括传感器和控制器,所述飞机振动环境预计装置将设备支架看做被动端,将支撑设备支架的飞机结构整体看做主动端,被动端目标点为设备安装位置,即为振动环境预计位置,所述传感器用于测量新工况下各个载荷识别位置的第一方向、第二方向、第三方向的应变,所述控制器包括存储器和处理器,所述存储器存储有计算机程序,所述程序被处理器执行时实现任一所述的飞机振动环境预计方法。
测试数据表明,基于应变传递路径模型的模型预计值与真实测试值基本保持一致,故可以获得结构的振动传递特性,不仅可以预测结构在不同激励下的响应,还可以研究结构本身并有针对性的对结构进行优化和改进,当结构存在损伤或裂纹时,相较于响应点到响应点构建的计算模型,基于应变的传递路径分析方法更为实用,更能直接反映出问题的本质。
附图说明
本发明的上述的以及其他的特征、性质和优势将通过下面结合附图和实施例的描述而变得更加明显,其中:
图1是根据一个或多个实施例的机舱的立体图。
图2是根据一个或多个实施例的机舱的后视图。
图3是沿图2中I-I方向的剖视图。
图4是新工况下目标点的航向振动环境预计曲线与测试响应曲线的对比图。
图5是新工况下目标点的垂向振动环境预计曲线与测试响应曲线的对比图。
具体实施方式
图1为某飞机典型的机舱1的立体图,该机舱1用于支撑设备支架2,设备支架2用于支撑设备或机载设备,在图示实施例中,设备支架2包括两支架梁。机舱1可以是飞机真实机舱,也可以试验场所的模拟机舱。在图1中设备模拟盒3模拟机载设备安装位置,设备模拟盒3可以替换为真实设备。
飞机飞行过程中,飞机飞行振动主要来自如下四个方面:(1)发动机振动对飞机结构的激励;(2)沿飞机结构外部的气动扰流;(3)飞机外挂装置的气动扰流,特别是开放式结构舱的气动扰流和共振;(4)飞机的机动飞行、气动抖振、着陆滑行等引起的振动。气动载荷及其它振动载荷作用在机舱1的蒙皮及主梁结构上,通过连接结构传递至设备支架2,再传递至设备3的安装位置,在此安装位置产生加速度响应,设备模拟盒3处选择为振动环境预计目标点。
后述实施例基于应变传递路径方法,进行飞机典型设备舱结构的振动环境预计,由于应变反映的是结构的受力与变形情况,因此,当结构存在损伤或裂纹时,基于应变的传递路径分析方法更为实用,更能直接反映出问题的本质。应变传递路径方法的计算模型如下:
为系统第k个振动响应输出;为系统第i个输入载荷;为第i个输入到第j个输出的传递函数,这里为频响函数;为频率;u为传递路径数量。参照公式(3)和(4)进行载荷识别,载荷识别位置通常为设备支架2的支架梁处。公式(3)为支架梁轴向力表达式,为与支架梁轴向弹性模量,支架梁截面,及轴向应变关联的系数。公式(4)为支架梁切向力表达式,为与支架梁切向弹性模量,支架梁截面,及切向应变关联的系数。其轴向应变和切向应变可以通过应变片测量得到。
具体而言,将设备支架2看做被动端,将支持设备支架的飞机结构整体的结构看做主动端,被动端的设备模拟盒3目标点为机载设备安装位置,即为振动环境预计位置。机舱1作为飞机结构整体的一个局部在图1中示出。如图2和图3所示,选取图1中机舱1的下层的中间的设备模拟盒3上响应点为目标点y,还选取设备模拟盒3的另外两个响应点为了提高模型预计精度的额外显示点u1、u2。目标点y振动响应由设备支架2振动传递叠加,选取F1、F2、F3、F4为载荷识别位置或者第i个载荷识别位置(i=1,2,3,4),载荷识别位置也可理解为路径输入点,这样包含了所有引起目标点y振动的传递路径。在设备支架2和设备模拟盒3(或者真实设备)为铆接或者类似铆接的连接方式时,以上所有响应及载荷均考虑两个方向,X方向和Z方向(飞机的航向和垂向),即载荷识别处的设备支架应变测量轴向和切向两个方向。在其他实施例中,根据连接结构的不同,可以考虑X方向、Y方向以及Z方向等三个方向的相应、载荷或应变测量。
在一实施例中,采用正弦扫频激励对所述主动端即进行激励,例如对机舱1的支架梁11进行激励,以构建计算模型。在又一实施例中,采用激振器对对所述主动端即进行激励,以验证所述计算模型。
在采用正弦扫频激励对所述主动端即进行激励的工况,目标点y和额外显示点u1、u2的响应可以表示为以下方程,该响应为加速度响应:
(2)
其中,为频率,、、分别为目标点y、额外显示点u1、u2的响应,i为整数变量,对应振动传递路径,、 、分别为第i个载荷识别位置轴向力到目标点y、额外显示点u1、u2的频响函数,、、分别为第i个载荷识别位置切向力到目标点y、额外显示点u1、u2的频响函数,、为第i个载荷识别位置处的轴向力和切向力,、分别为与第i个载荷识别位置的轴向、切向弹性模量以及应变、截面形状相关的系数,、第i个载荷识别位置轴向、切向的应变。
G iz、G iq、 G 1 iz、G 1 iq、 G 2 iz、G 2 iq为各频响函数与各应变的乘积,均为已知项。
方程(2)为超定方程组,采用SVD(奇异值分解方法)来求解。即采用前述工况数据求解出了参数和(i=1,...,4)。利用新工况的数据代入公式(5),获得计算模型,可对该工况下目标点的振动环境进行计算。
在其他实施例中,根据连接结构的不同,可以考虑X方向、Y方向以及Z方向等三个方向的相应、载荷或应变测量,相应的计算模型可表达如下:
(6),其中为目标点响应量,i为整数变量,u为传递路径数量,、、分别为与第i个载荷识别位置的第一方向、第二方向、第三方向弹性模量以及应变、截面形状相关的系数,、、分别为第i个载荷识别位置第一方向、第二方向、第三方向的载荷到目标点的频响函数,、、分别为第i个载荷识别位置第一方向、第二方向、第三方向的应变,为频率。
以图1至图3中飞机机舱结构为例,再进一步详细阐述基于应变传递路径的振动环境预计方法的过程。
步骤一:进行振动环境预计模型搭建的数据测试,需要测试的数据有以下两类:
a) 频响函数测试:
b) 工况测试量(以下各量均在运行工况下测试):
图4、图5为通过搭建的振动环境预计模型对新工况下的目标点进行预计,并与该工况下测试值进行对比,从图4和图5可得模型预计值与真实测试值在30Hz~2000Hz的频段中变化基本保持一致。比较新工况下预计值与测试值在30Hz~2000Hz频段内的均方根值,误差分别为12.6%和9.4%。
前述振动预计方法可以通过计算机程序执行,并将该计算机程序存储于计算机可读存储介质,所述程序被处理器执行时实现所述的飞机振动环境预计方法。
前述振动预计方法还可以用于飞机的实际飞行工况中,对飞行工况中不便于实时测量的目标点振动情况进行预计,这样的飞机振动环境预计装置包括传感器和控制器,所述飞机振动环境预计装置将设备支架看做被动端,将支撑设备支架的飞机结构整体看做主动端,被动端目标点为设备安装位置,即为振动环境预计位置,所述传感器用于测量新工况下各个载荷识别位置的第一方向、第二方向、第三方向的应变,所述控制器包括存储器和处理器,所述存储器存储有计算机程序,所述程序被处理器任一所述的飞机振动环境预计方法。
本发明虽然以较佳实施例公开如上,但其并不是用来限定本发明,任何本领域技术人员在不脱离本发明的精神和范围内,都可以做出可能的变动和修改。因此,凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何修改、等同变化及修饰,均落入本发明权利要求所界定的保护范围之内。
Claims (10)
1.飞机振动环境预计方法,其特征在于,包括如下步骤:
将设备支架看做被动端,将支撑设备支架的飞机结构整体看做主动端,被动端目标点为设备安装位置,即为振动环境预计位置;
获取计算模型:
其中为目标点响应量,i为整数变量,,u为传递路径数量,、、分别为与第i个载荷识别位置的第一方向、第二方向、第三方向弹性模量以及应变、截面形状相关的系数,、、分别为第i个载荷识别位置第一方向、第二方向、第三方向的载荷到目标点的频响函数,、、分别为第i个载荷识别位置第一方向、第二方向、第三方向的应变,为频率;
获取新工况下各个载荷识别位置的第一方向、第二方向、第三方向的应变,并代入所述计算模型,用所得结果进行目标点处的振动环境预计。
3.如权利要求2所述的飞机振动环境预计方法,其特征在于,获取所述计算模型的步骤中,采用正弦扫频激励对所述主动端进行激励,以进行各项所述测试。
4.如权利要求2或3所述的飞机振动环境预计方法,其特征在于,采用激振器对所述主动端进行激励,以对所述计算模型进行验证。
5.如权利要求1所述的飞机振动环境预计方法,其特征在于,所述主动端包括支架梁,所述传递路径包括所述支架梁向所述被动端的各个载荷传递路径。
6.如权利要求1所述的飞机振动环境预计方法,其特征在于,所述主动端和所述被动端通过紧固件连接。
7.如权利要求6所述的飞机振动环境预计方法,其特征在于,所述主动端和所述被动端铆接。
8.如权利要求1所述的飞机振动环境预计方法,其特征在于,所述第一方向为飞机的航向,所述第二方向为飞机的垂向,所述第三方向垂直于所述第一方向和所述第二方向。
9.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述程序被处理器执行时实现如权利要求1所述的飞机振动环境预计方法。
10.飞机振动环境预计装置,包括传感器和控制器,其特征在于,所述飞机振动环境预计装置将设备支架看做被动端,将支撑设备支架的飞机结构整体看做主动端,被动端目标点为设备安装位置,即为振动环境预计位置,所述传感器用于测量新工况下各个载荷识别位置的第一方向、第二方向、第三方向的应变,所述控制器包括存储器和处理器,所述存储器存储有计算机程序,所述程序被处理器执行时实现如权利要求1所述的飞机振动环境预计方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211083308.8A CN115169162B (zh) | 2022-09-06 | 2022-09-06 | 飞机振动环境预计方法、装置及计算机可读存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211083308.8A CN115169162B (zh) | 2022-09-06 | 2022-09-06 | 飞机振动环境预计方法、装置及计算机可读存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115169162A CN115169162A (zh) | 2022-10-11 |
CN115169162B true CN115169162B (zh) | 2023-02-03 |
Family
ID=83482233
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211083308.8A Active CN115169162B (zh) | 2022-09-06 | 2022-09-06 | 飞机振动环境预计方法、装置及计算机可读存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115169162B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115962887B (zh) * | 2023-03-16 | 2023-06-09 | 湖南云箭科技有限公司 | 试验类推法计算飞机外挂物挂机振动疲劳载荷谱的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103995918A (zh) * | 2014-04-17 | 2014-08-20 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种机翼变形和振动对飞机传递对准影响的分析方法 |
CN113408200A (zh) * | 2021-06-17 | 2021-09-17 | 中国航空综合技术研究所 | 航空装备振动环境分析和预测方法 |
CN114813005A (zh) * | 2022-06-29 | 2022-07-29 | 中国飞机强度研究所 | 一种飞机部件振动疲劳特性测试系统及其方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7338259B2 (en) * | 2004-03-02 | 2008-03-04 | United Technologies Corporation | High modulus metallic component for high vibratory operation |
-
2022
- 2022-09-06 CN CN202211083308.8A patent/CN115169162B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103995918A (zh) * | 2014-04-17 | 2014-08-20 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种机翼变形和振动对飞机传递对准影响的分析方法 |
CN113408200A (zh) * | 2021-06-17 | 2021-09-17 | 中国航空综合技术研究所 | 航空装备振动环境分析和预测方法 |
CN114813005A (zh) * | 2022-06-29 | 2022-07-29 | 中国飞机强度研究所 | 一种飞机部件振动疲劳特性测试系统及其方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115169162A (zh) | 2022-10-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8393206B1 (en) | Dry wind tunnel system | |
EP2904368B1 (en) | Turbine blade fatigue life analysis using non-contact measurement and dynamical response reconstruction techniques | |
JP4992084B2 (ja) | 構造物の損傷の診断システムおよび方法 | |
CN106525226B (zh) | 一种基于现场振动载荷识别的评估方法及系统 | |
US9183180B2 (en) | Method for in-flight assessment of freedom from flutter of an airplane | |
Van der Auweraer | Structural dynamics modeling using modal analysis: applications, trends and challenges | |
US20130297266A1 (en) | Method For Improving Determination Of Mode Shapes For A Mechanical Structure And Applications Hereof | |
Moncayo et al. | Identification and monitoring of modal parameters in aircraft structures using the natural excitation technique (NExT) combined with the eigensystem realization algorithm (ERA) | |
Mao et al. | The construction and comparison of damage detection index based on the nonlinear output frequency response function and experimental analysis | |
CN115169162B (zh) | 飞机振动环境预计方法、装置及计算机可读存储介质 | |
Özsoy et al. | An accelerated life test approach for aerospace structural components | |
EP2957883A2 (en) | Method for prognostics of an aircraft structure based on structural testing | |
Żółtowski et al. | The use of modal analysis in the evaluation of welded steel structures. | |
Pagani et al. | Static and dynamic testing of a full-composite VLA by using digital image correlation and output-only ground vibration testing | |
Luczak et al. | Identification of the test setup influence on the modal properties of a short wind turbine blade during fatigue test | |
Kwon | Uncertainty of bridge flutter velocity measured at wind tunnel tests | |
Derkevorkian et al. | Computational studies of a strain-based deformation shape prediction algorithm for control and monitoring applications | |
Luczak et al. | Static and dynamic testing of the full scale helicopter rotor blades | |
Medeiros et al. | Vibration-Based damage identification applied for composite plate: Experimental analyses | |
Martins et al. | Detecting damage in a UAV composite wing spar using distributed fiber optic strain sensors | |
Petrone et al. | An innovative health monitoring system for aircraft landing gears | |
Mehboob et al. | Numerical study for evaluation of a vibration based damage index for effective damage detection | |
Luczak et al. | Contact Versus Non‐Contact Measurement of a Helicopter Main Rotor Composite Blade | |
Mironov et al. | The demonstrator of helicopter Structural Health Monitoring technique | |
CN111380476B (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 |