陳嘉禾

（湖南大學(xué) 土木工程學(xué)院，湖南 長(zhǎng)沙 410082）

設(shè)計(jì)污水管道就必須對(duì)污水管道進(jìn)行水力計(jì)算，對(duì)此傳統(tǒng)的方法是根據(jù)各水力參數(shù)的限制范圍，在水力計(jì)算圖中查出一組可行的設(shè)計(jì)方案。這實(shí)際上是一項(xiàng)十分繁瑣的工作，并且精度不高，不利于設(shè)計(jì)方案的優(yōu)化。筆者在本文中將介紹一種新的方法——用MATLAB 來編程。

MATLAB 是目前常用的工程軟件，它集數(shù)值計(jì)算、矩陣計(jì)算和圖形繪制等功能于一體，廣泛應(yīng)用于各種數(shù)學(xué)及工程計(jì)算中。它采用數(shù)學(xué)語言編寫程序，更加通俗易懂，并且自帶矩陣和數(shù)組的運(yùn)算函數(shù)，減少了編程和調(diào)試的工作量，這就是MATLAB 相比于其他語言如Fortran、C++等的顯著優(yōu)點(diǎn)。

截至目前，已經(jīng)有不少學(xué)者成功應(yīng)用MATLAB 解決了一些工程計(jì)算問題，如進(jìn)行梯形明渠、供熱管網(wǎng)、燃?xì)夤芫W(wǎng)的水力計(jì)算；建立供水管網(wǎng)平差模型等。大大縮短了計(jì)算時(shí)間，提高了計(jì)算精度和工作效率。

1 污水管道水力計(jì)算的試算法及程序框圖

1.1 污水管道水力計(jì)算的試算法

對(duì)于無壓圓管流，其水力計(jì)算基本公式有：

對(duì)以上公式進(jìn)行整理可知，h/D和R均為θ 的函數(shù)，而v、I和θ 均同時(shí)為θ 和v的函數(shù)。因此，可先假定設(shè)計(jì)流速v的值，再算出其他參數(shù)的值。

此外，還應(yīng)考慮到在污水管道水力計(jì)算過程中可能涉及到的約束條件：

a. 管徑D對(duì)應(yīng)的最小流量Q（min）；

b. 管徑D對(duì)應(yīng)的最大充滿度h/D（max）；

c. 管徑D對(duì)應(yīng)的最大流速v（max）；

d. 管徑D對(duì)應(yīng)的最小水力坡度I（min）；

e. 管徑D對(duì)應(yīng)的最大水力坡度I（max）。

結(jié)合上述分析及約束條件，可整理得出污水管道水力計(jì)算的試算法，步驟如下：

（1）繪出計(jì)算草圖，對(duì)各點(diǎn)、管段進(jìn)行編號(hào)；

（2）計(jì)算各管段設(shè)計(jì)流量；

（3）根據(jù)約束條件a 和設(shè)計(jì)流量，求出各管段對(duì)應(yīng)的最大管徑；

（4）對(duì)每段管段的每種可能管徑，使設(shè)計(jì)流速v從最小值（0.6 m/s）以步長(zhǎng)（0.01 m/s）增大到相應(yīng)的最大值（約束條件c），計(jì)算相應(yīng)的θ、R、h/D、I；

（5）判斷求出的h/D和I是否同時(shí)滿足約束條件b、d、e；

（6）若滿足，則記錄這一組設(shè)計(jì)方案；若不滿足，則繼續(xù)試算下一個(gè)設(shè)計(jì)流速v，直至計(jì)算出所有管段的所有可行方案為止。

1.2 程序框圖

試算法計(jì)算程序框圖如圖1 所示。

圖1 試算法計(jì)算程序框圖

2 計(jì)算實(shí)例

2.1 例題

【例】用試算法計(jì)算排水管道，如圖2 所示。設(shè)節(jié)點(diǎn)流 量Q1=24.28，Q2=12.4，Q3=14.6，Q4=9.5，Q5=8.7，Q6=12.8，單位均為L(zhǎng)/s。

圖2 具有6個(gè)排出口的污水管道示意圖

2.2 示例程序

用MATLAB 編寫如下程序就能對(duì)這一管道進(jìn)行水力計(jì)算。

%水力約束矩陣,按列依次為D-Q(min)對(duì)應(yīng)關(guān)系,D-(h/D)(max)對(duì)應(yīng)關(guān)系,D-v(max)對(duì)應(yīng)關(guān)系,D-I(min)對(duì)應(yīng)關(guān)系,D-I(max)對(duì)應(yīng)關(guān)系

%標(biāo)準(zhǔn)管徑向量

注:在實(shí)際應(yīng)用中,技術(shù)人員應(yīng)根據(jù)工程實(shí)際情況適當(dāng)改變水力約束矩陣中的數(shù)值。(如街坊管道的最小管徑為200 mm,市政管道的最小管徑為300 mm;坡度閾值可參考實(shí)際地形坡度設(shè)定等)

2.3 計(jì)算結(jié)果

由于計(jì)算結(jié)果數(shù)目較大, 此處只列出第一段管段第一種管徑(200 mm)的計(jì)算結(jié)果:

設(shè)計(jì)流量計(jì)算結(jié)果(L/s):

管徑流速充滿度水力坡度

第1 段管段：當(dāng)管徑為200 mm 時(shí)共有118 種設(shè)計(jì)方案。

3 分析討論

3.1 與手工查表的比較

手工查表需要先在不同管徑對(duì)應(yīng)的水力計(jì)算表中找到流量所對(duì)應(yīng)的列，然后在考慮多種約束條件的情況下選擇一組合適的水力參數(shù)值。這一過程步驟繁瑣、用時(shí)較長(zhǎng)，加之水力計(jì)算表的數(shù)據(jù)密集性較高，給技術(shù)人員準(zhǔn)確查找數(shù)據(jù)帶來了不小的困擾。此外，第一次查出的數(shù)據(jù)可能無法滿足后期設(shè)計(jì)埋深的要求，導(dǎo)致技術(shù)人員需要重新查找，這也將耗費(fèi)大量時(shí)間?？傮w來說，這種方法過程繁瑣、效率較低。

在筆者的計(jì)算機(jī)上運(yùn)行整個(gè)程序，總用時(shí)為749 s（在不同的計(jì)算機(jī)上運(yùn)行結(jié)果不同，用時(shí)與計(jì)算機(jī)的速度有關(guān)），這比手工查表節(jié)省了許多時(shí)間。因此，與手工查表相比應(yīng)用MATLAB 進(jìn)行污水管道水力計(jì)算，速度快，效率高。

3.2 對(duì)設(shè)計(jì)方案優(yōu)化的意義

對(duì)程序計(jì)算結(jié)果進(jìn)行整理，見表1。

表1 各管段設(shè)計(jì)結(jié)果方案數(shù)統(tǒng)計(jì)

由表可知，MATLAB 能在保證準(zhǔn)確度的情況下給出多種設(shè)計(jì)方案，這不僅能避免繁冗的手工查找工作，而且有利于技術(shù)人員根據(jù)工程項(xiàng)目的特殊限制要求對(duì)計(jì)算結(jié)果加以篩選和比較，從而快速確定最佳設(shè)計(jì)方案。

4 結(jié)論

使用MATLAB 進(jìn)行水力計(jì)算，速度快、效率高，有利于設(shè)計(jì)方案的優(yōu)化，是給排水科學(xué)與工程從業(yè)者應(yīng)用計(jì)算機(jī)進(jìn)行污水管道水力計(jì)算的理想工具。