泥沙运动力学
泥沙颗粒跃移运动机理

泥沙颗粒跃移运动机理白玉川;陈有华;韩其为【摘要】Saltation of solid grains is an important form of bed load transport. Force of movement particle was studied and sediment saltation movement model neglecting particle rotation and turbulence effects was established on the basis of the results. Using numerical calculation method to solve the motion equation, get the statistical properties of saltation parameters (height of saltation, length of saltation, velocity of saltation), which agreed well with the previous experimental data. Through regression analysis, the mathematical expressions of saltation parameters and sediment transport rate were obtained, which is very important for studying the mechanics of bed load transport and bed load transport rate.%水流中泥沙颗粒的跳跃运动是推移质运动的重要形式,根据泥沙颗粒运动过程的受力分析,建立了忽略颗粒旋转和紊动影响的泥沙颗粒跃移运动模型,并采用数值计算方法,求解了跃移运动方程,得到了跃移运动参数(跃移长度、跃移高度、跃移速度)的统计特性,所得结论与实验资料吻合较好;通过对数值计算结果回归分析,得到了跃移运动参数以及输沙率的数学表达式,对于了解推移质运动机理有重要的意义.【期刊名称】《天津大学学报》【年(卷),期】2012(045)003【总页数】6页(P196-201)【关键词】推移质;跃移;输沙率【作者】白玉川;陈有华;韩其为【作者单位】天津大学水利工程仿真与安全国家重点实验室,天津300072;天津大学河流海岸工程泥沙研究所,天津300072;天津大学河流海岸工程泥沙研究所,天津300072;天津大学河流海岸工程泥沙研究所,天津300072【正文语种】中文【中图分类】TV142泥沙颗粒在水流作用下发生3种运动模式:跃移、滚动(滑动)和悬移.床面泥沙以滚动、滑动还是跃移的方式运动,取决于水流强度和床面特性,当床面切应力超过临界起动应力时,床沙开始以滑动或滚动形式运动[1-4].如果床面切应力继续增加,颗粒开始从床面跃起进入跳跃运动状态,这些以滑动、滚动和跳跃形式运动的泥沙称为推移质.关于推移质的运动形态,文献[5-8]认为跳跃是推移质的主要形式,而文献[2,9-10]则认为泥沙起动的初始阶段以及当水流强度较小时,床面剪切力较小,泥沙则主要以滚动、滑动的形式运动.笔者已对以滚动、滑动形式运动的接触型推移质进行了研究,并取得了一定的成果[11-12].这里针对跃移型推移质进行研究,以达到全面了解推移质运动的目的.对于颗粒的跳跃运动,许多学者都做了大量的研究工作,Einstein[13]基于水槽实验得出推移质层厚度大约是2D,跃高大约是100D;Van Rijn[3]建立了颗粒单步跃移模型,利用上述研究者的实验数据检验此模型,证明其可适用于低雷诺数水流条件;惠遇甲等[14]利用高速摄像技术测量了床面附近颗粒的运动状态,并给出了滑动、滚动、跃移以及悬移各占的比例;Nino等[4,15]利用高速摄像技术测量了跃移轨迹并建立了二维连续跳跃模型;文献[16-19]对颗粒的跃移运动进行了一系列系统的研究,利用实时图像显示技术和高速摄像技术测量了跃移轨迹和速度,分别建立了单步跳跃模型、多颗粒连续跳跃模型以及三维连续跳跃模型.笔者分别建立泥沙颗粒跳跃运动起步阶段和飞行阶段的运动模型,并在此基础上分析研究跃移运动参数的统计特性以及推移质输沙规律.1.1 泥沙受力模型当水流流过松散河床时,单颗粒床面泥沙在水流的作用下,主要承受上举力LF、拖曳力DF、重力W.考虑单向流运动情况,设泥沙在水流作用下移动的纵向和垂向速度分别为u、v,水流作用在床面沙粒上的流速为bU,则上举力LF、拖曳力DF、重力W的表达式分别为式中:DC为泥沙颗粒拖曳力阻力系数;LC为泥沙颗粒上举力系数;D为泥沙粒径;ρ为水的密度;γ、sγ分别为水和泥沙的容重;1A、2A和3A均为泥沙颗粒面积修正系数.1.2 跃移型推移质运动的动力学方程1.2.1 跃移起步段(初始段)在力学分析的基础上,根据牛顿第二定律,建立泥沙跃移起步阶段运动方程为将拖曳力、上举力、重力的表达式代入,并引入希尔兹参数、颗粒雷诺数,同时令[CDA1+CLA2μ]=G′,式(4)、式(5)可以化简为式(6)、式(7)的各项系数具有明确的物理含义,为此引入如下物理概念[12-13]:颗粒惯性参数阻力参数颗粒运动的临界参数可以看出A、B、E、F均为颗粒雷诺数和y的函数,C、H均为颗粒雷诺数、希尔兹参数以及y的函数,故颗粒运动方程可以写为1.2.2 跃移飞行段跃移型推移质跃移飞行阶段,摩擦力等于零,运动方程可写为式中:相对流速ur=[(Ub-u)2+v2]0.5;(Ub-u)/ur和v/ur分别表示阻力与水流方向夹角的余弦和正弦.2.1 临界起动状态当泥沙颗粒处于临界起动状态时,此时可以近似忽略其垂向运动,只考虑纵向运动,纵向运动速度U=0,加速度.由运动方程式(8)可以得出C=0,由此可以得到很容易看出,式(12)其实就是代表着泥沙起动的希尔兹(Sheilds)曲线.2.2 临界扬动状态当泥沙处于由静止扬起临界状态时,此时可以忽略其纵向运动,只考虑垂向运动,垂向运动速度v=0,加速度=0,由运动方程式(9)得出H=0,进一步可以得到以沙粒雷诺数*Re为横坐标、希尔兹参数Θ为纵坐标绘出泥沙颗粒处于临界状态时的关系曲线,如图1所示.3.1 方程数值求解跃移颗粒飞行阶段运动方程是不能直接求得解析解的,必须借助于数值计算方法进行求解.采用四阶龙格-库塔法对方程组计算可得出轨迹、跃高和跃长等参数,在计算过程中做了如下的假定,同时对一些参数的确定方法也加以说明.1)假定计算中的颗粒是粒径为D的圆球,且为均匀密度.2)初始速度、的确定关于跃移初始速度,文献[20]认为u0/u*和v0/u*近似等于2,文献[21]认为跃移初始速度在u*∼2.5u*之间,文献[3]采用了2u*,文献[17]的分析中采用了2.5u*.基于上述前人的研究成果,本文计算中采用=2.5u*,=2.5u*.3)初始位置床沙颗粒在水流作用下开始运动,只有其跃起一定高度以后才认为是跳跃运动.胡春宏等[22]认为,当颗粒在运动中的最高点离开床面超过1,D高度后,即为跃移运动;Van Rijn[3]的研究中跃移初始位置采用0.6,D,初始位置为颗粒重心距床面的距离,因此与文献[14]的研究结果是一致的,故本文中跃移初始位置亦采用0.6,D.4)阻力系数DC和上举力系数LC阻力系数采用上举力系数LC受颗粒大小、形状、旋转角速度以及水流条件影响,一般可当作常量处理[23],文献[4,16-19]的分析中均采用文献[23]中建议的L0.2C=,本文分析中也采用上述结果.5)采用对数流速分布式中:12α≈;χ为与沙粒雷诺数有关的函数.3.2 量纲分析为了对跃移运动参数进行统计特性分析,同时便于与别人的实验资料进行比较,需要对颗粒跃移运动所涉及的水流条件和颗粒参数无量纲化.反映颗粒跃移运动的参数主要包括水深H、坡降J、水的密度ρ、泥沙粒径D、跃移泥沙密度sρ、运动黏性系数μ、颗粒表面粗糙度sk、床沙粒径bD以及重力加速度g.利用π定理,可得无量纲方程[17],即式中:F为无量纲跃移运动参数,包括跃高、跃长以及跃移速度;τ 为无量纲切应力,为摩阻流速,;*Re为颗粒雷诺数,.由于H≫D,Db≫ks,同时假定D=Db,则可以略去这3项,将τ、*Re合并,略去*U,可得参数.用(U*-U*c)代替U*.用参数代替τ,因此式(16)可以简化为3.3 数值计算结果及分析数值计算所得几组典型跃移轨迹见图2和图3.将数值计算所得跃移特征(包括跃高、跃长、跃移速度)与文献[4,16-17]的实验数据进行比较,跃高、跃长、跃移速度与*T的关系分布见图4~图6.可以看出,无量纲跃高、跃长均随*D和*T的增加而增大,但是当*T较大时,随*D的增幅程度变得不太明显.这说明,当水流强度较大时,无量纲跃高、跃长与颗粒形状的关系不显著.文献[5,24]指出,跃长仅与颗粒形状和大小有关,与水流强度无关,而本文中的分析结果则说明此结论只适用于强水流情况,在水流强度较低时,跃长是与水流强度有关的.由图7可以看出,数值计算得到跃高点长(跃移颗粒达到最高点时对应的跃长)约占整个跃长的35%,文献[17]的实验结果约为40%,文献[14]的研究结果为30%.根据数值计算结果,得到跃移运动参数的回归方程为式中HS、LS和VS分别为跃高、跃长和跃移速度.3.4 推移质输沙率推移质输沙率表示为式中bC为推移质浓度.根据文献[17],bC可表示为将式(18)、式(20)、式(22)代入式(21)中得推移质单宽输沙率为由于文献[25-26]中的公式均是基于泥沙跳跃机理建立的,因此将式(23)与文献[25-26]中的2个公式进行比较,如图8所示.从图8中可以看出Yalin公式[25]计算结果明显偏小,式(23)在低强度下计算结果较文献[25-26]中的公式计算结果偏大,而随着水流强度增大,计算结果与文献[25-26]中的公式计算结果吻合程度较好.(1) 针对跃移型推移质泥沙运动,从泥沙颗粒典型受力出发,分别建立了床面层内跃移质运动起步阶段和飞行阶段的理论模型.(2) 根据跃移颗粒起步阶段运动方程特性,从理论上推导得到了泥沙由静止起动、扬动的函数表达式,从理论上解释了希尔兹曲线的物理意义.(3) 利用数值计算方法求得的跃移轨迹与前人实验资料吻合较好,并且从中可以看出:颗粒跃移轨迹形状类似弹道,颗粒的跃移高度、长度以及跃移速度均随水流强度的增大而增大,在相同的水流强度下,粒径大的颗粒的跃高、跃长以及跃移速度要比粒径小的大;当水流强度较低时,无量纲跃高、跃长均随*D的增大而增大,跃高、跃长都与水流强度有关;当水流强度较大时,无量纲跃高、跃长随*D的变化不太显著;跃高点大约处于1/3的跃长处.(4) 根据数值计算结果,回归得到了跃移运动参数(跃高、跃长、跃移速度)的数学表达式,并进一步得到了相应的推移质输沙率公式,将所得输沙率公式与其他公式对比,说明其可以反映实际输沙情况.【相关文献】[1]钱宁,万兆惠. 泥沙运动力学[M]. 北京:科学出版社,1983.Qian Ning,Wan Zhaohui.Kinematical Mechanics of Sediment[M]. Beijing:Science Press,1983(in Chinese).[2] Peng Gao. Transition between two bed-load transport regimes:Saltation and sheet flow [J]. Journal of Hydraulic Engineering,ASCE,2008,134(3):340-348.[3] Van Rijn L C. Sediment transport(Part I):Bed load transport[J]. Journal of Hydraulic Engineering,ASCE,1984,110(10):1431-1456.[4] Nino Y,Garcia M. Experiments on saltation of sand in water[J]. Journal of Hydraulic Engineering,ASCE,1998,124(10):1014-1025.[5] Einstein H A. Formula for the transportation of bedload[J]. Transactions of the American Society of Civil Engineers,1942,107:561-597.[6] Sekine M,Kikkawa H. M echanics of saltating grainsⅡ[J]. Journal of Hydraulic Engineering,ASCE,1992,118(4):536-558.[7] Wiberg P L,Smith J D. Calculation of the critical shear stress for motion of uniform and heterogeneous sediments[J]. Water Resources Research,1987,23(8):1471-1480. [8] Wiberg P L,Smith J D. Model for calculation bed load transport of sediment[J]. Journal of Hydraulic Engineering,ASCE,1989,115(1):101-123.[9]韩其为,何明民. 泥沙起动规律及起动流速[M]. 北京:科学出版社,1999.Han Qiwei,He Mingmin. The Incipient Disciplinarian and Incipient Velocity ofSediment[M]. Beijing:SciencePress,1999(in Chinese).[10] Mazumder B S,Anindita B,Satya P O. Near-bed particle motion due to turbulent flow using image-processing technique[J]. Journal of Flow Visualization and Image Processing,ASCE,2008,15(1):1-15.[11] Bai Yuchuan,Xu Haijue,Xu Dong,et al. Nonlinear dynamical characteristics of bed load motion[J]. Science in China (Series E):Technological Sciences,2006,49(3):365-384.[12] Xu Haijue,Bai Yuchuan,Ng Chiu-On. Nonlinear mechanism of bed load transport[J]. Transactions of Tianjin University,2009,15(2):126-129.[13] Einstein H A. The Bed-Load Function for Sediment Transport in Open Channel Flows[M]. Washington,USA:Department of Agriculture,Soil Conservation Service,1950.[14]惠遇甲,胡春宏. 水流中颗粒跃移的运动学特征[J].水利学报,1991,4(12):59-64.Hui Yujia,Hu Chunhong. Saltation characteristics of particle motions in water[J]. Journal of Hydraulic Engineering,1991,4(12):59-64(in Chinese).[15] Nino Y,Garcia M. Gravel saltation (2):Modeling[J]. Water Resources Research,1994,30(6):1915-1924.[16] Lee Hongyuan,Hsu Insong. Investigation of saltating particle motions[J]. Journal of Hydraulic Engineering,ASCE,1994,120(7):831-845.[17] Lee Hongyuan,Chen Yenhsu,You Jiingyun,et al. Investigations of continuous bed load saltating process[J]. Journal of Hydraulic Engineering,ASCE, 2000,126(9):691-700.[18] Lee Hongyuan,You Jiingyun,Lin Yingtien. Continuous saltating process of multiple sediment particles[J]. Journal of Hydraulic Engineering,ASCE,2002,128(4):443-450.[19] Lee Hong-Yuan,Lin Ying-Tien,You Jiingyun,et al. On three-dimensional continuous saltating process of sediment particles near the channel bed[J] Journal of Hydraulic Research,ASCE,2006,44(3):374-389.[20] Abbott J E,Francis J R C. Saltation and suspension trajectories of solid grains in a water stream[J]. Philosophical Transactions of the Royal Society,1977,284(1321):225-254.[21] White B R,Shulz J C. Magnus effect in saltation[J]. Journal of Fluid Mechanics,1977,81(3):497-512.[22] Hu Chunhong,Hui Yujia. Bed-load transport(Part I):Mechanical characteristics[J]. Journal of Hydraulic Engineering,ASCE,1996,122(5):245-261.[23] Wiberg P L,Smith J D. A theoretical model for saltating grains in water[J]. Journal of Geophysical Research,1985,90(4):7341-7354.[24] Einstein H A,E1-Sammi E A. Hydrodynamic forces on a rough wall[J]. Reviews of Modern Physics,1949,21(3):520-524.[25] Yalin M S. Mechanics of Sediment Transport [M]. New York:Pergamon Press,1972.[26] Engelund F,Fredsoe J. A sediment transport model for straight alluvial channels[J]. Nordic Hydrology,1976,7(5):293-306.。
第三节河流泥沙运动力学的理论与实践

第三节河流泥沙运动力学的理论与实践中国河流,特别是北方河流自古多沙。
多沙河流的水流运动有其特殊性。
这种区别于清水河流的特殊性使多沙河流的治理变得更加复杂。
古代对河流泥沙运动的理论认识起源于春秋战国时期,在两汉取得突出的进步,经过北宋的发展,至明代后期达到高峰。
这些大多处于定性阶段的理论认识在古代世界是领先的,其中一些还在治河实践中得到应用。
一泥沙运动力学的起源与张戎的贡献(一)对泥沙运动的观察水流有清浊之分,古人早有观察和记录。
《诗•小雅•谷风之什》曰:“相彼泉源,载清载浊。
”战国时人解释河水变浊的原因是:“夫水之性清,土者汨之,故不得清。
”②《尔雅•释水》并且具体解释黄河之所以含沙量高的原因是“河出昆仑墟,色白,所渠并千七百一川,色黄。
”晋代学者郭璞(273~324)对黄河之所以黄注解说:“潜流地中,汨漱沙壤,所受渠多,众水溷淆,宜其浊黄。
”即黄河之浊是众支流挟沙汇入所致。
水流含沙是水流的搬运作用。
北宋科学家沈括(1031~1095)曾指出,以黄河为首的华北诸河多浊流,河水挟带的泥沙都是上中游被冲蚀的土壤。
所以他的结论是:“所谓大陆者,皆浊泥所堙耳。
”③并且举例说,舜杀治水失败的鲧于羽山,这个羽山当时是在东海里,现在则在大陆,以为证明,正确地解释了华北平原的成因。
他还认为,浙江温州的雁荡山诸峰挺立,“原其理,当是为谷中大水冲激,沙土尽去,惟巨石岿然挺立耳”。
④那么,泥沙是否能被流水携带,还要看水流的速度与泥沙的粒径和比重。
东汉初年王充(27~97)指出:“湍濑之流,沙石转而大石不移,何者?大石重而沙石轻也。
”⑤而相反的情况则是大哲学家老子所概括的:“浑兮其若浊,孰能浊以止?静之徐清。
”⑥当挟沙水流静止下来的时候,便失去拖曳力和上举力,原来被流水挟带的泥沙随之沉降下来。
在春秋战国时期的水利建设中,已实际应用了水流冲淤的概念。
《考工记.匠人》中说:“凡沟,必因水势;防,必因地势。
善沟者,水漱之;善防者,水淫之。
泥沙运动力学

泥沙运动力学水力学基础持续介质:流体是由持续散布的流体质点组成的介质。
粘性力:由于存在内摩擦,一层流体对相对运动的另一层流体产生阻力。
牛顿内摩擦定律:牛顿流体粘性切应力与流体切边率成正比关系。
拉格朗日法:着眼于流体质点,跟从流体质点一路运动,记录流体质点在运动进程中各类物理量随所到位置和时刻的转变和规律。
欧拉法:着眼于空间点,把流体物理量表示为时空位置和时刻的函数。
流体无力量在不同时刻的时空散布。
迹线:流体质点在空间运动时所刻画出来的曲线叫做轨迹。
流线:流线是如此的曲线, 于某一固按时刻,该曲线上各点的速度方向与该点的切线方向一致。
理想流体:不可紧缩的、没有粘滞性称为理想流体。
持续方程: 不可紧缩流体:密度为常数,那么和时刻无关。
那么可紧缩定常流:变量与时刻无关。
那么雷诺数:Re=惯性力/粘性力弗洛伊德数:Fr=惯性力/重力伯努利方程:由于有粘滞力因此总能量必然是慢慢减少的。
泥沙特性1.泥沙来源泥沙:在流体中运动或受水流,风力,波浪,冰川和重力作用移动后沉积下来的固体颗粒碎屑。
泥沙的来源:岩石的风化是泥沙的重要的来源。
它包括机械的分离和化学的分解两个方面。
泥沙矿质的组成:长石,石英,云母石,高岭土,氧化铁泥沙特性 有泥沙颗粒的特性和泥沙群体的特性两种。
泥沙颗粒的特性要紧有:①重度,单位体积泥沙颗粒的重量,以千克/米3表示,其数值随泥沙的岩性不同而异,矿物成份主若是石英和长石,泥沙的重度一样约2650千克/米3。
②粒径,泥沙颗粒大小的一种量度,有不同方式表示。
经常使用的有等容粒径即体积与泥沙颗粒相等的球体的直径;筛径,即用具有不同孔径的标准筛,对泥沙进行分筛求出的粒径;沉降粒径,即依照粒径与沉降速度的关系算出的粒径等。
③沉速,指泥沙颗粒在无边界静水内的沉降速度,以米/秒或毫米/秒表示。
它也可作为泥沙颗粒大小的一种量度,故又称泥沙的水力粗度。
沉速综合反映颗粒和水的特性,因此是泥沙运动的一个重要参数。
七泥沙的推移质运动

1 概述
按照运动形态的不同,水流中运动的泥沙可分为悬移质和推移质两个部分。推移质通常沿河 床滚动、滑动、跳跃前进,动静间歇,与床面静止的泥沙时常交换,前进速度远较水流速度小。 而悬移质泥沙是在水流中浮游前进,其顺水运动速度与水流运动速度基本上相同。由于推移质沿 着床面运动,运动机理复杂,参数难以精确测量,不确定因素较多。因此,与悬移质相比,推移 质是泥沙研究中相对薄弱的课题。本调研将从泥沙起动、沙波运动、动床阻力、推移质输沙率和 推移质实(试)验等几个方面,对推移质的研究发展动态进行跟踪阐述。
究院,2010)。 对于床沙粒径比较均匀,用中值粒径或平均粒径来计算起动条件,能与实际基本相符。但对
于床面泥沙级配很不均匀的情况,不可能用一个代表粒径来计算起动条件,为此引进临界起动概
率的概念,与之相对应的泥沙颗粒粒径被称为临界起动粒径。根据 Meyer-Peter 和 Einstein 的推 移质输沙率公式,把输沙率很小(接近于零)作为泥沙起动的水力条件,由此得到的临界水力强
=
C
这里,Wc 是泥沙起动的水流功率,C 是常数。
对于某一河床比降和天然沙,式(3)可写为:
qc
=
K
d3/2 S
这里,qc 是泥沙起动的单宽流量,S 是河床比降。
(3) (4)
根据实验资料,Schoklitsch(1934)得到了与式(4)类似的泥沙起动公式,即 qc = 0.0194d / S 4 / 3 ,
动流速 Uc(垂线平均流速)可表示为:
Uc = η
ρs
− ρ
ρ
gd
⎜⎛ ⎝
h d
⎟⎞m ⎠
(1)
式中,ρs, ρ分别是泥沙和水的密度,g 是重力加速度,d 为泥沙粒径,h 是水深,m 是指数,η是 综合系数,可由起动流速实测资料反算。Shamov(1952)根据实验资料,求得η和 m 分别为 1.14 和 1/6。在这类公式中,李保如(1959)公式也有一定的代表性。若垂线流速分布采用对数形式, 则代表性的公式有 Goncharov(1962)和 Levy(1956)公式。
第三章 泥沙特性

泥沙颗粒)的双电层。
+
+ +
+
+吸
++ + +
附
层
扩散层
内 泥沙颗粒 层 外层 中性水
双电层
32
1、泥沙颗粒周围的双电层
2、双电层的外层 ①、吸附层 ②、扩散层
①、吸附层(固定 层,不活动层)
紧靠内层的反离子, 由于受静电引力大, 便与颗粒表面牢固 地结合在一起,称 吸附层。
++ + +
+
+
+ ++
c很缓,粒径变 化范围大,各组 粒径含量接近, 组成不均匀,级 配良好
请想想:砼搅拌中要求组成沙级配良好,对应均匀沙或非均匀沙?
第三章 泥沙特性
1.1 泥沙的几何特性
(二)沙样组成与粒配曲线
第三章 泥沙特性
1.1 泥沙的几何特性
(二)沙样组成与粒配曲线
1.1.2
第三章 泥沙特性
1.1 泥沙的几何特性
+ +
+ +
+ +
+
+ +
+
+吸
++ + +
附
层
内 泥沙颗粒 层 外层
双电层
扩散层 中性水
33
1、泥沙颗粒周围的双电层 2、双电层的外层 ①、吸附层 ②、扩散层
②、扩散层(活动层)
扩散层:距内层较远的 反离子与颗粒表面结合 的就不牢固,具有一定 的活动性,这一层叫做 扩散层。
++ + +
河流泥沙动力学习题及答案

河流泥沙动力学习题1.某河道悬移质沙样如下表所列。
要求:(1)用半对数坐标纸绘出粒径组的沙重百分数P 的分布图,绘出粒径的累积分布曲线,求出d 50、d pj 、ϕ(2575d d =)的数值。
(2)用对数概率坐标纸绘出粒径组的沙重百分数P 的分布图,绘出粒径的累积分布曲线,求出d 50、ϕ的数值。
(3)用方格纸绘出粒径组的沙重百分数P 的分布图,绘出粒径的累积分布曲线。
解:根据题意计算出小于某粒径之沙重百分数,列表如上。
(1)、半对数坐标纸上粒径组的沙重百分数P 的分布图及粒径的累积分布曲线,从下述半对数坐标纸上的粒配累计曲线上可查得中值粒径m m 054.050=d ,m m 075.075=d ,m m 041.025=d 。
平均粒径:069.01008675.6141141==∆∆=∑∑==i ii iipj pdp d , 非均匀系数:353.1041.0075.02575===d d ϕ。
半对数坐标纸上的沙重百分数p的分布图2468101214161820220.010.11粒径(mm)沙重百分数(%)半对数坐标纸上的粒配累积分布曲线1020304050607080901000.010.11粒径(mm)小于某粒径之沙重百分数(%)(2)、对数概率坐标纸上粒径组的沙重百分数P 的分布图及粒径的累积分布曲线 (3)、方格纸上粒径组的沙重百分数P 的分布图及粒径的累积分布曲线方格纸上的沙重百分数p的分布图2468101214161820220.050.10.150.20.250.3粒径(mm)沙重百分数(%)方格纸上的粒配累积分布曲线10203040506070809010000.050.10.150.20.250.3粒径(mm)小于某粒径之沙重百分数(%)2.已知泥沙沉降处于过渡区的动力平衡方程式为(ω可查表):223231)(ωρωρυγγd K d K d K s +=-令上式为 A=B+C要求计算并绘制d ~C B C +及d Re ~C B C +的关系曲线。
水流挟沙力计算公式研究综述

[5] [4 ]
U ω
0.23
1
Δ h h
0 .04
d2 5 d7 5
0.1 06
(9 )
赵龙保[1 1]根据椒江口的实测资料,引入前期含沙量的概 念,获得了较为理想的挟沙力经验公式: S f (e ) = a ( U f (e ) ) (Se ( f ) )
m n
(10 )
(2 )
ρ Uh m η
U
2
U
0 . 33
B h
0 .3 2 . 85 μ r
gh
ω
舒安平[35 ]从挟沙水 流紊动能量平衡时均方程理论出发, 建立了高含沙水流挟沙能力公式:
Sv = P f (μ ) r k2 fm 8
3/ 2 N
γ U3 m γ γ gR ω s m
曹如轩[ 3 6]利用南科所以及黄河等室内外实测资料,分粒 径组考虑,分别讨论了高含沙水流,建立了以下公式:
/( gR ω )
3
]
1.5 1.5
]
(2 1 )
S*i = k μ r
α
ω i
β
(Pi
+ P KDPbi )
U3
2.波、流共同作用下的挟沙力公式
刘家驹
[21 ]
γ γ s gRω i γ
(32 )
认为淤泥 质海岸的主要动力因素是以特征流速 吴伟明[3 1]运用床面泥沙交换机理对平衡输沙水流挟沙力 进行了探讨, 建立了平衡输沙时的非均匀沙分组挟沙力公式: (2 2 )
Z
第 11 卷
泉与窦国仁的研究,得出的长江口半潮平均挟沙力公式为:
S * = K1 U U + K2 gh ghω
河流动力学第四章 推移质运动

沙莫夫公式
=1.144
m=1/6 适用范围:
D>0.2mm
岗恰洛夫公式
对数流速分布 适用范围:
0.08-1.50mm
1
Uc 1.144
s
gD
(
h D
)
6
Uc
1.07 lg
8.8h D95
s
gD
§4.2. 泥沙的起动
三、无粘性均匀沙的起动拖曳力
起动拖曳力
0
hJ
U
2 *
Krammer方法:定性标准
最常用的方法 具体内容
♥ 无泥沙运动:静止 ♥ 轻微的泥沙运动:个别动,可数 ♥ 中等强度泥沙运动:少量动,不可计数 ♥ 普遍的泥沙运动:普遍动,床面变形
§4.2. 泥沙的起动
五、与泥沙起动有关的几个问题
泥沙起动具有随机性 泥沙条件 ♥ 大小、形状 ♥ 级配、密度:均匀沙,非均匀沙 ♥ 床面平整、颗粒排列 水流条件 ♥ 水流的紊动 ♥ 流速的大小
♥ 推移质运动达到一定规模,床面起伏 ♥ 泥沙颗粒在床面的集体运动
用途
♥ 推移质运动的一种主要形式 ♥ 构成河床地形的基本元素 ♥ 影响:水流结构,河道阻力,泥沙运动,河床演变
主要内容
沙波形态和运动状态 沙波的产生和消亡
§4.3.1. 沙波形态和运动状态
沙波介绍
名词:波峰、波谷、波长、波高 特点:迎水面:较为平坦、背水面:相对较陡
♥ 悬移质中的较粗部分 ♥ 推移质中的较细部分
同一泥沙组成:表现不同
♥ 水流较强时:悬移质 ♥ 水流较弱时:推移质
§4.1. 泥沙运动的形式
推移质与悬移质间的转换过程
悬移区 床面层 层移区 河床
(悬移质)
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
2018年博士生入学考试科目《泥沙运动力学》
考试大纲
一、要求了解和掌握的内容
C1、了解泥沙运动力学的研究内容,掌握泥沙问题在生产实践中的重要性。
C2.1、掌握泥沙的来源以及机械分离、化学分解的类型和原因,了解岩石风化的速度及产物。
C2.2、掌握泥沙群体特性、粒配曲线体现的特性以及平均粒径、中值粒径、等容粒径等计算确定方法,了解细颗粒泥沙的物理化学特性、泥沙的分类。
C2.3、了解浑水的一般性质,掌握浑水容重、含沙量的表示方法和计算公式。
C3、掌握泥沙的沉降过程和沉降状态,掌握球体颗粒泥沙在层流、紊流状态下的重力、阻力的计算方法和公式以及沉速的计算公式,了解天然泥沙颗粒沉速的影响因素和计算方法。
C4、掌握紊流的基本性质、紊动切应力的概念、水流切应力的垂线分布,了解紊流经典理论。
C5、掌握无粘性泥沙颗粒主要受力分析及一般表达式、泥沙运动的主要形式、划分推移质和悬移质的意义及其物理本质的区别、床沙质和冲泻质的概念及其划分的意义。
了解粘性泥沙的受力特点。
C6、掌握沙波的概念及其发展消长的阶段和过程、研究沙波运动意义,了解沙波形成的机理和原因、床面形态的判别方法。
C7、掌握近壁流区的重要性、河道阻力的组成单元、河床河岸及综合阻力的划分及计算方法、垂线流速对数分布公式(7.38)、垂线流速指数分布公式(7.38),了解水流能量及损失特性、沙粒及沙波阻力的确定方法。
C8.1、掌握泥沙起动的基本概念、起动的表达方式,了解泥沙起动的
判别标准。
C8.2、掌握希尔兹(Shields)泥沙起动曲线的特点、起动公式及其推导过程,了解斜坡上的泥沙起动
C8.3、掌握按滚动和滑动模式的起动流速受力特性和推导过程,掌握沙莫夫起动流速公式,了解其它流速起动公式以及粘性泥沙的起动。
C9.1、掌握推移质输沙规律的流派及其代表、拜格诺提出的推移质运动遵循的规律、爱因斯坦推移质输沙理论的主要认识。
C9.2、掌握拜格诺、爱因斯坦推移质输沙率公式的推导思路和过程,了解以流速为主要参数的推移质输沙率公式推导过程。
C9.3、掌握拜格诺、沙莫夫推移质输沙率公式以及无因次水流参数、无因次水流强度参数、无因次单宽输沙率的表达式、非均匀推移质输沙率的确定方法。
C10.1、掌握悬浮泥沙运动机理、泥沙扩散方程,掌握含沙量沿水深分布的扩散理论(10.16)的推导过程、各变量含义及其不足,了解重力理论的基本思路及其推导过程。
C10.2、掌握含沙量与悬移质输沙率之间的关系以及爱因斯坦、维利卡诺夫、拜格诺悬移质输沙率公式及其推导思路和原理,了解其推导过程。
C10.3、掌握各向泥沙扩散系数的确定方法。
C11、掌握水流挟沙力公式的主要类型、代表人物以及影响水流挟沙力的主要因素,了解水流挟沙力的确定方法。
C12、掌握泥沙的存在使水流紊动强度增强或减弱的原因、了解泥沙的存在对流速分布的影响。
C13、掌握高含沙水流的流动特性和泥沙运动特性,泥石流的类型、形成条件和运动特征。
C14、掌握异重流的基本概念、类型、研究意义、物理本质、形成条件、运动特性、泥沙淤积,掌握水电站取水口及船闸引航道泥沙淤积机理和防治技术。
二、考试形式与试卷结构
(一)考试形式
考试形式为笔试,考试时间为3小时,满分为100分。
(二)试卷结构
主要含论述和综合计算题等大题,总题量8~12题。
三、主要参考书目
钱宁、万兆惠,《泥沙运动力学》,科学出版社,2003年。