本发明公开了一种双叶片无堵塞离心叶轮水力设计方法，其中双叶片无堵塞离心叶轮水力结构包括前盖板、叶片和后盖板；1）基于速度系数法求解双叶片无堵塞离心叶轮主要水力几何参数；2）基于二次贝塞尔曲线控制算法对叶片轴面形状进行控制；3）基于四次贝塞尔曲线控制算法对叶片型线安放角进行控制；4）基于NSGA II遗传优化算法对双叶片离心叶轮进行优化设计；5）根据上述步骤1）—4）中所得到的数据作为双叶片无堵塞离心叶轮的尺寸参数，制作符合给定的流量、扬程、泵效率和转速值的双叶片无堵塞离心叶轮。本发明设计的双叶片无堵塞离心叶轮能够提高叶轮的过流和抗缠绕能力，改善叶轮内部流动状况，提高叶轮的运行稳定性和水力效率。
1.一种双叶片无堵塞离心叶轮水力设计方法，其中双叶片无堵塞离心叶轮水力结构包括前盖板(1)、叶片(2)和后盖板(3)；其特征在于，包括如下步骤：1)根据给定的流量、扬程、泵效率和转速值，基于速度系数法求解双叶片无堵塞离心叶轮水力几何参数，双叶片无堵塞离心叶轮水力几何参数包括：叶轮轮毂直径Dh、叶轮进口当量直径D1、叶片进口宽度b1、叶轮出口直径D2、叶片出口宽度b2、叶片进口安放角β1、叶片出口安放角β2、叶片包角φ、叶片前缘后掠角θ；2)基于二次贝塞尔曲线控制算法对叶片轴面形状进行控制，通过调整叶片前盖板(1)的轴面投影线和叶片后盖板(3)的轴面投影线对叶片的过流面积进行合理控制；3)基于四次贝塞尔曲线控制算法对叶片型线安放角进行控制，在叶片型线绘制过程中通过调整叶片前盖板型线叶片安放角控制线、叶片中间型线叶片安放角控制线和叶片后盖板型线叶片安放角控制线对叶片的型线进行合理控制，其中叶片中间型线指的是位于叶片前盖板和叶片后盖板中间的叶片型线；4)基于NSGA II遗传优化算法对双叶片离心叶轮进行优化设计，通过数值计算的方法或试验手段对叶轮的关键几何参数进行优化，其中所述的关键几何参数为：叶轮出口直径D2、叶片出口宽度b2、叶片出口安放角β2，提高双叶片无堵塞离心叶轮的整体性能；5)根据上述步骤1)—4)中所得到的数据作为双叶片无堵塞离心叶轮的尺寸参数，制作符合给定的流量、扬程、泵效率和转速值的双叶片无堵塞离心叶轮。
2.如权利要求1所述的一种双叶片无堵塞离心叶轮水力设计方法，其特征在于：所述的步骤1)中基于速度系数法求解双叶片无堵塞离心叶轮水力几何参数的方法为：通过给定的流量、扬程、泵效率和转速值，对双叶片无堵塞离心叶轮水力设计参数的计算方法如下：
其中式中：d为泵轴直径，m；Mn为计算扭矩，N·m；τ为材料许用切应力，Pa；Pc为计算轴功率，kW；n为叶轮给定转速，r/min；ρ为介质密度，kg/m3；g为重力加速度，N/kg；Q为给定流量，m3/s；H为给定扬程，m；ηp为给定泵效率；
其中式中：K0为叶轮进口当量直径修正系数，值为4.5～5.5；Q为给定流量，m3/s；n为叶轮给定转速，r/min；
Km1＝0.000534ns+0.09857 (8)
其中式中：K1为叶片进口宽度修正系数，值为2～2.5；v1m为叶片进口绝对速度轴向分量，m/s；Km1为叶片进口速度轴向分量修正系数；ns为比转速；D1m为叶片平均进口直径，m；Q为给定流量，m3/s；叶轮进口当量直径D1，m；
双叶片叶轮的直径比一般多叶片叶轮直径稍大，其计算公式如下：
其中式中：ns为比转速；g为重力加速度，N/kg；H为给定扬程，m；n为叶轮给定转速，r/min；
对于输送含有长纤维、大颗粒杂质的双叶片离心叶轮，其出口宽度应大于多叶片叶轮，计算公式如下：其中式中：ns为比转速；g为重力加速度，N/kg；H为给定扬程，m；n为叶轮给定转速，r/min；
其中式中：u1为叶片进口圆周速度，m/s；β1m为叶片中间流线进口安放角，°；v1m为叶片进口绝对速度轴向分量，m/s；D1m为叶片平均进口直径，m；n为叶轮给定转速，r/min；
对于双叶片离心叶轮，在叶型绘制中，当取三条型线加以控制时，将这三条型线定义为前盖板型线、中间型线和后盖板型线；其中叶片中间型线的进口安放角根据公式(13)求得，对于前盖板型线的进口安放角β1s和后盖板型线的进口安放角β1h通过公式(15)和(16)求得；
当选择更多的型线,即多于三条型线,对叶片加以控制时，则线性插值的方式求得其他型线的进口安放角；
由于双叶片离心叶轮出口宽度b2加大，故叶片出口安放角比多叶片叶轮要小；12°<β2<
25°时，双叶片叶轮将取得较优的特性曲线，比转速较小时，叶片出口安放角取大值；比转速较大时，取小值；
对于双叶片叶轮，叶片包角过大会增加叶片的摩擦面积，增加水力损失；而叶片包角过小时又降低叶片对流体的控制能力和流动的稳定性；因此为了获得更好的水力性能，对于双叶片离心叶轮的叶片包角取220°～280°；
对于输送含有长纤维、大颗粒杂质的双叶片离心叶轮，叶片前缘设计成后掠形式，便于颗粒和杂质污水污物向外输送，叶片前缘后掠角θ的值为30°～100°。
3.如权利要求1所述的一种双叶片无堵塞离心叶轮水力设计方法，其特征在于：所述的步骤2)中，基于二次贝塞尔曲线控制算法对叶片轴面形状进行控制的方法为：对于双叶片离心叶轮，其轴面由叶片前盖板轴面投影线、叶片中间截面轴面投影线、叶片后盖板轴面投影线、叶片出口边轴面投影线和叶片进口边轴面投影线组成；对于叶片前盖板轴面投影线和叶片后盖板轴面投影线采用二次贝塞尔曲线进行控制，二次贝塞尔曲线的函数见公式(17)；其中叶片前盖板轴面投影线的贝塞尔曲线由给定点P1s、P3s和P5s加以控制，其中给定点P1s、P3s和P5s是根据叶片轴面计算几何参数(D1，D2，ShroudZ2)进行确定，同时采用中介点P2s和P4s对贝塞尔曲线方程的t值加以控制，其中中介点P2s和P4s由t值进行确定；叶片后盖板轴面投影线的贝塞尔曲线由给定点P1h、P3h和P5h加以控制，其中给定点P1h、P3h和P5h是根据叶片轴面计算几何参数(d，D2，ShroudZ2,b2)进行确定，同时采用中介点P2h和P4h对贝塞尔曲线方程的t值加以控制，其中中介点P2h和P4h由t值进行确定；叶片中间截面轴面投影线基于中间线性插值算法根据叶片前盖板轴面投影线和叶片后盖板轴面投影线的具体位置进行自动调整；
B(t)＝(1-t)2P1+2t(1-t)P2+t2P3,t∈[0,1] (17)
式中：B(t)为贝塞尔控制曲线方程表示符号；P1，P2和P3为二次贝塞尔曲线的控制点，控制点根据叶片轴面计算几何参数进行确定；t为贝塞尔曲线控制值。
4.如权利要求1所述的一种双叶片无堵塞离心叶轮水力设计方法，其特征在于：所述的步骤3)中，基于四次贝塞尔曲线控制算法对叶片型线安放角进行控制的方法为：对于双叶片离心叶轮，其叶片型线由叶片前盖板型线、叶片中间型线和叶片后盖板型线组成；所述的叶片前盖板型线、叶片中间型线和叶片后盖板型线所对应的叶片安放角控制线分别为叶片前盖板型线叶片安放角控制线、叶片中间型线叶片安放角控制线和叶片后盖板型线叶片安放角控制线；所述的三条叶片安放角控制线采用四次贝塞尔曲线进行控制，四次贝塞尔曲线的函数见公式(18)；其中叶片前盖板型线叶片安放角控制线的贝塞尔曲线由给定点O1s、O2s、O3s、O4s和O5s加以控制，其中给定点O1s和O5s是根据叶片前盖板型线进出口安放角β1s和β2进行确定,而O2s、O3s和O4s则根据前盖板型线目标安放角的分布规律进行调整；叶片中间型线叶片安放角控制线的贝塞尔曲线由给定点O1m、O2m、O3m、O4m和O5m加以控制，其中给定点O1m和O5m是根据叶片前盖板型线进出口安放角β1m和β2进行确定,而O2m、O3m和O4m则根据中间型线目标安放角的分布规律进行调整；叶片后盖板型线叶片安放角控制线的贝塞尔曲线由给定点O1h、O2h、O3h、O4h和O5h加以控制，其中给定点O1h和O5h是根据叶片后盖板型线进出口安放角β1h和β2进行确定,而O2h、O3h和O4h则根据中间型线目标安放角的分布规律进行调整；
B(t)＝(1-t)4P1+4t(1-t)3P2+8t2(1-t)2P3+4t3(1-t)1P4+t4P5,t∈[0,1] (18)式中：B(t)为贝塞尔控制曲线方程表示符号；P1、P2、P3、P4和P5为四次贝塞尔曲线的控制点，控制点根据叶片进出口安放角进行确定；t为贝塞尔曲线控制值。
5.如权利要求1所述的一种双叶片无堵塞离心叶轮水力设计方法，其特征在于：所述的步骤4)中，基于NSGA II遗传优化算法对双叶片离心叶轮进行优化设计的方法为：对于所述的双叶片离心叶轮，为了提高叶轮的水力效率，通过采用NSGA II遗传优化算法对叶轮的关键几何参数进行优化，由公式(20)和公式(22)知，影响泵效率、扬程的因素为叶轮出口直径D2、叶片出口宽度b2、叶片出口安放角β2和叶片数Z由上式知：H～[D2,b2,Z,β2] (20)
式中：k1H和k2H为扬程系数；D2为叶轮出口直径，m；β2为叶片出口安放角，°；b2为叶片出口宽度，m；Z为叶片数；H为叶轮设计扬程，m；g为重力加速度，N/kg；ρ为介质密度，kg/m3；η为泵效率，％；P2为轴功率，W；
对于双叶片离心叶轮，由于其叶片数Z为2，故效率仅为D2、b2和β2的目标函数；在叶片优化过程中，将效率目标函数定义为F(X)，见公式(23)；
式中：X为变量数组；η为泵效率，％；D2为叶轮出口直径，m；β2为叶片出口安放角，°；b2为叶片出口宽度，m；
针对现有NSGA II在处理多模态问题时存在的局部收敛问题，采用精英控制策略和动态拥挤距离评价方法对原算法进行改进，以更好的对双叶片离心叶轮的目标函数进行求解；精英控制策略是对每一个非支配层的最大个体数目进行约束，来达到控制Pareto最优解层个数的目的，其函数表达式如下：其中式中：Nj为第j个非支配层的最大个体数；N为种群大小；r∈[0,1]，为衰减率；
针对原算法没有考虑每一非支配层上个体分布的均匀性，为改善个体分布的均匀性，采用一种动态的拥挤度评价方法，动态拥挤距离ddc,i表达式如下：式中：dc,j为拥挤距离；Vi表示相邻的两个体之间的差异；Nobj为目前种群个数；
分布表示第i+1和第i-1个体在第m个目标处适应度值；
根据非支配排序和拥挤度理论，对目标函数进行求解并得到控制变量的Pareto最优解集。
一种双叶片无堵塞离心叶轮水力设计方法
技术领域
[0001] 本发明属于离心泵水力设计方法领域，具体涉及一种双叶片无堵塞离心叶轮水力设计方法，主要用于快速有效的对双叶片无堵塞离心叶轮进行水力设计及优化，在保证颗粒通过能力前提下，以进一步提升水力模型的效率，实现节能减排的目的。
背景技术
[0002] 双叶片无堵塞离心叶轮具有良好的无堵塞性能和较高的效率，其被广泛地应用于市政污水处理、工业流程处理等诸多领域。目前市场上的双叶片无堵塞叶轮的效率普遍较低，其水力设计方法较为传统，主要根据经验对叶片进行绘型，无法快速有效的对叶片的轴面形状以及叶片型线的安放角进行控制，同时并没有采用相关优化算法对叶轮进行优化设计。目前，针对双叶片无堵塞离心叶轮水力设计方法的研究主要集中在：1)施卫东在其博士论文《高效无堵塞泵的研究开发及内部流场数值模型》中提出的基于经验系数法的双叶片无堵塞离心叶轮设计方法；在此基础上曹卫东等在其论文《污水泵水力设计方法及统计规律》和朱荣生等在其论文《高效叶片式污水泵叶轮的优化设计》中进一步发展了双叶片无堵塞离心叶轮的经验系数设计法；该经验系数法有效地得到叶轮的基本几何参数，但是设计人员很难对叶片轴面形状和叶片安放角进行快速调整，同时也缺少叶轮设计参数最优化设计的过程；2)在双叶片无堵塞离心叶轮水力优化设计方面，施卫东等在其论文《后掠式双叶片污水泵优化设计与试验》中采用了正交试验法对双叶片叶轮进行优化设计，该方法通过设计正交试验表，对不同参数组合的方案进行试验，通过对比分析得出较优参数组合解，以提高水力模型的效率；但该方法仍有一定的局限性，其最后选取的方案不一定是最优解；3)关于叶片式无堵塞叶轮设计方法公开的专利主要包括“泵叶轮”(专利号：CN 98120840.1)、“一种三流道无堵塞离心泵叶轮设计方法”(专利号：CN 201410146868.2)、“一种多叶片无堵塞泵叶轮及其设计方法”(专利号：CN 201410211176.1)和“一种双叶片冲压式无堵塞叶轮设计方法”(专利号：CN 200910234652.0)等，以上专利大都给出了无堵塞叶轮相关几何参数的计算方法，但是没有给出叶片轴面和叶片型线安放角的控制方法以及叶轮优化设计方法。
[0003] 针对上述存在的问题，本专利采用贝塞尔曲线算法对叶片轴面形状以及叶片型线安放角分布线进行控制，设计人员通过调整控制点对叶片轴面形状以及叶片型线进行控制，进而实现对流道过流面积的控制，以达到稳定叶轮内的流体流动状态和提升效率的目的；再则，通过采用NSGA II优化算法对叶轮关键几何参数进行最优化遍历，以得到效率最佳的参数组合解。因此，该方法具有重要的学术和工程应用价值。
[0004] 经检索，至今尚未见关于该方法的文献和申报专利。
发明内容
[0005] 本发明所要解决的技术问题是，已有双叶片无堵塞离心叶轮设计方法存在以下几类缺点：1)过分依赖设计人员的经验以及相关经验系数的选取；2)无法快速有效的对叶片轴面形状及叶片型线安放角进行控制；3)叶轮水力优化设计流程的缺失或具有一定的局限性。本发明的目的是提供一种双叶片无堵塞离心叶轮水力设计新方法，通过对叶片轴面形状和叶片型线安放角的控制以及结合NSGA II遗传优化算法对双叶片无堵塞叶轮进行水力优化设计。
[0006] 为达到上述目的，本发明的技术方案是：
[0007] 一种双叶片无堵塞离心叶轮水力设计方法，其中双叶片无堵塞离心叶轮水力结构包括前盖板、叶片和后盖板；包括如下步骤：1)根据给定的流量、扬程、泵效率和转速值，基于速度系数法求解双叶片无堵塞离心叶轮水力几何参数，双叶片无堵塞离心叶轮水力几何参数包括：叶轮轮毂直径Dh、叶轮进口当量直径D1、叶片进口宽度b1、叶轮出口直径D2、叶片出口宽度b2、叶片进口安放角β1、叶片出口安放角β2、叶片包角φ、叶片前缘后掠角θ；2)基于二次贝塞尔曲线控制算法对叶片轴面形状进行控制，通过调整叶片前盖板的轴面投影线和叶片后盖板的轴面投影线对叶片的过流面积进行合理控制，以优化叶轮的水力性能；3)基于四次贝塞尔曲线控制算法对叶片型线安放角进行控制，在叶片型线绘制过程中通过调整叶片前盖板型线叶片安放角控制线、叶片中间型线叶片安放角控制线和叶片后盖板型线叶片安放角控制线对叶片的型线进行合理控制,以优化叶轮的水力性能，其中叶片中间型线指的是位于叶片前盖板和叶片后盖板中间的叶片型线；4)基于NSGA II遗传优化算法对双叶片离心叶轮进行优化设计，通过数值计算的方法或试验手段对叶轮的关键几何参数进行优化，其中所述的关键几何参数为：叶轮出口直径D2、叶片出口宽度b2、叶片出口安放角β2，提高双叶片无堵塞离心叶轮的整体性能；5)根据上述步骤1)—4)中所得到的数据作为双叶片无堵塞离心叶轮的尺寸参数，制作符合给定的流量、扬程、泵效率和转速值的双叶片无堵塞离心叶轮。
[0008] 所述的步骤1)中基于速度系数法求解双叶片无堵塞离心叶轮水力几何参数的方法为：通过给定的流量、扬程、泵效率和转速值，对双叶片无堵塞离心叶轮水力设计参数的计算方法如下：
[0009] 1.1叶轮轮毂直径Dh
[0010] Dh＝1.5d (1)
[0011]
[0012]
[0013]
[0014] 其中式中：d为泵轴直径，mm；Mn为计算扭矩，N·m；τ为材料许用切应力，Pa；Pc为计算轴功率，kW；n为叶轮给定转速，r/min；ρ为介质密度，kg/m3；g为重力加速度，N/kg；Q为给定流量，m3/s；H为给定扬程，m；ηp为给定泵效率；
[0015] 1.2叶轮进口当量直径D1
[0016]
[0017] 其中式中：K0为叶轮进口当量直径修正系数，推荐值为4.5～5.5；Q为给定流量，m3/s；n为叶轮给定转速，r/min；
[0018] 1.3叶片进口宽度b1
[0019]
[0020]
[0021] Km1＝0.000534ns+0.09857 (8)
[0022]
[0023]
[0024] 其中式中：K1为叶片进口宽度修正系数，推荐值为2～2.5；v1m为叶片进口绝对速度轴向分量，m/s；Km1为叶片进口速度轴向分量修正系数；ns为比转速；D1m为叶片平均进口直径，m；Q为给定流量，m3/s；叶轮进口当量直径D1，m；
[0025] 1.4叶轮出口直径D2
[0026] 双叶片叶轮的直径比一般多叶片叶轮直径稍大，其计算公式如下：
[0027]
[0028] 其中式中：ns为比转速；g为重力加速度，N/kg；H为给定扬程，m；n为叶轮给定转速，r/min；
[0029] 1.5叶片出口宽度b2
[0030] 对于输送含有长纤维、大颗粒等杂质的双叶片离心叶轮，其出口宽度应大于多叶片叶轮，计算公式如下：
[0031]
[0032] 其中式中：ns为比转速；g为重力加速度，N/kg；H为给定扬程，m；n为叶轮给定转速，r/min；
[0033] 1.6叶片进口安放角β1
[0034]
[0035]
[0036] 其中式中：u1为叶片进口圆周速度，m/s；β1m为叶片中间流线进口安放角，°；v1m为叶片进口绝对速度轴向分量，m/s；D1m为叶片平均进口直径，m；n为叶轮给定转速，r/min；
[0037] 对于双叶片离心叶轮，在叶型绘制中，当取三条型线加以控制时，通常将这三条型线定义为前盖板型线、中间型线和后盖板型线；其中叶片中间型线的进口安放角根据公式(13)求得，对于前盖板型线的进口安放角β1s和后盖板型线的进口安放角β1h通过公式(15)和(16)求得；
[0038] β1s＝β1m-5 (15)
[0039] β1h＝β1m+5 (16)
[0040] 其中式中：β1m为叶片中间流线进口安放角，°；
[0041] 当选择更多的型线(即多于三条型线)对叶片加以控制时，则以通常线性插值的方式求得其他型线的进口安放角；
[0042] 1.7叶片出口安放角β2
[0043] 由于双叶片离心叶轮出口宽度b2加大，故叶片出口安放角比多叶片叶轮要小；通常12°<β2<25°时，双叶片叶轮将取得较优的特性曲线，比转速较小时，建议叶片出口安放角取大值；比转速较大时，建议取小值；
[0044] 1.8叶片包角φ
[0045] 对于双叶片叶轮，叶片包角过大会增加叶片的摩擦面积，增加水力损失；而叶片包角过小时又降低叶片对流体的控制能力和流动的稳定性；因此为了获得更好的水力性能，对于双叶片离心叶轮的叶片包角建议取220°～280°；
[0046] 1.9叶片前缘后掠角θ
[0047] 对于输送含有长纤维、大颗粒等杂质的双叶片离心叶轮，叶片前缘设计成后掠形式，便于颗粒和杂质等污水污物向外输送，叶片前缘后掠角θ的推荐值为30°～100°。
[0048] 所述的步骤2)中，基于二次贝塞尔曲线控制算法对叶片轴面形状进行控制的方法为：对于双叶片离心叶轮，其轴面由叶片前盖板轴面投影线、叶片中间截面轴面投影线、叶片后盖板轴面投影线、叶片出口边轴面投影线和叶片进口边轴面投影线组成；对于叶片前盖板轴面投影线和叶片后盖板轴面投影线采用二次贝塞尔曲线进行控制，二次贝塞尔曲线的函数见公式(17)；其中叶片前盖板轴面投影线的贝塞尔曲线由给定点P1s、P3s和P5s加以控制，其中给定点P1s、P3s和P5s是根据叶片轴面计算几何参数(D1，D2，ShroudZ2)进行确定，同时采用中介点P2s和P4s对贝塞尔曲线方程的t值加以控制，其中中介点P2s和P4s由t值进行确定；叶片后盖板轴面投影线的贝塞尔曲线由给定点P1h、P3h和P5h加以控制，其中给定点P1h、P3h和P5h是根据叶片轴面计算几何参数(d，D2，ShroudZ2,b2)进行确定，同时采用中介点P2h和P4h对贝塞尔曲线方程的t值加以控制，其中中介点P2h和P4h由t值进行确定；叶片中间截面轴面投影线基于中间线性插值算法根据叶片前盖板轴面投影线和叶片后盖板轴面投影线的具体位置进行自动调整。
[0049] B(t)＝(1-t)2P1+2t(1-t)P2+t2P3,t∈[0,1] (17)
[0050] 式中：B(t)为贝塞尔控制曲线方程表示符号；P1，P2和P3为二次贝塞尔曲线的控制点，控制点根据叶片轴面计算几何参数进行确定；t为贝塞尔曲线控制值。
[0051] 所述的步骤3)中，基于四次贝塞尔曲线控制算法对叶片型线安放角进行控制的方法为：对于双叶片离心叶轮，其叶片型线由叶片前盖板型线、叶片中间型线和叶片后盖板型线组成；叶片前盖板型线、叶片中间型线和叶片后盖板型线所对应的叶片安放角控制线分别为叶片前盖板型线叶片安放角控制线、叶片中间型线叶片安放角控制线和叶片后盖板型线叶片安放角控制线；三条叶片安放角控制线采用四次贝塞尔曲线进行控制，四次贝塞尔曲线的函数见公式(18)；其中叶片前盖板型线叶片安放角控制线的贝塞尔曲线由给定点O1s、O2s、O3s、O4s和O5s加以控制，其中给定点O1s和O5s是根据叶片前盖板型线进出口安放角β1s和β2进行确定,而O2s、O3s和O4s则根据前盖板型线目标安放角的分布规律进行调整；叶片中间型线叶片安放角控制线的贝塞尔曲线由给定点O1m、O2m、O3m、O4m和O5m加以控制，其中给定点O1m和O5m是根据叶片前盖板型线进出口安放角β1m和β2进行确定,而O2m、O3m和O4m则根据中间型线目标安放角的分布规律进行调整；叶片后盖板型线叶片安放角控制线的贝塞尔曲线由给定点O1h、O2h、O3h、O4h和O5h加以控制，其中给定点O1h和O5h是根据叶片后盖板型线进出口安放角β1h和β2进行确定,而O2h、O3h和O4h则根据中间型线目标安放角的分布规律进行调整；
[0052] B(t)＝(1-t)4P1+4t(1-t)3P2+8t2(1-t)2P3+4t3(1-t)1P4+t4P5,t∈[0,1] (18)[0053] 式中：B(t)为贝塞尔控制曲线方程表示符号；P1、P2、P3、P4和P5为四次贝塞尔曲线的控制点，控制点根据叶片进出口安放角进行确定；t为贝塞尔曲线控制值。
[0054] 所述的步骤4)中，基于NSGA II遗传优化算法对双叶片离心叶轮进行优化设计的方法为：对于双叶片离心叶轮，为了提高叶轮的水力效率，通过采用NSGA II遗传优化算法对叶轮的关键几何参数进行优化，由公式(20)和公式(22)知，影响泵效率、扬程的因素为叶轮出口直径D2、叶片出口宽度b2、叶片出口安放角β2和叶片数Z
[0055]
[0056] 由上式知：H～[D2,b2,Z,β2] (20)
[0057] ∵
[0058] ∴η～[D2,b2,Z,β2] (22)
[0059] 式中：k1H和k2H为扬程系数；D2为叶轮出口直径，m；β2为叶片出口安放角，°；b2为叶片出口宽度，m；Z为叶片数；H为叶轮设计扬程，m；g为重力加速度，N/kg；ρ为介质密度，kg/m3；η为泵效率，％；P2为轴功率，W；
[0060] 对于双叶片离心叶轮，由于其叶片数Z为2，故效率仅为D2、b2和β2的目标函数；在叶片优化过程中，将效率目标函数定义为F(X)，见公式(23)；
[0061]
[0062] 式中：X为变量数组；η为泵效率，％；D2为叶轮出口直径，m；β2为叶片出口安放角，°；b2为叶片出口宽度，m；
[0063] 针对现有NSGA II在处理多模态问题时存在的局部收敛问题，采用精英控制策略和动态拥挤距离评价方法对原算法进行改进，以更好的对双叶片离心叶轮的目标函数进行求解。精英控制策略是对每一个非支配层的最大个体数目进行约束，来达到控制Pareto最优解层个数的目的，其函数表达式如下：
[0064]
[0065] 其中式中：Nj为第j个非支配层的最大个体数；N为种群大小；r∈[0,1]，为衰减率；
[0066] 针对原算法没有考虑每一非支配层上个体分布的均匀性，为改善个体分布的均匀性，采用一种动态的拥挤度评价方法，动态拥挤距离ddc,i表达式如下：
[0067]
[0068]
[0069] 式中：dc,j为拥挤距离；Vi表示相邻的两个体之间的差异；Nobj为目前种群个数；分布表示第i+1和第i-1个体在第m个目标处适应度值。
[0070] 本发明的有益效果是：1)基于速度系数法快速有效的对双叶片无堵塞离心叶轮的基本几何参数进行计算，构建了基本叶轮轴面形状；2)采用贝塞尔曲线算法对叶片轴面形状和叶片型线安放角进行有效控制，便于设计人员快速的对叶片轴面形状和叶片安放角分布规律进行控制，通过对叶轮的过流面积分布的控制，有效提高了叶轮内流体流动稳定性和效率；3)采用NSGA II遗传优化算法对叶片关键几何参数进行最优化求解，进一步提高了叶轮水力效率。
附图说明
[0071] 图1为双叶片无堵塞叶轮的水力结构图；
[0072] 图2为双叶片无堵塞叶轮轴面几何参数图；
[0073] 图3为叶轮轴面形状贝塞尔曲线控制原理图；
[0074] 图4为叶轮进口到出口过流面积分布示意图；(横坐标为叶片轴面单位距离，纵坐标为叶片轴面过流断面面积)；
[0075] 图5为叶片型线分布示意图；(横坐标为叶片型线X方向圆周半径，纵坐标为叶片型线Y方向圆周半径)；
[0076] 图6为叶片型线安放角贝塞尔控制线分布图；(横坐标为叶片型线单位距离，纵坐标为叶片型线安放角)；
[0077] 图7为带有叶片厚度的双叶片无堵塞叶轮示意图；
[0078] 图8为双叶片无堵塞叶轮左视图。
具体实施方式
[0079] 实施例1
[0080] 下面将结合附图和实施例对本发明作进一步的详细描述。本实施例经验系数和推荐值是通过大量的试验统计数据修正获得的，也就是常规本领域技术人员都在使用的经验值；最终的设计值需通过数值计算结果进行校核。
[0081] 结合图1、图2、图3、图4、图5、图6、图7和图8，一种双叶片无堵塞离心叶轮水力设计方法，其中双叶片无堵塞离心叶轮水力结构包括前盖板1、叶片2和后盖板3，如图1、2所示；包括如下步骤：1)根据给定的流量、扬程、泵效率和转速值，基于速度系数法求解双叶片无堵塞离心叶轮主要水力几何参数，双叶片无堵塞离心叶轮水力几何参数包括：叶轮轮毂直径Dh、叶轮进口当量直径D1、叶片进口宽度b1、叶轮出口直径D2、叶片出口宽度b2、叶片进口安放角β1、叶片出口安放角β2、叶片包角φ、叶片前缘后掠角θ；
[0082] 本实施例所述的步骤1)中基于速度系数法求解双叶片无堵塞离心叶轮主要水力几何参数的方法为：给定的流量65m3/h、扬程40m、泵效率60％和转速值2900r/min，对双叶片无堵塞离心叶轮主要水力设计参数的计算方法如下：
[0083] 1.1轮毂直径Dh
[0084] Dh＝1.5d＝45 (1)
[0085]
[0086]
[0087]
[0088] 其中式中：d为泵轴直径，mm；Mn为计算扭矩，N·m；τ为材料许用切应力,轴材料为45号钢，对应的材料许用切应力为49MPa；Pc为计算轴功率，kW；n为叶轮给定转速，r/min；ρ为介质密度，介质为清水，密度为1000kg/m3；g为重力加速度，N/kg；Q为给定流量，m3/s；H为给定扬程，m；ηp为给定泵效率；
[0089] 1.2叶轮进口当量直径D1
[0090]
[0091] 其中式中：K0为叶轮进口当量直径修正系数，推荐值为4.5～5.5，本实施例K0＝3
4.5；Q为给定流量，m/s；n为叶轮给定转速，r/min；
[0092] 1.3叶片进口宽度b1
[0093]
[0094]
[0095] Km1＝0.000534ns+0.09857＝0.15 (8)
[0096]
[0097]
[0098] 其中式中：K1为叶片进口宽度修正系数，推荐值为2～2.5,本实施例K1＝2.3；v1m为叶片进口绝对速度轴向分量，m/s；Km1为叶片进口速度轴向分量修正系数；ns为比转速；D1m为叶片平均进口直径，m；Q为给定流量，m3/s；叶轮进口当量直径D1，m；
[0099] 1.4叶轮出口直径D2
[0100] 双叶片叶轮的直径通常比一般多叶片叶轮直径稍大，其计算公式如下：
[0101]
[0102] 其中式中：ns为比转速；g为重力加速度，N/kg；H为给定扬程，m；n为叶轮给定转速，r/min；
[0103] 1.5叶片出口宽度b2
[0104] 对于输送含有长纤维、大颗粒等杂质的双叶片离心叶轮，其出口宽度应大于多叶片叶轮，计算公式如下：
[0105]
[0106] 其中式中：ns为比转速；g为重力加速度，N/kg；H为给定扬程，m；n为叶轮给定转速，r/min；
[0107] 1.6叶片进口安放角β1
[0108]
[0109]
[0110] 其中式中：u1为叶片进口圆周速度，m/s；β1m为叶片中间流线进口安放角，°；v1m为叶片进口绝对速度轴向分量，m/s；D1m为叶片平均进口直径，m；n为叶轮给定转速，r/min；
[0111] 对于本实施例的双叶片离心叶轮，在叶型绘制中，当取三条型线加以控制时，通常将这三条型线定义为前盖板型线、中间型线(位于叶片前盖板型线和后盖板型线的中间位置)和后盖板型线；其中叶片中间型线的进口安放角根据公式(13)求得，对于前盖板型线的进口安放角β1s和后盖板型线的进口安放角β1h通过公式(15)和(16)求得；
[0112] β1s＝β1m-5＝20 (15)
[0113] β1h＝β1m+5＝30 (16)
[0114] 其中式中：β1m为叶片中间流线进口安放角，°；
[0115] 当选择更多的型线(即多于三条型线)对叶片加以控制，则通常线性插值的方式求得其他型线的进口安放角；其他型线包括介于前盖板型线和中间型线之间的型线，以及中间型线和后盖板型线之间的型线。
[0116] 1.7叶片出口安放角β2
[0117] 由于双叶片离心叶轮出口宽度b2加大，故叶片出口安放角比多叶片叶轮要小；通常推荐值为12°<β2<25°时，双叶片叶轮将取得较优的特性曲线，比转速较小时，建议叶片出口安放角取大值；比转速较大时，建议取小值；本实施例取18°；
[0118] 1.8叶片包角φ
[0119] 对于双叶片叶轮，叶片包角过大会增加叶片的摩擦面积，增加水力损失；而叶片包角过小时又降低叶片对流体的控制能力和流动的稳定性；因此为了获得更好的水力性能，对于双叶片离心叶轮的叶片包角建议取220°～280°，本实施例取250°；
[0120] 1.9叶片前缘后掠角θ
[0121] 对于输送含有长纤维、大颗粒等杂质的双叶片离心叶轮，叶片前缘设计成后掠形式，便于颗粒和杂质等污水污物向外输送，叶片前缘后掠角θ的推荐值为30°～100°；本实施例取30°。
[0122] 2)基于二次贝塞尔曲线控制算法对叶片轴面形状进行控制，通过调整叶片前盖板1的轴面投影线1-1和叶片后盖板3的轴面投影线1-3对叶片的过流面积进行合理控制，以优化叶轮的水力性能；如图3、4所示；
[0123] 本实施例的步骤2)中，基于二次贝塞尔曲线控制算法对叶片轴面形状进行控制的方法为：对于双叶片离心叶轮，其轴面由叶片前盖板轴面投影线1-1、叶片中间截面轴面投影线1-2(叶片中间指的是位于叶片前盖板和后盖板之间的中间位置)、叶片后盖板轴面投影线1-3、叶片出口边轴面投影线1-4和叶片进口边轴面投影线1-5组成；对于叶片前盖板轴面投影线1-1和叶片后盖板轴面投影线1-3采用二次贝塞尔曲线进行控制，二次贝塞尔曲线的函数见公式(17)；其中叶片前盖板轴面投影线的贝塞尔曲线由给定点P1s、P3s和P5s加以控制，其中给定点P1s、P3s和P5s是根据叶片轴面计算几何参数(D1，D2，ShroudZ2)进行确定，同时采用中介点P2s和P4s对贝塞尔曲线方程的t值加以控制，其中中介点P2s和P4s由t值进行确定；叶片后盖板轴面投影线的贝塞尔曲线由给定点P1h、P3h和P5h加以控制，其中给定点P1h、P3h和P5h是根据叶片轴面计算几何参数(d，D2，ShroudZ2,b2)进行确定，同时采用中介点P2h和P4h对贝塞尔曲线方程的t值加以控制，其中中介点P2h和P4h由t值进行确定；叶片中间截面轴面投影线1-2基于中间线性插值算法根据叶片前盖板轴面投影线1-1的控制点的坐标和叶片后盖板轴面投影线1-3的控制点的坐标进行自动调整。
[0124] B(t)＝(1-t)2P1+2t(1-t)P2+t2P3,t∈[0,1] (17)
[0125] 式中：B(t)为贝塞尔控制曲线方程表示符号；P1，P2和P3为二次贝塞尔曲线的控制点，控制点根据叶片轴面计算几何参数进行确定；t为贝塞尔曲线控制值，为了调整叶片轴面投影线，B(t)是一个动态的函数，其具体函数表达式随着控制点坐标以及t值进行变换。
[0126] 所述的步骤2)中，对于叶片轴面形状的前盖板轴面投影线和后盖板轴面投影线采用更高次的贝塞尔曲线进行控制；同时在对轴面过流面积调整时，应遵循面积均匀过渡的原则，避免过流面积出现局部突升或陡降。
[0127] 3)基于四次贝塞尔曲线控制算法对叶片型线安放角进行控制，在叶片型线绘制过程中通过调整叶片前盖板型线叶片安放角控制线3-1、叶片中间型线(叶片中间型线指的是位于叶片前盖板和后盖板中间位置的型线)叶片安放角控制线3-2和叶片后盖板型线叶片安放角控制线3-3对叶片的型线进行合理控制，以优化叶轮的水力性能；如图5、6所示；
[0128] 所述的步骤3)中，基于四次贝塞尔曲线控制算法对叶片型线安放角进行控制的方法为：对于双叶片离心叶轮，其叶片型线由叶片前盖板型线2-1、叶片中间型线2-2和叶片后盖板型线2-3组成；叶片前盖板型线2-1、叶片中间型线2-2和叶片后盖板型线2-3所对应的叶片安放角控制线分别为叶片前盖板型线叶片安放角控制线3-1、叶片中间型线叶片安放角控制线3-2和叶片后盖板型线叶片安放角控制线3-3；三条叶片安放角控制线采用四次贝塞尔曲线进行控制，四次贝塞尔曲线的函数见公式(18)；其中叶片前盖板型线叶片安放角控制线3-1的贝塞尔曲线由给定点O1s、O2s、O3s、O4s和O5s加以控制，其中给定点O1s和O5s是根据叶片前盖板型线进出口安放角β1s和β2进行确定,而O2s、O3s和O4s则根据前盖板型线目标安放角的分布规律进行调整；叶片中间型线叶片安放角控制线3-2的贝塞尔曲线由给定点O1m、O2m、O3m、O4m和O5m加以控制，其中给定点O1m和O5m是根据叶片前盖板型线进出口安放角β1m和β2进行确定,而O2m、O3m和O4m则根据中间型线目标安放角的分布规律进行调整；叶片后盖板型线叶片安放角控制线3-3的贝塞尔曲线由给定点O1h、O2h、O3h、O4h和O5h加以控制，其中给定点O1h和O5h是根据叶片后盖板型线进出口安放角β1h和β2进行确定,而O2h、O3h和O4h则根据中间型线目标安放角的分布规律进行调整；
[0129] B(t)＝(1-t)4P1+4t(1-t)3P2+8t2(1-t)2P3+4t3(1-t)1P4+t4P5,t∈[0,1] (18)[0130] 式中：B(t)为贝塞尔控制曲线方程表示符号；P1、P2、P3、P4和P5为四次贝塞尔曲线的控制点，控制点根据叶片进出口安放角进行确定；t为贝塞尔曲线控制值。
[0131] 所述的步骤3中，叶片型线叶片安放角控制线能多于三条，增加的叶片型线所对应的进出口安放角通过线性插值的方法获取，型线安放角控制线同样采用四次贝塞尔曲线进行控制；在调整叶片安放角控制线的时候，避免出现局部突升或陡降，以保证流体流动的稳定性，提高效率。
[0132] 4)基于NSGA II遗传优化算法对双叶片离心叶轮进行优化设计，通过数值计算的方法或试验手段对叶轮的关键几何参数进行优化，其中所述的关键几何参数为：叶轮出口直径D2、叶片出口宽度b2、叶片出口安放角β2，提高双叶片无堵塞离心叶轮的整体性能；
[0133] 所述的步骤4)中，基于NSGA II遗传优化算法对双叶片离心叶轮进行优化设计的方法为：对于双叶片离心叶轮，为了提高叶轮的水力效率，通过采用NSGA II遗传优化算法对叶轮的关键几何参数进行优化，由公式(20)和公式(22)知，影响泵效率、扬程的因素为叶轮出口直径D2、叶片出口宽度b2、叶片出口安放角β2和叶片数Z
[0134]
[0135] 由上式知：H～[D2,b2,Z,β2] (20)
[0136] ∵
[0137] ∴η～[D2,b2,Z,β2] (22)
[0138] 式中：k1H和k2H为扬程系数；D2为叶轮出口直径，m；β2为叶片出口安放角，°；b2为叶片出口宽度，m；Z为叶片数；H为叶轮设计扬程，m；g为重力加速度，N/kg；ρ为介质密度，kg/m3；η为泵效率，％；P2为轴功率，W；
[0139] 对于双叶片离心叶轮，由于其叶片数Z为2，故效率仅为D2、b2和β2的目标函数；在叶片优化过程中，将效率目标函数定义为F(X)，见公式(23)；
[0140]
[0141] 式中：X为变量数组；η为泵效率，％；D2为叶轮出口直径，m；β2为叶片出口安放角，°；b2为叶片出口宽度，m；
[0142] 针对现有NSGA II在处理多模态问题时存在的局部收敛问题，采用精英控制策略和动态拥挤距离评价方法对原算法进行改进，以更好的对双叶片离心叶轮的目标函数进行求解。精英控制策略是对每一个非支配层的最大个体数目进行约束，来达到控制Pareto最优解层个数的目的，其函数表达式如下：
[0143]
[0144] 其中式中：Nj为第j个非支配层的最大个体数；N为种群大小；r∈[0,1]，为衰减率；
[0145] 针对原算法没有考虑每一非支配层上个体分布的均匀性，为改善个体分布的均匀性，采用一种动态的拥挤度评价方法，动态拥挤距离ddc,i表达式如下：
[0146]
[0147]
[0148] 式中：dc,j为拥挤距离；Vi表示相邻的两个体之间的差异；Nobj为目前种群个数；分布表示第i+1和第i-1个体在第m个目标处适应度值；
[0149] 所述的步骤4中，在采用NSGA II算法对效率目标函数进行优化时，需要通过数值计算的方法或试验手段对不同变量的水力模型进行计算或试验，以建立效率与关键叶轮几何参数的近似目标函数；通过NSGA II优化后的双叶片叶轮模型最终需要采用试验手段对其进行验证，以保证优化过程中的准确性。
[0150] 5)根据上述步骤1)—4)中所得到的数据见下表，将上述计算得到的几何参数汇总入表1，将表1中的计算数据作为双叶片无堵塞离心叶轮的尺寸参数，制作符合给定的流量、扬程、泵效率和转速值的双叶片无堵塞离心叶轮。如图7和图8所示。
[0151] 表1
[0152]几何参数 设计值 几何参数 设计值
叶轮轮毂直径Dh 45mm 叶轮进口当量直径D1 83mm
叶片进口宽度b1 48mm 叶轮出口直径D2 190mm
叶片出口宽度b2 38mm 前盖板进口安放角β1s 20°
中间进口安放角β1m 25° 后盖板进口安放角β1h 30°
叶片出口安放角β2 18° 叶片包角φ 250°
叶片前缘后掠角θ 30°
