[发明专利]一种基于POD的非定常洞壁干扰修正方法

权利要求书

1.一种基于POD的非定常洞壁干扰修正方法，其特征在于，包括以下步骤：

步骤1，设定风洞试验条件，包括：缩减频率k、雷诺数Re、平均迎角α₀和振幅A；

步骤2，选取N种尺寸的试验模型，分别为：试验模型model₁,model₂,...,model_N，其模型尺寸，即弦长分别为：c₁,c₂,...,c_N；试验模型model₁,model₂,...,model_N之间为等比例放大缩小调节后的试验模型；

步骤3，试验模型model₁,model₂,...,model_N在相同的步骤1设定的风洞实验条件下，分别进行风洞实验，对应得到气动力系数Cx随迎角α变化的实验条件下的气动力系数变化曲线，分别表示为：Cx₁(α),Cx₂(α),...,Cx_N(α)；

步骤4，分别计算实验条件下的气动力系数变化曲线Cx₁(α),Cx₂(α),...,Cx_N(α)的气动力线斜率，对应表示为：

对进行外插计算，得到弦长c＝0时的气动力线斜率k_c＝0；

步骤5，对实验条件下的气动力系数变化曲线Cx₁(α),Cx₂(α),...,Cx_N(α)进行POD分析，获得主模态下的气动力系数变化曲线

步骤6，计算主模态下的气动力系数变化曲线的气动力线斜率k_m；

步骤7，根据下式，使用主模态对气动力系数进行重构，得到主模态重构后的气动力系数变化曲线Cx(c)：

其中：

a＝k_c＝0/k_m

a表示重构系数；

步骤8，在主模态重构后的气动力系数变化曲线Cx(c)中，定位到线性段；对线性段对应的气动力系数求均值，得到与主模态对应的线性段气动力系数均值

步骤9，对于实验条件下的气动力系数变化曲线Cx₁(α),Cx₂(α),...,Cx_N(α)中的每一种气动力系数变化曲线，表示为Cx_i(α)，i＝1,2,...,N，定位到对应的线性段，对定位到的线性段对应的气动力系数求均值，得到与模型尺寸c_i对应的线性段气动力系数均值

因此，对于实验条件下的气动力系数变化曲线Cx₁(α),Cx₂(α),...,Cx_N(α)，分别得到对应的线性段气动力系数均值

步骤10，对进行外插计算，得到c＝0时的线性段气动力系数均值

步骤11，采用下式，将c＝0时的线性段气动力系数均值与主模态对应的线性段气动力系数均值进行差量计算，得到线性段气动力系数差量△Cx：

步骤12，采用下式，对步骤7得到的主模态重构后的气动力系数变化曲线Cx(c)进行平移操作，平移距离为△Cx：

由此得到c＝0时的无洞壁干扰修正时的气动力系数变化曲线Cx(c＝0)。

2.根据权利要求1所述的基于POD的非定常洞壁干扰修正方法，其特征在于，步骤4中，对进行外插计算，得到c＝0时的气动力线斜率k_c＝0，具体为：

以模型尺寸c为横坐标，以气动力线斜率k为纵坐标，建立直角坐标系；在直角坐标系中，对进行拟合，得到气动力线斜率拟合曲线；

计算气动力线斜率拟合曲线中，c＝0时的气动力线斜率，即得到k_c＝0。