飞行器颤振分析网格模型沃尔什建模方法

摘要

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

权利要求

1.一种飞行器颤振分析网格模型沃尔什建模方法，其特征包括以下步骤：步骤1：以飞行器机体轴系 分析复杂颤振模型，在机体轴系选取 个网格点，坐标为: , 振动时第 个网格点坐标为时间 和其它两轴位置的函数，为了便于表达第个网格点在 轴的振动分量，以 为例，下标 为网格点标号，下标第二个字母 分别表示振动在机体轴系 的三个轴分量；为了简化问题，考虑第 个网格点在 轴方向振动时，考虑第 个网

格点在 轴方向振动时， ，，考虑第 个网格点在 轴方向振动时，；为了便于书写，将 、和 简写为 、 和 ；在网格点邻域内建立的近似模型为：公式1

式中， 为在机体轴系网格点 的邻域内 轴向振动函数， 、 为 轴向振动方程的结 构系数函数 ， 分别为在机 体轴系网格 点，处 轴向振动时对应于网格点 的变化值；

为在机体轴系网格点 的邻域内 轴向振动函数，、 为 轴向振动方程的结构系数函数， 分别为在机体轴系网格点 ，处轴向振动时对应于网格点 的变化值； 为在机体轴系网格点 的邻域内 轴向振动函数， 、为 轴向振动方程的结构系数函数，分别为在机体轴系网格点 处 轴向振动时对应于网格点的变化值； 为在网格点 的等效激励函数， 为时间；

为参数向量， 表示网格点 的温度， 为飞行高度， 为马赫数， 为网格点 的气流环境影响， 为大气密度；

步骤2：对应步骤1的机体轴系网格点 ，安装微型温度传感器，、 、 三个轴向的气流和位置以及振动传感器，在机翼上下方和所有舵面两边安装微型 、 、 三个轴向的气流和位置以及振动传感器，同时在机身加装大于1000帧/秒的图像传感器记录观测机翼翼尖、所有舵面的振动幅值和频率；飞机机载传感器记录时间、飞行高度、马赫数，大气密度；

步骤3：将飞行器到达给定高度和马赫数后颤振试验的过程表达成有效颤振飞行试验，有效颤振飞行试验数据采样时间为 ， 为正整数， 为记录数据的采样周期， 为有效颤振飞行试验的总采样次数；通过颤振飞行试验获得机体轴系网格点 ，在采样时间 时刻的测量值、 、 和 测量值；

步骤4：根据机体轴系网格点 ，安装微型 、 、 轴向气流传感器，在机翼上下方和所有舵面两边安装微型 、 、 轴向气流传感器，确定时刻机体轴系网格点 ，的激励函数；

对 、 、

分别采用给定函数逼近，得到：且 关于 连续可导，关于 连续可导， 关于 连续可导；这样，可得：以及 ;

步骤5：令：

以及

可将公式1描述成：

公式2

令 、 、 ，

式中：

， ， 对应的沃尔什展开系数；

、 、 为对应于 的沃尔什展开的阶次；

为沃尔

什函数； ， 为0或 的二进制表示式的二进制数值，为二进制值位数，表示符号函数， 表示 绝对值的最大值；可得式中：

式中： 、 ， ，分别为 单位矩阵及零矩阵；

令

，

式中：

、 , ， 、 , , 、 , 为对应的沃尔什级数的系数，可得

或写成

公式3

以公式3第一项为例，对两边

求 偏导数，可得

根据步骤3和步骤4得到的 、 、 和 ，和 的测试值，可得：                                   公式4式中，

进 而 可

得： ，代入

可以

按照下式及最小二乘估计得出 公式5。

说明书

飞行器颤振分析网格模型沃尔什建模方法

技术领域

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

背景技术

[0002] 颤振是弹性结构在均匀气流中受到空气动力、弹性力和惯性力的耦合作用而发生的一种大幅度振动现象。对于飞机而言，在飞行中受到不确定扰动后会发生振动。此时，由于气流的作用，飞机的弹性结构如机翼、尾翼或操纵面将会产生附加气动力；作为一种激振力，附加气动力将加剧结构的振动。同时空气对飞机结构的阻尼力又试图减弱振动；在低速飞行时，由于阻尼力占优，扰动后的振动逐渐消失；当达到某一飞行速度即颤振临界速度颤振边界后，激振力占优，平衡位置失稳，将产生大幅度振动，导致飞机在数秒内解体，酿成灾难性后果；可以说，从航空工业起步的那一天起，颤振就一直是航空界研究的热门问题。

[0003] 为避免颤振事故发生，新机研制必须经历颤振试验环节，以确定不发生飞行颤振的稳定飞行包线；开展颤振问题研究主要有两类途径，一是数值计算：这需要对分析对象进行数学建模，此过程需要在结构、气动等方面引入一定的假设，难以考虑真实存在的各种非线性因素和建模误差的影响，分析结果具有一定的参考价值，但可能与实际情况有较大的偏差；二是试验手段：与颤振有关的试验主要有风洞试验和飞行试验。风洞试验可以考虑气动力影响，但此方法要求将试验对象进行缩比设计，缩比模型与真实结构存在一定的差别，且由于风洞洞壁与支架的干扰气动力难免失真；此外对于高速、热环境等情况，风洞试验模拟费用昂贵且实施困难。飞行试验可以完全模拟试验对象的真实工作环境，但试验的条件受限、费用高且风险大，飞机一旦在空中发生颤振，会在几秒甚至更短的时间内解体，飞行员几乎没有处置时间，逃脱概率基本为零。

[0004] 地面颤振模拟试验就是一种可以有效弥补传统试验不足的、极具生命力的颤振研究方法。地面试验以飞行器地面颤振试验系统为研究对象，以多学科设计优化理论研究为核心，密切结合飞行器地面颤振试验系统的工程特点，突破等效试验建模方法、多点分布式气动力建模与控制、颤振试验一体化检测方法等关键技术，着力解决飞行器颤振气动力模型难实现、多点激振力无法精确控制、颤振试验结果无法反复回放等问题，提高总体设计水平。

[0005] 航空界、力学界虽然较早对避免颤振的问题进行了研究，但目前的研究还是初级阶段，没有形成一个系统的理论方法体系；现有的方法缺乏飞行器等价地面颤振试验方法和评价；特别是现有技术方法难以描述飞行器在不同飞行速度、大气密度、气流环境、不同温度等气动力和强度变化影响下的复杂颤振模型，使得颤振地面试验研究难以有工程化进展。

发明内容

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

[0007] 本发明解决其技术问题采用的技术方案是，一种飞行器颤振分析网格模型沃尔什建模方法，其特征包括以下步骤：

[0008] 步骤1：以飞行器机体轴系 分析复杂颤振模型，在机体轴系选取 个网格点，坐标为: ,  振动时 第 个网格点坐标为时间 和其它两轴位置的函数，为了便于表达第
个网格点在 轴的振动分量，以 为例，下标 为网格点标号，下
标第二个字母 分别表示振动在机体轴系 的三个轴分量；为了简化问题，考虑第 个网格点在 轴方向振动时，
考虑第 个网
格点在 轴方向振动时， ，
，考虑第 个网格点在 轴方向振动时，
；为了便于书写，将 、
和 简写为 、 和 ；在网格点邻域内建立的近似模型为：

[0009]

[0010] 公式1

[0011] 式中， 为在机体轴系网格点 的邻域内轴向振动函数， 、 为
轴向振动方程的结构系数函数， 分别为在机体轴系网格点
，处 轴向振动时对应于网格点 的变化值；
为在机体轴系网格点 的邻域内 轴向振动函数，
、 为 轴向振动方程的结
构系数函数， 分别为在机体轴系网格点 ，处
轴向振动时对应于网格点 的变化值； 为在机体
轴系网格点 的邻域内 轴向振动函数， 、
为 轴向振动方程的结构系数函数，
分别为在机体轴系网格点 处 轴向振动时对应于网格点
的变化值； 为在网格点 的等效激励函数， 为时间；
为参数向量， 表示网格点 的温度， 为飞
行高度， 为马赫数， 为网格点 的气流环境影响， 为大气密度；

[0012] 步骤2：对应步骤1的机体轴系网格点 ，安装微型温度传感器， 、 、 三个轴向的气流和位置以及振动传感器，在机翼上下方和所有舵面两边安装微型 、 、 三个轴向的气流和位置以及振动传感器，同时在机身加装大于1000帧/秒的图像传感器记录观测机翼翼尖、所有舵面的振动幅值和频率；飞机机载传感器记录时间、飞行高度、马赫数，大气密度；

[0013] 步骤3：将飞行器到达给定高度和马赫数后颤振试验的过程表达成有效颤振飞行试验，有效颤振飞行试验数据采样时间为 ，  为正整数，  为记录数据的采样周期， 为有效颤振飞行试验的总采样次数；通过颤振飞行试验获得机体轴系网格点 ，在采样时间 时刻的测量值
、 、 和 测量值；

[0014] 步骤4：根据机体轴系网格点 ，安装微型 、 、 轴向气流传感器，在机翼上下方和所有舵面两边安装微型 、 、 轴向气流传感器，确定时刻机体轴系网格点 ，的激励函数
；

[0015] 对 、 、分别采用给定函数逼近，得到：

[0016]

[0017] 且 关于 连续可导，关于 连续可导， 关于
连续可导；这样，可得：

[0018]

[0019] 以及 ;

[0020] 步骤5： 令：

[0021]

[0022] 以及 

[0023] 可将公式1描述成：

[0024] 公式2

[0025] 令 、 、 ，

[0026] 式中：

[0027]

[0028] ， ， 对应的沃尔什展开系数；

[0029]

[0030]

[0031] 、 、 为对应于 的沃尔什展开的阶次；

[0032] 为沃尔什函数； ， 为0或 的二进制表示式的二进制数值，为二进制值位数， 表示符号函数；， 表示 绝对值的最大值；可得

[0033]

[0034]

[0035] 式中：

[0036]

[0037] 式中： 、 ， ，分别为 单位矩阵及零矩阵；

[0038] 令 

[0039] ，

[0040] 式中：

[0041]

[0042]

[0043] 、 , ， 、 , , 、 ,为对应的沃尔什级数的系数，

[0044] 可得

[0045]

[0046] 或写成

[0047]

[0048] 公式3

[0049] 以公式3第一项为例，对

[0050]两边求 偏导数，可得

[0051] 根据步骤3和步骤4得到的 、 、 和 ，和 的测试值，可得：

[0052]                                    公式4

[0053] 式中，进 而 可
得： ，代入

[0054]可以按照下式及最小二乘估计得出

[0055]  公式5。

[0056] 本发明的有益结果是：在飞行器机体轴系选择多个网格点，考虑不同飞行速度、大气密度、气流环境、不同温度等气动力和强度变化影响下按照机体轴系分解方法表示复杂颤振网格模型，根据建立该模型的要求提出安装传感器和数据、图像记录要求，通过有效颤振飞行试验获取数据，通过气流传感器测量值获得激励函数，通过气流传感器测量值获得激励函数，采用沃尔什函数对振动变量进行逼近和等效描述，按照辨识方法同时确定了机体轴系坐标网格点处三个轴向振动方程求解，从而给出了完整的复杂颤振模型网格模型建模技术方案，解决了现有技术不能有效表达气动力和强度变化影响下复杂颤振模型的技术问题。

[0057] 下面结合具体实例对本发明作详细说明。

具体实施方式

[0058] 步骤1：以飞行器机体轴系 分析复杂颤振模型，在机体轴系选取 个网格点 ，坐 标为: ,  振 动时第 个网 格点坐 标为时间 和其它两轴位置的函数，为了便于表达第
个网格点在 轴的振动分量，以 为例，下标 为网格点标号，下
标第二个字母 分别表示振动在机体轴系 的三个轴分量；为了简化问题，考虑第 个网格点在 轴方向振动时，
考虑第 个网
格点在 轴方向振动时， ，
，考虑第 个网格点在 轴方向振动时，
；为了便于书写，将 、
和 简写为 、 和 ；在网格点邻域内建立的近似模型为：

[0059]

[0060] 公式1

[0061] 式中， 为在机体轴系网格点 的邻域内轴向振动函数， 、 为
轴向振动方程的结构系数函数， 分别为在机体轴系网格点
，处 轴向振动时对应于网格点 的变化值；
为在机体轴系网格点 的邻域内 轴向振动函数，
、 为 轴向振动方程的结
构系数函数， 分别为在机体轴系网格点 ，处
轴向振动时对应于网格点 的变化值； 为在机体
轴系网格点 的邻域内 轴向振动函数， 、
为 轴向振动方程的结构系数函数，
分别为在机体轴系网格点 处 轴向振动时对应于网格点
的变化值； 为在网格点 的等效激励函数， 为时间；
为参数向量， 表示网格点 的温度， 为飞
行高度， 为马赫数， 为网格点 的气流环境影响， 为大气密度；

[0062] 步骤2：对应步骤1的机体轴系网格点 ，安装微型温度传感器， 、 、 三个轴向的气流和位置以及振动传感器，在机翼上下方和所有舵面两边安装微型 、 、 三个轴向的气流和位置以及振动传感器，同时在机身加装大于1000帧/秒的图像传感器记录观测机翼翼尖、所有舵面的振动幅值和频率；飞机机载传感器记录时间、飞行高度、马赫数，大气密度；

[0063] 步骤3：将飞行器到达给定高度和马赫数后颤振试验的过程表达成有效颤振飞行试验，有效颤振飞行试验数据采样时间为 ，  为正整数，  为记录数据的采样周期， 为有效颤振飞行试验的总采样次数；通过颤振飞行试验获得机体轴系网格点 ，在采样时间 时刻的测量值
、 、 和 测量值；

[0064] 步骤4：根据机体轴系网格点 ，安装微型 、 、 轴向气流传感器，在机翼上下方和所有舵面两边安装微型 、 、 轴向气流传感器，确定时刻机体轴系网格点 ，的激励函数
；

[0065] 对 、 、分别采用给定函数逼近，得到：

[0066]

[0067] 且 关于 连续可导，关于 连续可导， 关于
连续可导；这样，可得：

[0068]

[0069] 以及 ;

[0070] 步骤5： 令：

[0071]

[0072] 以及 

[0073] 可将公式1描述成：

[0074] 公式2

[0075] 令 、 、 ，

[0076] 式中：

[0077]

[0078] ， ， 对应的沃尔什展开系数；

[0079]

[0080]

[0081] 、 、 为对应于 的沃尔什展开的阶次；

[0082] 为沃尔什函数； ， 为0或 的二进制表示式的二进制数值，为二进制值位数， 表示符号函数；， 表示 绝对值的最大值；可得

[0083]

[0084]

[0085] 式中：

[0086]

[0087] 式中： 、 ， ，分别为 单位矩阵及零矩阵；

[0088] 令 

[0089] ，

[0090] 式中：

[0091]

[0092]

[0093] 、 , ， 、 , , 、 ,为对应的沃尔什级数的系数，

[0094] 可得

[0095]

[0096] 或写成

[0097]

[0098] 公式3

[0099] 以公式3第一项为例，对

[0100]两边求 偏导数，可得

[0101] 根据步骤3和步骤4得到的 、 、 和 ，和 的测试值，可得：

[0102]                                    公式4

[0103] 式中，进而可得：
，代入

[0104]可以按照下式及最小二乘估计得出

[0105]  公式5。