Impact force to a rigid obstruction from a granular mass sliding down a smooth incline

The impact force to a rigid obstruction from a granular mass sliding down a smooth incline provides insights into the solid-like and fluid-like behaviors of granular avalanches and useful information for risk assessment and engineering design against landslides. In this study, a series of 2-D flume tests were performed to systematically investigate the effects of inclination angle, sliding distance, and initial relative density on the flow front velocity and impact force on a rigid obstruction. The experimental results show that for inclination angles smaller than the critical state friction angle of sand, an increase in the sliding distance and/or initial relative density results in smaller impact forces; for higher inclination angles, the trend is reversed. Based on the experimental results, an analytical equation is proposed to estimate the flow front velocity and an empirical approach is presented to estimate the maximum impact force based on elastic solid and hydrodynamic methods. The proposed equations are found to provide more accurate predictions for the maximum impact force than similar equations in the literature.     

Experimental study on impulse waves generated by a viscoplastic material at laboratory scale

This experimental study provides insight into impulse waves generated by a viscoplastic material. The viscoplastic material chosen is a stable polymeric gel called Carbopol Ultrez 10, which is approximately modeled as Herschel–Bulkley model. As observed from high-speed cameras, the viscoplastic material such as Carbopol moves as a long and thin train of material along the slope, and only a fraction of the sliding mass is engaged in generating the leading wave. Therefore, our primary objective is to study how much of the initial slide mass is able to contribute to the leading wave formation. For the sake of distinguishing the actual slide mass acting on the leading wave formation with the initial mass, we define the submarine slide mass when the leading wave reaches its maximum wave height as “effective mass”. In this work, we held the still water depth and slope angle constant, and varied the initial slide mass and slope length. Then, we measured the slide velocity, slide thickness, and slide mass at impact, as well as the wave amplitude and wave height. The results indicate that the effective mass is dependent on both the initial slide mass and the slope length. The ratio of the effective mass to the initial slide mass is less than 20% in our experimental range, and the ratio increases with larger initial mass. In addition, we also examined our experimental data with previous empirical equations developed from granular slides. By considering the effective mass instead of the initial slide mass, the prediction of impulse waves generated by viscoplastic material is significantly improved.     

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

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

Study of the interaction between dry granular flows and rigid barriers with an SPH model

This study uses an incompressible smoothed‐particle hydrodynamics model to investigate the interaction between dry granular material flows and rigid barriers. The primary aim is to summarise some practical guidelines for the design of debris‐resisting barriers. The granular materials are modelled as a rigid‐perfectly plastic material where the plastic flow corresponds to the critical state. The coupled continuity equation and momentum equation are solved by a semi‐implicit algorithm. Compared with flows in controlled flume experiments, the model adequately reproduces both the kinetic of the flows and the impact force under various conditions. Then the numerical simulations are used to study the detailed interaction process. It is illustrated quantitatively that the interaction force consists of two parts, ie, the earth pressure force caused by the weight of the soil and a dynamic force caused by the internal deformation (flowing mass on top of a dead zone). For the estimation of impact load, this study suggests that an increased earth pressure coefficient depending on the Froude number should be incorporated into the hydrostatic model.     

二维条件下颗粒柱体崩塌入水堆积过程

在我国西南地区,以塔柱状危岩体崩塌为扰动源的涌浪次生灾害时有发生。涌浪特征与其扰动源的初始形状和破坏模式关系密切,塔柱状危岩体由下至上、空中崩解的破坏模式与颗粒柱体崩塌相似。区别于以往多试验采用的刚体材料和利用运动距离获得初始动力状态的涌浪试验,本研究设计并进行了不同高度和宽度的颗粒柱体在不同水深下自然崩塌的物理模型试验。运用PIV技术分析了颗粒体及水体的运动特征,按颗粒体与水体的相互作用关系将整个过程大致分为了3个阶段,并研究了颗粒柱体初始形态对试验结果的影响。结果表明,颗粒柱体崩塌后的残余体积受颗粒体初始形态影响较大,水深对其影响较小。涌浪产生早期,首浪位置与颗粒运动位置有着良好的对应关系;而后,两者的位置逐渐拉开。颗粒柱体高度越高,首浪最大波高出现时的位置与崩塌颗粒群前端位置的差值越大。这可能与流固相互作用后期水体流速大于颗粒流动速度有关。试验产生的涌浪大多具有完整波形,这与柱状颗粒体的散体性及其覆盖堆积特征对水体的作用呈不均匀性有关。涌浪传播过程中,水体在z方向上抬升幅度具有一定的渐变性。试验分析结果揭示了柱状危岩体引发涌浪的机制特点,为涌浪灾害防治提供了基础资料。     

Shearing of a narrow granular layer with polar quantities

The paper deals with numerical investigations of the behaviour of granular bodies during shearing. Shearing of a narrow layer of sand between two very rough boundaries under constant vertical pressure is numerically modelled with a finite element method using a hypoplastic constitutive relation within a polar (Cosserat) continuum. The constitutive relation was obtained through an extension of a non‐polar one by polar quantities, viz. rotations, curvatures, couple stresses using the mean grain diameter as a characteristic length. This relation can reproduce the essential features of granular bodies during shear localization. The material constants can be easily determined from element test results and can be estimated from granulometric properties. The attention is laid on the influence of the initial void ratio, pressure level, mean grain diameter and grain roughness on the thickness of shear zones. The results of shearing are also compared to solutions without the polar extensions. The FE‐calculations demonstrate that polar effects manifested by the appearance of grain rotations and couple stresses are significant in the shear zone, and its thickness is sensitive to the initial void ratio, mean grain diameter and layer height. The effect of the pressure level is rather low within the considered range. Copyright © 2001 John Wiley & Sons, Ltd.     

刚性块体压入颗粒体系时的受力及力链演变

颗粒体系由大量离散颗粒组成,普遍存在于自然界中,比如砂土地基、泥石流及滑坡体等。外荷载通常沿着准直形的路径在体系内传递,形成区别于单个颗粒和整个体系的细观结构--力链,亦即颗粒物质具有典型的多尺度特征。采用颗粒离散元方法模拟刚性块体压入颗粒体系的过程,计算得到了刚性块体底部所受阻力与压入深度的幂率关系;给出力链的判断准则和搜索强力链程序,构建颗粒体系强力链网络,并分析了压入试验过程中力链的演变及其长度的分布规律。     

侧限条件下充填散体压缩承载和传力特性

充填散体的压缩承载和传力特性对防控充填开采矿山的地表变形具有重要意义。通过河卵石、碎石、尾矿、河砂、尾砂5种散体材料和7组不同粒径河卵石的侧限压缩试验,以及河卵石数值模拟压缩试验,分析了其压实过程和粒径对其变形、承载的影响,揭示了散体压缩过程中的传力特性。结果表明：散体压缩过程可分为空隙压密、颗粒破碎压密和压实固结3个阶段,且散体类型和粒径对其非线性压实过程基本没有影响,但对其内部结构重组过程有较大影响;压缩过程中,存在一个临界密实度使不同粒径河卵石的承载能力次序发生反转,且河卵石散体平均变形模量与颗粒特征因子呈负指数关系;基于河卵石数值模拟压缩试验得到的散体内部力链强度分布和“压力拱”,很好地证实了室内压缩试验中散体垂直应力沿加载方向呈递减的趋势,以及侧向应力在整个散体高度上先增后减这一试验结果。     

Elastic analysis of laterally loaded pile in multi-layered soil

An analysis is developed to determine the response of laterally loaded piles in layered elastic media. The differential equations governing pile deflections in different layers due to a concentrated static force and/or moment acting at the pile head are obtained using the principle of minimum potential energy and calculus of variations. The differential equations are solved analytically using the method of initial parameters. Pile deflection, slope of the deformed axis of the pile, bending moment and shear force can be reliably obtained by this method for the entire pile length. The input parameters needed for the analysis are the pile geometry and the elastic constants of the soil and pile. It is observed that soil layering has a definite impact on pile response and must be taken into account for proper analysis and design. The analysis forms the basis for future formulations that can consider stress–strain nonlinearity.     

基于三维有限元法的巷道影响特征

在矿井瞬变电磁法探测中,巷道空间对探测结果存在影响,不利于探测资料的解释。采用三维有限元法,通过建立地质模型模拟了有无巷道时瞬变电磁场的分布情况,绘制了10 us和0.1 ms时刻有无巷道以及巷道不同长度和高度下XOZ平面磁场强度等值线分布图。结果表明:磁场在高阻区扩散速度比低阻区慢;巷道长度的变化对巷道长度方向磁场分布影响较大;巷道高度的变化对巷道上下方磁场分布影响较大。巷道空间在早期对磁场分布影响较大,晚期则影响不大,这为探测资料异常的解释提供了理论依据。      

The influence of a fixed axial top load on the dynamic response of a single pile

The influence of the gravity and inertia of a concentrated mass attached to the top of a single pile on the pile’s dynamic response was investigated in this study. Based on Muki’s fictitious pile method, the static initial axial force of the pile caused by the gravity of the concentrated mass was calculated first. Using the obtained initial axial force of the pile and taking into account the influence of the inertia of the concentrated mass, the second kind of Fredholm integral equations describing the dynamic response of the preloaded single pile were established via the fictitious pile method. Using the obtained integral equations, the dynamic response of the preloaded single pile to an incident Rayleigh wave was calculated. The influence of the pile–soil Young’s modulus ratio, the frequency of the incident Rayleigh wave and the pile length on the dynamic response of the preloaded pile was examined. The numerical results of this study show that the concentrated mass attached to the top of the pile will affect both the dynamic axial and shear forces of the pile. Generally, the concentrated mass attached to the pile top has a greater influence on the shear force of the pile than on the axial force.     

Experimental Study on the Mobility of Channelized Granular Mass Flow

Granular mass flows (e.g., debris flows/avalanches) in landslide-prone areas are of great concern because they often cause catastrophic disasters as a result of their long run-out distances and large impact forces. To investigate the factors influencing granular mass flow mobility, experimental tests were conducted in a flume model. Granular materials consisting of homogeneous sand and non-homogeneous sandy soil were used for studying particle size effects. Run-out tests with variable flow masses, water contents, and sloping channel confinement parameters were conducted as well. The results indicated that granular mass flow mobility was significantly influenced by the initial water content; a critical water content corresponding to the smallest flow mobility exists for different granular materials. An increase in the total flow mass generally induced a reduction in the travel angle (an increase in flow mobility). Consistent with field observations, the travel angles for different granular materials decreased roughly in proportion to the logarithm of mass. The flume model tests illustrate that the measured travel angles increase as the proportion of fine particles increases. Interestingly, natural terrain possesses critical confinement characteristics for different granular mass flows.     

泥石流冲击弧形拦挡坝动力响应研究

Dam foundation is subjected to a larger impact force when debris flow runs up, causing stress concentration and local impact failure. To address this problem, in this study the vertical structures are optimized into arc-shaped dams. Based on the principle of momentum and energy conservation, the theoretical calculations of the impact process of debris flow and arc-shaped dam are carried out, and the formulas of impact force and maximum run-up height of debris flow are deduced. The theoretical formulas are verified through a series of physical model tests of debris flow impact arc-shaped dam. The results show that the results of the physical model are highly consistent with those of the theoretical calculations, indicating that the proposed theoretical formulas are applicable in the calculation of the impact of debris flow on arc-shaped dam. The debris velocity, impact force and the maximum run-up height are proportional to the flume slope of debris flow. The impact force and the maximum run-up height are mainly controlled by Froude number(Fr), flume slope(?), and arc-shaped radius(R). Both the impact force and the maximum run-up height have quadratic relationships with the Froude number, and are inversely proportional to the cosine of the flume slope. Compared with the rigid vertical structures, the arc-shaped dams have no signicicant influence on the maximum run-up height, but it can reduce the normal impact force on the dam considerably, and the structure strength can also be enhanced by the strengthening of local structure. This study provides a theoretical and technical support for the dam structure design.     

竖向直排冻结条件水平冻胀力试验研究

竖向直排冻结方式常用于软弱地层中斜井掘进和深基坑止水帷幕的冻结法施工中,目前国内外对该冻结方式下水平冻胀力的分布研究鲜见报导。以某斜井冻结工程为背景,推导出温度场、应力场和水分场等相似模拟准则,设计了竖向直排冻结模型试验并细化具体试验方案,利用大型三维模拟冻结试验系统,对竖向直排冻结条件下不同深度土体的水平冻胀力分布特性进行了试验研究。结果表明,竖向直排冻结过程中埋深对水平冻胀力的影响明显,在其他条件相近的情况下同一水平面中水平冻胀力的大小与温度密切相关;冻胀力的变化率主要受冻结锋面位置的影响,距离冻结锋面较近时冻胀力的变化率逐步达到峰值;已冻土体的冻胀力随冻结锋面向外发展而趋于平稳。     

嵌岩桩抗滑特性的物理模型试验研究

进行嵌岩桩与滑坡体相互作用的框架式滑坡物理模型试验,根据微型土压力盒和电阻应变计的监测数据,分别研究滑坡推力作用下模型桩的受力特征、桩身弯矩分布规律及模型变形破坏模式。试验结果表明,嵌岩桩加固后的滑坡,桩后推力随深度的增加呈抛物线型分布形式,合力作用点约在滑动面以上模型桩自由段的1/2处;嵌岩模型桩有明显的抗滑特性,承担了大部分桩后推力,传递至桩前土压力值较小且稳定;模型桩的桩身弯矩分布形式不同于普通抗滑桩弯矩分布形式,自由段埋深0~15 cm范围内为主要弯矩承受区域,最大弯矩截面位于滑面上模型桩自由段1/3处,滑动面处桩身弯矩绝对值较小;滑坡模型在沿滑面推力加载作用下发生桩后滑体越桩滑动破坏。该试验成果为嵌岩桩的抗滑特性研究提供了科学依据,可为该类型抗滑桩设计提供一定的指导性建议。     

管缝式锚杆在拉拔荷载作用下受力分析模型

建立了管缝式锚杆在拉拔荷载作用下的受力分析模型,该模型能够反映锚杆破坏的渐进过程。分析表明,在拉拔荷载作用下管缝式锚杆在荷载作用点最先发生相对滑移,且随着荷载的增加滑移段的长度逐渐增加,直到整根锚杆与围岩发生相对滑移。根据能量守恒原理得到了管缝式锚杆在安装入钻孔后对围岩径向作用力的计算公式。对管缝式锚杆进行了非局部摩擦分析,得到了基于非局部摩擦模型的管缝式锚杆界面剪应力分布。研究了管缝式锚杆直径对锚杆界面剪应力分布的影响,同时得到了锚杆最大抗拉拔力与锚杆长度之间的关系。将理论模型预测值与试验结果对比,验证了理论模型的有效性。     

微型抗滑桩单桩设计计算模型及算法研究

考虑微型桩与周围岩土体间摩擦力的作用,提出了一种新的微型抗滑桩单桩设计计算模型,并给出了具体算法。由微型桩加固滑坡体的变形特点,分析微型桩与岩土体之间的相互作用机制,将微型桩与周围岩土体的摩擦作用引入其受力分析中;根据微型桩上各部分受力特点的不同,将微型桩分成上部摩擦受拉段、中部滑坡推力作用段和下部锚固段3段进行分析,推导了微型桩总的变形控制方程及各分段的变形控制方程;采用初参数法对控制方程进行求解,得到了微型抗滑桩上的内力分布及变形规律。计算结果表明,在滑坡推力作用下,微型桩的变形主要发生在滑面附近及以上桩段,滑面附近桩段将产生较大的内力和弯曲变形,受拉段弯曲变形较小,近似水平移动;桩与岩土体间的摩擦力是微型桩与周围岩土体相互作用的重要组成部分,摩擦力的作用能显著减小微型桩的弯曲变形,有效控制滑坡体的位移。     

Experimental Study of Dry Granular Flow and Impact Behavior Against a Rigid Retaining Wall

Shallow slope failure in mountainous regions is a common and emergent hazard in terms of its damage to important traffic routes and local communities. The impact of dry granular flows consisting of rock fragments and other particles resulting from shallow slope failures on retaining structures has yet to be systematically researched and is not covered by current design codes. As a preliminary study of the impact caused by dry granular flows, a series of dry granular impact experiments were carried out for one model of a retaining wall. It was indirectly verified that the total normal force exerted on a retaining wall consists of a drag force (F _d), a gravitational and frictional force (F _gf), and a passive earth force (F _p), and that the calculation of F _d can be based on the empirical formula defined in NF EN Eurocode 1990 (Eurocode structuraux. Base de calcul des structures, AFNOR La plaine Saint Denis, 2003). It was also indirectly verified that, for flow with Froude number from 6 to 11, the drag coefficient (C _d) can be estimated using the previously proposed empirical parameters.     

层流区Shields曲线分布规律试验

为探明层流条件下Shields曲线的分布规律,运用电荷耦合组件(Charge Coupled Device,CCD)成像技术结合激光扫描及计算机图像处理技术,在水和甘油混合液中对无黏性均匀玻璃颗粒的起动拖曳力进行测量,并根据层流绕流理论建立了无黏性均匀颗粒的滚动起动模型,提出层流条件下无黏性均匀颗粒的量纲一起动拖曳力公式.试验观测发现在层流范围内希尔兹数呈规则的带状分布,流体作用引起颗粒床表面粗化并导致床面颗粒突起减小,使颗粒起动拖曳力增大1倍以上.结果表明颗粒床表面结构性状直接影响着颗粒起动拖曳力的大小,在层流区Shields曲线具有带状分布特性.     

颗粒形态对离散介质剪切强度的影响

自然条件下,颗粒介质大多以非规则单元形态存在。非规则几何形态对颗粒介质的宏观力学性能有很大影响。针对颗粒单元的不同几何形态,采用团颗粒单元对离散介质的直剪试验过程进行了离散元数值计算,详细地讨论了颗粒形态对离散介质剪切强度的影响。该非规则颗粒由不同形态、不同数目、镶嵌尺寸、组合方位和颗粒大小的球形颗粒进行随机构造,其在局部与整体坐标之间的转动、力矩和方位关系通过4元素方法进行确定,基本球体颗粒之间的作用力采用具有Mohr-Coulomb摩擦定侓的Hertz-Mindlin非线性接触模型,并考虑了非线性法向粘滞力的影响。通过构造7种具有相同的质量概率分布的不同形态的团颗粒,在不同法向应力下,对团颗粒的直剪试验进行了离散元模拟,分析了不同形态颗粒的剪切强度。通过对不同形态颗粒介质剪切强度的数值分析,进一步揭示了非规则颗粒间的咬合互锁效应,为分析非规则颗粒的宏观动力特性提供了依据。