首页 / 专利库 / 飞机类型 / 飞行器 / 微型飞行器 / 飞行器颤振分析网格模型傅里埃建模方法

飞行器颤振分析网格模型傅里埃建模方法

阅读:510发布:2020-05-14

专利汇可以提供飞行器颤振分析网格模型傅里埃建模方法专利检索,专利查询,专利分析的服务。并且为了克服 现有技术 不能有效表达 气动 力 和强度变化影响下复杂 颤振 模型的问题,本 发明 提供了一种 飞行器 颤振分析 网格模型 傅里埃建模方法,该方法在飞行器 机体 轴系选择多个网格点,在不同飞行速度、大气 密度 、气流环境、不同 温度 等气动力和强度变化影响下按照机体轴系分解方法表示复杂颤振网格模型,根据建立该模型的要求提出安装 传感器 和数据、图像记录要求,通过有效颤振飞行试验获取数据,通过气流传感器测量值获得激励函数,采用傅里埃函数对振动变量进行逼近和等效描述,按照辨识方法同时确定了机体轴系坐标网格点处三个轴向振动方程求解,解决了现有技术不能有效表达气动力和强度变化影响下复杂颤振模型的技术问题。,下面是飞行器颤振分析网格模型傅里埃建模方法专利的具体信息内容。

1.一种飞行器颤振分析网格模型傅里埃建模方法,其特征包括以下步骤:
步骤1:以飞行器机体轴系 分析复杂颤振模型,在机体轴系选取 个网格点,坐标为: , 振动时第 个网格点坐标
为时间 和其它两轴位置的函数,为了便于表达第
个网格点在 轴的振动分量,以 为例,下标 为网格点标号,下
标第二个字母 分别表示振动在机体轴系 的三个轴分量;为了简化问题,考虑第 个网格点在 轴方向振动时,
考虑第 个网
格点在 轴方向振动时, ,
,考虑第 个网格点在 轴方向振动时,
;为了便于书写,将 、
和 简写为 、 和 ;
在网格点邻域内建立的近似模型为:
式中,
为在机体轴系网格点 的邻域内 轴向振动函
数, 、 为 轴向振动方程的
结构系数函数, 分别为在机体轴系网格点 ,
处 轴向振动时对应于网格点 的变化值; 为在
机体轴系网格点 的邻域内 轴向振动函数, 、
为 轴向振动方程的结构系数函数,
分别为在机体轴系网格点 ,处 轴向振动时
对应于网格点 的变化值; 为在机体轴系网格点
的邻域内 轴向振动函数, 、
为 轴向振动方程的结构系数函数,
分别为在机体轴系网格点 处 轴向振动时对应于网格点
的变化值; 为在网格点 的等效激励函数, 为时间;
为参数向量, 表示网格点 的温度, 为飞
行高度, 为赫数, 为网格点 的气流环境影响, 为大气密度
步骤2:对应步骤1的机体轴系网格点 ,安装微型温度传感器
、 、 三个轴向的气流和位置以及振动传感器,在机翼上下方和所有面两边安装微型 、 、 三个轴向的气流和位置以及振动传感器,同时在机身加装大于1000/秒的图像传感器记录观测机翼翼尖、所有舵面的振动幅值和频率;飞机机载传感器记录时间、飞行高度、马赫数,大气密度;
步骤3:将飞行器到达给定高度和马赫数后颤振试验的过程表达成有效颤振飞行试验,有效颤振飞行试验数据采样时间为 ,  为正整数,  为记录数
据的采样周期, 为有效颤振飞行试验的总采样次数;通过颤振飞行试验获得机体轴系网格点 ,在采样时间 时刻的测量值
、 、 和 测量值;
步骤4:根据机体轴系网格点 ,安装微型 、 、 轴向气流传
感器,在机翼上下方和所有舵面两边安装微型 、 、 轴向气流传感器,确定时刻机体轴系网格点 ,的激励函数

对 、 、
分别采用给定函数逼近,得到:
且 关于 连续可导,
关于 连续可导, 关于 连续可导;这样,可得:
以及 ;
步骤5: 令:
以及 
可将(1)式描述成:
  (2)
令  、


式中: 、 、 为常数项;
, , ,对应
的傅里埃展开系数;

, 为频率, 、 、
为对应于 的傅里埃展开的阶次;
式中:

令 
式中: 为常数,
、 、 、 , , 、 、 、 , ,
、 、 、 , ,为对应的傅里埃级数的系数,
可得:
(3)
以(3)式第一项为例,令 ,

  (4)
两边求 偏导数,可得
根据步骤3和步骤4得到的 、 、 和 ,
和 的测试值,可得:
  (5)
式中: ,


进而
可得: ,带入(4)式中,得
进一步可写成:
  (6)
式中:
将(6)式写成: ,有
可以按照最小二乘估计出 。

说明书全文

飞行器颤振分析网格模型傅里埃建模方法

技术领域

[0001] 本发明涉及民用飞机、战斗机、无人机等飞行器飞行安全地面综合试验方法,特别涉及飞行器颤振分析网格模型傅里埃建模方法,属于航空航天与信息技术领域。

背景技术

[0002] 颤振是弹性结构在均匀气流中受到空气动、弹性力和惯性力的耦合作用而发生的一种大幅度振动现象。对于飞机而言,在飞行中受到不确定扰动后会发生振动。此时,由于气流的作用,飞机的弹性结构如机翼、尾翼或操纵面将会产生附加气动力;作为一种激振力,附加气动力将加剧结构的振动。同时空气对飞机结构的阻尼力又试图减弱振动;在低速飞行时,由于阻尼力占优,扰动后的振动逐渐消失;当达到某一飞行速度即颤振临界速度颤振边界后,激振力占优,平衡位置失稳,将产生大幅度振动,导致飞机在数秒内解体,酿成灾难性后果;可以说,从航空工业起步的那一天起,颤振就一直是航空界研究的热问题。
[0003] 为避免颤振事故发生,新机研制必须经历颤振试验环节,以确定不发生飞行颤振的稳定飞行包线;开展颤振问题研究主要有两类途径,一是数值计算:这需要对分析对象进行数学建模,此过程需要在结构、气动等方面引入一定的假设,难以考虑真实存在的各种非线性因素和建模误差的影响,分析结果具有一定的参考价值,但可能与实际情况有较大的偏差;二是试验手段:与颤振有关的试验主要有洞试验和飞行试验。风洞试验可以考虑气动力影响,但此方法要求将试验对象进行缩比设计,缩比模型与真实结构存在一定的差别,且由于风洞洞壁与支架的干扰气动力难免失真;此外对于高速、热环境等情况,风洞试验模拟费用昂贵且实施困难。飞行试验可以完全模拟试验对象的真实工作环境,但试验的条件受限、费用高且风险大,飞机一旦在空中发生颤振,会在几秒甚至更短的时间内解体,飞行员几乎没有处置时间,逃脱概率基本为零。
[0004] 地面颤振模拟试验就是一种可以有效弥补传统试验不足的、极具生命力的颤振研究方法。地面试验以飞行器地面颤振试验系统为研究对象,以多学科设计优化理论研究为核心,密切结合飞行器地面颤振试验系统的工程特点,突破等效试验建模方法、多点分布式气动力建模与控制、颤振试验一体化检测方法等关键技术,着力解决飞行器颤振气动力模型难实现、多点激振力无法精确控制、颤振试验结果无法反复回放等问题,提高总体设计平。
[0005] 航空界、力学界虽然较早对避免颤振的问题进行了研究,但目前的研究还是初级阶段,没有形成一个系统的理论方法体系;现有的方法缺乏飞行器等价地面颤振试验方法和评价;现有技术方法难以描述飞行器在不同飞行速度、大气密度、气流环境、不同温度等气动力和强度变化影响下的复杂颤振模型,使得颤振地面试验研究难以有工程化进展。

发明内容

[0006] 为了克服现有技术不能有效表达气动力和强度变化影响下复杂颤振模型的问题,本发明提供了一种飞行器颤振分析网格模型傅里埃建模方法,该方法在飞行器机体轴系选择多个网格点,在不同飞行速度、大气密度、气流环境、不同温度等气动力和强度变化影响下按照机体轴系分解方法表示复杂颤振网格模型,根据建立该模型的要求提出安装传感器和数据、图像记录要求,通过有效颤振飞行试验获取数据,通过气流传感器测量值获得激励函数,采用傅里埃函数对振动变量进行逼近和等效描述,按照辨识方法同时确定了机体轴系网格点处三个轴向振动方程求解,解决了现有技术不能有效表达气动力和强度变化影响下复杂颤振模型的技术问题。
[0007] 本发明解决其技术问题采用的技术方案是,一种飞行器颤振分析网格模型傅里埃建模方法,其特征包括以下步骤:
[0008] 步骤1:以飞行器机体轴系 分析复杂颤振模型,在机体轴系选取 个网格点,坐标为: ,  振动时第 个网格点坐标为时间 和其它两轴位置的函数,为了便于表达第
个网格点在 轴的振动分量,以 为例,下标 为网格点标号,下
标第二个字母 分别表示振动在机体轴系 的三个轴分量;为了简化问题,考虑第 个网格点在 轴方向振动时,
考虑第 个网
格点在 轴方向振动时, ,
,考虑第 个网格点在 轴方向振动时,
;为了便于书写,将 、
和 简写为 、 和 ;
[0009] 在网格点邻域内建立的近似模型为:
[0010]式中, 为在机体轴系网格点 的邻域内 轴向振
动函数, 、 为 轴向振动方
程的结构系数函数, 分别为在机体轴系网格点
,处 轴向振动时对应于网格点 的变化值;
为在机体轴系网格点 的邻域内 轴向振动函数,
、 为 轴向振动方程的结
构系数函数, 分别为在机体轴系网格点 ,处
轴向振动时对应于网格点 的变化值; 为在机体
轴系网格点 的邻域内 轴向振动函数, 、
为 轴向振动方程的结构系数函数,
分别为在机体轴系网格点 处 轴向振动时对应于网格点
的变化值; 为在网格点 的等效激励函数, 为时间;
为参数向量, 表示网格点 的温度, 为飞
行高度, 为赫数, 为网格点 的气流环境影响, 为大气密度;
[0011] 步骤2:对应步骤1的机体轴系网格点 ,安装微型温度传感器, 、 、 三个轴向的气流和位置以及振动传感器,在机翼上下方和所有面两边安装微型 、 、 三个轴向的气流和位置以及振动传感器,同时在机身加装大于1000/秒的图像传感器记录观测机翼翼尖、所有舵面的振动幅值和频率;飞机机载传感器记录时间、飞行高度、马赫数,大气密度;
[0012] 步骤3:将飞行器到达给定高度和马赫数后颤振试验的过程表达成有效颤振飞行试验,有效颤振飞行试验数据采样时间为 ,  为正整数,  为记录数据的采样周期, 为有效颤振飞行试验的总采样次数;通过颤振飞行试验获得机体轴系网格点 ,在采样时间 时刻的测量值
、 、 和 测量值;
[0013] 步骤4:根据机体轴系网格点 ,安装微型 、 、 轴向气流传感器,在机翼上下方和所有舵面两边安装微型 、 、 轴向气流传感器,确定时刻机体轴系网格点 ,的激励函数

[0014] 对 、 、分别采用给定函数逼近,得到:
[0015]
[0016] 且 关于 连续可导,关于 连续可导, 关于
连续可导;这样,可得:
[0017]
[0018] 以及 ;
[0019] 步骤5: 令:
[0020]
[0021] 以及 
[0022] 可将(1)式描述成:
[0023]   (2)
[0024] 令 、
[0025] 、,
[0026] 式中: 、 、 为常数项;
[0027]
[0028] , , ,对应的傅里埃展开系数;
[0029]
[0030]
[0031]
[0032]
[0033]
[0034] ,
[0035] , 为频率, 、 、为对应于 的傅里埃展开的阶次;
[0036]
[0037]
[0038] 式中:
[0039] ;
[0040] 令 
[0041] 式中: 为常数,
[0042]
[0043]
[0044] 、 、 、 , , 、 、 、 ,, 、 、 、 , ,为对应的傅里埃级数的系数,
[0045] 可得:
[0046] (3)
[0047] 以(3)式第一项为例,令,对
[0048]   (4)
[0049] 两边求 偏导数,可得
[0050]
[0051] 根据步骤3和步骤4得到的 、 、 和 ,和 的测试值,可得:
[0052]   (5)
[0053] 式中: ,
[0054] ,,
[0055]进而可得: ,带入(4)式中,得
[0056]
[0057] 进一步可写成:
[0058]   (6)
[0059] 式中:
[0060] 将(6)式写成: ,有
[0061]
[0062] 可以按照最小二乘估计出 。
[0063] 本发明的有益结果是:在飞行器机体轴系选择多个网格点,考虑不同飞行速度、大气密度、气流环境、不同温度等气动力和强度变化影响下按照机体轴系分解方法表示复杂颤振网格模型,根据建立该模型的要求提出安装传感器和数据、图像记录要求,通过有效颤振飞行试验获取数据,通过气流传感器测量值获得激励函数,通过气流传感器测量值获得激励函数,采用傅里埃函数对振动变量进行逼近和等效描述,按照辨识方法同时确定了机体轴系网格点处三个轴向振动方程求解,从而给出了完整的复杂颤振模型网格模型建模技术方案,解决了现有技术不能有效表达气动力和强度变化影响下复杂颤振模型的技术问题。
[0064] 下面结合具体实例对本发明作详细说明。

具体实施方式

[0065] 步骤1:以飞行器机体轴系 分析复杂颤振模型,在机体轴系选取 个网格点,坐标为: ,  振动时第 个网格点坐标为时间 和其它两轴位置的函数,为了便于表达第
个网格点在 轴的振动分量,以 为例,下标 为网格点标号,下
标第二个字母 分别表示振动在机体轴系 的三个轴分量;为了简化问题,考虑第 个网格点在 轴方向振动时,
考虑第 个网
格点在 轴方向振动时, ,
,考虑第 个网格点在 轴方向振动时,
;为了便于书写,将 、
和 简写为 、 和 ;
[0066] 在网格点邻域内建立的近似模型为:
[0067]式中, 为在机体轴系网格点 的邻域内 轴向振
动函数, 、 为 轴向振动方
程的结构系数函数, 分别为在机体轴系网格点
,处 轴向振动时对应于网格点 的变化值;
为在机体轴系网格点 的邻域内 轴向振动函数,
、 为 轴向振动方程的结
构系数函数, 分别为在机体轴系网格点 ,处
轴向振动时对应于网格点 的变化值; 为在机体
轴系网格点 的邻域内 轴向振动函数, 、
为 轴向振动方程的结构系数函数,
分别为在机体轴系网格点 处 轴向振动时对应于网格点
的变化值; 为在网格点 的等效激励函数, 为时间;
为参数向量, 表示网格点 的温度, 为飞
行高度, 为马赫数, 为网格点 的气流环境影响, 为大气密度;
[0068] 步骤2:对应步骤1的机体轴系网格点 ,安装微型温度传感器, 、 、 三个轴向的气流和位置以及振动传感器,在机翼上下方和所有舵面两边安装微型 、 、 三个轴向的气流和位置以及振动传感器,同时在机身加装大于1000帧/秒的图像传感器记录观测机翼翼尖、所有舵面的振动幅值和频率;飞机机载传感器记录时间、飞行高度、马赫数,大气密度;
[0069] 步骤3:将飞行器到达给定高度和马赫数后颤振试验的过程表达成有效颤振飞行试验,有效颤振飞行试验数据采样时间为 ,  为正整数,  为记录数据的采样周期, 为有效颤振飞行试验的总采样次数;通过颤振飞行试验获得机体轴系网格点 ,在采样时间 时刻的测量值
、 、 和 测量值;
[0070] 步骤4:根据机体轴系网格点 ,安装微型 、 、 轴向气流传感器,在机翼上下方和所有舵面两边安装微型 、 、 轴向气流传感器,确定时刻机体轴系网格点 ,的激励函数

[0071] 对 、 、分别采用给定函数逼近,得到:
[0072]
[0073] 且 关于 连续可导,关于 连续可导, 关于
连续可导;这样,可得:
[0074]
[0075] 以及 ;
[0076] 步骤5: 令: 
[0077] 以及 
[0078] 可将(1)式描述成:
[0079]   (2)
[0080] 令 、
[0081] 、,
[0082] 式中: 、 、 为常数项;
[0083]
[0084] , , ,对应的傅里埃展开系数;
[0085]
[0086]
[0087]
[0088]
[0089]
[0090] ,
[0091] , 为角频率, 、 、为对应于 的傅里埃展开的阶次;
[0092]
[0093]
[0094] 式中:
[0095] ;
[0096] 令 
[0097] 式中: 为常数,
[0098]
[0099]
[0100] 、 、 、 , , 、 、 、 ,, 、 、 、 , ,为对应的傅里埃级数的系数,
[0101] 可得:
[0102] (3)
[0103] 以(3)式第一项为例,令 ,对
[0104]   (4)
[0105] 两边求 偏导数,可得
[0106]
[0107] 根据步骤3和步骤4得到的 、 、 和 ,和 的测试值,可得:
[0108]   (5)
[0109] 式中: ,
[0110] , ,
[0111]进而可得: ,带入(4)式中,得
[0112]
[0113] 进一步可写成:
[0114]   (6)
[0115] 式中:
[0116] 将(6)式写成: ,有
[0117]
[0118] 可以按照最小二乘估计出 。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈