CN116611175B - 一种大展弦比飞机体自由度颤振的预测方法 - Google Patents
一种大展弦比飞机体自由度颤振的预测方法 Download PDFInfo
- Publication number
- CN116611175B CN116611175B CN202310875463.1A CN202310875463A CN116611175B CN 116611175 B CN116611175 B CN 116611175B CN 202310875463 A CN202310875463 A CN 202310875463A CN 116611175 B CN116611175 B CN 116611175B
- Authority
- CN
- China
- Prior art keywords
- aspect ratio
- ratio aircraft
- aircraft
- high aspect
- degree
- 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
- 238000000034 method Methods 0.000 title claims abstract description 38
- 230000008878 coupling Effects 0.000 claims abstract description 13
- 238000010168 coupling process Methods 0.000 claims abstract description 13
- 238000005859 coupling reaction Methods 0.000 claims abstract description 13
- 238000005452 bending Methods 0.000 claims description 27
- 239000011159 matrix material Substances 0.000 claims description 22
- 230000010355 oscillation Effects 0.000 claims description 11
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 238000002474 experimental method Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 238000009434 installation Methods 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000001687 destabilization Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- 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
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/28—Fuselage, exterior or interior
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- 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)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Mathematical Optimization (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Evolutionary Computation (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Automation & Control Theory (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明涉及一种大展弦比飞机体自由度颤振的预测方法,属于飞行器气动弹性技术领域,解决了现有技术中不能获得大展弦比飞机体自由度颤振产生的具体条件、进而不能对大展弦比飞机体自由度颤振进行定量预测的问题。本发明通过构建大展弦比飞机刚柔耦合模型并进行求解,预测大展弦比飞机体自由度颤振;明确了在给定大展弦比飞机的迎角的情况下,准确地预测获得大展弦比飞机体自由度颤振临界速度;方法简洁清晰,便于航空工程应用;可塑性强;能够在航空工程应用中准确地、便捷地预测获得大展弦比飞机体自由度颤振临界速度,有效避免了大展弦比飞机飞行过程中体自由度颤振的发生,保障飞行员和大展弦比飞机的安全。
Description
技术领域
本发明涉及飞行器气动弹性技术领域,具体涉及一种大展弦比飞机体自由度颤振的预测方法。
背景技术
气动弹性力学的动不稳定性问题主要表现为颤振问题,现有大量关于飞行器颤振的研究,比如中国专利:CN104182560B_飞行器颤振预测分析方法和装置、CN104615863B_一种带控制面的三自由度机翼的颤振边界预测方法、CN112861348B_连续变动压条件下的颤振临界动压预测方法、系统及介质和CN114444332B_一种风洞试验颤振稳定性参数预测方法及装置。近年来,大展弦比飞机在工程应用领域有了越来越多的表现。大展弦比飞机大柔性的特点让其在实际应用中出现了许多常规飞行器未曾出现的问题,体自由度颤振就是其中之一。所以大展弦比飞机需要有能预测体自由度颤振的方法以进行应对。
现有针对大展弦比飞机体自由度颤振的分析,张兰钦等于2020年在《大展弦比飞机体自由度颤振的非线性时域分析方法研究》中存在大展弦比飞机体自由度颤振的非线性时域分析方法,该方法通过刚弹耦合的稳定性分析,给出了大展弦比飞机各结构参数对体自由度颤振的影响,但是对体自由度颤振在什么条件下会发生并无具体说明。由于体自由度颤振对大展弦比飞机的飞行安全造成严重影响,甚至会破坏飞机本身结构,因此对大展弦比飞机体自由度颤振给出定量预测是一个亟待解决的问题。
综上,现有技术中存在不能获得大展弦比飞机体自由度颤振产生的具体条件、进而不能对大展弦比飞机体自由度颤振进行定量预测的问题。
发明内容
鉴于上述问题,本发明提供了一种大展弦比飞机体自由度颤振的预测方法,通过构建大展弦比飞机刚柔耦合模型,联立大展弦比飞机的刚性短周期模态振荡频率与大展弦比飞机的柔性一阶弯曲振动频率并进行求解,能够预测大展弦比飞机体自由度颤振;本发明方法解决了现有技术中不能获得大展弦比飞机体自由度颤振产生的具体条件、进而不能对大展弦比飞机体自由度颤振进行定量预测的问题。
本发明提供了一种大展弦比飞机体自由度颤振的预测方法,包括如下步骤:
步骤S1.针对大展弦比飞机建立大展弦比飞机刚性纵向小扰动方程;其中,该大展弦比飞机刚性纵向小扰动方程包括大展弦比飞机纵向小扰动系统矩阵;
步骤S2.对大展弦比飞机纵向小扰动系统矩阵进行简化处理,得到大展弦比飞机纵向小扰动简化系统矩阵;
求解大展弦比飞机纵向小扰动简化系统矩阵的特征根,计算得到大展弦比飞机的刚性短周期模态振荡频率;
步骤S3.针对大展弦比飞机建立大展弦比飞机机翼结构柔性弯曲振动方程;
步骤S4.求解大展弦比飞机机翼结构柔性弯曲振动方程,计算得到大展弦比飞机的柔性一阶弯曲振动频率;
步骤S5.联立大展弦比飞机的刚性短周期模态振荡频率与大展弦比飞机的柔性一阶弯曲振动频率,得到大展弦比飞机刚柔耦合模型;
获取大展弦比飞机的物理参数和与大展弦比飞机飞行速度相关的性能参数,代入大展弦比飞机刚柔耦合模型,得到大展弦比飞机体自由度颤振的产生条件;
根据大展弦比飞机体自由度颤振的产生条件,对大展弦比飞机体自由度颤振进行预测,获得大展弦比飞机体自由度颤振临界速度。
进一步地,大展弦比飞机纵向小扰动简化系统矩阵表示为:
;
其中,;/>;/>;;/>表示大展弦比飞机的发动机推力;/>表示大展弦比飞机的迎角;/>表示大展弦比飞机的发动机推力方向与大展弦比飞机的机体轴形成的安装角;/>表示大展弦比飞机的升力系数对大展弦比飞机的迎角的偏导数;/>表示动压;/>表示大展弦比飞机的机翼面积;/>表示大展弦比飞机的质量;/>表示大展弦比飞机的俯仰力矩系数对大展弦比飞机的迎角的偏导数;/>表示大展弦比飞机的平均气动弦长;/>表示沿垂直于大展弦比飞机对称平面的大展弦比飞机惯性轴的惯性矩;/>表示大展弦比飞机的俯仰力矩系数对大展弦比飞机的迎角关于时间的变化率的偏导数;/>表示大展弦比飞机的飞行速度;/>表示大展弦比飞机的俯仰力矩系数对大展弦比飞机的俯仰角速度的偏导数。
进一步地,大展弦比飞机的刚性短周期模态振荡频率表示为:
。
进一步地,大展弦比飞机机翼结构柔性弯曲振动方程表示为:
;
该大展弦比飞机机翼结构柔性弯曲振动方程的边界条件为:
;
其中,表示单位展长机翼偏离平衡位置的弯心的位移;/>表示单位展长机翼的弯曲刚度;/>表示机翼单位展长质量;/>表示作用在单位展长机翼的升力;以上参数均是的函数;/>表示机翼的展向坐标;/>表示时间;/>表示机翼展长的一半。
进一步地,大展弦比飞机的柔性一阶弯曲振动频率表示为:
;
其中,为以下方程数值最小的解:
。
进一步地,步骤S5中的所述获取大展弦比飞机的物理参数,具体指根据大展弦比飞机的自身属性,通过计算或测量得到物理参数的具体数值;所述物理参数包括单位展长机翼的弯曲刚度、机翼单位展长质量/>、机翼展长的一半/>、大展弦比飞机的机翼面积、大展弦比飞机的平均气动弦长/>、沿垂直于大展弦比飞机对称平面的大展弦比飞机惯性轴的惯性矩/>、大展弦比飞机的发动机推力方向与大展弦比飞机的机体轴形成的安装角/>和大展弦比飞机的质量/>。
进一步地,步骤S5中的所述获取与大展弦比飞机飞行速度相关的性能参数,具体指通过试验测得性能参数与大展弦比飞机的飞行速度的关系曲线后,再对关系曲线拟合处理得到性能参数与大展弦比飞机飞行速度/>的函数关系式,最后基于函数关系式利用大展弦比飞机飞行速度/>来表示性能参数;性能参数具体包括动压/>、大展弦比飞机的俯仰力矩系数对大展弦比飞机的迎角的偏导数/>、大展弦比飞机的俯仰力矩系数对大展弦比飞机的俯仰角速度的偏导数/>、大展弦比飞机的发动机推力/>和大展弦比飞机的升力系数对大展弦比飞机的迎角的偏导数/>。
进一步地,步骤S5中的所述大展弦比飞机刚柔耦合模型表示为:
。
进一步地,步骤S5中的所述大展弦比飞机体自由度颤振的产生条件,具体指大展弦比飞机的飞行速度和大展弦比飞机的迎角/>之间满足的函数关系。
进一步地,步骤S5中的所述根据大展弦比飞机体自由度颤振的产生条件,对大展弦比飞机体自由度颤振进行预测,获得大展弦比飞机体自由度颤振临界速度,具体指先给定大展弦比飞机的迎角,然后计算得到大展弦比飞机的飞行速度/>,该大展弦比飞机的飞行速度/>就是大展弦比飞机的体自由度颤振临界速度。
与现有技术相比,本发明至少具有现如下有益效果:
(1)本发明的大展弦比飞机体自由度颤振的预测方法,考虑大展弦比飞机独特构型的影响,明确了在给定大展弦比飞机的迎角的情况下,能够准确地预测获得大展弦比飞机体自由度颤振临界速度。
(2)本发明的大展弦比飞机体自由度颤振的预测方法,简洁清晰,便于航空工程应用。
(3)本发明的大展弦比飞机体自由度颤振的预测方法,可塑性强,后续需要更精确的大展弦比飞机的体自由度颤振临界速度解时,只需要用更合适的理论模型替换第一步或第三步的理论模型重新求解即可。
(4)本发明的大展弦比飞机体自由度颤振的预测方法,基于上述有益效果,能够在航空工程应用中准确地、便捷地预测获得大展弦比飞机体自由度颤振临界速度,有效避免了大展弦比飞机飞行过程中体自由度颤振的发生,保障飞行员和大展弦比飞机的安全。
附图说明
附图仅用于示出具体实施例的目的,而并不认为是对本发明的限制。
图1为本发明公开的大展弦比飞机体自由度颤振的预测方法的流程图。
具体实施方式
为了能够更清楚地理解本发明的上述目的、特征和优点,下面结合附图和具体实施方式对本发明进行进一步的详细描述。需要说明的是,在不冲突的情况下,本发明的实施例及实施例中的特征可以相互组合。另外,本发明还可以采用其他不同于在此描述的其他方式来实施,因此,本发明的保护范围并不受下面公开的具体实施例的限制。
当结构弹性振动频率接近刚体运动频率时,在一定情况下会使飞机出现一种有别于经典弯扭耦合颤振的失稳现象,称为体自由度颤振;对于对称构型飞行器,其表现为刚体短周期模态与机翼对称一弯模态耦合诱发的颤振;对于非对称构型飞行器,其表现为刚体滚转模态与机翼反对称一弯模态耦合诱发的颤振。
本发明的一个具体实施例,如图1所示,公开了一种大展弦比飞机体自由度颤振的预测方法,包括如下步骤:
步骤S1.针对大展弦比飞机建立大展弦比飞机刚性纵向小扰动方程;其中,该大展弦比飞机刚性纵向小扰动方程包括大展弦比飞机纵向小扰动系统矩阵:
其中,向量是大展弦比飞机纵向状态矢量,是一个四阶列向量;/>表示大展弦比飞机纵向状态矢量/>对时间的导数;/>为大展弦比飞机纵向控制矢量,是一个二阶列向量;为大展弦比飞机纵向小扰动系统矩阵,是一个4阶方阵;/>为大展弦比飞机纵向小扰动控制矩阵,是一个四行二列的矩阵。
由于后续步骤只会用到大展弦比飞机纵向小扰动系统矩阵,故详细给出/>,表示为:
其中,
;
;
;
;
;
;
;
;
其中,表示大展弦比飞机的飞行速度;/>表示大展弦比飞机的机翼面积;/>表示空气密度;/>表示大展弦比飞机的飞行马赫数;/>表示大展弦比飞机质量;/>表示大展弦比飞机的发动机推力方向与大展弦比飞机的机体轴形成的安装角;/>表示大展弦比飞机的迎角;/>表示大展弦比飞机的发动机推力系数对大展弦比飞机的飞行马赫数的偏导数;/>表示大展弦比飞机的发动机推力系数;/>表示动压;/>表示大展弦比飞机的阻力系数对大展弦比飞机的飞行速度的偏导数;/>表示大展弦比飞机的阻力系数;/>表示大展弦比飞机的发动机推力;/>表示大展弦比飞机的阻力系数对大展弦比飞机的迎角的偏导数;/>表示重力加速度;/>表示大展弦比飞机的升力系数对大展弦比飞机的飞行速度的偏导数;/>表示大展弦比飞机的升力系数;/>表示大展弦比飞机的升力系数对大展弦比飞机的迎角的偏导数;/>表示大展弦比飞机的平均气动弦长;/>表示沿垂直于大展弦比飞机对称平面的大展弦比飞机惯性轴的惯性矩;/>表示大展弦比飞机的俯仰力矩系数对大展弦比飞机的速度的偏导数;/>表示大展弦比飞机的俯仰力矩系数对大展弦比飞机的迎角的偏导数;/>表示大展弦比飞机的俯仰力矩系数对大展弦比飞机的迎角关于时间的变化率的偏导数;/>表示大展弦比飞机的俯仰力矩系数对大展弦比飞机的俯仰角速度的偏导数。
步骤S2.对大展弦比飞机纵向小扰动系统矩阵进行简化处理,得到大展弦比飞机纵向小扰动简化系统矩阵。
由于大展弦比飞机的体自由度颤振仅与大展弦比飞机的刚性短周期模态相关,故将大展弦比飞机纵向小扰动系统矩阵简化为如下二阶方阵/>:
;
其中,是大展弦比飞机纵向小扰动简化系统矩阵。
求解大展弦比飞机纵向小扰动简化系统矩阵的特征根,计算得到大展弦比飞机的刚性短周期模态振荡频率。
具体来说,求矩阵的特征根/>,即求如下特征方程的解:
;
将上述特征方程展开后为关于的一元二次代数方程:
;
该一元二次代数方程为二阶系统振动方程,基于振动系统运动理论,可由上述方程计算得到大展弦比飞机的刚性短周期模态振荡频率;其中,该刚性短周期模态振荡频率同时也是大展弦比飞机的飞行速度/>的函数:
。
步骤S3.基于工程梁理论,针对大展弦比飞机建立大展弦比飞机机翼结构柔性弯曲振动方程:
;
该大展弦比飞机机翼结构柔性弯曲振动方程的边界条件为:
;
其中,表示单位展长机翼偏离平衡位置的弯心的位移;/>表示单位展长机翼的弯曲刚度;/>表示机翼单位展长质量;/>表示作用在单位展长机翼的升力;以上参数均是的函数;/>表示机翼的展向坐标;/>表示时间;/>表示机翼展长的一半。
步骤S4.求解大展弦比飞机机翼结构柔性弯曲振动方程,计算得到大展弦比飞机的柔性一阶弯曲振动频率:
;
其中,为以下方程数值最小的解:
;
为大展弦比飞机的柔性一阶弯曲振动频率,该柔性一阶弯曲振动频率为一常数。
步骤S5.联立大展弦比飞机的刚性短周期模态振荡频率与大展弦比飞机的柔性一阶弯曲振动频率,得到大展弦比飞机刚柔耦合模型:
。
获取大展弦比飞机的物理参数和与大展弦比飞机飞行速度相关的性能参数,代入大展弦比飞机刚柔耦合模型,得到大展弦比飞机体自由度颤振的产生条件。
获取大展弦比飞机的物理参数,具体指根据大展弦比飞机的自身属性,通过计算或测量得到物理参数的具体数值;所述物理参数包括单位展长机翼的弯曲刚度、机翼单位展长质量/>、机翼展长的一半/>、大展弦比飞机的机翼面积/>、大展弦比飞机的平均气动弦长/>、沿垂直于大展弦比飞机对称平面的大展弦比飞机惯性轴的惯性矩/>、大展弦比飞机的发动机推力方向与大展弦比飞机的机体轴形成的安装角/>和大展弦比飞机的质量。
获取与大展弦比飞机飞行速度相关的性能参数,具体指通过试验测得性能参数与大展弦比飞机的飞行速度的关系曲线后,再对关系曲线拟合处理得到性能参数与大展弦比飞机飞行速度/>的函数关系式,最后基于函数关系式利用大展弦比飞机飞行速度/>来表示性能参数;性能参数具体包括动压/>、大展弦比飞机的俯仰力矩系数对大展弦比飞机的迎角的偏导数/>、大展弦比飞机的俯仰力矩系数对大展弦比飞机的俯仰角速度的偏导数/>、大展弦比飞机的发动机推力/>和大展弦比飞机的升力系数对大展弦比飞机的迎角的偏导数/>。
大展弦比飞机体自由度颤振的产生条件,具体指大展弦比飞机的飞行速度和大展弦比飞机的迎角/>之间满足的函数关系。
根据大展弦比飞机体自由度颤振的产生条件,对大展弦比飞机体自由度颤振进行预测,获得大展弦比飞机体自由度颤振临界速度。
具体来说,先给定大展弦比飞机的迎角,然后计算得到大展弦比飞机的飞行速度/>,该大展弦比飞机的飞行速度/>就是大展弦比飞机的体自由度颤振临界速度。
与现有技术相比,本发明的大展弦比飞机体自由度颤振的预测方法,考虑大展弦比飞机独特构型的影响,明确了在给定大展弦比飞机的迎角的情况下,能够准确地预测获得大展弦比飞机体自由度颤振临界速度;方法简洁清晰,便于航空工程应用;可塑性强,后续需要更精确的大展弦比飞机的体自由度颤振临界速度解时,只需要用更合适的理论模型替换第一步或第三步的理论模型重新求解即可;基于上述有益效果,能够在航空工程应用中准确地、便捷地预测获得大展弦比飞机体自由度颤振临界速度,有效避免了大展弦比飞机飞行过程中体自由度颤振的发生,保障飞行员和大展弦比飞机的安全。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。
Claims (10)
1.一种大展弦比飞机体自由度颤振的预测方法,其特征在于,包括如下步骤:
步骤S1.针对大展弦比飞机建立大展弦比飞机刚性纵向小扰动方程;其中,该大展弦比飞机刚性纵向小扰动方程包括大展弦比飞机纵向小扰动系统矩阵;
步骤S2.对大展弦比飞机纵向小扰动系统矩阵进行简化处理,得到大展弦比飞机纵向小扰动简化系统矩阵;
求解大展弦比飞机纵向小扰动简化系统矩阵的特征根,计算得到大展弦比飞机的刚性短周期模态振荡频率;
步骤S3.针对大展弦比飞机建立大展弦比飞机机翼结构柔性弯曲振动方程;
步骤S4.求解大展弦比飞机机翼结构柔性弯曲振动方程,计算得到大展弦比飞机的柔性一阶弯曲振动频率;
步骤S5.联立大展弦比飞机的刚性短周期模态振荡频率与大展弦比飞机的柔性一阶弯曲振动频率,得到大展弦比飞机刚柔耦合模型;
获取大展弦比飞机的物理参数和与大展弦比飞机飞行速度相关的性能参数,代入大展弦比飞机刚柔耦合模型,得到大展弦比飞机体自由度颤振的产生条件;
根据大展弦比飞机体自由度颤振的产生条件,对大展弦比飞机体自由度颤振进行预测,获得大展弦比飞机体自由度颤振临界速度。
2.根据权利要求1所述的大展弦比飞机体自由度颤振的预测方法,其特征在于,大展弦比飞机纵向小扰动简化系统矩阵表示为:
;
其中,;/>;/>;/>;/>表示大展弦比飞机的发动机推力;/>表示大展弦比飞机的迎角;/>表示大展弦比飞机的发动机推力方向与大展弦比飞机的机体轴形成的安装角;/>表示大展弦比飞机的升力系数对大展弦比飞机的迎角的偏导数;/>表示动压;/>表示大展弦比飞机的机翼面积;/>表示大展弦比飞机的质量;/>表示大展弦比飞机的俯仰力矩系数对大展弦比飞机的迎角的偏导数;/>表示大展弦比飞机的平均气动弦长;/>表示沿垂直于大展弦比飞机对称平面的大展弦比飞机惯性轴的惯性矩;/>表示大展弦比飞机的俯仰力矩系数对大展弦比飞机的迎角关于时间的变化率的偏导数;/>表示大展弦比飞机的飞行速度;/>表示大展弦比飞机的俯仰力矩系数对大展弦比飞机的俯仰角速度的偏导数。
3.根据权利要求2所述的大展弦比飞机体自由度颤振的预测方法,其特征在于,大展弦比飞机的刚性短周期模态振荡频率表示为:
。
4.根据权利要求3所述的大展弦比飞机体自由度颤振的预测方法,其特征在于,大展弦比飞机机翼结构柔性弯曲振动方程表示为:
;
该大展弦比飞机机翼结构柔性弯曲振动方程的边界条件为:
;
其中,表示单位展长机翼偏离平衡位置的弯心的位移;/>表示单位展长机翼的弯曲刚度;/>表示机翼单位展长质量;/>表示作用在单位展长机翼的升力;以上参数均是/>的函数;/>表示机翼的展向坐标;/>表示时间;/>表示机翼展长的一半。
5.根据权利要求4所述的大展弦比飞机体自由度颤振的预测方法,其特征在于,大展弦比飞机的柔性一阶弯曲振动频率表示为:
;
其中,为以下方程数值最小的解:
。
6.根据权利要求5所述的大展弦比飞机体自由度颤振的预测方法,其特征在于,步骤S5中的所述获取大展弦比飞机的物理参数,具体指根据大展弦比飞机的自身属性,通过计算或测量得到物理参数的具体数值;所述物理参数包括单位展长机翼的弯曲刚度、机翼单位展长质量/>、机翼展长的一半/>、大展弦比飞机的机翼面积/>、大展弦比飞机的平均气动弦长/>、沿垂直于大展弦比飞机对称平面的大展弦比飞机惯性轴的惯性矩/>、大展弦比飞机的发动机推力方向与大展弦比飞机的机体轴形成的安装角/>和大展弦比飞机的质量。
7.根据权利要求6所述的大展弦比飞机体自由度颤振的预测方法,其特征在于,步骤S5中的所述获取与大展弦比飞机飞行速度相关的性能参数,具体指通过试验测得性能参数与大展弦比飞机的飞行速度的关系曲线后,再对关系曲线拟合处理得到性能参数与大展弦比飞机飞行速度/>的函数关系式,最后基于函数关系式利用大展弦比飞机飞行速度/>来表示性能参数;性能参数具体包括动压/>、大展弦比飞机的俯仰力矩系数对大展弦比飞机的迎角的偏导数/>、大展弦比飞机的俯仰力矩系数对大展弦比飞机的俯仰角速度的偏导数/>、大展弦比飞机的发动机推力/>和大展弦比飞机的升力系数对大展弦比飞机的迎角的偏导数/>。
8.根据权利要求7所述的大展弦比飞机体自由度颤振的预测方法,其特征在于,步骤S5中的所述大展弦比飞机刚柔耦合模型表示为:
。
9.根据权利要求8所述的大展弦比飞机体自由度颤振的预测方法,其特征在于,步骤S5中的所述大展弦比飞机体自由度颤振的产生条件,具体指大展弦比飞机的飞行速度和大展弦比飞机的迎角/>之间满足的函数关系。
10.根据权利要求9所述的大展弦比飞机体自由度颤振的预测方法,其特征在于,步骤S5中的所述根据大展弦比飞机体自由度颤振的产生条件,对大展弦比飞机体自由度颤振进行预测,获得大展弦比飞机体自由度颤振临界速度,具体指先给定大展弦比飞机的迎角,然后计算得到大展弦比飞机的飞行速度/>,该大展弦比飞机的飞行速度/>就是大展弦比飞机的体自由度颤振临界速度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310875463.1A CN116611175B (zh) | 2023-07-18 | 2023-07-18 | 一种大展弦比飞机体自由度颤振的预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310875463.1A CN116611175B (zh) | 2023-07-18 | 2023-07-18 | 一种大展弦比飞机体自由度颤振的预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116611175A CN116611175A (zh) | 2023-08-18 |
CN116611175B true CN116611175B (zh) | 2023-09-12 |
Family
ID=87676726
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310875463.1A Active CN116611175B (zh) | 2023-07-18 | 2023-07-18 | 一种大展弦比飞机体自由度颤振的预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116611175B (zh) |
Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1992013758A1 (en) * | 1991-02-04 | 1992-08-20 | Flex Foil Technology, Inc. | Flexible tailored elastic airfoil section |
CN103745066A (zh) * | 2014-01-21 | 2014-04-23 | 北京航空航天大学 | 一种大展弦比机翼结构刚度指标的确定方法 |
CN103853890A (zh) * | 2014-03-12 | 2014-06-11 | 温州职业技术学院 | 一种高超声速飞行器气动弹性剪裁方法 |
CN104615863A (zh) * | 2015-01-14 | 2015-05-13 | 南京航空航天大学 | 一种带控制面的三自由度机翼的颤振边界预测方法 |
CN104965991A (zh) * | 2015-07-08 | 2015-10-07 | 中国人民解放军军械工程学院 | 基于传递函数的机翼颤振速度确定方法 |
CN104978485A (zh) * | 2015-06-23 | 2015-10-14 | 中国航空工业集团公司西安飞机设计研究所 | 一种大展弦比飞机机翼弯曲刚度的计算方法 |
CN106927022A (zh) * | 2017-03-23 | 2017-07-07 | 清华大学 | 基于自展开折叠翼技术的超大展弦比飞机 |
CN110309579A (zh) * | 2019-06-27 | 2019-10-08 | 复旦大学 | 一种针对弹性飞机阵风响应的仿真分析方法和系统 |
CN113111430A (zh) * | 2021-03-06 | 2021-07-13 | 北京航空航天大学 | 基于非线性气动力降阶的弹性飞机飞行动力学建模方法 |
CN113753257A (zh) * | 2021-09-19 | 2021-12-07 | 中国航空工业集团公司西安飞机设计研究所 | 一种桁架支撑机翼构型飞机展弦比优化方法 |
CN115186531A (zh) * | 2022-06-15 | 2022-10-14 | 北京理工大学 | 基于位姿特性的机器人加工稳定性预测方法及装置 |
CN115186376A (zh) * | 2022-06-21 | 2022-10-14 | 中国航空工业集团公司哈尔滨空气动力研究所 | 一种小展弦比飞翼飞行器滚转力矩模型及其参数辨识方法 |
CN115186444A (zh) * | 2022-06-15 | 2022-10-14 | 北京理工大学 | 基于机器人铣削稳定性提升的工艺调控方法及装置 |
CN115659523A (zh) * | 2022-12-27 | 2023-01-31 | 北京航空航天大学 | 一种大展弦比无人机刚柔耦合建模分析方法 |
CN115655635A (zh) * | 2022-12-14 | 2023-01-31 | 中国空气动力研究与发展中心高速空气动力研究所 | 用于全机体自由度颤振或阵风试验的两自由度支撑系统 |
CN115924079A (zh) * | 2022-12-15 | 2023-04-07 | 上海夏佩尔信息技术有限公司 | 一种大展弦比太阳能飞机快速爬升方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8448893B2 (en) * | 2009-10-26 | 2013-05-28 | Aerion Corporation | Laminar flow wing optimized for transonic cruise aircraft |
US8272594B2 (en) * | 2009-10-26 | 2012-09-25 | Aerion Corporation | Laminar flow wing optimized for supersonic cruise aircraft |
-
2023
- 2023-07-18 CN CN202310875463.1A patent/CN116611175B/zh active Active
Patent Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1992013758A1 (en) * | 1991-02-04 | 1992-08-20 | Flex Foil Technology, Inc. | Flexible tailored elastic airfoil section |
CN103745066A (zh) * | 2014-01-21 | 2014-04-23 | 北京航空航天大学 | 一种大展弦比机翼结构刚度指标的确定方法 |
CN103853890A (zh) * | 2014-03-12 | 2014-06-11 | 温州职业技术学院 | 一种高超声速飞行器气动弹性剪裁方法 |
CN104615863A (zh) * | 2015-01-14 | 2015-05-13 | 南京航空航天大学 | 一种带控制面的三自由度机翼的颤振边界预测方法 |
CN104978485A (zh) * | 2015-06-23 | 2015-10-14 | 中国航空工业集团公司西安飞机设计研究所 | 一种大展弦比飞机机翼弯曲刚度的计算方法 |
CN104965991A (zh) * | 2015-07-08 | 2015-10-07 | 中国人民解放军军械工程学院 | 基于传递函数的机翼颤振速度确定方法 |
CN106927022A (zh) * | 2017-03-23 | 2017-07-07 | 清华大学 | 基于自展开折叠翼技术的超大展弦比飞机 |
CN110309579A (zh) * | 2019-06-27 | 2019-10-08 | 复旦大学 | 一种针对弹性飞机阵风响应的仿真分析方法和系统 |
CN113111430A (zh) * | 2021-03-06 | 2021-07-13 | 北京航空航天大学 | 基于非线性气动力降阶的弹性飞机飞行动力学建模方法 |
CN113753257A (zh) * | 2021-09-19 | 2021-12-07 | 中国航空工业集团公司西安飞机设计研究所 | 一种桁架支撑机翼构型飞机展弦比优化方法 |
CN115186531A (zh) * | 2022-06-15 | 2022-10-14 | 北京理工大学 | 基于位姿特性的机器人加工稳定性预测方法及装置 |
CN115186444A (zh) * | 2022-06-15 | 2022-10-14 | 北京理工大学 | 基于机器人铣削稳定性提升的工艺调控方法及装置 |
CN115186376A (zh) * | 2022-06-21 | 2022-10-14 | 中国航空工业集团公司哈尔滨空气动力研究所 | 一种小展弦比飞翼飞行器滚转力矩模型及其参数辨识方法 |
CN115655635A (zh) * | 2022-12-14 | 2023-01-31 | 中国空气动力研究与发展中心高速空气动力研究所 | 用于全机体自由度颤振或阵风试验的两自由度支撑系统 |
CN115924079A (zh) * | 2022-12-15 | 2023-04-07 | 上海夏佩尔信息技术有限公司 | 一种大展弦比太阳能飞机快速爬升方法 |
CN115659523A (zh) * | 2022-12-27 | 2023-01-31 | 北京航空航天大学 | 一种大展弦比无人机刚柔耦合建模分析方法 |
Non-Patent Citations (1)
Title |
---|
柔性翼飞行器刚柔耦合动态特性研究;陈洋;北京理工大学学报;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN116611175A (zh) | 2023-08-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Stenfelt et al. | Yaw control of a tailless aircraft configuration | |
CN113868771B (zh) | 一种考虑结构和气动非线性的飞行动力学建模方法 | |
Castellani et al. | Flight loads prediction of high aspect ratio wing aircraft using multibody dynamics | |
Danowsky et al. | Control oriented aeroservoelastic modeling of a small flexible aircraft using computational fluid dynamics and computational structural dynamics-invited | |
Koreanschi et al. | Flutter analysis of a morphing wing technology demonstrator: numerical simulation and wind tunnel testing | |
CN116611175B (zh) | 一种大展弦比飞机体自由度颤振的预测方法 | |
Tong et al. | Dynamic response analysis under atmospheric disturbances for helicopters based on elastic blades | |
Cooper et al. | Wind tunnel testing of a high aspect ratio wing model | |
Fouda et al. | Effect of Wing Planform on Airplane Stability and Control Authority in Stall | |
Li et al. | A comparison of the fixed-axes and the mean-axes modeling methods for flexible aircraft simulation | |
Rizzi et al. | Virtual aircraft design and control of transcruiser-A canard configuration | |
Tal et al. | Unsteady aeroservoelastic modeling of flexible wing generic transport aircraft with variable camber continuous trailing edge flap | |
Bansal et al. | Stochastic variations in aerodynamic influence coefficients (aics) on flutter prediction of a generic wing | |
GILBERT et al. | The effects of aeroelastic deformation on the unaugmented stopped-rotor dynamics of an X-Wing aircraft | |
Chinmaya et al. | Aeroelasticity-in general and flutter phenomenon | |
CN115438603B (zh) | 一种移动风场环境中弹性飞机风场动响应的确定方法 | |
Wang et al. | Longitudinal flight dynamics and control of highly flexible solar UAV | |
Szymanski et al. | Experimental Investigation of the Static and Dynamic Stability Derivatives Using a Novel 3-DoF Mechanism | |
Murua et al. | A discrete-time state-space model with wake interference for stability analysis of flexible aircraft | |
CN114625159B (zh) | 一种基于被控变量的结冰飞机控制方法 | |
Meng et al. | Static Aeroelastic Effects on Rudder Efficiency of Flying Wing Aircraft | |
Murua et al. | T-tail flutter: potential-flow modelling and experimental validation | |
Li et al. | On the modeling of asymmetric disturbance effect and rejection control for fixed-wing aircraft | |
Güner et al. | New aerodynamic modeling methodology for transonic flutter calculations | |
Ullah et al. | Effect of Reduced Frequency and Mach Number on Damping Characteristics of SACCON UCAV |
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 |