水利工程论文-考虑滑移的离心泵叶片S2流面反问题计算方法.doc

水利工程论文-考虑滑移的离心泵叶片S2流面反问题计算方法摘要：针对离心泵叶轮叶片数较少的特点，将有限叶片数带来的滑移的影响引入S2流面反问题计算过程之中，提出了用滑移率修正在假定流动轴对称条件下确定的叶片安放角的做法.从而使S2流面反问题计算方法更接近于实际情况.文中给出了计算实例，并与现有的S2流面反问题计算结果进行了对比.关键词：离心泵叶轮叶片滑移S2流面反问题计算S2流面反问题计算，就是根据流动参数确定叶轮叶片几何形状的过程.反问题计算与水泵的设计密切相关.目前已有的S2流面反问题计算方法，一般是建立在叶片无穷多，即流动轴对称前提条件下的.这种假定对一般的叶轮机械(如透平式压缩机，混流式水轮机等)来讲是可以接受的，但对离心泵而言，由于叶片数一般只有3至6片，这样，轴对称假定与叶轮内部的实际流动状况有较大差异，造成目前已有的S2流面反问题计算方法在离心泵上的应用效果不理想，以至于离心泵叶轮的水力设计方法还是主要采用以经验系数为主的一元设计理论.这样，很难从根本上解决水泵设计中存在的本质问题.因此，对S2流面反问题计算方法进行研究，就显得尤为重要.为此，本文在已有研究结果的基础上，将滑移理论引入到S2流面的反问题计算过程之中，提出了一种具有工程实用价值的S2流面反问题计算方法.一、S2流面反问题计算基本公式1.1基本运动方程取S2流面上流线的轴面投影为正交曲线坐标的q1轴(见图1)，取与q1垂直的方向为q2轴，圆周切线方向为q3轴.同时假定流体为理想液体，在轴对称流动的情况下，可写出沿q-2方向的运动方程式1：(1)图1研究叶轮的准正交曲线坐标这里，下标i(i=1,2,3)对应q1,q2,q3三个坐标分量，H1为拉梅系数.E1为叶片进口处单位质量液体的能量，E1=p1/+v21/2.K1为进口处速度矩，K1=(vur)1.为叶轮旋转角速度.q1，q2是指流面偏导数，就是流场参数沿相对流面S2变化时的偏导数，按复合函数求导法则确定.考虑下式：(2)式(1)可写成：并引入F(q1,q2)从而，沿q2方向的运动方程可写成：(以下依次为(3),(4),(5)(3)(4)(5)式(5)是求解S2流面反问题的基本公式.1.2滑移率的计算公式在离心泵叶轮内，因叶片数相对较少，因此产生轴向旋涡.轴向旋涡使流线偏离叶片骨线2，这就是一般文献中所说的离心泵的滑移.对滑移的研究，目前主要集中在叶轮出口处，即通过滑移系数计算叶轮出口处的绝对速度圆周分量的差值，而对叶片进口处和叶片中部的滑移，研究得不多.为此，笔者在文献3中提出了一种确定叶片中部滑移量的方法.在该方法中，首先引入了“滑移率”的概念，以表示，其物理意义为：在叶片表面(骨线)上某一P点处，叶片安放角(K)p，与液流角()p的正切值之差为滑移率()p，即：()p=tg(K)p-tg()p(6)在进行反问题计算时，由于K是未知的，因此，不可能通过式(6)来计算()p，而根据笔者在文献3中的推导，有：()p=(mp-m)/(m-m)(7)式中，mp为轴面流线在叶片P点处的长度值；m和m分别为叶片进口处和叶片出口处的轴面流线长度值；为叶片出口处的滑移率.在叶片S2流面反问题计算过程，可根据泵的扬程H、滑移系数，以及其它参数，通过迭代计算得出：=(vm)/(u-(vu+u)-(vm)/(u-vu)(8)式中，为按Weisner方法计算得到的滑移系数，1.3叶片型线的计算根据在水力机械中广泛使用的叶片型线微分方程式，有3：(d)/(dm)=-(1)/(rtgk)(9)将式(6)代入式(9)中，并根据速度三角形，将tg用vm/(r-v3)代替，有：(d)/(dm)=(v3r-r2)/(vmr2-(v3r)r+r)(10)考虑到dm=H1dq1，式(5)/H1q1可近似用式(9)中的d/dm来代替，同时注意到在=常数的特征线上有d=0，这样，式(5)可写为：(d(v3r)/(dq2)=F(q1,q2)(vmr2)/(v3r-r2-r)(11)从而，通过求解式(10)和式(11)便可完成S2流面反问题的计算.二、S2流面反问题计算的基本方法2.1F(q1,q2)的确定从式(11)中可以看出，F(q1,q2)与轴面流场有关，主要取决于vur沿q2方向的变化率.在反问题计算时，有两种方法可以用于确定F(q1,q2)：一是直接给定F(q1,q2)的分布，二是通过正问题的求解结果计算得出.在第一种确定F(q1,q2)的方法中，一般是先给定叶片进口边处的F(q1,q2)和叶片出口处的F(q1,q2)，然后给定分布函数G，从而确定叶片上任一点处的F(q1,q2)4：F(q1,q2)=F(q1,q2)+FII(q1,q2)-F1(q1,q2)G.(12)由于F(q1,q2)是一个二维函数，所以无论是对函数G，还是对F(q1,q2)和F(q1,q2)，都是很难给定的，因此，这种方法较少使用.在第二种方法中，一般是先假定(或给定)轴面流场，然后通过下式计算F(q1,q2)1：F(q1,q2)=vm()/(m)(1)/(cos)-(d)/(dl)(1)/(cos)+(sin)/(r)+(lnk)/(m)sin-(dvm)/(dl