Identification of deformation and failure characteristics in cataclinal slopes using physical modeling

This study investigates the deformation characteristics of cataclinal slopes in central Taiwan prior to landslide failure. Field surveys and physical model tests were performed to explain the gravitational deformation characteristics of cataclinal slopes under various conditions and to derive the deformation process and failure characteristics. The results show that the distribution of erosion gullies (different length of the slope mass), the extent of erosion (different thickness of the slope mass), the foliation dip angle, and the geological material critically affect the deformation of cataclinal slope masses in the study area. The results of physical model tests indicate that increasing the foliation dip angle, the thickness and the length of sliding mass, particle size (spacing between foliations) increases the depth of slope deformation. Foliation dip angle is the most critical factor that controls the deformation of slate slopes. When the cataclinal slopes reached maximum deformation, a shear failure and translational slide occurred within a short period. The deformation zone exhibited significant cracking at the scarp and the bulging of the slope toe, which facilitated the infiltration of surface water and groundwater, accelerating the deformation to failure.     

Dynamic behavior of railway embankment slope subjected to seismic excitation

To reveal the dynamic behavior of a railway embankment slope subjected to seismic excitation, a shaking table model test was performed on a 1:8 scale embankment slope. Different types of seismic wave of differing amplitudes were applied to study the dynamic behavior of the embankment slope, and white noise excitations were interspersed among the seismic waves to observe the changes of dynamic characteristics of the embankment slope. Residual deformation behaviors of the embankment slope were also investigated. The results of the tests show that the natural frequency of the embankment slope exhibits a decreasing trend and that the damping ratio exhibits an increasing trend. The embankment slope exhibits a significant amplification effect on the input acceleration, and the acceleration response differs greatly when subjected to different seismic excitations of differing spectrum characteristics. The filler of the embankment slope affects the changes of the spectrum characteristics of the seismic wave. The filler performs a filtering effect on high-frequency seismic waves and amplifies the energy of low-frequency seismic waves, especially when the frequency is close to the natural frequency of the embankment slope. A bidirectional excitation creates a greater acceleration response than a unidirectional excitation does. The seismic residual deformation of the embankment slope occurs under the seismic subsidence mode.     

小坡度海底土层地震液化诱发滑移分析方法

地震可使海底砂质、粉质土层液化并导致上部土层的滑移。基于有效应力有限元动力分析方法和Newmark刚性滑块理论,提出了一种计算海底小坡度（≤5o）土层地震液化引起侧向滑移的简化方法。该方法将波浪荷载简化为海底恒定的上覆压力和初始孔压,忽略了海水粘性对海底土层地震反应的影响,利用改进的Seed孔压模型进行动力分析和液化判别,用Newmark滑块理论计算了土层侧向滑移。通过算例和对比分析,研究了海水深度和土层坡度对侧向滑移的影响,表明该方法的有效性,可为近海工程场地地震地质灾害评价提供参考数据。     

土质边坡动力稳定性试验研究

以云贵高原某典型土质边坡为原型,采用了4种加速度震动波输入模式,设计完成了比例为1:6的小型振动台模型试验,结合FLAC3D数值模拟对边坡动力响应特性和边坡变形破坏规律进行分析。结果表明:当输入加速度低于某个临界值时,整个坡体的加速度响应基本保持一致,各部位放大效应增加不明显,当输入加速度逐渐增加,高于临界值时,坡体的卓越频率得到充分激励,各部位加速度响应大幅增加,此时边坡最易发生变形破坏,且加速度响应沿着坡高方向有显著的放大效应;剪应变增量时程曲线反映出在边坡震动破坏过程中,滑体后缘以张拉为主,中部及下部以剪切为主,而且剪出口剪应变增量的变化尤为关键,其增幅速度直接导致该部位抗剪强度降低速度增快;边坡震动变形破坏模式为崩塌-剪切滑移破坏,变形破坏过程可分为4个阶段。      

Analysis on mechanism of landslides under ground shaking: a typical landslide in the Wenchuan earthquake

A high steep rock hill with two-side slopes near National Road 213 is used as a prototype in this paper. The full process from initial deformation to sliding of the slope during ground shaking is simulated by a new discrete element method—continuum-based discrete element method. Then, the seismic responses of a high steep rock hill with two-side slopes are researched from the base of time, frequency and joint time–frequency domain using Hilbert–Huang transform and Fourier Transform. The findings are: first, the stress concentration phenomenon occurs at the top of the sliding mass, and then some tension and shear failure points appear, which expand from the top toward the toe of the sliding mass along the structural plane. At the same time, the number of tension failure points gradually increases. Then the toe of the sliding mass fails, and shears out from its toe which results in the landslide. If the material parameters are under the same conditions, the landslide in the middle of the slope occurs before that at the foot of slope, and the starting time of landslide and the arrival time of the peak ground acceleration are synchronous or the former slightly lags behind the latter. The difference of distribution and dissipation of earthquake energy in the sliding body and sliding bed is the major influence factor to induce the landslide. When the accelerations are small, the instantaneous frequency of accelerations between sliding bed and sliding body is generally consistent, the energy transmittance coefficients of the sliding structural plane and the controlled frequency band of the energy all range in a limitation; with the increase of the seismic intensity, the instantaneous frequency and the energy transmittance coefficients gradually decrease, and then they are steady within the lower limitation. At the same time, the controlled frequency band also shifts gradually from high frequency band to the lower one. Based on the input seismic wave, the peak acceleration amplifies as the increase of elevation, regardless of the monitoring points on the steep slope, gentle slope side or inside of the slope. Generally speaking, amplification of the vertical peak acceleration is stronger than that of the horizontal peak acceleration, and amplification of the peak acceleration on a steep slope is stronger than that on a gentle slope, and that of inside of the slope is the weakest amplification.     

层状反倾-顺倾边坡倾倒变形形成条件及发育规模特征

以往研究中倾倒变形研究在反倾边坡较多而在顺倾边坡中很少,对两种不同结构的边坡形成倾倒体的异同点更是存在较多空白。在总结大量倾倒体实例基础上,对层状反倾和顺倾边坡倾倒变形形成条件及发育规模特征进行了详细研究。结果表明,顺倾边坡若发生倾倒,通常表现为坡高 100 m,边坡坡角 35°,岩层倾角 60°,岩性为薄层或薄层与中层互层的软岩、软硬相间的岩石;反倾边坡当坡角 30°及岩层倾角 30°就可能发生倾倒,其岩性为薄层～中厚层状的软岩、硬岩及具有似层状结构的坚硬岩石均可;提出层状岩质边坡“倾倒临界倾角 ”的概念,对于顺倾边坡, 60°,当 时边坡将可能产生倾倒破坏,当 时边坡通常产生顺层面的“滑移-弯曲”或“滑移-拉裂”型破坏;反倾边坡 30°,当 时边坡才可能演化成明显倾倒变形,当 时边坡不会倾倒或倾倒不明显;对于地质条件基本相同而坡体结构不同的两种边坡,反倾边坡形成的倾倒体无论是发育分布面积还是倾倒深度通常是大于顺倾边坡的,而且一旦形成倾倒体,二者的规模通常是深层的、大型或特大型的;倾倒体的分布面积和倾倒深度均呈现一种随坡高的增加而增加的趋势,坡高 250 m时其分布面积和倾倒深度表现为陡然增加,并且反倾边坡增加幅度大于顺倾边坡。     

高陡边坡组合式支挡结构大型振动台试验研究

为研究上部采用锚索框架结构、下部设置为桩板墙的组合式支挡结构的抗震性能,开展了大型振动台模型试验。通过输入不同强度大瑞波测试了组合结构的加速度响应和动土压力响应,同时将动土压力的实测值与规范法、Mononobe-Okabe法（M-O法）的计算值进行了对比研究。研究表明：（1）组合式支挡结构在不同强度地震动激励下水平向和竖直向加速度沿坡高有不同程度的放大,且输入的激振加速度幅值越大,放大效应越明显;（2）组合式支挡结构下部桩板墙墙后的动土压力强度随激振加速度幅值的增大而增大,沿着墙高呈现出上部小、下部大的分布特性;（3）动土压力强度的计算,在低烈度区使用规范法和M-O法是合理的,但在高烈度区需要对规范法和M-O法进行修正。     

地震作用下堆积体边坡的坡面变形与失稳机制

地震作用下堆积体边坡的动力响应特性十分复杂,单一抗震安全系数不足以评价其动力稳定性。通过大型振动台试验,研究了连续多级地震荷载作用下,地震波的类型、卓越频率及峰值加速度对堆积体边坡坡面永久位移的影响,并初步分析其失稳机制。试验结果表明,相同峰值加速度下振动型地震波比冲击型地震波更容易产生坡面永久位移,地震波卓越频率对坡面永久位移也有重要影响;堆积体边坡在峰值加速度apeak=0.2g时开始有大颗粒石砾滚落,对应的坡面永久位移在apeak=0.2g～0.3g之间开始产生并显著增大,另外利用考虑坡面几何形态变化的改进Newmark法对坡顶的永久位移进行了估算。通过坡面永久位移评价堆积体边坡的动力稳定性有一定合理性。     

含顺层断续节理岩质边坡地震作用下的破坏模式与动力响应研究

基于二维颗粒流软件PFC2D的人工合成岩体技术（SRM）,研究了岩桥倾角和节理间距不同组合形式的含顺层断续节理岩质边坡在地震作用下的破坏模式与动力响应规律。研究结果显示：在地震动力作用下,含单潜在滑动面的顺层断续节理岩质边坡呈现出滑移-倾倒的混合破坏特征,含多潜在滑动面的顺层断续节理岩质边坡则主要发生倾倒破坏;由顺层断续节理以及岩桥交替连接所组成的潜在滑动面是控制边坡动力稳定性的关键因素。在地震动力作用下,最靠近坡脚的岩桥段首先萌生翼裂纹,使得拉应力得到释放,随后各节理相继萌生裂纹并扩展、贯通,最终导致坡体发生阶梯状整体失稳。裂纹扩展受顺层断续节理控制,萌生裂纹中以张拉裂纹为主,且裂纹数量与输入地震波的加速度曲线具有同步性。另一方面,节理面的存在对边坡动力响应产生明显影响,沿坡表以及沿水平方向上的峰值速度、峰值位移随着岩桥倾角的增大、节理间距的减小而增大,同时节理间距和岩桥倾角对于峰值加速度（PGA）放大系数的影响范围主要集中在坡表、坡肩;沿竖直方向上,峰值位移随着岩桥倾角、节理间距的增大而减小,PGA放大系数曲线随高程变化总体呈现U型分布特征。     

Influence of strong motion characteristics on permanent displacement of slopes

Landslides triggered by moderate to major earthquakes are a recognized seismic hazard. Arias Intensity (I _a) is a key intensity measure of the ground motion, but significant duration is widely used to define strong motion duration. We calculate Newmark’s displacements using earthquake records bracketing a broad range of Arias Intensity and significant duration employing Newmark’s rigid block method and a number of yield accelerations. Total landslide displacement increases with the increase in the energy content of the ground motion (I _a) above a threshold. Such threshold may be expressed as a function of yield acceleration of the slope regardless of the ground motion characteristics. Newmark displacement decreases with increasing duration for earthquakes with similar energy content. The wide scatter in the results converges when using formal dimensional analysis. Self-similar symmetry may facilitate the assessment of the performance of slopes during earthquakes. The mathematical framework for probabilistic determination of landslides displacement may be a useful aid to estimate the likelihood of landslide hazards provided that the geotechnical properties of the slopes are known.     

Insight into the seismic response of earth dams with an emphasis on seismic coefficient estimation

This paper presents a numerical investigation of the seismic response of earth dams by employing results from 110 nonlinear two-dimensional (2D) dynamic analyses of four different cross-sections with heights ranging from 20 to 120 m. The analyses were of a parametric nature, considering the effects of seismic excitation characteristics (intensity and frequency content), foundation soil stiffness, and the existence of typical stabilising berms and/or an impounded reservoir. The results of these studies indicated that the predominant period of a dam’s vibration was strongly affected by its height and the input motion characteristics. The results also indicated that the peak acceleration at the dam’s crest was governed by its height, the input motion characteristics, and the stiffness of the foundation soil, but not by the other parameters. These same analyses yielded results on pseudo-static seismic coefficients for a total of 1084 potential sliding masses within the analysed cross-sections, demonstrating that the seismic coefficients decreased as the sliding mass became deeper and bulkier, increased if the mass was located upstream rather than downstream, and were strongly affected by the seismic excitation characteristics and stiffness of the foundation soil. Moreover, these results allowed for a thorough evaluation of existing methodologies for seismic coefficient estimation, quantifying their accuracy and depicting their limitations. This evaluation process also illustrated the fact that there is currently no methodology accounting for all significant problem parameters.     

Applicability of Newmark method at regional, sub-regional and site scales: seismically induced Bullas and La Paca rock-slide cases (Murcia, SE Spain)

In this paper, the applicability of the Newmark method at regional, sub-regional and site scales has been investigated in the Lorca Basin (Murcia). This basin is located in one of the most seismically active regions of Spain. The area is very interesting for studying earthquake-induced slope instabilities as there are well-known cases associated with specific earthquakes. For the regional and sub-regional scales, a geographic information system has been used to develop an implementation of Newmark sliding rigid block method. Soil and topographic amplification effects have been particularly considered. Subsequently, ‘Newmark displacement’ maps for deterministic seismic scenarios have been produced. Some specific studies have also been performed using limit equilibrium methods to estimate the safety factor and the critical acceleration of certain slope instabilities at a site scale. These instabilities were the rock slides related to recent seismic series at the Lorca Basin: 2002 Bullas (M _w = 5.0) and 2005 La Paca (M _w = 4.8). Finally, the safety factor, critical acceleration and Newmark displacement values estimated at different scales have been compared to determine which scale is most suitable for the Newmark method.     

穿越不同密实度饱和砂土地层的盾构隧道地震响应三维数值分析

地震液化对隧道结构有重大威胁,且位于不同抗液化能力地层交界处的盾构隧道段更易发生严重的地震破坏。采用三维数值方法研究穿越不同密实度状态饱和砂土地层的盾构隧道的地震响应规律。饱和砂土用一种描述不同密实度砂土液化行为的边界面模型进行模拟,首先通过隧道液化上浮的振动台试验结果验证该本构模型的合理性。其次,应用多自由度连接弹簧表征管片环间相互作用,采用文献中的拼装管片的逐级加载试验结果验证该方法的可行性。最后,建立穿越两种不同密实度饱和砂土地层的盾构隧道三维数值模型,研究相对密实度、输入加速度峰值和交界面倾角对砂土地层-盾构隧道系统动力响应的影响。结果表明,可液化地层中隧道结构位移模式是水平地震激励下产生的水平位移与由于液化上浮效应产生的竖向位移的耦合作用,加之隧道在不同土层中变形存在差异,从而导致隧道呈现扭转的变形形态。在靠近交界面处,隧道整体上浮量急剧变化且该处结构上浮量随着交界面倾角增大而增大,同时管片结构弯矩出现突变,接头螺栓的环间剪切和拉伸位移也显著增加。分析结果进一步印证地震作用下盾构隧道在不同性质饱和砂土地层交界面处更易破坏,在设计阶段应予以重点关注。     

地震作用下预应力锚索边坡的稳定性分析

稳定性分析对边坡工程具有重要的意义。为探究地震对预应力锚固边坡的影响,基于预应力锚索具有安全储备的特点,分析了地震作用下边坡的分阶段模型。在传统Newmark法的基础上,采用拟动力计算,推导了顺层岩质边坡的位移、锚固力和安全系数计算公式。结合工程两个算例,探究了不同锚索模型和地震累积作用下的边坡变形规律。由结果可知,当锚索锚固力取固定值时,边坡计算位移偏大,随时间呈线性增长。当考虑锚索的安全储备时,位移增长速率随时间逐渐减小,最后趋近于定值。不同锚索模型对位移影响较大,锚索储备作用不可忽略。同时地震会对预应力锚固边坡产生永久性影响。当再次发生扰动时,其表现为临界滑动阈值的提升,位移增长率与受震经历的加速度幅值呈反指数变化。以露天矿区边坡为背景,在抗震设防标准下确定最佳支护方案。     

A simplified approach for studying the influence of vertical accelerations on seismic response of slopes

ABSTRACT

A simplified approach is presented for estimating permanent displacements in slopes as a result of both vertical and horizontal seismic accelerations. A study of 52 earthquake records showed that the time difference between maximum horizontal and vertical accelerations varied between 0 and 10.3 s. The approach is illustrated for an earth dam embankment by analysing the effects of five of the above earthquake records. The approach combines a pseudo-static slope stability analysis for estimation of the critical (or yield) horizontal-vertical acceleration combinations, and a Newmark type displacement analysis. Guidelines are presented for conservative choice of soil strength parameters of saturated clays for use in the stability analysis. While permanent displacements of up to 40 cm were predicted without considering the vertical acceleration component, no additional displacement above 3.5 cm resulted when this component was included. The predicted additional displacement was consistently less than 10%, and in 50% of the analyses, vertical acceleration led to smaller predicted displacements. The simple approach may be applied in analysis for any slope using real earthquake records. Using existing, empirical expressions for permanent displacement, based only on horizontal accelerations, the effect of the vertical accelerations may be conservatively estimated by increasing the displacement by 10%.     

重庆武隆鸡尾山滑坡视向滑动机制分析

重庆武隆鸡尾山山体为典型的斜倾厚层灰岩山体,其破坏模式不同于常见侧向崩塌-堆积层滑坡,属于真倾向滑移变形转为视向整体滑动的特殊失稳模式。在现场地质调查的基础上,从地层岩性、岩体结构、岩溶及地下水作用、软弱夹层等影响因素分析重庆武隆鸡尾山滑坡形成原因;根据滑坡破坏机制,基于关键块体控制理论,对鸡尾山滑坡进行三维稳定性极限平衡分析;利用三维离散元软件模拟鸡尾山滑坡的初始变形破坏过程,分析了鸡尾山滑坡不同影响因素条件下的视向滑动形成机制和变形破坏特征,并探讨了节理化和溶蚀岩体的参数取值。结果分析认为,在重力的长期作用下,鸡尾山山体初始沿真倾向方向滑移,沿岩溶发育的陡倾节理裂隙逐渐产生后缘及侧向裂缝,形成块状后部驱动块体,由于地下水等因素使软弱夹层软化,驱动块体下滑力增大,前缘阻滑关键块体内部应力积累,并最终沿强度较低的岩溶发育带发生剪切破坏,从而导致整体滑动;在进行滑坡稳定性极限平衡分析时,考虑实际的滑坡机制,将滑体分为驱动块体和关键块体分别进行力的解析,并以关键块体的安全系数代表滑坡的安全系数更加合理。数值模拟显示,软弱夹层强度降低、岩溶发育带剪切破坏后,滑体进入大变形阶段,表明关键块体控制和阻滑作用明显,软弱夹层强度降低是滑坡发生的关键因素。采矿形成的采空区对山体的影响主要是使上覆岩体压应力增大,但对滑体的变形无影响。     

考虑土体抗拉强度的边坡永久位移极限分析

通过野外观测与室内试验发现,边坡后缘往往存在拉应力区。拉应力区的存在会影响边坡的稳定性,而地震荷载的存在会放大这种影响。分析拉应力区对边坡稳定性的影响,当前主要采用的方式为:对强度准则中抗拉强度进行折减（即张拉截断）。文章通过极限分析上限原理和拟静力法,推导出边坡临界加速度计算方程。以边坡在不同参数组合下的位移系数为基础,输入实测地震波,采用改进的Newmark法对边坡进行位移分析。文章算例的结果表明:拉应力区的存在会大大降低边坡临界加速度,土体在完全张拉截断下的临界加速度对边坡可能会产生超过50%的折减。拉应力区的存在也可以使永久位移达到传统的摩尔库伦理论计算值的2倍之多。文中所有的结果皆以图表形式展示,非常便于理解以及读取数据。      

西安金盆水库放水塔附近滑坡特征及成因分析

放水塔附近滑坡是金盆水库右岸原1号滑坡体下游边界段的残留体,受短期强降水因素的诱发导致坡体变形、失稳滑动。滑坡体具顺层牵引滑动特征,且平面上,不同区段变形破坏程度不同,其表现形式与坡体平面旋转有一定的相似性。坡体滑动主要受4组结构面的控制,其中的软弱片理结构面（产状1501703555）与产状为2302903555的另一组结构面构成滑坡的滑动控制面。基岩内发育的大量的软弱片理结构面,大气强降水,滑坡上部相对平缓的地貌及人类工程对地表植被环境的破坏,滑坡体下部喷护层的存在等因素的综合影响,导致了坡体的失稳滑动。     

黄土基覆层边坡动力破坏特征的大型振动台试验研究

为了探讨黄土基覆层边坡动力破坏特征与加速度的响应关系,采用1∶20振动台试验,设置输入地震波幅值逐级增大,实时监测边坡裂缝试验全周期内发育规律,应用MATLAB动力破坏特征检测系统获取地表裂缝的基本信息,提出边坡表面动力破坏特征,并结合各工况监测点加速度峰值变化规律进行分析。得到如下结论:(1)输入加速度峰值0.6g时,坡顶和坡脚的裂缝宽度、裂缝面积均有跳跃性增长,表明土岩接触面部位土体已发生累进剪切破坏,滑面和坡面上的加速度峰值响应突变,振幅突变较大,表明边坡已经发生破坏。(2)沿坡面和滑面加速度放大系数均呈非线性增加,而沿坡面的放大系数在各工况下明显比沿滑面的大,说明加速度沿高程放大效应明显。(3)拉裂缝和剪切裂缝面积突增是边坡破坏的重要特征,贯通裂缝产生和加速度响应突变可以作为边坡动力破坏的依据。     

地震作用下顺层岩质边坡动力响应和破坏模式大型振动台试验研究

5•12汶川大地震触发了大量的顺层岩质滑坡,对其进行研究很有必要。根据动力模型试验的相似关系,设计制作了1个坡角大于岩层倾角的尺寸(高×长×宽)为1.6 m×1.75 m×0.8 m的顺层模型边坡,并完成了大型振动台试验。试验结果表明,在坡体表面和内部竖直方向上,加速度放大系数随着坡体高程增加而增大,并且随着高程增加,加速度放大系数增大的速度加快;在坡体内同一高程上,坡面处的加速度放大系数大于一定水平深度坡体内部的加速度放大系数,表现出趋表效应;地震波输入频率对坡体动力响应有明显影响,随着频率的增加,越接近坡体的自振频率,加速度放大效应越显著;加速度放大系数随着输入波振幅的增加,总体上表现为递减趋势;通过和均质边坡振动台试验加速度监测数据对比,发现坡体结构对坡体加速度放大系数也有一定的影响,结构面对地震波的反射和折射作用加大了坡体加速度的放大效应。,对试验过程中坡体破坏特征的描述和分析发现,边坡的破坏模式为地震诱发－坡肩拉裂张开－坡面中部出现裂缝－裂缝贯通－发生高位滑坡－转化为碎屑流－堆积坡脚。研究成果对地震灾区滑坡形成机制的认识和减灾防灾有一定的价值。