输水管数学模型.docVIP

输水管线的数学模型孟玲玲 俞剑明林云龙 摘要：本问题为一个水力学问题，是一个长输水管道输水的问题，先将长输水管道的数据进行简化，使复杂的管道数据更加简单直观化，便于分析。 首先管道建成后会进行小流量灌水，将小流量灌水问题分成了跌水过程，小流量上升过程，和拐点损失过程这三个过程。对于跌水过程，利用非满流的巴甫洛夫斯基公式计算水头损失，再通过能量守恒公式进行小间隔的递推分析从而解出速度方程，对于小流量上升过程，可看成连通器问题进行分析。通过以上局部求解的联系，得到当初始流量为0.6时，通完全管的时间为11.17天。而对于运动曲线的模拟，利用水利学连续性方程，动量方程和水管的边界条件，简化成等宽明渠的问题进行分析。 其次对于通过两次通水使通水量达到所需通水量，我们在设计通水量时就要考虑到每个管子拐点处在增加通水量时所增加的压力值平均，这样才能使增加水量时压力的增加量的最大值最小。其次我们就要计算拐点速度时间和压力分布，这里对于上升过程，由于水量的增加就不能看成连通器问题了，我们将其看成满贯上升问题，利用海曾－威廉公式和能量守恒公式进行分析。所以第一次用54分钟开启水库处入口阀门，时流量增至2.9，再用37分钟缓慢开启水库处入口阀门，是流量达到5.208。最后利用速度分布图得到水压分布。 关键字： 小流量 跌水过程 上升过程 拐点损失过程 正文： ·1 问题重述 某输水管线是利用114米地面落差有压、重力流输水工程。管径2.2米，管线长度176公里，地面高程变化很大。 当输水管线建成后，首次通水时，为了防止水流速波动产生水击破坏管线，只能以每秒0.6立方米的流速由某水库向管线内灌水。 当被供水城市水厂出流稳定后，再逐步提高入流速度达到每天45万吨（一期工程单线）的额定工作状态。 利用data1.txt 中提供的输水管线信息数据来解决下面的运动界面追踪问题： 1．第一阶段：假设以每秒Q（=0.6）立方米入流速度向管线内灌水，模拟出水头前沿面的运动过程，计算出水头前沿面顺序到达管线高程谷点和峰点的时间。 2．第二阶段：再分两阶段提高入流速度达到每天45万吨的额定工作状态，设计出通水方案，计算出由一种恒定流达到另一种稳定流后的管线水头压力分布和需要的时间。 ·2 问题分析 在第一阶段即首次通水时，水库以每秒0.6立方米的流速向管线内灌水，则对于管道而言，水流由峰点流向谷点时为非满流状态，由谷点流向峰点时为满流状态。利用水力学理论知识，利用巴普罗夫斯基公式和能量守恒可以得到不用点的速度的关系，而由于水头损失是一个和速度有关的量，利用微积分的思想，对于小间隔管长进行分析，求出小段前和末端的速度，这样就能得到离散的速度值，这就相当于在求解一个微分方程。对于上升过程就转化成一个简单的联通器进行考虑。 在第二阶段由于流量变大对于上升过程就不能当成是连通器，由于水流速过大时水流流入谷底时，基本已经把拐点处的空间已经填满，上升过程就能够当成是满流的上升，这样利用海曾-威廉公式和能量守恒就能类比上述方法得到求解。这样就能模拟出水流的运动过程了。 ·3 问题假设 为便于分析解决该问题，我们作出以下合理假设： 在第一阶段，水流由水库进入接管井时为满流状态，即水刚进入接管井时充满管的截面，由此可以根据流量Q就算出初速度。 在模型建立时，我们暂不考虑局部水头损失，水头损失以沿程水头损失算。 假设第一阶段为小流量灌水，上升过程相当于向连通器冲水。 假设第二阶段水流足够大，上升过程为一个满流流动的过程。 ·4 符号说明 ：水头损失 ：水断面面积 ：湿周 i:水利坡降 ：巴甫洛夫斯基公式粗糙系数 Q：流量 g: 重力加速度 D：管直径 n:摩擦系数 ·5 模型建立与求解 问题一为一个小流量灌水的问题，对于小流量灌水，我们不仅要通过水头损失，波努力方程等计算满流非满流的时管道各点的速度，时间，由于有些前面的一些峰值点低于后面的峰值点，前面的要被淹没。淹没的情况又有完全、非完全淹没的可能。所以我们分为以下几个阶段来考虑这个问题。 一：跌水过程 水头前沿面达到谷点的时间计算问题，可以根据能量守恒近似地计算出实际水体流动速度，跌水路程除以速度，得到水体到达谷点的时间。 由于水下流时是一个非满管的过程，我们利用巴甫洛夫斯基公式进行计算水头损失： 式中： 巴甫洛夫斯基公式粗糙系数，取n=0.0109 带入 得：