第二章 离心通风机的工作原理
§1 通风机的基本方程
由于叶轮是离心通风机的主要部件,当叶轮被原动机拖动以等速ω(s-1)旋转时,叶轮对气体作功,使气体获得了静压Ps和动能。下面主要讨论叶轮对气体作功的大小,即Euler方程。
一. 进出口速度三角形
先讨论一个后向叶轮。假设叶轮为理想叶轮:即有无限多个叶片的叶轮,叶轮内的流动损失不计,并且叶片为无限薄。如图2-1所示形成进出口速度三角形。
图2-1 后向叶片叶轮的速度三角形
在叶片进口1处,假设气流以相对速度无撞击进入叶轮。与叶轮的牵连运动 的反向夹角为β1,那么绝对速度 为:
与 的夹角为α1 。在无限多叶片的假设前提下,叶轮叶片出口2点的相对速度的方向与叶轮叶片的方向一致,与反向夹角为出口处叶片几何角β2A,令
与 的夹角为α2′,(上标′表示无限多叶片的假设条件下的出口值※),那1和2点的速度三角形可以得出:
(2-1)
式中下标m表示轴面速度,在径向叶轮中为径向速度,C1u 常称为叶轮进口的预旋。
二.Euler方程
在上述假设条件下,并令流体是不可压缩的,叶轮中的流动是稳定的,不随时间变化,根据动量矩定律,叶轮所获得的转矩M为气体在叶轮进口处动量矩M1=ρQT C1u r1和出口的动量矩M2=ρQT C’2u r2之差:
M=ρQT( r2 C’2u - r1 C1u) (2-2)
式中QT为叶轮中理论容积流量(M3/·S-1),Mω为原动机传给叶轮功率的大小,而Mω/ρQT表示叶轮对1Kg气体在一秒钟内所做的功:
Mω/ρQT =( r2 C2u′- r1 C1u)
那么,叶轮在无限多叶片假设前提下的理论全压升△P th∞为(简称无限多叶片理论全压):
△P th∞= Mω/ QT =ρ(r2 C2u ′- r1 C1u) (N/ m2) (2-3)
(2-3)式即为欧拉方程。
当为径向进口是, C1u =0,则
△P th∞=ρu2 C2u′=ρu22(1-C2m/u2 ctg β2A) (2-4)
如果通风机的入口有导流器,改变导流器角度,即可改变C1u,从而可以改变风机的全压。
(2-3)式如果用压头H th∞表示为:
H th∞=△P th∞ /γ=1/g(u2 C’2u -u1 C1u) (2-5)
上式表示压头与重度γ等无关。
二. 欧拉方程的物理概念
将进出口速度三角形的关系:
W´22 = u22 + C´22 - 2 u2 C´2 cosα2′= u22 + C´22 - 2 u2 C´2u
W12 = u22 + C12 - 2 u2 C1 cosα1= u12 + C12 - 2 u1 C1u (2-6)
代入方程(2-3)式
△P th∞=ρ/2(u22 - u12)+ ρ/2(W12 – W2´2)+
ρ/2(C2′2 - C12)
(2-7)
很显然,第三项表示气体流经叶轮后动能的增加,即动压升,即用△P th∞,d 表示:
△P th∞,d =ρ/2(C2′2 - C12) (2-8a)
那么,在通风机设计中,应力争△Pd在蜗壳式扩压器等元件中加以回收,将动能转化为静压。在回收静压时,有较大的损失,为此要求在叶轮设计中获得较大的静压。而第一,二项表示叶轮中气体的静压的增加,称为静压升△P th∞,s
△P th∞,s =ρ/2(u22 - u12)+ ρ/2(W12 – W2´2) (2-8b)
其中第一项表示由于离心力的作用使气体静压的增加,该项与圆周速度的平方成正比。在轴流风机中,由于u1=u2 ,这一项不存在。第二项表示叶道截面的变化使气体相对速度降低而转化为静压的增加。
§2 叶轮叶片型式和反作用度
由式(2-6)可知叶轮的理论全压△P th∞由静压和动压组成,为了表示叶轮中动压和静压之间的比值,引入叶轮的反作用度。
一. 叶轮的反作用度
定义叶轮的理论静压升△P th∞,s与全压△P th∞之比,称为叶轮的反作用度(反应度,反动度),用Ω表示(用于在理论,无限多叶片假设下,也可以加下标Ωth∞):
Ω(Ωth∞)=(u22 - u12+ W12 – W2´2)/2(u2 C’2u -u1 C1u)
当C1u =0,时
△P th∞,s =ρ/2(u22– W2´2+ C1m2)
Ω=(u22– W2´2+ C1m2)/2 u2
进而假设C1m =C2m ,则
Ω=(2u2 C´2u–C´2u2)/2 u2 C2u =1- C´2u /2 u2 (2-10)
把速度分量的关系(2-1)式代入:
Ω=1- (u2-C2mCtgβ2A )/2 u2 (2-11a)
Ω=1/2+1/2φ2m Ctgβ2A (2-11b)
式中φ2m=C2m/ u2,可见当β2A一定时,Ω与φ2m成正比。还可以令:
τ= C´2u /u2 =tgβ2A /(tgα2´+tgβ2A ) (2-12a)
Ω= 1-1/2τ (C1u =0, C1m =C2m) (2-12b)
如果叶轮的叶片出口角β2A不同,Ω也就不同。为了分析方便,引入无量纲的压力系数ψ,压力系数为压力与出口圆周速度动压ρ/2u22的比值,例如总压系数ψth∞,tat,或静压系数ψth∞,st, 当径向进口时为:
ψth∞,tat=2△P th∞ /ρu22 =2 C´2u /u2=2τ (2-13a)
ψth∞,st=2△P th∞,s /ρu22 =ψth∞,tatΩ (2-13b)
=2τ-τ2
图2-2表示,当β2A加大时,△P th∞加大,但是反作用Ω减少,即:(由2-10式)
ψth∞,tat=4(1-Ω) (2-13C)
图2-2
图2-2可见ψth∞,s是一条抛物线,在C´2u =0和C´2u =2u2时与横坐标相交。其最大值为1,其时C´2u =2u2,τ=1,β2A=90°,总压的一半是静压,一半是动压。在C´2u =2u2时,ψth∞,s =0,ψth∞,tat=4,即τ=2,只产生动能。当τ=0,ψth∞,tat=0,所以叶片不再有三种形式。
二. 叶片的型式
图2-3为由出口角β2A所决定的三种叶片型式。
图2-3
1. 前向叶片:β2A>90°为前向叶片。分为一般前向叶片和多叶前向叶片两种。
2. 后向叶片: β2A<90°为前向叶片,分为曲线型和直线型两种。
3. 径向叶片: 一般分径向直叶片和径向出口叶片,其β2A=90°。
由图2-3中的速度出口的三角形可以看出:前向叶片β2A>90°,C´2u >u2;后向叶片: β2A<90°, C´2u <u2。对于径向叶片C´2u =u2。根据式(2-4),对于径向进口的欧拉方程:
△ P th∞=ρu22(1-C2m/u2 ctg β2A)
所以在其他条件相同时,前向叶片输出的能量最高,后向叶片最低,径向叶片居中。上式用流量Q,及出口直径D2和叶轮宽度表示:
△P th∞=ρu22-ρQu2 /лD2 b2tg β2A (2-14a)
这一关系表明△P th∞=f(Q)是一直线关系,如图2-4(a)表示。前向叶片当Q增加时,压力△P th∞也增加;径向叶片△P th∞保持不变;而后向叶片,随Q的增加,△P th∞减少。从反作用角度看:后向叶片反作用度最大,效率高,而前向叶片反作用度小,动能大,且流动转折大,损失大,效率低。但是前向叶片能量大,即在尺寸,转速相同时,前向叶片的压力(全压)大。
(a)不同叶片角时的理论特征曲线 (b)功率曲线
图2-4 无限多叶片时的理论特征曲线
关于叶轮的功率N th为:
N th=△P th∞Q =ρu22 Q -ρQ2 u2 /(лD2 b2tgβ2A ) (2-14b)
N t∽∽Q的曲线由图2-4(b)上表示。β2A =90°时为一条直线。β2A>90°是前向叶片, N th∽∽Q的抛物线高于直线;β2A<90°的后向叶片的N th∽∽Q抛物线低于直线。
图2-4(a)中的静压△P th∞,s曲线,在径向进口和C1m =C2m的条件下,可用下式表示:
△P th∞,s=ρ/2(u22– W2´2+C2m2)
由于W2´=C2m/sinβ2A
△P th∞,s=ρ/2(u22–C2m2 /tg2β2A)
=ρ/2{u22–Q 2/(л2D22 b22tg 2β2A)} (2-14c)
即△P th∞,s~Q的曲线仍然是一条抛物线,无论前向或后向叶轮,均为β2A=90°的直线△P th∞,s=ρu22/2下方的抛物线。
§3. 气体在叶轮中的实际流动
这一节将研究有限叶片数对叶轮中气流的影响和叶轮中的流动方程。
一. 作用在叶轮中气体的作用力
图2-5中为流道中的一个微元ds×dn×b所受的力,其中ds为流动方向的微元长度,dn为垂直于ds方向的长度,b为叶轮的宽度。
a 后向叶片
b 前向叶片
图2-5 叶片流道中流体所受的力
作用在流体上的力有以下几种:曲率半径为R的弯曲流道会产生一个垂直于叶片流道的离心力W2/Rdm,由于叶轮的转动产生一个径向的离心力rω2dm,该力的分量rω2dmCosβ垂直于流动方向,还有Corriolis力垂直于流动方向,大小为2ωWdm。
垂直于流动方向的力会产生此方向的压力梯度dp/dn:对于后向叶片:
dp/dn=ρ(W2/R+rω2Cosβ-2ωW )
对于前向叶片:
dp/dn=ρ(W2/R+rω2Cosβ+2ωW )
沿流动方向s的作用力只有离心力的分量ρrω2sinβ,此外还有沿s方向的压力梯度dp/ds,那么沿s方向的动量方程:
WdW/ds= -dp/ρds + rω2sinβ (2-17a)
取sinβds=dr
ρWdW -ρrω2dr +dp =0 (2-17b)
沿s方向积分后:
W2/2 + p/ρ-u2/2 = H´=Const (2-18)
上式为沿叶轮中相对流线s方向的伯努利方程,即沿相对流线上伯努利常数不变。
式(2-18)微分后与(2-16)式比较,得出速度分布的微分方程,对于后向叶轮:
dw/dn = 2ω-W/R (2-19a)
对于前向叶轮:
dw/dn = -(2ω+W/R) (2-19b)
上述各式是分析叶轮中流动的基础。
例如对于后向直叶片,R=∞,上式积分后:
W=Wp + 2nω
式中Wp 为叶片工作面上的速度,n为相邻叶片间的距离。当Wp =0时,Wmax=2nω,流量为Q=n2ωb。令叶片间的速度为等速分布;即dw/dn =0,那么:
2ω-W/R =0 R=W/2ω=r1W/2u1
当W沿流动方向为常数时,R为常数,这就是简单的圆弧叶片,但是这种叶片强烈地向后弯曲,只能产生很小的压力升。
对于一种等压分布的流道,即dp/dn=0的流道:
W2/R + rω2Cosβ-2ωW =0
R= W2/ω(2W-uCosβ).
如果假设W的变化,就可以计算出R,然后可以把叶片造型出来,当2W-uCosβ=0, R趋于无穷大,为直线,从这一点开始,如W=常数,那么uCosβ=常数,R仍为无穷大,那从这一点开始就是直线叶片,这种叶片为Grun对于压气机所提出的一种叶片。
二. 有限叶片数的影响
1. 相对涡流
在有限叶片的情况下,无摩擦的流动假设会引起一 “相对涡流”如图2-6a所示。把流道封闭,当流道旋转时,由于流道内空气的惯性,对于一个与叶轮一起转动的观察者会看到流道中的空气做与叶轮旋转方向相反的旋转,这种流动称为相对涡流。当打开流道后,由于气流与涡流的叠加,在压力面上的速度变低,吸力面上的速度加大。当流量减少时,会在压力面上产生“回流”。
涡流的存在,使叶轮出口处的相对速度发生偏转,而不是叶片切线方向流出,即实际出气角β2小于出口叶片角β2A。
a)封闭流道中的轴向涡流 b)开口叶道中的相对流动
图2-6气体在叶道中的流动
2. 有限叶片数对理论扬程的影响
图2-7为各种型式叶片出口由于有限叶片数的影响而形成的速度三角形用虚线表示,其下标用2表示。可以看出:
1).各种型式的叶轮, 出口气流角β2<β2A。
2).平均出口绝对流速C2< C2´,α2<α2´。
3).对于相对速度W2, 后向叶片W2> W2´而前向叶片
W2> W2´;
4).各种叶轮C2u< C2u´,C2u- C2u´= ΔC2u.
图2-7 有限叶片数对出口速度三角形的影响
由于C2u的减少, 导致了可能达到的无限多叶片理论压力ΔP th∞的降低。有限叶片时的理论压力为ΔP th∞,那么ΔP th <ΔP th∞此时仍然假设流体中没有摩擦力的作用。ΔP th与ΔP th∞的比例,称为减功系数,用ε表示:
ε=ΔP th/ΔP th∞=H th/H th∞=u2 C2u/ u2 C2u´= C2u/C2u´
=1-ΔC2u/C2u (2-20)
ΔC2u也称为是滑移速度.
ε的计算多用半经验公式。
(1).Stodola计算公式:
对于后向叶轮:
ε=1-πsinβ2A /z(1-C2m/u2Ctgβ2A) (2-21a)
其中z为叶片数.
对于径向叶片:
ε=1-π/z (2-21b)
由于stodola在推导时仅考虑了叶片出口的影响,忽略了叶道曲率和径向深度的影响,具有近似性,对于叶道长,叶片数多的后向叶轮较为合适, 而Bruno Eek公式计及了叶道曲率的影响.
(2)Eek公式:
对于前后盘平行的叶轮b1=b2
ε=1/{1+sinβ2Aπ/z[1-(D1/D2)2]} (2-22a)
对于Cm为常数的叶轮(即rb=常数的叶轮)
ε=1/{1+sinβ2Aπ/2z[1-(D1/D2)2]} (2-22b)
对于β2A为20~120的一般通风机:
ε=1/{1+(1.5+1.1sinβ2A)π/z[1-(D1/D2)2]} (2-22c)
ε的计算还有很多,就不一一列举了.
三.考虑压缩性影响的理论全压△P th
在考虑压缩性的影响时, △P th与变化过程有关.在通风机中的等熵过程为参数标准:
△P th=(k/k-1)Rρ1T1※ {(P2tat/P1tat)k-1/k-1}
=(k/k-1)P1,s{(P2tat/P1tat)k-1/k-1} + ρ1/2(C22-C12) (2-23)
令P2t - P1t=△P, H th可以用下式近似表示:
H th=△P/r1(1-△P/2.8 P1tat)
四.气流进口冲角
当气流进入叶轮叶道时,相对速度的方向角β1与叶片进口安装角β1A不一致时,其差角i为进口气流冲角,如图2-8所示:
如果β1A>β1, i =β1A-β1>0为正冲角
如果β1A<β1, i =β1A-β1<0为负冲角
如果在设计工况, β1A=β1 ,即气流无冲击进入,在小流量工况, β1<β1A, i >0。如图2-8所示,(图2-8中下标d的变量为设计工况).在小流量工况,气流在非工作面上会出现分离,而且分离区会向下游发展。
图2-8 进口冲角
在大流量工况,β1A<β1,i<0,此时在工作面也出现分离,在工作面上的分离一般不扩大,这是因为在工作面后缘,附近为气流加速区域,呈顺压梯度,减弱了分离区的发展。
§4.通风机的无因次系数
根据叶轮机械的相似理论,要保证气流经过了两个叶轮后的流动过程相似必须满足几何相似,运动相似和动力相似.在通风机叶轮中由于空气的重度小可以忽略重力的作用.同时由于速度低可以忽略气体弹性力的作用。只考虑粘性力,压力和惯性力,这三个力组成一个力的三角形,当雷诺数相等时,粘性力和惯性力比例相等,且方向相同时,力的三角形相似。所以通风机中只要雷诺数相等就会满足动力相似的条件,可以证明只要满足以上的相似条件,通风机的无因次数就一定会相等。※
(1) 压力系数(P)ψ表示压力与叶轮出口动能ρu22/2之比。压力系数实际上就是欧拉数。例如:
ψth∞=△P th∞/ρu22/2 = 2 C2u′/u2 (2-24a)
ψth=△P th/ρu22/2 = 2 C2u/u2=2εC2u′/u2 (2-24b)
ψ=△P/ρu22/2= 2ηhεC2u′/u2
式中ηh是水力效率。
(2) 流量系数Ф,表示速度三角形相似。同时也表示斯特鲁哈数相等。由于风机的转速是稳定的,进口速度三角形相似是St数必然相等。对于离心式风机:
Ф=Q/u2F2=C2mπD 2b2 /(u 2π/4D22)=4C2m b2/D2u2 (2-25a)
对于轴流式风机:
Ф′=C2m /u2 (2-25b)
其中г=D1/D2,那么Ф′=Ф/(1-г2)
目前研究成果认为,流量系数的定义应利用进口面积更合理.即:
Ф1=Q/ u 1(π/4)D12 (2-25C)
这是因为叶栅的流量由进口面积所决定的.但是目前实用的仍然是Ф。
(3)功率系数λ。对于驱动功率(轴功率)NS:
NS=N e/η=△P Q/1000η (KW)
式中N e为有效功率, η是全压效率.功率系数λ为:
λ=1000 NS/(π/4ρD22u23)=ψφ/η (2-26)
根据无因次量的定义代入u2 =πD2n/60可以得出
△P=0.00137ψρD22n2
Q =0.04108φD23n
N =0.5635×10-7λρD25n3
(4) 比转速ns,转速系数σ和直径系数δ
在实际应用中已发现单靠ψ和φ来定量地估计离心叶轮的主要性能是不够的.还需要有一些与通风机尺寸和转速有关的系数.
比转速ns是一台“单位透平机械”的转速。“单位透平机械”是一台与给定机械完全相似的透平机械, 其尺寸产生1m的压头和1m3/s的流量时使用1单位功率,其转速就是给定机械的比转速ns,根据相似定理:相似通风机只要给定两个物理量的比值就可以求出其他量的比值。例如两个通风机1和2相似,同时ρ1=ρ2,那么给定全压比及流量比Q1/Q2就可以表示为其他量D1/D2,N1/N2,n1/n2等的比值:
D2/ D1 =( Q2/ Q1)1/2 (ΔP1/ΔP2) 1/4 (ρ2/ρ1) 1/4
n2/ n1 =( Q1/ Q2)1/2 (ΔP2 (ρ1/ρ2) 3/4
N2/ N1 =( Q2/ Q1) (ΔP2/ΔP1)
对于单位通风机2,Q2=1,ΔP2=1,n2=
那么n s =Q1/2/ΔP 3/4 (2-27a)
上式为比转速的定义式。
相似通风机的比转速相等,但n s相等的通风机,不一定相似。n s是有量纲的量。如果考虑了ρ的影响,应为:
n s =Q1/2/(ΔP/ρ) 3/4
其单位为1/S(m3/s)1/2(N/m2)-3/4(Kg/m3)3/4=1为无因次量。在工业上使用米——公斤力——秒(MKfs)制,即SI制的n s数值比按MKfs制表示的小5.54倍[(9.087)]
n s用无量纲量表示为:
n s =30Φ1/2/π1/2ρ3/4(Ψ/2)3/4 (2-27b)
对于进口标准状态ρ=1.2Kg/m3
n s =24.827Φ1/2/Ψ3/4 (2-27C)
比转速n s可以用来划分通风机的类型。n s大的通风机流量大,全压低多为轴流式。n s小的通风机流量小,压力高,多为离心式。例如:
n s =2.7~12(15~65)* 为前弯型离心通风机
n s =3.6~16.6(20~90) 为后弯型离心通风机
n s >16.6~17.6 (90~95) 单级双进气式并联离心通风机 ns'=1.414 ns
18~36(100~200) 为轴流式通风机
n s<1.8~2.7(10~15) 罗茨式其他型式的通风机
此外比转速n s可以反映风机的几何形状并用于相似设计。
*注()内表示的数据为Mkfs制中的n s数值。
另外两个用途较广的无因次系数转速系数σ和直径系数δ:
σ=n/nm =Φ1/2/Ψ3/4 (2-28a)
δ=D2/D2M=Φ1/4/Ψ1/2 (2-28b)
式中σ,δ分别表示通风机的转速式直径与一标准风机(Φ=1,Ψ=1风机)的转速nm 式或直径D2M之比。同时从(2-28a)可以看出σ与成ns正比,即
ns=24.827σ
同时也可以用σ,δ表示Φ与Ψ:
Φ=1/σ2δ2 (2-29a)
Ψ=1/σδ3 (2-29b)
各种叶型的通风机无量纲数范围见表2-1。
图2-9表示叶轮的最佳σ~δ曲线,在曲线上的叶轮有最高的效率(但在大σ的范围,即在轴流通风机的情况不太可靠。)
表2-1 各种型式通风机的无因次数
Φ |
Ψ |
ΦΨ |
σ |
δ |
ns=24.83σ |
1.0 |
2~4 |
2~4 |
0.35~0.6 |
1.14~1.19 |
8.65~14.8 |
1 |
2~3 |
2~3 |
0.438~0.592 |
1.19~1.32 |
10.8~14.6 |
0.3 |
0.75 |
0.225 |
0.68 |
1.7 |
16.8 |
0.2 |
0.6 |
0.12 |
0.657 |
1.965 |
16.25 |
0.13 |
1.0 |
0.13 |
0.361 |
2.72 |
9.0 |
0.03 |
1.1 |
0.033 |
0.162 |
5.92 |
4.0 |
0.00185 |
1.1 |
0.00203 |
0.04 |
24.4 |
1.0 |
0.1~0.2 |
0.05~0.01 |
0.005~0.02 |
1.6~3.8 |
1.0~1.78 |
39.6~94 |
0.3 |
0.5 |
0.15 |
0.924 |
1.535 |
22.85 |
0.3 |
0.7 |
0.21 |
0.715 |
1.62 |
17.68 |
图2-9 σ~δ图上最高效率曲线
§5. 通风机的损失和效率
离心通风机中的损失可以分成流动损失,泄漏损失,轮阻损失和机械损失,其中流动损失引起通风机压力的降低,泄漏损失引起流量的减少,轮阻损失和外部机械损失则必多耗功。
一. 流动损失
整个流动损失的计算目前尚欠完善,一般以水力学中ξρC2/2的公式为基础,按单个部件分别计算各种水力效率,其中阻力系数ξ由经验数据或实验测定。
1. 进口损失
气体在进入叶轮以前要通过进气道和90度弯头,产生的压力损失为:ΔP1
ΔP1=ξ1ρC12/2 (2-33a)
式中ξ1=0.15~0.25。
2. 叶轮进口处及扩压器进口处的冲击损失
若流量与额定流量不同,则进口处相对流动的方向将不再和叶片几何方向相一致,产生叶片进口冲击损失。如图2-8和2-10的进口速度三角形可以看出:当流量从额定值增加到Qx时,Cm增至C1mx, 进口气流角从增加至β′,w′分解为一个与w1平行的分量和另一个分量Cs, Cs的大小决定了冲击损失ΔP2:
ΔP2=ρCs2 (2-33b)
式中=0.7~0.9,而Cs为:
Cs=u1(C1MX /CMX -1)= u1 (QX /Q -1)
将 u1 =u2D1/D2代入
Cs= u2D1/D2(QX /Q -1)
ΔP2=ρu22D12/D22 (QX /Q -1)2 (2-33C)
那相对值为:
ΔP2/ΔP =μ1/ψD12/D22 (QX /Q -1)2 (2-33d)
可以看出D1/D2对于冲击损失有决定性影响,为了减少冲击损失,D1/D2比值应减小,在ψ值大的时候,冲击损失将减少。
扩压器导叶与叶轮中的冲击损失一样计算:
ΔP3/ΔP=μ1/ψ(QX /Q -1)2 (2-33f)
图2-10 进口速度三角形
3. 叶轮中的损失
目前还没有准确的方法分析叶轮中的损失,因为除了摩擦损失,还有流动的分离,二次流和涡流所造成的损失。在工程上多用水力学的方法计算叶轮中的损失,为此或简化成管道流,或简化成平板流动所造成的损失,目前后一种简化更好地被掌握了。例如平板的阻力W′
在叶道中的损失分成摩擦损失ΔP'和扩压损失ΔP''两种来计算,ΔP'就用平板阻力W′的方法计算。设Ftat表示叶轮两个侧壁面和叶片两面的总面积之和,W为叶轮中的总摩擦阻力,那么:
W所做的功L为:
假设叶片为有平均叶片角的的对数螺线。那么叶片的长度:
叶片的面积
其中:Z为叶片数,,而两侧面积为,因为:
那么些
(2-33h)
式中:
叶轮中扩压损失计算利用局部损失公式:
(2-33g)
那么叶轮中损失为:
除以全压:
(2-33i)
对于高效风机,为叶轮中效率。
4).蜗壳和出口扩压器中的损失
叶轮排出的动能ρC2/2必须在叶轮后扩压件中转换为压能,这个转换称为“静压恢复”, ρC2/2转换为压力能的百分数为扩压器的恢复系数。大多数离心式通风机,这一转换在蜗壳中进行,这一过程通常引起相当大的损失,因此性能优良的通风机,总要求扩压元件的扩散度不大,气流不会分离,则压力损失表示为:
ΔP4=(0.2~0.3)ρ/2(C22-C32)
式中C3为通风机出口的速度,一般它与进口速度相差不大,为了简单可只用ρC22一项代替:
ΔP4=(0.15~0.25)ρ/2C22 (2-33j)
也可以用C2u代替C2也不产生过大误差,这是因为C2与u之间的夹角很小,即α2<20°,而且C2u/u2=Ψ/2ηhyd
ΔP4/ΔP=Ψ/4ηhyd (2-33j)
式中ηhyd为水力效率。那么总的流动损失为:
ΔPh=ΔP1 + ΔP2 + ΔP3 +ΔP4 N/m2 (2-33k)
二. 泄漏损失
叶轮与机壳件装有密封环,其间隙尽可能最小。由于这一间隙的存在,使叶轮中的总流量中,有一小部分不由蜗壳排出。这部分的流量的压力升至排出压力值之间所做的功被损失掉。
图2-11所示:气体从通风机转轴与蜗壳之间的Δ间隙处的泄漏为外泄漏。此外蜗壳靠近前盘的气流,经叶轮进口与进气口之间的间隙δ,流回叶轮进口的低压区而引起的损失,称为内泄漏损失。一般外泄漏较小。
图2-11
为了计算通过间隙的流量,必须知道间隙前后的压力差。通风机的全压力为,ΔP=Ψρ/2u22, 而在叶轮出口周围的压力,由于在蜗壳内的静压恢复而小于全压。对于后向叶片而言,叶轮出口的静压近似等于2/3ΔP。此外由于空气的旋转,压力是向中心减少。所以在密封间隙处的压力略小于2/3ΔP,再者,蜗壳的形状也有影响。这样不考虑这些影响,设间隙处的压力为2/3ΔP,计算泄漏损失会大于实际数值。对于前向叶片Ψ=2~3,但反作用度只有后向叶片的一半,这样用Ψ=1,一般可以包括大多数情况,间隙处的压力可以按ΔPst'=2/3ρ/2u22Ψ计算,通过间隙处的流量Q'为:
Q'=FμC=μπD1δi u2(2Ψ/3)1/2
一般令Ψ=1.0,μ=0.7,令δi=D2/u
Q'/Q=4μD1(2Ψ/3)1/2/(D2uФ) (2-34b)
那么理论流量Qth为:
Qth=Q+ Q' (2-34C)
为了减少内泄漏,间隙的形式可以采用图2-9所示的套口形式,而不采用对口形式。一般间隙大小为(1/100~1200)D2,当δi/ D2从0.005减小到0.0005时,可以使通风机的效率提高3~4%。
三. 叶轮的轮阻损失
叶轮的轮阻损失功率用表示,叶轮的两个侧面与机壳的前后盘间,以及轮缘的空气摩擦所造成的。
Nr=ρu23×10-3 kw (2-35)
式中β为轮阻损失计算系数。β与雷诺数Re,轮盘与壳体间的间隙B/ D2,以及圆盘外侧的粗糙度有关,根据Stodola给定β=0.81~0.89,令叶轮的总有效功Ne=ΔPQ
ΔP=ΔPthηhydth Q=Qth-Q'
Ne=ΨΦπρ/2 u22D22/4 u2
Nr/Ne=8β/1000ΨΦπ (2-36a)
对于离心式叶轮Φ=4b2 C2m u2
取β=0.88
ρ=Nr/Ne=1/(1765b2 C2m u2) (2-36b)
令ψ=1
ρ=1/262(D1/D2)3
表(2-2)给出了计算数值:
表(2-2)
D1/D2 0.15 0.3 0.5 0.7
ξ 1.135 0.142 0.0305 0.011
由此可以看到,很小的直径比,叶轮轮缘占很大比重,在D1/D2=0.15,Nr>Ne,这时叶轮轮阻损失非常重要。
四. 轴承损失:
轴承损失也是一种机械损失,不过很小。令G为转子的重量,u'为轴径的圆周速度。那么一个滑动轴承克服的摩擦所做的功为轴承损失NmeCh:
NmeCh=μG u'×10-3 (kw)
μ的数据如下:滚动轴承:0.0015,套筒和滑动轴承0.005,米切尔自位轴承0.003。
五. 通风机的功率及效率:
(一)通风机功率
1. 有效功率Ne,离心通风机的有效功率Ne为有效容积流量Q的气体通过通风机后增加的总能量为ΔPQ,即
Ne=ΔPQ (w) (2-37a)
式中ΔP为全压,单位为[w]=[N/m2][m3/s]
2.内功率Ni
实际消耗于气体上的功率为内功率Ni,
Ni=(ΔP+ΔPh)(Q+Q)+Nr (2-37b)
3.轴功率
通风机输入的功率为轴功率Ns
N= Ni + Nmeeh (2-37e)
(二)通风机效率
通风机的各种效率分别定义如下:
1.流动(水力)效率ηh
流体在通风机叶轮中要产生实际的全压升ΔP的话,必须要达到理论全压升ΔPth(无冲击进口市),扣除流动损失,才能实现:
ΔPth =ΔP+ΔP1+ΔPimp+ΔP4
那么流动效率ηh为:
ηh=ΔP/ΔP=(ΔPth -ΔPh)/ΔPth (2-38a)
ΔPh=ΔP1+ΔPimp+ΔP4
2.容积效率ηC
ηC =Q/Q th =Q/Q+Q' (2-38b)
3.机械效率ηm
ηm= Ns =Ni /Ns =(Ns- Nm )/Ns (2-38c)
Nm为机械损失功率,由轴承损失和传递损失功率所组成的。
4.轮阻效率ηr
ηr=(Ni-Nr)/ Nr
Nr=(1-ηr) Ni (2-38d)
5.通风机的内效率ηi
ηi = Ne/Ni=ΔP Q/(ΔP+ΔPth)(Q+Q')+ Nr
ηi = 1/(1/ηhηC+(1-ηC)/ηi)
故ηi =ηhηCηr (2-38e)
6.全效率:
η= Ne/Ns=ηiηm =ηhηCηrηm (2-38f)
7.静压效率ηst
离心通风机最佳工况中,出口动压是全压的10~20%左右,当偏离最佳工况到最大流量工况区间工作时,还要大一些。如果动能不利用的话,就意味着损失的加大,为此引入静压效率ηst衡量通风机的性能:
ηst=ΔPst/ΔP η (2-38g)
这样静压内效率ηsti为:
ηsti=ΔPst/ΔPηi (2-38h)
§6 离心通风机的性能曲线
图2-12
将通风机的压力ΔP,功率N和效率η等随通风机的理论流量Q不同而变化的关系绘成曲线称为通风机的性能曲线。如图2-12所示。由流动损失等在不同工况下的数值难以精确计算,则理论曲线很不准,一般都用实验方法测出,用来检验统计参数与实测参数之间的一致程度如何,另外可判定通风机的适应性,例如要求通风机的效率曲线尽可能平坦,高效率区间大一些,以适应工况的变化,使通风机在较佳的工况下工作。
一. 不计任何损失时的性能曲线
先讨论无限多叶片数,而且进口又无预旋时的通风机理论压力,由式(2-4)
ΔP th∞=ρu22 {1-QT/(πD2b2u2 ctgβ2)}
于一定的叶轮,直径D2,转速n及β2不变,则有:
ΔP th∞=A±B
因式中A,B为常数,则ΔP th∞=f(QT)是一条直线关系,如图2-4所示。图2-13是以后向叶轮为例的。再计入不限叶片数的影响。ΔP th=ΔP th∞ε(ε<0,而且近似为一常数)所以ΔP th=f(QT)的曲线仍近似为一直线,在Ⅰ曲线之下的Ⅱ曲线便是。
二. 计入损失后的性能曲线:
叶轮对气体做功的同时,一方面使气体获得能量,另外还耗功,克服各种损失。无冲击时流动损失为:
ΔP h′=ΔP 1+ΔP imp+ΔP ψ
因为C i与流量QT成正比。
ΔP h′∝QT 2
任意工况之下ΔP h※=ΔP h′ (QT ※/QT) 2,如图2-13所示。扣除此流动损失,图2-10的ΔP~f(QT)曲线就如Ⅲ所示。
冲角损失ΔP 2的大小与(QT -Q0) 2成正比;如图2-13中ΔP 2~(QT -Q0) 2的曲线。
至于效率,一般在设计流量Q0时,流动情况较好,分离,冲击损失较小,效率较高。当流量加大时,流动损失都与流量的平方成正比增加,则效率下降,而小流量时,冲击损失存在,分离损失很大,此外泄漏和轮阻损失所占相对比例增大,使效率下降,这就构成图2-12中的曲线。
图2-13
图2-14
从图2-14可以看出:
前向叶轮的离心通风机压力曲线呈驼峰状,效率曲线比后向叶轮稍低一些,功率曲线N=f(Q)一直上升,称为可过载通风机。
反向叶轮压力曲线,单调下降的。效率曲线一般高一些,当超过Q0以后的功率N几乎不增加,故存在不过载的优点。
此外,还可以得静压曲线ΔP st=f(Q)。
当转速变化时,用下列关系可以确定曲线。
Q/Q1=n/n1, ΔP/ΔP1= (n/n1)2,
N/N1=(n/n1)3,
以上曲线为有因次参数的曲线。一般可以用于新产品的初步性能试验的使用。
三. 无因次性能曲线
通风机的主要无因次参数φ,ψ,λ,η之间的关系也可以用曲线表示。由于同一系列的通风机,其相应点的无因次参数都相同,故无因次参数性能曲线代表一系列通风机的性能,用起来很简便。这些曲线可以由试验直接得出,可以从有因次曲线计算转换而来。
四. 通风机的空气动力学略图
在通风机系列产品中,由于用相似设计方法设计的;它们的几何形状相似。无因次性能曲线也相同。通风机流道的各部分几何尺寸,可以用=100为基准。画出其相对于D2的尺寸,这样通流部件图再加上无因次性能曲线就组成了通风机的空气动力学简图式空气动力学略图。
五. 系列产品的综合性能曲线
把一个系列产品的性能参数,全压,流量,转速,蜗壳出口速度,效率和叶轮直径等关系综合在一起的曲线称为综合特性曲线,见图2-15的蜗壳出口速度Ca为横坐标,以ΔP,Q为纵坐标的曲线,它们的关系可简述如下:
设Fa为蜗壳出口面积,每个机ξ(对应一个叶轮直径)的流量为Q=Ca Fa,那么流量Q与Ca成正比,不同机ξ的流量不同,机ξ大的(直径大),蜗壳大,Q也大。
ΔP与Ca的关系按相似换算为:
ΔP=1/2ψρu22
因为Ca=φπD22 u2/4 Fa
所以Ca∝u2,且ΔP∝Ca2。
图2-15
六. 系列产品的对数坐标曲线
对数坐标曲线是把在标准进口状态下对应于无因次性能曲线上的一个工况点的所有同系列通风机的叶轮直径,转速n,圆周速度u2以及相应的流量Q,全压ΔP,功率N全部表示出来,即对数坐标曲线表示同系列通风机主要参数之间的关系。这种曲线可供使用部门方便地选择所需的通风机。制造厂可以用这种曲线合理地确定该型号系列产品的性能型谱,即确定机ξ和转速,为产品的标准化,通用化,系列化创造一系列条件。
图2-16
七. 对数坐标曲线的原理:
对数曲线以Q为横坐标,ΔP为纵坐标,u2也为纵坐标,那么等Q线,等ΔP 线,等u2线都平行于坐标轴,如图2-16所示。
1.等转速线,等直径线,等功率线的斜率有特点;
等速度线与流量坐标Q的夹角为33.7°;等直径线与流量坐标轴的夹角为63.5°;功率为常数也为一直线方程,斜率为α3,即等功率线与Q轴夹角为135°,等功率线是近似的,因为多条性能曲线上的全压效率不一定相等。
2. ξ参数坐标长度的比例关系:
根据上面的方程,可以得出以流量坐标mQ为基准,压力坐标mp=mQ。当ΔP=常数时,转速坐标mn=2mQ (方向相反), mD=2mQ,功率坐标mN=mQ,圆周速度坐标mu=2mQ。关于对数坐标曲线见图2-17。
图2-17 (a) 图2-17 (b)