碎屑流冲击桥墩的最大冲击力影响因素试验研究.pdf

第 63 卷 第 8 期2023 年8 月铁道建筑Railway EngineeringVol.63 No.8August 2023文章编号：10031995（2023）08007606碎屑流冲击桥墩的最大冲击力影响因素试验研究陈勇 姚昌荣 朱永发 周迅 安浩然西南交通大学 土木工程学院，成都 610031摘要 通过物理模型试验探究碎屑流对桥墩的冲击效应，研究碎屑流冲击距离、体量对桥墩最大冲击力的影响规律。试验再现了碎屑流起动加速、持续冲击、稳定堆积三阶段过程。试验结果表明：最大冲击力随冲击距离的增加呈指数函数负相关，随碎屑流体量的增加而呈指数函数正相关。对比冲击力试验值与半经验理论公式的预测值发现，基于流体动力学模型的理论公式能够较好地预测碎屑流冲击桥墩时的最大冲击力，理论公式预测值比试验值平均高18%，最大预测误差23%。研究结果可为山区桥梁碎屑流灾害防治提供参考。关键词 山区桥梁；桥墩；影响因素；试验研究；冲击力；碎屑流；运动过程中图分类号 U445.7 文献标识码 A DOI：10.3969/j.issn.10031995.2023.08.15引用格式：陈勇，姚昌荣，朱永发，等.碎屑流冲击桥墩的最大冲击力影响因素试验研究 J.铁道建筑，2023，63（8）：7681.碎屑流是指在崩塌或滑坡的运动过程中形成的碎屑流体，其质量大、速度快，具有很强的破坏性，严重威胁着交通基础设施和人民生命财产的安全1-2。桥梁是交通网络中的重要节点，破坏后抢修难度大，造成的损失十分严重。近年来，碎屑流冲毁桥梁的事故屡有发生。2016年，雅康高速公路天全段部分山体滑坡，碎屑流冲毁了在建的8根桥墩；2020年，雅西高速“920”高位山体崩塌引发碎屑流，导致姚河坝特大桥损毁。不同于单个滚石作用下的冲击，碎屑流冲击作用机理更复杂，且具有摩擦碰撞、时滞效应、分布作用等特性。冲击力是冲击效应最直接的衡量指标3。实际地质灾害中，碎屑流的冲击工况复杂多变，在设计过程中须定性分析。因此，有必要对碎屑流冲击桥墩的演化过程及最大冲击力影响因素进行研究，为工程设计提供参考。关于碎屑流冲击结构物的动力学研究主要包括理论方法、数值模拟方法和模型试验方法。理论研究方面，碎屑流细颗粒具有高速远程特性，其运动特性与泥石流浆体类似4。关于浆体计算模型，大多研究主要基于动力学模型变形，通过推导得到不同适用条件的半经验理论计算模型。Armanini等5提出适用于碎屑流冲击拦挡结构后发生反弹的冲击力计算模型。Delannay 等6建立了碎屑流处于过渡态时，由重力分力和惯性力组成的冲击力计算模型。此外，中国台湾、加拿大、日本等许多设计规范都是基于动力学模型进行制定的。由于公式中修正系数的变化范围较大，对计算结果会造成一定误差。数值模拟方面，碎屑流是由大量的离散颗粒组成的集合体，离散元（Discrete Element Method，DEM）作为一种处理非连续介质的模拟方法，已被诸多学者用于碎屑流研究。张睿骁等7采用DEM方法，研究了刚性挡板距离和坡度对碎屑流冲击力的影响规律。Liu 等8采用离散元和有限元耦合的数值方法研究了碎屑流对柔性障碍物的动力学作用。物理模型试验的试验条件可控性较高，是研究碎屑流的重要方法之一，国内外学者开展了大量碎屑流冲击结构物的模型试验。Jiang等9进行了碎屑流冲击挡墙模型试验，基于试验结果和已有冲击力模型，重点关注碎屑流的法向冲击力及作用点，提出了一种新的半经验冲击力模型。Wang等10通过物理模型试验研究了挡墙结构对碎屑流速度和堆积形式的影响。眭静等11根据碎屑流冲击刚性挡墙的试验结果，提出了一种冲击力力学模型，用于计算法向、切向冲击力及冲击力作用高度。收稿日期：20230404；修回日期：20230615基金项目：国家自然科学基金（51478400）；四川省科技计划项目（2021YJ0052）第一作者：陈勇（1999 ），男，硕士研究生。E-mail：通信作者：姚昌荣（1974），男，副教授，博士。E-mail：第 8 期陈勇等：碎屑流冲击桥墩的最大冲击力影响因素试验研究冷晓玉等12开展碎屑流模型试验，研究了冲击坡度和颗粒级配对碎屑流运动过程的影响规律。Kim 等13进行了碎屑流冲击圆柱形挡流板的模型试验，研究了挡流板高度和行数对碎屑流能量耗散的影响规律。Ng等14开展了碎屑流冲击双刚性拦挡结构的模型试验，研究了拦挡结构高度和距离对碎屑流冲击作用的影响。综上，已有研究主要集中于碎屑流对挡墙、拦挡结构等构筑物的冲击行为，关于碎屑流冲击桥墩结构的研究相对较少，缺乏相关的设计指南，较少关注碎屑流冲击桥墩时的最大冲击力。因此，本文通过室内模型试验模拟碎屑流冲击桥墩的运动过程，研究碎屑流对桥墩的冲击效应，以及碎屑流冲击距离、体量对桥墩最大冲击力的影响规律，为山区桥梁碎屑流灾害防治提供一定参考。1 试验概况 1.1试验装置试验装置主要包括斜槽、物料区、反力架、废料箱等部分，见图1。斜槽尺寸为2.0 m（长）0.3 m（宽）0.5 m（高），坡度为20。斜槽侧面采用透明亚克力板材质，表面覆有网格尺寸为2 cm的网格纸，以便统计碎屑流的运动距离。将前后挡板粘贴于斜槽中形成物料区，临时固定前挡板以便抽离，使碎屑流在重力作用下滑动。反力架由铝合金焊接而成，用于支撑固定六分力传感器及桥墩模型，其底部嵌固在地上。1.2桥墩模型及测量系统以实际工程山区高架桥为原型，跨径为25 m，桥墩墩高为10 m，直径为1.5 m，按1 50的比例缩尺制作双柱式桥墩模型。桥墩模型由树脂材料采用3D打印制成，直径为0.03 m。定义来流方向一侧的墩柱称为近柱，另一墩柱为远柱。测量系统包括六维力传感器和高速相机。两个六维力传感器与两个墩柱之间通过连接装置用螺栓连接，六维力传感器用螺丝固定于反力架上，见图2。传感器采样频率为50 Hz。在桥墩顶部、侧面布置高速相机，用于记录碎屑流的运动形态。1.3试验材料及方案设计自然界中的碎屑流体通常由形状大小各异的岩石块组成，且粒径级配分布十分广泛，若在现场测量粒径分布、形状及其他参数，或在室内模型试验中再现真实的级配比例及形状是非常困难的15。因此，采用模型试验或数值模拟的方式研究碎屑流体的一般冲击特性时，主要采用粒径分布均匀、尺寸范围较窄的碎石来模拟碎屑流体颗粒16-17。已有研究表明颗粒平均粒径与水槽宽度的比值在1/90 1/20时，可忽略水槽宽度对碎屑流体自由滑动限制的影响15。考虑试验中桥墩模型的直径及斜槽装置的宽度，试验材料选择无黏性的砾石。砾石初始粒径分布均匀，粒径分布范围为 8 16 mm，中值粒径（D50）与水槽宽度比值为1/23，满足试验要求。砾石材料参数见表1。本文以碎屑流冲击距离（L）、体量（V）为控制参数，通过改变初始堆积体高度（H）来实现 V 的变化。采用控制变量法研究碎屑流冲击距离、体量对其运动形态及冲击桥墩时的冲击力沿程演化规律。试验工况见表2，每组工况重复三次试验。图1碎屑流物理模型试验装置图2桥墩模型传感器布置（单位：mm）表1砾石材料参数参数干密度/(kgm-3)D50/mm不均匀系数休止角/（）数值1 690131.416.8表2试验工况工况A1A2A3B1B2B3C1C2C3L/m0.400.650.90H/m0.300.400.500.300.400.500.300.400.50V/m0.022 50.030 00.037 50.022 50.030 00.037 50.022 50.030 00.037 577铁道建筑第 63 卷2 碎屑流运动形态分析 2.1运动过程以试验组B2为例分析碎屑流运动过程。为便于观察比较，将图片旋转20并按水平排列，见图3（a）。可知：t=0.24 s时，前挡板完全提起后，位于碎屑流体底部的颗粒率先下滑移动，带动前缘颗粒移动的同时上部颗粒迅速下落；t=0.60 s时，前缘颗粒开始冲击近柱桥墩；t=0.77 s时，颗粒持续冲击桥墩，出现颗粒碰撞飞溅的现象，少量颗粒开始冲击远柱桥墩；t=1.60 s时，由于颗粒之间、颗粒与斜槽之间的摩擦作用，碎屑流体稳定堆积在斜槽与桥墩之间。综上，碎屑流运动过程可分为起动加速、持续冲击、稳定堆积三阶段过程。其他工况的试验现象与试验组B2类似，但有些许不同，见图3（b）。可知，当冲击距离增加后，碎屑流颗粒冲击桥墩时的碰撞飞溅程度较弱，颗粒飞溅高度更低。这主要是由于碎屑流冲击过程中，重力势能转换为动能和摩擦内能，而摩擦内能是最主要的耗能方式17。随着冲击距离的增加，颗粒之间、颗粒与斜槽之间由于摩擦消耗的能量更多，颗粒的冲击速度变小，冲击作用变弱。2.2运动速度定义前挡板提起瞬时为初始时刻（t=0），碎屑流前端颗粒自初始时刻到冲击近柱时刻的平均速度为v 1、冲击近柱后到滑出斜槽的平均速度为v 2。两段平均速度为对应的运动距离和运动时间之比。运动距离借助斜槽侧面网格纸进行统计，运动时间通过回放试验录像获得，绘制前端颗粒的平均运动速度，从而分析碎屑流体量、冲击距离对前端颗粒平均速度的影响规律，见图4。由图 4可知：当碎屑流冲击距离相同时（工况A1、A2、A3），v 1和v 2均随碎屑流体量的增大而增大。碎屑流体底部的颗粒先下滑移动，导致上部颗粒迅速下落。本文试验中碎屑流体量越大，初始堆积体的高度越高，且上部颗粒具有更大的重力势能，因此下滑过程中产生的动能越大，速度越高。当碎屑流体量相同时（A1、B1、C1），v 1随着冲击距离的增大而增大，v 2随着冲击距离的增大而减小。这主要是由于碎屑流运动过程中，重力势能转换为动能，冲击距离越远初始重力势能越大，自初始时刻到冲击桥墩时速度越大。而撞击桥墩后，因颗粒之间、颗粒与斜槽、颗粒与桥墩之间的摩擦作用，导致大部分颗粒堆积于斜槽和桥墩间，这种现象随着冲击距离的增加愈发明显，导致后部颗粒的能量不能有效传递到前端颗粒，因此前端颗粒的运动速度减小。3 冲击力特征 3.1冲击距离、体量对碎屑流冲击力的影响分析冲击力试验数据后发现近柱受到的冲击力远大于远柱，本文旨在探究碎屑流冲击桥墩时最大冲击力的影响因素规律，因此仅分析近柱的冲击力规律。不同冲击距离下碎屑流冲击力时程演化曲线，见图5。可知，冲击力可分为碎屑流体运动时的冲击力波动段和稳定堆积后的冲击力稳定段。随着碎屑流冲击距离的增加，碎屑流自起动加速到冲击桥墩的时间增加，桥墩受到的最大冲击力减小，桥墩受到碎屑流的堆积作用也变小，相同碎屑流体量下的试验工况规律相同。原因是当碎屑流冲击距离增加时，沿程颗粒之间、颗粒与滑槽之间的碰撞、摩擦所消耗的能量增加，更多的颗粒堆积在斜槽中，导致冲击桥墩时的颗粒数减少，所受冲击力峰值也减小。不同碎屑流体量下冲击力时程演化曲线见图6。可知，体量的增加对碎屑流开始冲击桥墩的时间影响不大。随着碎屑流体量的增加，桥墩所受到的最大冲图3碎屑流运动过程和颗粒飞溅效果对比图4各试验工况的前端颗粒平均速度78第 8 期陈勇等：碎屑流冲击桥墩的最大冲击力影响因素试验研究击力增加，颗粒稳定时的堆积作用越大。原因是当碎屑流体量增加时，前端颗粒开始冲击桥墩的时间不变，意味着速度变化不大，但同一时刻冲击作用于桥墩的碎屑流颗粒数量增加，导致最大冲击力和堆积作用力增加。3.2最大冲击力影响因素分析冲击距离、体量是影响碎屑流最大冲击力的重要因素，提取各工况下的最大冲击力数据进行曲线拟合，见图7。图中，R2为决定系数。可知：初始堆积体高度为 0.30、0.40、0.50 m 时，决定系数分别为0.983、0.998、0.991，回归拟合度较好。当碎屑流体量一致时，即初始堆积体高度相同时，最大冲击力随着冲击距离的增加呈指数函数负相关，与文献 18 规律基本一致。原因是当碎屑流冲击距离越大，颗粒碰撞飞溅的高度越低 参见图4（b），即冲击过程中碎屑流的运动速度越小，碎屑流颗粒之间碰撞所消耗的动能越小，导致最大冲击力下降的幅度较小。冲击距离为 0.40、0.65、0.90 m 时，决定系数分别为 0.869、0.998、0.973，回归拟合度较好。当冲击距离相同时，最大冲击力随着碎屑流体量的增加呈指数函数正相关。原因是碎屑流体量越大，冲击过程中作用于桥墩的碎屑流颗粒也越多，且碎屑流的运动速度相对较大，从而导致最大冲击力上升的幅度较大。综上，参照文献 19 对碎屑流最大冲击力（F）与L、V之间的关系进行拟合，拟合关系式为F=0.004 4V2+22.24L2-0.772 1VL+0.825 1V-23.42L-0.193 0，R2=0.967（1）最大冲击力与冲击距离、体量之间的关系见图8。可知：碎屑流冲击距离为0.90 m，体量为0.022 5 m时，最大冲击力为2.199 N，是所有工况中的最小值；冲击距离为0.40 m，体量为0.037 5 m时，最大冲击力为18.908 N，是所有工况中的最大值。碎屑流最大冲击力与冲击距离成反比，与碎屑流体量成正比。在山区桥梁碎屑流灾害防治中，应优先考虑碎屑流体量更大且冲击距离更近的冲击工况。4 试验值与半经验公式预测值对比 文献 4 研究表明，具有高速远程特征的碎屑流细颗粒，其运动特性与泥石流“浆体”类似4。因此，采用流体动力学模型对比本文最大冲击力试验结果与常见的半经验理论公式预测结果20。模型表达式为图6不同碎屑流体量下冲击力时程演化曲线图7最大冲击力与影响因素的拟合曲线图8最大冲击力与冲击距离、体量之间的关系图5不同冲击距离下碎屑流冲击力时程演化曲线79铁道建筑第 63 卷P=muv2mu（2）式中：P、mu、vmu分别为流体冲击压力、密度、运动速度；为动力修正系数。参考T/CAGHP 0062018 泥石流灾害防治工程勘察规范（试行）21中工程截面为圆形时的情况，取1.0。选取最大冲击力对应时刻的运动速度作为vmu，获取方式参考文献 22，即通过相邻几帧图像上碎屑流前端颗粒的运动距离变化，计算出相应的前缘速度。同时考虑碎屑流冲击桥墩的作用面积，桥墩直径取0.03 m；通过回放试验录像获得各工况下碎屑流在桥墩上的作用高度（0.02 0.09 m）。半经验公式为F=muv2mubh（3）式中：b为桥墩直径，m；h为碎屑流作用高度，m。对比不同工况下桥墩最大冲击力的试验值与半经验公式预测值，见图9。可知，半经验公式预测值与试验值最大误差为23，各工况下两者的数值差异较小，可见半经验公式能够较好地预测碎屑流冲击桥墩时的最大冲击力。不同工况下的半经验公式预测值比试验值平均高18，可见采用半经验公式计算的最大冲击力偏保守。受试验条件的限制，未能模拟出碎屑流冲击桥墩过程中实际的地质环境，如本文斜槽底部平整、无明显起伏，碎屑流冲击距离、体量变化范围小等；碎屑流的颗粒粒径考虑较为单一，没有考虑含水量、块状黏性物质等，因此试验中未能体现出碎屑流运动过程中的粘黏、破碎、筛分等现象。此外，本文未考虑试验装置的尺寸效应。试验条件的限制和未考虑试验装置尺寸效应对碎屑流体的运动形态、速度、冲击力时程演化等具有重要影响，后续可开展不同尺度的模型试验进行研究。5 结论 1）试验过程中再现了碎屑流起动加速、持续冲击、稳定堆积的三阶段过程。在撞击桥墩前，碎屑流的运动速度随冲击距离和体量的增大而增大；撞击桥墩后，运动速度随冲击距离的增大而减小。2）碎屑流冲击距离、体量是影响碎屑流最大冲击力的重要因素。最大冲击力随冲击距离的增加呈指数函数负相关，最大冲击力随碎屑流体量的增加呈指数函数正相关。3）基于流体动力学模型的半经验公式能够较好地预测碎屑流冲击桥墩时的最大冲击力，预测值比试验值平均高18，最大预测误差为23。采用半经验公式预测的最大冲击力偏于保守。参考文献1 叶振南，陈宗良，高幼龙，等.甘肃陇南泻流坡滑坡-碎屑流变形特征与成因机制 J.科学技术与工程，2021，21（22）：9250-9257.2 ZHANG M，YIN Y P，WU S R，et al.Dynamics of the Niumiangou Creek Rock Avalanche Triggered by 2008 Ms 8.0 Wenchuan Earthquake，Sichuan，China J.Landslides，2011，8（3）：363-371.3 NICOLETTI P G，SORRISO-VALVO M.Geomorphic Controls of the Shape and Mobility of Rock Avalanches J.Geological Society of America Bulletin，1991，103（10）：1365-1373.4 WIEGHARDT K.Forces in Granular FlowJ.Mechanics Research Communications，1974，1（1）：3-7.5 ARMANINI A，LARCHER M，ODORIZZI M.Dynamic Impact of a Debris Flow Front Against a Vertical WallC/Proceedings of the 5th International Conference on Debris-Flow Hazards Mitigation：Mechanics，Prediction and Assess-ment.Padua：Journal of Hydraulic Research，2011：1041-1049.6 DELANNAY R，Valance A，MANGENEY A，et al.Granular and Particle-Laden Flows：from Laboratory Experiments to Field Observations J.Journal of Physics D：Applied Physics，2017，50（5）：1-39.7 张睿骁，樊晓一，姜元俊.滑坡碎屑流冲击力影响因素的离散元模拟 J.铁道建筑，2019，59（4）：108-111.8 LIU C，YU Z X，ZHAO S C.Quantifying the Impact of a Debris Avalanche Against a Flexible Barrier by Coupled DEM-FEM Analyses J.Landslides，2020，17（1）：33-47.9 JIANG Y J，FAN X Y，SU L J，et al.Experimental Validation of a New Semi-empirical Impact Force Model of the Dry Granular Flow Impact Against a Rigid Barrier J.Landslides，2021，18：1387-1402.10 WANG Z F，HE M，WANG Y，et al.Dynamic Process of Rock Avalanche and Dynamic Interaction Between Rock Avalanche and Hazard-bearing Bodies J.Arabian Journal of Geosciences，2021，14（20）：1-11.11 眭静，姜元俊，樊晓一，等.碎屑流冲击刚性挡墙的力学模型研究 J.岩石力学与工程学报，2019，38（1）：121-132.12 冷晓玉，樊晓一.颗粒级配与场地条件对滑坡运动特征的图9最大冲击力试验值与半经验公式预测值对比80第 8 期陈勇等：碎屑流冲击桥墩的最大冲击力影响因素试验研究影响模型试验研究 J.铁道建筑，2016，56（5）：139-143.13 KIM B J，YUNE C Y.Flume Investigation of Cylindrical Baffles on Landslide Debris Energy Dissipation J.Landslides，2022，19（12）：3043-3060.14 NG C W W，CHOI C E，KOO R C H，et al.Dry Granular Flow Interaction with Dual-barrier SystemsJ.Gotechnique，2018，68（5）：386-399.15 JIANG Y J，TOWHATA I.Experimental Study of Dry Granular Flow and Impact Behavior Against a Rigid Retaining WallJ.Rock Mechanics and Rock Engineering，2013，46：713-729.16 GONG S，ZHAO T，ZHAO J，et al.Discrete Element Analysis of Dry Granular Flow Impact on Slit DamsJ.Landslides，2021，18（3）：1143-1152.17 柳春，余志祥，刘宇鹏，等.碎屑流冲击棚洞的最大冲击力影响因素研究 J.振动与冲击，2020，39（2）：195-203，222.18 陈果，钮志林，樊晓一，等.高速远程滑坡沿程速度演化与冲击力分布研究以三溪村滑坡为例 J.自然灾害学报，2022，31（3）：232-241.19 陈果.滑坡碎屑流动力参数与冲击力沿程演化研究 D.绵阳：西南科技大学，2022.20 刘洋，游勇，王海帆，等.颗粒流冲击力研究现状及讨论 J.防灾减灾工程学报，2020，40（5）：714-723.21 中国地质灾害防治工程行业协会.泥石流灾害防治工程勘查规范（试行）：T/CAGHP 0062018 S.武汉：中国地质大学出版社有限责任公司，2018.22 杨海龙，樊晓一，赵运会，等.偏转角度对滑坡-碎屑流运动影响的模型试验 J.山地学报，2017，35（3）：316-322.Experimental Study on Influencing Factors of Maximum Impact Force of Debris Flow Impacting Bridge PierCHEN Yong，YAO Changrong，ZHU Yongfa，ZHOU Xun，AN HaoranSchool of Civil Engineering，Southwest Jiaotong University，Chengdu 610031，ChinaAbstract Through physical model tests，the impact effect of debris flow on piers was explored，and the influence of debris flow impact distance and volume on the maximum impact force of piers was studied.The test reproduces the three-stage process of debris flow initiation acceleration，continuous impact and stable accumulation.The results show that the maximum impact force decreases exponentially with the increase of impact distance，and increases exponentially with the increase of impact volume.The comparison between the impact force test value and the predicted value of the semi-empirical theoretical formula shows that the theoretical formula based on the hydrodynamic model can better predict the maximum impact force when the debris flow impacts the pier.The predicted value of the theoretical formula is 18%higher than the experimental value，and the maximum prediction error is 23%.The relevant research results can provide some reference for the prevention and control of debris flow disasters in mountain bridges.Key words mountain bridge；pier；influencing factor；experimental research；impact force；debris flow；movement processCitation format：CHEN Yong，YAO Changrong，ZHU Yongfa，et al.Experimental Study on Influencing Factors of Maximum Impact Force of Debris Flow Impacting Bridge Pier J.Railway Engineering，2023，63（8）：7681.（编辑：郑冰 校对：苗蕾）81