离心泵的课程设计-第七组-胡极诸

1、 Zhejiang Sci-Tech University课 程 设 计 论 文 论文题目: 过 程 流 体 机 械 课 程 设 计 学科专业： 过 程 装 备 与 控 制 工 程 作者姓名： 胡 极 诸 指导教师： 崔 宝 玲 递交日期： 2012年7月3日 21摘 要离心泵广泛应用于疏浚、冶金、化工、水利、土建和环保等行业，随着上述行业的迅速发展，对离心泵的需求量不断增加。因此，提高泵的效率，扩大其运行范围，提高运行安全和可靠性，降低振动和噪声，不仅具有巨大的经济效益，而且对提高我国泵工业发展水平具有重要意义。这就需要对离心泵内部的流场结构和能量损失的发生机理进行深入研究。认识离心泵内部流

2、场的各种方法中，试验测试是最基本和最可信的研究方法。为了了解叶轮内部流动特性，自行设计了叶轮，采用1片普通圆柱形叶片进行内部流动特性数值模拟，得到两种工况下的叶轮内部流场分布，分析了不同流量下的流量扬程和流量效率关系。计算结果表明，扬程随着流量增大而减小，效率随着流量增大而减小，到达最佳工况点之后，效率开始下降。关键词：离心泵；水力设计；内部流场；数值模拟目录第一章 叶轮水力设计11.1设计参数11.1.1性能参数11.1.2设计内容11.1.3设计成果11.2泵的设计参数和结构形式计算和确定11.2.1泵形式及级数的确定11.2.2 泵效率计算21.2.3轴功率确定21.3泵进、出口直径的确

3、定21.3.1泵进口直径选择21.3.2 泵出口直径选择31.3.3. 计算泵进出口流速31.4叶轮主要几何参数的计算和确定31.4.1泵轴传递的扭矩31.4.2泵的最小轴径31.4.3叶轮轮毂直径41.4.4叶轮进口直径Dj的初步计算41.3.5叶片进口直径D1的初步计算41.3.6叶片出口直径的D1初步计算41.3.7叶片进口宽度b1 的初步计算51.3.7叶片进口宽度b2 的初步计算51.3.9 叶片包角 61.3.10叶轮叶片进口安放角61.3.11 叶片进口速度61.4精算叶轮外径71.4.1叶片出口排挤系数71.4.2理论扬程71.4.3静矩71.4.4有限叶片理论扬程修正系数71

4、.4.5无穷叶片理论扬程81.4.6出口轴面速度81.4.7出口圆周速度81.4.8出口直径81.4.9误差比较8第二章 欧拉方程性能计算102.1离心泵的理论扬程102.2 考虑有限叶片数的影响10第三章 计算流体力学方法及数值模拟与分析123.1计算模型的建立123.1.1叶轮模型的建立123.2网格划分133.3 流场计算边界条件153.4数值模拟基本参数163.5数值模拟结果分析163.5.1 fluent 计算163.5.2 叶轮内部流场的压力分布173.5.2 速度分布173.5.3水泵扬程和水力效率的预测173.6总结19参考文献21第一章 叶轮水力设计1.1设计参数1.1.1性

5、能参数流量Q=5m3/h 扬程H=50m转速n=1450r/min出口安装角ß2A=30°叶片数Z=41.1.2设计内容离心泵叶轮的水力计算欧拉方程性能计算CFD性能预测内部特性分析1.1.3设计成果二维图三维图共同得到的性能曲线计算的理论值针对两个工况进行计算分析1.2泵的设计参数和结构形式计算和确定1.2.1泵形式及级数的确定 泵的比转速【1】 （1-1）故泵的水力方案为：单级单吸式离心泵。1.2.2 泵效率计算（1）水力效率【1】 （1-2）（2）机械效率【1】 （1-3） （3）容积效率【1】 （1-4）（4）总效率【1】 （1-5）1.2.3轴功率确定 【1】 （

6、1-6）1.3泵进、出口直径的确定1.3.1泵进口直径选择查阅现代泵技术手册可知泵的进口流速一般取值为3m/s，但是从制造经济性和抗汽蚀性能考虑，在吸入口径小于250mm时，可取吸入口流速vs=1.01.8m/s,所以我们不妨取vs=1.0/s。 （1-7）取进口直径=45mm。1.3.2 泵出口直径选择【1】为减小泵体积和排出管路直径，Dd<Ds,一般有Dd=0.71.0Ds （1-8）取=38mm。1.3.3. 计算泵进出口流速【1】泵进口流速： （1-9）泵出口流速： （1-10）1.4叶轮主要几何参数的计算和确定1.4.1泵轴传递的扭矩【1】 （1-11）其中计算功率 【1】 （

7、1-12）1.4.2泵的最小轴【1】对于45号钢，经调质处理之后HB=241286，=(441539)×105 N/m2,取=441×105 N/m2 ，则最小轴颈【1】 （1-13）取1.4.3叶轮轮毂直径查阅现代泵技术手册可知叶轮轮毂直径必须保证轴孔在开键槽以后有一定的厚度，使轮毂具有一定的强度，一般dh=1.21.4d，在这里我们取K=1.3 【1】 （1-14）取1.4.4叶轮进口直径Dj的初步计算D0可由系数速度系数法求得 【1】 （1-15）K0系数，根据统计资料选取主要考虑效率 K0 =3.54.0兼顾效率和气蚀 K0 =4.04.5主要考虑气蚀 K0 =4.

8、55.5对于穿轴叶轮： 【1】 （1-16）1.3.5叶片进口直径D1的初步计算查阅现代泵技术手册可知K1=(0.71.0)，由于泵的比转速为10.49比较小，故K1应取较大值1.0【1】 （1-17）1.3.6叶片出口直径的D1初步计算【1】 （1-18） 【1】 （1-19）由现代泵技术手册可知，由于比转速小于60，上式计算的应乘以修正系数K，取K为1.3，则【1】 （1-20）【1】 （1-21）1.3.7叶片进口宽度b1 的初步计算 【1】 （1-22）取 。其中为加速系数，其取值范围为，考虑效率为主时取较大值，考虑气蚀性能为主时可取较小值，这里取0.8。1.3.7叶片进口宽度b2 的

9、初步计算【1】 （1-23）其中 【1】 （1-24）但根据上述该公式计算得到的低比转速高速离心泵的叶轮出口宽度b2, 较小，故采用以下经验公式【2】 （1-25）取1.3.9 叶片包角 一般可取85°110°，少数可达150°。如果叶片数z大，则应小一些；如果叶片数z小，则应取大一些。故取130°。1.3.10叶轮叶片进口安放角对于低比转速离心泵，C. 普费莱德纳尔提出了如下的计算公式【1】 （1-26）由此公式推导出叶片进口角【1】（1-27）叶片进口角取28°。1.3.11 叶片进口速度【1】 （1-28）计算变工况条件下的泵进口流速:当

10、=Q=5m3/h时, （1-29）当=1.2Q=6m3/h时， （1-30）1.4精算叶轮外径1.4.1叶片出口排挤系数【1】 （1-31）2 叶片出口真是厚度，通常取2 =2-4mm2 叶轮出口轴面截线与流线的夹角，通常取70°-90°取=4mm =90°。1.4.2理论扬程【1】 （1-32）1.4.3静矩S为叶片轴面投影图中线对旋轴的静距【1】 （1-33）1.4.4有限叶片理论扬程修正系数【1】 （1-34）为修正系数a与泵结构形式有关的经验系数，取a=0.6。 对导叶式压水室 a=0.6 对蜗壳式压水室 a=0.60.85对于环形压水室 a=0.851.

11、0 1.4.5无穷叶片理论扬程【1】 （1-36）1.4.6出口轴面速度【1】 （1-37）1.4.7出口圆周速度【1】（1-38）1.4.8出口直径【1】 （1-39）1.4.9误差比较【1】 （1-40）依上述步骤进行第二次精算可得 符合要求。表1-1 最终参数参数名称计算结果最终取值 泵的最小直径d （mm）1215叶轮轮毂直径dh（mm）19.520叶轮进口直径Dj（mm）4545叶片进口直径D1（mm）4545叶片出口直径D2（mm）376384叶轮进口宽度b1（mm）88叶轮出口宽度b2（mm）5.56 叶片出口安放角ß2A（°）（3030 叶片进口安放角

12、23;1A（°）27.928叶片数Z44叶片进口速度（m/s）0.87，1.05 叶片包角（°）130130第二章 欧拉方程性能计算2.1离心泵的理论扬程【1】假定叶片为无穷多，由方程式 （2-1）式中：F2叶轮出口有效过流面积， （2-2）2 ：在计算时可假定2=0.80.9；R2叶轮外圆半径；u2圆周速度， （2-3）D2出口处研究处点的直径；2A叶片出口安放角30°；计算结果： （2-4）2.2 考虑有限叶片数的影响经过分析与推证，对低比转速叶轮，普氏将滑移系数的计算公式【1】具体总结（2-5）a与泵体结构有关的经验系数；对导叶压液室a=0.6;对涡壳压液室

13、a=0.650.85;对环形压液室a=0.851.O；F封闭周线所围面积， （2-6）T叶片节距， （2-7）计算结果: （2-8）第三章 计算流体力学方法及数值模拟与分析3.1计算模型的建立叶轮流道的实体造型是进行流场数值模拟的基础，模型的精确与否直接影响了计算结果的可信度。本文利用自行设计的叶轮，结合UG软件进行叶片、叶轮的精确三维实体造型，为后续的数值模拟奠定了良好的基础。3.1.1叶轮模型的建立根据叶轮的结构参数，利用UG软件进行叶片、叶轮流道的精确三维实体造型5，生成叶轮的三维实体模型，图3-1是主视图及子午视图，图3-2是叶轮三维造型，图3-3是叶轮流到实体模型图3.1 主视图及子

14、午视图 图3.2 半开式叶轮三维造型实体 A 流道正面 B 流道背面图3.3 叶轮流道实体建模 3.2网格划分网格是CFD模型的几何表达式，也是模拟和分析的载体。网格质量对CFD计算精度和计算效率有重要的影响。对于复杂的模型，网格生成是一个漫长而枯燥的过程，经常需要进行大量的试验才能取得成功。GAMBIT是作为FLUENT求解器的专用前处理软件包，用来为数值模拟生成网格模型。GAMBIT提供了混合网格、结构网格和非结构网格等多种类型的网格。本文设计的叶轮比较简单，进行网格划分时只需要把进口壁和叶片部分划的网格间距不一样即可。利用UG软件生成几何模型保存成STEP格式文件后，导入GAMBIT中进

15、行混合网格划分，其步骤如下：首先把整个模型根据不同的层次分成两个体，得到分别为Z1进口部分和Z2叶片部分两个体。设第一个体Z1,进口部分的网格间距为2.5毫米，得到网格数为57524。模型见图3.4。 图3.4 进口部分网格图设第二个体Z2部分的网格间隔为2.5毫米，得到网格数为416777。模型见图3.5。. 图3.5 叶片部分流道网格图最后得到计算网格总数为474301。得到叶轮网格划分示意图3.6。图3.6 叶轮实体模型网格划分示意图3.3 流场计算边界条件 下面对圆柱形叶片流道叶片的叶轮进行数值计算。其中，叶片数Z=4，设计工况点为转速n=1450r/min,流量Q=5m3/h。对叶轮

16、内部流动而言，主要存在进口边界、出口边界和固壁边界。（1）进口边界对于不可压流动，进口边界取为速度入口边界条件(velocity-inlet)。在计算出轴向速度时，假设进口边界上无切向速度与径向速度，并且轴向速度是均匀分布的，因此根据流量即可计算出进口速度。（2）出口边界在计算区域的出口边界上，流动情况完全由区域内部外推得到，可以认为流动己充分发展，即对上游的流动参数没有影响，故取自由出流边界条件(pressure-out)。（3）固壁边界对于固壁边界条件，因不存在壁面有抽吸或喷射等情况，所以固壁表面取为无滑移边界条件。对壁面的k、，可根据标准壁面函数法进行取值在旋转坐标系中，计算采用segr

17、egated方法，湍流模型k-,模型中的系数均采用默认值。数值模拟及分析3.4数值模拟基本参数进行数值模拟的基本参数如表3-1所示： 表3-1 数值模拟的基本参数参数值单位流量Q5,6m3/h转速n1450r/min出口角30°叶片数Z4 3.5数值模拟结果分析3.5.1 fluent 计算运用fluent软件，导入mash文件，定义边界条件，计算收敛曲线如图在Q=Q，Q=1.2Q的工况下，收敛曲线图如图所示图3.7 Q=5m3/h 中收敛图 图3.8 Q=6m3/h 中收敛图3.5.2 叶轮内部流场的压力分布图3.7和3.8示出叶轮中间截面流场的压力分布图。从图3.9和图3.10可

18、以看出，从进口到出口，示一个静压不断增大的过程。液体水进入泵体后，通过扩散逐步把动能转化为压能。由于粘性流体和惯性力的作用，随着半径增大，压力逐渐增大。由于转速较高，在进口范围内形成低压区，在叶轮出口存在压力较大的区域。 图3.9 Q=5m3/h 中截面静压图 图3.10 Q=6m3/h 中截面静压图3.5.2 速度分布 图3.9,3.10所示水泵中间截面的速度分布图。由图3.8，3.10可以看出，随着叶轮圆周直径的增大，圆周速度增大，则速度随之增大 图3.11 Q=5m3/h 中截面速度分布图 图3.12 Q=6m3/h 中截面速度分布图3.5.3水泵扬程和水力效率的预测3.5.3.1 水泵

19、扬程预测水泵叶轮带着液体旋转时把力矩传给液体，使液体的运动状态发生变化，从而完成了能量的转换。叶轮就是通过叶片把力矩传给液体，使液体的能量增加。在模拟计算三维流场的基础上，通过计算水泵进口面的环量差和转子的轴功，可求得水泵的实际扬程和水泵的效率【6-7】， 在水泵中，规定叶轮的扬程出口总水头与进口总水头之差，叶轮的扬程，如式所示【1】 (3-1)式中：pD , pS泵出口、进口处液体的静压力;vD、 vS泵出口、进口处液体的速度;zD 、zS泵出口、进口到任选的测量基准面的距离。3.5.3.2 水泵的效率预测叶轮的水力效率，如下式所示 【8】 (3-2)式中：h水力效率；Q水泵的流量M叶片正、

20、背面表面受到的绕Z轴的力矩之和，N·m计算结果见表3.2表3.2各工况点计算结果工况点Q=5 m3/hQ=6 m3/h进口压力pD pa332378354298出口压力pS pa484414460135 进口速度vD m/s0.871.05 出口速度vS m/s12.1412.18 扭矩 N·m15.8815.07总效率 %4115通过无限叶片和有限叶片公式（第二章）分别计算得出两个工况点下的理论流量，如表3.3所示。扬程流量、效率流量曲线图如图3.13所示。表3.3扬程计算工况点Q=5m3/hQ=6m3/h实际扬程H m3021无限叶片公式计算结果 m3121有限叶片公式计算结果 m2815图3.13 扬程、效率流量曲线图3.6总结 1）通过对多个工况点数值模拟，就可以得到泵的外特性和内特性，如果计算更过工况点，就可以绘制出水泵扬程曲线图，如果计算多个工况点，就可以绘制出水泵的扬程流量曲线和流量效率曲