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

Hua-Xiang Zhu  Zhen-Yu Yin  Qiang Zhang 《国际地质力学数值与分析法杂志》2020,44(3):389-404

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.  相似文献   

Revisiting the existence of an effective stress for wet granular soils with micromechanics          下载免费PDF全文

Jérôme Duriez  Richard Wan  Mehdi Pouragha  Félix Darve 《国际地质力学数值与分析法杂志》2018,42(8):959-978

A possible effective stress variable for wet granular materials is numerically investigated based on an adapted discrete element method (DEM) model for an ideal three‐phase system. The DEM simulations consider granular materials made of nearly monodisperse spherical particles, in the pendular regime with the pore fluid mixture consisting of distinct water menisci bridging particle pairs. The contact force‐related stress contribution to the total stresses is isolated and tested as the effective stress candidate for dense or loose systems. It is first recalled that this contact stress tensor is indeed an adequate effective stress that describes stress limit states of wet samples with the same Mohr‐Coulomb criterion associated with their dry counterparts. As for constitutive relationships, it is demonstrated that the contact stress tensor used in conjunction with dry constitutive relations does describe the strains of wet samples during an initial strain regime but not beyond. Outside this so‐called quasi‐static strain regime, whose extent is much greater for dense than loose materials, dramatic changes in the contact network prevent macroscale contact stress‐strain relationships to apply in the same manner to dry and unsaturated conditions. The presented numerical results also reveal unexpected constitutive bifurcations for the loose material, related to stick‐slip macrobehavior.  相似文献   

用双三次样条函数和GPS资料反演现今中国大陆构造形变场   总被引：38，自引：20，他引：18  

杨少敏  游新兆  杜瑞林  乔学军  王琪 《大地测量与地球动力学》2002,22(1):68-75

将中国大陆现今构造变动视为一种连续的地壳变形,利用双三次样条函数模拟了近期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.  相似文献   

On the equivalence of two methods for computing partial derivatives of seismic waveforms—II. Laterally homogeneous initial model     

Robert J. Geller  Tatsuhiko Hara  Seiji Tsuboi 《Geophysical Journal International》1990,102(2):499-502

  相似文献   

A continuum‐discrete hydromechanical analysis of granular deposit liquefaction     

M. Zeghal  U. El Shamy 《国际地质力学数值与分析法杂志》2004,28(14):1361-1383

A coupled continuum‐discrete hydromechanical model was employed to analyse the liquefaction of a saturated loose deposit of cohesionless particles when subjected to a dynamic base excitation. The pore fluid flow was idealized using averaged Navier–Stokes equations and the discrete element method was employed to model the solid phase particles. A well established semi‐empirical relationship was utilized to quantify the fluid–particle interactions. The conducted simulations revealed a number of salient micro‐mechanical mechanisms and response patterns associated with the deposit liquefaction. Space and time variation of porosity was a major factor which affected the coupled response of the solid and fluid phases. Pore fluid flow was within Darcy's regime. The predicted response exhibited macroscopic patterns consistent with experimental results and case histories of the liquefaction of granular soil deposits. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

An explicit representation of steady state response of a beam on an elastic foundation to moving harmonic line loads     

Lu Sun 《国际地质力学数值与分析法杂志》2003,27(1):69-84

A closed‐form deflection response of a beam rest is presented in this paper using the integral transform method. The theory of linear partial differential equations is used to represent the deflection of beam subjected to a moving harmonic line load in integration form. The solution is finally carried out using the inverse Fourier transform. To evaluate the integration analytically, poles of the integrand are identified with the help of algebraic equation theory. Residue theorem is then utilized to represent the integration as a contour integral in the complex plane. Closed‐form deflections and numerical results are provided for different combinations of load frequency and velocity. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

小波分析理论在GPS技术中的应用     

ZHENG Zuoya HUANG Cheng LU Xiushan 《中国科学院上海天文台年刊》2003,(1)

在GPS数据处理中 ,存在着误差影响、影响波的干扰、周跳和数据量大等问题。误差影响和影响波的干扰实质是在接收卫星信号时受到其它因素的影响 ;周跳是由于卫星信号的失锁而造成信号的不连续 ;数据量大是因为GPS观测需要采样间隔小又连续观测所致。由于小波理论具有时频分析、波形分解、特征提取和快速小波变换等特性 ,应用小波变换和波形分解可以解决误差影响和影响波的干扰的问题 ;应用特征提取可以解决周跳检测问题 ;应用快速小波变换可进行数据压缩  相似文献   

Bond rolling resistance and its effect on yielding of bonded granulates by DEM analyses     

M. J. Jiang  H. S. Yu  D. Harris 《国际地质力学数值与分析法杂志》2006,30(8):723-761

A discrete element modelling of bonded granulates and investigation on the bond effect on their behaviour are very important to geomechanics. This paper presents a two‐dimensional (2‐D) discrete element theory for bonded granulates with bond rolling resistance and provides a numerical investigation into the effect of bond rolling resistance on the yielding of bonded granulates. The model consists of mechanical contact models and equations governing the motion of bonded particles. The key point of the theory is that the assumption in the original bond contact model previously proposed by the authors (55th CSCE‐ASCE Conference, Hamilton, Ont., Canada, 2002; 313–320; J. Eng. Mech. (ASCE) 2005; 131 (11):1209–1213) that bonded particles are in contact at discrete points, is here replaced by a more reliable assumption that bonded particles are in contact over a width. By making the idealization that the bond contact width is continuously distributed with the normal/tangential basic elements (BE) (each BE is composed of spring, dashpot, bond, slider or divider), we establish a bond rolling contact model together with bond normal/tangential contact models, and also relate the governing equations to local equilibrium. Only one physical parameter β needs to be introduced in the theory in comparison to the original bond discrete element model. The model has been implemented into a 2‐D distinct element method code, NS2D. Using the NS2D, a total of 86 1‐D, constant stress ratio, and biaxial compressions tests have been carried out on the bonded granular samples of different densities, bonding strengths and rolling resistances. The numerical results show that: (i) the new theory predicts a larger internal friction angle, a larger yielding stress, more brittle behaviour and larger final broken contact ratio than the original bond model; (ii) the yielding stress increases nonlinearly with the increasing value of β, and (iii) the first‐yield curve (initiation of bond breakage), which define a zone of none bond breakage and which shape and size are affected by the material density, is amplified by the bond rolling resistance in analogous to that predicted by the original bond model. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

富氏谱分析法模拟地壳垂直运动的地学依据及其策应力模型     

张勤  范一中 《大地测量与地球动力学》2001,21(4):76-80

应用地壳波浪与镶嵌构造学说对富氏谱分析法提取地壳垂直形变信息的科学性做了地质学意义上的阐释 ,并提出了根据多期形变资料提取特定波段上构造策应力的数学模型  相似文献   

The similitude equation in magnetotelluric inversion     

S. E. Dosso  D. W. Oldenburg 《Geophysical Journal International》1991,106(2):507-509

  相似文献