走滑断层相关源-汇体系演化的砂箱物理模拟实验

以走滑断层为边界的隆升剥蚀区和山前沉积区组合是盆山体系的重要类型之一。由于走滑断层两盘不断发生相对滑动,因此走滑断层两侧的源区和沉积区之间的空间关系随走滑运动不断发生变化。走滑断层两侧的源-汇系统之间的耦合机制是一个有待研究的问题。针对这一问题,本文设计了一个简化的构造地貌砂箱实验,用来揭示源区特定母岩的沉积物分布特征以及两者之间的时空关系。实验结果表明,特定母岩在沉积区的沉积信号分布范围受包括山前洪积扇形态、岩体在流域内的分布等多种因素影响。由于洪积扇向沉积区呈扇状展开的形态特征,沉积信号在离走滑断层更远的位置分布范围更宽。当一个岩体分布在不同流域内时,其沉积物分布在两个洪积扇内,从而使其沉积信号分布范围更宽。相反,规模较小的洪积扇会被相邻较大的洪积扇限制形态从而使得沉积范围变窄。由于上述原因,在利用沉积信号和岩体之间的错开距离计算获得的走滑速率,可能会低估或者高估真实的走滑速率。本文的实验为研究以走滑断层为边界的盆山耦合系统及其源-汇体系提供了重要的模型参考。     

Constitutive modeling of cyclically loaded clays with entrapped gas bubbles under undrained and unexhausted conditions

Acta Geotechnica - Fine-grained gassy marine sediments are widely distributed around the world, and the cyclic behavior of this type of sediments has a great influence on the foundations of...     

Sliding element simulating the response of slip surfaces and its application for the prediction of earthquake-induced landslide movement using one-dimensional dynamic analyses

Landslides - Earthquake-induced landslides involve excessive movement of slopes, usually along slip surfaces. This seismic movement of slopes may depend crucially on (a) the soil response along the...     

Analogue and numerical modeling of normal fault patterns produced due to slip along a detachment zone

In this paper, we investigate normal fault patterns produced by the sliding motion along a gently dipping normal fault by using analogue model tests and numerical modeling. The motivation for this study was offered by microseismic test data that indicate the existence of an active low-angle shear zone at a depth of 9–11 km in the extensional region of high seismic activity of the Gulf of Corinth (Greece). Both modeling techniques seem to support the hypothesis that the system of high-angle normal faults that are responsible for the final asymmetrical graben formation initiate at the tip of the active basal detachment nearest to the free surface. The normal faults propagate upwards with progressive sliding of the inclined basal plane, resulting in a first phase of symmetrical graben configuration that is delimited by a main synthetic fault and an antithetic fault forming a Rankine zone. Subsequent sliding on the inclined base induces a family of secondary antithetic normal faults, which are responsible for the asymmetry of the failure pattern and the diffusive character of deformation in that area. Shear deformation is more intense and localized along the synthetic normal fault than along the antithetic faults. Elaboration on the analogue test results has led to the phenomenological relations among four main parameters that describe the geometry of grabens, namely, (i) the width and (ii) the maximum subsidence of the graben, (iii) the dip angles of the conjugate normal faults, and (iv) the amount of sliding along the low-angle normal fault. However, analogue models do not produce the system of synthetic faults that is observed in the Gulf of Corinth. The effects of both friction angle variation along the detachment base and of the constitutive behavior of the model material on the configuration of the final structural pattern were also studied with a series of numerical continuum models. It was found that (a) the fault pattern of the Gulf of Corinth may be reproduced with either a strain-softening material with low elastic modulus or a constant strength material, and (b) two consecutive grabens, such as those of Gulfs of Corinth and Evia, may also be reproduced by an appropriate combination of variation of dip and frictional properties along the hypothesized detachment zone.     

爆破方法识别滑带试验研究

以往对滑坡体中软弱层与滑带的识别主要依赖于钻孔并辅之以工程师的经验。采用在钻孔中进行爆破的方法,通过对地表接收信号的波形分析即可十分方便地识别出软弱层与滑带。简要阐述了该方法的理论依据,通过对现场爆破实验测量信号的分析可以发现,该方法可以减少钻孔数量以及提高滑带识别的精度与效率。     

Constitutive modeling of unsaturated aggregated soils

A coupled water retention–mechanical constitutive model for unsaturated aggregated soils is presented here. Based on the multi‐scale experimental results, the model incorporates the inter‐particle bonding, fabric and partial saturation effects in a single framework. It is formulated within the framework of hardening elasto‐plasticity and is based on the critical state concept. Prior to model validation, we evaluate the model parameters and propose determination procedures for the main new parameters. Finally, the model is examined for its capability in simulating the experimental results of aggregated and bonded soils. Results of these simulations show that the model addresses the most features arising from the combined effects of soil structure and partial saturation. Copyright © 2010 John Wiley & Sons, Ltd.     

基于整体稳定性分析法的边坡临界滑动面搜索

临界滑动面搜索是边坡稳定分析中很重要的内容,针对现有方法的不足,提出了一种新的搜索方法。基于边坡稳定性的整体分析法,建立了确定临界滑面的非线性优化模型。该模型将滑面离散点的纵坐标和安全系数都视为独立变量,目标函数取为安全系数本身,约束条件是平衡方程组和滑面凸性。由于目标函数是线性函数,约束条件至多是二次多项式,非线性程度较低,可采用经典优化算法和常见的非线性优化工具求解,比如Matlab。通过实例分析并与传统的滑面搜索方法进行对比表明,所建议的方法在数值稳定性及收敛性方面均具有优势。     

A numerical procedure for predicting rainfall‐induced movements of active landslides along pre‐existing slip surfaces

A numerical model to predict landslide movements along pre‐existing slip surfaces from rainfall data is presented. The model comprises: a transient seepage finite‐element analysis to compute the variations of pore water pressures due to rainfall; a limit equilibrium stability analysis to compute the factors of safety along the slip surface associated with transient pore pressure conditions; an empirical relationship between the factor of safety and the rate of displacement of the slide along the slip surface; an optimization algorithm for the calibration of analyses and relationships based on available monitoring data. The model is validated with reference to a well‐monitored active slide in central Italy, characterized by very slow movements occurring within a narrow band of weathered bedrock overlaid by a clayey silt colluvial cover. The model is conveniently divided and presented in two parts: a groundwater model and a kinematic model. In the first part, monthly recorded rainfall data are used as time‐dependent flow boundary conditions of the transient seepage analysis, while piezometric levels are used to calibrate the analysis by minimizing the errors between monitoring data and computed pore pressures. In the second part, measured inclinometric movements are used to calibrate the empirical relationship between the rate of displacement along the slip surface and the factor of safety, whose variation with time is computed by a time‐dependent stability analysis. Copyright © 2007 John Wiley & Sons, Ltd.     

The effects of method of generating circular slip surfaces on determining the critical slip surface by particle swarm optimization

The main objectives of slope stability analysis are evaluating factor of safety for a given slip surface and determining the critical slip surface for a given slope. Factor of safety is usually calculated by limit equilibrium method. The main steps to determine the critical slip surface are generating trial slip surfaces as probable solutions and searching among them to determine the one with the lowest factor of safety. Although the process of searching the critical slip surface received much attention between researchers, the significance of method of generating slip surfaces is seldom addressed in the literature. The authors believe that this ignorance can affect the accuracy of the results of slope stability analysis even in the simplest problems with circular slip surfaces. Consequently, this paper focused on the method of generating circular trial slip surfaces as the simplest mechanism of sliding and considered its effect on determining the critical slip surface. A new method of generating circular slip surface was presented, which is more efficient and less restricted than the conventional method. A computer program was also developed to determine the critical slip surface of slopes by using particle swarm optimization. The performances of the proposed method and developed computer program were verified during comparative studies and sensitivity analysis. Based on the results, the effect of method of generating circular slip surfaces on determining the critical slip surface was confirmed successfully. In all considered problems, the proposed method of generating circular slip surfaces led to the lower values of factor of safety compare with the conventional method.     

Constitutive modeling of three-dimensional non-coaxial characteristics of clay

Acta Geotechnica - A simple critical state-based constitutive model is proposed to represent the non-coaxial plastic deformations of saturated clay subjected to monotonic shearing with fixed...     

不同应力路径下水土耦合超固结黏土分叉分析

基于超固结黏土三维弹塑性本构模型,导出水-土耦合部分排水条件下饱和黏土应变局部化的萌生条件,求出在不同应力路径下,单相和水-土耦合条件下超固结土分叉的三维解析解和数值解。理论分析表明,应力罗德角 为-30°、15°和30°时,无论是单相条件,还是水土耦合两相介质条件,该模型均无分叉现象产生;应力罗德角为0°时,单相条件下该模型有分叉现象产生,水土耦合不排水条件下该模型无分叉现象出现,水-土耦合部分排水条件下分叉与否则与参数选取有关;应力罗德角为-15°时,无论是单相条件,还是水-土耦合条件,该模型均有分叉现象产生。利用嵌入了上述本构模型的有限元软件ABAQUS,对单相和水-土耦合条件下的多单元立方体应变局部化分叉现象进行数值分析。分析表明,在应力罗德角为 30°时,单相条件和水土耦合条件下都无分叉现象产生;应力罗德角为15°时,单相条件下有分叉现象产生,水-土耦合不排水条件下无分叉现象出现,水-土耦合部分排水条件下何时分叉则与渗透系数大小有关;应力罗德角为-15°和0°时,单相条件和水土耦合条件下该模型均有分叉现象产生。同一应力路径下数值解相对于理论解更易出现分叉。     

考虑土拱效应的倾斜滑移面间竖向应力研究

假定两滑移面相互平行,且与水平面呈一定角度?,对滑移面间土体沿竖向取水平薄层作为微分单元体,通过作用在单元体上的水平力和竖向力的平衡条件,基于土体主应力轴旋转理论,得到考虑土拱效应及倾斜角度的倾斜滑移面间竖向应力的理论公式,及对应不同倾斜角度及滑移面-土摩擦角的土侧压力系数。将得到的理论公式与Handy等公式进行比较,验证了该公式的合理性。研究结果表明：对于无黏性土,竖向应力沿深度先近似线性增大,后增加缓慢并逐渐趋于定值。竖向应力随着倾斜角度和滑移面间距的增大而增大,随滑移面-土摩擦角的增大而减小;土侧压力系数随着滑移面-土摩擦角的增大而增大,而随内摩擦角的增大而减小。     

Cauchy problem of three‐dimensional critical slip surfaces of slopes

Using the continuum‐based finite element approach in conjunction with the strength reduction technique, one can easily find out the factor of safety of a slope with rather high precision, but will encounter some obstacles in the accurate location of the corresponding three‐dimensional critical slip surface (CSS). On the basis of the Mohr–Coulomb failure criterion and the stress field in the limit equilibrium state of the slope, it is deduced that the three‐dimensional CSS is the solution of the Cauchy problem of a quasi‐linear partial differential equation (PDE) of first order. By means of the method of characteristics for the problem, the three‐dimensional CSS that may take any shape can be determined without the specification of its geometry. Copyright © 2010 John Wiley & Sons, Ltd.     

Constitutive modeling of dense gravelly soils subjected to cyclic loading

A constitutive model for dense gravelly soils was developed to reproduce their responses under cyclic loadings. Its application aims at nonlinear dynamic analyses of earth structures involving gravelly soils, such as rockfill dams and railroad ballasts. The framework of generalized plasticity was modified to incorporate the concept of stress distance for better simulation of unloading and reloading responses. It was then combined with the theory of critical state soil mechanics to develop the constitutive model. The model has the following important features: unified simulation of particle breakage through translating critical state line, smooth transition from unloading to reloading in the stress space, and proper modeling of cyclic hysteresis, cyclic densification, and cyclic hardening of dense gravelly soils. Most of the model parameters can be obtained through simple calculation using conventional triaxial test results, and their calibration process was discussed. The model was used to simulate the cyclic responses of three gravelly soils with satisfactory accuracy. Copyright © 2014 John Wiley & Sons, Ltd.     

Geodimeter measurements of slip and strain accumulation along the san andreas fault

The U.S. Geological Survey conducts repeated geodimeter surveys of trilateration networks in central California in order to study the processes of slip and strain accumulation along the San Andreas fault. The precision of distance measurement is described by a standard deviation where a = 3mm, b = 2 · 10⁻⁷, and L is the line length. Within the precision of measurement, no anomalous strain episodes preceding earthquakes or even strain discontinuities at the time of earthquakes were detected from repeated measurements of lines near the epicenters of small (magnitude 4.5–5.1) earthquakes. Annual measurements of small (5-km aperture) strain polygons near the San Andreas fault have not proved strain accumulation in a 3-year period. Repeated measurements of longer lines over periods of 8 to 14 years indicate changes that cannot be attributed to fault slip and must represent strain accumulation at the level of a few parts in 10⁷ per year.     

Constitutive modeling of the destructuration and anisotropy of natural soft clay

The mechanical behavior of natural clays is affected by their inherent anisotropy and metastable soil structure. A simple hierarchical model that considers initial anisotropy and destructuration was formulated within the framework of critical state soil mechanics. In the proposed model, stress sensitivity and a destructuration index were introduced to account for the degree of bonding and the rate of destructuration, respectively. An inclined yield surface was used to incorporate the effect of the initial anisotropy. The proposed model can be degenerated to the Modified Cam Clay model by setting the initial stress sensitivity equal to unity and using a horizontal yield surface. Reasonable agreement between the model simulations and the experimental results on a variety of stress paths demonstrated that the proposed model can capture well the deformation behavior of natural clay and reconstituted soil. The model was implemented into the finite element program for the numerical analysis of an embankment on soft clay improved with prefabricated vertical drains. The numerical predictions were compared with the field-measured data in terms of embankment settlement. Additionally, the numerical simulations were analyzed in terms of horizontal displacements, excess pore water pressure, mean effective stress and volumetric strain. All of the simulations and comparisons indicate the importance of considering the effects of plastic anisotropy, interparticle bonding and destructuration caused by loading beyond yield stress and field disturbance in analyzing the behavior of an embankment on natural soft clay.     

Identification of representative slip surfaces for reliability analysis of soil slopes based on shear strength reduction

Although a slope may have numerous potential slip surfaces, its failure probability is often governed by several representative slip surfaces (RSSs). Previous efforts mainly focus on the identification of circular RSSs based on limit equilibrium methods. In this paper, a method is suggested to identify RSSs of arbitrary shape based on the shear strength reduction method. Monte Carlo simulation is used to generate a large number potential slip surfaces. The RSSs are identified through analyzing the failure domains represented by these samples. A kriging-based response surface model is employed to enhance the computational efficiency. These examples shows that the RSSs may not always be circular, and that the suggested method can effectively locate the RSSs without making prior assumptions about the shape of the slip surfaces. For the examples investigated, the system failure probabilities computed based on the shear strength reduction method are comparable to, but not the same as those computed based on the limit equilibrium methods. The suggested method significantly extends our capability for identifying non-circular RSSs and hence probabilistic slope stability analysis involving non-circular slip surfaces.     

基于颗粒物质力学与连续介质热力学的岩土本构模型初探

首先简要介绍颗粒物质力学与模拟岩土材料本构特性的热力学方法,其次对力链及其对应的强弱网络的形成、力学特性与能量耗散特点与机制进行深入地分析,随后在Collins提出的土体热力学模型的基础上,考虑强弱网络结构的应力应变特征,引入合理的假设,探讨建立符合热力学原理的宏细观结合的岩土本构模型的思路与步骤     

A practical procedure for searching critical slip surfaces of slopes based on the strength reduction technique

Once an earth slope is led to the limit equilibrium state by means of the finite element–strength reduction technique, a plastic zone will go through the slope from the toe to the top. Based on the observation that the critical slip surface will be within the plastic zone and be composed of the points at which the equivalent plastic strain arrives at the maximum in the vertical direction. First, some vertical lines are arranged along the horizontal direction. Then, along each of the vertical lines the point with the maximum equivalent plastic strain is found out. All such points take on a shape of wave and form the functional data. Finally, the functional data are smoothed by the least squares and the position of the critical slip surface is determined. The proposed procedure is verified through comparison with other procedures.