泥沙沉速公式

泥沙沉速公式研究泥沙沉速时很多方法都是从球体出发，但天然泥沙并非球体，它在下沉时受到的阻力比球体大，其阻力系数通常根据实验确定，关于泥沙的沉速，中外学者提出不少计算公式。岗恰洛夫公式⑴层流区(D<0.15皿)：3=1丫sTgD224F—v(2)紊流区(D(2)紊流区(D>1.5rn):gD⑶过渡区(0.15<D<1.5m)：考虑层流区和紊流区沉速公式的结构形式，认为在过渡区来说，几个主要变量的次方，应该介于层流区与紊流区之间。考虑量纲法则得到过渡区沉速公式pEV1/3B反映粒径和温度变化改变粘滞性影响的一个附加因素。D0=0.15cm,计算时D应与D0的单位一致。1/12

沙玉清公式(1)层流区(D<0.1的)：①=1丫s7gD224f-v(2)紊流区(D>2胸)：3=1.143=1.14gD(3)过渡区(0.1<D<2胸)：引进两个无因次判数，一个名为沉速判数Sa,是沙粒雷诺数sD/v及阻力系数CD的函数。A1/3g1/3v1/3g1/3)另一个名为粒径判数①，是沙粒雷诺数与沉速判数的函数。3DVSTag1/3。丫3DVSTaV2/3适用条件：当粒径为0.062到2.0mm时张瑞瑾公式泥沙下沉时的有效重力w=K(Y -Y)D31SK1为泥沙体积系数泥沙下沉时颗粒所受阻力F=KpuD3+KpD232K2和K3都是无因次系数由W=F得到Ku3=— 2Ku3=— 2 +KD3I2KTD' 3 'KY-Y+ 1SKY3gD并通过大量的泥沙实测资料得出:gD-C1=13.95,C2=1.09则TOC\o"1-5"\h\z\o"CurrentDocument"7 _u\2V-y -u①二13.95 +1.09's'gD-13.95\o"CurrentDocument"* D） ~1~~ D此公式可满足层流区、紊流区和过渡区的要求。（1）层流区：1y-yd23= sg25~6y-o-（2）紊流区：3=1.044，yS―ygD窦国仁公式①既然在过渡区，绕流阻力有表面阻力和形状阻力两部分组成，所以过渡区总阻力F总二F表+F形。②绕流属层流状态时，因为泥沙颗粒背后不发生分离现象，在泥沙的表面全部处于滞性状态泥沙所受的阻力是F表。当进入过渡状态时，泥沙表面一部分面积仍然受到表面阻力F表而剩下的那一部分的面积由于水流的分离现象，产生了形状阻力。根据实验观察，这个分离区域的大小，可以用分离角。来表示。所以F表、F形〜。有关。

③雷诺数ReD表示惯性力与粘滞阻力的对比关系，而过渡区绕流阻力的变化过程正是粘滞阻力和惯性阻力的互相消长过程，因此这个分离角。应该与雷诺数有关。=f(ReD)。d9_A

dRR9随着ReD的增大而增大，但增长率随ReD的增大而减小。用微分方程表示即 eDReD9=Aln-eD(A比例系数一一常数)积分后 。(C积分常数)9=lg4ReD当ReD三0.25时，泥沙周围没有分离现象，这时9=0,绕流呈层流状态。当ReD叁350时，泥沙周围布满着涡体，此时9二n,绕流呈紊流状态。把这两个边界条件代入得到：C=0.25；A=0.434.•.9=2.3x0.4341nRD=lg(4R)0.25eD……(9以弧度计)但规定ReD<0.25,取9=0；ReD>350,取9二n。④建立F表、F形〜9之间的关系3F=K3兀目加(1+ R)•①表面阻力K一颗粒形状的修正系数。球体K=1；泥沙K=4/3。少一关于泥沙表面受粘滞阻力作用的修正系数，实际上是分离区外的球冠面积与整个球面积的比值。球冠面积=2球冠面积=2兀rh=兀DfD+Dcos9〕122J咚G+cos9)过渡区阻力系数：32ReD(1+31一R)_(1+cos9)+1.2sm2916eD，2兀D3/ 、W二JL)令W=F总，解出:适用条件：ReD=0.25〜350;9=0〜n。当丁二15℃,相应泥沙直径口=0.2~2.0咂。当ReD<0.25,取9=0；当ReD>350,取9二n。泥沙的起动流速公式泥沙颗粒由静止状态转为运动状态时的临界水流条件,称起动的条件。起动流速：当用水流垂线平均流速来表示起动条件时,该流速称为临界起动流速。(1)无粘性均匀沙的起动流速公式考虑的条件：泥沙粒径是均匀的；颗粒之间没有粘结力促使颗粒运动的力是水流推移力FD和上举力FL,抗拒其运动的力是泥沙的有效重力W。u2F=CaD2yDD1 2gu2F=CaD2y_o_LL2 2g量纲为kgm/t2,a1=n/4；a2=n/4。CD,CL推移力及上举力系数。al,a2垂直于水流方向及铅直方向的沙粒面积系数。

w=a3(ys-y)D3a3沙粒的体积系数，a3=n/6，yS：沙粒的容重。如沙粒采取滚动形式起动，则起动临界条件的动力平衡方程为K1DFD+K2DFL=K3DW式中K1D、K2D、K3D为FD、FL、W的相应力臂，得到：K1DCa1D2yK1DCa1D2yU2 0 F2gDC2La2D2yU2 o—2gDa33(yK3a3(ys—K3a3(ys—y)DKCay—0—+KCay——1D1 2g 2L2 2g2K3a3(ys-y)gD(K1C：a1+12CLa2)y2Ka。 /y「y0CK1CDa1+3K32CLa2' $Y0CU0C表示起底流速（或临界底流速）用垂线平均流速来表达代替底流速，采用指数流速分布公式有U一垂线上某点（y点）的速度（时均值）y—该点离开床面的距离U—y=h处即水流表面流速h—水深 m—指数

则垂线平均流速则Uh=(1+m)U，代入前式有U=(1+m)U如取y=aD处的流速作为作用在沙粒上流速的平均值，有U=(1+m)amU0=nn为综合系数综合系数n只有通过起动流速的试验资料反求。沙莫夫公式根据试验求得n=1.144,m=1/6。=1.144gD对于天然泥沙平均重率为2.65t/m3,上式可以简化为U=4.6D1/3h1/6当床面大量起动时

U=6D1/3hi/6C适用范围D>0.2的岗恰洛夫公式取对数形式流速分布得到岗恰洛夫公式取对数形式流速分布得到8.8h=1.07lg8.8h=1.07lg-D~~95适用D=0.08〜1.5的(2)粘性泥沙的起动流速公式仁公式U=仁公式U=0441g11于加0.19个广公式单位为g、cm、s制。8=0.213x10-4cm3为薄膜水参数，可取

及为粘结力参数'k及为粘结力参数'ks为河床粗糙度，=2.56cm3/s2d<0.5mm, =0.5mm=0.05cm/>0.5mm,Ks=d张瑞瑾公式\l/2丫一丫 710+*\l/2— 17.6— d+6.05x10-7 -0.72JCl=l.34,C2=0.00000496,s=0.72,m=0.14,参考粒径d=lmni,将以上参数代入式（3-u）,并将公式中各变量均采用国际单位（kg、m、s制），得U=1Z6‘S7+6.05x10-7c⑺I~1~ 印")唐存本公式数m对于天然河道为1/6,4j7[hJ对于平整床面（如实验室水槽或d<0.01nmi的天然河道）为24j7[hJps：考虑空隙在内的床面泥沙的实际干密度;p；。：床面泥沙达到密实后的稳定干容重，其值为1.6t/m3oC：粘结力系数,C=2.842x1042V/mC：粘结力系数,C=2.842x1042V/m沙玉清公式0.267y/4+