[发明专利]一种基于主成分分析的翼型网格快速变形方法

权利要求书

1.一种基于主成分分析的翼型网格快速变形方法，其特征在于，具体步骤包括：

步骤1：选择n个不同的翼型，生成对应的流场网格；

步骤2：提取每个翼型流场网格的网格点坐标，将每个网格单元与网格坐标点的映射关系写入矩阵E，将网格点坐标写入矩阵P；

步骤3：采用主成分分析法，计算矩阵P的主成分系数矩阵C、主成分得分矩阵S、解释度向量e和坐标均值向量μ；

步骤4：计算解释度大于99.9％对应的主成分维数m；

步骤5：根据主成分系数矩阵C和坐标均值向量μ，计算初始翼型网格点坐标向量N的前m维主成分得分向量Sm，将Sm作为初始翼型的外形参数；

步骤6：改变Sm各维度的数值，生成新的主成分得分向量Sm’；

步骤7：Sm’与矩阵C前m列的转置相乘，然后与坐标均值向量μ相加，得到变形后的翼型网格坐标向量Nt；

步骤8：根据向量Nt和矩阵E按照特定网格文件格式重新生成新的翼型网格，实现翼型网格变形。

2.如权利要求1所述的基于主成分分析的翼型网格快速变形方法，其特征在于，所述步骤1中，生成n个不同翼型对应的流场网格时，采用批处理方式，保证n个流场网格的网格坐标点具有相同的拓扑结构，即具有相同的矩阵E。

3.如权利要求1所述的基于主成分分析的翼型网格快速变形方法，其特征在于，所述步骤1中，为保证网格变形的精度，n应不小于200。

4.如权利要求1所述的基于主成分分析的翼型网格快速变形方法，其特征在于，所述步骤2中，矩阵E的行向量存储每个网格单元对应的网格节点编号。

5.如权利要求1所述的基于主成分分析的翼型网格快速变形方法，其特征在于，所述步骤2中，将每个翼型的网格点坐标重整成一维行向量，存入矩阵P的行向量中。矩阵P的每行存储一个翼型的所有网格点坐标向量。

6.如权利要求1所述的基于主成分分析的翼型网格快速变形方法，其特征在于，所述步骤5中，Sm的具体计算方法为：Sm＝(N-μ)·Cm。