钢筋石笼起动流速试验与流场结构数模分析.docx

1、第47卷第2期2014年2月Vol.47No.2Feb.2014天津大学学报(自然科学与工程技术版)JournalofTianjinUniversity(ScienceandTechnology)DOklO.ll784/tdxbz201206059钢筋石笼起动流速试验与流场结构数模分析叶恩立I,周宜红心，肖焕雄I,任磊3(1.武汉大学水资源与水电工程科学国家取点实验室.武汉430072;2. 三峡大学水利与环境学院.宜昌443002;3.爱尔兰国立大学工程与信息学院，高威)摘要：介绍了钢筋石笼起动流速的研究现状，针对其不足，用石茏实际边长代誉化引直径，经过理论推导.提出了以拖拽力系数为主要待定

2、参数的钢筋石笼起动流速简化公式，并在此基独上进行了axaxc形石茏的起动流速水槽试验.通过试验，分析了扃度与稳定性的关系，得到了不同扁度所时应的拖拽力系教值并讨论了底部犀擦力系数，给钢筋石茏建定性所带来的重要影响.此外，为更好地了解钢筋石笼周图流场的特桂，构定了一个透水性绕流数值模型，并以大涡模拟为计算手段.对透水性绕流场进行了三维数值模拟.通过将模拟站果与非透水性绕流场从流场结构、涡体脱落、静压强分布以及拖拽力系数等方面进行定性、定量的对比分析，得到了透水性对流场的具体影响结果.关键词：钢筋石笼；起动流速；水槽试验；方钝体透水性绕流；大涡模拟中图分类号：TV13I.61文献标志码：A文章编号

3、：0493-2137(2014)02-0108-08CriticalVelocityTestandFlow-Field-StructureNumericalAnalysisofRock-FilledSteelCageUnderWaterYeEnli',ZhouYihong1,2,XiaoHuanxiong1,RenLei3(I.StateKeyLaboratoryofWaterResourcesandHydropowerEngineeringScience,WuhanUniversity,Wuhan430072.China；2.CollegeofHydraulicandEnvironm

4、entalEngineering,ThreeGorgesUniversity,Yichang443002,China；CollegeofEngineeringandInformatics.NationalUniversityofIreland,Galway,Ireland)Abstract:Anintroductiontothestatusofresearchinthefieldofthecriticalvelocityofrock-filledsteelcageisdescribed.Consideringthedeficienciesofstudiesinthepast,aconcisef

5、ormulainwhichthedragcoeftlcicntisthemainundeterminedparameterisproposedtocalculatethecriticalvelocityofrock-filledsteelcage,intheformula,cagelength,insteadofthetransformeddiameter,isemployed.Aflumetestisalsocarriedoutforthecageintheshapeofaxaxc.Accordingtotheexperimentalresults,therelationshipbetwee

6、ntheflatnessandthestabilityofthecageunderwaterisanalyzed,andthevalueofdragcoefficientunderdiflerenttlatnessisobtained.Thegreatinfluenceoffrictioncoefficientonthestabilityofrock-fillcdsteelcageisalsodiscussed.Besides,forabetterunderstandingoftheflowfieldaroundthecage,alargeeddysimulationisemployedtos

7、imulatethethreedimensionalstructureoftheflowpassingaperviouscubic-bluntbodyatRe=2.2x104.Thus,thecharacteristicsofthiskindofflowfieldanditsconcreteefleetsareconcludedfromtheaspectsofflowfieldstructure,vortexshedding,dragcoefficientandstaticpressuredistribution.Keywords:rock-filledsteelcage；criticalve

8、locity；flumetest；flowpassingaperviouscubic-bluntbody；largeeddysimulation如何减轻截流难度，一直是水利水电工程设计、从而减轻截流难度的角度出发，钢筋石笼尤其是钢筋施工中颇受关注的问题.从增大抛投料自身稳定性石笼串在立堵截流工程中得到了越来越多的应用.收稿日期：2012-06-26;修回日期：2012-09-04.基金项目：国家“卜五”科技支撑计划项口(2008BAB29B02-2-2).作者简介：叶恩V.(1987-)，男，博士研究生,yecnli2OO5.通讯作者：周宜红，zyhwhu2OO3.钢筋石笼是将工程开挖料或现场

9、易于获得的天然中小石料装入巳编制好的钢筋笼而制成的，具有相对较强的透水性和柔性.葛洲坝、三峡以及溪洛渡等工程的截流实践表明，其稳定性优于同吨位的混凝土块体且不逊于同吨位的特大石.同时，其材料的获取和制作的过程却要比特大石和混凝土块体方便得多，因此，受到越来越多的关注和应用.但到目前为止，国内外有关钢筋石笼稳定性的研究成果很少，使得钢筋石笼(串)在工程应用中缺少相应的理论指导，给截流设计工作带来了困难.因此，对钢筋石笼稳定性机理的研究具有实际意义和学术价值.针对钢筋石笼起动流速现行研究中的不足，本文从石笼的特性出发，用实际边长代替化引直径，对单个钢筋石笼的起动流速进行了理论推导，提出了以拖拽力系

10、数为主要待定参数的钢筋石笼起动流速公式.结合。xxc形石笼的起动流速水槽试验，分析了扁度与稳定性的关系，得到了不同扁度所对应的拖拽力系数值，并指出了底部摩擦力系数，对钢筋石笼稳定性的影响.此外，由于传统试验手段很难对钢筋石笼周围流场结构进行无干扰精确观测，所以为了更好地了解钢筋石笼周围流场的特性.构建了一个透水性绕流数值模型，并以大涡模拟为计算手段，对透水性绕流场进行了三维数值模拟.通过将模拟结果与非透水性绕流场从流场结构、涡体脱落、静压强分布以及拖拽力系数等方面进行定性、定量的对比分析，得到了透水性对流场的具体影响结果，为研究钢筋石笼稳定性机理提供了参考.1研究现状目前，关于截流抛投材料的起

11、动流速研究有很多成果，也出现了许多起动流速公式，其中大多数均脱胎于经典的伊兹巴斯公式，即式中v、K、O、匕、尸分别为作用在截流材料上的来流平均速度、稳定系数、抛投块体化引直径、抛投块体容重以及水的容重.伊兹巴斯公式及其类似公式(以下简称伊类公式)中的未知自变量少，应用起来相对简单，但其将影响材料稳定性的接触面粗糙程度、拖拽力系数、上举力系数、流速分布等因素全部用稳定系数K来体现和替代，不免显得过于笼统，K值的物理意义也不够明确，最终导致其值波动较大.汪定扬根据试验和收集到的资料，发现K值的变化范围是0.42.7,这样就给在实际应用伊类公式时K值的准确选取带来了困难.为此汪定扬通过多年的试验研究

12、，利用日本学者岩恒雄一的摩擦系数公式对K值进行了分解，得到了块体稳定实用计算公式，其基本关系为v=(£广京+B(枷2g矢梧(2)式中：X、8、为待定参数，根据平堵、立堵等形式的不同由试验确定；H为截流龙门水深；H/D为相对水深；夕。为相对糙度；。为断面流速分布系数，根据不同情况由试验确定.以上所述是针对一般截流块体所提出来的，而有关钢筋石笼稳定性的研究成果却很少.2010年，李学海对钢筋石笼的稳定性进行了一系列卓有成效的试验研究，引入了泥沙研究中所常用的扁度系数2，并结合块体稳定实用计算公式，提出了钢筋石笼抗冲稳定计算公式，其基本关系为v=H+8(新叫2g笠福(3)式中：、c分别为钢

13、筋石笼长、中、短轴的长度；A、8的值分起动和止动两种情况，通过试验拟合得到.在文献3中，经过验证，式具有相当的精度.但以上研究仍然存在以下两方面不足.(1) 用化引直径。来简化与球体形状偏离较大的方块体，则化引过程(一般是等体积化引)缺乏清晰的物理和数学根据；并且，块体的摆放方式也不能得到体现，例如，长轴横水向摆放与顺水向摆放时的起动流速显然不同，而采用化引直径则认为与摆放方式无关了(因为此时不存在长、短轴了).虽然式(3)有扁度系数的矫正，而旦在应用时也能够得到工程允许的结果，但从学术研究的角度，用3个几何尺度来定义钢筋石笼的形状，则能更清楚地研究影响其稳定性的因素，有助于理解其稳定性机理.

14、(2) 式(2)和式(3)都没能很好地处理块体底部摩擦力的问题.例如，假定块体与底部接触面十分光滑(接近于没有摩擦力并忽略分子间的作用力)，那么根据经典力学，块体即使受到很小的力(接近于零)的作用也会发生起动，即起动流速接近于零.这一点从式(2)和式(3)中并不能得到体现，说明此类公式并没能全面地反映出起动流速关系.本文正是针对以上不足作为研究出发点.F,.，其计算公式为2问题分析钢筋石笼在水流的作用下主要受5个力的作用：拖拽力席、上举力弓、重力G、浮力已和阻力其受力分析如图1所示.从图中可以看出，钢筋石笼与其他截流材料相比，在受力个数及性质方面都一样.导致钢筋石笼与其他截流材料稳定性的不同主

15、要是与其他截流材料相比，钢筋石笼有3方面的差异：一是透水性不同；二是与水体接触产生的表面摩擦力不同；三是与河床接触产生的底部摩擦力不同.研究钢筋石笼的起动流速需将重点放在这3点不同上，其中前两方面的不同是通过拖拽力和上举力体现出来的.图1钢筋石笼受力分析示意Fig.lSchematicdiagramofforceanalysisofrock-filledsteelcage浸没在水中的物体在与水流发生相对运动的情况下，由于水的黏性，在其运动方向上会产生拖拽力%.拖拽力由压差阻力与e和摩擦阻力&组成，两者的比例关系取决于物体的形状.但在计算拖拽力时，并没有将两者分开来，其计算式为2式中：乌

16、为拖拽力系数；4为特征面面积.特征面的选取需根据块体的形状而定，分为下列3种情况.(1) 迎水而.适合迎水面相对较大、短粗类钝体，例如球体、圆柱体、导弹、鱼雷等形体.(2) 俯视面.即为来流方向上的投影面，适合扁平形体，例如机翼、平板等.(3) 湿面.主要适用于船只等浮体.此处，研究钢筋石笼时，应该选择第1种，即选择迎水而为特征面来计算拖拽力.钢筋石笼在流水中，由于顶部和底部的流速不相等，根据伯努利方程，会产生压力差，表现为上举力尤=。/：4式中：q为上举力系数；4为钢筋石笼在来流方向的投影面积.上述Cd和CL并不是常数，而是根据块体形状、雷诺数、攻角等的不同而发生变化.对于它们之间变化的规律

17、，有很多经验性的公式和成果，但还远远不够，不能给出成熟的计算公式.对于不同的问题，往往还是通过试验，根据式和式反求Cd和CL.针对本文要研究的六面体类的块体，乌只和块体形状有关，q只和来流雷诺数有关.以滑动为条件，列出钢筋石笼起动时的力学平衡方程，即=甲，可得式中/为底部摩擦力系数.如上所述，钢筋石笼与相同形状的六面体混凝土块体相比有3方面的差别，前两方面的差别在式(6)中是通过Cd和q体现出来的.所以当底部情况一定时，研究钢筋石笼起动流速的关键在于确定式(6)中的Cd和q.当然，不同的钢筋石笼还涉及到密实度的问题，这在式中口J以通过/；来体现.3模型试验及分析钢筋石笼起动流速试验是在1:50

18、的水槽中进行的.水槽长22m,宽32cm,深30cm,槽身、槽底均采用有机玻璃制造.本试验以1:50的比尺、按不同的质量、分不同形状制作了试验所需的钢筋石笼，具体规格见表1.其中，由于工程中用立方体和长条形(包括扁形)的钢筋石笼居多，所以与韩其为【可在泥沙研究中所定义的扁度系数不同，这里重新定义了扁(a、b、c所示尺寸见图(a、b、c所示尺寸见图1).另外，为了排除迎水面摆放形式(即将。边垂直于底板放置还是平行于底板放置)给问题分析所带来的影响，本次试验只制作了axaxc形(即的钢筋石笼，后文得到的上与&的关系也仅限于oxoxc形的钢筋石笼.起动试验程序：预先将试验钢筋石笼放置在试验段

19、底板巳为减小位置随机性所带来的影响，每次将石笼放置在同一位置上；然后调节上游闸门的开度.直至石笼开始滑动(滑动速度较慢)为止，取出石笼，用小旋桨流速仪测量钢筋石笼迎水面中心点处的流速，每次起动后记录15次流速，分析时取算数平均.为降低脉动随机性带来的影响，以2min为调节上游开度的时间间隔.需要说明的是，由于本文研究的流速不是沿垂线的平均流速，并且不用像研究泥沙细颗粒时需要考虑薄膜水附加下压力I"),所以在这里不用考虑每次试验的水深.表1试验钢筋石笼规格Tab.lSizeofrock-filledsteelcagesfortest&aCrj(tm')质量原型/m模型/

20、cm原型/m模型/cm原型/t模型/g1.01.2502.501.252.502.0484.0321.51.1002.201.653.302.0484.0322.01.0002.002.004.002.0484.0322.50.9251.852.304.602.0484.0323.00.8751.752.605.202.0484.0321.01.7503.501.753.501.91310.0851.51.5003.002.254.501.91310.0852.01.4002.802.755.501.91310.0852.51.3002.603.206.401.91310.0853.01.20

21、02.403.657301.91310.0851.02.3004.602.304.601.79821.81751.52.0004.003.006.001.79821.81752.01.8253.653.657.301.7982L81752.51.6903.384.258.501.79821.81753.01.6003.204.809.601.79821.81751.03.0006.003.006.001.76947.73821.52.6205.243.937.861.76947.73822.02.3804.764.759.501.76947.73822.52.2104.425.5011.001

22、.76947.73823.02.0804.166.2412.481.76947.7382经过多组反复试验，经过算数平均整理过的试验结果见表2和图2.表2流速试验结果Tab.2Resultsofvelocitytests质量九速度/(m广)石=1.00=1.504=2.004=2.504=3.004.02.743.353.534.403.6610.03.203.854.204.983.9521.83.334.264.655.254.5047.73.654.404.905.884.95从图2中可以发现的值从1.0到2.5变化时，钢筋石笼的起动流速有逐渐增大的趋势，并且在4=2.5时起动流速达到最大

23、.究其原因，从拖拽力来看，随着&值的增大，迎水面的面积在减小，假定拖拽力系数不变，则拖拽力越小.起动流速就越大.但在4=3.0时有所减小，则很口J能是由于底部面积增大，上举力的影响超过了拖拽力，导致了起动流速的减小.-o-X,=1.0-A-=1.54=2.0-O-丈,=2.5图2不同人对应的质和速度关系Fig.2Relationshipbetweenmassandvelocityunderdifferent&由于本次试验条件的限制，缺乏单独针对Cd和/而进行的试验，并旦.考虑到石笼的透水性，在底部面积不大时，上举力的影响远小于拖拽力对稳定性的影响，所以为了分析方便，将上举力和拖

24、拽力的影响都通过拖拽力系数来体现，为此,式(6)简化为射=低施匕-尸)(7)式(7)就是钢筋石笼起动流速简化计算公式.式(7)从参数来看，仿佛与石笼的迎水面边长。、力没关系，其实，当c一定时，3对石笼稳定性的影响可以通过Cd来体现，因为正如上文所述，CD是形状的函数.需要说明的是，对于钢筋石笼来讲，其匕(体现了透水性强弱)也会对Cd产生影响，但考虑到本次试验钢筋石笼(实际工程中也一样)的。变化不大，对Cd产生的影响有限，所以本文在分析Cd时，就忽略了匕的变化对Cd产生的影响.根据试验所得流速数据以及实际测得的/值(/=0.21),代入式(7),就可以得到XQXC形钢筋石笼不同人值所对应的拖拽力

25、系数Cd(见表3).表3不同人对应的Cd值Tab.3ValuesofCDunderdifferent44Cd1.00.681.50.562.00.572.50.493.00.79为了验证式(7)在钢筋石笼为axoxc形时的准确性，现以李学海在光滑水泥面上所进行的立面体钢筋石笼起动流速试验数据为基础，对&=1.0情况下(缺少其他情况资料)这两次试验的数据和由公式计算出的数值进行了对比，结果见图3,其中，根据对钢筋石笼与光滑水泥面的摩擦系数所进行的试验，摩擦系数取为0.48.数据Fig.3Relationshipbetweenmassweightandvelocityaccord-ingt

26、ocomputationandtestunderdifferentfwhen人=1.0由图3可知，本文公式其有一定的精度和适用性.另外可以看出，摩擦系数，对钢筋石笼的稳定性有很大的影响.例如当/=0.21时，47.7t的钢筋石笼的起动流速为3.65m/s,而当f=0.48时，2t的钢筋石笼就已经可以抵抗4m/s的流速了.需要说明的是，底部摩擦系数对稳定性的这种影响不是钢筋石笼所特有的，对所有块体的稳定性研究都需十分重视摩擦系数的影响.现阶段由于缺乏对摩擦系数实质性的研究(事实上，关于摩擦力的本质目前尚未有定论)，所以在实际应用中还不能够由函数关系来对摩擦系数做出确定，往往还需要用试验来确定，有

27、时英至不直接考虑摩擦系数的影响.因此，怎样实现和评估试验环境中的摩擦系数与实际工程中的摩擦系数之间的转换值得深入研究，否则，由于摩擦系数的重要影响，当把室内试验所获得的结论向实际情况转化的过程中就会存在较大的差别，从而使室内试验失去了其指导实际工程的作用.4数值模拟及分析用传统试验手段很难对精细流场结构进行清晰的观察，并且任何接触性测量都不能对反映流场性质的物理量:进行无干扰地测量，所以本文以发展较快的CFD技术为研究手段，建立了一个透水性方钝体绕流模型，以更好地观察钢筋石笼周围流场的特性，继而从另一个角度来帮助理解透水性对钢筋石笼稳定性的影响.4.1模型建立4.1.1湍流模型选取钝体绕流，特

28、别是像方钝体这样带有锐缘(sharpedge)的绕流，周围流场结构相当复杂，其中涉及到碰撞、分离、回流以及非定常涡体脱落等复杂的湍流现象，具有高度的三维特性.目前，在计算流体力学领域，一般有3种方法用于模拟湍流：直接数值模拟(DNS)、雷诺平均法(RANS)和大涡模拟(LES).DNS在计算复杂的湍流时，必须采用很小的时间和空间步长，网格划分需要达到科尔莫戈洛夫微尺度(Kolmogorovmicroscale)的精度.据计算叫当雷诺数达到1x106时，应用DNS计算三维问题需要的网格节点数约为1x10”，所以，用DNS模拟湍流汁算代价过高，很难实现.RANS模型中，主要包括标准k-e模型、修正

29、Ay模型以及雷诺应力(RSM)模型.20世纪80年代,Mochida等°】首次发现了标准卜£模型在模拟像方钝体绕流这样带有碰撞现象的三维流动时，在碰撞区域会产生湍动能&过大的失真现象，如图4所示.之后,Rodi(m通过与试验对比,发现标准k-E模型对方钝体顶部流动的模拟与试验结果有较大的出入；之后，经过修正的*七模型，例如Kato等提出的lkA-e模型，由于较好地克服了湍动能过大现象,其模拟结果相对于标准k-e模型来讲，与试验结果符合得较好.然而，对于方钝体后部回流区的模拟，所有RANS模型的模拟结果都存在回流区长度过长而失真的现象，而大涡模拟则能够更好地处理以上问

30、题.所以本文以大涡模拟方法为研究手段(并不对大涡模拟数值方法本身进行讨论)对问题展开研究.(a)风洞试验(b)标准农£模型图4标准A-e模型模拟时产生的湍动能过大现象(取中垂面)Fig.4Over-productionofturbulentenergybythestandardk-emodel(centersectionofcube)4.1.2计算条件本文来流雷诺数屁=2.2x10。(依据方钝体边长和进口平均流速得到)，参考前人计算经验，规定本文模拟的计算条件，见表4.文中采用标准Smagorinsky亚格子格式，G取条件设定值计算区域进口条件出口条件顶部和侧边条件边壁条件网格划分上

31、游3.5.下游18H,侧边距离6H,高度2H（见图5）采用均匀来流.速度进11边界条件满足零梯度条件采用周期边界条件栗用无滑移边喂条件采用结构何格.最小间隔0.05.共I25（x）x59（v）x30（2）个缁元体（见图6）表4计算条件Tab.4Computationalconditions0.12,用二阶中心差分离散空间导数，用Adams-Bashfbrth格式离散对流项,Crank-Nicolson格式离散扩散项.采用变时间步长，'，迭代收敛性较好且趋于密定时，采用的时间步长&=0.0002s.文中进行了5x105个时间步的计算，以达到稳定的周期条件，随后又进一步进行了5x1

32、05个时间步的计算以便在做数理统计处理时能得到较为可靠的平均值.为了使方钝体具有透水性，在方钝体迎水面、背水面、侧面及顶面开【1,具体尺寸见图5,计算域网格如图6所示.图5计算模型示意Fig.5Aschematicdiagramofcomputationalmodel图6计算域网格划分Fig.6Meshesofcomputationaldomain4.2计算结果及分析4.2.1流场结构观察图7显示了在x=6.5H的截面上，透水性和非透水性绕流在不同时刻的速度场云图.从图中可以看出，无论是透水性还是非透水性绕流.在方钝体后都会产生涡体脱落现象.但由于透水性的影响,与非透水性绕流相比.在此横而上透

33、水性绕流场主要有4个方面的不同：方钝体后部涡体产生能力强,涡体形态更加丰富，流态更加发杂；透水体顶部没有产生明显的回流;绕流场最大流速较非透水性绕流下降r约12.5%;涡体脱落过程没有明显的周期性.（a）IE透水性/=60.5s/=6l.0s/=6l.5s/=62.0s/=62.5s1=63.0s速度/(m-s')3.23.02.82.62.42.22.01.81.61.41.21.00.80.60.40.2速度/(ms)2.82.62.42.22.01.81.61.41.21.00.80.60.40.2(b)透水性图7x=6.5H截面各时刻非透水性和透水性绕流速度场云图Fig.7Co

34、ntoursofvelocityonthex=6.5A/planeofflowpassinganon-perviousbodyandaperviousbodyatdifferenttimes以上特性是基于对x=6.5H截面上各时刻的速度场云图进行观察而得出的，但绕流场具有高度的三维特性，进一步分析则需要更多不同的观察角度.图8显示了在t=60.5s以及/=6l.0s时刻各截面的速度场云图.从图中凹以发现，以x=6.5/7截面为界，方体两侧流场结构各自相似.涡体的产生能力，在/=60.5s时刻，从x=7.0/7截面至x=6.0H截面呈递增趋势，而在/=61.0s时刻呈递减趋势.在其他时刻，这种递

35、增、递减现象交替出现.呈现出一定的周期性.此种现象，很可能是由于当流体经过透水体两侧时产生的涡体不对称脱落所造成的（类似于卡口涡街）.4.2.2静压强分布从方.钝体周围的静压强方面，透水性绕流和非透水性绕流也有着显著的差别，如图9所示.从图中可以发现：在迎水面的透水部位，静压强较非透水体相同部位要小得多；在顶部，透水体净压强分布比较平稳，儿乎没有波动，且其绝对值也要小于非透水体相应部位的静压强值；另外，透水体下游侧所受静压x=7.0Hx=6.87/x=65Hx=6.2Hx=6.0Wx=7.0Hx=6.87/x=65Hx=6.2Hx=6.0W.80.4.2.0.8.6.4.2.0.8.6.4.2

36、2.2.2.2.2.1.Ll,Ll.o.o.o.o.强，除了从过渡到CQ时，在C点至透水部位有一个较大的波动外，其余部位压强分布较为平稳，同时，从图7和图9口J以发现，。点附近也正是透水体背水面产生涡体脱落以及流体分离的主要部位.x=7.0Wx=6.8/x=6.5/fx=6.2Hx=6.0Hx=7.0Wx=6.8/x=6.5/fx=6.2Hx=6.0H(b)/-61.0s图8纵向各截面速度场云图Fig.8Contoursofvelocityinseveralsections通过以上对方钝体边壁t：静压强分布的分析，可以推断出透水体受力，无论是压差阻力还是上举力，都较非透水体要小.这也从一定程度

37、上解释了在水流作用下，透水性块体的稳定性优于相同形状相同质量的非透水性块体.-o-非透水A-透水CD位置图9沿x=6.5H截面轮廓线平均静压强分布曲线Fig.9Distributionoftime-averagedvalueofstaticpressurealongtheoutlineofx=6.5Hplane4.2.3拖拽力系数计算拖拽力系数C。定义为%(8)(8)扣4A=-pcosedj式中：扁和日分别为方钝体表面微元面积必上的切应力和压应力；。为表面微元上的法线与流速方向的夹角.根据式(8)结合数值模拟可提取的结果，可以算出在本文计算条件卜.，透水体的拖拽力系数约为0.63,非透水体的拖

38、拽力系数约为1.02.其中，透水情况下的拖拽力系数与试验部分相同情况下所测得的拖拽力系数值0.68(见表3)符合得较好.5结论本文用石笼实际边长代替化引直径，提出了钢筋石笼起动流速简化公式(式(7).式中以拖拽力系数为主要待定参数，相对于稳定系数K,选取更加方便且有依据(根据扁度系数选取).通过将公式应用于不同工况，验证了公式的准确性和适用性.通过对xaxc形石笼进行水槽试验，以及对透水性绕流场的三维数值模拟，本文还得出了以下结论：(1) axaxc形石笼稳定性随&值的增加先提高后减小，并在义=2.5时达到最稳状态；(2) 透水体顶部没有明显回流，流场最大流速较非透水性绕流流场F降了约

39、12.5%;涡体主要在下游网流区尾部产生且脱落过程随机性较强,没有明显的周期性，而两侧的涡体交替产生，呈现出一定的周期性；(3) 通过压强分析.发现透水体所受压差阻力和上举力都要远小于非透水体；(4) 由数值模拟计算出透水体的拖拽力系数为0.63,与试验结果相差不大.限于试验条件，本文没有将上举力系数分离出来进行研究，并且部分结论也只适用于单个石笼的特定形状，更夏杂的情况可做进一步研究.另外，摩擦系数对截流块体稳定性的研究影响显著，值得引起注意.参考文献：1 肖焕椎.施工水力学M.北京：水利电力出版社，1992.XiaoHuanxiong.HydraulicsinConstruction.Be

40、ijing：ChinaWaterpowerPress,1992(inChinese).2 汪定扬.立堵截流实用水力学计算JL水利学报，1983(9)：11-19.WangDingyang.PracticalhydrauliccalculationinenddumpriverclosureJ.JournalofHydraulicEngineering,1983(9):ll-19(inChinese).3 李学海.深段覆盖层河床截流若干关键技术问题研究D.武汉：武汉大学水资源与水电工程科学国家重点实验室.2010.LiXuchai.SeveralKeyTechnologyStudyfbrDeep-

41、ThickCoveringLayerRiverbedClosureD.Wuhan:StateKeyLaboratoryofWaterResourcesandHydropowerEngineeringScience,WuhanUniversity,2010(inChinese).4 韩其为.川江卵石形态与堆积特征调查LG/全国重点水库泥沙观测研究协作组.推移质泥沙测验技术文件汇编.1980.HanQiwei.ResearchontheshapeandthecharacteristicofgravelinChinaGJ/CooperationGroupofResearchontheSediment

42、inImportantReservoirofChina.TechnicalTestDocumentationofBedloadSediment.1980(inChinese).5WhiteFM.FluidMechanicsM.NewYork：McGraw-HillBookCompanies,1994.6EvettJB.FundamentalsofFluidMechanicsM.NewYork：McGraw-HillBookCompanies,1987.7 窦国仁论泥沙起动流速J.水利学报，1960：44-60.DouGuoren.StudiesonthecriticalvelocityofsedimentJJ.JournalofHydraulicEngineering,1960(4):44-60(inChinese).8 张瑞瑾.谢鉴衡，陈文彪.河流动力学M.武汉：武汉大学出版社，2007.ZhangRuijin,XieJianhe