基于平均集合卡尔曼滤波的叶栅风洞实验工况参数修正方法

本发明属于航空航天实验测量领域，涉及一种基于流场数值模拟和平均集合卡尔曼滤波的叶栅风洞实验工况参数修正方法。

背景技术：

1、叶栅风洞实验是测量叶轮机械叶型性能的常用方法，其工况参数主要指来流攻角和来流马赫数：一般通过转动叶栅安装转盘到对应角度以获得来流攻角；来流马赫数则通过调节风洞气源的阀门或者风机转速，使实验件前部的马赫数传感器显示到给定的马赫数。然而在实际的风洞实验过程中，风洞中的气流受到风洞壁面影响，气流在流经马赫数传感器后存在损失和折转，真实吹到叶栅区域的气流马赫数和攻角并非给定的工况参数。在大攻角、高马赫数时，工况偏差一般更会为明显，严重影响了叶栅风洞实验的效果。

2、现有的叶栅风洞实验工况参数修正方法使用流场数值模拟并利用了叶片表面全部的测量数据，但是在叶片吸力面后50％区域一般会发生流动分离，但通常的湍流模型对分离区预测准确性较差，校正结果会被分离区预测偏差污染，影响工况参数校正结果的准确性。

3、所以目前叶栅风洞实验领域缺乏一种准确、可靠的工况参数修正方法。

技术实现思路

1、本发明所要解决的技术问题是：

2、为了避免现有修正方法的样本数据存在流动分离导致修正错误的问题，本发明提供一种基于流场数值模拟和平均集合卡尔曼滤波的叶栅风洞实验工况参数修正方法。

3、为了解决上述技术问题，本发明采用的技术方案为：

4、一种基于流场数值模拟和平均集合卡尔曼滤波的叶栅风洞实验工况参数修正方法，其特征在于，包括：

5、步骤1：浮动实验工况参数，构建多组来流工况参数样本，对每组工况参数进行流场数值模拟；

6、步骤2：将浮动工况参数及数值流场的叶片弦长前50％区域的测点位置的相应物理量组成状态向量；

7、步骤3：计算样本协方差矩阵，生成子样本噪声向量及噪声协方差矩阵；

8、步骤4：计算子样本的卡尔曼滤波矩阵；

9、步骤5：对每个子样本进行卡尔曼滤波；

10、步骤6：对所有样本的子样本的状态向量进行平均，从平均后的状态向量中提取修正后的工况参数。

11、本发明进一步的技术方案：步骤1具体如下：根据来流参数的调控精度设置浮动范围，采用拉丁超立方方法进行随机抽样，构建来流参数的样本，在湍流模型对每个样本进行叶栅流场数值计算。

12、本发明进一步的技术方案：其特征在于，步骤2具体如下：

13、将浮动的来流参数及其数值流场结果组建样本状态向量xf,i：

14、

15、式中，main,i和αi分别为来流参数浮动后第i个样本的来流马赫数和来流攻角，表示第i个样本的流场中位于实验中前50％叶片弦长区域第j个测点位置的对应被测物理量，j＝1,2,...,m，m表示验中前50％叶片弦长区域测点的数目。

16、本发明进一步的技术方案：所述的被测物理量包括表面等熵马赫数、静压系数、静压。

17、本发明进一步的技术方案：步骤3具体如下：

18、步骤3-1：计算样本状态向量平均值

19、

20、步骤3-2：计算样本协方差矩阵p：

21、

22、

23、步骤3-3：生成噪声向量wi,l和噪声矩阵wl：其中l＝1,…,l，表示第l个子样本，l为子样本的个数，噪声向量的长度于测点数目相同，即为一组零均值的高斯分布的随机数，其分布方差可参照测量误差，将噪声向量组合成噪声矩阵：

24、wl＝[w1,l,w2,l,…,wn,l]

25、步骤3-4：计算子样本的噪声协方差矩阵：

26、

27、本发明进一步的技术方案：步骤4所述的卡尔曼滤波矩阵计算式如下：

28、kl＝pht(hpht+rl)-1

29、式中，h为观测矩阵：

30、h=[0m×2,im×m]

31、用于将状态向量映射为相应的实验测量值，其中0m×2表示m×2维的零矩阵，im×m表示m×m维的单位矩阵。

32、本发明进一步的技术方案：步骤5具体如下：

33、对集合样本的子样本进行卡尔曼滤波，得到每个子样本校正后结果xa,i,l：

34、xa,i,l＝xf,i+kl(yexp-hxf,i+wi,l)

35、式中，yexp为叶栅风洞实验中前50％叶片弦长区域测得的表面等熵马赫数组成得列向量。

36、本发明进一步的技术方案：步骤6具体如下：

37、对所有滤波后样本进行平均，获得平均后的状态向量

38、

39、向量前两个数据即为修正后来流马赫数和来流攻角。

40、一种计算机系统，其特征在于包括：一个或多个处理器，计算机可读存储介质，用于存储一个或多个程序，其中，当所述一个或多个程序被所述一个或多个处理器执行时，使得所述一个或多个处理器实现上述的方法。

41、一种计算机可读存储介质，其特征在于存储有计算机可执行指令，所述指令在被执行时用于实现上述的方法。

42、本发明的有益效果在于：

43、本发明提供的一种基于流场数值模拟和平均集合卡尔曼滤波的叶栅风洞实验工况参数修正方法，利用流场数值模拟和平均集合卡尔曼滤波方法修正了叶栅风洞的来流马赫数和攻角，修正结果精度高，即便对偏差较大的参数也能够较好的修正，计算量小；本发明提出的平均集合卡尔曼滤波方法相比于已有的集合卡尔曼滤波方法，更简单，无需迭代，鲁棒性更好。具体如下：

44、1、因为在叶片后50％区域一般会发生流动分离，目前的湍流模型预测流动分离都不十分准确，故后50％的数据被湍流模型预测偏差污染，会带来成修正错误。因此，本发明采用前50％弦长区域的数据而没有采用全部数据，而前50％区域是没有流动分离的附着流动，湍流模型一般能够非常准确地预测。

45、2、本发明采用平均集合卡尔曼滤波方法进行数据同化，提出了子样本的概念，即生成合成噪声的时候，对每个样本生成若干个噪声向量，对每个样本构成了若干个子样本，对每个子样本都进行滤波运算，再进行总体平均。相对于现有技术中的集合卡尔曼滤波方法鲁棒性更好。

46、3、采用平均集合卡尔曼滤波方法，计算结果稳定随机性小，鲁棒性好，结果准确，其需要的样本少，计算量少。在集合卡尔曼滤波方法中有两部分随机性，一个是样本的生成，一个是随机噪声，传统方法仅进行一次平均，无法完全消除随机性，为了结果稳定，只能使用较多的样本，平均集合卡尔曼滤波使用了子样本，进行了两个尺度的平均，更好地消除了结果的随机性，所以稳定性更好，较少的样本就能获得的稳定的结果。

技术特征：

1.一种基于流场数值模拟和平均集合卡尔曼滤波的叶栅风洞实验工况参数修正方法，其特征在于，包括：

2.根据权利要求1所述的一种基于流场数值模拟和平均集合卡尔曼滤波的叶栅风洞实验工况参数修正方法，其特征在于，步骤1具体如下：根据来流参数的调控精度设置浮动范围，采用拉丁超立方方法进行随机抽样，构建来流参数的样本，在湍流模型对每个样本进行叶栅流场数值计算。

3.根据权利要求2所述的一种基于流场数值模拟和平均集合卡尔曼滤波的叶栅风洞实验工况参数修正方法，其特征在于，步骤2具体如下：

4.根据权利要求3所述的一种基于流场数值模拟和平均集合卡尔曼滤波的叶栅风洞实验工况参数修正方法，其特征在于，所述的被测物理量包括表面等熵马赫数、静压系数、静压。

5.根据权利要求4所述的一种基于流场数值模拟和平均集合卡尔曼滤波的叶栅风洞实验工况参数修正方法，其特征在于，步骤3具体如下：

6.根据权利要求5所述的一种基于流场数值模拟和平均集合卡尔曼滤波的叶栅风洞实验工况参数修正方法，其特征在于，步骤4所述的卡尔曼滤波矩阵计算式如下：

7.根据权利要求6所述的一种基于流场数值模拟和平均集合卡尔曼滤波的叶栅风洞实验工况参数修正方法，其特征在于，步骤5具体如下：

8.根据权利要求7所述的一种基于流场数值模拟和平均集合卡尔曼滤波的叶栅风洞实验工况参数修正方法，其特征在于，步骤6具体如下：

9.一种计算机系统，其特征在于包括：一个或多个处理器，计算机可读存储介质，用于存储一个或多个程序，其中，当所述一个或多个程序被所述一个或多个处理器执行时，使得所述一个或多个处理器实现权利要求1所述的方法。

10.一种计算机可读存储介质，其特征在于存储有计算机可执行指令，所述指令在被执行时用于实现权利要求1所述的方法。

技术总结
本发明涉及一种平均集合卡尔曼滤波的叶栅风洞实验工况参数修正方法，属于航空航天实验测量域。本发明的基于流场数值模拟和平均集合卡尔曼滤波的叶栅风洞实验工况参数修正方法主要步骤包括：对欲修正实验工况的来流马赫数和攻角进行浮动，构建多组来流参数样本并进行叶栅流场数值计算；使用平均集合卡尔曼滤波方法，修正来实验测量时的流马赫数和来流攻角。本发明能够有效修正叶栅风洞实验的来流边界条件，能有效克服湍流模型产生的分离区预测偏差，修正能力强，无需特意选取湍流模型，计算量小，无需迭代，特别适合于叶栅风洞实验参数的快速修正。

技术研发人员：高丽敏,刘锬韬,但玥,李泽锐
受保护的技术使用者：西北工业大学
技术研发日：
技术公布日：2024/1/13