水錘計算方法

水錘計算方法PAGEPAGE82——第一節(jié)概述一、水電站的不穩(wěn)定工況機組在穩(wěn)定運行時，水輪機的出力與負荷相互平衡，這時機組轉(zhuǎn)速不變，水電站有壓引水系統(tǒng)(壓力隧洞、壓力管道、蝸殼及尾水管)中水流處于恒定流狀態(tài)。在實際運行過程中，電力系統(tǒng)的負荷有時會發(fā)生突然變化(如因事故突然丟棄負荷，或在較短的時間內(nèi)啟動機組或增加負荷)，破壞了水輪機與發(fā)電機負荷之間的平衡，機組轉(zhuǎn)速就會發(fā)生變化。此時水電站的自動調(diào)速器迅速調(diào)節(jié)導葉開度，改變水輪機的引用流量，使水輪機的出力與發(fā)電機負荷達到新的平衡，機組轉(zhuǎn)速恢復到原來的額定轉(zhuǎn)速。由于負荷的變化而引起導水葉開度、水輪機流量、水電站水頭、機組轉(zhuǎn)速的變化，稱為水電站的不穩(wěn)定工況。其主要表現(xiàn)為：(1)引起機組轉(zhuǎn)速的較大變化由于發(fā)電機負荷的變化是瞬時發(fā)生的，而導葉的啟閉需要一定時間，水輪機出力不能及時地發(fā)生相應變化，因而破壞了水輪機出力和發(fā)電機負荷之間的平衡，導致了機組轉(zhuǎn)速的變化。丟棄負荷時，水輪機在導葉關(guān)閉過程中產(chǎn)生的剩余能量將轉(zhuǎn)化為機組轉(zhuǎn)動部分的動能，從而使機組轉(zhuǎn)速升高。反之增加負荷時機組轉(zhuǎn)速降低。(2)在有壓引水管道中發(fā)生“水錘”現(xiàn)象當水輪機流量發(fā)生變化時，管道中的流量和流速也要發(fā)生急劇變化，由于水流慣性的影響，流速的突然變化使壓力水管、蝸殼及尾水管中的壓力隨之變化，即產(chǎn)生水錘。導葉關(guān)閉時，在壓力管道和蝸殼中將引起壓力上升，尾水管中則造成壓力下降。反之導葉開啟時，在壓力管道和蝸殼內(nèi)引起壓力下降，而在尾水管中引起壓力上升。(3)在無壓引水系統(tǒng)(渠道、壓力前池)中產(chǎn)生水位波動現(xiàn)象。無壓引水系統(tǒng)中產(chǎn)生的水位波動計算在第八章已介紹。二、調(diào)節(jié)保證計算的任務水錘壓力和機組轉(zhuǎn)速變化的計算，一般稱為H0+△H，而左端水庫保持不變?yōu)镠，因此“邊界”處的水體不能保持平衡，管道中的水體在△H壓差作用下將逆流向水庫。在t=L/a后的dt1時段內(nèi)，首先是緊靠水庫dxn管段內(nèi)發(fā)生變化，流速將由0變?yōu)?V0，壓力由H0+△H變?yōu)镠0；管壁及水體隨著水錘壓力的消失恢復至原狀。同理接再經(jīng)dt2、dt3、…時段，在相應dxn-1、dxn-2、dxn-3…管段中將發(fā)生同樣的變化，如圖9-1(b)。直到t=2L/a時刻，整個管道中的壓力、流速、管徑及水的密度均恢復到初始狀態(tài)。這說明，水錘波在水庫處要發(fā)生反射，反射特點是“等值異號”反射，即反向波與入射波的數(shù)值相同，均為，但符號相反，升壓波反射為降壓波。(3)t=2L/a~3L/a。當t=2L/a時，水錘波傳播到閥門處A點，由于閥門已關(guān)閉，加之水流的慣性作用，管道中的水繼續(xù)流向水庫。在t=2L/a~3L/a時段內(nèi)，首先是緊近閥門dx1管段內(nèi)發(fā)生變化，依次傳到dx2、dx3…管段，到3L／a時刻，流速將由-V0變?yōu)?，壓力由H0變?yōu)镠0－△H，管徑為D－△D，水的密度變?yōu)棣眩鳓选．旈y門全關(guān)閉時，水錘波在閥門處的反射特點是“等值同號”反射，即反向波與入射波的數(shù)值和符號不變，從水庫傳來降壓波仍反射為降壓波。(4)t=3L/a~4L/a。當t=3L/a時，水錘波又回到水庫處D點，由于管道壓力比水庫低△H，則D點壓力不能維持平衡，因此水庫的水又向閥方向流動，這時水庫將閥門反射回來的降壓波又反射為升壓波，到t=4L/a時，管道流速將由0變?yōu)閂0，壓力由H0－△H0變?yōu)镠0T=4L/a稱為水錘波的“周期”。每經(jīng)一個周期，水錘現(xiàn)象就重復一次上述過程。水錘波在管中傳播一個來回的時間tr=2L/a，稱之為“相”，兩個相為一個周期T=2tr。閥門突然開啟時，水錘現(xiàn)象與上述情況相反。如果不存在水力摩阻，則上述的水錘過程將無休止地反復下去，但由于水力摩阻的存在，水錘過程不可能無休止地振蕩下去，壓力波因摩擦損失而逐漸衰減，在一定時段內(nèi)逐漸消失。綜上所述，我們可以初步得出以下幾點結(jié)論：(1)水錘壓力實際上是由于水流速度變化而產(chǎn)生的慣性力。當突然啟閉閥門時，由于啟閉時間短、流量變化快，因而水錘壓力往往較大，而且整個變化過程是較快的。(2)由于管壁具有彈性和水體的壓縮性，水錘壓力將以彈性波的形式沿管道傳播。(3)水錘波同其它彈性波一樣，在波的傳播過程中，外部條件發(fā)生變化處(即邊界處)均要發(fā)生波的反射。其反射特性(指反射波的數(shù)值及方向)決定于邊界處的物理特性。二、水錘波的傳播速度在水錘過程的分析與計算中，波速是一個重要的參數(shù)。它的大小與管壁材料、厚度、管徑、管道的支承方式以及水的彈性模量等有關(guān)。由水流的連續(xù)方程并考慮水體和管壁的彈性后，可導出水錘波的傳播速度為(9-1)式中K——水的體積彈性模量，一般為2.06×103MPa；E——管壁材料的縱向彈性模量(鋼村E=2.06×105MPa，鑄鐵E=0.98×105MPa，混凝土E=2.06×104MPa)；g——重力加速度；D——管道內(nèi)徑；δ——管壁厚度。為聲波在水中的傳播速度，隨水溫度和壓力的升高而加大，一般可取為1435m/s。在缺乏資料的情況下，露天鋼管的水錘波速可近似地取為1000m/s，埋藏式鋼管可近似取為1200m/s，鋼筋混凝土管可取900m/s~1200m/s。第三節(jié)水錘基本方程及邊界條件為求解水錘壓力升高問題，需要建立基本方程。基本方程與相應的邊界條件聯(lián)立，用解析方法或數(shù)值計算方法求解水錘值及其變化過程。一、水錘基本方程(一)、基本方程對有壓管道而言，不論在何種情況下都應滿足水流的運動方程及連續(xù)方程。當水管材料、厚度及直徑沿管長不變時，其運動方程為：(9-2)將管道材料及水體當作彈性體考慮，其連續(xù)方程為：(9-3)式中H——壓力水頭；V——管道中的流速，向下游為正；a——水錘波傳播速度；f——水流摩擦阻力系數(shù)；D——管道直徑；x——距離，其正方向與流速取為一致；t——時間。上面二式中，因流速V與波速a相比數(shù)量較小，故可忽略和項。另外，為了簡化計算，使方程線性化，忽略摩擦阻力的影響。當x軸改為取閥門端為原點，向上游為正時，如圖9-2，方程(9-2)、(9-3)可簡化為：(9-4)(9-5)式(9-4)和式(9-5)為一組雙曲線型偏微分方程，其通解為：(9-6)(9-7)圖9-2水擊計算示意圖式中H0和V0為初始水頭和流速；F和f分別為兩個波函數(shù)，其量綱與水頭H相同，故可視為壓力波。F(t-x/a)表示以波速a沿x軸負方向傳播的壓力波，即逆水流方向移動的壓力波，稱為逆流波；f(t+x/a)表示以波速a沿圖9-2水擊計算示意圖任何斷面任何時刻的水錘壓力值等于兩個方向相反的壓力波之和，而流速值為兩個壓力波之差再乘以－g/a。如果知道了t時刻在x位置處的水錘波函數(shù)F(t-x/a)，則當時間變?yōu)閠1=t+Δt，研究x1=x+aΔt處的逆流波函數(shù)=，其值不變，證明了F(t-x/a)沿逆水流方向的傳播特性。反之研究t1=t+Δt時刻在位置x1=x-aΔt處的順流波函數(shù)，可以證明f(t+x/a)沿順水流方向的傳播特性。(二)、水錘計算的連鎖方程若已知斷面A(見圖9-2)在時刻t的壓力為，流速為，由(9-6)和(9-7)消去f后，得：同理可寫出時刻后B點的壓力和流速的關(guān)系：由于，由上述二式得(9-8)同理：(9-9)方程(9-8)和(9-9)為水錘連鎖方程。連鎖方程給出了水錘波在一段時間內(nèi)通過兩個斷面的壓力和流速的關(guān)系。但前提應滿足水管的材料、管壁厚度、直徑沿管長不變。水錘連鎖方程(9-8)和(9-9)用相對值來表示為：(9-10)(9-11)式中稱為管道特性系數(shù)；，稱為水錘壓力相對值；為管道相對流速。二、水錘的邊界條件應用水錘基本方程計算壓力管道中水錘時，首先要確定其起始條件和邊界條件。(一)、初始條件當管道中水流由恒定流變?yōu)榉呛愣鲿r，把恒定流的終了時刻看作為非恒定流的開始時刻。即當t=0時，管道中任何斷面的流速V=V0；如不計水頭損失，水頭H=H0。(二)、邊界條件1．管道進口。管道進口處一般指水庫或壓力前池。水庫水位變化比較慢，在水錘計算中不計風浪的影響，認為水庫水位為不變的常數(shù)是足夠精確的。壓力前池的水位變化情況與渠道的調(diào)節(jié)類型有關(guān)。自動調(diào)節(jié)渠道的前池水位變化雖大，但與管道中水錘計算時間相比，變化還是緩慢的。非自動調(diào)節(jié)渠道，水位變化較小，一般只有幾米，在水錘計算中也認為前池水位不變。所以管道進口邊界條件為：Hp=H02．分岔管。分岔管的水頭應該相同，即Hp1=Hp2=Hp3=…=Hp分岔處的流量應符合連續(xù)條件，即ΣQ=03．分岔管的封閉端。在不穩(wěn)定流的過程中，當某一機組的導葉全部關(guān)閉，或某一機組尚未裝機，而岔管端部用悶頭封死，其邊界條件為：Qp=04．調(diào)壓室。把調(diào)壓室作為斷面較大的分岔管，其邊界條件為:調(diào)壓室內(nèi)有自由水面，而隧洞、調(diào)壓室與壓力管道的交點和分岔管相同。5、水輪機。水電站壓力管道出口邊界為水輪機，水輪機分沖擊式和反擊式，兩種型式的水輪機對水錘的影響不同。(1)沖擊式水輪機。沖擊式水輪機的噴嘴是一個帶針閥的孔口。水輪機轉(zhuǎn)速變化對孔口出流沒有影響，對沖擊式水輪機，噴嘴全開時斷面積為ωmax，流量系數(shù)為φ0，根據(jù)《水力學》的孔口出流規(guī)律，過流量為：當孔口關(guān)至ωi時一般假定：φ0=φ，均為流量系數(shù)，所以式中：，稱為相對開度，為任意時刻水錘壓力相對值。而,所以(9-12)這是沖擊式水輪機噴嘴的出流規(guī)律，也即閥門處A點的邊界條件。(2)反擊式水輪機。反擊式水輪機有如下特點：(i)反擊式水輪機有蝸殼、尾水管及導水葉，過流特性與孔口出流不完全相同。(ii)反擊式水輪機的轉(zhuǎn)速與水輪機的流量互相影響。(iii)流量突然改變時，不僅在壓力管道中，而且在蝸殼、尾水管中也發(fā)生水錘。尾水管中發(fā)生的水錘現(xiàn)象與蝸殼相反，即導水葉關(guān)閉時發(fā)生負水錘，開啟時發(fā)生正水錘。蝸殼、尾水管中的水錘影響水輪機的流量，繼而又對水錘產(chǎn)生影響。由此可見，反擊式水輪機的過水能力與水頭H、導葉開度a和轉(zhuǎn)速n有關(guān)。即Q=Q(H,a,n)，需要綜合運用管道水錘計算方程、水輪機運轉(zhuǎn)特性曲線、水輪機組轉(zhuǎn)速方程等進行求解，因此增加了問題的復雜性。為了簡化計算，常假定壓力管道出口邊界條件為沖擊式水輪機，然后再加以修正。圖9-3水輪機開度變化規(guī)律(三)開度按直線規(guī)律變化圖9-3水輪機開度變化規(guī)律水輪機導葉和閥門的關(guān)閉規(guī)律與調(diào)速系統(tǒng)的特性有關(guān)，實際的關(guān)閉規(guī)律如圖9-3所示。從全開(τ0=1.0)到全關(guān)(τ=0)的全部歷時為TZ，曲線開始一段接近水平，關(guān)閉的速度極慢，這是由于調(diào)節(jié)機構(gòu)的慣性所決定的，在這段過程中，引起的水錘壓力很小，對水錘計算沒有多大實際意義。在接近關(guān)閉終了時，閥門的關(guān)閉速度又逐漸減慢，曲線向后延伸，這種現(xiàn)象只對閥門關(guān)閉接近終了時的水錘壓力有影響。因此為了簡化計算，常取閥門關(guān)閉過程的直線段加以適當延長，得到Ts，Ts稱為有效關(guān)閉時間。在缺乏資料的情況下，可近似取Ts=(0.6~0.95)TZ。直線規(guī)律關(guān)閉時，一個相長，一個相的開度變化，負號表示閥門關(guān)閉；正號表示閥門開啟。第四節(jié)簡單管水錘的解析計算簡單管是指壓力管道的管徑、管壁材料和厚度沿管長不變。解析法的要點是采用數(shù)學解析的方法，引入一些符合實際的假定，直接建立最大水錘壓力的計算公式。簡單易行，物理概念清楚，可直接得出結(jié)果。直接水錘和間接水錘水錘有兩種類型：直接水錘和間接水錘。(一)直接水錘水錘波在管道中傳播一個來回的時間為2L/a，稱為“相”。當水輪機開度的調(diào)節(jié)時間TS≤2L/由于水管末端未受水庫反射波的影響，因此基本方程(9-6)和(9-7)中的波函數(shù)f(t+x/c)=0，然后從二式中消去F(t-x/c)得直接水錘公式(9-13)公式(9-13)只適用于TS≤2L/a(1)當閥門關(guān)閉時，管內(nèi)流速減小，V-V0<0為負值，△H為正，產(chǎn)生正水錘；反之當開啟閥門時，即V-V0>0，△H為負，產(chǎn)生負水錘。(2)直接水錘壓力值的大小只與流速變化(V-V0)的絕對值和水管的水錘波速a有關(guān)，而與開度變化的速度、變化規(guī)律和水管長度無關(guān)。當管道中起始流速V0＝4m/s，a＝1000m/s，終了流速V＝0時，壓力升高值為：=-1000(0-4)/9.81=407.7m，因此在水電站中應當避免直接水錘。(二)間接水錘若水輪機開度的調(diào)節(jié)時間TS＞2L/a，當閥門關(guān)閉過程結(jié)束前，水庫異號反射回來的降壓波已經(jīng)到達閥門處，因此水管末端的水錘壓力是由向上游傳播的水錘波F和反射回來的水錘波f發(fā)生間接水錘時，水錘壓力波的消減、增加過程是十分復雜的。間接水錘是水電站中經(jīng)常發(fā)生的水錘現(xiàn)象，也是要研究的主要對象。工程中最關(guān)心的是最大水錘壓力。由于水錘壓力產(chǎn)生于閥門處，從上游反射回來的降壓波也是最后才達到閥門，因此最大水錘壓力總是發(fā)生在緊鄰閥門的斷面上。下面應用前面的水錘連鎖方程(9-10)和(9-11)及管道邊界條件，推求閥門處各相水錘壓力的計算公式。二、計算水管末端各相水錘壓力的公式(一)第一相末的水錘壓力設閥門為A點，水庫為B點，水錘波從A到B點的連鎖方程為：邊界和初始條件：t=0時，；在水庫進水口B點，所以：將A點邊界條件代入上式：(2)水錘波從B到A的連鎖方程：因，和A點邊界條件，上式變?yōu)橐?t=tr=2L/a(9-14)(二)第二相末的水錘壓力(l)寫出水錘波從A→B的連鎖方程式：由B點的邊界條件得＝0，上式可改寫成：所以：(2)寫出水錘波從B→A的連鎖方程式：把，，和代人上式，并用2代替4t表示第二相末，得：(9-15)(三)第n相末的水錘壓力用同樣原理可以得出以后任意n相末的水錘壓力計算公式，其一般公式為：(9-16)利用式(9-14)~(9-16)，可以依次求出各相末閥門處的水錘壓力，得出水錘壓力隨時間的變化關(guān)系。上面是閥門關(guān)閉情況，當閥門或?qū)~開啟時，管道中的流速增加，壓力降低，產(chǎn)生負水錘，其相對值用y表示，用同樣的方法可求出各相末計算公式。此時，求出的y本身為負值。(9-17)……(9-18)上述水錘壓力計算公式的條件：(1)沒有考慮管道摩阻的影響，因此只適用于不計摩阻(如水頭較高、管道較短等)的情況；(2)采用了孔口出流的過流特性，只適用于沖擊式水輪機，對反擊式水輪機必須另作修改；(3)這些公式在任意開關(guān)規(guī)律下都是正確的，可以用來分析非直線開關(guān)規(guī)律對水錘壓力的影響。三、水錘波在水管特性變化處的反射水錘發(fā)生后，水錘波在水管末端和水管特性變化處(水管進口、分岔、變徑段、閥門等)都要發(fā)生反射。當入射波到達水管特性變化處之后，一部分以反射波的形式折回，一部分以透射波的形式繼續(xù)向前傳播。反射波與入射波的比值稱反射系數(shù)，以r表示。透射波與入射波的比值稱透射系數(shù)，以s表示，兩者的關(guān)系為(9-19)(一)水錘波在水管末端的反射水錘波在水管末端的反射特性取決于水管末端的出流規(guī)律。對于水斗式水輪機，其噴嘴的出流規(guī)律為，當時，可近似地取為。在入射波未達到的時刻，，。設有一入射波傳到閥門后發(fā)生反射，產(chǎn)生一反射波折回，由方程(9-7)得：閥門處的水錘壓力為入射波與反射波的疊加結(jié)果，根據(jù)式(9-6)以上二式消去，簡化后得閥門的反射系數(shù)為(9-20)根據(jù)水錘常數(shù)和任意時刻的開度，可利用式(9-20)確定閥門在任意時刻的反射系數(shù)。當閥門完全關(guān)閉時，，閥門處發(fā)生同號等值反射。(二)水錘波在管徑變化處的反射如圖9-4所示的變徑管，入射波F1從1管傳來，在變徑處發(fā)生反射。反射波為f1，透射波為F2，由方程(9-6)和(9-7)及水流在變徑處的連續(xù)性，可推導出反射系數(shù)(9-21)圖9-4變徑管式中，。圖9-4變徑管r為正表示反射是同號的，其結(jié)果是使管1中水錘壓力的絕對值增大；反之，r為負表示反射是異號的，其結(jié)果是使水管1中水錘壓力的絕對值減小。若管2斷面趨近于零，則ρ2→∞，r=1，為同號等值反射，這相當于水管末端閥門完全關(guān)閉情況。若管2斷面為無限大，則V2=0，ρ2=0，r=-1，為異號等值反射，這相當于水庫處的情況。(三)水錘波在分岔處的反射如圖9-5所示，入射波F1從1管傳來，在分岔處發(fā)生反射，反射波為f1，透射波為F2和F3，根據(jù)基本方程(9-6)和(9-7)及此處水流的連續(xù)性，導出反射系數(shù)為圖9-5分岔管(9-22)圖9-5分岔管式中，Q為總管流量，A為水管斷面積。四、開度依直線變化的水錘進行水錘計算，最重要的是求出其最大值。在開度依直線規(guī)律變化情況下，不必用連鎖方程求出各相末水錘，再從中找出最大值，可用簡化方法直接求出。(一)開度依直線變化的水錘類型當閥門開度依直線規(guī)律變化時，根據(jù)最大壓力出現(xiàn)的時間可歸納為兩種類型：(1)最大水錘壓力出現(xiàn)在第一相末，，如圖9-6(a)，稱為第一相水錘。(2)最大水錘壓力出現(xiàn)在第一相以后的某一相，其特點是最大水錘壓力接近極限值，即>，如圖9-6(b)，稱為極限水錘。產(chǎn)生這兩種水錘現(xiàn)象的原因是由于閥門的反射特性不同造成的，閥門處的反射特性可由其反射系數(shù)確定。1．第一相水錘根據(jù)式(9-20)，當<1時，r為正，水錘波在閥門處的反射為同號。在閥門關(guān)閉過程中，閥門處任意時刻的水錘壓力由三部分組成：閥門不斷關(guān)閉所產(chǎn)生的升壓波、經(jīng)水庫反射回來壓力波、經(jīng)閥門反射向上游的壓力波。(1)第一相中，根據(jù)水庫異號反射的特性，升壓波到達水庫后反射回的降壓波還未到達閥門處，因此該處水錘壓力即是閥門關(guān)閉所產(chǎn)生的升壓波，在第一相末達到。圖9-6開度為直線關(guān)閉時的水擊類型(2)第二相末，水庫傳來的降壓波到達閥門處，如果此時閥門處具有同號反射的特性，則在該處反射仍為降壓波，兩個降壓波之和將超過第二相中由于閥門關(guān)閉所產(chǎn)生的升壓波，因而第二相末的水錘壓力<。圖9-6開度為直線關(guān)閉時的水擊類型(3)第三相末，由于第二相中閥門同號反射回去的降壓波，經(jīng)水庫異號反射為升壓波，這兩個升壓波共同作用，又使閥門處的水錘壓力開始升高，>。根據(jù)閥門同號反射的規(guī)律，水錘壓力將環(huán)繞某一值上下波動，最后趨于。由于最大水錘壓力出現(xiàn)在第一相末，>，故稱為第一相水錘。2．極限水錘根據(jù)式(9-20)，當>1時，r為負，水錘波在閥門處的反射為異號。在閥門關(guān)閉過程中，閥門處任意時刻的水錘壓力仍由上述三部分組成。第一相末，水庫反射回的降壓波還未到達閥門處，該處水錘壓力只是閥門關(guān)閉所產(chǎn)生的升壓波，即。第二相末，水庫傳來的降壓波到達閥門處，因閥門處為異號反射，則在該處反射為升壓波，它和在第二相中閥門繼續(xù)關(guān)閉產(chǎn)生的升壓波共同作用，使第二相中閥門處的水錘壓力繼續(xù)升高，使>。在以后各相，閥門處水錘壓力逐漸增加，趨近某一極限值。由于最大水錘壓力為，>，故稱為極限水錘。(二)開度依直線變化時水錘的簡化計算當調(diào)節(jié)閥門按直線規(guī)律啟閉，與的關(guān)系為：當閥門關(guān)閉時(9-23)當閥門開啟時(9-24)1．第一相水錘計算的簡化公式當<0.5時，，則(9-14)可簡化為：令，稱為水錘特性常數(shù)，關(guān)閉時用正值，開啟時為負值?？紤]到和的關(guān)系，代入上式可解得第一相末水錘壓力值為：關(guān)閉閥門時(9-25)開啟閥門時(9-26)發(fā)生第一相水錘的條件是<1，對于丟棄負荷情況，=1，有。若a=1000m/s，Vmax=5m/s，則H0>250m，故在丟棄負荷的情況下，只有高水頭電站才有可能出現(xiàn)第一相水錘。(2)極限水錘計算簡化公式根據(jù)(9-15)，第n相和第n+1相末的水錘壓力計算公式為：上二式相減，得：如果水錘波傳播的相數(shù)n足夠多，可認為，上式可以簡化為設，上式可寫為：解得：(9-27)當水錘壓力≤0.5時，，可得到更為簡化的近似公式：(9-28)(9-29)(3)間接水錘類型的判別條件僅用大于還是小于1作為判別水錘類型的條件是近似的。水錘的類型除與有關(guān)，還與有關(guān)。很明顯，這兩種情況的分界條件必須是。將式(9-14)的值用代替，得將式代入上式，則以、代人上式得將上式代入式中，解得值為(9-30)如果公式(9-30)滿足，則。公式(9-30)代表一根曲線，如圖9-7所示。圖9-7水錘類型圖中同時繪出了的直線。曲線表示極限水錘和第一相水錘的分界線，直線表示第一相水錘和直接水錘的分界線。共有五個分區(qū)：I區(qū)為極限正水錘；II區(qū)為第一相正水錘；III區(qū)為直接水錘；IV區(qū)為極限負水錘；V區(qū)為第一相負水錘。簡單判別方法：<1.0時，常發(fā)生第一相水錘；>1.5時，常發(fā)生極限水錘；1.0<<1.5時，則隨值的不同而發(fā)生第一相或極限水錘，個別情況下發(fā)生直接水錘。此時按圖9-7判別。最后，為了方便水錘壓力的計算，將計算公式匯總于表10-1。五、起始開度對水錘的影響水電站可能在各種不同的負荷情況下運行，當機組滿負荷運行時，起始開度＝1；當機組只擔任部分負荷運行時，＜l。因此機組由于事故丟棄負荷時的起始開度可能有各種數(shù)值。從前面的水錘壓力計算公式可以繪制出圖9-8。圖中的曲線和分界點說明了起始開度對水錘壓力的影響。表10-1水錘壓力計算公式匯總表開度計算公式近似公式起始終了關(guān)閉直接水擊010間接水擊01開啟直接水擊101間接水擊11011圖9-8起始開度對水錘壓力的影響由極限水錘只與有關(guān)，而與無關(guān)，圖中是一根平行于軸的水平線。對第一相水錘，隨著的減小而增大，所以在圖中表示為一根曲線。對直接水錘，，為一通過坐標軸原點的直線，其斜率為2。圖中三條曲線的交點為：直接水錘和第一相水錘：令和相等，可以解出：第一相水錘和末相水錘令和相等，可以解出：因此可得出以下結(jié)論：(l)當起始開度，>1時，，最大水錘壓力發(fā)生在閥門關(guān)閉的終了，即極限水錘；(2)當起始開度時，最大水錘壓力發(fā)生在第一相末；(3)當起始開度時，發(fā)生直接水錘。但由于直接水錘壓力的大小與初始開度成正比，所以不一定是最大的水錘值；(4)當閥門起始開度為臨界開度時，發(fā)生最大直接水錘，由得：水輪機存在空轉(zhuǎn)流量Qxx、相應的空轉(zhuǎn)開度為τxx、水輪機在該開度下運行，不能輸出功率，能量僅消耗于克服摩阻。因此，機組不可能在小于τxx開度下運行。如果τxx＞σ/ρ，說明該機組不可能發(fā)生直接水錘。另外，閥門實際關(guān)閉規(guī)律并非直線，根據(jù)水輪機調(diào)運器特性，關(guān)閉終了時存在延緩現(xiàn)象，小初始開度時的實際關(guān)閉時間要長于τ0Ts，水錘壓力比計算值要小，一般不起控制作用。六、開度變化規(guī)律對水錘壓力的影響前面有關(guān)第一相或極限水錘的一些概念及計算公式是在假定閥門開度按直線變化條件求得的。在水電站運行實踐中，閥門的啟閉不完全是按直線而往往采用非直線的規(guī)律。圖9-9繪出了三種不同的關(guān)閉規(guī)律，三種規(guī)律都具有相同的關(guān)閉時間，同時繪出了與之相應的三種水錘壓力變化過程線。由圖可以看出，開度的變化規(guī)律不同，水錘壓力的變化過程也不同。圖9-9開度變化規(guī)律對水錘壓力的影響曲線Ⅱ表示開始階段關(guān)閉速度較快，因此水錘壓力迅速上升到最大值，而后關(guān)閉速度減慢，水錘壓力逐漸減??；曲線Ⅲ的規(guī)律與曲線Ⅱ相反，關(guān)閉速度是先慢后快，而水錘壓力是先小后大。水錘壓力的上升速度隨閥門的關(guān)閉速度的加快而加快，最大壓力出現(xiàn)在關(guān)閉速度較快的那一時段末尾。從圖中可以看出，關(guān)閉規(guī)律Ⅰ較為合理，最不利的是規(guī)律Ⅲ。由此可見，通過調(diào)速器或針閥等設備，采取比較合理的啟閉規(guī)律，可以作為減小水錘壓力和解決調(diào)節(jié)保證問題的措施之一。在高水頭電站中常發(fā)生第一相水錘，可以采取先慢后快的非直線關(guān)閉規(guī)律，以降低第一相水錘值；在低水頭水電站中常發(fā)生極限水錘，可采取先快后慢的非直線關(guān)閉規(guī)律，以降低末相水錘值。七、水錘壓力沿管長的分布以上討論的都是水管末端A點(閥門或?qū)~處)的水錘壓力。在進行壓力管道強度設計時，不僅需要計算管道末端的壓力，而且需要管道沿線各點的最大正水錘壓力和最大負水錘壓力的分布情況，以便進行管道的強度設計及檢驗管道內(nèi)部是否有發(fā)生真空的可能。圖9-10水錘壓力沿管道的分布第一相水錘和極限水錘沿管長的分布規(guī)律是不同的，下面分別予以討論。(一)極限水錘壓力的分布規(guī)律理論研究證明，極限水錘無論是正、負水錘，管道沿線線的最大水錘壓力均按直線規(guī)律分布，如圖9-10中實線所示。若管道末端A點的最大水錘為和，則任意點C點的最大水錘為(9-32)和(9-33)(二)第一相水錘壓力的分布規(guī)律研究證明，第一相水錘壓力沿管線不按直線規(guī)律分布，正水錘壓力分布曲線是向上凸的，負水錘壓力分布曲線是往下凹的，如圖9-10中的虛線所示。任意點C的最大水錘升壓值發(fā)生在A點的最大水錘升壓傳到C點時，即比A點出現(xiàn)最大水錘升壓滯后(L-l)/a，其值為(9-34)式中為第一相末A點的水錘壓力，即可直接用簡化公式求得；為第一相終了前2l/a秒時A點的水錘壓力，可用第一相水錘簡化公式求得，只需用代替式中的即可。式(9-34)的近似表達式為(9-35)式中；上面的兩式可以看出，等號右端的第一項為管長為L時A點第一相末的水錘壓力，第二項為管長為L-l(相當于水庫移至C點)時A點第一相末的水錘壓力，C點最大水錘壓力為兩者之差。對于第一相負水錘，任意點C的最大水錘降壓為(9-36)式中為閥門開啟2l/a時A點的負水錘，可用表10-1中的公式求解，用代替式中的即可。相當于管長為l(即閥門移至C點)時第一相水錘。式(9-36)可近似表示為(9-37)式中繪制水錘壓力沿管線分布圖時，應根據(jù)管線的布置情況，選擇幾個代表性的斷面，求出各斷面上的最在正、負水錘壓力。當丟棄負荷時可不計管路的水損失，在上游最高靜水位上繪制水錘壓力分布圖；當增加負荷時，必須計算開啟終了時管路的水頭損失與流速水頭，在上游最低水位線以下，考慮水頭損失、流速水頭與負水錘壓力，繪制水錘壓力分布圖。第五節(jié)復雜管道水錘計算前面所討論的是簡單管道的水錘問題。簡單管的直徑、管壁厚度和管材料均不隨管長而變化，因此整根水管的特性是不變的。在實際工程中，這種簡單管是不多見的，常見的是復雜管路系統(tǒng)，共有三種類型：(1)管壁厚度、直徑和材料隨水頭增加自上而下逐段改變，這種復雜管稱為串聯(lián)管。(2)分岔管，這在分組供水和聯(lián)合供水中經(jīng)常遇到。(3)裝有反擊式水輪機的管道系統(tǒng)，應考慮蝸殼和尾水管的影響，而且其過流特性與孔口出流不一樣，流量不僅與作用水頭有關(guān)，而且與水輪機的機型和轉(zhuǎn)速有關(guān)。一、串聯(lián)管水錘的簡化計算圖9-11串聯(lián)管示意圖由于串聯(lián)管各管段的V0和a不同(如圖9-11)，因此表示水管特性的系數(shù)和各異。在實用中常把串聯(lián)管轉(zhuǎn)化為等價的簡單管來計算。所謂等價就是將串聯(lián)管轉(zhuǎn)化為簡單管后應滿足管長、相長和管中水體動能等與原管相同的原則。這種簡化計算方法稱為“等價水管法”。圖9-11串聯(lián)管示意圖設一根串聯(lián)管的管道特性為：L1，V1，a1；L2，V2，a2；……；Ln，Vn，an，等價管的總長為：。根據(jù)管中水體動能不變的要求，則LVm=L1V1+L2V2+……LnVn=∑LiVi，由此可得加權(quán)平均流速：(9-39)根據(jù)相長不變的要求，水錘波按平均波速由斷面閥門A傳到水庫斷面D所需的時間等于水錘波在各段傳播時間的總和，即由此可得波速的加權(quán)平均值：(9-40)對于間接水錘，管道的平均特性常數(shù)為(9-41)(9-42)(9-43)求出管道平均特性常數(shù)后，可按簡單管的間接水錘計算公式求出復雜管道的間接水錘值。二、分岔管的水錘壓力計算如圖9-12所示，分岔管除了管徑和管壁厚度沿管軸線變化外，同時還增加了分岔，其水錘壓力計算比串聯(lián)管更復雜。分岔管的水錘計算方法之一是截肢法。這種方法的特點是：當機組同時關(guān)閉時，選取總長度最大的一根支管，如圖9-12(a)中的支管2，將其余的支管截掉，變成圖(b)所示的串聯(lián)管道，然后用各管段中實際流量求出各管段的流速，再用加權(quán)平均的方法求出串聯(lián)管中的平均流速和平均波速，最后采用串聯(lián)管的簡化公式相應地求出水錘值。當壓力水管的主管較長、支管較短(例如支管長度為主管的10％以內(nèi))的情況下，計算結(jié)果誤差不大，否則誤差較大。(a)(b)圖9-12分岔管的截肢法三、蝸殼、尾水管水錘壓力計算反擊式水輪機的過流部件包含有蝸殼和尾水管。蝸殼和尾水管中的水流現(xiàn)象十分復雜，水錘基本方程主要假定之一是水流為一元流，這一假定對蝸殼和尾水管是不合適的，因此蝸殼和尾水管中的水錘計算一般只能用近似方法。首先將蝸殼視作壓力水管的延續(xù)部分，并假想把導葉移至蝸殼的末端，尾水管也作為壓力管道的一部分。這樣把壓力管道、蝸殼和尾水管組合視為一串聯(lián)管，再將該串聯(lián)管簡化為等價簡單管進行計算。設壓力水管、蝸殼及尾水管長度、平均流速和水錘波速分別為LT、VT、aT；Lc、Vc、ac；Lb、Vb、ab，則L=LT+Lc+LbVm=(LTVT+LcVc+LbVb)/L于是可求出等價管和特性系數(shù)、，求出管道末端最大水錘壓力值。然后以管道、蝸殼、尾水管三部分水體動能為權(quán)，將水錘壓力值進行分配，求出壓力管道、蝸殼末端和尾水管進口的水錘壓力。壓力水管末端最大壓力上升相對值為：(9-44)蝸殼末端最大水錘壓力上升相對值：(9-45)尾水管在導葉或閥門之后，水錘現(xiàn)象與壓力管道相反，其進口處壓力下降相對值為：(9-46)求出尾水管的負水錘后，應校核尾水管進口處的真空度Hr，以防水流中斷。(9-47)式中Hs—水輪機的吸出高度Vb—尾水管進口斷面在出現(xiàn)yb時的流速。第六節(jié)水錘計算的計算機方法根據(jù)簡化水錘方程即數(shù)學物理中的波動方程導出的水錘計算連鎖方程曾廣泛用于計算管道水錘壓力，其缺點是不能用于分析復雜管路和復雜邊界的水錘，并且不能計入管道摩擦阻力的影響。計算機的飛速發(fā)展和應用研究帶來了計算上的革命，Gray和Streeter合作首先介紹了用計算機計算管道水錘的特征線法，隨后Streeter出版了瞬變流專著奠定了用計算機分析管道水錘的基礎。用特征線法計算水錘可分析復雜管路也可處理復雜的邊界條件，也可以計入摩擦阻力的影響(在低水頭水電站中摩擦阻力的影響較大)，下面主要介紹特征線法的計算機算法。一、特征線方程特征線方法是將偏微分方程轉(zhuǎn)化為全微分方程的型式，再對全微分方程進行積分，得到有限差分方程進行數(shù)值計算。首先對本章第三節(jié)介紹的水錘基本方程(9-2)和(9-3)進行適當簡化。假設管道是水平的，且沿管道長度引水管的直徑不變。另外，水錘的發(fā)生和衰減過程是在很短的時間內(nèi)完成的，所以在式(9-2)中，式(9-3)中。這樣可以得到簡化以后的水錘基本方程，分別命名為L1和L2：(9-48)(9-49)引入特征值λ，將上面的兩個方程進行線性組合，得：(9-50)將其整理為：(9-51)特征線方法就是選擇兩個不同的實數(shù)特征值和，使得方程(9-51)成為一組全微分方程，并與方程(9-48)和(9-49)完全等價。設方程(9-51)的解為和，則：，(9-52)對比方程(9-51)和(9-52)，假如下面的關(guān)系成立：(9-53)則方程(9-51)可以轉(zhuǎn)化為全微分方程：(9-54)并且由式(9-53)可以得出：(9-55)及(9-56)式(9-56)說明，壓力管道中的水壓力以波的型式傳播，其傳播速度為a。當其取正值時，水錘壓力波向水庫方向傳播，取負值時水錘波向水輪機方向傳播，壓力管道中的水錘壓力就等于這兩種波的疊加。在發(fā)生水錘的過程中，壓力管道中的水壓力分布不僅與時間有關(guān)，而且與位置有關(guān)，這是由于水錘波在管道中來回傳播，管壁的阻力可以使水錘波逐漸減弱，而波的傳播與疊加使得不同位置的壓力也不盡相同。當特征值λ分別取正值和負值時，將其代入方程(9-54)，可以得到兩組方程，分別用C+和C-來命名，即：(9-57a)(9-57b)將上述方程的解在x-t平面上展開，就不難對它加以形象化說明。因為對于一個給定的管道，a通常是常數(shù)，于是方程(9-57a)在x-t平面上畫出來是一根直線AP；同樣，方程(9-57b)在x-t平面上是另一根直線BP，見圖9-13。我們將這些x-t平面上斜率為±1/a的直線分別稱為正特征線和負特征線。沿C+特征線，方程(9-57a)成立；沿C-特征線，方程(9-57圖9-13特征線圖9-14單一管道求解的x-t網(wǎng)絡圖二、基本求解方法為了用有限差分法求解常微分方程(9-57)，首先將管道在長度方向離散成N等份，每一等份的長度為，每隔時間計算一次水錘壓力的分布，則在長度方向和時間方向的離散可以形成一個計算網(wǎng)格，如圖9-14所示。如果計算的時間步長取為，則網(wǎng)格的對角線斜率為+1/a或-1/a，即滿足方程(9-57a)或(9-57b)中的第二個方程。如果A點的變量v和H是已知的，那么沿著C+方向的特征線從A到P進行積分，同時注意到adt/g=dx/g，并用流量Q代替流速v，可以得到：(9-58)同理，對方程組C-，-adt/g=dx/g，可得：(9-59)式中所包含的Q|Q|是與x有關(guān)的變量，在近似計算中可用A點或B點的值表示。若取得足夠小，其一次近似即可滿足要求，代入A、B值得：(9-60a(9-60b)在上面兩個方程中，A點和B點的變量值是已知的，而未知量只有HP和QP，兩個方程聯(lián)立可以求解之。發(fā)生水錘過程之前或發(fā)生之初(即t=0時)，管道中的水流呈穩(wěn)定流狀態(tài)，各點的H、V是已知的。在時刻，管道中任一點的流動狀態(tài)可由式(9-60a)和式(9-60b)解出，進而可以再對時刻的流動狀態(tài)進行計算。但需要注意的是，對管道兩端的邊界點，由于只能利用式(9-60a)和式(9-60b)中的一個方程，所以還必須應用管道的邊界條件才能求解。引入流量與流速的關(guān)系Q=VA(A為管道斷面積)，并將計算過程中與管道特性有關(guān)的常數(shù)進行簡化，令：則方程(9-60)可以寫成:(9-61a)(9-61b)對于特征網(wǎng)格上的任意截面i點，上述兩個方程可以改寫為：(9-62a)(9-62b)令(9-63a)(9-63b)代入(9-62)得：(9-64a(9-64b)因此可以求解出：(9-65a(9-65b)或(9-65c)觀察圖9-14中的網(wǎng)格，可以看到，系統(tǒng)中的兩個端點，從第一時步以后，開始影響內(nèi)部的點。所以，為了求得任意時刻的解，必須引入相應的邊界條件。水電站有壓引水系統(tǒng)的邊界條件見本章第三節(jié)。水及計算的步驟總結(jié)為：(1)確定計算時間步長Δt。由于采用矩形網(wǎng)格進行計算，故一般取Δt=Δx/a?？紤]到水錘波速a是確定的，所以關(guān)鍵在于選定Δx。通?？筛鶕?jù)管道布置及精度要求將整個管路系統(tǒng)分成很多管段，各管段的兩端或為內(nèi)點，或為邊界點。由于波速隨管道特性而變化，而Δt又是常數(shù)，所以不同管道的管段長Δx是不相同的。另外，從數(shù)學上可以證明，只有當時，差分計算格式才是穩(wěn)定的。(2)計算各節(jié)點在恒定流狀態(tài)下(即起始狀態(tài))的水壓力分布和流量值。(3)增加一個Δt，按上述所列的公式計算該時刻管道各內(nèi)部節(jié)點處的水頭和流量。(4)計算同一時刻水輪機處的水頭，流量。三、計算程序根據(jù)上面的計算原理，編制了一個簡單管道的水錘壓力計算程序，主要目的是說明計算程序的編制方法。其主要變量說明：NUMDX：在管道長度方向的節(jié)點數(shù)；NUMDT：時間方向的分步數(shù)；MTYPE1：水庫端的邊界條件代碼，1為恒定水位，2為正弦變化水位，3為半正弦流量變化；MTYPE2：水輪機端邊界條件代碼，1為孔口出流，2為封閉端，3為反擊式水輪機(待增加)；TCLOSE：水輪機關(guān)閉時間；DX：管道長度離散的步長；PIPEL：管道長度；PIPED：管道直徑；QBEGN：初始流量；QEND：關(guān)閉結(jié)束后的流量；HEAD0：水輪機的凈水頭(以水輪機出流孔口處為0)；COFLOS：管道的糙率系數(shù)；WSPEED：水錘波傳播速度；QPUMP：模擬半正弦流量變化Q=Q0+ΔQ|sinωt|中的Q0；DELTQ：模擬半正弦流量變化Q=Q0+ΔQ|sinωt|中的ΔQ；DELTH：模擬正弦變化水位的ΔH；XCOOR(50)：管道長度方向的節(jié)點坐標；HEADS(50,1000)：管道內(nèi)的水頭分布和隨時間的變化；FLOWQ(50,1000)：管道內(nèi)的流量分布和隨時間的變化。源程序如下：C###################################################C簡單壓力管道水錘計算程序C特征線方法C2003年5月C###################################################CDIMENSIONHEADS(50,1000),FLOWQ(50,1000),XCOOR(50)CHARACTER*15AFILE,BFILE,TITLE*80COMMON/PARAM1/NUMDX,NUMDT,MTYPE1,MTYPE2COMMON/PARAM2/TCLOSE,DT,PIPEL,PIPED,QBEGN,QENDCOMMON/PARAM3/HEAD0,COFLOS,WSPEEDCOMMON/PARAM4/GRAVC,SECTA,VINIT,HLOSSCOMMON/PARAM5/QPUMP,DELTH,OMEGA,TBEGN,TEND,DELTQDATAGRAVC/9.81/PI/3.1415926/CC打開輸入和輸出數(shù)據(jù)文件CWRITE(*,2350)READ(*,'(A15)')AFILEOPEN(10,FILE=AFILE)WRITE(*,2360)READ(*,'(A15)')BFILEOPEN(20,FILE=BFILE)2350FORMAT(/'PleasetypetheDATAfilename....',$)2360FORMAT(/'Pleasetypetheresultfilename....',$)CC基本數(shù)據(jù)的輸入和輸出CREAD(10,'(A80)')TITLEREAD(10,1000)NUMDX,NUMDT,MTYPE1,MTYPE2READ(10,1010)TCLOSE,PIPEL,PIPED,QBEGN,QENDIF(MTYPE1.EQ.1)READ(10,1010)HEAD0,COFLOS,WSPEEDIF(MTYPE1.EQ.2)READ(10,1010)HEAD0,COFLOS,WSPEED,DELTHIF(MTYPE1.EQ.3)READ(10,1010)HEAD0,COFLOS,WSPEED,QPUMP,DELTQ1000FORMAT(10I5)1010FORMAT(7F10.0)CWRITE(20,2000)TITLEWRITE(20,2010)NUMDX,NUMDT,MTYPE1,MTYPE2WRITE(20,2020)TCLOSE,PIPEL,PIPED,QBEGN,QENDIF(MTYPE1.EQ.1)WRITE(20,2030)HEAD0,COFLOS,WSPEEDIF(MTYPE1.EQ.2)WRITE(20,2040)HEAD0,COFLOS,WSPEED,DELTHIF(MTYPE1.EQ.3)WRITE(20,2050)HEAD0,COFLOS,WSPEED,QPUMP,DELTQ2000FORMAT(//A80)2010FORMAT(/'NUMDX=',I5/'NUMDT=',I5/'MTYPE1=',I5/#'MTYPE2=',I5)2020FORMAT(/'TCLOSE=',F12.3/#'PIPEL=',F12.3/'PIPED=',F12.3/#'QBEGN=',F12.3/'QEND=',F12.3)2030FORMAT(/'HEAD0=',F12.3/'COFLOS=',F12.3/#'WSPEED=',F12.3)2040FORMAT(/'HEAD0=',F12.3/'COFLOS=',F12.3/#'WSPEED=',F12.3/'DELTH=',F12.3)2050FORMAT(/'HEAD0=',F12.3/'COFLOS=',F12.3/#'WSPEED=',F12.3/'QPUMP=',F12.3/#'DELTQ=',F12.3)CC幾個計算常數(shù)CCLOSS=PIPEL*COFLOS*COFLOS/(PIPED/4.0)**(4.0/3.0)SECTA=(PI*PIPED*PIPED)/4.0VINIT=QBEGN/SECTAHLOSS=CLOSS*VINIT*VINITDX=PIPEL/NUMDXNUMDX=NUMDX+1DT=DX/WSPEEDDOI=1,NUMDXX=(I-1)*DXXCOOR(I)=XENDDOWRITE(20,'(A10,F12.3)')'DT=',DTWRITE(20,2060)(XCOOR(I),I=1,NUMDX)2060FORMAT(/'節(jié)點坐標'/(6F12.3))CC初始狀態(tài)下的流量和水頭分布CX0=XCOOR(1)DO10I=1,NUMDXX=XCOOR(I)DX=X-X0HEADS(I,1)=HEAD0-HLOSS*DX/PIPELFLOWQ(I,1)=QBEGN10CONTINUECC開始計算CDO100IT=2,NUMDTTIME=(IT-1)*DTDO90IX=2,NUMDX-1QL=FLOWQ(IX-1,IT-1)QR=FLOWQ(IX+1,IT-1)HL=HEADS(IX-1,IT-1)HR=HEADS(IX+1,IT-1)B=WSPEED/(GRAVC*SECTA)DXL=XCOOR(IX)-XCOOR(IX-1)R1=COFLOS*DXL/(2*GRAVC*PIPED*SECTA*SECTA)DXR=XCOOR(IX+1)-XCOOR(IX)R2=COFLOS*DXR/(2*GRAVC*PIPED*SECTA*SECTA)CPI=HL+B*QL-R1*QL*ABS(QL)CMI=HR-B*QR+R2*QR*ABS(QR)HI=(CPI+CMI)/2.0QI=(HI-CMI)/BHEADS(IX,IT)=HIFLOWQ(IX,IT)=QI90CONTINUECC邊界點的計算CCALLBOUNDL(IT,TIME,HEADS,FLOWQ,XCOOR)CALLBOUNDR(IT,TIME,HEADS,FLOWQ,XCOOR)C100CONTINUECC輸出結(jié)果CWRITE(20,3000)DO200IT=1,NUMDTTIME=IT*DTWRITE(20,3010)IT,TIMEWRITE(20,3020)(HEADS(I,IT),I=1,NUMDX)WRITE(20,3030)IT,TIMEWRITE(20,3020)(FLOWQ(I,IT),I=1,NUMDX)200CONTINUECWRITE(20,3100)DO300IX=1,NUMDXWRITE(20,3110)IX,XCOOR(IX)WRITE(20,3020)(HEADS(IX,I),I=1,NUMDT)WRITE(20,3120)IX,XCOOR(IX)WRITE(20,3020)(FLOWQ(IX,I),I=1,NUMDT)300CONTINUESTOPC3000FORMAT(/3X,'水錘計算結(jié)果'//3X,'每個時段沿管長的分布')3010FORMAT(/3X,'時步(IT)=',I4,5X,'時間(TIME)=',F12.4//3X,'水頭分布')3020FORMAT(6F13.4)3030FORMAT(/3X,'時步(IT)=',I4,5X,'時間(TIME)=',F12.4//3X,'流量分布')3100FORMAT(/3X,'每個節(jié)點的變化時程曲線')3110FORMAT(/3X,'節(jié)點號(IX)=',I5,5X,'坐標(X)=',F12.3//3X,'水頭分布')3120FORMAT(/3X,'節(jié)點號(IX)=',I5,5X,'坐標(X)=',F12.3//3X,'流量分布')ENDCSUBROUTINEBOUNDL(IT,TIME,HEADS,FLOWQ,XCOOR)CC左端(水庫端)邊界條件CDIMENSIONHEADS(50,1),FLOWQ(50,1),XCOOR(1)COMMON/PARAM1/NUMDX,NUMDT,MTYPE1,MTYPE2COMMON/PARAM2/TCLOSE,DT,PIPEL,PIPED,QBEGN,QENDCOMMON/PARAM3/HEAD0,COFLOS,WSPEEDCOMMON/PARAM4/GRAVC,SECTA,VINIT,HLOSSCOMMON/PARAM5/QPUMP,DELTH,OMEGA,TBEGN,TEND,DELTQCCCONSTANTSCB=WSPEED/(GRAVC*SECTA)DXR=XCOOR(2)-XCOOR(1)QR=FLOWQ(2,IT-1)HR=HEADS(2,IT-1)R2=COFLOS*DXR/(2*GRAVC*PIPED*SECTA*SECTA)CMI=HR-B*QR+R2*QR*ABS(QR)CGOTO(100,200,300),MTYPE1CCH=H0C100HEAD=HEAD0HEADS(1,IT)=HEADFLOWQ(1,IT)=(HEAD-CMI)/BRETURNCCH=H0+DH*SIN(Wt)C200HEAD=HEAD0+DELTH*SIN(OMEGA*TIME)HEADS(1,IT)=HEADFLOWQ(1,IT)=(HEAD-CMI)/BRETURNCCQ=Q0+DQ*|SIN(Wt)|C300DELTQ=DELTHFLOWQ(1,IT)=QPUMP+DELTQ*ABS(SIN(OMEGA*TIME))HEADS(1,IT)=CMI+B*FLOWQ(1,IT)RETURNENDCSUBROUTINEBOUNDR(IT,TIME,HEADS,FLOWQ,XCOOR)CC右端(水輪機端)邊界條件CDIMENSIONHEADS(50,1),FLOWQ(50,1),XCOOR(1)COMMON/PARAM1/NUMDX,NUMDT,MTYPE1,MTYPE2COMMON/PARAM2/TCLOSE,DT,PIPEL,PIPED,QBEGN,QENDCOMMON/PARAM3/HEAD0,COFLOS,WSPEEDCOMMON/PARAM4/GRAVC,SECTA,VINIT,HLOSSCOMMON/PARAM5/QPUMP,DELTH,OMEGA,TBEGN,TEND,DELTQCB=WSPEED/(GRAVC*SECTA)DXL=XCOOR(NUMDX)-XCOOR(NUMDX-1)R1=COFLOS*DXL/(2*GRAVC*PIPED*SECTA*SECTA)CGOTO(100,200,300),MTYPE2CC孔口出流C100IF(TIME.LT.TCLOSE)TAO=1.0-TIME/TCLOSEIF(TIME.GE.TCLOSE)TAO=0CV=(QBEGN*TAO)**2/(2.0*HEAD0)QL=FLOWQ(NUMDX-1,IT-1)HL=HEADS(NUMDX-1,IT-1)CPI=HL+B*QL-R1*QL*ABS(QL)FLOWQ(NUMDX,IT)=-B*CV+SQRT((B*CV)**2+2.0*CV*CPI)HEADS(NUMDX,IT)=CPI-B*FLOWQ(NUMDX,IT)RETURNCC封閉端C200FLOWQ(NUMDX,IT)=0QP=0NL=NUMDX-1QL=FLOWQ(NL,IT-1)HL=HEADS(NL,IT-1)HEADS(NUMDX,IT)=HL-B*(QP-QL)-R1*QL*ABS(QL)RETURNCC反擊式水輪機C300RETURNEND四、計算實例有一長400m的水輪機管道，直接從水庫引水。水輪機閥門在全開狀態(tài)時，管道內(nèi)水流流量56.55m3/s，凈水頭H0=120m。管道直徑4m，其糙率系數(shù)為0.012，水錘波速為1200m數(shù)據(jù)文件為：41900112.40.01400.04.0056.550.0120.00.0121200.00.010.020.030.040.050.060.070.080.090.0100.0110.0120.0130.0140.0150.0160.0170.0180.0190.0200.0210.0220.0230.0240.0250.0260.0270.0280.0290.0300.0310.0320.0330.0340.0350.0360.0370.0380.0390.0400.0計算結(jié)果：閥門處最大水錘壓力225.7691m，最小水錘壓力15.6038m，相當于水頭壓力升高105.7691圖9-15閥門處水壓力變化過程第六節(jié)機組轉(zhuǎn)速變化計算機組與電力系統(tǒng)解列后負荷變?yōu)榱悖藭r多余的能量轉(zhuǎn)化為機械能，使機組轉(zhuǎn)速上升，水輪機調(diào)節(jié)機構(gòu)開始關(guān)閉導葉，水輪機的引用流量逐漸減小，機組出力逐漸下降，同時在引水系統(tǒng)產(chǎn)生水錘壓力。當關(guān)閉到空轉(zhuǎn)開度時出力變?yōu)榱?。導葉關(guān)閉過程中所產(chǎn)生的能量，完全被機組轉(zhuǎn)動部分所消耗，造成機組轉(zhuǎn)速的升高。在機組調(diào)節(jié)過程中，轉(zhuǎn)速變化通常以相對值表示，稱為轉(zhuǎn)速變化率β，又稱為暫態(tài)不均衡率。丟棄負荷(9-66)增加負荷(9-67)式中：n0——機組額定轉(zhuǎn)速；nmax——丟棄負荷后的最高轉(zhuǎn)速。nmin——增加負荷后的最低轉(zhuǎn)速。一、機組運動方程機組作為剛體繞主軸旋轉(zhuǎn)，其運動方程為(9-68)式中——機組角速度；Mt——作用于機組上的動力矩；Mg——作用于機組上的阻力矩；M——不平衡力矩；J——機組轉(zhuǎn)動部分的慣性矩。負荷不變時，Mt=Mg，M=0，，轉(zhuǎn)速不變；機組丟棄負荷時，Mt>Mg，M>0，，轉(zhuǎn)速上升；增加負荷時，Mt<Mg，M<0，，轉(zhuǎn)速下降。二、機組轉(zhuǎn)速變化率計算近似公式圖9-16圖9-16丟棄負荷時出力與時間關(guān)系當丟棄負荷時，假定導葉按直線規(guī)律關(guān)閉，忽略了轉(zhuǎn)速變化、效率變化等因素的影響，且先不考慮水錘壓力升高對出力的影響，則水輪機的出力與時間成直線變化，如圖9-16所示。當丟棄負荷后在關(guān)閉時間Ts1內(nèi)產(chǎn)生多余的能量為：(9-69)式中Ts1——導葉關(guān)閉至空轉(zhuǎn)的時間；對于沖擊式和混流式水輪機Ts1=0.9Ts；對于軸流式水輪機Ts1=0.7Ts；102——單位轉(zhuǎn)換系數(shù)，由出力N0的千瓦數(shù)變?yōu)閗g.m/s；N0——機組丟棄負荷前的出力，以kW計。未被發(fā)電機輸出的能量E完全轉(zhuǎn)化為機組轉(zhuǎn)動部分的動能，使其轉(zhuǎn)速增加，其關(guān)系式為：(9-70)式中J——慣性矩，，G為轉(zhuǎn)動部分重量(t)，D是轉(zhuǎn)動部分慣性直徑(m)，如果以kg計，?！獧C組額定角速度，，n0是機組每分鐘的額定轉(zhuǎn)速?！獧C組最大角速度，將上述已知值和代入式(9-70)，得解上式并去掉不合理的根得：(9-71)在調(diào)節(jié)過程中由于水錘壓力影響，使水輪機出力增加或減少，因此需要考慮水錘修正系數(shù)f，其值根據(jù)圖9-17給出的曲線查得。丟棄負荷時：(9-72)增加負荷時：