CN108446462B - 飞行器颤振分析网格模型埃米特建模方法 - Google Patents
飞行器颤振分析网格模型埃米特建模方法 Download PDFInfo
- Publication number
- CN108446462B CN108446462B CN201810172987.3A CN201810172987A CN108446462B CN 108446462 B CN108446462 B CN 108446462B CN 201810172987 A CN201810172987 A CN 201810172987A CN 108446462 B CN108446462 B CN 108446462B
- Authority
- CN
- China
- Prior art keywords
- mesh point
- formula
- flutter
- body shafting
- function
- 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
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
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- 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
- 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，以飞行器机体轴系分析复杂颤振模型，在机体轴系选取个网格点，坐标为:,振动时第个网格点坐标 为时间和其它两轴位置的函数，为了便于表达第个网格点在轴的振动分量，以为例，下标为网格点标号，下标第二个字母分别表示振动在机体轴系的三个轴分量；为了简化问题，考虑第个网格点在轴方向振动时，考虑第个网格点在轴方向振动时， ，，考虑第个网格点在轴方向振动时，；为了便于书写，将、和简写为、和；
在网格点邻域内建立的近似模型为：
公式1
式中，为在机体轴系网格点的邻域内轴向振动函数，、为轴向振动方程的结构系数函数，分别为在机体轴系网格点，处轴向振动时对应于网格点的变化值；为在机体轴系网格点的邻域内轴向振动函数，、为轴向振动方程的结构系数函数，分别为在机体轴系网格点，处轴向振动时对应于网格点的变化值；为在机体轴系网格点的邻域内轴向振动函数，、为轴向振动方程的结构系数函数，分别为在机体轴系网格点处轴向振动时对应于网格点的变化值；为在网格点的等效激励函数，为时间；为参数向量，表示网格点的温度，为飞行高度，为马赫数，为网格点的气流环境影响，为大气密度；
步骤2，对应步骤1的机体轴系网格点，安装微型温度传感器，、、三个轴向的气流和位置以及振动传感器，在机翼上下方和所有舵面两边安装微型、、三个轴向的气流和位置以及振动传感器，同时在机身加装大于1000帧/秒的图像传感器记录观测机翼翼尖、所有舵面的振动幅值和频率；飞机机载传感器记录时间、飞行高度、马赫数，大气密度；
步骤3，将飞行器到达给定高度和马赫数后颤振试验的过程表达成有效颤振飞行试验，有效颤振飞行试验数据采样时间为， 为正整数， 为记录数据的采样周期，为有效颤振飞行试验的总采样次数；通过颤振飞行试验获得机体轴系网格点，在采样时间时刻的测量值、、和测量值；
步骤4，根据机体轴系网格点，安装微型、、轴向气流传感器，在机翼上下方和所有舵面两边安装微型、、轴向气流传感器，确定时刻机体轴系网格点，的激励函数；
对、、分别采用给定函数逼近，得到，
且关于连续可导，关于连续可导，关于连续可导；这样，可得：
以及;
步骤5， 令：
以及
可将公式1描述成：
公式2
令、、，
式中，
，，对应的埃米特展开系数；
、、为对应于的埃米特展开的阶次；
为阶埃米特正交多项式的递推形式，可得
式中，
其余;
令
，
式中：
、,，、,,、,为对应的埃米特级数的系数，
可得
或写成
公式3
以公式3第一项为例，对
两边求偏导数，可得
根据步骤3和步骤4得到的、、和，和的测试值，可得：
公式4
式中，
进而可得，代入
可以按照下式及最小二乘估计得出
公式5。
本发明的有益结果是：在飞行器机体轴系选择多个网格点，考虑不同飞行速度、大气密度、气流环境、不同温度等气动力和强度变化影响下按照机体轴系分解方法表示复杂颤振网格模型，根据建立该模型的要求提出安装传感器和数据、图像记录要求，通过有效颤振飞行试验获取数据，通过气流传感器测量值获得激励函数，通过气流传感器测量值获得激励函数，采用埃米特函数对振动变量进行逼近和等效描述，按照辨识方法同时确定了机体轴系网格点处三个轴向振动方程求解，从而给出了完整的复杂颤振模型网格模型建模技术方案，解决了现有技术不能有效表达气动力和强度变化影响下复杂颤振模型的技术问题。
下面结合具体实例对本发明作详细说明。
具体实施方式
步骤1，以飞行器机体轴系分析复杂颤振模型，在机体轴系选取个网格点，坐标为:, 振动时第个网格点坐标 为时间和其它两轴位置的函数，为了便于表达第个网格点在轴的振动分量，以为例，下标为网格点标号，下标第二个字母分别表示振动在机体轴系的三个轴分量；为了简化问题，考虑第个网格点在轴方向振动时，考虑第个网格点在轴方向振动时， ，，考虑第个网格点在轴方向振动时，；为了便于书写，将、和简写为、和；
在网格点邻域内建立的近似模型为：
公式1
式中，为在机体轴系网格点的邻域内轴向振动函数，、为轴向振动方程的结构系数函数，分别为在机体轴系网格点，处轴向振动时对应于网格点的变化值；为在机体轴系网格点的邻域内轴向振动函数，、为轴向振动方程的结构系数函数，分别为在机体轴系网格点，处轴向振动时对应于网格点的变化值；为在机体轴系网格点的邻域内轴向振动函数，、为轴向振动方程的结构系数函数，分别为在机体轴系网格点处轴向振动时对应于网格点的变化值；为在网格点的等效激励函数，为时间；为参数向量，表示网格点的温度，为飞行高度，为马赫数，为网格点的气流环境影响，为大气密度；
步骤2，对应步骤1的机体轴系网格点，安装微型温度传感器，、、三个轴向的气流和位置以及振动传感器，在机翼上下方和所有舵面两边安装微型、、三个轴向的气流和位置以及振动传感器，同时在机身加装大于1000帧/秒的图像传感器记录观测机翼翼尖、所有舵面的振动幅值和频率；飞机机载传感器记录时间、飞行高度、马赫数，大气密度；
步骤3，将飞行器到达给定高度和马赫数后颤振试验的过程表达成有效颤振飞行试验，有效颤振飞行试验数据采样时间为， 为正整数， 为记录数据的采样周期，为有效颤振飞行试验的总采样次数；通过颤振飞行试验获得机体轴系网格点，在采样时间时刻的测量值、、和测量值；
步骤4，根据机体轴系网格点，安装微型、、轴向气流传感器，在机翼上下方和所有舵面两边安装微型、、轴向气流传感器，确定时刻机体轴系网格点，的激励函数；
对、、分别采用给定函数逼近，得到
且关于连续可导，关于连续可导，关于连续可导；这样，可得：
以及;
步骤5，令
以及
可将公式1描述成：
公式2
令、、，
式中：
，，对应的埃米特展开系数；
、、为对应于的埃米特展开的阶次；
为阶埃米特正交多项式的递推形式，可得
式中，
其余;
令
，
式中：
、,，、,,、,为对应的埃米特级数的系数，
可得
或写成
公式3
以公式3第一项为例，对
两边求偏导数，可得
根据步骤3和步骤4得到的、、和，和的测试值，可得：
公式4
式中，
进而可得，，代入
可以按照下式及最小二乘估计得出
公式5。
Claims (1)
1.一种飞行器颤振分析网格模型埃米特建模方法，其特征包括以下步骤：
步骤1：以飞行器机体轴系分析复杂颤振模型，在机体轴系选取个网格点，坐标为:, 振动时第个网格点坐标 为时间和其它两轴位置的函数，为了便于表达第个网格点在轴的振动分量，以为例，下标为网格点标号，下标第二个字母分别表示振动在机体轴系的三个轴分量；为了简化问题，考虑第个网格点在轴方向振动时，考虑第个网格点在轴方向振动时， ，，考虑第个网格点在轴方向振动时，；为了便于书写，将、和简写为、和；
在网格点邻域内建立的近似模型为：
公式1
式中，为在机体轴系网格点的邻域内轴向振动函数，、为轴向振动方程的结构系数函数，分别为在机体轴系网格点，处轴向振动时对应于网格点的变化值；为在机体轴系网格点的邻域内轴向振动函数，、为轴向振动方程的结构系数函数，分别为在机体轴系网格点，处轴向振动时对应于网格点的变化值；为在机体轴系网格点的邻域内轴向振动函数，、为轴向振动方程的结构系数函数，分别为在机体轴系网格点处轴向振动时对应于网格点的变化值；为在网格点的等效激励函数，为时间；为参数向量，表示网格点的温度，为飞行高度，为马赫数，为网格点的气流环境影响，为大气密度；
步骤2：对应步骤1的机体轴系网格点，安装微型温度传感器，、、三个轴向的气流和位置以及振动传感器，在机翼上下方和所有舵面两边安装微型、、三个轴向的气流和位置以及振动传感器，同时在机身加装大于1000帧/秒的图像传感器记录观测机翼翼尖、所有舵面的振动幅值和频率；飞机机载传感器记录时间、飞行高度、马赫数，大气密度；
步骤3：将飞行器到达给定高度和马赫数后颤振试验的过程表达成有效颤振飞行试验，有效颤振飞行试验数据采样时间为， 为正整数， 为记录数据的采样周期，为有效颤振飞行试验的总采样次数；通过颤振飞行试验获得机体轴系网格点，在采样时间时刻的测量值、、和测量值；
步骤4：根据机体轴系网格点，安装微型、、轴向气流传感器，在机翼上下方和所有舵面两边安装微型、、轴向气流传感器，确定时刻机体轴系网格点，的激励函数；
对、、分别采用给定函数逼近，得到：
且关于连续可导，关于连续可导，关于连续可导；这样，可得：
以及;
步骤5： 令：
以及
可将公式1描述成：
公式2
令、、，
式中：
，，对应的埃米特展开系数；
、、为对应于的埃米特展开的阶次；
为阶埃米特正交多项式的递推形式，可得
式中，
其余;
令
，
式中：
、,，、,,、,为对应的埃米特级数的系数，
可得
或写成
公式3
以公式3第一项为例，对
两边求偏导数，可得
根据步骤3和步骤4得到的、、和，和的测试值，可得：
公式4
式中，
进而可得：，代入
可以按照下式及最小二乘估计得出
公式5。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810172987.3A CN108446462B (zh) | 2018-03-02 | 2018-03-02 | 飞行器颤振分析网格模型埃米特建模方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810172987.3A CN108446462B (zh) | 2018-03-02 | 2018-03-02 | 飞行器颤振分析网格模型埃米特建模方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108446462A CN108446462A (zh) | 2018-08-24 |
| CN108446462B true CN108446462B (zh) | 2019-11-19 |
Family
ID=63192870
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810172987.3A Active CN108446462B (zh) | 2018-03-02 | 2018-03-02 | 飞行器颤振分析网格模型埃米特建模方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108446462B (zh) |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2235172A1 (en) * | 2007-12-12 | 2010-10-06 | Aarhus Universitet (University Of Aarhus) | Crystal structure of a type iic p-type atpase |
| JP5237192B2 (ja) * | 2009-05-19 | 2013-07-17 | 川崎重工業株式会社 | 風洞用飛散物回収装置 |
| CN102103036A (zh) * | 2009-12-18 | 2011-06-22 | 西安费斯达自动化工程有限公司 | 基于振动和音频信息的汽车发动机故障诊断系统 |
| JP6303260B2 (ja) * | 2012-12-06 | 2018-04-04 | 株式会社村田製作所 | 正極活物質およびその製造方法、正極、電池、電池パック、電子機器、電動車両、蓄電装置ならびに電力システム |
| CN103809448B (zh) * | 2014-02-28 | 2016-05-25 | 西安费斯达自动化工程有限公司 | 飞行器多回路模型簇颤振抑制复合根轨迹补偿鲁棒控制器设计方法 |
| CN107491616B (zh) * | 2017-08-24 | 2020-09-18 | 北京航空航天大学 | 一种适用于格栅构型舵面的结构有限元参数化建模方法 |
-
2018
- 2018-03-02 CN CN201810172987.3A patent/CN108446462B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN108446462A (zh) | 2018-08-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| González-Rocha et al. | Sensing wind from quadrotor motion | |
| Foster et al. | Dynamics modeling and simulation of large transport airplanes in upset conditions | |
| CN108132134A (zh) | 基于风洞自由飞试验的气动导数辨识方法和系统 | |
| CN109238636A (zh) | 一种级间分离风洞自由飞试验模型分离轨迹预估方法 | |
| CN103471803A (zh) | 一种模型自由飞试验的气动参数确定方法 | |
| RU2615220C1 (ru) | Способ определения управляющего сигнала по углу крена модели гиперзвукового летательного аппарата (ГЛА) для контроля аэродинамической идентичности по числам Рейнольдса траекторий полёта модели и натурного изделия при проведении опережающих лётных исследований | |
| KR101842966B1 (ko) | 자기력을 활용한 비행체의 동안정 미계수 추출방법 | |
| Luo et al. | Simulate the aerodynamic olfactory effects of gas-sensitive UAVs: A numerical model and its parallel implementation | |
| CN108318203A (zh) | 飞行器等价地面颤振综合试验方法 | |
| CN102364477A (zh) | 一种无附加气动阻尼的飞行器颤振特性分析方法 | |
| CN108333945B (zh) | 飞机颤振试验分布式充分激励输入信号设计方法 | |
| CN108388741B (zh) | 飞行器颤振分析网格模型沃尔什建模方法 | |
| CN108446461A (zh) | 飞行器颤振分析网格模型建模方法 | |
| CN108446462B (zh) | 飞行器颤振分析网格模型埃米特建模方法 | |
| CN108399297B (zh) | 飞行器颤振分析网格模型勒让德建模方法 | |
| CN108398231B (zh) | 飞行器颤振分析网格模型Hartley建模方法 | |
| CN108387360B (zh) | 飞行器颤振分析网格模型切比雪夫建模方法 | |
| CN108303897A (zh) | 飞行器颤振分析网格模型拉盖尔建模方法 | |
| Reed et al. | Coaxial Rotor Wake and Prop Induction Impact on a Horizontal Tail Using HPCMP CREATE™-AV Helios | |
| CN108341073A (zh) | 飞机颤振高温试验图像测试方法 | |
| CN108387359B (zh) | 飞行器颤振分析网格模型傅里埃建模方法 | |
| Perez Perez | Forward flight rotor performance at martian atmospheric densities and sensitivity to low Reynolds numbers | |
| Lin et al. | Wake Vortex Encounter Modeling and Simulation for Small Fixed-Wing UAS with Inner Loop Attitude Controller | |
| Weishaeupl et al. | Flight Test Driven Development of Low Cost UAVs-Pitfalls and Opportunities | |
| Sundresan et al. | Review of aeroelasticity testing technology |
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 |