Fully coupled elastoplastic hydro-mechanical analysis of unsaturated porous media using a meshfree method     

Omid Ghaffaripour  Golnaz A. Esgandani  Arman Khoshghalb  Babak Shahbodaghkhan 《国际地质力学数值与分析法杂志》2019,43(11):1919-1955

This paper presents the first application of an advanced meshfree method, ie, the edge-based smoothed point interpolation method (ESPIM), in simulation of the coupled hydro-mechanical behaviour of unsaturated porous media. In the proposed technique, the problem domain is spatially discretised using a triangular background mesh, and the polynomial point interpolation method combined with a simple node selection scheme is adopted for creating nodal shape functions. Smoothing domains are formed on top of the background mesh, and a constant smoothed strain, created by applying the smoothing operation over the smoothing domains, is assigned to each smoothing domain. The deformation and flow models are developed based on the equilibrium equation of the mixture, and linear momentum and mass balance equations of the fluid phases, respectively. The effective stress approach is followed to account for the coupling between the flow and deformation models. Further coupling among the phases is captured through a hysteretic soil water retention model that evolves with changes in void ratio. An advanced elastoplastic constitutive model within the context of the bounding surface plasticity theory is employed for predicting the nonlinear behaviour of soil skeleton. Time discretisation is performed by adopting a three-point discretisation method with growing time steps to avoid temporal instabilities. A modified Newton-Raphson framework is designed for dealing with nonlinearities of the discretised system of equations. The performance of the numerical model is examined through a number of numerical examples. The state-of-the-art computational scheme developed is useful for simulation of geotechnical engineering problems involving unsaturated soils.  相似文献   

Variations in glacier volume and snow cover and their impact on lake storage in the Paiku Co Basin,in the Central Himalayas     

Xingxing Zhang  Rui Wang  Zhijun Yao  Zhaofei Liu 《水文研究》2020,34(8):1920-1933

Glaciers and snow cover are important constituents of the surface of the Tibetan Plateau. The responses of these phenomena to global environmental changes are sensitive, rapid and intensive due to the high altitudes and arid cold climate of the Tibetan Plateau. Based on multisource remote sensing data, including Landsat images, MOD10A2 snow product, ICESat, Cryosat-2 altimetry data and long-term ground climate observations, we analysed the dynamic changes of glaciers, snow melting and lake in the Paiku Co basin using extraction methods for glaciers and lake, the degree-day model and the ice and lake volume method. The interaction among the climate, ice-snow and the hydrological elements in Paiku Co is revealed. From 2000 to 2018, the basin tended to be drier, and rainfall decreased at a rate of −3.07 mm/a. The seasonal temperature difference in the basin increased, the maximum temperature increased at a rate of 0.02°C/a and the minimum temperature decreased at a rate of −0.06°C/a, which accelerated the melting from glaciers and snow at rates of 0.55 × 10⁷ m³/a and 0.29 × 10⁷ m³/a, respectively. The rate of contribution to the lake from rainfall, snow and glacier melted water was 55.6, 27.7 and 16.7%, respectively. In the past 18 years, the warmer and drier climate has caused the lake to shrink. The water level of the lake continued to decline at a rate of −0.02 m/a, and the lake water volume decreased by 4.85 × 10⁸ m³ at a rate of −0.27 × 10⁸ m³/a from 2000 to 2018. This evaluation is important for understanding how the snow and ice melting in the central Himalayas affect the regional water cycle.  相似文献   

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

Three-dimensional finite element simulation of fracture propagation in rock specimens with pre-existing fissure(s) under compression and their strength analysis     

Ramin Pakzad  Shanyong Wang  Scott W. Sloan 《国际地质力学数值与分析法杂志》2020,44(10):1472-1494

A FORTRAN program, consistent with the commercially available finite element (FE) code ABAQUS, is developed based on a three-dimensional (3D) linear elastic brittle damage constitutive model with two damage criteria. To consider the heterogeneity of rock, the developed FORTRAN program is used to set the stiffness and strength properties of each element of the FE model following a Weibull distribution function. The reliability of the program is assessed against available experimental results for granite cylindrical specimens with a throughgoing, flat and inclined fissure. The calibration procedure of the material parameters is explained in detail, and it is shown that the compressive to tensile strength ratio can have a substantial influence on the failure response of the specimens. Numerical simulations are conducted for models with different levels of heterogeneity. The results show a smaller load bearing capacity for models with less homogeneity, representing gradual coalescence of fully damaged elements forming throughout the models during loading. The maximum load bearing capacity is studied for various combinations of inclination angles of two centrally aligned, throughgoing and flat fissures of equal length embedded in cylindrical models under uniaxial and multiaxial loading conditions. The key role of the compressive to tensile strength ratio is highlighted by repeating certain simulations with a lower compressive to tensile strength ratio. It is proven that the peak loads of the rock models with sufficiently small compressive to tensile strength ratios containing two throughgoing fissures of equal length are similar, provided that the minimum inclination angles of the models are the same. The results are presented and discussed with respect to the existing experimental findings in the literature, suggesting that the numerical model applied in this study can provide useful insight into the failure behaviour of rock-like materials.  相似文献   

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

Implementation of a coupled plastic damage distinct lattice spring model for dynamic crack propagation in geomaterials          下载免费PDF全文

Chao Jiang  Gao‐Feng Zhao 《国际地质力学数值与分析法杂志》2018,42(4):674-693

This paper studies dynamic crack propagation by employing the distinct lattice spring model (DLSM) and 3‐dimensional (3D) printing technique. A damage‐plasticity model was developed and implemented in a 2D DLSM. Applicability of the damage‐plasticity DLSM was verified against analytical elastic solutions and experimental results for crack propagation. As a physical analogy, dynamic fracturing tests were conducted on 3D printed specimens using the split Hopkinson pressure bar. The dynamic stress intensity factors were recorded, and crack paths were captured by a high‐speed camera. A parametric study was conducted to find the influences of the parameters on cracking behaviors, including initial and peak fracture toughness, crack speed, and crack patterns. Finally, selection of parameters for the damage‐plasticity model was determined through the comparison of numerical predictions and the experimentally observed cracking features.  相似文献   

基于噪声谱比法的内蒙古西部地区台站场地响应研究     

张珂  王鑫  王树波  张红才  宋晓燕 《大地测量与地球动力学》2020,40(6):646-651

利用内蒙古西部12个台站的地脉动噪声数据，采用噪声谱比法研究台站的场地响应情况。分析表明，台站场地响应按曲线形态可分为3类，且可能受地形地貌、局部构造和台基状况等条件影响。对比分析噪声谱比法与Moya方法的场地响应结果发现，场地响应曲线形态基本一致，只有极少数台站存在明显差异。结果表明，内蒙古西部区域大部分台站的台基状况总体较好，场地响应曲线较为平坦，无明显频率放大点。  相似文献   

河外H2O超脉泽源及与其相关的分子谱线的观测和研究     

俞志尧 《天文学进展》2002,20(3):256-264

自从在活动星系核NGC4945的视线方向上发现第一个河外H2O超脉泽源以来，迄今为止已发现了19个河外H2O超脉泽源，对与活动星系核成协和河外H2O超脉泽源及分子谱线的观测和研究是探测和研究活动星系核核区中央源，拱核气体和尘埃环性质的非常有效的工具，主要评述对河外H2O超脉泽源及与其相关分子谱线的搜索，观测和理论研究现状。  相似文献   

Muon density measurements with the KASCADE central detector     

《Astroparticle Physics》2002,16(4):183-386

Frequency distributions of local muon densities in high-energy extensive air showers (EAS) are presented as signature of the primary cosmic ray energy spectrum in the knee region. Together with the gross shower variables like shower core position, angle of incidence, and the shower sizes, the KASCADE experiment is able to measure local muon densities for two different muon energy thresholds. The spectra have been reconstructed for various core distances, as well as for particular subsamples, classified on the basis of the shower size ratio N_μ/N_e. The measured density spectra of the total sample exhibit clear kinks reflecting the knee of the primary energy spectrum. While relatively sharp changes of the slopes are observed in the spectrum of EAS with small values of the shower size ratio, no such feature is detected at EAS of large N_μ/N_e ratio in the energy range of 1–10 PeV. Comparing the spectra for various thresholds and core distances with detailed Monte Carlo simulations the validity of EAS simulations is discussed.  相似文献   

On the accretion geometry of Cyg X-1 in the low/hard state     

F. Eugenio Barrio  Chris Done  Sergei Nayakshin 《Monthly notices of the Royal Astronomical Society》2003,342(2):557-563

  相似文献