一种多能互补条件下水轮机转轮的多工况优化设计方法

1.本发明属于水力机械旋转部件优化领域，具体为一种多能互补条件下水轮机转轮的多工况优化设计方法。
技术背景
2.目前新能源发电发展迅速，新能源正朝着绿色低碳型方向转变。风电、光伏作为可再生清洁能源，来源丰富且容易获取，风光水多能互补发电作为能源转型利用的一种新模式，也是未来能源发展的重要方向之一。风光水多能互补系统中，风力、光伏发电本身强烈的随机性、灵活性要求水电机组作为补偿出力，致使水电机组频繁运行在偏工况区。而转轮作为混流式水轮机的核心部件之一，对混流式水轮机的性能、安全稳定运行起着至关重要的作用。
3.现有的转轮大多数在高负荷区可以稳定运行，当混流式水轮机处于多能互补条件下，现有转轮难以满足风光水多能互补发电的运行要求。


技术实现要素：

4.针对现有技术中存在的问题，本发明提供一种多能互补条件下水轮机转轮的多工况优化设计方法，通过对混流式水轮机在不同工况下的进口安放角、出口安放角进行优化设计，能够快速并准确的提高风光水多能互补条件下混流式水轮机转轮的设计优化效率，使水轮机在风光水多能互补条件下能够稳定高效运行。
5.本发明是通过以下技术方案来实现：
6.一种多能互补条件下水轮机转轮的多工况优化设计方法，包括如下步骤：
7.步骤1，建立混流式水轮机转轮初始叶片的三维模型，调整b
é
zier曲线的阶数，使b
é
zier曲线表达出所述三维模型的安放角几何分布特征，得到参数化后的优化设计变量；
8.步骤2，将参数化后的优化设计变量进行建模，得到初步优化设计后的混流式水轮机转轮叶片的三维模型，将所述三维模型的进口安放角和出口安放角的变化范围进行约束，得到约束后优化设计变量；
9.步骤3，根据参数化后的优化设计变量和约束后优化设计变量，构建优化设计的样本空间，之后对其中的每个混流式水轮机转轮模型进行水力性能数值计算，得到每个混流式水轮机分别在超低负荷工况、低负荷工况和额定负荷区运行下的运行效率，建立每个混流式水轮机转轮在超低负荷工况、低负荷工况和额定负荷区运行下的优化目标函数，最后构建约束后优化设计变量与优化目标函数间的响应面模型；
10.步骤4，当响应面模型满足要求时进行步骤5，当响应面模型不满足要求，增加步骤3中样本空间的样本数后重新构建所述的响应面模型，直至得到的响应面模型满足要求，然后进行步骤5；
11.步骤5，对步骤4得到的响应面模型进行全局寻优，得到最优解集，从最优解集中选择混流式水轮机及其转轮叶片加权平均效率均提高的解集，所述的解集对应多能互补条件
下多工况优化后的水轮机转轮。
12.优选的，步骤1将混流式水轮机转轮初始叶片的三维模型沿叶高方向均分为4个三维截面，使b
é
zier曲线表达出4个三维截面的安放角几何分布特征。
13.进一步，步骤1所使用的b
é
zier曲线为四点三次b
é
zier曲线。
14.进一步，步骤2通过约束四点三次b
é
zier曲线的四个控制点对所述的进口安放角和出口安放角的变化范围进行约束，具体将四点三次b
é
zier曲线的四个控制点的纵坐标变化范围分别约束在(-4.5
°
～+4.5
°
)、(-3
°
～+3
°
)、(-3.5
°
～+3.5
°
)、(-5.5
°
～+5.5
°
)。
15.优选的，步骤3中，基于拉丁超立方抽样实验设计方法，生成若干组样本点，形成样本空间，之后形成若干个混流式水轮机转轮模型。
16.优选的，步骤3采用如下过程对样本空间内每个混流式水轮机转轮模型的水力性能数值进行计算：
17.对样本空间中的每个混流式水轮机转轮叶片模型进行批量网格划分，同时对每个混流式水轮机剩余部件进行几何建模并生成网格，将得到的所有混流式水轮机超低负荷工况、低负荷工况和额定负荷区工况下的全流道数值模拟离散模型，对包含水轮机所有过流部件的离散模型进行数值计算，得到每个混流式水轮机在超低负荷工况、低负荷工况和额定负荷区工况下的运行效率。
18.优选的，步骤3基于几何平均超传递近似法，确定优化目标函数中超低负荷工况、低负荷工况和额定负荷区下的效率权重系数，优化目标函数为η
max
＝ω1η1+ω2η2+ω3η3，其中η
max
为加权平均效率的最大值，ω1，ω2,ω3分别为超低负荷工况、低负荷工况和额定负荷区下的效率权重系数。
19.优选的，步骤3利用多元多项式构建约束后优化设计变量与优化目标函数间的响应面模型，得到kriging模型。
20.优选的，步骤4按照拉丁超立方抽样方法生成测试样本，将响应面模型预测所得的优化目标函数值与测试样本的数值计算结果进行对比，若两者的偏差小于或等于5％，则响应面模型精度满足要求；若两者的偏差大于5％，则响应面模型精度不满足要求。
21.优选的，步骤5采用多目标遗传算法对步骤4得到的响应面模型进行全局寻优。
22.与现有技术相比，本发明具有以下有益的技术效果：
23.本发明一种多能互补条件下水轮机转轮的多工况优化设计方法，先使b
é
zier曲线能准确表达出三维模型的安放角几何分布特征，之后便可通过改变翼型安放角即可修改转轮叶片截面上的翼型形状，实现对进口安放角和出口安放角的优化，进而修改转轮叶片的形状。为了更好的控制混流式水轮机转轮叶片进、出口安放角的变化，并且避免叶片变形过大，对各个设计变量的变化范围进行约束，之后根据参数化后的优化设计变量和约束后优化设计变量，构建了优化设计的样本空间，便可对每个混流式水轮机转轮模型进行水力性能数值计算，得到每个混流式水轮机分别在超低负荷工况、低负荷工况和额定负荷区运行下的运行效率，以便建立每个混流式水轮机转轮在超低负荷工况、低负荷工况和额定负荷区运行下的优化目标函数，最后得以构建约束后优化设计变量与优化目标函数间的响应面模型，该响应面模型能预测优化目标函数值，可判断响应面模型是否满足要求，最后对响应面模型进行全局寻优，可得到最优解集，从这些最优解集中可选择混流式水轮机及其转轮叶片加权平均效率均提高的解集，这些解集对应多能互补条件下多工况优化后的水轮机转
轮。本发明在得到多工况的优化水轮机转轮模型时，结合风、光、水多能互补条件下水轮机在超低负荷工况区、低负荷工况区的运行时间要求，各个工况设置不同的权重因子，提高超低负荷工况区、低负荷工况区的占比，满足风、光、水多能互补条件下水电机组频繁运行在超低负荷工况区、低负荷工况区的要求，提高了其在多能互补条件下混流式水轮机的运行效率，达到多工况优化设计的目的。本发明在设置优化参数范围时考虑工程实际需求，有效的避免了在多工况优化设计过程中产生变形过大的叶片。
24.进一步的，本发明在优化设计过程中使用响应面模型和多目标遗传算法相结合的方法，对混流式水轮机转轮叶片进行全局寻优，形成pareto最优，避免陷入局部最优，具有快速、准确的优点。
附图说明
25.图1为本发明所述的多工况优化设计方法中的设计方法流程图。
26.图2为本发明所述的多工况优化设计方法中的转轮优化前叶片三维翼型示意图。
27.图3为本发明所述的多工况优化设计方法中的基于bezier曲线的叶片安放角拟合示例图。
28.图4为本发明所述的多工况优化设计方法中的优化前后转轮叶片模型对比的俯视图。
29.图5为本发明所述的多工况优化设计方法中的优化前后转轮叶片模型对比的主视图。
30.图中：1-初始叶片，2-优化叶片。
具体实施方式
31.下面结合附图及实施例对本发明的技术方案作进一步的描述说明，所述是对本发明的解释而不是限定。
32.本发明一种风光水多能互补条件下混流式水轮机转轮的多工况优化设计方法，以超低负荷工况、低负荷工况、额定负荷区下的水轮机加权平均效率为优化目标，叶片进口安放角、出口安放角为待优化的几何参数，并结合工程实际，不产生几何变形过大的叶片，如图1所示，具体包括如下步骤：
33.步骤一：转轮叶片的参数化表达
34.转轮叶片的参数化表达是开展优化设计的基础，其作用是利用较少的设计参数来准确的描述转轮叶片的几何形状，以便于能够便捷的通过改变设计参数得到新的转轮叶片。采用b
é
zier曲线实现转轮叶片的参数化表达，通过改变翼型安放角即可修改转轮叶片截面上的翼型形状，实现对进口安放角和出口安放角的优化，进而修改转轮叶片的形状。
35.先建立混流式水轮机转轮初始叶片的三维模型，将该叶片的三维模型沿叶高方向均分为4个三维截面，之后调整b
é
zier曲线的阶数，使b
é
zier曲线能准确表达出该三维模型的4个三维截面的安放角几何分布特征，得到参数化后的优化设计变量，具体通过以下方式实现：
36.先编写混流式水轮机转轮初始叶片的三维造型文件，利用该文件生成混流式水轮机转轮初始叶片的三维模型。将该叶片的三维模型沿叶高方向均分为如图2所示的4个三维
截面，基于b
é
zier曲线对混流式水轮机转轮初始叶片的不同截面安放角进行参数化控制，通过调整b
é
zier曲线的阶数，能够保证准确的表达出混流式水轮机转轮叶片4个三维截面安放角的几何分布特征。描述任意一个截面安放角的b
é
zier曲线控制点个数决定了混流式水轮机转轮叶片的优化设计变量，即p＝(p1,l,pn)，其中n为优化设计变量的总数。
37.具体在实施时，采用四点三次b
é
zier曲线对各个截面安放角几何分布进行参数化，每个三维截面的控制点各4个优化设计变量，因此混流式水轮机转轮叶片上的优化设计变量总计为16个优化设计变量，其中如图3所示的任意一条四点三次b
é
zier曲线两端的控制点直接决定了进口安放角和出口安放角的大小。
38.步骤二，逆向建模，设计变量的约束
39.通过逆向建模对转轮叶片参数化的进度进行校验，为了更好的控制转轮叶片进口安放角、出口安放角的变化，并且避免叶片变形过大，对各个优化设计变量的变化范围进行约束。
40.具体地，将步骤1得到的优化设计变量重新进行建模，比较步骤1的初始叶片和该叶片的拟合程度，经校验，参数化没有让叶片形状失真。
41.进口安放角和出口安放角的变化范围通过约束b
é
zier曲线的四个控制点p0、p1、p2、p3的纵坐标来实现。变化范围为p0±
4.5
°
、p1±3°
、p2±
3.5
°
、p3±
5.5
°
。在混流式水轮机转轮中，叶片进口安放角、出口安放角会对转轮的性能产生重要影响。因此将进口安放角和出口安放角的优化设计变量约束范围适当增加。
42.步骤三：不同工况权重因子确定、优化目标的建立
43.以不同运行工况的水轮机转轮加权平均效率为优化目标进行样本空间内每个混流式水轮机转轮的水力性能评价。而加权平均效率的计算最重要的步骤是权重因子的确定，因此需要确定合理的不同工况的权重因子。
44.以混流式水轮机转轮不同工况的加权平均效率为优化目标，即η＝(η1,l,ηm)，其中，m为优化运行工况数。转轮优化运行工况为超低负荷工况、低负荷工况、额定负荷区，即m＝3。确定优化目标函数为：η
max
＝ω1η1+ω2η2+l+ωmηm。其中η
max
为不同工况加权平均效率的最大值，ω1，ω2lωm为各个运行工况下的效率权重系数。
45.基于几何平均超传递近似法，确定混流式水轮机在各个工况的效率权重系数，具体通过以下方式实现：
46.首先在各工况间通过两两比较进行不同工况的相对重要性评估，生成二元互补矩阵；利用二元互补矩阵构造出超传递近似矩阵；最后利用特征向量法求超传递近似矩阵的最大特征值，该最大特征值所对应的特征向量即各个工况的权重系数。
47.在具体实施时，一个具体的水电机组长时间运行在超低负荷工况、低负荷区、额定负荷区，计算得到ω1＝0.48，ω2＝0.35，ω3＝0.17。
48.步骤四：建立样本空间
49.根据步骤1得到的优化设计变量及步骤2得到的优化设计变量的约束范围构建优化设计的样本空间。
50.根据步骤一得到的优化设计变量和步骤二得到的变量约束条件，基于拉丁超立方抽样实验设计方法，生成m组样本点，形成样本空间，得到m个混流式水轮机转轮叶片模型。
51.具体在实施时，考虑计算时间以及精度，混流式水轮机转轮优化设计的样本点个
数设置为200个，具体如表1～表5。
52.表1混流式水轮机转轮生成的1-40个样本空间，其中0为初始。
[0053][0054][0055]
表2混流式水轮机转轮生成的41-80个样本空间。
[0056] x1x2x3x4x5x6x7x8x9x10x11x12x13x14x15x164140.3032.1638.6446.7835.4333.0043.8462.7339.9435.0455.6268.5734.8533.1469.1976.284241.3527.6235.7048.9234.9437.7640.1766.9036.5537.2762.4267.0335.5938.3267.4577.904338.5230.8438.0552.8233.6039.5637.8464.0443.8136.4154.5070.8736.0433.2875.6683.284443.0931.7834.7454.0637.2737.5444.1257.6440.0238.4259.8264.0437.8337.6265.5070.294543.9427.4533.8247.0838.7236.4840.4960.8735.9534.3764.4972.0538.2132.8274.3378.904642.1228.7036.5152.2733.5338.5744.2460.0642.8539.5763.9071.0136.9937.7674.9682.594740.7029.1438.4252.4232.8236.6736.2757.8237.1536.5655.2774.0132.0534.4764.8775.904839.4127.3633.4651.8736.2535.4642.9965.2836.1135.5658.4070.5232.8632.3076.0772.134943.3727.5339.7446.7332.4735.1340.3356.1537.3534.8259.7064.8734.4732.6866.3373.985044.3027.0639.9648.3237.5536.1938.2856.0238.8633.5560.2462.7832.1933.5972.3279.905140.3429.4434.0847.4335.1533.1942.8764.4743.6134.4853.3263.8337.9335.1769.8879.055237.7128.6533.0953.1636.5338.8638.4059.8141.2534.2260.5371.9133.9836.8869.7475.595339.9427.4237.5048.4734.4135.5341.7464.1037.2735.9356.2174.2136.8832.7563.6280.055441.8827.7433.7153.1135.8634.4337.0764.6040.8539.7656.5171.4336.0836.8165.3669.445537.1127.6537.1749.0731.8336.5240.0558.9439.7436.3457.2866.2737.5134.1972.7371.985643.9026.9834.8951.9236.3533.6643.1564.7242.7734.5656.9268.0135.8331.7466.2080.90
5740.5028.6238.1252.6237.0936.6339.6562.9838.4639.1660.4172.8932.7536.3275.5280.365839.7331.2533.4951.6233.9534.8737.6766.3439.4240.1364.0871.6433.5233.7363.0680.745941.8327.0936.1447.7833.6740.2240.4558.6336.1534.9658.5267.9436.4333.5667.0369.756039.3329.0633.1346.9336.4233.5937.3559.1338.7435.8961.9576.0338.2532.4073.5076.136141.1928.1238.6054.3632.8633.8542.7959.8739.1037.4563.3162.9237.6532.1969.1275.216239.1330.4934.1945.0336.6335.2737.1961.3741.4136.5256.6969.5435.7635.2467.9376.216340.9531.2839.5949.6232.4336.0841.9056.7036.2333.5958.0568.9932.4034.6865.2271.216438.7628.7635.9253.7137.8036.8140.7456.6437.3137.9459.2375.2633.0735.9065.9981.056537.1929.7038.4945.4335.9339.6338.2462.8636.8736.0462.3165.9936.2936.1566.7578.446636.9931.7535.4848.2737.6938.7538.1259.0041.2934.0460.7774.3533.3536.6474.0673.296740.8726.5434.3046.9832.5438.4239.2567.5240.9336.6761.3670.0336.9237.3071.6980.516836.8329.9033.1647.9337.6634.7241.2268.0837.7038.9459.6475.6834.3631.9872.8070.066942.0031.8938.9348.5732.8935.3143.2767.2139.5834.8955.0968.5033.2836.2566.2672.217040.1828.4740.0353.8636.1435.3537.5560.9336.3138.1264.0275.7537.6233.3863.7677.447143.2132.1037.2445.6832.1233.5242.9165.9641.5339.4261.8371.5038.1434.2271.9780.677241.4331.4337.0648.1336.6034.9142.7556.8337.9438.9062.6665.4335.5232.5463.1374.757343.0529.0334.1254.2635.5440.0342.2364.6642.1340.2457.7575.1235.0337.3464.9472.527440.2227.5638.0954.7637.4134.2140.3761.5542.0539.0963.9664.9431.7031.6769.8178.367538.4828.5636.1050.9236.3235.6439.0465.0340.4639.7257.6970.1731.7733.2466.1376.827639.7731.0540.3753.9135.2939.5238.5267.4643.1739.2064.9768.6433.3133.9876.0078.057741.2330.2037.6850.2732.6439.8542.3560.1242.4537.1263.0768.2237.7932.8965.7170.607841.9630.1733.6048.9735.2236.4539.4161.3042.3337.7554.9165.3634.6437.0263.4874.907942.8427.4734.6751.7232.7535.9341.3067.7737.0735.4555.5668.4333.2435.6664.7373.448037.7928.0038.6745.2838.4038.3940.8657.9541.1739.5353.9773.5238.1836.4372.3980.97
[0057]
表3混流式水轮机转轮生成的81-120个样本空间。
[0058] x1x2x3x4x5x6x7x8x9x10x11x12x13x14x15x168137.7531.6637.8353.8132.2635.0237.6364.5341.0533.7053.7974.6337.4835.1368.7781.668238.0029.9934.4945.5332.6140.1843.1965.3436.5940.8354.2673.6636.5334.4362.6581.438342.5228.8840.1850.6734.6638.1743.3156.7740.7035.0856.7568.9235.3131.7063.9075.288444.2629.4139.0454.5135.5837.3236.8767.3336.4734.4459.8873.1034.8234.2667.1777.368542.6831.6937.2054.0137.2438.9738.0057.8937.5433.8154.2073.7337.0635.8764.1176.908638.3626.5136.2154.1132.3635.9043.5566.4040.7738.5356.8664.1134.7138.1474.2782.288741.9229.6734.7849.8233.1033.4142.4359.9436.3535.1556.1070.5936.5735.5273.9280.288842.4427.1238.3149.7234.8336.9238.0462.5541.4939.0560.5967.5233.1735.7667.5273.068940.1031.9534.8553.2638.2936.7439.5360.5041.5733.7457.4063.3436.1534.0166.6182.669036.4226.4839.2351.4235.6835.9742.8365.9040.3835.8656.9875.4734.4337.7970.1678.679139.2127.3934.2748.6736.2838.9039.2959.6341.6538.3163.6767.1734.6838.2165.1570.759238.4026.4538.7151.9732.0437.4337.4366.8339.2239.8761.6074.7736.1835.4863.0073.219342.7630.0839.8153.6638.4739.4838.1660.5639.9837.5354.1465.2233.4936.0170.3081.829437.3129.2337.5746.4337.2034.6139.9763.6041.7737.4261.1866.5532.5138.1865.6482.059542.2830.9639.8547.3833.0339.2636.5568.2036.2736.2364.2066.4833.5936.5072.6772.759644.0230.3138.9747.8336.6737.5137.3160.6241.7336.2658.5872.4734.0132.5867.6681.979742.9729.8538.0147.1834.7639.8136.3562.2437.2336.7554.6264.8035.3438.4670.7972.679842.6027.1835.2648.0838.2234.5436.8357.1438.6238.0153.6169.6833.6636.0471.3484.519937.8429.7635.6645.1332.0137.2137.0357.0839.9037.1959.4171.9835.6236.7173.9974.9810043.6930.9335.6248.8233.7440.2537.7660.8136.0337.7958.6469.8933.0334.0570.9374.2910143.8129.2037.9053.2132.7234.9442.3157.2641.9336.9764.7374.4934.0838.3565.7880.2010240.9929.5240.2952.6735.4736.1239.8561.2443.5334.6356.2768.1536.9536.2966.8271.7510339.2527.9138.5350.3736.8834.1436.9164.2942.2140.4661.0664.1832.9633.8475.9478.1310440.3830.1436.9845.3833.7033.3043.4761.4938.9036.4963.4376.5234.9233.8068.7074.5210542.0426.7433.7547.7333.4634.1038.8858.3839.7034.4158.2864.4535.4538.0463.8379.4410642.3231.1336.3247.9834.6239.3044.1667.9539.2639.3858.8870.3838.3532.5174.8975.1310738.8528.7337.3245.9835.3634.8341.1062.6740.5840.2060.8376.4534.9934.6474.4769.6010843.9831.9236.6245.3334.8036.9939.1663.6643.2135.7162.9063.0637.0936.5367.5970.9010942.3627.8539.6351.8233.7735.0941.7867.6438.3040.5461.6566.2032.4737.2071.6269.3711041.6730.7236.8045.7835.1235.7942.5163.5441.3336.1264.8565.2938.2837.1670.5170.37
11137.5531.1933.2746.5337.1633.0438.8466.9639.3438.7962.7863.7636.2234.1263.2074.1311242.6428.8238.8250.4732.9337.6539.1258.0140.6239.2757.3467.8735.7332.0968.8472.0611336.5030.0239.1546.1835.7534.2937.7158.5141.0937.0856.8166.6232.1236.0866.9674.6711437.6329.5837.6553.0638.4435.3841.9456.3336.4337.6062.2568.2937.3737.9067.8075.8211540.0630.5234.2350.8232.0835.8238.3262.6142.8136.1954.5676.5934.3332.7267.7372.4411641.4731.4838.7950.2235.0836.8839.9358.2042.4938.4663.1374.1534.1934.7566.4781.3611741.7128.5037.7947.5336.5639.9237.8058.3237.7436.8656.0466.0631.9832.0271.5575.5211843.6128.2634.5250.3234.8737.0342.0659.3139.3039.6158.3473.9437.2033.0374.2084.1211940.6226.4240.2249.6734.3039.1238.5666.2740.1837.9063.1969.1334.5032.2669.5377.8212042.4826.6536.0354.6635.3334.9838.6464.4142.6935.6359.1172.1236.0131.8864.3283.12
[0059]
表4混流式水轮机转轮生成的121-160个样本空间。
[0060][0061][0062]
表5混流式水轮机转轮生成的161-200个样本空间。
[0063][0064][0065]
步骤五：样本空间内混流式水轮机转轮水力性能数值分析
[0066]
水轮机转轮水力性能的数值分析包含两个基本过程。1、对计算域进行离散，即对计算域进行网格划分。2、水轮机转轮内流场的数值计算，即采用计算流体动力学理论对整个水轮机转轮进行流场计算。
[0067]
对步骤四生成的样本空间中的每个混流式水轮机转轮进行水力性能的数值计算，具体通过以下方式实现：
[0068]
首先编写网格划分脚本文件，对步骤四生成的样本空间中的每个混流式水轮机转轮叶片模型进行批量网格划分，同时对每个混流式水轮机的剩余部件(包括蜗壳、尾水管、活动导叶、固定导叶等)进行几何建模并生成网格，将得到的所有混流式水轮机各个工况下的全流道数值模拟离散模型，对不同工况下包含水轮机所有过流部件的离散模型进行数值计算，得到每个混流式水轮机多工况下的运行效率。
[0069]
步骤六：建立响应面模型
[0070]
构建步骤二约束后的优化设计变量与步骤三的优化目标函数间的响应面模型，获得响应面特征参数，并对步骤三中的优化目标函数进行求解分析并判断该响应面模型精度是否符合要求，若符合要求则进行下一步；若不符合，则重复步骤二～步骤六直至其获得的
响应面模型精度满足要求，具体过程说明如下。
[0071]
构建响应模型的具体措施如下：
[0072]
利用多元多项式构建步骤一得到的优化设计变量与步骤三中的优化目标函数间的响应模型符合kriging(克里金)模型，再根据步骤五中的样本空间的数值计算结果，得到响应面特征参数，即各个工况下的水轮机的水力效率。
[0073]
按照拉丁超立方抽样方法生成10组测试样本，利用该测试样本检验上一步所建立的响应面模型的精度，将响应面模型预测所得的优化目标函数值与测试样本的数值计算结果进行对比，若两者结果偏差大于5％，则响应面模型精度不满足要求，需增加样本数至步骤四的样本空间，重新进行数值计算并建立新的响应面模型直至其精度满足优化要求；若精度满足要求则输出响应面模型。
[0074]
本发明所列举的200个样本空间所形成的响应面模型预测所得的优化目标函数值与测试样本的数值计算结果进行对比，两者结果偏差小于5％。
[0075]
步骤六：遗传全局寻优算法优化
[0076]
采用多目标遗传算法在构建的响应面模型上进行全局寻优，确保最终解为全局最优，并且其解的形式为pareto最优解集，可根据实际需要在pareto最优解集中进行择优。
[0077]
具体地，选取moga遗传算法进行优化，以设置pareto解百分比范围作为判定moga遗传算法是否收敛的准则，包括如下步骤：
[0078]
步骤a.设置种群大小并随机产生初始种群，即初始样本数；
[0079]
步骤b.设置迭代样本数，计算各个样本的目标函数值；
[0080]
步骤c.根据各个样本点的适应度进行选择、交叉及变异，生成新的群体；
[0081]
步骤d.检查是否收敛，若没有，转至步骤b，若收敛，则进入步骤e；
[0082]
步骤e.输出候选解集中的pareto最优解集及相应的优化目标函数值；
[0083]
步骤f.根据实际情况，从pareto解集中挑选出最适合该优化方案的解；
[0084]
moga遗传算法的具体参数设置见表6。
[0085]
表6 moga遗传算法的设置参数
[0086][0087][0088]
步骤七：选取最优结果
[0089]
根据步骤六得到的最优叶片的解集中选取一组相较于初始叶片，超低负荷工况和低负荷工况效率有所提高且额定负荷区效率变化较小的叶片作为最终优化叶片。
[0090]
选取的最终优化叶片为步骤4中编号为71的叶片，优化前后加权平均效率对比如
表7所示。优化前叶片的三维翼型如图2所示，优化前后的转轮叶片模型对比图如图4和图5所示。从表7可以看出优化后转轮叶片相较于原始叶片，转轮加权平均效率提高0.32％，整机运行加权平均效率提高2.95％。
[0091]
表7优化前后加权平均效率对比
[0092]