水利工程论文-充沙长管袋抛掷沉降运动的力学研究.doc

水利工程论文-充沙长管袋抛掷沉降运动的力学研究摘要：本文对充沙长管袋这种新型的筑坝、抢险整治工程材料进行了抛掷沉降运动的力学分析研究与模拟计算。根据充沙长管袋的几何、材料特性，本文分析了管袋沉降机理与沉降速度；在多个小块柔性连接体的假设下，对抛掷在河流中的管袋进行了力学分析与动力学过程的研究，建立了管袋沉降的连续体数学模型。通过合理简化，导出了实用的管袋位移与转角的计算公式。根据管袋抛掷的观测资料，确定了沉降计算的有关参数；针对黄河水流、泥沙特点及现实施工条件进行了多种管袋抛掷沉降的模拟计算，提出了适应黄河丁坝防护的管袋尺寸及位移范围值。关键词：丁坝充沙长管袋沉降速度位移与转角近年河务部门尝试采用将填充泥沙的土工织物长管袋抛掷在丁坝周围，利用其柔性、贴体、防冲的特点，抵抗洪水对丁坝坝体及附近河床的冲刷，达到抢险护坝的目的。有关柔性长管袋抛掷后的力学特性以及管袋在河流中的位移、稳定定位的确定，是河道坝岸工程设计、抢险施工技术遇到的新问题。本文拟就长管袋抛掷沉降的力学过程进行初步的理论分析与模拟计算。1管袋抛掷的力学模型1.1长管袋的水平作用力为了便于理论分析，下面将实际工程中采用的长管袋概化为直径d长度L的柔性圆柱体。按水平抛掷管袋法线x与水流方向夹角为，见图1.此时洪水来流流速可以分解为垂直于管袋轴线的x方向和平行于管袋轴线的y方向的分速度x、y：x=cos,y=sin(1)水流对管袋的作用力F也可分解为垂直和平行于管袋轴线方向的分力Fx、Fy，分别进行研究。在河流中沿管袋轴线的流速分布不均匀，所以取单位长管段研究垂直管袋轴线的水流作用力：fx=CdmAd(x)2/2g(2)式中：fl洪水对单位长管段的水平作用力；Cd管袋绕流阻力系数；m浑水容重；Ad单位长管段的绕流特征面积；x水流与管袋间的相对速度。由式(2)，可将垂直于管袋轴线方向合力表示为：(3)1.2管袋在水中的沉降特性充沙长管袋在河流中的水力特性是指长管袋在水中的沉降规律，常以长管袋的沉降速度表征。抛掷在水中的管袋受到重力作用下沉，初始向下作加速运动，同时也产生绕流阻力。当管袋的有效重力与管袋所受的绕流阻力相平衡时，管袋将以匀速继续下沉。管袋在水体中匀速下沉的速度称为管袋沉速，用表示。管袋沉速的大小主要与管袋的尺寸、泥沙密实度和绕流状态有关。实际河流洪水流速较高、Re较大，观测表明管袋沉降时引起周围水体强烈的扰动，属于强紊动沉降状态。长管袋形状不很规则，沉降规律比较复杂。如果设管袋与河水的容重分别为s和m，则管袋在水中沉降时受到的有效重力为：W=(s-m)d2/4L(4)管袋下沉时受到水体的绕流阻力为：Fd=CdmAd2/2g=CdmdL2/2g(5)式中：管袋沉速；Cd垂向绕流阻力系数；Ad绕流特征面积。当管袋匀速下沉时，必定有W=Fd，由上述力学分析则可导出管袋沉降速度：(6)上式是管袋沉降规律的一般形式，式中绕流阻力系数与管袋沉降状态有关。实际管袋的形状不规则，在沉降时还受水流紊动的影响，因此实际管袋在河流中的沉速应写成如下形式：(7)式中：K与管袋形状、尺寸有关的沉降系数。根据相关实测资料分析，其值一般可取0.821.15之间。浑水容重m与含沙量有关：m=+0.623S；S浑水含沙量，单位为kg/m3.长管袋抛投在水中沉降，沉降时间与抛投点处的水深H和管袋的沉速有关。如果将动水的影响以系数c修正，则长管袋在水中的沉降时间t可表示为：t=cH/(8)图1丁坝前抛掷管袋示意图2水平面内管袋抛掷受力示意1.3管袋沉降的连续体数学模型充沙长管袋可视为多个小块体的柔性连接体，假定各小块体为刚体，它们的运动由作用在其上的合力及合力矩所决定，或平移或围绕其形心的转动。从简化考虑，分别研究沿x-y和y-z两个平面运动，用增量法逐步计算出小块体的位移。一单个块体在变化的外力Fk作用下产生的运动，可用牛顿第二运动定律描述：(9)式中：小块体的形心运动速度；Fk小块体形心处所受到的合力，m小块体的质量；t时间。上式左边采用中心差分格式，在时间t可表达为：(10)将以上两式整理为：(11)上式中半个时步点的速度可以用小块体形心位移u的形式写出：(12)由于力的产生依赖于位移，所以力位移的计算在同一时步内同时进行。块体在多个力及重力作用下，其运动速度方程变为：(13)(14)式中：块体对于其形心的角速度；I块体的惯性矩；gi块体的重力加速度分量；i块体形心的速度分量；(i=x,y).根据以上两式得出新速度，可以由下式确定出块体的新位置。(15)(16)式中：i块体绕其形心的转动量；xi块体形心的坐标分量；(i=x,y).这样每一迭代过程中产生块体的新位置都导致新的接触力。合力与合力矩产生线加速度和角加速度，块体的速度和位移可由整个时步增量的积分获得。用这种方法可以计算出管袋每一时刻的位置，比较精确地描述管袋的实际运动过程。1.4管袋沉降的计算方法这里主要考虑在水平面内，确定管袋顺丁坝轴线抛掷的位移及转角。为了简化计算，假定管袋运动速度很小，水流对管袋的作用力因管袋的移动仅改变很小，这样可将水流与管袋的相对速度用水流速度代替。为了正确反映水流对管袋的作用力，应该对流速加以修正：=C(17)式中：C流速修正系数。修正系数C随水流速度的不同而变化。将管袋的运动分为水平运动和竖直向下运动两部分。对于水平运动，抛掷管袋受力示意如图2所示。其中Fx为作用于与管袋轴线垂直方向分力，是Fx偏离管袋重心的系数。根据牛顿第二运动定律：ax=Fx/m(18)考虑到管袋质量m=s/4gd2L,有单位长管袋绕流面积为：Ad=d1=d(19)将式(1)、(2)、(3)、(17)、(19)代入式(18)得：(20)管袋所受到的合力矩M为：M=LFx，则管袋的转动惯量为：I=m/48(3d2+4L2)(21)而管袋绕其形心的角加速度为：=M/I(22)将式(1)、(2)、(3)、(17)、(19)、(21)代入式(22)得：(23)管袋重心的水平位移为：u=1/2axt2(24)将式(20)、(8)、(7)代入式(24)得：(25)因为相对管袋原始位置的转角为：=1/2t2(26)将式(23)、(8)、(7)代入式(26)得：(27)式中为：(28)若将管袋顺水流方向抛掷(管袋轴线平行水流方向)，推导过程与上述类似故从略，则位移计算公式为：u=CdC2c2/4K2H22/Lgd2m/s(s-m)(29)2管袋抛掷沉降的模拟计算2.1计算参数的确定在管袋运动理论分析中出现的有关参数，如绕流阻力系数Cd、管袋形状影响系数K、沉降综合影响系数c、相对运动影响的流速修正系数Cv等，都是影响管袋沉降物理过程的重要因素。但这些系数目前都缺乏可供计算的公式，因此可靠的参数只有依靠大量管袋抛掷实测资料来推求、校验。目前在黄河堤防工程中如丁坝防护，已有使用管袋的实例。根据收集到的管袋抛掷现场观测资料，一般水深H在2m14m之间，流速v在0.7m/s1.5m/s之间，大河含沙量S=60kg/m3左右；管袋直径D=0.8m1.0m，管袋长L=7m10m，管袋填充度即管袋内泥浆或土体容重s=1200kg/m32100kg/m3.对这些实测资料进行分析，初步率定了有关参数见表1、表2.应用有关参数进行了管袋位移及转角的模拟计算，计算值基本与实测值吻合，见图3.这表明所选参数基本合理，可以用来预测典型设计条件下的管袋沉降问题。表1顺水流方向抛掷管袋参数参数名称CdKcCv参数数值0.651.13.60.550.65表2垂