非均一性斜坡落石运动轨迹

现有的落石运动轨迹研究均将斜坡视为半无限均匀的坡面,而实际的工况下,坡面由不均一材料组成。根据落石不同的运动模式,将坡面简化为多层岩土体材料的结构形式,在滚动阶段,基于Hertz接触理论得到落石在不同坡面条件下运动特征,并给出坡面分层条件下落石切向摩擦系数的计算公式,在碰撞阶段,基于能量守恒定理,采用准静态接触力学理论得到法向恢复系数公式,理论推导得到任意层厚岩土体材料坡面条件下落石法向和切向恢复系数的解析解,然后求解得到碰撞后落石运动参数,并根据碰撞后落石的速度及回弹量关系,给出碰撞后落石运动模式转换的判别条件,最后结合运动学定理,得到3种运动模式任意分层坡面条件落石的运动轨迹,将理论公式应用于某山区落石运动轨迹预测,验证了理论公式的适用性和有效性。     

边坡滚石运动轨迹分段循环算法

边坡滚石是在基本建设中常遇到的一种工程灾害。其防护结构的设计依据主要是边坡滚石的运动轨迹,因而要对这种灾害进行防治,首先要确定滚石的运动轨迹。根据落石运动全过程的三个阶段,即滚动(滑动)运动阶段、抛物运动阶段和与坡面碰撞阶段,利用分段循环算法分别得到了三个运动阶段的运动速度计算公式。以某岩质边坡构建相似模型,将该理论公式应用于落石运动轨迹预测,验证了公式的可行性和有效性。与现有方法相比,分段循环算法具有计算原理简单明确、易于使用等优点。该方法可用于斜坡区滚石运动特征的预测,并可作为斜坡区滚石灾害防治的依据。     

隧道仰坡落石冲击模型试验研究与机制探讨

以在建的成都－兰州（成兰）铁路某隧道出口仰坡为研究对象,研发了大型落石冲击边坡模型试验系统,该系统由模型试验台架、落石三维空间释放装置、高速摄影系统等组成,结构强度高并能满足试验过程的可视化需求。根据现场岩样实测物理力学参数,采用凡士林、硅油、水泥、细沙、重晶石粉、石膏以及滑石粉等原料,研制出现场千枚岩和砂岩的相似材料,并用于模拟边坡和落石。通过开展物理模型试验研究了近长方体落石初始滑移角度、形状、质量等3种因素对碰撞恢复系数这一落石运动轨迹模拟中关键参数的影响,同时在落石下落冲击的过程中,利用高速摄影仪和坐标卡尺记录落石完整的运动轨迹,后期通过视频处理软件分析落石与边坡碰撞前后的速度变化,根据公式计算出落石法向和切向的碰撞恢复系数。试验结果表明：随着初始滑移角度的增加,板状试件的法向、切向恢复系数分别呈减小和增大的趋势,且切向恢复系数变化更加明显;落石法向恢复系数随质量的增加呈明显减小的趋势,而切向恢复系数值变化不大;对于近长方体落石,立方体试件的法向和切向恢复系数大于厚板状试件,条形试件的两向恢复系数值受碰撞接触形态的影响,呈现两极化的趋势。最后,结合试验现象以及前人的研究成果对落石运动冲击机制进行了初步探讨。     

明棚洞落石冲击荷载计算方法研究

根据落石冲击回填土明棚洞过程中的基本力学和运动学规律,采用Laplace变换,从落石冲击半无限土体、落石冲击有限厚度土体、考虑回填土与下部结构相互作用3个方面推导了落石冲击荷载的理论计算公式,通过数值模拟和相关试验数据对理论公式进行了验证,并与目前常见落石冲击荷载计算公式进行了对比研究。研究表明：所提出的公式理论计算结果与数值模拟结果和相关试验结果均有较稳定的规律,理论计算值比数值模拟值大6%～41%,与Pichler现场试验的95%分位值差值仅在6%以内;尽管所提出的公式计算值比目前所有公式计算结果均偏大,但与数值模拟和试验数据更吻合,更能反映出真实的落石冲击荷载;所提出的计算公式可以反映缓冲层厚度和缓冲层性质的影响,缓冲层厚度越小,落石冲击荷载越大,缓冲层厚度影响系数可根据h/r查表确定;缓冲层和结构的动力相互作用导致结构承受冲击荷载相较缓冲层顶部冲击荷载增大,结构动力放大系数与缓冲层厚度相关,缓冲层厚度越大,结构动力放大系数越小。所提出的计算公式可以反映落石大小、落石形状、落石冲击能量、缓冲层厚度、缓冲层性质等因素对冲击荷载的影响,计算荷载值可直接用于明棚洞结构设计。     

边坡落石灾害的颗粒流模拟方法

落石是边坡的一种常见的动力地质灾害。基于颗粒流理论,提出了落石灾害的二维的可视化模拟方法。通过颗粒流中不同参数的ball单元的组合模拟了各种形状和特性的落石体;通过wall单元的组合模拟坡面;选择ball-ball和ball-wall的接触刚度模型、粘结模型和滑移模型等接触本构模型来反映落石运动过程中落石岩块内部各部分和落石与坡面作用行为;实现了边坡落石灾害的二维可视化模拟。模型模拟出落石的自由落体、弹跳、滑动和滚动等运动形式、运动过程岩块的崩解;得到了落石的运动速度、位移等运动学参数。工程实例初步验证了该方法的可行性、合理性。     

落石冲击混凝土棚洞力学特性研究

落石冲击棚洞结构作用过程复杂,缺乏统一的落石冲击力表达式。首先,将落石简化为刚性球体,基于Hertz接触理论,推导得到落石冲击力半正弦算法的理论表达式,考虑落石冲击下棚洞的非弹性特征,根据落石与材料碰撞过程中落石加速度曲线特征,采用函数拟合法推导得到落石法向冲击下其冲击力的理论计算方法;然后,基于ANSYS/LS-DYNA软件建立落石冲击棚洞数值计算模型,研究不同冲击速度下落石冲击棚洞动力特征;最后,与现存常见的多种方法进行对比,得出以下结论：Hertz半正弦法得到的落石冲击力远大于函数拟合法和数值法,而函数拟合法和数值法得到的落石冲击力时程曲线相接近,表明函数拟合法更能反映落石与棚洞接触碰撞动力关系;对比其他计算方法可以得到,Hertz算法适用于分析无能量损失下的弹性碰撞问题,而Logistic算法适用于材料大塑性变形的情况,弹塑性接触理论结果和动力有限元结果存在差异,而采用函数拟合推导的计算方法得到的落石最大冲击力和落石冲击作用时间与动力有限元法更接近,更能反映落石冲击棚洞动力响应特征,推导的落石冲击力计算方法可为工程实践中棚洞防护设计提供理论参考。     

鄂西山区崩塌落石运动特征及危险性分析

以远安县瓦坡崩塌区为例,通过地质调查、野外测绘、无人机航拍,建立了瓦坡崩塌区三维模型,基于Rockfall Analyst（RA）分析软件,实现了瓦坡崩塌区大量崩塌落石三维空间下运动路径、高度、能量等要素模拟,探究了崩塌区落石的三维运动特征,开展了崩塌区落石危险性评估,为崩塌落石的防治提供科学参考。研究结果表明:崩塌区危岩主要破坏形式为倾倒式,目前处于欠稳定状态;模拟落石三维运动轨迹与已有落石点基本重合,说明此次模拟结果与实际情况较为吻合;落石运动过程中以碰撞弹跳、自由飞落为主,落石主要集中在崩塌区下部冲沟及公路内,部分落石达到居民区,在崩塌落石区下部公路、冲沟及崩塌区右侧危险性较高;崩塌落石防治工程建议采用危石孤石清除+被动防护网,在公路内侧、斜坡下方分别设置5 m高和3 m高抗冲击力2 000 kJ的被动防护网,可有效拦截落石。     

北川羌族自治县都坝乡政府后山崩塌落石

针对四川省绵阳市北川羌族自治县都坝乡政府后山崩塌的落石,结合影响落石运动的各种条件因素,分析并确定落石的运动形式,选取合适的计算公式,初步得出了落石运动的速度及其对坡下危害目标的冲击能量。     

SNS被动防护系统在崩塌地质灾害治理中的应用

朝阳市位于辽宁西部低山丘陵区,近年来,崩塌落石地质灾害频繁发生,瓦房子镇团山子为典型的软硬岩接触导致的陡崖崩塌落石地质灾害。朝阳地区常用的崩塌地质灾害防护结构多为圬工结构,圬工结构存在材料用量大搬运困难、施工强度大、施工进度缓慢、原始生态环境破坏大、不美观等缺点。综合考虑瓦房子镇团山子崩塌的地质环境条件,分析崩塌地质灾害特征,通过多种治理方案比选,确定采用SNS被动防护系统治理。该系统克服了圬工结构的缺点,能够有效地拦截落石,在受力上具有以柔克刚的优点。通过对瓦房子崩塌落石运动轨迹计算,给出了一整套计算边坡落石运动轨迹的公式,为类似崩塌地质灾害治理方案选择及落石运动轨迹计算提供了参考。     

基于运动特征分析的高陡边坡崩塌落石防治措施建议

崩塌落石具有突发性、随机性、隐蔽性及致灾严重性等特征。为对其进行有效防治,以某煤矿内崩塌落石为研究对象,采用RocFall软件对已崩落石的运动轨迹进行了模拟,反演拟合了下垫面参数,在此基础对拟崩落石区域进行数值模拟和理论计算得出其运动特征,为崩塌落石防治措施选择提供了合理依据。     

Rockfall hazard and risk analysis for Monte da Lua,Sintra, Portugal

The prediction of rockfall trajectories below a rock cliff is essential in susceptibility, hazard and risk maps, particularity close to populated areas. The Monte da Lua hill area in Portugal, a tourist destination close to the historic city of Sintra (UNESCO World Heritage), is a typical granite boulder chaos landscape where from time to time rockfalls occur, the last such event having occurred on 29 January 2002. This area is therefore suitable to develop a rockfall study in order to provide hazard and risk maps a basis for mitigation measures. A preliminary investigation of the area leads to the identification of 188 potentially dangerous boulders. Detailed locations and geotechnical characteristics in terms of geometry, strength and context were sampled for each boulder. Digital elevations at 1 × 1 m resolution, known rockfall trajectory and building locations are provided in a GIS project for the study together with the spatial database of boulder characteristics. The modelling approach was conducted in two steps: (1) discrimination of the boulders in terms of static and dynamic mobility behaviour with multivariate analysis; (2) stochastic simulation of rockfall trajectories. The rockfall trajectory algorithm proposed is straightforward and is only dependent on elevation data, initial location of boulders and a friction angle. Due to the slope of the area, it assumes that rockfall is always of the rolling or sliding type. The friction angle was calibrated on the basis of the rockfall travel distance recorded on 29 January 2002 and generates simulated “realistic” trajectories. A smaller friction angle increases all simulated trajectories, leading to more “pessimistic” scenarios. The combined analysis of trajectories and potential damage to buildings and discrimination in terms of static and dynamic behaviour provides a final table in which all 188 sampled boulders are classified in one of the five risk grades.     

Numerical Simulation in Rockfall Analysis: A Close Comparison of 2-D and 3-D DDA

Accurate estimation of rockfall trajectory and motion behaviors is essential for rockfall risk assessment and the design and performance evaluation of preventive structures. Numerical simulation using discontinuous deformation analysis (DDA) is effective and helpful in rockfall analysis. Up to now, there have been many reports on application of two-dimensional (2-D) DDA programs. In this paper, the major advantages of rockfall analysis using 2-D and extensions to three-dimensional (3-D) analysis are presented. A practical 3-D DDA code is demonstrated to be capable of simulating free falling, rolling, sliding, and bouncing with high accuracy. Because rockfall trajectories and motion behaviors can be described as combinations of these four types, this demonstration indicates that the implemented code is capable of providing reliable rockfall analysis. Finally, specific tests are conducted to compare 2-D and 3-D DDA rockfall analysis in predicting trajectory and dynamic behavior. The results indicate that 3-D DDA simulations are more appropriate for rough tree-laden inclined slopes in providing detailed spatial distribution, whereas 2-D DDA simulations have better efficiency for slopes dominated by valleys and ravines. These results can help in selecting the appropriate DDA simulation for rockfall analysis.     

地震作用下康定市郭达山危岩带运动特征

康定地区多为深切河谷地貌,山坡陡峻,基岩裸露,崩塌落石多发,地震频发。为开展地震作用下崩塌运动特征和规律研究,以康定市郭达山危岩带为研究对象,采用颗粒流离散元软件(particle flow code in 2 dimension,PFC^2D)模拟危岩带不同部位崩塌源(坡顶孤石、坡体上部碎裂岩体、坡体中部老崩塌堆积体、坡体下部块状危岩)在面波震级(surface ware magnitude,Ms) 8.0地震作用下的运动特征和破坏过程。研究结果表明: ①坡顶孤石质量越小越易启动,孤石越接近球形越易发生倾覆和滚动,沿临空面飞出后运动类型以坠落、碰撞、滚动为主,坡顶孤石的运动速率最大时达11.8 m/s; ②坡体上部碎裂岩体破坏过程可分为裂隙延伸贯通—启动坠落—碰撞解体—滚动堆积4个阶段,位于碎裂岩体上部的块石运动距离最远,达269 m; ③老崩塌堆积体块石自前向后形成碎屑流,沿坡面运动类型以滚动、碰撞为主; ④下部块状危岩运动特征为启动—滑动—挤压—解体—再滑动—再挤压—堆积; ⑤不同部位的崩塌在运动过程中块石会发生碰撞、摩擦、挤压、解体,快速消耗自身动能,导致运动距离和速率骤降。采用离散元模拟能够更全面更精细化的认识深切河谷区的崩塌,可为崩塌灾害的工程治理和山区城镇的防灾减灾提供科学依据。     

A simplified approach for rockfall ground penetration and impact stress calculations

This paper presents a simple procedure for calculating rockfall penetration and ground impact stress for design of pipeline burial depth. Rockfall may damage pipelines by penetrating through cover soils and damaging the pipeline by puncture or excessive impact stress. Design engineers may face a challenging question regarding what burial depth is required for adequate pipe protection. Greater burial depth may provide better protection; however, there are usually site and economic constraints. This paper discusses a simple procedure for calculating the maximum rock penetration and ground impact stress due to rockfall. The procedure may be used to estimate the minimum burial depth required to protect pipelines from rockfall damage. A case study involving a large boulder that fell down a slope and landed on a buried high-pressure gas pipeline is presented to demonstrate and verify the procedure. A calculation was performed before the removal of the boulder to evaluate whether any damage was made to the pipeline. The calculated results agreed well with the actual conditions observed after the boulder was removed.     

基于三维非连续变形分析的巨石崩塌运动研究

巨石作为一种大体积、大质量的滚石,其崩塌失稳及高速、高能远程运动往往导致沿途建筑物和交通线路的毁灭性灾难。以西藏自治区G318国道K4580典型滑坡为工程背景,基于三维非连续变形分析（three-dimensional discontinuous deformation analysis,简称3D-DDA）方法研究巨石崩塌失稳及运动全过程的特征与现象。分别建立该边坡未滑坡、浅层滑坡后和深层滑坡后3种坡形的巨石崩塌3D-DDA数值模型。采用滚石运动横向偏移经验模型,验证3D-DDA巨石运动模拟的准确性。在此基础上,分析巨石崩塌失稳机制及破坏后沿3种不同坡形边坡的运动轨迹和动能演进等运动特征。结果表明,3D-DDA能够有效模拟巨石崩塌失稳、运动发展、剧烈冲击碰撞直至最终静止等整个动力学过程。巨石崩塌表现为滑动→倾倒?滑动→倾倒→翻转?下落的失稳模式转换;巨石运动表现为碰撞、弹跳、飞跃、滚动、滑动等多种运动形式以及横向偏移、侧向偏转等三维空间运动特征,经过道路并与高架桥发生碰撞,引发巨石灾害。不同坡面几何特征下的巨石运动偏移量、弹跳高度、运动至坡底碰撞时间、最终稳定时间等均随着未滑坡、浅层滑坡后、深层滑坡后3种坡形变化而减小。通过3D-DDA巨石崩塌运动分析,预测巨石运动全过程、影响范围、冲击能量、停积位置等,可为巨石防灾减灾对策或措施制定提供依据。     

Overturning stability of L-shaped rigid barriers subjected to rockfall impacts

Reinforced concrete barriers are commonly used as defence measures in hilly areas to contain falling boulders and landslide debris. These barriers are conventionally designed to satisfy the conditions of force and momentum equilibrium with a factor of safety. A major limitation of this approach is that the inertial resistance of the barrier is neglected such that the design could be over-conservative. This paper presents a novel displacement-based approach for the assessment of overturning stability of rigid L-shaped barriers subjected to rockfall impacts. Analytical solutions, which are derived based on conservation of momentum and energy, are used to take into account the contributions of the self-weight and, thus, the inertial resistance of the barrier in resisting an impact. The actual amount of energy transferred from the impacting boulder to the barrier is considered by including the coefficient of restitution between the two objects. The accuracy of the analytical solutions has been confirmed by laboratory impact experiments. Numerical assessments conducted using the new solutions indicate that a reasonably sized rigid barrier, due to its own inertial resistance, may adequately withstand the impact action of a heavy boulder rolling down a hillslope without relying on any anchorage to its support. A range of geometric design of the barriers with L-shaped cross sections also has been considered and analysed. The new approach presented in this paper is easy to apply in practice and will be useful for engineers designing concrete barriers as passive rockfall mitigation measures.     

地震荷载下危岩运动特征的模拟研究

西南地区地形复杂,位于地震高发区域。为了进一步研究该区域地震作用下落石的运动特征,以我国九寨沟景区内荷叶寨后山危岩带2为例。首先研究危岩带2中单一危岩体在自重荷载下的失稳情况,分别采用RocFall和非连续变形分析方法(DDA)进行模拟,结果表明:DDA和RocFall的模拟结果相近。DDA中落石水平位移模拟结果为210.6 m,RocFall中大部分落石水平位移分布在205 m附近。鉴于RocFall的局限性,采用DDA方法继续研究地震作用下危岩体的密度、高度、体积和地震波峰值加速度等因素对落石运动特征如运动轨迹、水平位移和速度的影响规律,并监测落石在坡脚附近的最大弹跳高度。研究结果表明:(1)危岩带2所处边坡自然坡角范围为65°~75°,落石在坡面上加速,在高程约2292 m(与居民区相对高差约为132 m)处发生抛射,运动方式以滚动、碰撞、弹跳为主;(2)岩体软弱结构面的贯通时间随落石密度、体积以及地震波峰值加速度增大而减小;(3)峰值速度随边坡高度、体积及地震波峰值加速度增大而增大;(4)水平位移随落石密度、边坡高度及地震波峰值加速度增大而增大;(5)坡脚处落石最大弹跳高度随各因素变量增大而增大。继续研究内部存在节理裂隙的岩体(3条平行节理将危岩分为K1~K4部分)在地震荷载下的运动特征:(1)在地震荷载作用下,危岩沿节理面发生倾倒破坏,落石水平位移模拟结果为191.5~243.8 m;(2)落石在坡脚附近的最大弹跳高度范围为0~13.8 m。     

The Son Poc rockfall (Mallorca,Spain) on the 6th of March 2013: 3D simulation

The Son Poc rockfall took place on the 6th of March 2013 in the municipality of Bunyola, on the southern side of the Tramuntana Range (Mallorca) and after a rainy and cold period on the region. A volume of rock of 4.000 m³ was detached from the cliff crowning the peak falling down by toppling. The impact of the boulder caused its fragmentation, and numerous boulders bounced and rolled downslope with volumes from 1 to 35 m³, following two trajectories: southwest (SW) and southeast (SE). The SE trajectory, with a larger runout (376 m), reached an urban area, where some of the boulders hit the roofs and walls of nearby houses, stopping others in their gardening areas. Fortunately, no fatalities occurred despite of the presence of some people at that moment, but the event caused great concern in a region which lives from and for tourism. The Son Poc rockfall has been simulated using RocPro3D software which uses GIS technology to produce 3D rockfall trajectories lines, estimated velocity and energy of falling blocks, as well as bounce heights, impacts, and stopping points. The results are in agreement with field observations and with a very good accuracy between real and modeled outcomes.     

基于线黏弹性接触的落石冲击地面参量理论研究

落石与坡面或者防护结构接触过程中冲击参量随时间变化特征是描述落石碰撞过程的重要指标,对于揭示落石与坡面相互作用机制以及采取合理的防护措施具有重要的意义。在现有的相关设计规范中,并没有给出关于落石冲击力时程关系的计算方法,仅参照有关规范或经验方法确定一个落石最大冲击力值。为此,首先基于线黏弹性接触理论,建立落石冲击地面力学模型,根据位移-速度组合初始条件以及速度-加速度组合初始条件分别推导得到两种落石冲击特征参量理论解析解;然后基于 ANSYS/LS-DYNA 非线性动力学软件,建立落石冲击地面三维数值模型,研究球体落石冲击地面力学特点;最后将理论结果对比室内试验和已有的研究成果,得出以下结论：（1）Hertz 弹性接触理论结果中各参量变化在加载阶段和恢复阶段均呈现对称的趋势,速度、加速度初始条件下的落石冲击特征参量和动力有限元法非常接近,而位移-速度初始条件组合并不适用于研究落石冲击下的动力特征;（2）不同速度和下落高度下,落石最大冲击力值随落石下落高度和冲击速度的增大而增加,而落石冲击作用时间随下落高度和冲击速度的增加而减小;（3）计算结果得到的落石最大冲击力以及落石冲击作用时间与室内试验结果和已有成果相接近,相比室内试验和有限元结果的震荡性,此结果更能体现落石冲击力变化规律; （4）多种冲击速度下,对比不同方法得到的落石最大冲击力,可知计算结果均在各种冲击力计算结果的范围内,具有很好的可靠性。考虑到现有研究理论的不足,难以求解落石冲击力时程关系,求解结果丰富了落石碰撞理论,可以指导工程有关落石灾害的防护设计。