P′条件下,如果上、下游水位高差Z很大,则水舌具有很大动能,容易把下游水体推开一段距离,发生远驱水跃,临近堰壁的下游处,水深仍小于堰顶,则发生自由出流。试验表明,薄壁堰发生淹没出流的充分条件是:同时满足下列两条件图7.5图7.6h下>P′及z/P′<0.7·198·
第7章堰流薄壁堰淹没出流的流量Q′可按下式近似计算H2n0.385Q′=Q[1-()](7.6)H1式中:Q表示堰前水头为H1时的自由出流流量;H2是下游水位高出堰顶的高度(见图7.5);n为指数,对于矩形堰n=3/2,对于三角堰n=5/2。2)直角三角形薄壁堰3当所需测量的流量较小,如Q<0.1m/s时,采用矩形薄壁堰则因水头过小,测量水头的相对误差增大,一般改用直角三角形薄壁堰,其计算公式为5/2Q=C0H式中:C0为三角堰的流量系数。当H=0.05~0.25m时,可用下列公式计算2.47Q=0.0154Hl/s式中H———水头,以厘米计。(2)实顶堰实用堰是水利工程中用来挡水同时又能泄水的建筑物,它的剖面形式是随着生产的发展而不断改进的。如采用不便加工成曲线的条石或其他当地材料修建的中、低溢流堰,堰顶剖面常作成折线形,称为折线形实用堰(图7.7(a)、(b))。如用混凝土修筑的中、高溢流堰,堰顶制成适合水流情况的曲线形,称为曲线型实用堰(图7.7(c)、(d))。图7.7曲线型实用堰又可分为非真空堰和真空堰两类。如果堰的剖面曲线基本上与薄壁堰的水舌下缘外形相符,水流作用在堰面上的压强仍近似为大气压强,称为非真空堰(图7.7(c))。若堰的剖面曲线低于薄壁堰的水舌的下缘,溢流水舌局部地脱离堰面,脱离处的空气被水流带走而形成真空区,这种堰称为真空堰(图7.7(d))。真空堰由于堰面上真空区的存在,与管嘴的水力性质相似,增加了堰的过流能力,即增大了流量系数。但是,由于真空区的存在,水流不稳定而引起水工构筑物的振动,且易使堰面发生空蚀破坏。为了防止这种情况发生,多将实用堰剖面外形稍稍伸入薄壁堰溢流水舌下缘以内,如图7.7(c)所示。实用堰的流量计算公式,可按分析薄壁堰的方法得出与式(7.2)相同的形式。若同时考虑·199·
水力学侧收缩和淹没出流的影响,可将(7.2)式右边乘以系数ε≤1及σ≤1,便得堰流流量计算的普遍公式3/2Q=εσmb2gH0(7.7)式中:ε为侧收缩系数,当无侧收缩时ε=1.0,当有侧收缩时ε<1.0并可用下式计算H0ε=1-α(7.8)b+H0其中,α为考虑坝墩头部形状影响的系数,矩形坝墩α=0.20,半圆形或尖形坝墩α=0.11,流线形尖墩α=0.06(图7.8)。图7.8图7.9σ为淹没系数,当自由出流时σ=1.0,当淹没出流时σ<1.0,其值与H2/H(见图7.9)有关,可根据表7.1查得。实用堰形成淹没出流的条件与薄壁堰相同。表7.1H2/HσH2/HσH2/HσH2/Hσ0.050.9970.520.930.740.8310.920.5700.100.9950.540.9250.760.8140.930.5400.150.9900.560.9190.780.7960.940.5060.200.9850.580.9130.800.7760.950.4700.250.9800.600.9060.820.7500.960.4210.300.9720.60.8970.840.7240.970.3570.360.9640.640.8880.860.6950.980.2740.400.9570.660.8790.880.6630.990.1700.460.9450.700.8560.900.6210.9950.100.500.9350.720.8440.910.5961.000.00(3)宽顶堰水利工程中,宽顶堰溢流的现象和形式是很常见的。如进水闸,不论有底坎或无坎(平底),其水流均属宽顶堰溢流;流经隧洞、涵管进口、小桥孔及施工围堰的水流等,水流流经这些构筑物时,因断面变小,流速增大,自由水面发生降落,在2.5<δ/Η<10范围内,具有宽顶堰的水力性质。1)流量公式及流量系数对堰前断面0-0至堰顶上收缩断面C-C(图7.10)的水流应用能量方程和连续条件按照前述方法可推导出宽顶堰流量公式为式(7.2)。将侧收缩及淹没出流的影响纳入公式中便得式(7.9):·200·
第7章堰流(7.9)式中,B为堰顶过水宽度。宽顶堰的流量系数m取决于堰的进口形式和堰的相对高度P/H,可按下列经验公式进行计算:堰坎进口为直角形(图7.11(a))图7.10m=0.32+0.013-P/H0.46+0.75P/H(7.10)图7.11堰坎进口为圆弧形(图7.11(b))3-P/Hm=0.36+0.01(7.11)P1.20+1.50H由以上两式可知,当P/H=3时,m为最小值;当P=0时,m=0.385,为最大值。m值的范围对于直角进口为m=0.32~0.385,当P/H>3时,m=0.32;对圆弧形进口为m=0.36~0.385,当P/H>3时,m=0.36。2)淹没出流条件及淹没系数σ由实验得知:宽顶堰的淹没条件为H2≥0.8H0(7.12)式中,H0为堰前水头,H2为下游水面高出堰顶的高度(图7.12)图7.12图7.13宽顶堰的淹没系数σ值主要随相对淹没度H2/H的增大而减小,可查表7.2。表7.2宽顶堰的淹没系数σH20.800.810.820.830.840.850.860.870.880.89H0σ1.000.9950.990.980.970.960.950.930.900.87H20.900.910.920.930.940.950.960.970.98H0σ0.840.820.780.740.700.640.590.500.40·201·
水力学3)侧收缩系数ε可按以下经验公式计算ε值:34Hbbε=1-α0·(1-)(7.13)0.2H+PBB式中,α0为堰的边墩或中墩的形状系数,矩形墩α0=0.19,圆形墩α0=0.1;b为溢流孔净宽;B为上游引渠宽;P为上游堰顶高(如图7.13)。例7.1有一矩形无侧收缩薄壁堰。已知堰宽B=0.5m,上、下游堰高P=P′=0.5m,堰前水头H=0.2m,求下游水深分别为h下=0.4m及h下=0.6m时通过薄壁堰的流量。解①求h下=0.4m时的流量因h下P′,下游水面高于堰顶。又上下游水位差z=P+H-h下=0.5+0.2-0.6=0.1mz0.1==0.2<0.7P′0.5满足薄壁堰淹没出流的两个条件,故为淹没出流。薄壁堰淹没出流的流量Q′可按式(7.6)计算H2n0.385Q′=Q[1-()]H1式中:n=1.5;H2=h下-P′=0.6-0.5=0.1m;H1=H=0.2m;Q为堰前水头H1=H=30.2m时的自由出流流量,即Q=0.0866m/s。所以0.11.50.3853Q′=0.0866[1-()]=0.0732m/s0.2例7.2一直角进口无侧收缩宽顶堰,堰宽b=3.0m,堰坎高P=P′=0.38m,堰前水头H=0.6m。求当下游水深分别为0.6m及0.9m时通过此堰的流量。解①求流量系数m因为P/H=0.38/0.6=0.633<3,所以流量系数m由式(7.10)计算:3-P/Hm=0.32+0.01×=0.3450.46+0.75P/H因为流量待求,行近流速v0和堰前总水头H0均为未知,故需要试算。②第一次近似计算:先假设不计行近流速水头并取α=1,即设H01=H=0.6m,H2=h下-P′=0.6-0.38=0.22m,·202·
第7章堰流0.8H01=0.8×0.6=0.48m,因H2<0.8H01,故按自由出流计算,无侧收缩,ε=1.0。3/23/23Q1=mb2gH01=0.345×3.0×2×9.8×0.6=2.130m/sQ12.130v01===0.724m/sb(H+P)3(0.6+0.38)③第二次近似计算用v0=v01计算堰前总水头:22v010.724H02=H+=0.6+=0.627m2g19.6因H2<0.8H02=0.8×0.627=0.502m,故仍为自由出流,其流量为3/23/23Q2=mb2gH02=0.345×3.0×19.6×0.627=2.273m/s2.273v02==0.773m/s3(0.38+0.6)则应有22v020.773H03=H+=0.6+=0.630m2g19.6④第三次近似计算3/2Q3=0.345×319.6×0.63=2.291mQ3-Q22.291-2.273||=||=0.8%<5%Q32.2913满足精度要求,则当下游水深为0.6m时通过此堰的流量Q≈Q3=2.291m/s。⑤当下游水深h下=0.9m时,先假设堰上总水头H01=H=0.6m,此时H2=h下-P′=0.9-0.38=0.52m,H2>0.8H01,故为淹没出流。由H2/H01=0.87查表7.2,得σ1=0.93,第一次近似计算流量得3/23/23Q1=σ1mb2gH01=0.93×0.345×3.019.6×0.6=1.98m/sQ1.98v01===0.674m/sb(P+H)3(0.38+0.6)22v010.674H02=H+=0.6+=0.623m2g19.6第二次近似计算流量H2/H02=0.52/0.623=0.83>0.8,故为淹没出流,查表7-2得σ2=0.98,所以3/23Q2=0.98×0.345×3.019.6×0.623=2.208m/s2.208v02==0.751m/s3(0.38+0.6)22v020.751H03=H+=0.6+=0.629m2g19.6第三次近似计算流量H2/H03=0.52/0.629=0.83>0.8,为淹没出流,查表7.2得σ3=0.98,所以3/2Q3=0.98×345×3.019.6×0.629=2.24m·203·
水力学Q3-Q22.24-2.208||=||=1.4%<5%Q32.243故所求流量为Q≈Q3=2.24m/s。思考题7.1简述堰流的水力计算特点。7.2堰有几种类型,如何判别?7.3简述宽顶堰的水流特点。7.4宽顶堰实现淹没出流的充要条件是什么?7.5堰流流量计算公式是如何推导出来的?习题7.1一直角进口无侧收缩宽顶堰,宽度b=2m,堰高P=0.5m,堰前水头为1.8m,设为自由出流,求通过堰的流量。7.2一矩形进口宽顶堰,堰宽b=2m,堰高P=P′=1m,堰前水头H=2m,上游渠宽B=3m,边墩为矩形。下游水深h下=2.8m,求过堰流量(设行近流速v0可忽略不计)。7.3一无侧收缩矩形薄壁堰,堰宽b=0.5m,堰高P=P′=0.4m,堰前水头H=0.6m,下游水深h=0.6m,求通过的流量。7.4一圆角进口无侧收缩宽顶堰,堰高P=P′=3.5m,堰顶水头H限制为0.85m,通过3堰顶的流量Q=20m/s,求堰宽b及不发生淹没出流的下游最大水深。7.5直角三角形薄壁堰,堰前水头H=0.2m,求通过此堰的流量。若流量增加一倍,问水头变化如何?7.6图示潜水坝,厚度δ=2m,坝高P=P′=1m,上游水位高出坝顶0.6m,下游水位高出坝顶0.1m,求通过坝顶的单宽流量。7.7设有一取水闸,堰坎系矩形进口宽顶堰,坎高P=P′=1m,堰前水头H=2m,堰下游水深h下=1.0m,题7.6图堰宽b=2m,引水渠宽B=3m。求取水闸通过的流量。7.8在上题中,如果下游水深h下=2.8m,其他条件均不改变,问通过流量是多少,设v0可忽略不计。37.9有一曲线形滚水坝(流线型墩),已知堰上水头H=2.4m,流量Q=80m/s,自由出流,流量系数m=0.48,试求该滚水坝的宽度(设行近流速水头可忽略不计)。·204·
第8章地下水动力学基础8.1概述液体在孔隙介质中的流动称为渗流。水在土壤空隙和岩石裂缝中的流动,是渗流的一个重要部分,又称为地下水运动。水在土壤中的存在状态有几种不同的类型,以水蒸气的形式散逸于土壤空隙中的水称为气态水;由于分子力的作用而聚集于土壤颗粒周围,其厚度小于最小分子层厚度的水称为吸着水;厚度在分子作用半径以内的水层称为薄膜水;由于表面张力作用而聚集于土壤颗粒周围的水称为毛细水;当孔隙介质中含水量甚大,受重力作用而运动的水称为重力水。地下水动力学研究的主要对象是重力水的运动。渗流运动的特性与孔隙介质的粒径、级配、均匀性、排列情况以及孔隙的大小、形状及孔隙系数等因素密切相关。从渗流的角度可将土壤分为均质土壤与非均质土壤。均质土壤是指其渗透性质不随空间位置而变化的土壤,否则为非均质土壤。均质土壤又可分为各向同性的和非各向同性的土壤,均质各向同性的土壤是指其渗透性质与渗流方向无关的土壤,例如,均质砂土就是均质各向同性土壤,而黄土和各个方向上有不同裂缝的岩石就是均质非各向同性土壤。本章只研究均质各向同性土壤中的重力水的恒定流。在土木建筑工程中有许多问题都涉及到渗流运动。例如地下水水源利用,它涉及到水井和集水廊道等集水建筑物的设计、产水量的计算等问题;堰、坝、渠道侧坡的修建涉及构筑物的稳定性问题及渗漏损失问题,在土木建筑施工中涉及基坑排水、路基排水等问题。本章所讨论的内容有:渗流模型和渗流基本定律,地下水的均匀流与非均匀流,集水廊道和井的水力计算。8.2渗流基本定律(1)渗流模型天然土壤中的颗粒形状及粒径大小各不相同,颗粒间的孔隙形状、大小及分布无一定规则。水在孔隙中的渗流运动是很复杂的,按实际情况进行分析将十分困难。因此,在研究渗流运动时,人们将孔隙介质所占据的空间模型化,认为该空间内没有土壤的颗粒(骨架)存在,只·205·
水力学有水充满全部空间,并沿主流方向作为连续介质而运动,这个空间中所通过的流量、断面上的压力以及流动阻力(水头损失)均与实际渗流相等,这样的空间称为渗流理论的简化模型,或简称渗流模型。设ω是渗流模型的过水断面面积,Q为通过该过水断面的流量,则定义:v=Q/ω为渗流的断面平均流速。渗流模型的过水断面面积ω不等于真实渗流的过水断面面积ω′(ω′是孔隙介质断面上的孔隙面积),ω′<ω。因此,上述定义的渗流断面平均流速v的值比真实的渗流平均流速v′小。设n=ω′/ω并称为土壤的孔隙率,则v=nv′。各种土壤的孔隙率大致如表8.1所示。表8.1土壤种类粘土粉砂中粗混合砂均匀砂孔隙率0.45~0.550.40~0.500.35~0.400.30~0.40土壤种类细、中混合沙砾石砾石和砂砂岩孔隙率0.30~0.350.30~0.400.20~0.350.10~0.20(2)达西定律液体在孔隙介质中流动时,由于液体粘滞性的作用,必然有能量损失。早在1852~1855年,法国学者达西(H·Darcy)对砂质土壤进行了大量渗流实验研究,总结出渗流水头损失与渗流速度之间的基本关系,即渗流达西定律。达西实验装置如图8.1所示,一上端开口的直立圆筒,内装颗粒均匀的砂土,上部由供水管A供水,并用溢流管B以恒定水位,渗透过砂体的水通过底部滤水网C流入容器D,并由此测定渗透流量。在筒壁上接通两测压管相距l,以测量1-1和2-2断面上的渗透压强。由于达西实验中的渗流流速很小,渗流流态为层流,所以可以忽略流速水头,因此1-1和2-2断面的测压管水头差ΔΗ就是渗流在l长度的水头损失hω,从而水力坡度J为:J=hω/l=(h1-h2)/l(8.1)达西以不同尺寸的圆筒和不同类型的土壤进行了图8.1大量的实验,观测,发现在不同尺寸的圆筒和不同类型的土壤渗流中所通过的渗流流量Q与圆筒的横断面积ω和水力坡度J成正比,并与土壤的渗透性质有关。可以表示为:Q=kwJ(8.2)Q或v==kJ(8.3)w式中:k为反映土壤渗透性质的系数,称为渗流系数,其单位为速度的单位,v为渗流流速。式(8.3)称为渗流达西定律的表达式。渗流达西定律描述了当渗流流速很小时,渗流能量损失与渗流流速之间的基本关系,揭示了渗流层流的基本规律:渗流层流的断面平均流速与水·206·
第8章地下水动力学基础力坡度的一次方成正比。(3)渗流系数渗流系数是反映孔隙介质渗透特性综合指标的重要参数。渗流系数的大小主要取决于土壤颗粒的形状、大小、不均匀系数及水温等等。要精确测定渗流系数的数值较为困难,通常采用经验公式法、实验室测定法、现场观测法等多种方法测算渗流系数的概值。本书仅大概介绍实验室测定法和现场测定法。1)实验室测定法实验室测定法通常使用类似于达西渗流实验所采用的装置在实验室进行实验,测出Q,h1及h2,采用式(8.2)计算渗流系数。此法简便易测,若选取的土壤是实际的未扰动土壤,并有足够数量有代表性的土壤进行实验,其结果是可靠的。2)现场观测法在现场利用钻井或原有井作抽水或灌水试验,然后根据井的公式(在后面的小节中讨论)计算渗流系数k。这种方法是可靠的测定方法,且实用意义大,可以取得大面积平均渗流系数值,但经济耗费大。-1渗流系数k的量纲为〔LT〕,常用cm/s或m/d表示,其中d表示天。作近似计算时,可以采用表8.2给出的水在土壤中渗流系数的概值。表8.2渗透系数k土名m/dcm/s-6<0.005<6×10粘土-6-10.005~0.16×10~1×10-4-4轻亚粘土0.1~0.51×10~6×10-4-4黄土0.25~0.53×10~6×10-4-3粉砂0.5~1.06×10~1×10-3-3细砂1.0~5.01×10~6×10-3-2中砂5.0~20.06×10~2×10-2-2均质中砂35~504×10~6×10-2-2均质粗砂60~707×10~8×10-2-2圆砾50~1006×10~8×10-1-1卵石100~5001×10~6×10-1无填充物卵石500~10006×10~1×10-2-2稍有裂隙岩石20~602×10~7×10-2裂隙多的岩石>60>7×10(4)达西定律的适用范围和非线性渗流定律达西定律表明渗流的沿程水头损失与流速的一次方成正比,即水头损失与断面平均流速成线性关系。凡符合这种规律的渗流,称为层流渗流或线性渗流。达西渗流定律又称为线性·207·
水力学渗流定律。当渗流流速较大(如在重粗颗粒土壤中或堆石中的渗流)时,水头损失与流速之间不再呈线性关系,当流速达到一定数值后,水头损失与流速的平方成正比,这种渗流称为非线性渗流。土壤的渗透性质十分复杂,难以找到线性渗流与非线性渗流的确切的判别准则。有人建议直接引用土壤粒径,这种方法过于粗略,大多数人建议如同管渠流动一样采用雷诺数。这方面有多种研究成果。下面仅介绍两个实验研究成果。一种是直接采用雷诺数的通常表达式:vdRe=(8.4)ν式中v———渗流断面平均流速,以cm/s计;d———骨架或土壤的特征粒径,通常采用d10,即筛分时占10%的重量的土粒所通过的筛孔直径,以cm计;2ν———的运动粘滞系数,以cm/s计。一般可取Re≤1~10作为线性定律的上限值。另一种是考虑土壤孔隙率n的雷诺表达式:1vdRe=(8.5)0.75n+0.23ν当实际土壤的雷诺数Re<7~9,则为线性渗流。工程上所遇到的较多渗流问题属线性渗流,但渗水路堤,堆石坝等的渗流则不符合线性渗流定律。颗粒极细的粘土,能否运用渗流达西定律进行计算也尚待研究。1901年福希海梅(Forchheimer)提出渗流水头损失的一般表达式为:2J=au+bu(8.6)式中:a和b为待定系数,由实验测定。当b=0时,即为线性渗流定律,当渗流进入紊流阻尼平方区时,a=0,即水头损失与流速的平方成正比;若a和b都不等于零,则为一般的非线性渗流定律。一些实验结果表明,渗流紊流开始于Re=60~150,达西定律在Re≥1~10时已不适2用了。因此Re≈10~150间的层流区,也有bu项出现。本章仅限于讨论符合达西定律的渗流。8.3地下水的均匀流与非均匀流采用渗流模型研究,即认为地下水的运动是连续的,因此可以应用研究地表明渠水流的方法将渗流分为均匀渗流和非均匀渗流。服从达西定律的渗流具有某些地表明渠流所没有的特点。(1)恒定均匀流与渐变流流速分布1)均匀渗流在均匀渗流中,若视不透水层的顶坡为渠道的底坡,设其坡度为i,地下水水面线平行于2v渠底,服从线性律的渗流流速很小,因此可以略去流速水头不计,总水头线与测压管水头线2g重合。而地下水水面线就是测压管水头线。均匀渗流任一断面的测压管水头线坡度(或水力·208·
第8章地下水动力学基础坡度)都相同,因此有:J=Jp=i=常数(8.7)即各个断面水力坡度都等于底坡。均匀渗流每一过水断面上的压强分布都与静水压强分布相p同,即服从(z+)=常数,故在断面上各个点的水力坡度都相等。根据达西定律渗流流场中γ某点的渗流流速为:u=kJ=ki(8.8)即均匀渗流区域中任一点的渗流流速都相等。因此均匀渗流过水断面上的流速分布图为矩形,且断面上流速分布图沿程不变,全渗流区各点渗流流速相等,如图8.2所示。图8.2图8.32)渐变渗流如图8.3所示渐变渗流,任取相距为dL的过水断面1-1和2-2。在渐变渗流断面上压强分布近似服从静水压强的分布规律,因此,1-1断面上各点的测压管水头都是H1,断面2-2上各点的测压管水头均为H2。断面1-1与2-2之间任一流线的水头损失相同,为:dH=H2-H1又,由于渐变流流线间的夹角小,流线的曲率小,流线族几乎为平行直线,因此可以认为1-1和2-2断面间的所有流线长度均近似为dL,故过水断面上各点的水力坡度相等为:dHJ=-=常数dL根据达西定律,过水断面上各点的渗流流速u都相等,断面平均流速就等于点的渗流流速,即:v=u=kJ(8.9)该式称为裘皮幼公式。该公式表明:渐变渗流过水断面上各点的渗流流速相等,因而断面平均渗流流速等于断面上任一点的渗流流速,不同的过水断面上的流速不相等,因为水力坡度是沿程变化的。(2)渐变渗流基本方程下面讨论渐变渗流的水力要素沿程变化的规律。如图8.4所示无压恒定渐变渗流。不透水层顶坡为渗流底坡i,任取一过水断面1-1,其含水层水深为h,测压管水头为H,渠底距基准面的高度为z,则有:H=h+z该断面上各点的水力坡度为:dHd(h+z)J=-=-dLdL·209·
水力学dzdhdh=--=i-dLdLdL由于水深h沿程变化,因此,渐变渗流的不同断面具有不同的测压管水头线坡度。将测压管水头线坡度J的表达式代入裘皮幼公式(8.8)中可得断面平均渗流流速和通过该断面的渐变渗流流量:dhv=k(i-)dL图8.4dhQ=kw(i-)(8.10)dL式(8.9)为适用于各种底坡渐变渗流的基本微分方程,也是分析和绘制渐变渗流水面线(称为浸润曲线)的理论依据。(3)渐变渗流浸润曲线的类型无压渗流的地下水水面称为浸润面,在流动的纵剖面上它是一条曲线,称为浸润曲线。对于渐变渗流,当流速水头可忽略时,总水头线与测压管水头线相重合,浸润曲线既是测压管水头线又是总水头线。由于渗流流动过程中必然存在水头损失,总水头线总是沿程下降的,因此,浸润曲线也只能是沿程下降的,不可能是水平线,亦不可能沿程上升。前面已经得出渐变渗流微分方程,积分之,可得浸润曲线方程。在分析地表明渠水面曲线时,正常水深和临界水深起着重要的作用,这里沿用地表明渠流的概念,讨论渗流问题时将均匀渗流的水深h0称为正常水深;将不透水层顶坡作为渗流底坡,按其坡度是否大于零,依次分为顺坡渗流:i>0;平坡渗流:i=0;逆坡渗流:i<0。dh均均流的水深沿程不变:=0,决定均匀渗流水力要素的基本方程为:dLQ=kω0i(8.11)式中:ω0为相应于正常水深h0的过水断面面积ω0=bh0。由于可以略去流速水头,断面单位2αv能(比能)Es=h+实际上就近似等于水深h,因此就不存在临界水深、缓流、急流等概念,在2g分析浸润线时只需用实际水深与特征水深(正常水深)相比较,故渐变渗流浸润曲线类型及其位置的分区比地表明渠水面曲线少,在三种坡度情况下总共只有四种浸润曲线类型。下面分别对三种不同的坡度进行讨论。1)顺坡(i>0)渗流:均匀渗流只发生于顺坡。将式(8.10)中的渗流流量用式(8.11)代替可得:kω0i=kω(i-dh),由此可得出顺坡浸润曲线的微分方程为:dLdhw0h0=i(1-)=i(1-)(8.12)dLwh记正常水深线为N-N线,N-N线把渗流划分为(a)和(b)两个区域,如图(8.5)所示。dh(a)区:h>h0,ω>ω0,因此>0,即:浸润线为壅水曲线。在浸润线上游端,当h→h0时,dLdhdh有=0,浸润线以N-N线为渐近线,在曲线下游端,当h→∞时,→i,曲线以水平线为渐dLdL·210·
第8章地下水动力学基础近线。dh(b)区:h0,为顺坡渗流,又因h2>h1,故浸润曲线为顺坡壅水曲线。由顺坡渗流基本方程(8.11)·212·
第8章地下水动力学基础图8.9dhh0=i(1-)(a)dLhh式中h0为相应于均匀渗流的水深。令:=ηh0则:dh=h0dη代入(a)式中得:ηdηi=dLη-1h0等式两端同时积分:η2η+1-1iη2-1iL∫dη=L,得:η2-η1+ln=η1η-1h0η1-1h0改为常用对数表达式:iLη2-1=η2-η1+2.30lg(b)h0η1-1h1h2(b)式中:η1=,η2=h0h0h11h21.9由题给条件i=0.02,L=180m,η1==,η2==,代入(b)式后并化简得:h0h0h0h01.9-h01h0lg=(0.02×180-1.9+1.0)=1.1741.0-h02.3采用试算法解得:h0=0.945m。故每米长渠道所渗出的流量为:3q=kh0i=0.005×0.945×100×0.02=0.00945cm/s·cm为绘制浸润曲线,将i=0.02,h0=0.945m,h1=1m代入(b)式可得:h2h2-0.945L=47.25(-1.058+2.30lg)(c)0.9450.0548分别假设h2=1.2m,1.4m,1.7m,1.9m,(因为渗流总长只有180m,才可以这样计算。否则应该分段求和,特别是当曲线变化较快时,分段应取很短。)代入(c)式相应的L为80.6m,117.7m,156.7m,180m,连接这些坐标点,即可绘出浸润曲线如图8.9所示。·213·
水力学8.4集水廊道的渗流计算集水廊道是建造于无压含水层中用以集取地下水源或降低地下水位的集水建筑物。设一位于不透水层上的矩形断面集水廊道如图8.10所示,底坡i=0。廊道中不断抽水时,地下水流向廊道,水面沿程下降,在其两侧形成对称于廊道轴线的浸润曲面。刚开始抽水时,这种渗流属于非恒定流。若含水层很大,廊道很长,渗流持续一段时间后,可近似地形成无压渐变渗流,廊道中保持某一恒定水深h,两侧浸润曲线的形状位置基本不变,在所有垂直于廊道轴线的剖面上,渗流情况相同,可作为平面渗流问题讨论。图8.10取廊道右侧单位长度研究,设q为集水廊道单位长度上自一侧渗入的单宽流量,由式(8.9):dhQ=kω(i-)dLqdhq得:dL=-hdh或:=-kdLkhdhdz建立坐标系xoz如图8.10所示,x坐标与流向相反,故:=-,代入上式得渐变渗流的基dLdx本微分方程为:dzq=kzdx将该式分离变量并积分,代入边界条件:x=0时,z=h1,可得到集水廊道浸润曲线方程:222qz-h1=x(8.14)k可见,浸润曲线是抛物型曲线,当x越大,地下水位的降落就越小。设在x=L处,地下水位降落趋近于零,z等于含水层厚度H,L为集水廊道的影响范围,将这一边界条件代入式(8.14)中可得集水廊道单位长度上每侧的产水量公式为:22k(H-h1)q=(8.15)2LH-h式中,k与地质条件有关,由抽水试验决定。若引入浸润曲线的平均坡度:J=这一概念,L则式(8.15)可以改写为:kq=(H+h)J(8.16)2式(8.16)可用以估算q,J的数值可根据土壤性质确定,由表8.3选取。·214·
第8章地下水动力学基础表8.3浸润曲线的平均坡度土壤类别J值粗砂及冰川沉积土0.003~0.005砂土0.005~0.015微弱粘性砂土0.03亚粘土0.05~0.10粘土0.158.5单井的水力计算井在给水工程中是吸取地下水的建筑物,应用很广。从井中抽水可使井附近的天然地下水位降落,可起到排水或降低地下水位的作用,也可向井中输水,使地下水水位提高。根据水文地质条件,可将井分为潜水井和承压井两种基本类型。潜水井指在具有自由水面的潜水层中开凿的井,又称为无压井,若井底达不透水层,则称为完全井;若井底未达到不透水层,则称为不完全井;承压井指在两个不透水层之间的含水层中凿井,含水层压强大于大气压强,承压井又称为自流井。当地下水开采量较大,补给来源不足或者需要精确测定水文地质参数时,应按非恒定流考虑。但在地下水来源充沛,开采量远小于天然补给量的情况下,经一段时间抽水后,可按恒定流分析井的渗流情况。严格地讲,井的渗流属于三维渗流,求解是非常复杂的,但若忽略运动要素沿z轴方向的变化,并采用轴对称假设,即可采用一维渐变渗流的裘皮幼公式进行分析。(1)潜水井(无压井)具有自由水面的地下水称为无压地下水或潜水。在潜水中修建的井称为潜水井或无压井。井的断面通常为圆形,水由井壁渗入井中。如图8.11所示完全潜水井,井底为不透水层,含水层厚度为H,井的半径为r0。抽水时地下水从四周径向对称流入井内,形成对井中心垂直轴线对称的漏斗形浸润曲面。设抽水量不变地连续抽水,且含水层体积很大,抽水过程中不致使含水层厚度H有所改变,则流向水井的地下水渗流为恒定渗流,浸润面的位置不变,井中水深h也保持不变。流向水井的过水断面是一系列圆柱面,各径向剖面的渗流状况相同。图8.11可以运用裘皮幼公式计算断面平均流速。距井轴为r的过水断面,其高度为z,面积为dzdz2πrz,此圆柱面上各点的水力坡度皆为,故断面平均流速为:v=k,由Q=ωv得渗流流drdr·215·
水力学zrdzQdr量:Q=2πrzk,将该式分离变量,并从r0至r及h至z取定积分,即∫:zdz=∫由drhr02πkr此可求得完全潜水井的浸润曲线方程:22Qr0.732Qrz-h=ln=lg(8.17)πkr0kr0设在半径r=R的过水断面上,潜水水深z=H,即该处天然地下水位已不受井抽水影响,则距离R称为井的影响半径。将r=R时,z=H这一边界条件代入式(8.17)中可得:22k(H-h)Q=1.366(8.18)Rlgr0式(8.18)即为完全潜水井的产水量公式,称为裘皮幼产水量公式。对于一定的产水量Q,地下水面相应的最大降落深度为:S=H-h称为水位降深。从而有:22SH-h=(H+h)(H-h)=2H(1-)S2HS当H远大于S时,n1,可略去之,则(8.18)式可以简化为:2HkHSQ=2.73(8.19)Rlgr0式(8.19)表明:产水量Q与k、H及S成正比,而Q是随lgR而变化,故R值对Q的影响较小。一般,影响半径由抽水试验测定。在估算中,R值可按经验酌情选用:粗粒土壤R=700~1000m;中粗粒土壤R=250~700m;细粒土壤R=100~200m。也可以采用以下经验公式估算R值:R=3000Sk(8.20)式中:S、R以m计,k以m/s计。不完全井的产水量不仅来自井壁四周,而且还来自井底,其产水量公式一般由经验公式测定,此处不详述。例8.4一完全潜水井,井的半径r0=0.5m,含水层厚度H=8m,土壤的渗流系数k=0.0015m/s,抽水时井中水深h=5m,试估计井的产水量。解水位降深:S=H-h=8-5=3m由经验公式估算井的影响半径R的值:R=3000Sk=3000×3×0.0015=342.6m取影响半径R=350m,可求得井的产水量为:2222k(H-h)0.0015×(8-5)3Q=1.366=1.366×=0.028m/sR350lglgr00.5(2)自流井(承压井)如图8.12所示,含水层位于两不透水层之间,含水层的压强大于大气压强,这样的含水层称为承压含水层。凿井穿过上面的不透水层,从含水层中取水,这样的井称为自流井或承压·216·
第8章地下水动力学基础井。若井底直达下部不透水层的表面,则为完全自流井,图8.12所示为完全自流井。在本书中仅讨论一种最为简单的情况,即如图8.12所示,下面的不透水层水平,并且含水层厚度t为定值的完全自流井。未抽水时,在含水层压力作用下水深上升到H,H即为含水层的天然总水头,井中水面的延长面为地下水天然水头面,它高于t,若含水层压力较大,还有可能高出地面。若从井中抽水,井中水深由H降至h。若在上部不透水层钻若干小井作为测压管用,则可观测到抽水井外的测压管水头线沿渗流方向沿程下降,当抽水一段时间图8.12后,可近似地形成一个对称于井轴的漏斗形水头降落曲面。如图8.12所示。承压井渗流的过水断面为一系列高度为t的圆柱面,各径向剖面的渗流情况相同,除井周附近的区域外,测压管水头线的曲率很小,恒定抽水时,可作为恒定渐变渗流分析。dz对于半径为r的圆柱面过水断面,断面平均流速为:v=k,由Q=vω可得平坡渗流微分dr方程(或渗流流量公式):dzQ=wv=2πrtkdr式中z为半径等于r的过水断面的测压管水头。将上式分离变量,并从r0到r,h到z积分得:Qrz-h=0.366lg(8.21)ktr0式(8.21)即为自流井的测压管水头线方程。若同样引入影响半径的概念,当z=H时,r=R,(当r>R以后,测压管水头高度保持为H),可得:kt(H-h)ktsQ=2.73=2.73(8.22a)RRlglgr0r0RQlgr0或S=(8.22b)2.73kt影响半径R也可按照完全潜水井的方法确定。例8.5如图8.13所示,一完全自流井的半径r0=0.1m,含水层厚度t=5m,在离井中心r1=10m处钻一观测孔。在未抽水前,测得地下水的天然总水头3H=12m。现抽水流量Q=30m/h,井中水位降深S0=2m,观测孔中水位降深S1=1m,试求含水层的渗流系数k及影响半径R。解由题给条件知:S1=H-h1,所以h1=H-S1=12-1=11m,S0=H-h0,所以h0=10m,由式(8.20)图8.13·217·
水力学Qrz-h=0.366lgktr0将r=r1,z=h1,h=h0各条件代入得:3010h1-h0=0.366×lg=13600×k×50.1于是可解得渗流系数:k=0.00122m/skt(H-h)再由式(8.21):Q=2.73中解得影响半径R:Rlgr02.73kt(H-h)2.73×0.00122×5×(12-10)×3600lgR=+lgr0=+lg0.1Q30=3.997-1≈33所以影响半径:R=10=1000m。(3)大口井与基坑排水大口井是用以集取浅层地下水的一种井,井径较大,大致为2~10m或者更大,这种井类似于一个很大的坑。基坑排水与大口井集水相似,其计算方法基本相同。大口井可以是完全井,也可以是不完全井,但一般都是不完全井。井壁可以是透水的,也可以是不透水的。井底进水量往往很大,常为总产水量的主要部分。对于井壁与井底同时进水的大口井,其分析十分复杂。本章讨论假设井壁不透水,而只有井底进水的大口井的渗流。设有一大口井,井壁四周为不透水层,井底为半球形,紧接下层深度为无穷大的含水层。供水是由井底的渗流提供的。如图8.14所示。图8.14图8.152半球底大口井的渗流流线是径向的,过水断面为与井底同心的半球面,Q=ωv=2πr·dzk,分离变量积分:drrzdr∫Q2=2π∫kdzrrH-S0当r=R时,Z=H,且Rmr0,故得:Q=2πkr0S(8.23)·218·
第8章地下水动力学基础式(8.23)为半球底大口井的产水量公式。对于平底的大口井,其过水断面近似为椭圆,流线是双曲线,如图8.15所示。其产水量公式为:Q=4kr0S(8.24)式(8.23)和式(8.24)两式的计算结果相差甚大。当含水层比井的半径大8~10倍时,采用(8.23)式为好。8.6井群的水力计算井群是指多个井同时工作,井与井之间的距离小于一个井的影响半径的多个井的组合,如图8.16所示。抽水时,各井之间相互影响,渗流区地下水流比较复杂,其浸润面的形状也十分复杂,因此,井群的水力计算也比单井复杂得多。(1)完全潜水井井群的浸润曲面方程先讨论完全井渗流运动的连续性微分方程:将xoy坐标平面建立在渗流流场的不透水层上如图8.17所示,则浸润曲面的方程可表示为:z=f(x,y)。在渗流流场中自不透水层至浸润面取一底面积为dxdy的微小柱体,其高度为z,渗流通过该微小柱体的质量守恒。从adeh(ρQx)面流入柱体的质量流量为:ρQx=ρωxvx,从bcgf面流出柱体的质量流量为:ρQx+dx从x侧面abcd流入柱体的质量流量为:ρQy=ρωyvy图8.16图8.17(ρQy)从efgh面流出柱体的质量流量为:ρQy+dyy根据质量守恒定律得:(ρQx)(ρQy)(ρQx+dx-ρQx)+(ρQy+dy-ρQy)=0xy于是可得连续性微分方程为:(ρQx)(ρQy)dx+dy=0(8.25)xy此式对于完全潜水井及完全承压井均适用。·219·
水力学zQx=zdy·k·x对于完全潜水井,根据达西渗流定律有:zQy=zdx·ky将其代入式(8.25)中并考虑到在不可压缩流体条件下,ρ=常数,于是可得:2222(z)(z)2+2=0(8.26)xy2式(8.26)为完全潜水井浸润面z所应满足的微分方程。由该方程可知,式中z是满足线性方2程(即拉普拉斯方程)的函数,因此,函数f(z)=z可以叠加。即:当井群的所有井共同工作22时,所形成的z函数为井群中各井(记为第i个井)单独工作时的zi之和,即:2222z=z1+z2+⋯=∑Zi(8.27)i设井群中的第i个井的抽水量为Qi,井中水深为hi,井的半径为roi,由式(8.17)知:20.732Qiri2zi=lg+hi(8.28)kroi将上式代入式(8.27)中,得:nn220.732Qiri2z=∑zi=∑(lg+hi)(8.29)i=1i=1kroiQ0若各井产水量相同,即:Q1=Q2=⋯⋯=Qn=,Q0为n个井的总产水量,则:n20.732Q02z=〔lg(r1r2⋯⋯rn)-lg(r01r02⋯⋯ron)〕+∑hi(8.30)k设井群的影响半径为R,在影响半径上取一点A,A点距各井很远,即:r1≈r2≈⋯⋯rn=R,而z=H,代入式(8.30)中得:20.732Q2H=〔nlgR-lg(r01r02⋯⋯r0n〕+∑hi(8.31)k将式(8.30)与式(8.31)相减得:220.732Q01z=H-〔lgR-lg(r1r2⋯⋯rn)〕(8.32)kn式(8.32)为完全潜水井井群的浸润曲面方程。式中影响半径R可采用下式计算:R=575SHk(8.33)式中S———井群中心的水位降深,以m计;H———含水层厚度,以m计。(2)完全潜水井群产水量公式由式(8.32)完全潜水井井群的浸润曲面方程可以解得当各井产水量相等时完全潜水井群产水量公式为:221.366k(H-z)Q0=(8.34)1〔lgR-lg(r1r2⋯⋯rn)〕n(3)自流井井群的测压管水头方程用分析完全潜水井井群的方法去分析完全自流井井群,对于承压含水层的厚度t为常数的情况可得:·220·
第8章地下水动力学基础22zz2+2=0xy即:完全承压井的测压管水头函数z(x,y)满足拉普拉斯方程,具有可叠加性。于是,完全承压井井群的测压管水头方程为:n0.366Q01z=∑zi=H-〔lgR-lg(r1r2⋯⋯rn)〕(8.35)i=1ktn井群的产水量为:2.732kt(H-z)Q0=(8.36)1〔lgR-lgr1r2⋯⋯rn〕n因为第i个单自流井测压管水头方程为:0.366Qirizi-hi=lg(a)ktroi当z=H时,r=R,代入(a)式中得:0.366QiRH-hi=lg(b)ktroi将(a)、(b)两式相减可得单井的水头降深:0.366QiRSi=H-zi=lgktri当井群中各井抽水量相等时,总产水量Q0=nQi,则由(8.34)式得井群的水头降落:nn0.36610.366QiRS=H-z=nQi[lgR-lg(r1r2⋯rn)]=∑lg=∑Si(8.37)ktni=1ktrii=1式(8.37)说明自流井井群同时均匀地抽水时,任一点A的水头降落等于各井单独抽水时A点的水头降落之和。这就是自流井井群的水头降落叠加原理。例8.6一如图8.18所示的无压完全井井群,用以降低基坑中的地下水位。已知a=50m,b=20m,各井的抽水量相等,其总的抽水流量Q0=6l/s,各井的半径均为r0=0.2m,含水层厚度H=10m,土壤为粗砂,其渗流系数k=0.01cm/s,取影响半径R=800m,试求:B点和G点的地下水位降低值SB和SG。解总抽水量:3Q0=6l/s=0.006m/s对于G点:a2b2rA=rC=rF=rD=()+()=26.93m22rB=rE=10m对于B点:rA=rC=25m,rB=0.2m,图8.18rF=rD=30.02m,rE=20m由完全潜水井群的浸润线方程:22Q01z=H-0.732〔lgR-lg(r1r2r3r4r5r6)〕kn可得G点:·221·
水力学220.006142z=10-0.732〔lg800-lg(26.93·10)〕0.000160.006=100-0.732×〔2.90-1.28〕=100-71.15=28.850.0001z=5.37m因此G点水位降SB=H-z=10-5.37=4.63m对于B点:20.006122z2=10-0.732〔lg800-lg(25·0.2·30.02·20)〕0.000162=10-43.92〔2.9-1.06〕=19.13z=4.37因此B点的水位降深:SB=H-z=10-4.37=5.63m思考题8.1什么是渗流模型?为什么要引入这一概念?8.2渗流流速指的是什么?它与真实渗流的流速有什么区别?8.3试比较达西渗流定律的表达式与裘皮幼公式有何异同?各自的应用条件是什么?8.4影响渗流系数的因素有哪些?8.5地表上棱柱形渠道的水面曲线有12条,为什么渐变渗流的浸润曲线只有4条?它们都是些什么类型?8.6何为潜水层?何为自流层?8.7什么是完全井与不完全井?8.8影响潜水井渗流流量的主要因素有哪些?影响自流井渗流流量的主要因素有哪些?8.9什么是大口井?8.10什么是井的影响半径?自流井有没有影响的半径?8.11什么是井群?8.12如何求完全承压井井群的水头面方程?习题8.1在实验室中,根据达西渗流定律测定某土壤的渗流系数时,将土壤装在直径D=20cm的圆筒中,在40cm的水头差作用下,经过一昼夜测得渗透水量为15l,两测压管间的距离为30cm,试求:该土壤的渗流系数k。8.2圆柱形滤水器,其直径d=1.2m,滤层高1.2m,渗流系数k=0.01cm/s,求H=0.6m时的渗流流量Q。8.3已知渐变流浸润曲线在某一过水断面上的坡度为0.005,渗流系数为0.004cm/s,·222·
第8章地下水动力学基础求过水断面上的点渗流流速及断面平均渗流流速。8.4厚度t=15m的承压含水层,有两个观测井,距离l=200m,测得观测井1中水位为64.22m,观测井2中的水位为63.44m,如图所示。含水层由粗砂组成,已知渗流系数k=45m/d,试求该含水层单位宽度(每米)的渗流量q。8.5如图所示两水池A、B,中间为一水平不透水层上的砂壤土山丘,已知A池水位为15m,B池水位为10m,不透水层高程为5m,砂壤土的渗流系数k=0.0005cm/s,两水池间的距离l=题8.2500m,试求:1)单宽流量q;2)浸润曲线坐标y=f(x)。题8.4图题8.5图8.6为了查明地下水储藏情况,在含水层土壤中相距S=500m处打两钻孔1和2,测得两个钻孔中水深分别为h1=3m,h2=2m,不透水层的底坡i=0.0025,渗流系数k=0.005cm/s,试求:1)渗流单宽流量q;2)两钻孔中间断面C-C处的地下水深度hc。8.7一水平不透水层上的渗流层,宽800m,渗流系数为0.003m/s,在沿渗流方向相距1000m的两个观测井中,分别测得水深为8m及6m,求渗流流量Q。题8.6图8.8如图所示不透水层上的排水廊道,已知:垂直于低面方向长100m,廊道中水深h0=2m,天然含水层水深H=4m,土壤的渗流系数k=0.001cm/s,廊道的影响半径R=200m。试求:1)廊道的排水流量Q;2)距廊道100m处C点的地下水深hc。8.9在公路沿线建造一条排水明沟(如图)以降低地下水位。含水层厚度H为1.2m,土壤渗流系数k为0.012cm/s,浸润曲线的平均坡度J为0.03,沟长L为100m,试求从两侧流向排水明沟的流量,并绘制浸润曲线。8.10某工地以潜水为给水水源,钻探测知含水层为沙夹卵石层,含水层厚度H=6m,渗流系数k=0.0012m/s,现打一完全井,井的半径r0=0.15m,影响半径R=300m,求井中水位·223·
水力学题8.8图题8.9图降深S=3m时的产水量。8.11完全潜水井,直径为80cm,含水层厚为6m,渗流系数k=3.6cm/min,井中水位降落为2m,水位恒定,要求抽水量。38.12在潜水井中进行抽水试验。测得恒定的产水量为92m/h,在距井轴8m处设观测井,测得水位降深为58cm。距井轴25m处设观测井,测得水位降深为46cm,并测得未抽水前潜水层厚度为12.60m,试确定含水层渗流系数k。8.13承压井的现场测定,井半径为20cm,距井30m和10m处各设有一个观测孔,孔中水位降落分别为20cm和42cm,含水层厚度为6m,井中水位为34.8m,Q=24m/h,H=9m,要求渗流系数k和井壁内外水位差。8.14完全自流井中h3时,则有(n-3)个指数需用其他指数值的函数来表示。将所求得的各αi值,代回式(9.3)即可得诸因素间的函数关系式。例9.1由第四章讨论知,流动有两种形态:层流和紊流,流态相互转变时的流速称临界流速。实验表明,恒定有压管流的下临界流速υc与管径d,液体密度ρ,动力粘性系数μ有关。试用量纲分析法求出它们的函数关系。解按瑞利法解本题。首先将关系式写成指数关系:αααυc=kd1ρ2μ3其中k为无量纲量。-1-3-1-1各量的量纲分别为:dimυc=LT,dimd=L,dimρ=ML,dimμ=MLT。将上式指数方程写成量纲方程:-1α-3α-1-1αLT=(L)1(ML)2(MLT)3则有·228·
第9章量纲分析和相似原理L:1=α1-3α2-α3T:-1=-α3M:0=α2+α3解得α3=1,α2=-1,α1=-1。将各指数代入原式,得μνυc=k=kρdd若将上式化为无量纲形式,有υcdk=ν式中无量纲数k称临界雷诺数,以Rec表示,即υcdRec=ν根据雷诺实验,该值在恒定有压管流动中为2000,可以用来判别层流和紊流。例9.2根据观察,实验与理论分析,认为圆管流动中管壁切应力τ0与液体的密度ρ,动力粘性系数μ,断面平均流速υ,管径d及管壁粗糙凸出高度Δ有关。试用瑞利法求τ0的表达式。解根据上述影响因素,将关系式写成指数关系:ααααατ0=kρ1μ2υ3d4Δ5其中k为无量纲数。写出量纲关系式为-1-2-3α-1-1α-1αααMLT=(ML)1(MLT)2(LT)3L4L5则有:M:1=α1+α2L:-1=-3α1-α2+α3+α4+α5T:-2=-α2-α3这是5个未知数3个方程的方程组,以α1,α5为待定指数,分别求出α2,α3,α4为α2=1-α1,α3=1+α1,α4=-1+α1-α5因此α1-α1+α-1+α-αατ0=kρ1μ1υ1d15Δ5ρυ·dαμΔα2=k()1()()5ρυμρυ·ddρυ·d其中k,α1,α5可取任何数值都不会影响上式的量纲和谐性。极易证明为无量纲量,水力μ学中称为雷诺数,记为Re,则5α-1Δα2Δ2τ10=k(Re)()ρυ=f(Re,)ρυddΔλλ2令f(Re,)=,式中λ为沿程阻力系数,由实验确定。所以τ0=ρυ。d88从以上二例可知,在独立影响因素不多余3时,用瑞利法很易求得表达某一物理过程的方程。反之,则出现待定指数,分析起来较为困难。此种情况下可采用π定理方法。(2)π定理π定理是量纲分析更为普遍的定理。π定理指出·229·
水力学任何一个物理过程,如果包含有n个物理量x1,x2,⋯,xn,则这个物理过程可用一完整的函数关系表示为f(x1,x2,⋯,xn)=0(9.5)其中m个物理量在量纲上是互相独立的,其余(n-m)个物理量是非独立的,则此物理过程可用(n-m)个无量纲数π表示的函数关系来描述,即F(π1,π2,⋯,πn-m)=0(9.6)因无量纲数是以符号π表示,所以称之为π定理。π定理在1915年由布金汉首先提出,又称为布金汉π定理。现在介绍应用π定理建立一物理过程的物理方程的步骤:①确定影响某一物理过程的物理量,并写成式(9.5)。②从n个物理量中选取m个在量纲上互相独立的物理量,称为基本物理量。对于不可压缩流体的运动,m一般为3。在实践中,常分别选几何学的量(如管径d、水头H等),运动学的量(如速度v、加速度a等)和动力学的量(如密度ρ、动力粘性系数μ等)各一个,作为独立物理量。③三个基本物理量依次与其余物理量组合成一个无量纲的π数,这样一共可写出(n-3)个π项:αβγπ1=x11x21x31x4αβγπ2222=x1x2x3x5⋯⋯⋯⋯⋯αβγπn-3n-3n-3n-3=x1x2x3xn式中αi,βi,γi为各π项的待定指数。④每个π项是无量纲数,可根据量纲和谐原理,求出各π项的指数αi,βi,γi。⑤写出描述此物理过程的无量纲关系式F(π1,π2,⋯,πn-3)=0例9.3用π定理求解例9.2中τ0的表达式。解①拟定函数关系式f(τ0,ρ,μ,υ,d,Δ)=0②从各物理量中选取d(几何量),υ(运动量),ρ(动力量)为基本物理量。③写出n-3=6-3=3个无量纲π项:αβγπ1=d1υ1ρ1τ0(1)αβγπ2222=dυρμ(2)αβγπ3=d3υ3ρ3Δ(3)④根据量纲和谐原理,各π项的指数分别确定如下:对(1)式,其量纲式为α-1β-3γ-1-2dimπ1111=L(LT)(ML)(MLT)则:L:0=α1+β1-3γ1-1T:0=-β1-2M:0=γ1+1联立以上三式求解得α1=0,β1=-2,γ1=-1。则可得到-1-2π1=τ0ρυ·230·
第9章量纲分析和相似原理同理,求得-1-1-1-1π2=μdυρ=(Re)-1π3=Δd⑤将各π项代入式(9.6)得无量纲数方程为τ01ΔF(2,,)=0ρυRedτ0Δ或写成2=f1(Re,)ρυdΔλ令f1(Re,)=,则d8λ2τ0=ρυ8量纲分析法在水力学研究中很有用处。但量纲分析毕竟是一种数学分析方法,有一定的局限性。它要求正确选择与物理过程有关的影响因素,正确选择基本物理量应具有独立的量纲等,最后在确定函数关系式的具体形式时,还必须依靠理论分析和实验的成果。9.3流动相似的概念模型实验的结果能够应用于原型的重要条件是,模型与原型应保证流动相似,即两个流动的对应点上所有表征流动状况的同名物理量之间(如速度、压强、各种力)相互平行,并都维持各自的固定比例关系,则这两个流动是相似的。这要求两个流动应满足几何相似、运动相似和动力相似以及初始条件和边界条件的相似。在下面的讨论中,原型(prototype)中的物理量标以下标p,模型(model)中的物理量标以下标m。(1)几何相似几何相似是指原型和模型两个流场的几何形状相似。要求两个流场中所有相应长度都维持一定的比例关系,对应角度相等。即lpλl=(9.7)lmθp=θm(9.8)式中lp和θp代表原型某一部位的长度和某两线段的夹角;lm和θm代表模型相应部位的长度和相应两线段的夹角;而λl称为长度比尺。几何相似的结果必然使任何两个相应的面积ω和体积V也都维持一定的比例关系,即ωp2λω==λl(9.9)ωmVp3λV==λl(9.10)Vm由上可知,长度比尺λl是几何相似的基本比尺,其他比尺均可通过长度比尺λl来表示。λl视实验场地与实验要求不同而取不同的值,在水工模型实验中,通常λl在10与100范围内取值。·231·
水力学(2)运动相似运动相似是指原型和模型两个流场对应点上同名的运动学量成比例,指原型和模型两个流动的空间对应点(包括边界上的点)处,质点流动在相应瞬间的速度和加速度分别方向相同,大小维持一定的比例。tpλt=(9.11)tmυpupλlλυ=λu===(9.12)υmumλtapλlλa==2(9.13)amλt式中λt称为时间比尺,λv或λu称为流速比尺,λa称为加速度比尺。(3)动力相似动力相似是指原型和模型两个流场对应点上各同名作用力方向分别相互平行,大小维持一定的比例关系。所谓同名作用力是指具有同一物理性质的力。作用在液体上的力通常有重I力G,粘性力Fμ,压力P,弹性力FE,表面张力Fσ,惯性力F等。两个流动动力相似,力的比尺λF为IGpFμpPpFEpFσpFpλF======I(9.14)GmFμmPmFEmFσmFm(4)边界条件相似与初始条件相似边界条件相似是指原型和模型两个流场的边界性质相同。边界条件可分为几何的、运动的和动力的几个方面,例如,原型中为固体边壁,模型中也应为固体边壁,原型中液面为自由液面,模型中液面也应为自由液面。对于非恒定流动,原型和模型两个流场还应满足初始条件相似。但在恒定流中,初始条件则失去实际意义。初始条件和边界条件的相似是保证两流动相似的必要条件。9.4相似准则几何相似是运动相似和动力相似的前提和依据,动力相似是决定两个流动运动相似的主导因素,运动相似是几何相似和动力相似的表现,在几何相似的前提下,要保证流动相似,则要实现动力相似。作用于液体上的重力、粘性力、压力、弹性力、表面张力等总是企图改变液体的流动状态,而惯性力却企图维持液体原有运动状态。液体运动的变化和发展就是惯性力和其他各种作用力相互作用的结果。I设作用在液体上外力的合力为F,液体质量为m,产生的加速度为a,惯性力为F=-ma,则由动力相似有:I22FpFpmpapρpVpap322ρplpυpλF==I===λρλlλa=λρλlλυ=22(9.15)FmFmmmamρmVmamρmlmυm也可写为·232·
第9章量纲分析和相似原理FpFm22=22(9.16)ρplpυpρmlmυmF式中22为一无量纲数,称为牛顿数。以Ne表示。即ρlυFNe=22(9.17)ρlυ式(9.16)可用牛顿数表示为(Ne)p=(Ne)m(9.18)上式表明,两个流动的动力相似,必须两个流动各相应处的牛顿数相等。这是流动相似的重要标志和判据,也称为牛顿相似准则。式(9.16)亦可写成λF22=1λρλlλυλFλFλl/λυλFλt式中22=3=λρλlλυλρλlλυλmλυ即λFλt=1(9.19)λmλυλFλt式中λm是原型和模型流动的质量比尺。称为相似判据。它和牛顿数一样,是用来判别相λmλυ似现象的重要标志。所以得出结论:对动力相似的流动,其相似判据为1,或相似流动的牛顿数相等。要使流动完全满足牛顿相似准则———两个流动的牛顿数相等,就要求作用在相应点上各种同名力均有同一的力的比尺。但由于各种力的性质不同,影响它们的物理因素不同,实际上很难做到这一点。在某一具体流动中占主导地位的力往往只有一种,因此模型实验中只要让这种力满足相似条件即可。这种相似虽是近似的,但实践证明,由此得到的结果也是令人满意的。下面分别介绍只考虑一种主要作用力的相似准则。(1)重力相似准则重力是液流现象中最常遇到的一种作用力,凡有自由液面并且允许液面上下自由变动的各种流动,如明渠流动,堰坝溢流,闸孔出流等都是重力起主要作用的流动。根据式(9.15),力的比尺λF可写为322GpρplpgpρplpυpλF==3=22Gmρmlmgmρmlmυm22υpυm化简后得=gplpgmlm开方后有υpυm=(9.20)gplpgmlm式中υ/gl为一无量纲数,称佛汝德(Froude)数,以Fr表示。即υFr=(9.21)gl·233·
水力学它表征了惯性力与重力的对比关系。式(9.20)用佛汝德数表示为(Fr)p=(Fr)m(9.22)上式表明,若作用在液体上的外力主要是重力,要使两个流动动力相似,必须两个流动相应处的佛汝德数相等,称为重力相似准则或称为佛汝德相似准则。式(9.20)也可写成比尺形式有2λυ=1(9.23)λgλl(2)粘滞力相似准则管道中的有压流动、污水处理池中的流动及潜体绕流问题等,主要是受水流阻力作用。重力对这种流动的机理不起作用,阻力又主要与粘滞力的作用有关,所以这类流动的相似就要求粘滞力作用相似。根据式(9.15),力的比尺λF可写为22FμpμplpυpρplpυpλF===22Fμmμmlmυmρmlmυm化简后得υplpυmlm=(9.24)νpνm式中υ·l/ν为一无量纲数,是前已介绍过的雷诺数Re,式中l为断面的特征几何尺寸,可为管径d或水力半径R。雷诺数Re表征了惯性力与粘滞力之间的对比关系。式(9.24)可用雷诺数表示为(Re)p=(Re)m(9.25)上式表明:若作用在液体上的外力主要是粘滞力,要使两个流动动力相似,必须两个流动相应处的雷诺数相等,称为粘滞力相似准则或称为雷诺相似准则。式(9.24)也可写成比尺形式有λυλl=1(9.26)λν(3)压力相似准则若作用于液体上的外力主要是动水压力P时,根据式(9.15),力的比尺λF可写为222PppplpρplpυpλF==2=22Pmpmlmρmlmυm化简后得:pppm2=2(9.27)ρpυpρmυm2式中p/ρυ为一无量纲数,称欧拉(Euler)数,以Eu表示,即pEu=2(9.28)ρυ它表征了惯性力与压力的对比关系。式(9.27)可用欧拉数表示为(Eu)p=(Eu)m(9.29)上式表明:若作用在液体上的外力主要是动水压力时,要使两个流动动力相似,必须两个流动相应处的欧拉数相等,称为压力相似准则或称为欧拉相似准则。式(9.27)也可写成比尺关系有λp2=1(9.30)λρλυ在不可压缩液体的有压流动中,起作用的是压强差Δp,而不是压强的绝对值。所以欧拉数也·234·
第9章量纲分析和相似原理可表示为ΔpEu=2(9.31)ρυ欧拉准则不是独立的准则,当雷诺准则和佛汝德准则得到满足时欧拉准则自动满足。作用在相似液流上的同名作用力不止以上三类,因此,还有另外一些需要满足的准则,如弹性力相似准则、表面张力相似准则等。限于篇幅,不再赘述。9.5模型实验模型实验依据相似原理,制成与原型相似的小尺度模型进行实验研究,以实验结果判断或预测原型的流动现象。进行模型实验首先需解决模型律的选择及模型设计两个问题。(1)模型律的选择模型律的选择应依据流动相似准则。为尽可能使模型与原型完全相似,除首先考虑几何相似以外,各独立的相似准则应同时满足,事实上,这是很困难甚至不可能的,只能对所研究的流动问题进行深入分析,找出影响该流动的主要作用力并选用相应的相似准则。例如,当粘滞力起主导作用时,则选用雷诺准则设计模型,称雷诺模型;当重力起主导作用时,则选用佛汝德准则设计模型,称佛汝德模型。下面分析这二种模型。1)雷诺模型在雷诺模型中,由于原型与模型的雷诺数相等,可根据式(9.26)来确定长度比尺与其他比尺的关系。流速比尺由下式确定:λυ=λν/λl(9.32)这就是说λυ取决于λν与λl之比,不能任意选择。例如当模型与原型中的液体相同,且温度也相同时,运动粘性系数比尺λν=1,则流速比尺与长度比尺之间为倒数关系,即λυ=1/λl(9.33)上式表明,雷诺模型尺度越小,模型中流速越快,即模型中流速将远大于原型中流速,这是雷诺模型的一个特点,雷诺模型中其他比尺亦可导出,例如流量比尺λQ2λQ=λυλω=λυλl=λl(9.34)时间比尺λt2λt=λl/λυ=λl(9.35)2)佛汝德模型在佛汝德模型中,由于原型和模型的佛汝德数相等,可根据式(9.23)来确定长度比尺λl与其他比尺的关系。当模型与原型流动均在地球上时,λg=1,所以流速比尺λυ由下式表示λυ=λl(9.36)佛汝德模型中其他比尺亦可导出由长度比尺λl表示,例如,流量比尺λQ与时间比尺λt分别为25/2λQ=λυλl=λl(9.37)λt=λl/λυ=λl(9.38)若要满足粘滞力与重力同时相似,即要保证模型与原型流动中的雷诺数和佛汝德数一一·235·
水力学对应相等。如果模型与原型采用同一种介质,由雷诺数相等条件,有λυ=1/λl由佛汝德数相等条件,有λυ=λl显然,只有λl=1,才能同时满足以上条件,即模型不能缩小,失去了模型实验价值。如果模型与原型采用不同的介质,有λυ=λν/λl=λl或3/2λν=λl(9.39)即实现流动相似有两个条件:一是模型流的流速为原型流流速的1/λl倍;二是必须按λν=3/2λl来选择运动粘性系数的比尺,但通常这一条件难以实现。从上述分析可见,一般情况下同时满足两个或两个以上作用力相似是难以实现的。实际中,往往仅考虑满足一个影响流动的主要作用力的相似,而忽略其他次要力的相似。(2)模型设计模型设计的步骤如下:①根据实验场地、经费、模型的制作条件和仪器、设备的量测条件确定几何比尺λl。一般情况下,按λl=10∽100选定。②根据几何比尺缩小原型的几何尺寸,得出模型的几何边界尺寸。③根据作用在原型流动上的主要作用力,选择模型律。如佛汝德模型,雷诺模型等。④按所选择的模型律推算各物理量的比尺,例如速度比尺、流量比尺、时间比尺等。从而由模型测得的各物理量的数据,推算原型液流中各物理量的相关数值。进而对设计方案进行优化,确定出安全可靠且经济的原型设计方案。3例9.4有一直径为15cm的输油管,管长5m,管中要通过的流量为0.18m/s,现用水来作2模型实验,当模型管径与原型一样,水温为10℃(原型用油的运动粘性系数νp=0.13cm/s),问水的模型流量应为多少才能达到相似?若测得5m长模型水管两端的压强水头差为3cm,试求在100m长的输油管两端的压强差应为多少(用油柱高表示)?解①因为圆管中流动主要受粘滞力作用,所以相似条件应满足雷诺准则,即(Re)p=(Re)m或υpdpυmdm=νpνm因为dp=dm,即λl=1,则上式可化简为:υpυm=νpνm又因Q=ωυ,而ωp=ωm,所以上式又可写成QpQm=νpνm22将已知油的νp=0.13cm/s,水的νm=0.0131cm/s代入上式,可得水的模型流量为νm0.01313Qm=Qp=0.18×=0.0181m/sνp0.13·236·
第9章量纲分析和相似原理②研究压强问题,须按欧拉准则,才能保证原型与模型压强相似,即(Eu)p=(Eu)mΔppΔpm或2=2ρpυpρmυmΔppΔpm或2=2γpυp/gpγmυm/gmΔpm因gp=gm,且已知模型测得压强水头差=3cm,则原型输油管两端的压强差(油柱)为γm2ΔppΔpmυp=·2γpγmυmπ2π2已知υp=Qp/dp=0.18/×0.15=10.19m/s44π2π2υm=Qm/dm=0.0181/×0.15=1.027m/s44所以lp=lm=5m长输油管的压差油柱为2Δpp10.19hp==0.03×2=2.95mγp1.027则在100m长的输油管两端的压强差为2.95×100=59m(油柱高)。53例9.5溢流坝的最大下泄流量为1000m/s,用缩小比尺λl=60的模型进行实验,试求模型中最大流量为多少?如在模型中测得坝上水头Hm为8cm,测得模型坝脚处收缩断面流速υm=1m/s,试求原型情况下相应的坝上水头和收缩断面流速各为多少?解为了使模型水流能与原型水流相似,首先必须做例9.5图到几何相似。由于溢流现象中起主要作用的是重力,其他作用力如粘滞阻力和表面张力等均可忽略,故要使模型与原型相似,必须满足佛汝德准则。根据佛汝德模型,流量比尺为5/25/2λQ=λl=60=27885则模型中流量3Qm=Qp/λQ=1000/27885=0.0359m/s=35.9l/s长度比尺为λl=60,则原型坝上水头为Hp=λlHm=60×8=480cm=4.8m流速比尺为λυ=λl=60=7.75则收缩断面处原型流速为υp=λυυm=7.75×1=7.75m/s·237·
水力学思考题9.1什么是物理量的量纲和单位?它们有何区别?9.2如何保证基本量的量纲是独立的?9.3简述瑞利法和π定理。9.4什么是几何相似,运动相似,动力相似?三者的关系如何?9.5试分别按雷诺准则和佛汝德准则导出下列各物理量的比尺(用长度比尺表示)。速度、加速度、流量、时间、力、压强、功、功率等。习题9.1用基本量纲L、M、T推导出力偶矩M,动能T,动量K及转动惯量J的量纲。9.2整理下列各组物理量为无量纲数:(1)τ,υ,ρ;(2)ν,l,υ;(3)F,l,υ,ρ;(4)σ,l,υ,ρ;(5)υ,g,l。9.3水泵单位时间抽送重度为γ的液体体积是Q,单位重量液体由水泵内获得的总能量为H(单位:米液柱高)。试用瑞利法证明水泵输出功率为P=kγQH。9.4实验观察与理论分析指出,水平等直径恒定有压管流的压强损失Δp与管长l,直径d,管壁粗糙度Δ,运动粘性系数ν,密度ρ,流速υ等因素有关。试用π定理求出计算压强损失的公式及沿程水头损失hf的公式。9.5如图所示的孔口出流,实验知道,孔口出流时,孔口断面流速υ与下列因素有关:孔口作用水头Η,孔口直径d,重力加速度g,液体密度ρ,动力粘性系数μ及表面张力系数σ。试用π定理推求孔口流量公式。9.6水流围绕一桥墩流动时,将产生绕流阻力F,该阻力与桥墩的宽度b(或桥墩直径d),水流速度υ,水的密度ρ,动力粘性系数μ及重力加速题9.5图度g有关。见图示。试用π定理推导绕流阻力表达式。9.7有一管径为200mm的输油管道,油的运动-52粘性系数ν=4.0×10m/s,管道内通过的流量是30.12m/s。若用直径为50mm的管道并以20℃的水题9.6图做模型实验,试求在流动相似时模型管内应通过的流量。若测得1m长模型输水管两端压强水头差为5mm,试求在100m长输油管两端压强差应为多少(用油柱高表示)?9.8有一处理废水的稳定池,池的宽度为25m,池长100m,池中水深2m,池中水温为20℃,水力停留时间15d(水力停留时间定义为池的容积与流量之比),成缓慢均匀流。设制作模型的长度比尺λl=20,在同种介质中实验,求模型尺寸及模型中的水力停留时间。(提示:按雷诺模型进行设计。)·238·
第9章量纲分析和相似原理9.9一桥墩长lp=24m,墩宽bp=4.3m,水深hp=8.2m,河中水流平均流速υp=2.3m/s,两桥墩间的距离Bp=90m,试取λl=50来制作模型,确定模型尺寸及其中的平均流速υm和流量Qm。9.10采用长度比尺λl=25的模型来研究闸下出流情况(如图所示),重力为流动的主要作用力。试求:①当原型闸门前水深Hp=14m时,模型中相应水深Hm=?②若模型实验测得闸下出口断面平均流速υm=3.1m/s,流量Qm=56l/s,由此推算出原型相应流速υp=?,流量Qp=?。③若模型中水流作用于闸门的力Fm=124N,问原型闸门所受的力Fp=?。9.11一溢流坝(参见例9.5附图)泄水流量为3150m/s,现按重力相似准则设计模型,如实验室供题9.10图3水量仅有0.08m/s,为这个模型选取几何比尺;原型坝高H0p=20m,坝顶水头Hp=4m,问模型最高为多少(H0m+Hm)?·239·
习题参考答案第1章92-821.1k=1.568×10N/m,β=6.38×10m/N1.2Δp=20000个大气压221.3v=0.01062m/s,μ=0.104N·s/m21.4τ=1150N/m21.5μ=0.054N·s/m31.6ρ=816kg/m31.7约16.30m1.83mm,1mm第2章2.1(1)p0=14.70kPa,(2)p0=11.03kPa2.23=14cm2.3pAabs=107.80kPa,pA=9.8kPa2.4(1)p1=19.7kPa,(2)pabs/γ=17m水柱,(3)pv=29.5kPa322.5γ1=6.86kN/m,pabs=106.33kN/m2.6H=1.14m22.7p0=-4.9kN/m,hv=0.5m水柱2.8h1=5m水柱,h2=0.38m水银柱22.9pv=19.6kN/m2.10hE=12.5m,hF=12.2m,hG=10.6m,hp=0.6m2.11(1)hA=4.6m水柱,(2)hA=0.81m水柱,(3)hA=0.3m水柱222.12pvm=27.07kN/m,p0=-38.05kN/m22.13p0=-27.9kN/m22.14(1)0.19m水柱,(2)0.784kN/m222.15(1)2.04kN/m,(2)2.94kN/m2.16F=27kN22.17p5=263.4kN/m2.1914.48m水柱,4.48m水柱,17.78m水柱,4.48m水柱·240·
习题参考答案32.20γ0=7.84kN/m22.21p0abs=6.664kN/m,z2=0.68m2.22H=0.4m2.23P=45.74kN,h0=2.03m2.24x=0.796m2.25F=32.67kN2.26M=933.4kN·m2.27PA=4.82kN,PB=7.05kN,PC=3.21kN2.28(1)T≥32.2kN,(2)T≥28kN3d3π2.29P=γ,yD=d12322.30Px=29.23kN,Pz=2.56kN2.31Pz=1.2kN方向向上2.32H=3.05r2.33P=2663.33kN,hD=6.91m2.34Pz=71.84kN2.35ΔH=2.52m,Px=19.4kN过球心向左2.36Px=82.05kN方向向右,Pz=79.92kN方向向上2.37Px=26.1kN,Pz=41.5kN2.38Pz=73.5N方向向上第3章3.5d0=25mm3.6υ2=3.18m/s3.7υ=7.5cm/s3.8Q2=19.4l/sγQ2=0.194kN/s33.9Q=1m/sυ1=2.5m/s3.10υ1=2cm/s3.11Q=102l/sυA=3.85m/s3.12p3abs=9mHO23.13υB=2.73m/sυA=φυB=2.68m/s3.14HA=7.816m>HB=5.051m水流由A向Bhw=2.765m3.15Q=51.18l/s3.16d1=9.8cm3.17Q=173l/s3.18H=5.186m3.19p1=-6731Pap2=16955Pa3.20Q=17.5l/spB=11.27Pa3.21p1=154.25kPad2=27mm·241·
水力学3.22(1)pB=29.4kPa(2)υ1=1.98m/s,υ2=3.96m/s33.23Q=0.55m/s3.24p2=10.06kPaQ=260l/s3.25d=75.5mmpB=-53.89kPa3.26hs=4.5-0.5-1=3m3.27H1=2.57m>H2=1.74m流向由A到Bhw=0.83m3.28h=7m3.29Q=67.3l/sp=79.2kPa3.30Q=27.1l/s3.31p2=44.1kPa3.32H=1.23mp23.33h=-=0.24mγ33.34Q=5.98m/s3.35F=100N(→)3.36R60°=252NR90°=504NR180°=1008N3.37θ=30°RX=456.5NRY3.38RX=3.815kNRY=3.415kNtanθ==0.895θ=41°81′RX3.39RX=0.462kN3.40R=1.70kN3.41RX=-1240.5kN(→)RY=-1750.9kN(↑)3.42R=153255.48N第4章4.1Qmax=0.471l/s4.2Re1/Re2=24.3Re=1929.6<2000层流4.4Re=7640>Rec=500紊流24.5λ=0.0332τ0=1.81N/m24.6(1)hf=15.02m(2)τr=0=0(3)τr=100=36.8N/m4.7d=0.15m4.8d=6.4mm4.9d=2.25mm4.10紊流光滑管-524.11ν=4.498×10m/s4.12紊流过渡区4.13(1)Qmax=12l/s(2)Qmin=174l/s4.14(1)Q=77.6l/s(2)hf=40.47mH2O·242·
习题参考答案4.15λ1=0.0279λ2=0.0339λ3=0.02914.16hf1=39.08mH2Ohf2=53.97mH2Ohf3=107.94mH2O4.17公式法hf=12.82mH2O查图法hf=10.1mH2O4.18(1)紊流过渡区hf=0.045m(2)层流hf=0.13m(3)紊流粗糙区hf=10.9m4.19λ=0.0352hf=8.19mH2O4.20Q=38.31l/s4.21t=10℃Δp=100.6kPaλ=0.0301hf=10.26mH2Ot=15℃Δp=104.45kPaλ=0.0315hf=10.6mH2O34.22(1)v=0.80m/s(2)Q=5.10m/s4.23莫迪图J=0.0070曼宁公式J=0.00894.24pmin=342kPa4.25H=43.9m4.26Q=25.4l/s34.27Q=0.159m/sh1=2.916mh2=7.084m4.28F=206N(→)4.29λ=0.04n=0.0114.30(1)Q=63.1l/s(2)Fx=0.783kN(→)Fy=0.130kN(↓)F合=0.794kN2(3)τ0=0.908N/m4.31dmax=0.868m0.54.32C=62.8m/sλ=0.0199v*=0.183m/shf=0.517mH2O4.33Q=45.8l/s第5章-23-23-235.1(1)Q1=3.73×10m/s;(2)Q2=2.16×10m/s;(3)Q3=2.366×10m/s5.2Q=16.1l/s5.3Q=193l/s35.4pv/γ=4.5m,Q=4.28m/s25.5p0=44.2kN/m5.6Q2=36.2l/s,H2=1.896m5.7ΔΗ=2.1cm5.8Q=3.11l/s5.9t=17.87s2D1.5+λl/d5.10t=2(H-H/2)d2g5.11t1=177.78s,t2=84.63s35.12Q=0.375m/s,pv=47.67kPa5.13d=1.0m·243·
水力学5.14Q=49.3l/s5.15Q=48.2l/s5.16Q=513l/s,hv=2.83m水柱5.17Np=35.84kW5.18H=21.8m,N=13.4kW5.19Hs=5.07m5.20d=0.5m,hs=4.28m,N=82.2kW5.21Q=109l/s5.22(1)Q=735l/s,(2)H=29.12m水柱5.23pA/γ=21.18m水柱5.24H=26.14m5.25H=10.42m5.261.2655.27hfAB=13.57m5.28Q1=57.6l/s,Q2=42.4l/s,hf=9.2m5.29H=19m5.30d=600mm5.31d1-2=150mm,d2-3=d3-5=100mm,d2-4=75mm5.32Ht=14m第6章6.1b=3.0m,i=0.00126.3h=3.92m,b=3.24m6.4i=0.0005436.5Q=33.66m/s,v=1.22m/s6.6h=1.49m6.7h0=0.6m,i=0.0020.56.8C=48.74m/s6.9h0=1.66m36.10Q=1.72m/s,v=0.512m/s6.11i=0.00033,v=1.63m/s6.12h=0.42m,b=4.15m6.13n=0.0326.14i=0.000256.15Δh=0.7m6.16i=0.00216.17h=2.47m,b=4.94m36.18ΔQ=6.01m/s6.19h=0.84m·244·
习题参考答案6.20hk=0.754m6.23hk=1.07m6.25hk=0.907m6.26ik=0.0056.27陡坡渠道6.28缓流6.29缓流、缓坡6.30急流6.31Fr<1,缓流6.33h″=1.59m,lj=8.9m,Δhw=1.12m6.34h″=1.584m6.37h上=1.72m第7章37.1Q=2.29m/s37.2Q=6.94m/s37.3Q=0.504m/s7.4b=17.2m,hmax=4.09m37.5Q=0.025m/s,H=0.263m37.6q=0.69m/(s·m)37.7Q=8.2m/s37.8Q=7.0m/s7.9b=9.9m第8章-48.1k=4.16×10cm/s-538.2Q=5.65×10m/s-58.3v=u=2×10cm/s228.4q=2.53米/日(m/d)-728.5q=3.75×10m/s-728.6q=5.71×10m/shc=2.60m8.7Q=67.2l/s-538.8Q=6.0×10m/shc=3.16m-438.9Q=4.32×10m/s8.10Q=13.4l/s8.11Q=6.4l/s8.12k=0.0032m/s-48.13k=8.82×10m/s,水位差=3m8.15R=184m·245·
水力学338.16Q=340m/d,Q=216m/d第9章2τυ·lFσυ9.2(1)2;(2);(3)22;(4)2;(5)ρυνρlυρlυgl2lΔν2Δνlυ9.4Δp=f(,,)·ρυ;hf=f(,)··ddυddυdd2g2πddμσ9.5Q=f(,,2)·2gH4HρHgHρgH22μgb9.6F=ρbυf(,2)ρbυυΔp9.7Qm=0.76m/s;()p=12.254mγ9.8Bm=1.25m;Lm=5m;Hm=0.1m;tm=54min39.9υm=0.325m/s;Qm=0.0914m/s39.10Hm=0.56;υp=15.5m/s;Qp=175m/s;Fp=1937.5kN9.11λ1=20.38;H0m+Hm=1.18m·246·
参考文献[1]清华大学水力学教研室编.水力学.上、下册.北京:人民教育出版社,1980[2]蒋觉先.水力学.北京:高等教育出版社,1993[3]成都科技大学水力学教研室编.水力学.上、下册.北京:人民教育出版社,1979[4]西南交通大学水力教研室编.水力学(第三版).北京:高等教育出版社,1983[5]刘鹤年编.水力学.北京:中国建筑工业出版社,1998[6]闻德荪,魏亚东,李兆年,王世和编.工程流体力学(水力学)(第一版).北京:高等教育出版社,1991[7]禹华谦主编.工程流体力学(水力学)(第一版).成都:西南交通大学出版社,1999[8]李士豪主编.流体力学(第一版).北京:高等教育出版社,1990[9]周谟仁主编.流体力学泵与风机(第二版).北京:中国建筑工业出版社,1985[10]黄儒钦主编.水力学教程(第二版).成都:西南交通大学出版社,1998[11]MelvynKay.PracticalHydraulics.E&FNSpon.animprintofRoutledge11NewFetterLane.LondonEC4P4EE,1998·247·'