本发明属于应用于叶轮水力性能的优化技术，具体是混流泵叶轮的水力优化设计方法，以提高混流泵性能在实际运行中的稳定性，降低其对不确定因素的敏感度。
背景技术：
：混流泵因其具有使用流量、扬程变化范围大、高效区宽、在小流量高扬程环境下能对性能曲线中的驼峰现象有较好的改善、结构紧凑、易气动等优点，被广泛地用于电站冷却水循环系统、海水淡化装置、农田灌溉、市政给排水等国民经济领域。其次，混流泵是一种兼具离心泵高扬程和轴流泵大流量两者优点的一种泵型，其可以很好地吸收两者的优点，并逐渐被应用到传统的离心泵和轴流泵领域。随着计算机技术的高速发展，在混流泵的叶轮设计中用到最多的是将计算流体力学(简称cfd)中的分析和智能优化算法相结合的水力优化设计方法。但是，在已有的对混流泵的水力优化设计中，大多数是只考虑其设计工况的单点优化和考虑多个离散工况点的多点优化，且优化过程中往往直接选取叶轮的几何参数作为优化变量。例如中国专利申请号cn201510686219.6的文献中公开了一种超低比转速离心泵叶轮多工况多目标水力优化方法，以叶轮的进口直径、出口直径、出口安放角、出口宽度、包角为优化变量，对三个工况下的加权平均效率和叶轮加权平均径向力进行了优化；又例如中国专利申请号cn201010520561.6的文献中同样也以叶轮的15个几何参数为优化变量，公开了一种基于cfd的离心泵多工况水力优化方法；中国专利申请号cn201310638000.x的文献中给出了一种只考虑设计工况点的双流道多目标优化方法。这些传统优化方法中的环境参数、优化变量及目标函数均是确定性的，不能反映客观存在的不确定因素对设计方案的影响，且通过该些传统优化设计方法得到的最优解的性能往往在某个特定工况条件下是最优的，一旦使用条件发生变化，混流泵的性能就会发生急剧的恶化、运行的稳定性就会变差，甚至会出现不能正常工作的情况。经检索，尚未发现有关将流量工况参数的随机不确定性波动考虑在内的鲁棒优化设计方法。技术实现要素：本发明旨在提供一种混流泵叶轮的鲁棒优化设计方法，考虑到流量工况参数的随机不确定性波动，求解混流泵叶轮鲁棒优化的多目标问题，进而得到控制叶片形状的最优参数组合，能增强混流泵水力性能对流量工况参数随机不确定波动的抗干扰能力和运行的稳定性。为达到以上目的，本发明采用的技术方案是包括以下步骤：(a)选取叶片速度矩分布函数f(x)以及函数f(x)中的n个参数x1,x2,…,xn作为优化变量x＝[x1,x2,…,xn]；用均匀试验设计方法生成具有n个参数、m个水平数的均匀试验表，绘制出m个叶轮三维实体图；(b)在流量工况参数的变动区间[ql，qu]内选取q个流量工况点，组成m×q组数据，ql是流量工况参数变动区间的下限值，qu是流量工况参数变动区间的上限值；(c)建立鲁棒性衡量函数在鲁棒性衡量函数r(x)的基础上将鲁棒优化问题分解为鲁棒子优化函数r(x)和最优子优化函数o(x)：minimizer(x)＝(r1(x),r2(x),…,rk(x))，maximizeo(x)＝(f1(x),f2(x),…,fk(x))，h为在干扰邻域bδ内进行随机抽样的规模数；ξj为抽样规模h中的第j个抽样样本；norm(ξj)为关于ξj的无穷向量范数；x∈ω，ω为可行解空间，minimize为求鲁棒子优化函数r(x)的最小值；rk(x)为第k个优化指标在解空间中变化区域的大小，maximize为o(x)的最大值；fk(x)为第k个优化指标的函数值；(d)对(r1(x),r2(x),…,rk(x))和(f1(x),f2(x),…,fk(x))分别进行加权处理，得到对应的鲁棒性函数robust(x)和最优性函数optimal(x)，建立鲁棒优化数学模型minimize(robust(x),optimal(x))；(e)针对步骤(b)中m×q组数据，依据鲁棒优化数学模型minimize(robust(x),optimal(x))得到每一组数据所对应的robust(x)和optimal(x)，建立优化变量x＝[x1,x2,…,xn]与robust(x)和optimal(x)之间的函数关系式并求解该关系式，得到同时使robust(x)和optimal(x)的函数值达到最小化时所对应的一组最优解本发明采用上述技术方案后具有的有益效果为：(1)本发明基于混流泵叶片二元设计理论，通过采用cfd分析技术、多元回归分析和多目标遗传算法(简称“moga”)，将鲁棒优化问题分解为鲁棒子优化问题和最优子优化问题来求解，使整个优化结果更接近于鲁棒优化指标。(2)在建立鲁棒优化数学模型的过程中，充分考虑了混流泵在实际运行过程中流量工况参数随机不确定性波动的干扰因素，并将其考虑在建立优化变量与优化目标函数关系式的过程中，使优化后混流泵的水力性能得到提升的同时也提高了该性能在流量工况参数发生随机不确定波动环境中的稳定性。(3)在鲁棒优化数学模型的建立过程中对所有的数据进行了标准化处理，以此消除了数量级之间的差异，提高了计算结果的准确性。(4)经过本发明的鲁棒优化后，混流泵的水力性能及该性能对实际环境中不确定因素的抗干扰能力得到了提升，这对泵在节能降耗、节约成本方面具有重要的意义。附图说明图1为混流泵叶轮的叶片轴面投影图；图2为混流泵叶轮的三维叶片及其骨线示意图；图3为混流泵叶轮的三维实体图；图4为混流泵叶轮非结构网格划分示意图；图5为本发明设计流程图。图中，6：叶片进口边；7：叶片出口边；1-1＇：叶轮轮缘处流线；5-5＇：叶轮轮毂处流线；2-2＇、3-3＇、4-4＇：轮缘与轮毂间的中间流线；s0：叶片进口边到出口边之间任意一条流线的弧线总长度；s：从任意一条流线的进口到该流线上任意一点的弧线长度；8：三维叶片的进口；9：三维叶片的出口；10：三维叶片的轮缘；11：三维叶片的轮毂；12：三维叶片上的叶片骨线；13：混流泵叶轮进口；14：混流泵叶轮前盖板；15：混流泵三维叶片。具体实施方式本发明首先以混流泵叶轮叶片二元设计理论为基础，选取控制速度矩变化规律的任一种速度矩分布函数里面的n(n≥1)个参数作为优化变量来间接控制叶片的骨线形状；其次，根据鲁棒性的定义，即解的抗干扰能力或解的稳定性，建立针对混流泵水力性能的鲁棒优化数学模型；最后，根据建立的鲁棒优化数学模型，用多元回归分析方法来建立优化目标与优化变量之间的函数关系式，并结合多目标遗传算法(moga)对上述函数关系式进行求解，最终，得到一组满足使混流泵水力性能达到最优化的同时还可以使该水力性能在不确定流量工况参数扰动下的稳定性达到最优化的n个参数的组合。具体步骤如下：(1)依据《现代泵理论与设计》(关醒凡编著，中国宇航出版社，2011)中混流泵叶片的二元设计理论选取速度矩分布函数。如图1和图2所示，混流泵叶轮的叶片的进口边6到出口边7之间有5条流线，分别是叶轮轮缘处流线1-1＇、叶轮轮毂处流线5-5＇以及轮缘与轮毂间的中间流线2-2＇、3-3＇、4-4＇。以其中任意一条流线上任意一点k为例，k点处的速度矩vur有9种分布形式，且通常用分段函数法、指数法、经验公式法、多项式法、反对称二次抛物线旋绕分步法等速度矩分布函数来控制速度矩vur的分布规律，以此来间接控制如图2所示的叶片上的叶片骨线12。图2中，三维叶片由进口8、出口9、轮缘10和轮毂11组成，在进口8、出口9之间是叶片骨线12。速度矩vur的分布规律可以描述为式(1)所示：vur＝vu1r1+f(x)*δ(vur)(1)式中，vu1r1为叶片进口边6处的速度矩；f(x)为选取的vur分布规律的速度矩分布函数中的任意一种函数，x表示任意一条流线的相对长度，x＝s/s0；以附图1中的中间流线3-3＇为例，s表示沿着中间流线3-3＇从叶片轴面投影的进口边6到该流线上任意一点k之间的弧线长度，s0表示沿着中间流线3-3＇从叶片轴面投影的进口边6到叶片轴面投影的出口边7之间的弧线总长度；△(vur)为叶片出口边7与叶片进口边6之间的速度矩差值，△(vur)可以表示为式(2)所示：δ(vur)＝vu2r2-vu1r1(2)式中，vu2r2为叶片出口边7处的速度矩，vu1r1为叶片进口边6处的速度矩。(2)选取函数f(x)中的n(n≥1)个参数x1,x2,…,xn作为优化变量x＝[x1,x2,…,xn]，并依据优化过程的实际要求给定优化变量x中每个参数x1,x2,…,xn的取值范围：其中，xn表示所选取的第n个参数，x1、分别表示参数x1取值范围的下限值和上限值，x2、分别表示参数x2取值范围的下限值和上限值，xn、分别表示参数xn取值范围的下限值和上限值。(3)根据确定的优化变量x＝[x1,x2,…,xn]及其包含的每个参数x1,x2,…,xn的取值范围，用常规的均匀试验设计的方法生成具有n个参数(因素)、m个水平数的均匀试验表，m≥1，即表示该均匀试验表有n+1列、m+1行。均匀试验表中的因素的个数与所选参数的个数n是相等的，水平数m是由均匀试验设计方法中均匀试验表的使用表确定的(参见方开泰编著的《正交均匀试验设计》(科学出版社，1994))。在生成的均匀试验表中，除去第一行和第一列外其余每一行中的n个数值可以用三维建模软件pro/engineer间接绘制出一个叶轮，如此可以得到m个叶轮，叶轮的三维实体图如图3所示，叶轮包含有混流泵叶轮进口13、混流泵叶轮前盖板14和混流泵三维叶片15。以n＝3、m＝5为例，用均匀试验设计方法生成的均匀试验表的式样如下表1所示，该均匀试验表式样有3+1列、5+1行，其中第2列、第2-6行中的j11、j13、j15、j12、j14表示在参数x1的取值范围中的取值，但j11、j13、j15、j12、j14彼此之间各不相等，即但j11≠j13≠j15≠j12≠j14；同理，第3列、第2-6行中的j25、j22、j24、j21、j23表示在参数x2的取值范围中的取值，但j25、j22、j24、j21、j23彼此之间各不相等，即但j25≠j22≠j24≠j21≠j23；第4列、第2-6行中的x33、x35、x31、x34、x32表示在参数x3的取值范围中的取值，但j33、j35、j31、j34、j32彼此之间各不相等，即因此，在表l第2-6行、2-4列中每一行中的三个数值(j11，j25，j33)、(j13，j22，j35)、(j15，j24，j31)、(j12，j21，x34)、(j14，j23，j32)可以用三维建模软件pro/engineer间接绘制出一个叶轮，如此得到5个不同的叶轮。表1均匀试验表式样mx1x2x31j11j25j332j13j22j353j15j24j314j12j21j345j14j23j32(4)采用网格划分软件icem对步骤(3)中得到的m个叶轮分别进行如图4所示的非结构网格划分。针对每一个叶轮，在流量工况参数的变动区间[ql，qu](ql表示流量工况参数变动区间的下限值，qu表示流量工况参数变动区间的上限值)内用蒙特卡罗取样法选取q个流量工况点，如此组成m×q组数据；用ansyscfx软件分别对该m×q组数据进行数值模拟计算，得到m×q组数值模拟结果。以q＝3(即取q1、q2、q3三个流量工况点，且满足q1∈[ql，qu]、q2∈[ql，qu]、q3∈[ql,qu])为例，结合上表1均匀试验表，得到下表2。如表2中第1-5列、2-6行数据所示，每一个叶轮对应三个流量工况点q1、q2、q3，如此得到5×3＝15组数据，用ansyscfx软件分别对该15组数据进行数值模拟计算，得到15组数值模拟结果。表2m(5)×q(3)＝15组数据(5)用水力性能在解空间中变化区域的函数r(x)的大小来衡量其鲁棒性的好坏。r(x)越小表示水力性能的鲁棒性越好，r(x)越大表示水力性能的鲁棒性越差。r(x)的表达式如式(3)：其中，x＝[x1,x2,…,xn]表示步骤(2)中的优化变量；f(x)表示混流泵水力性能优化指标的函数值，水力性能的优化指标可以是泵的扬程h、效率η、高效区范围he等，选择何种优化指标要根据设计要求而定；ω表示可行解空间；δ＝(δ1,δ2,…,δn)表示干扰向量；bδ为以干扰向量δ为半径的干扰邻域；ξ为bδ中的子向量。bδ的定义如式(4)所示：bδ＝{ξ|ξ＝(ξ1,ξ2,…,ξn),ξn∈[-δn,δn]}(4)为了利于鲁棒性衡量函数r(x)的计算，用蒙特卡罗积分法(montecarlointegral)对上述式(3)进行近似得到式(5)：式中，h为在干扰邻域bδ内进行蒙特卡罗随机抽样的规模数；ξj为抽样规模h中的第j个抽样样本；norm(ξj)为关于ξj的无穷向量范数，norm(ξj)的求解(根据方保镕编著的《矩阵论》(清华大学出版社，2004)进行)。(6)在步骤(5)给出混流泵水力性能鲁棒性衡量函数r(x)的基础上，将混流泵水力性能的鲁棒优化问题分解为两个子优化问题，即鲁棒子优化函数r(x)和最优子优化函数o(x)。鲁棒子优化函数r(x)的定义式如式(6)所示：minimizer(x)＝(r1(x),r2(x),…,rk(x))，x∈ω(6)式中，minimize表示求鲁棒子优化函数r(x)的最小值；k表示混流泵水力性能优化指标的个数，优化指标可以是泵的扬程h、效率η、高效区范围he等，选择何种优化指标要根据设计目的而定；x＝[x1,x2,…,xn]表示步骤(2)中的优化变量；ω表示可行解空间；rk(x)表示第k个优化指标在解空间中变化区域的大小，其值由式(5)计算得到。最优子优化函数o(x)的定义式如式(7)所示：maximizeo(x)＝(f1(x),f2(x),…,fk(x))，x∈ω(7)式中，maximize表示求最优子优化函数o(x)的最大值；k表示混流泵水力性能优化指标的个数，优化指标可以是泵的扬程h、效率η、高效区范围he等，选择何种优化指标要根据设计目的而定；x＝[x1,x2,…,xn]表示步骤(2)中的优化变量；ω表示可行解空间；fk(x)表示第k个优化指标的函数值。(7)将两个子优化函数r(x)和o(x)中的(r1(x),r2(x),…,rk(x))和(f1(x),f2(x),…,fk(x))分别进行加权处理，分别得到一个鲁棒性函数robust(x)和一个最优性函数optimal(x)。robust(x)和optimal(x)的表达式如式(8)所示：式中，ωi和μi分别代表第i个优化指标的fi'(x)和ri'(x)所对应的权重因子，ωi和μi可由模糊判断矩阵理论确定或者根据经验给出，且两者满足和的条件。fi'(x)和ri'(x)是分别对fi(x)和ri(x)进行标准化处理得到的，标准化的公式如式(9)所示：式中，fi,max和fi,min分别是fi(x)在可行解空间中的最大值和最小值，ri,max和ri,min分别是ri(x)在可行解空间中的最大值和最小值。最终，由鲁棒性函数robust(x)和最优性函数optimal(x)得到要建立的鲁棒优化数学模型，如式(10)所示：minimize(robust(x),-optimal(x))(10)式中，x＝[x1,x2,…,xn]表示步骤(2)中的优化变量，minimize表示求两个函数robust(x)和-optimal(x)的最小值。(8)依据步骤(7)中所建立的鲁棒优化数学模型minimize(robust(x)，-optimal(x))，分别针对步骤(4)中得到的m×q组数值模拟结果进行分析，得到每一组所对应的robust(x)和-optimal(x)，如此构成m×q组初始样本空间，m×q组初始样本空间与m×q数值的区别在于m×q组初始样本空间还包含了m×q所对应的robust(x)和-optimal(x)。以表2为例，对5×3＝15组数值模拟结果进行分析得到每一组所对应的robust(x)和-optimal(x)，如此形成下表3，下表3中的t1-t15、w1-w15是对15组数值模拟结果分析得到的具体数值：表3初始样本空间(9)根据步骤(8)中构建的m×q组初始样本空间，选取m×q组初始样本空间中的(m×q－10)组数据来用多元回归分析方法分别建立优化变量x＝[x1,x2,…,xn]与robust(x)和-optimal(x)之间的函数关系式，用剩余的10组数据来检验所建函数关系式的预测精度，得到预测结果robust(预测)和-optimal(预测)。多元回归分析方法是指将一个变量视为因变量，将一个或多个变量视为自变量，进而建立起因变量与自变量之间线性或非线性的函数关系式的一种方法。在本发明中，因变量为robust(x)和-optimal(x)，自变量即为作为优化变量的x＝[x1,x2,…,xn]。在多元回归分析中，以表3中的数据为例，分别选中因变量robust(x)在第6列、2-6行的数据和优化变量x＝[x1,x2,…,xn]在第3-5列、2-6行的数据，就可以得到因变量robust(x)与优化变量x之间的函数关系式；同理，分别选中因变量-optimal(x)在第7列、2-6行的数据和优化变量x＝[x1,x2,…,xn]在第3-5列、2-6行的数据，就可以得到因变量-optimal(x)与优化变量x之间的函数关系式。最终，用多元回归分析法所建立的优化变量x与robust(x)和-optimal(x)之间的函数关系式分别如式(11)、式(12)所示：式中，n为步骤(2)中所选取的作为优化变量的参数的个数，h0、ht、htt、htw及l0、lt、ltt、ltw为两个函数关系式中的系数，该些系数的值可从用excle软件中回归分析法得到的系数表中读取得到(10)根据步骤(9)两个函数关系式(11)和(12)，采用多目标遗传算法(moga)对函数关系式进行求解。在使用moga工具箱的过程中，在软件的特定位置处输入优化变量x＝[x1,x2,…,xn]中每个参数的下限值lowlimits＝[x1,x2,...,xn]和上限值采用二进制的编码方式在每个参数的下限值lowlimits＝[x1,x2,...,xn]和上限值之间随机产生200个个体，并将其作为父代种群；通过选择、变异和复制算子产生子代种群，选择算法采用tournament、交叉算法采用intermediate；交叉和变异概率分别设置为0.8和0.2，最优前端个体系数pareto-fraction设置为0.3，最大进化代数设置为500。在求解过程中，对比用步骤(9)中所建立两个函数关系式的预测结果robust(预测)、-optimal(预测)与用cfd的计算结果robust(cfd)、-optimal(cfd)间的相对误差当ε1和ε2均满足≤5％时，输出当前的寻优结果；若ε1和ε2不满足≤5％的条件时，将当前的结果添加到matlab寻优系统的数据空间中，按照步骤(9)重新建立化变量x＝[x1,x2,…,xn]与robust(x)和-optimal(x)之间的函数关系式，直至满足上述≤5％的收敛标准。最终，得到同时使式(11)和式(12)的函数值达到最小化时所对应的一组解即一组满足使混流泵水力性能达到最优化的同时还可以使该水力性能在不确定流量工况参数扰动下的稳定性达到最优化的n个参数的组合。(11)根据步骤(10)寻优得到的一组解即一组满足使混流泵水力性能达到最优化的同时还可以使该水力性能在不确定流量工况参数扰动下的稳定性达到最优化的n个参数的组合，用三维建模软件pro/engineer间接对其建模得到新的混流泵叶轮。以下提供本发明的一个实施例：实施例本实施例中所采用的混流泵的主要设计参数为：比转速ns＝520.17，额定流量qd＝4500m3/h，额定扬程h＝10.34m，转速n＝735r/min。且在本具体实施方式中，取泵的扬程h和泵的效率η作为混流泵水力性能鲁棒优化问题中的优化指标，即本具体实施方式中的优化指标的个数m＝2。具体实施步骤为：(1)选取一四阶无量纲的速度矩分布函数f(x)来控制速度矩vur的分布规律，f(x)的表达式如式(13)所示：f(x)＝ax4+bx3+cx2+dx+e(13)式中，x表示混流泵叶片轴面投影图中任意一条流线的相对长度，x＝s/s0；以图1中中间流线3-3＇为例，s表示沿着流线3-3＇从叶片进口边6到该流线上任意一点k之间的弧线长度，s0表示沿着流线3-3＇从叶片进口边6到叶片出口边7之间的弧线总长度；a、b、c、d、e为多项式的系数，可以根据给定的条件进行求解，给定的求解条件如式(14)：(2)选步骤(1)中的3个参数x1、x2、x3作为混流泵水力性能鲁棒优化问题中的优化变量x＝[x1,x2,x3]，依据本实施例在优化过程中的实际需求给定每个参数x1、x2、x3的取值范围：0≤x1≤1.2、0.2≤x2≤0.8、0≤x3≤1.2，即x1∈[0,1.2]、x2∈[0.2,0.8]、x3∈[0,1.2]。(3)根据步骤(2)中确定的优化变量x及其所包含的3个参数x1、x2和x3的取值范围，用均匀试验设计法生成具有3个因素、37个水平数的均匀试验表，即表示该均匀试验表有3+1列、37+1行。在生成的均匀试验表中，除去第一行和第一列外其余每一行中的3个数值可以用三维建模软件pro/engineer间接绘制出一个叶轮，如此可以得到37个叶轮，叶轮的三维实体图如图3所示。(4)采用网格划分软件icem对步骤(3)中得到的37个叶轮进行如图4所示的非结构网格划分，针对每一个叶轮，在流量工况参数的变动区间[0.8qd，1.2qd](0.8qd表示流量工况参数变动区间的下限值，1.2qd表示流量工况参数变动区间的上限值)内用蒙特卡罗取样法选取4个流量工况点，如此组成37×4＝148组数据；用ansyscfx软件分别对该148组数据进行数值模拟计算，得到148组数值模拟结果。(5)根据前述鲁棒性的定义，给出衡量混流泵水力性能在流量工况参数发生随机不确定性波动环境中鲁棒性好坏的标准——用水力性能在解空间中变化区域的大小r(x)来衡量其鲁棒性的好坏。r(x)越小表示水力性能的鲁棒性越好，r(x)越大表示水力性能的鲁棒性越差。r(x)的表达式如式(15)所示：式中，h为在干扰邻域bδ内进行蒙特卡罗随机抽样的规模数；f(x)表示优化变量x与混流泵水力性能优化指标之间的函数关系式，本实施例中水力性能的优化指标取为泵的扬程h和效率η；ξj为抽样规模h中的第j个样本；norm(ξj)为关于ξj的无穷向量范数，norm(ξj)的求解(根据方保镕编著的《矩阵论》(清华大学出版社，2004)进行)。(6)在步骤(5)给出混流泵水力性能鲁棒性衡量标准r(x)的基础上，将本实施例中选取的关于泵的扬程h和效率η的鲁棒优化问题分解为两个子优化问题，即鲁棒子优化函数r(x)和最优子优化函数o(x)。泵的扬程h和效率η的鲁棒子优化函数r(x)的定义式如式(16)所示：minimizer(x)＝(rh(x),rη(x))(16)式中，minimize表示求泵的扬程h和效率η的鲁棒子优化问题r(x)的最小值；x＝[x1,x2,x3]表示本实施例步骤(2)中的优化变量；rh(x)表示混流泵扬程h在解空间中变化区域的大小，rη(x)表示混流泵效率η在解空间中变化区域的大小，rh(x)和rη(x)的值由式(15)计算得到。泵的扬程h和效率η的最优子优化函数o(x)的定义式如式(17)所示：maximizeo(x)＝(h(x),η(x))(17)式中，maximize表示求泵的扬程h和效率η的最优子优化问题o(x)的最大值；x＝[x1,x2,x3]表示本实施例步骤(2)中的优化变量；h(x)表示混流泵的扬程值，η(x)表示混流泵的效率值。(7)将泵的扬程h和效率η的鲁棒子优化函数r(x)中的rh(x)、rη(x)及泵的扬程h和效率η的最优子优化函数o(x)中的h(x)、η(x)分别进行加权处理，得到一个鲁棒性函数robust(x)和一个最优性函数optimal(x)。robust(x)和optimal(x)的表达式分别为：式中，μh、μη分别表示r′h(x)和r′η(x)所对应的权重因子，ωh、ωη分别表示h'(x)和η'(x)所对应的权重因子，在本实施例中取ωh＝ωη＝0.5、μh＝μη＝0.5。h'(x)和η'(x)是分别对h(x)和η(x)进行标准化处理得到的，r′h(x)和r′η(x)是分别对rh(x)和rη(x)进行标准化处理得到的。标准化的公式如式(19)和(20)：式中，hmax、hmin和ηmax、ηmin分别表示h(x)和η(x)在可行解空间中的最大值和最小值，rh,max、rh,min和rη,max、rη,min分别表示rh(x)和rη(x)在可行解空间中的最大值和最小值。最终，得到本实施例要建立的鲁棒优化数学模型，如式(21)所示：minimize[robust(x),-optimal(x)](21)式中，x＝[x1,x2,x3]表示本实施例步骤(2)中的优化变量，minimize表示求两个函数robust(x)和-optimal(x)的最小值。(8)依据本实施例步骤(7)中所建立的鲁棒优化数学模型minimize[robust(x),-optimal(x)]，分别对本实施例步骤(4)中得到的148组数值模拟结果进行分析，得到每一组所对应的robust(x)和-optimal(x)，如此构成148组初始样本空间。(9)根据本实施例步骤(8)中构建的148组初始样本空间，选取其中的137组数据来用多元回归分析方法分别建立优化变量x＝[x1,x2,x3]与robust(x)和-optimal(x)之间的函数关系式，用剩余的10组数据来检验所建函数关系式的预测精度。多元回归分析方法是指将一个变量视为因变量，将一个或多个变量视为自变量，进而建立起因变量与自变量之间线性或非线性的函数关系式的一种方法。在本实施例中，因变量为robust(x)和-optimal(x)，自变量即为作为优化变量的x＝[x1,x2,x3]。最终，用多元回归分析法所建立的优化变量x＝[x1,x2,x3]与robust(x)和-optimal(x)之间的函数关系式如式(22)和(23)：式中，h0、ht、htt、htw及l0、lt、ltt、ltw为两个函数关系式中的系数，该些系数的值可从用excle软件中回归分析法得到的系数表中读取得到。最终，将建立好的两个函数关系式(22)和(23)保存到特定的工作目录中。(10)根据步骤(9)中保存到特定的工作目录中的两个函数关系式(22)和(23)，用matlab软件自带的moga工具箱对其进行求解。在使用moga工具箱的过程中，在软件的特定位置处输入优化变量x＝[x1,x2,x3]中每个参数的下限值lowlimits＝[0,0.2,0]和上限值uplimits＝[1.2,0.8,1.2]，采用二进制的编码方式在每个参数的下限值lowlimits＝[0,0.2,0]和上限值uplimits＝[1.2,0.8,1.2]之间随机产生200个个体，并将其作为父代种群；通过选择、变异和复制算子产生子代种群，选择算法采用tournament、交叉算法采用intermediate；交叉和变异概率分别设置为0.8和0.2，最优前端个体系数pareto-fraction设置为0.3，最大进化代数设置为500。其次，在求解过程中，由预测精度得到预测结果计算相对误差，根据相对误差判断是否输出最优解。对比用本实施例步骤(9)中所建立两个函数关系式的预测结果robust(预测)、-optimal(预测)与用cfd的计算结果robust(cfd)、-optimal(cfd)间的相对误差：当ε1和ε2均满足≤5％时，输出当前的寻优结果，即最优解；若ε1和ε2不满足≤5％的条件时，将当前的结果添加到matlab寻优系统的数据空间中，按照本实施例步骤(9)重新建立化变量x＝[x1,x2,x3]与robust(x)和-optimal(x)之间的函数关系，直至ε1和ε2均满足≤5％的收敛标准。最终，得到同时使式(22)和式(23)的函数值达到最小化时所对应的一组解即一组满足使混流泵混流泵扬程h和效率η达到最优化的同时还可以使扬程h和效率η在不确定流量工况参数扰动下的稳定性达到最优化的x1、x2、x3参数的组合。(11)根据步骤(11)寻优得到的一组解即一组满足使混流泵混流泵扬程h和效率η达到最优化的同时还可以使扬程h和效率η在不确定流量工况参数扰动下的稳定性达到最优化的x1、x2、x3参数的组合，用三维建模软件pro/engineer间接对其建模得到新的混流泵叶轮。当前第1页12