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.     

Developing constitutive models from EPR‐based self‐learning finite element analysis

A constitutive model that captures the material behavior under a wide range of loading conditions is essential for simulating complex boundary value problems. In recent years, some attempts have been made to develop constitutive models for finite element analysis using self‐learning simulation (SelfSim). Self‐learning simulation is an inverse analysis technique that extracts material behavior from some boundary measurements (eg, load and displacement). In the heart of the self‐learning framework is a neural network which is used to train and develop a constitutive model that represents the material behavior. It is generally known that neural networks suffer from a number of drawbacks. This paper utilizes evolutionary polynomial regression (EPR) in the framework of SelfSim within an automation process which is coded in Matlab environment. EPR is a hybrid data mining technique that uses a combination of a genetic algorithm and the least square method to search for mathematical equations to represent the behavior of a system. Two strategies of material modeling have been considered in the SelfSim‐based finite element analysis. These include a total stress‐strain strategy applied to analysis of a truss structure using synthetic measurement data and an incremental stress‐strain strategy applied to simulation of triaxial tests using experimental data. The results show that effective and accurate constitutive models can be developed from the proposed EPR‐based self‐learning finite element method. The EPR‐based self‐learning FEM can provide accurate predictions to engineering problems. The main advantages of using EPR over neural network are highlighted.     

RECENT TECTONIC DEFORMATION ANOMALY AND EARTHQUAKES IN GANSU-NINGXIA-QINGHAI AREA

1　IntroductionAccordingtogeologicalstructure ,theGansu Ningxia Qinghaiareabelongstothenortheastmar ginofQinghai Tibetblock .Thisareahasbeenpaidmuchattentionby geo specialistsinChinaandabroadbecauseofitssignificanttectonicmovement,itsintensiveseismicity ,anditsimportanceinearth quakehazardmitigation .IntheDevelopmentPro gramonNationalKeyBasicResearchesundertheProject“MechanismandPredictionofContinentalStrongEarthquakes”,themechanismsofcontinen talstrongearthquakesarestudied ,usinghypo…     

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

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

用双三次样条函数和GPS资料反演现今中国大陆构造形变场

将中国大陆现今构造变动视为一种连续的地壳变形,利用双三次样条函数模拟了近期GPS测定的大陆内部及周边地区412个测站速率,反演大陆地区自洽的构造变动速度场和应变率场.模拟结果显示:印度板块与欧亚板块的碰撞、挤压是构成中国大陆内部岩石层水平形变的主要驱动力.印度板块在东喜马拉雅构造结深深插入青藏高原,造成地壳大规模的缩短和抬升.青藏高原东南部的喜马拉雅带、拉萨和羌塘地块以及青藏高原东南边的川滇地区,内部构造活动强烈,其内部的构造变形包含地壳碎片的冲断、褶皱和侧向逃逸.大陆地壳(或岩石圈)的增厚,尤其是喜马拉雅山脉南北向的快速缩短和青藏高原东西向的缓慢拉张,大约吸收了印欧板块会聚量的85%,西藏中东地区东西向的拉张速率达到了(16±2.0)mm/a,且顺时针方向扭转明显.印度板块相对欧亚板块运动的欧拉极为(29.7°N, 19.3°E, 0.392°/Ma);华南地块相对于欧亚大陆向东(102°±7.4°)南的运动速率是(11±1.54)mm/a,华南块体相对欧亚板块运动的欧拉极为(62.25°N, 126.56°E, 0.141°/Ma);塔里木地块相对较稳定,其西部运动速度高于东部运动速度,作顺时针方向旋转.总体上讲,中国大陆运动方向为北偏东呈辐射状,从西部近南北方向的运动转向东部地区东南方向的运动,绕东喜马拉雅构造结有一顺时针方向的旋转.横穿喜马拉雅构造带及青藏内部的南北向压缩速率为(19±2.0)mm/a,横穿西天山构造带的南北向压缩平均速率为(13±1.5)mm/a,横穿东天山构造带的南北向压缩平均速率为(6.0±1.4)mm/a.阿尔金断裂带的左旋走滑速率为(6±1.2)mm/a.     

Application of micropolar plasticity to post failure analysis in geomechanics

A micropolar elastoplastic model for soils is formulated and a series of finite element analyses are employed to demonstrate the use of a micropolar continuum in overcoming the numerical difficulties encountered in application of finite element method in standard Cauchy–Boltzmann continuum. Three examples of failure analysis involving a deep excavation, shallow foundation, and a retaining wall are presented. In all these cases, it is observed that the length scale introduced in the polar continuum regularizes the incremental boundary value problem and allows the numerical simulation to be continued until a clear collapse mechanism is achieved. The issue of grain size effect is also discussed. Copyright © 2004 John Wiley & Sons, Ltd.     

青藏块体东北缘断层形变与中强地震

对祁连山－海原断裂带近期断层形变特征进行了初步研究，发现多场地，大范围的断层活动异常是中等强度地震发生的显著背景，并且往往与大陆地震活动的阶段性总体状况相呼应；区域形变存在明显的特征量，包括特征形态和特征时间，同一场地在不同地震前的异常特征具有重复性，但会受到背景差异显著的不同地震的影响，目前形变状况表明研究区仍具有发生中强地震的地壳运动背景。     

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

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

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

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

描述地壳垂直运动过程的一种函数形式

本文在已有数据处理方法的基础上，利用近代数值逼近理论，给出了从时空域角度描述地壳垂直运动过程的一种具体的函数解析形式。最后给出了一个实际算例。