基于分段精细寻优策略的叶轮多工况气动优化方法及装置

1.本发明涉及离心压气机叶轮的气动设计技术领域，特别是指一种基于分段精细寻优策略的叶轮多工况气动优化方法及装置。


背景技术：

2.离心压气机是保障国防安全和促进国民经济发展的重要动力装备，已广泛应用于航空航天、船舶、化工和新能源等战略需求领域。根据全国能源基础与标准化委员会的有关统计资料，工业压气机系统年耗电量约占全国总发电量的6％-9％左右。随着中国在第七十五届联合国大会上提出“碳达峰、碳中和”的目标承诺及推进，提高离心压气机的气动性能对节能减排具有积极意义。
3.在离心压气机叶轮的复杂应用场景中需同时兼顾多个不同工况的气动综合性能水平，会面临设计变量多、“黑箱”搜索和优化效率低等难题，增加了气动设计优化的难度。


技术实现要素：

4.为了减少设计变量、缩小设计空间、避免盲目搜索、提高计算效率和改善寻优质量，高效求解离心压气机叶轮的最佳几何构型，提出了本发明。
5.为解决上述技术问题，本发明提供如下技术方案：
6.一方面，本发明提供了一种基于分段精细寻优策略的叶轮多工况气动优化方法，该方法由电子设备实现，该方法包括：
7.s1、获取待优化的离心压气机叶轮的几何构型。
8.s2、将待优化的离心压气机叶轮的几何构型输入到构建好的叶轮多工况气动优化模型；其中，叶轮多工况气动优化模型包括全局多工况气动优化模型以及局部多工况气动优化模型。
9.s3、基于待优化的离心压气机叶轮的几何构型、全局多工况气动优化模型以及局部多工况气动优化模型，得到多工况场景下离心压气机叶轮优化后的最佳气动性能及叶片几何外形。
10.可选地，s3中的基于待优化的离心压气机叶轮的几何构型、全局多工况气动优化模型以及局部多工况气动优化模型，得到多工况场景下离心压气机叶轮优化后的最佳气动性能及叶片几何外形包括：
11.s31、对待优化的离心压气机叶轮的几何构型进行全局参数化建模，得到叶轮中叶片吸力面与单位样条曲面的映射模型以及叶片压力面与单位样条曲面的映射模型。
12.s32、设定全局多工况气动优化模型的目标函数和约束条件，基于全局多工况气动优化模型的目标函数和约束条件、叶轮中叶片吸力面与单位样条曲面的映射模型以及叶片压力面与单位样条曲面的映射模型，得到全局优化后的叶片几何构型和气动综合性能。
13.s33、基于全局优化后的叶片几何构型和气动综合性能，得到多个叶片的局部几何区域。
14.s34、基于多个叶片的局部几何区域，对全局优化后的叶片几何构型进行局部参数化建模，建立叶片局部几何构型与空间控制体的映射模型。
15.s35、基于全局优化后的叶片几何构型和气动综合性能，设定局部多工况气动优化模型的目标函数和约束条件，基于局部多工况气动优化模型的目标函数和约束条件以及叶片局部几何构型与空间控制体的映射模型，得到多工况场景下离心压气机叶轮优化后的最佳气动性能及叶片几何外形。
16.可选地，s31中的对待优化的离心压气机叶轮的几何构型进行全局参数化建模，得到叶轮中叶片吸力面与单位样条曲面的映射模型以及叶片压力面与单位样条曲面的映射模型包括：
17.s311、获取待优化的离心压气机叶轮的几何构型中叶片各截面型线。
18.s312、构建单位映射样条曲面。
19.s313、根据叶片各截面型线以及单位映射样条曲面，得到叶片吸力面与单位样条曲面的映射模型以及叶片压力面与单位样条曲面的映射模型。
20.可选地，s312中的单位映射样条曲面为伯恩斯坦bernstein基函数贝塞尔bezier曲面。
21.可选地，s312中的单位映射样条曲面的数学表达式如下式(1)所示：
[0022][0023]
式中，是贝塞尔bezier曲面上的点坐标，p
k,l
是bezier曲面控制顶点，控制点总数为(m+1)
×
(n+1)，m为横轴布局点个数，n为纵轴布局点个数；v和u是两个变化范围为[0,1]的自变量；和是伯恩斯坦bernstein基函数。
[0024]
可选地，s32中的基于全局多工况气动优化模型的目标函数和约束条件、叶轮中叶片吸力面与单位样条曲面的映射模型以及叶片压力面与单位样条曲面的映射模型，得到全局优化后的叶片几何构型和气动综合性能包括：
[0025]
基于全局多工况气动优化模型的目标函数和约束条件、叶轮中叶片吸力面与单位样条曲面的映射模型以及叶片压力面与单位样条曲面的映射模型，采用多目标优化算法，得到全局优化后的叶片几何构型和气动综合性能。
[0026]
可选地，多目标优化算法为元启发式优化算法或超启发式优化算法。
[0027]
可选地，s34中的对全局优化后的叶片几何构型进行局部参数化建模，建立叶片局部几何构型与空间控制体的映射模型包括：
[0028]
基于b样条基函数的自由曲面变形ffd方法，对全局优化后的叶片几何构型进行局部参数化建模，建立叶片局部几何构型与空间控制体的映射模型。
[0029]
可选地，s34中的叶片局部几何构型与空间控制体的映射模型的数学表达式如下式(2)所示：
[0030][0031]
式中，为叶片表面坐标，为ffd控制框架上的控制顶点；(s,t,u)为在控制框架内的局部坐标；i,j,k为ffd控制框架三个方向的标号；l,m,n为ffd控制框架在
三个方向上的划分数；n
i,d
(s)，n
j,e
(t)，n
k,f
(u)分别对应于d,e,f阶的b样条基函数。
[0032]
另一方面，本发明提供了一种基于分段精细寻优策略的叶轮多工况气动优化装置，该装置应用于实现基于分段精细寻优策略的叶轮多工况气动优化方法，该装置包括：
[0033]
获取模块，用于获取待优化的离心压气机叶轮的几何构型。
[0034]
输入模块，用于将待优化的离心压气机叶轮的几何构型输入到构建好的叶轮多工况气动优化模型；其中，叶轮多工况气动优化模型包括全局多工况气动优化模型以及局部多工况气动优化模型。
[0035]
输出模块，用于基于待优化的离心压气机叶轮的几何构型、全局多工况气动优化模型以及局部多工况气动优化模型，得到多工况场景下离心压气机叶轮优化后的最佳气动性能及叶片几何外形。
[0036]
可选地，输出模块，进一步用于：
[0037]
s31、对待优化的离心压气机叶轮的几何构型进行全局参数化建模，得到叶轮中叶片吸力面与单位样条曲面的映射模型以及叶片压力面与单位样条曲面的映射模型。
[0038]
s32、设定全局多工况气动优化模型的目标函数和约束条件，基于全局多工况气动优化模型的目标函数和约束条件、叶轮中叶片吸力面与单位样条曲面的映射模型以及叶片压力面与单位样条曲面的映射模型，得到全局优化后的叶片几何构型和气动综合性能。
[0039]
s33、基于全局优化后的叶片几何构型和气动综合性能，得到气动性能改善潜力较大的多个叶片的局部几何区域。
[0040]
s34、基于多个叶片的局部几何区域，对全局优化后的叶片几何构型进行局部参数化建模，建立叶片局部几何构型与空间控制体的映射模型。
[0041]
s35、基于全局优化后的叶片几何构型和气动综合性能，设定局部多工况气动优化模型的目标函数和约束条件，基于局部多工况气动优化模型的目标函数和约束条件以及叶片局部几何构型与空间控制体的映射模型，得到多工况场景下离心压气机叶轮优化后的最佳气动性能及叶片几何外形。
[0042]
可选地，输出模块，进一步用于：
[0043]
s311、获取待优化的离心压气机叶轮的几何构型中叶片各截面型线。
[0044]
s312、构建单位映射样条曲面。
[0045]
s313、根据叶片各截面型线以及单位映射样条曲面，得到叶片吸力面与单位样条曲面的映射模型以及叶片压力面与单位样条曲面的映射模型。
[0046]
可选地，单位映射样条曲面为伯恩斯坦bernstein基函数贝塞尔bezier曲面。
[0047]
可选地，单位映射样条曲面的数学表达式如下式(1)所示：
[0048][0049]
式中，是贝塞尔bezier曲面上的点坐标，p
k,l
是bezier曲面控制顶点，控制点总数为(m+1)
×
(n+1)，m为横轴布局点个数，n为纵轴布局点个数；v和u是两个变化范围为[0,1]的自变量；和是伯恩斯坦bernstein基函数。
[0050]
可选地，输出模块，进一步用于：
[0051]
基于全局多工况气动优化模型的目标函数和约束条件、叶轮中叶片吸力面与单位样条曲面的映射模型以及叶片压力面与单位样条曲面的映射模型，采用多目标优化算法，
得到全局优化后的叶片几何构型和气动综合性能。
[0052]
可选地，多目标优化算法为元启发式优化算法或超启发式优化算法。
[0053]
可选地，输出模块，进一步用于：
[0054]
基于b样条基函数的自由曲面变形ffd方法，对全局优化后的叶片几何构型进行局部参数化建模，建立叶片局部几何构型与空间控制体的映射模型。
[0055]
可选地，输出模块，进一步用于：
[0056][0057]
式中，为叶片表面坐标，为ffd控制框架上的控制顶点；(s,t,u)为在控制框架内的局部坐标；i,j,k为ffd控制框架三个方向的标号；l,m,n为ffd控制框架在三个方向上的划分数；n
i,d
(s)，n
j,e
(t)，n
k,f
(u)分别对应于d,e,f阶的b样条基函数。
[0058]
一方面，提供了一种电子设备，所述电子设备包括处理器和存储器，所述存储器中存储有至少一条指令，所述至少一条指令由所述处理器加载并执行以实现上述基于分段精细寻优策略的叶轮多工况气动优化方法。
[0059]
一方面，提供了一种计算机可读存储介质，所述存储介质中存储有至少一条指令，所述至少一条指令由处理器加载并执行以实现上述基于分段精细寻优策略的叶轮多工况气动优化方法。
[0060]
本发明实施例提供的技术方案带来的有益效果至少包括：
[0061]
上述方案中，提出了基于分段精细寻优策略的离心叶轮多工况气动优化方法，建立了两个优化阶段：基于全局优化系统的第一阶段和基于局部优化系统的第二阶段，研究了空间复杂曲面参数化方法、多目标优化算法及寻优策略的最佳匹配关系，求解了等熵效率最大值，通过两次寻优实现了复杂曲面叶片精细化高效构型，最大限度每个阶段的开发潜力，实现了寻优过程中减少设计变量、缩小设计空间、加快收敛速度、增强探索能力及形性优化控制的多重目的，具有一定的推广应用价值。
附图说明
[0062]
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0063]
图1是本发明实施例提供的基于分段精细寻优策略的叶轮多工况气动优化方法流程示意图；
[0064]
图2是本发明实施例提供的基于分段精细寻优策略的离心叶轮多工况气动优化方法流程示意图；
[0065]
图3是本发明实施例提供的型面设计顶点分布局图；
[0066]
图4是本发明实施例提供的主叶片前缘、中部优化变量设置图；
[0067]
图5是本发明实施例提供的主叶片后缘优化变量设置图；
[0068]
图6是本发明实施例提供的分流叶片前缘、中部优化变量分布图；
[0069]
图7是本发明实施例提供的分流叶片后缘优化变量分布图；
[0070]
图8优化前后主叶片前缘、中部的ffd框架控制顶点变形图；
[0071]
图9优化前后主叶片尾缘ffd框架控制顶点变形图；
[0072]
图10优化前后分流叶片前缘、中部的ffd框架控制顶点变形图；
[0073]
图11优化前后主叶片尾缘ffd框架控制顶点变形图；
[0074]
图12是本发明实施例提供的额定工况和常用工况的流量-绝热效率性能曲线图；
[0075]
图13是本发明实施例提供的额定工况和常用工况的总压比-绝热效率性能曲线图；
[0076]
图14是本发明实施例提供的基于分段精细寻优策略的叶轮多工况气动优化装置框图；
[0077]
图15是本发明实施例提供的一种电子设备的结构示意图。
具体实施方式
[0078]
为使本发明要解决的技术问题、技术方案和优点更加清楚，下面将结合附图及具体实施例进行详细描述。
[0079]
如图1所示，本发明实施例提供了一种基于分段精细寻优策略的叶轮多工况气动优化方法，该方法可以由电子设备实现。如图1所示的基于分段精细寻优策略的叶轮多工况气动优化方法流程图，该方法的处理流程可以包括如下的步骤：
[0080]
s1、获取待优化的离心压气机叶轮的几何构型。
[0081]
s2、将待优化的离心压气机叶轮的几何构型输入到构建好的叶轮多工况气动优化模型。
[0082]
其中，叶轮多工况气动优化模型包括全局多工况气动优化模型以及局部多工况气动优化模型。
[0083]
s3、基于待优化的离心压气机叶轮的几何构型、全局多工况气动优化模型以及局部多工况气动优化模型，得到多工况场景下离心压气机叶轮优化后的最佳气动性能及叶片几何外形。
[0084]
可选地，s3中的基于待优化的离心压气机叶轮的几何构型、全局多工况气动优化模型以及局部多工况气动优化模型，得到多工况场景下离心压气机叶轮优化后的最佳气动性能及叶片几何外形包括：
[0085]
s31、对待优化的离心压气机叶轮的几何构型进行全局参数化建模，得到叶轮中叶片吸力面与单位样条曲面的映射模型以及叶片压力面与单位样条曲面的映射模型。
[0086]
一种可行的实施方式中，如图2所示，构建第一阶段全局优化系统，可以采用具有全局特性的型面映射参数化方法，对原始离心压气机叶轮的叶片几何构型进行参数化建模，创建原始叶片吸力面与单位样条曲面的映射方程、压力面与单位样条曲面的映射方程。
[0087]
可选地，s31中的对待优化的离心压气机叶轮的几何构型进行全局参数化建模，得到叶轮中叶片吸力面与单位样条曲面的映射模型以及叶片压力面与单位样条曲面的映射模型包括：
[0088]
s311、获取待优化的离心压气机叶轮的几何构型中叶片各截面型线。
[0089]
一种可行的实施方式中，获取叶片各截面型线，对各截面型线进行归一化，归一化
数学表达式如下式(1)(2)所示：
[0090][0091][0092]
式中，ξ
i,j
和η
i,j
分别是型线归一化后的横坐标和纵坐标，i∈(1,n
p
)，n
p
指的是各个径向截面的点数。j∈(1,ns),ns指的是总的径向截面数。lc指的是径向第j个截面第c段型线长度，lj是第j个截面上各段型线的总长；lr是轴向第i个截面上第r段型线长度，而li是i截面上各段型线的总长。
[0093]
s312、构建单位映射样条曲面。
[0094]
可选地，s312中的单位映射样条曲面为伯恩斯坦bernstein基函数贝塞尔bezier曲面。
[0095]
可选地，s312中的单位映射样条曲面的数学表达式如下式(3)所示：
[0096][0097]
式中，是贝塞尔bezier曲面上的点坐标，p
k,l
是bezier曲面控制顶点，构建m
×
n阶连续的bezier曲面，需要(m+1)
×
(n+1)个控制顶点，即横轴布局m个点，纵轴布局n个点；v和u是两个变化范围为[0,1]的自变量；和是伯恩斯坦bernstein基函数，数学表达式如下式(4)(5)所示：
[0098][0099][0100][0101][0102]
建立待优化的叶片各截面型线与bernstein基bezier曲面点一一对应，通过扰动bezier曲面实现叶片几何构型。
[0103]
s313、根据叶片各截面型线以及单位映射样条曲面，得到叶片吸力面与单位样条曲面的映射模型以及叶片压力面与单位样条曲面的映射模型。
[0104]
一种可行的实施方式中，利用两个6
×
3阶bezier曲面对主叶片和分流叶片分别进行曲面参数化，设计变量布局如图3所示，每个曲面在ξ方向上设7个控制顶点(0、0.1、0.3、0.5、0.7、0.9、1.0)，在η方向上设4个点(0、0.4、0.7、1.0)。为了保证前缘的一阶连续性，固定ξ1、ξ2和ξ'2不动。同时设定吸力面控制顶点的变化与对应压力面控制顶点在构型过程中的变化保持一致，可以确保叶片厚度不会变薄或者畸形，以压力面控制顶点的沿垂直平面
的变化量作为优化变量。
[0105]
s32、设定全局多工况气动优化模型的目标函数和约束条件，基于全局多工况气动优化模型的目标函数和约束条件、叶轮中叶片吸力面与单位样条曲面的映射模型以及叶片压力面与单位样条曲面的映射模型，得到全局优化后的叶片几何构型和气动综合性能。
[0106]
可选地，s32中的基于全局多工况气动优化模型的目标函数和约束条件、叶轮中叶片吸力面与单位样条曲面的映射模型以及叶片压力面与单位样条曲面的映射模型，得到全局优化后的叶片几何构型和气动综合性能包括：
[0107]
基于全局多工况气动优化模型的目标函数和约束条件、叶轮中叶片吸力面与单位样条曲面的映射模型以及叶片压力面与单位样条曲面的映射模型，采用多目标优化算法，得到全局优化后的叶片几何构型和气动综合性能。
[0108]
可选地，多目标优化算法可以是元启发式优化算法或超启发式优化算法等。
[0109]
一种可行的实施方式中，第一阶段多工况全局优化目标函数的数学表达式如下式(6)(7)所示：
[0110]
maxη
roc
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0111]
maxη
noc
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(7)
[0112]
式中，η
noc
是原始叶轮常用工况的绝热效率，η
roc
是原始叶轮额定工况的绝热效率。
[0113]
约束条件的数学表达式如下式(8)(9)所示：
[0114]
π
noc_opt-π
noc_ori
≥0
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(8)
[0115]
π
roc_opt-π
roc_ori
≥0
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(9)
[0116]
式中，π
noc_ori
和π
noc_opt
分别是优化前和优化后常用工况的总压比。π
roc_ori
和π
roc_opt
分别是优化前和优化后额定工况的总压比。
[0117]
一种可行的实施方式中，步骤s32的具体过程可以是：
[0118]
s321、布局得到单位样条曲面的设计顶点参数和单位样条曲面的设计空间参数，并采用拉丁超立方抽样方法初始化样本种群p
t
。
[0119]
s322、基于叶片吸力面与单位样条曲面的映射模型、叶片压力面与单位样条曲面的映射模型以及样本种群p
t
，得到叶片吸力面的变化量以及叶片压力面的变化量。
[0120]
s323、基于叶片吸力面的变化量以及叶片压力面的变化量，得到新的叶片几何构型。
[0121]
s324、基于待优化叶片几何构型生成网格模板文件，根据网格模板文件对新的叶片几何构型进行网格划分，得到新的叶片网格模型。
[0122]
s325、对新的叶片网格模型进行多工况的定常数值模拟计算，得到样本种群p
t
中每一个个体的多工况气动性能参数。
[0123]
s326、设定多工况下新的叶片几何构型优化流程的目标函数和约束条件，基于目标函数和约束条件对样本种群p
t
进行遗传算子操作，得到子代种群p
t*
，进而得到子代种群p
t*
中每一个个体的多工况气动性能参数。
[0124]
s327、基于样本种群p
t
以及子代种群p
t*
组成新种群，并对新种群进行分组。
[0125]
s328、基于分组后的新种群以及空间密度算子模型，生成下一代的父代种群p
t+1
，转去执行步骤s326，直到达到预设迭代次数，输出原始离心压气机叶轮优化后的多工况最佳叶片几何构型。
[0126]
其中，多目标优化算法具体流程：随机生成初始种群pg，开始进化；依据优化耗时总时间指定进化过程的迭代次数或者收敛精度；对pg进行二元锦标赛、二进制交叉和多项式变异，产生新的后代适应度评估，得到每个个体的多目标值；pg和中的所有个体根据非支配排序为f1,f2,f3…fn
；f1,f2,f3…fn
被划分为m1,m2,m3三组。根据空间密度算子模型排序，生成下一个种群；返回第三步，直到满足最大迭代次数或者收敛精度退出寻优。
[0127]
空间密度算子模型排序的具体流程包括：m1和m2被分为一组q；找到空间密度最小的两个个体，其中至少一个个体属于m2；如果一个个体属m1，另一个属于m2，直接从q中删除属于m2的个体；如果两个个体都属于m2，则删除与q中其它个体空间密度最小的个体；返回第二步，直到m1和m2总数达到种群规模。
[0128]
两个个体空间密度的数学表达式如下式(10)所示：
[0129][0130]
和是具有n维决策变量的两个个体。
[0131]
数值计算利用numeca的euranus求解器计算三维稳态reynolds平均navier-stokes方程得到离心叶轮稳态流场，湍流模型采用s-a一方程模型，时间进行采用四阶显式runge-kutta(龙格-库塔)模型，采用具有二阶和四阶人工粘性项的有限体积中心差分格式控制空间离散化过程中的伪数值振荡，利用多重网格、局部时间步长和隐残差加快算法的收敛速度。
[0132]
s33、基于全局优化后的叶片几何构型和气动综合性能，得到气动性能改善潜力较大的多个叶片的局部几何区域。
[0133]
一种可行的实施方式中，基于“人在回路”思想，通过分析步骤s32的寻优结果，得到可进一步开发且潜力大的局部几何区域。在监控优化时可根据离心叶轮气动设计经验对优化流程做出引导和改进，减少盲目性，提高探索成功率和优化效率。
[0134]
通过分析第一阶段的全局优化结果，全局优化方法对叶顶处优化效果明显，而叶根处变化较小，叶顶处发生变形的区域面积较大，而叶根处发生变形的区域面积较小，所以叶根处的构型对性能提升的潜力较大；从全局优化方法可知，为了保证两个曲面在前缘处的连续性和避免网格生成时叶根和叶尖不相交的情况，全局优化方案限制了前缘控制点变化，前缘处有较大的改善空间；根据离心叶轮几何造型对气动性能的影响经验可知，沿型线方向，叶片前缘、中部、尾缘对气动性能影响较大。上述区域仍有较大改善的潜能，为下一步局部优化提供了理论依据。
[0135]
s34、对全局优化后的叶片几何构型进行局部参数化建模，建立叶片局部几何构型与空间控制体的映射模型。
[0136]
一种可行的实施方式中，构建第二阶段局部优化系统，具体地，可以是创建离心叶轮气动优化设计局部参数化方法，建立叶片局部几何构型与空间控制体的映射模型，基于步骤s33得到的几何区域，布局设计变量和设计空间。其中，优化控制变量选定了前缘和中部的叶根区域以及尾缘，共计18
×
2＝36个设计变量，设计空间变量范围为优化顶点与周向相邻控制顶点的连线矢量的一半，其设计变量的局部坐标和变化范围如表1，主叶片设计变量布局如图4和图5，分流叶片设计变量布局如图6和图7。
[0137]
表1
[0138]
前缘中部尾缘(1，2，1)；[-0.5,0.5](6，1，1)；[-0.5,0.5](9，1，1)；[-0.5,0.5](1，2，2)；[-0.5,0.5](6，2，1)；[-0.5,0.5](9，1，2)；[-0.5,0.5] (6，3，1)；[-0.5,0.5](9，1，3)；[-0.5,0.5] (6，1，2)；[-0.5,0.5](9，1，4)；[-0.5,0.5] (6，2，2)；[-0.5,0.5](9，1，5)；[-0.5,0.5] (6，3，2)；[-0.5,0.5](9，3，1)；[-0.5,0.5]
ꢀꢀ
(9，3，2)；[-0.5,0.5]
ꢀꢀ
(9，3，3)；[-0.5,0.5]
ꢀꢀ
(9，3，4)；[-0.5,0.5]
ꢀꢀ
(9，3，5)；[-0.5,0.5]
[0139]
可选地，s34中的对全局优化后的叶片几何构型进行局部参数化建模，建立叶片局部几何构型与空间控制体的映射模型包括：
[0140]
基于b样条基函数的自由曲面变形ffd(free form deformation，自由曲面变形)方法，对全局优化后的叶片几何构型进行局部参数化建模，建立叶片局部几何构型与空间控制体的映射模型。
[0141]
一种可行的实施方式中，首先，利用b样条基函数的ffd方法建立叶片几何构型与空间控制网格的映射模型，将几何构型置于网格控制网格内，通过顶点的位移实现控制体的变形，内置的构型几何跟随控制体发生弹性变形，实现参数化构型。
[0142]
其中，b样条基函数的ffd模型数学表达式如下式(11)所示：
[0143][0144]
式中，为叶片表面坐标，为ffd控制框架上的控制顶点；(s,t,u)为在控制框架内的局部坐标；i,j,k为ffd控制框架三个方向的标号；l,m,n为ffd控制框架在三个方向上的划分数；n
i,d
(s)，n
j,e
(t)，n
k,f
(u)分别对应于d,e,f阶的b样条基函数。n
i,d
(s)，n
j,e
(t)，n
k,f
(u)数学定义按照de boor-cox递推，如下式(12)、(13)所示：
[0145][0146][0147]nj,e
(t)和n
k,f
(u)的数学定义与n
i,d
(s)原理相同。
[0148]
其次，原始控制顶点已知，经过得到新的控制顶点及变形后的控制网格，进而引起叶片表面变形。
[0149]
最后，结合和局部坐标(s,t,u)即可得到表面变形量，局部坐标(s,t,u)可由蒙特卡洛算法求出，则变形后物体表面坐标的数学表达如下式(14)所示：
[0150][0151]
其中，蒙特卡洛算法流程如下：
[0152]
首先，建立映射函数与真实叶片数据点的误差模型，数学表达式如下式(15)所示：
[0153][0154]
式中，s,t,u为映射参数，q为映射值与真实值的误差为，a
real
为真实坐标,是样条曲面控制顶点坐标，i,j,k是ffd控制框架三个方向的标号；l,m,n是ffd控制框架在三个方向上的划分数；n
i,d
(s)，n
j,e
(t)，n
k,f
(u)分别对应于d,e,f阶的b样条基函数，其中s、t和u是映射参数。
[0155]
其次，初始化局部坐标(s,t,u)＝(s0,t0,u0)，计算q0，选定一正数t。
[0156]
再次，在区间[-t,t]上生成随机数向量n，计算q1＝q0(s0+ns,t0+n
t
,u0+nu)。当q1＜q0，(s,t,u)＝(s0+ns,t0+n
t
,u0+nu)，q0＝q1。若随机生成的多组随机向量仍不满足q1＜q0，则令t＝t/2。
[0157]
最后，返回步骤2，直到q0＜ε，(s,t,u)＝(s
best
,t
best
,u
best
)，求出局部坐标。
[0158]
s35、基于全局优化后的叶片几何构型和气动综合性能，设定局部多工况气动优化模型的目标函数和约束条件，基于局部多工况气动优化模型的目标函数和约束条件以及叶片局部几何构型与空间控制体的映射模型，得到多工况场景下离心压气机叶轮优化后的最佳气动性能及叶片几何外形。
[0159]
一种可行的实施方式中，基于第一阶段全局优化结果，设定第二阶段局部优化的目标函数和约束条件，结合步骤s34创建的局部参数化方法、多目标优化算法和数值计算寻优第二阶段最佳叶片几何构型和气动综合性能，完成多工况场景下离心叶轮复杂曲面叶片气动构型的分段精细高效优化。
[0160]
步骤s35的过程可以是：首先，基于叶轮几何构型生成的网格模板文件，对新的叶轮几何构型进行网格划分，得到新的叶轮网格模型。网格主拓扑采用h&i，叶尖间隙拓扑采用ho，采用fine/turbo的autogrid5模块生成.trb文件的网格模板。
[0161]
其次，对新的叶轮网格模型进行数值计算，得到叶轮的多工况气动性能参数。
[0162]
其中，数值计算利用numeca的euranus求解器计算三维稳态reynolds平均navier-stokes方程得到离心叶轮稳态流场，湍流模型采用一方程模型，时间进行采用四阶显式runge-kutta(龙格-库塔)模型，采用具有二阶和四阶人工粘性项的有限体积中心差分格式控制空间离散化过程中的伪数值振荡，利用多重网格、局部时间步长和隐残差加快算法的收敛速度。
[0163]
最后，设定离心压气机叶轮多工况气动优化流程的目标函数和约束条件，并根据多工况气动性能参数，得到ffd控制框架上控制顶点的最优解，主叶片控制顶点变形如图8和图9，分流叶片控制顶点变形如图10和图11。
[0164]
其中，第二阶段多工况局部优化目标函数的数学表达式如下式(16)(17)所示：
[0165]
maxη
roc_baseline
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(16)
[0166]
maxη
noc_baseline
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(17)
[0167]
式中，η
noc_baseline
是第一阶段优化后得到基准叶轮常用工况的绝热效率，η
roc_baseline
是第一阶段优化后得到基准叶轮额定工况的绝热效率。
[0168]
约束条件的数学表达式如下式(18)(19)所示：
[0169]
π
noc_opt-π
noc_ori
≥0
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(18)
[0170]
π
roc_opt-π
roc_ori
≥0
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(19)
[0171]
式中，π
noc_ori
和π
noc_opt
分别是原始叶轮和第二阶段优化后常用工况的总压比。π
roc_ori
和π
roc_opt
分别是原始叶轮和第二阶段优化后额定工况的总压比。
[0172]
多目标优化算法具体流程：随机生成初始种群pg，开始进化；依据优化耗时总时间指定进化过程的迭代次数或者收敛精度；对pg进行二元锦标赛、二进制交叉和多项式变异，产生新的后代适应度评估，得到每个个体的多目标值；pg和中的所有个体根据非支配排序为f1,f2,f3…fn
；f1,f2,f3…fn
被划分为m1,m2,m3三组。根据空间密度算子模型排序，生成下一个种群；返回第三步，直到满足最大迭代次数或者收敛精度退出寻优。
[0173]
空间密度算子模型排序的具体流程包括：m1和m2被分为一组q；找到空间密度最小的两个个体，其中至少一个个体属于m2；如果一个个体属m1，另一个属于m2，直接从q中删除属于m2的个体；如果两个个体都属于m2，则删除与q中其它个体空间密度最小的个体；返回第二步，直到m1和m2总数达到种群规模。
[0174]
两个个体空间密度的数学表达式如下式(20)所示：
[0175][0176]
和是具有n维决策变量的两个个体。
[0177]
基于分段精细寻优策略的离心叶轮多工况气动优化方法，最大限度每个阶段的开发潜力，具有全局特性的bernstein基bezier曲面参数化方法适用第一阶段，实现全局快速构型；具有局部强支撑特性的b样条基函数ffd参数化方法适用第二阶段，实现局部自由灵活构型。同时，在几何控制参数上增加的离心约束有利于光滑叶片的生成。通过本发明的优化方法，获得了离心压气机叶轮的最佳叶片几何构型，提高了多工况的综合气动性能，性能参数改善情况如表2应用实例优化前后气动性能对比(额定工况)和表3应用实例优化前后气动性能对比(常用工况)：
[0178]
表2
[0179][0180]
表3
[0181][0182]
研究结果表明优化后气动性能曲线明显整体上移，额定工况和常用工况的流量-绝热效率性能曲线如图12，额定工况和常用工况的总压比-绝热效率性能曲线如图13，经过两阶段的气动优化，离心压气机叶轮的综合性能得到大幅改善：额定工况的绝热效率提高1.77％，总压比提高0.18％，流量增加了10.4％，裕度提高了7.8％；常用工况的绝热效率提高1.6％，总压比不变，流量增加了10.2％，裕度提高了11.8％。
[0183]
通过应用案例可知，本发明利用合理的优化策略，充分发挥了第一阶段和第二阶段的各自优势、扬长避短，实现了高效、快速、有效地求解了复杂多工况下离心式压气机叶片最佳几何构型和多工况气动性能，实现了形性优化的目的，同时验证了该方法的可行性和普适性。
[0184]
本发明建立了离心叶轮多工况气动优化设计全局参数化模型，构建全局优化系统，利用多目标优化算法求解多工况绝热效率最大值，进而得到第一阶段全局寻优后的最佳叶片几何构型和气动综合性能；基于第一阶段的优化结果和“人在回路”思想，建立叶片气动优化设计局部参数化模型，构建局部优化系统，利用多目标优化算法进一步求解多工况绝热效率最大值，进而得到第二阶段局部寻优后的最佳叶片几何构型和气动综合性能，通过两阶段寻优实现了多工况场景下的离心叶轮复杂曲面叶片精细化高效构型。
[0185]
本发明实施例中，提出了基于分段精细寻优策略的离心叶轮多工况气动优化方法，建立了两个优化阶段：基于全局优化系统的第一阶段和基于局部优化系统的第二阶段，研究了空间复杂曲面参数化方法、多目标优化算法及寻优策略的最佳匹配关系，求解了等熵效率最大值，通过两次寻优实现了复杂曲面叶片精细化高效构型，最大限度每个阶段的开发潜力，实现了寻优过程中减少设计变量、缩小设计空间、加快收敛速度、增强探索能力及形性优化控制的多重目的，具有一定的推广应用价值。
[0186]
如图14所示，本发明实施例提供了一种基于分段精细寻优策略的叶轮多工况气动优化装置1400，该装置1400应用于实现基于分段精细寻优策略的叶轮多工况气动优化方法，该装置1400包括：
[0187]
获取模块1410，用于获取待优化的离心压气机叶轮的几何构型。
[0188]
输入模块1420，用于将待优化的离心压气机叶轮的几何构型输入到构建好的叶轮多工况气动优化模型；其中，叶轮多工况气动优化模型包括全局多工况气动优化模型以及局部多工况气动优化模型。
[0189]
输出模块1430，用于基于待优化的离心压气机叶轮的几何构型、全局多工况气动优化模型以及局部多工况气动优化模型，得到多工况场景下离心压气机叶轮优化后的最佳气动性能及叶片几何外形。
[0190]
可选地，输出模块1430，进一步用于：
[0191]
s31、对待优化的离心压气机叶轮的几何构型进行全局参数化建模，得到叶轮中叶片吸力面与单位样条曲面的映射模型以及叶片压力面与单位样条曲面的映射模型。
[0192]
s32、设定全局多工况气动优化模型的目标函数和约束条件，基于全局多工况气动优化模型的目标函数和约束条件、叶轮中叶片吸力面与单位样条曲面的映射模型以及叶片压力面与单位样条曲面的映射模型，得到全局优化后的叶片几何构型和气动综合性能。
[0193]
s33、基于全局优化后的叶片几何构型和气动综合性能，得到多个叶片的局部几何区域。
[0194]
s34、基于多个叶片的局部几何区域，对全局优化后的叶片几何构型进行局部参数化建模，建立叶片局部几何构型与空间控制体的映射模型。
[0195]
s35、基于全局优化后的叶片几何构型和气动综合性能，设定局部多工况气动优化模型的目标函数和约束条件，基于局部多工况气动优化模型的目标函数和约束条件以及叶片局部几何构型与空间控制体的映射模型，得到多工况场景下离心压气机叶轮优化后的最佳气动性能及叶片几何外形。
[0196]
可选地，输出模块1430，进一步用于：
[0197]
s311、获取待优化的离心压气机叶轮的几何构型中叶片各截面型线。
[0198]
s312、构建单位映射样条曲面。
[0199]
s313、根据叶片各截面型线以及单位映射样条曲面，得到叶片吸力面与单位样条曲面的映射模型以及叶片压力面与单位样条曲面的映射模型。
[0200]
可选地，单位映射样条曲面为伯恩斯坦bernstein基函数贝塞尔bezier曲面。
[0201]
可选地，单位映射样条曲面的数学表达式如下式(1)所示：
[0202][0203]
式中，是贝塞尔bezier曲面上的点坐标，p
k,l
是bezier曲面控制顶点，控制点总数为(m+1)
×
(n+1)，m为横轴布局点个数，n为纵轴布局点个数；v和u是两个变化范围为[0,1]的自变量；和是伯恩斯坦bernstein基函数。
[0204]
可选地，输出模块1430，进一步用于：
[0205]
基于全局多工况气动优化模型的目标函数和约束条件、叶轮中叶片吸力面与单位样条曲面的映射模型以及叶片压力面与单位样条曲面的映射模型，采用多目标优化算法，得到全局优化后的叶片几何构型和气动综合性能。
[0206]
可选地，多目标优化算法为元启发式优化算法或超启发式优化算法。
[0207]
可选地，输出模块1430，进一步用于：
[0208]
基于b样条基函数的自由曲面变形ffd方法，对全局优化后的叶片几何构型进行局部参数化建模，建立叶片局部几何构型与空间控制体的映射模型。
[0209]
可选地，输出模块，进一步用于：
[0210][0211]
式中，为叶片表面坐标，为ffd控制框架上的控制顶点；(s,t,u)为
在控制框架内的局部坐标；i,j,k为ffd控制框架三个方向的标号；l,m,n为ffd控制框架在三个方向上的划分数；n
i,d
(s)，n
j,e
(t)，n
k,f
(u)分别对应于d,e,f阶的b样条基函数。
[0212]
本发明实施例中，提出了基于分段精细寻优策略的离心叶轮多工况气动优化方法，建立了两个优化阶段：基于全局优化系统的第一阶段和基于局部优化系统的第二阶段，研究了空间复杂曲面参数化方法、多目标优化算法及寻优策略的最佳匹配关系，求解了等熵效率最大值，通过两次寻优实现了复杂曲面叶片精细化高效构型，最大限度每个阶段的开发潜力，实现了寻优过程中减少设计变量、缩小设计空间、加快收敛速度、增强探索能力及形性优化控制的多重目的，具有一定的推广应用价值。
[0213]
图15是本发明实施例提供的一种电子设备1500的结构示意图，该电子设备1500可因配置或性能不同而产生比较大的差异，可以包括一个或一个以上处理器(central processing units，cpu)1501和一个或一个以上的存储器1502，其中，存储器1502中存储有至少一条指令，至少一条指令由处理器1501加载并执行以实现下述基于分段精细寻优策略的叶轮多工况气动优化方法：
[0214]
s1、获取待优化的离心压气机叶轮的几何构型。
[0215]
s2、将待优化的离心压气机叶轮的几何构型输入到构建好的叶轮多工况气动优化模型；其中，叶轮多工况气动优化模型包括全局多工况气动优化模型以及局部多工况气动优化模型。
[0216]
s3、基于待优化的离心压气机叶轮的几何构型、全局多工况气动优化模型以及局部多工况气动优化模型，得到多工况场景下离心压气机叶轮优化后的最佳气动性能及叶片几何外形。
[0217]
在示例性实施例中，还提供了一种计算机可读存储介质，例如包括指令的存储器，上述指令可由终端中的处理器执行以完成上述基于分段精细寻优策略的叶轮多工况气动优化方法。例如，计算机可读存储介质可以是rom、随机存取存储器(ram)、cd-rom、磁带、软盘和光数据存储设备等。
[0218]
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成，也可以通过程序来指令相关的硬件完成，所述的程序可以存储于一种计算机可读存储介质中，上述提到的存储介质可以是只读存储器，磁盘或光盘等。
[0219]
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。