水利工程论文-河工模型相似设计及特殊处理技术.doc

水利工程论文-河工模型相似设计及特殊处理技术摘要：本文从河工模型相似原理出发，由水流对床面产生的拖曳力和床面泥沙可动指标之比，导出适用于底沙的冲积类河工变态模型的变率e，是由模型砂密度来决定的关系式：=e2，进而确定了冲积河流的模型设计的变率e=L/h和=(sp-1)/(sm-1)之间的关系式来进行底沙模型设计。另外对一些特殊模型设计的处理技术作了较详细的介绍，对河工模型设计和试验研究者具有参考意义。关键词：河工模型相似律河床是河道水流与其底部泥沙相互作用的产物，在相互作用过程中逐渐形成适应水流的河床形态，水流的变化影响河床的冲刷和淤积变化，河床的冲淤变化也反过来影响水流变化，对于挟沙水流来说，水流和河床是一对矛盾的统一体，相互制约和作用，这种演变过程会使河床发生纵向和横向的变化。在平原冲积河流上，因水流的周期性变化使河床的这种变化过程也带有周期性的重复特点，其长时期观测是有规律可循的，自然现象是有一定相似准则可以模拟的，因次分析和相似原理是经常联合使用的，当研究的物理现象复杂，纯理论不足以解决这些复杂现象问题时，不得不借助于模型试验手段，模型试验就有相似问题，意味着模型与原型的无因次量必须相等。1理论关系式河床演变过程的基本相似，必须在模型上同时满足三个主要条件：a河床形态相似；b原型和模型的水流水力学相似；c原型和模型上的水流与河床泥沙相互作用过程相似。a要求是重要的，原、模型几何比尺是定值，b要求为原模型水力学动力相似，应遵循无因次数Fr相同，即Fr=v2/gh相等。v2p/ghp=v2m/ghm(1)式中下标p为原型值，m为模型值，这个条件就是重力相似条件，由此得v=h1/2(2)对于c要求是河床泥沙有两种相互保持不变的力，一种使河床泥沙颗粒脱离河底的力(上举力和悬浮力)，另一种使泥沙颗粒保持在河底的力(重力)，这两种力相互作用过程，表现出河床泥沙颗粒脱离河床而被水流挟带到新的位置，或保持不动。河床泥沙颗粒可动指标：以水流在床面上产生的拖曳力(剪应力)和泥沙起动指标来表征1。Kg=U2*/(s-)/gd=hJ/d(3)式中(s-)/=为泥沙的水下容重比，与科巴斯2的修正的颗粒弗劳德数=0/(s-)d是相一致的。因此有hpJp/pdp=hmJm/mdm(4)有d=hJ=2h/L(5)由U=C得v=c1/2h1/2J=1/2h(6)c=1/1/2J=(L/h)1/2由曼宁-斯崔克勒公式的经验公式1/n=26/d1/690，得c=-1/6d1/6h=(L/h)1/2(7)由此得d=4h/3L(8)将(8)代入(5)得=(L/h)2=e2(9)式中e为模型的变率。(9)式则是由模型砂的比重来决定的模型变率e，反过来也可以由模型变率e来确定模型砂材料的比重m。例：某河工模型，L=200,h=40,e=5,应用(9)式，=e2=25即(s-)p/(s-)m=25，(2.65-1)/(-1)=25,=1.066t/m3此模型砂以选用聚苯乙烯塑料砂=1.056t/m3为佳。2悬沙模型的相似律前面给出以底沙为主的模型设计简捷计算方法，但对于河工模型相似律中，由于挟沙水流的运动规律十分复杂，各家考虑的思路，及选用的起动流速及输沙率公式的不同，导致推导出的相似律也不尽相同，往往会存在诸多矛盾。但是，试验人员都得根据研究的课题抓住主要矛盾，保证主要方面相似，而对次要矛盾作些简化处理。对于悬移质泥沙运动相似，其相似律有两种思路，多数是从沉降相似来考虑，沉降公式按泥沙粒径d的Rew数来分个区域来表征：Rew0.2=(1/18v)(s-)/gd2为Stokes公式Rew103=1.720.2Rew103=(4/3)(1/Cd)(s-/)gd式中Cd为阻力系数，Cd=24/Red(1+3/16)Red由泥沙粒径d的不同，由不同的沉速公式导出的相似律也不尽相同。由于Cd系数的计算很繁杂，窦国仁3则将Cd简化为指数关系：Cd=/(Red)认为：Red0.25，(即d850,(d2mm时)，=0.43；=0；0.25Red850,(0.1mmd2mm)值在10之间变化，需要通过试算来定。由此可以导出：沉速比尺与粒径比尺关系：Re0.2=(s-/)2dRe103=1/2(s-/)1/2d0.2Re0.5mm=dvf=16.73(s-)/gd2/51/5h1/5(13)用(12)或(13)式来校验vf=v=1/2h,至于模型中其他相似律，不谈自明。3河工模型设计中特殊处理技术详细研究河工模型相似律之后，有许多特殊处理技术，可以达到事半功倍的效果。3.1采用原型砂粒径的模型相似律某些研究项目，由于原型砂很细，即便使用轻质砂，可以作到模型砂粒径变粗些，可能还是很困难，根据模型设计的相似原理，模型砂的粒径在一定程度上是可以变动的，也就是可以使模型砂粒径等于原型砂，即d=1。将其代入(8)式，d=4h/3L=1则有h=3/4L(14)v=1/2h=3/8L(15)变率e=L/h=1/4L(16)Q=Lhv=17/8L(17)t=5/8L(18)=/1/2L至于模型砂的材料，根据(9)和(16)式则有=e2=1/2L(19)这套处理方法往往会对这类模型设计带来许多方便简捷之处。3.2采用原型砂作为模型砂的相似律有时模型砂的选择还会有困难，采用原型砂在模型上作试验研究，即：d=1,=1(=1)由(5)式得(14)h=1/2L(14)v=1/2h=1/4L(19)E=L/h=1/4L(16)Q=Lhv=7/4L(20)t=Lv=3/4L(21)=(Lh)/p底沙输沙率采用以密实体积来计算的公式P=(1/200)d(v-vk)(v3/v3k)-1(d0.24mm适用)由上式得：p=dv=(L1/2L)dv=(3/2L)1/4L=5/4L(22)若采用窦国仁单宽输沙率公式pg=(K0/C0)(s)/(s-)(v-vk)(v/g)则有=/()4v/2C0因=1,=1,所以，=(h/L)v=1/2L；C0=(L/h)1/2=L(23)表明选用不同的输沙率公式，得到的冲淤时间比尺t1也不同，因而在应用中应注意选择适合当地输沙率公式来计算。以上均为以底沙为主的，采用原型砂粒径及原型砂的河工模型相似律。3.3“人工小河”自然河工模型生产中有时会遇到各种复杂课题，而且又受到环境及条件的限制，变态模型相似律表明可以使模型砂变粗些，这给研究者带来方便。当不允许作变态模型时，“人工小河”自然河工模型就可以达到上述的要求，其实质是略掉一些模型砂的相似律的条件，以模型量与原型量的比例来确定模型的其他比尺量，这方面笔者1987年曾承担漳卫新河一段险工段模型试验研究项目4，该险工段长度为4km，床沙很细d50=0.03mm，主要目的研究险工段丁坝冲刷问题，必须作成正态模型，按常规进行模型设计很困难，在L=h=80正态模型上，以水流相似原理施放一次洪水过后资料进行验证，采用“人工小河”自然河工模型原理，不管模型砂粒径相似，用浸泡饱和的木屑sm=1.16t/m3,d50m=0.50.6mm，来铺地形，施放该次洪水，试验过程中不断观测地形变化和冲刷深度，返复几次，达到与原型的地形和深度一致时，从中确定冲淤时间比尺t1和含沙量比尺s，以此决定了正式试验中的各类比尺，取得了满意的成果，该试验成果经第二年实际过水实测资料得以证实表明成果可靠，这种方法与黄委会李保如，张红武5自然河工模型的方法相类似。4系列模型延伸法沙玉清6早期发表了系列模型延伸理论，是以原型砂在不同几何比尺的模型上进行系列试验，然后将试验成果依比尺关系外延到原型(L=1)。所以系列模型至少要作两个以上的比尺模型，费时费力，这种方法比较适用于单体的局部冲刷和淤积的模型，其成果可信度是高的。在系列模型延伸法的基础上，为了克服要作几个不同比尺的系列模型的缺点，发展出系列模型砂延伸法。就是在同一个几何比尺模型上改变几次模型砂的粒径d，等于变更d，最后要求得以物理量和d点绘关系向外延伸到d=1时为最后结果，作者曾用此法作过几个求冲刷深度的模型试验，结果令人满意。这类方法适用于局部冲刷等模型试验研究项目，其优点是更换几次模型砂较作几个模型经济得多。近几年来又有学者发表了水沙模拟系列试验二次延伸法，其特点为：分两步处理的试验方法，一是按模型相似律设计计算出有关因素的模拟量，试验时不要求全部模拟量都满足，该模型为虚设的；二是以虚设模型作为模拟对象，以不满足相似要求的量值为系列进行系列试验，将试验结果先按指数关系延伸到虚设模型上，再按模型相似律延伸到原型，这种方法为二次延伸法。例如：某火电厂的专用供水水库（）人工围成长3km，宽1.5km，最大蓄水深7.43m，库容3000万m3，水库不间断向电厂输水1.5m3/s，每年两次向水库充水16.4天，入库Q=12m/s，Qs=0.5kg/m，而供水Qs=0.05kg/m，充水过程中泥沙在库内沉积以满足电厂供水含沙量要求，模型采用变态，以泥沙沉降相似为主，L=200,由于采用原型砂作为模型砂，=1=vhL,求得垂直比尺为34.2，这样绝大多数水、沙相似都能满足，唯模型流态达不到紊流。按重力相似，入库Q=12m/s的模型流量为0.31l/s；只有模型流态达到1.5l/s以上时才能达到紊流流态。试验中的固定进库含沙量为0.5kg/m3，分别以1.5、3.0、6.0l/s流量进行试验，分别量测出库含沙量，将三个不同的入库流量比尺与出库含沙量在对数纸上点绘关系，延伸得到虚拟模型(