A novel coupled FDM-DEM modelling method for flexible membrane boundary in laboratory tests

Difficulties are involved in discrete element method (DEM) modelling of the flexible boundary, that is, the membranes covering the soil sample, which can be commonly found in contemporary laboratory soil tests. In this paper, a novel method is proposed wherein the finite difference method (FDM) and DEM are coupled to simulate the rubber membrane and soil body, respectively. Numerical plane strain and triaxial tests, served by the flexible membrane, are implemented and analysed later. The effect of the membrane modulus on the measurement accuracy is considered, with analytical formulae derived to judge the significance of this effect. Based on an analysis of stress-strain responses and the grain rotation field, the mechanical performances produced by the flexible and rigid lateral boundaries are compared for the plane strain test. The results show that (1) the effect of the membrane on the test result becomes more significant at larger strain level because the membrane applies additional lateral confining pressure to the soil body; (2) the tested models reproduce typical stress and volumetric paths for specimens with shear bands; (3) for the plane strain test, the rigid lateral boundary derives a much higher peak strength and larger bulk dilatation, but a similar residual strength, compared with the flexible boundary. The latter produces a more uniform (or ‘diffuse') rotation field and more mobilised local kinematics than does the former. All simulations show that the proposed FDM-DEM coupling method is able to simulate laboratory tests with a flexible boundary membrane.     

Dynamic modelling of retrogressive landslides with emphasis on the role of clay sensitivity

This paper presents a detailed numerical study of the retrogressive failure of landslides in sensitive clays. The dynamic modelling of the landslides is carried out using a novel continuum approach, the particle finite element method, complemented with an elastoviscoplastic constitutive model. The multiwedge failure mode in the collapse is captured successfully, and the multiple retrogressive failures that have been widely observed in landslides in sensitive clays are reproduced with the failure mechanism, the kinematics, and the deposition being discussed in detail. Special attention has been paid to the role of the clay sensitivity on each retrogressive failure as well as on the final retrogression distance and the final run‐out distance via parametric studies. Moreover, the effects of the viscosity of sensitive clays on the failure are also investigated for different clay sensitivities.     

对青藏东北缘现今块体划分、运动及变形的初步研究

利用2维非连续变形分析方法（DDA），以位移代替围压作为边界约束力，研究青藏东北缘现今块体划分及其运动变形。根据该地区地质构造及地震活动，以GPS点测量位移作为模拟结果约束点，得出了较合理的块体划分模型和随时间演化的主应变分布图，并把应变高值区与近几年来发生的5级以上地震作对比，得出了研究区内地震危险性可能较大的区域。另外，对模拟的甘青块体与阿拉善块体的边缘带断裂左旋运动做了大概计算。     

莱州湾南岸潍河下游地区咸水入侵灾害成因及特征

通过对莱州湾南岸咸水入侵较严重地潍河下游地区晚更新世以来沉积特征及现代自然环境条件变化的分析，探讨了沉积相对咸水入侵产生及空间范围特征的环境机理。晚更新世以来的三次海平面升降变化造成了潍河下游地区海陆沉积环境交替，形成了巨厚的海陆交互相沉积层。海进时期，大面积的滨海平原被淹没，在近海平原洼地滞留的海水经过蒸发、浓缩变为卤水，成为咸水入侵的物源；海退后陆源碎屑在滨海地区沉积形成了巨厚的古河道砂层。20世纪70年代末期以来，随着对地下淡水的过度开采，淡咸水水头压力差减小．卤水通过古河道砂层快速南侵。通过对潍河下游地区100余个地质钻孔水化学连续监测资料分析，阐明了咸水入侵的特征，有针对性地提出了咸水入侵的防治措施。     

P‐wave transmission across fractures with nonlinear deformational behaviour

Stress wave attenuation across fractured rock masses is a great concern of underground structure safety. When the wave amplitude is large, fractures experience nonlinear deformation during the wave propagation. This paper presents a study on normal transmission of P‐wave across parallel fractures with nonlinear deformational behaviour (static Barton–Bandis model). The results show that the magnitude of transmission coefficient is a function of incident wave amplitude, nondimensional fracture spacing and number of fractures. Two important indices of nondimensional fracture spacing are identified, and they divide the area of nondimensional fracture spacing into three parts (individual fracture area, transition area and small spacing area). In the different areas, the magnitude of transmission coefficient has different trends with nondimensional fracture spacing and number of fractures. In addition, the study reveals that under some circumstances, the magnitude of transmission coefficient increases with increasing number of fractures, and is larger than 1. Copyright © 2006 John Wiley & Sons, Ltd.     

Smeared crack approach: back to the original track

This paper briefly reviews the formulations used over the last 40 years for the solution of problems involving tensile cracking, with both the discrete and the smeared crack approaches. The paper focuses on the smeared approach, identifying as its main drawbacks the observed mesh‐size and mesh‐bias spurious dependence when the method is applied ‘straightly’. A simple isotropic local damage constitutive model is considered, and the (exponential) softening modulus is regularized according to the material fracture energy and the element size. The continuum and discrete mechanical problems corresponding to both the weak discontinuity (smeared cracks) and the strong discontinuity (discrete cracks) approaches are analysed and the question of propagation of the strain localization band (crack) is identified as the main difficulty to be overcome in the numerical procedure. A tracking technique is used to ensure stability of the solution, attaining the necessary convergence properties of the corresponding discrete finite element formulation. Numerical examples show that the formulation derived is stable and remarkably robust. As a consequence, the results obtained do not suffer from spurious mesh‐size or mesh‐bias dependence, comparing very favourably with those obtained with other fracture and continuum mechanics approaches. Copyright © 2006 John Wiley & Sons, Ltd.     

昆仑山8.1级地震前中国大陆的构造应变背景

利用“网络工程”1998～2001年累积的1181个测站的GPS重复观测资料，采用双三次样条函数模型建立中国大陆水平运动模型速度场，用大地坐标在椭球面上计算各类应变场，详细分析了2001年昆仑山8．1级地震前中国大陆水平构造应变场空间分布特征。各类构造应变场的最高值都出现在喜马拉雅构造带与昆仑山地块内(地震断裂带南侧)，鲜水河—安宁河断裂带次之。分析表明，昆仑山8．1级地震正好发生在张性面膨胀应变率的高值区，第一、第二和最大剪应变率高值区边缘的突变区和最大、最小主应变率的高值区。     

Strain localization in sand: an overview of the experimental results obtained in Grenoble using stereophotogrammetry

Experimental results are presented from the extensive program of drained plane strain compression tests on sand carried out in Grenoble over the last two decades. Systematic analysis of photographs of the deforming specimen allowed for measuring deformations and determining strain fields throughout the test, that is: prior to, at, and after the onset of strain localization. The principles, details and accuracy of the procedure are described, as well as its suitability to properly depict the patterns of deformation. Findings concerning the occurrence and progression of strain localization are discussed. The issues of shear band orientation and thickness are addressed, as well as temporary and persistent complex localization patterns, and the volumetric behaviour inside a band after its formation. The influence of such variables as initial state of the sand (effective stress and relative density), specimen size and slenderness, as well as grain size, is discussed. Copyright © 2004 John Wiley & Sons, Ltd     

地震前后垂直形变场动态演化的量化指标

提出一种表达垂直形变场动态演化过程的量化指标——区域应变率、应变集中度。在此基础上，对南北地震带各水准监测区近30年的垂直形变资料进行了实际计算．并结合具体震例进行了对比研究。结果表明：量化指标在一定程度上能够反映地震前后形变场的动态演化过程和地壳运动状态，对地震的中长期预报有一定的积极意义。     

广西老厂铅锌矿床地质特征及其成因探讨

本文是在专题研究基础上,探讨广西老厂铅锌矿床地质特征和矿床成因。该矿床位于海洋山复式背斜南段的老厂穹窿构造核部,铅锌矿化受构造、地层和岩性控制。矿体呈脉状,产于NE向组张扭性断裂构造中;围岩为寒武系白云岩,板岩、砂岩;矿石类型以锌为主,矿物成分简单,以闪锌矿、方铅矿与石英、白云石为主,具充填、交代矿石结构、构造,矿石平均品位含Zn9.94%、Pb1.50%;围岩蚀变主要为硅化、绢云母化和碳酸盐化。据地层岩石元素丰度值、稀土元素配分和稳定同位素数据、资料,探讨了成矿物质来源于围岩,热液主要来自大气降水。因此,认为本矿床实属产于寒武系轻变质岩中的层控矿型大气降水热液铅锌矿床。