热水供热系统的水锤计算.doc

《热水供热系统的水锤计算.doc》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、热水供热系统的水锤计算清华大学热能系蔡启林北京建筑工程学院李锐摘要　热水供热系统发生水锤是集中供热发展进程中出现的实际问题，一旦发生，破坏力大，会导致财、物的重大损失，甚至危及人身安全．本文叙述了热水供热系统水锤的概念，建立了基本方程组，处理了常见的各种复杂边界条件，考虑了出现蒸汽穴的影响．所编制的软件对实例的运行结果表明：物理上看合理、正确，计算结果与某国咨询公司计算结果吻合较好．一、水锤的概念　　　　当压力管道中的流体因某些原因而产生流速的急剧变化时，由于流体的惯性作用而引起管道的流体压力急剧变化，这种现象作为水锤现象或流体瞬变过程．　　水锤是管道瞬变流动中的一种压力波．下面以图1所示的简

2、单的供热系统模型为例来说明水锤波的形成机理及其传播、反射和叠加过程．为便于分析，假设循环水泵a产生的压头与流量无关，水泵的吸入端和压出端连接膨胀水箱e和f，使管网在任何水力工况下，水泵吸入端6和压出端1的总压头都保持固定不变。　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　图1　热水供热网模型图　　　　　　　　　　　　　a--循环水泵；b--止回阀；c--干线供水管上的调节阀；　　　　　　　　　　　　　　　　d--热用户；e--膨胀水箱；f--膨胀水箱　　　　当阀C全开、供热系统处于稳态时．管中流速为VQ，倘若由于某种原因，阀C突然瞬间关闭，则点3流速突然滞止为零，点2出现

3、突然压力升高，从而使此部分流体的密度增大、流体压缩和管道膨胀；点3出现突然的压力降低，使这部分流体的密度减小、流体膨胀和管壁压缩．阀门C前后的压力以弹性波的形式由阀门C迅速传向全部热水网管道。　　当快速关闭阀门C时，上游管道中的后续流体仍以原来流速+V0继续向阀门C方向运动，与首先停留的流体相毗邻的流体一起填充了前一部分流体压缩和管壁膨胀而扩大了的容积之后，其流速也相继从+V0骤降为零，随即产生升压，流体压缩和管壁膨胀．这样依次由阀门C向循环水泵a延伸．当这种升压运动传播到假设压力恒定点1时，1一2管道中的全部流体速度都变为零，在这一瞬间，管道中流体的动能全部转化为流体的弹性压缩能和管壁的变形

4、能，1-2管在正压力作用下呈现膨胀状态。　　当水锤波达到点1时，由于点1压力不变，所以在压力波达到的瞬间，点1处于一种不平衡状态。点1处流体密度首先恢复到原来状态，升压随即消失。由这一瞬间开始，压力波向反方向运动，管段1-2中流体产生与原来相反的流动，这种降压以同样的传播速度向阀门C处传播，当其返回阀门C处时，管道1-2中的压力全部为H1，即与点1的压力相等。　　由于管道1-2中流体向循环水泵a流动，而阀门处于完全关阀状态，无流体补充以维持这种流动，从而引起该处压力下降，流体密度减小和管壁的收缩，这种降压运动又以同样传播速度向a处传播。当到达点1时，管道中的流速都变为零。在这一瞬间，管道中全部

5、流体的动能变为流体的弹性膨胀能和管壁的收缩变形能，管道中压力降低。管道1-2呈收缩状态。　　虽然管道1-2中流体压力降低，但由于点1的水头高于管道1-2中的压力，于是使点1处的流体产生向阀门C的流动，管道1-2中的压力开始升高，这种压力波以同要的速度向阀门C处传播，当传播到阀门C处后，压力回升到稳定值H1。管道1-2呈收缩状态。　　由于磨擦以及流体和管壁的非完全弹性等作用，使振荡发生阻尼，压力波逐渐衰减，从而最终使流动消失。　　在阀门C的另一例，点3将出现类似情景。但开始时是减压波向阀门C下游管道传播，液体汽化现象称为气穴或空穴。由于气穴处是气体而不是液体．而破坏水流的连续性．发生这种边界条件

6、变化，通常会使上游的液柱加速而下游液柱减速，从而使上游的液柱赶上下游的波柱，导致气泡破裂、溃灭，或者当气穴中的气体被带到下游压强较大的区域时，发主气穴的条件消失，也发生空泡的溃灭．在空泡溃灭的过程中，空泡溃灭水锤极大，可以大到足以破坏管子。　　在管道中发生水锤现象后，任何一点的流速和压力不仅与该点的位置有关，而且与时间有关，因而是非稳定流．在水锤现象发生的过程中，由于液体受到巨大压强的作用，在计算中必须考虑液体的压缩性和管壁的弹性变形，否则，将会导致与实际不符合的结果。　　二、基本方程　　　　热水供热网的水力瞬变可以用下面的一些对拟线性双曲编微方程描述：　　　　　　　　　　　　　　　　（1）　

7、　连续方程：　　　　　　　　　　　　　　　　　　　　　　　　（2）　　式中：H--测压管水头（mH2O）;t--时间，（S）;　　　　　a--水锤波传播速度，（m/s）；g--重力加速度，（m/s2）;　　　　　A--管道横截面积，（m2）；Q--管道流量，（m3/s）　　　　　x---计算点距起始点的坐标蹁，（m）f-管道沿程磨擦阻力系数；　　　　　D--管内径，（m）　　应用特征线解法可以将方