共查询到20条相似文献，搜索用时 46 毫秒
1.
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. 相似文献
2.
This paper describes the mechanical behavior of artificially cemented sands with strong, intermediate, and weak bond strengths, using experimentation and 3D discrete element method (DEM) simulation. The focus is on the features of bond breakage and the associated influences on the stress–strain responses. Under triaxial shearing, the acoustic emission rate captured in the experiment and the bond breakage rate recorded in the simulations show resemblance to the stress–strain response, especially for strongly and intermediately cemented samples, where a strain softening response is observed. The simulations further reveal the shear band formation coincides with the development of bond breakage locations due to the local weakness caused by the bond breakages. Strain softening and volumetric dilation are observed inside the shear band, while the region outside the shear band undergoes elastic unloading. The weakly cemented sample exhibits a strain hardening response instead; bond breakages and the associated local weaknesses are always randomly formed such that no persistent shear band is observed. Note that in the DEM simulation, the flexible membrane boundary is established by a network of bonded membrane particles; the membrane particle network is further partitioned into finite triangular elements. The associated algorithm can accurately distribute the applied confining pressure onto the membrane particles and determine the sample volume. 相似文献
3.
This article presents a fundamental study on the role of particle breakage on the shear behavior of granular soils using the three‐dimensional (3‐D) discrete element method. The effects of particle breakage on the stress ratio, volumetric strain, plastic deformation, and shear failure behavior of dense crushable specimens undergoing plane strain shearing conditions are thoroughly investigated through a variety of micromechanical analyses and mechanism demonstrations. The simulation of a granular specimen is based on the effective modeling of realistic fracture behavior of single soil particles, which is demonstrated by the qualitative agreement between the results from platen compression simulations and those from physical laboratory tests. The simulation results show that the major effects of particle breakage include the reduction of volumetric dilation and peak stress ratio and more importantly the plastic deformation mechanisms and the shear failure modes vary as a function of soil crushability. Consistent macro‐ and micromechanical evidence demonstrates that shear banding and massive volumetric contraction depict the two end failure modes of a dense specimen, which is dominated by particle rearrangement–induced dilation and particle crushing–induced compression, respectively, with a more general case being the combination and competition of the two failure modes in the medium range of soil crushability and confining stress. However, it is further shown that a highly crushable specimen will eventually develop a shear band at a large strain because of the continuous decay of particle breakage. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
4.
深部岩石在工程中具有高应力、大变形等典型特点,因此,高围压下考虑岩石残余强度的三轴试验对于分析深部岩石力学特性具有重要意义。离散单元法是分析岩石力学特性的重要数值方法,但是长期以来采用离散单元法定量模拟岩石的三轴试验一直存在诸多挑战,即数值模拟与室内试验得到的应力-应变全过程曲线难以定量匹配。采用改进的三维胶结抗弯-扭模型对深部砂岩考虑残余强度时的三轴试验进行了定量模拟,实现了数值模拟与室内试验应力-应变全过程曲线的定量匹配,获得了岩石较大的峰值/残余内摩擦角及非线性强度包线,克服了经典BPM模型存在的3个突出问题。通过参数分析,研究了峰值/残余内摩擦角及黏聚力与离散元微观参数之间的关系,同时这些大量的算例也证明了该模型具有较高的计算效率,可以满足模拟三维室内常规试验的要求。 相似文献
5.
粗粒料是一定级配的岩石颗粒集合体,具有独特的物理力学特性。以粗粒料室内三轴固结排水试验成果为基础,基于离散元颗粒流理论,从细观角度出发,以PFC3D为工具,通过自编程及二次开发,得到按级配生成的粗粒料三轴试验数值模型。引入clump颗粒考虑颗粒形状对粗粒料强度及变形的影响,分析剪胀、颗粒形状、颗粒重排的关系。结果表明：颗粒形状是影响粗粒料强度与变形的主要因素,在其他细观参数一定的情况下,改变颗粒形状,可以显著影响粗粒料的力学行为;BPM模型的应力-应变关系只在低围压下与试验值吻合,随着围压的增大,偏差越来越大;而引入clump颗粒的PFC3D数值模型能很好地模拟粗粒料室内三轴固结排水试验的应力-应变特性,但由于BPM及clump都是刚性颗粒,没有考虑颗粒变形及破碎,造成应变剪胀偏大。 相似文献
6.
A simple three‐dimensional distinct element modeling of the mechanical behavior of bonded sands 
下载免费PDF全文
下载免费PDF全文 This paper presents a simple three‐dimensional (3D) Distinct Element Method (DEM) for numerical simulation of the mechanical behavior of bonded sands. First, a series of micro‐mechanical tests on a pair of aluminum rods glued together by cement with different bond sizes were performed to obtain the contact mechanical responses of ideally bonded granular material. Second, a 3D bond contact model, which takes into account the influences of bond sizes, was established by extending the obtained 2D experimental results to 3D case. Then, a DEM incorporating the new contact model was employed to perform a set of drained triaxial compression tests on the DEM bonded specimens with different cement contents under different confining pressures. Finally, the mechanical behavior of the bonded specimens was compared with the available experimental results. The results show that the DEM incorporating the simple 3D bond contact model is able to capture the main mechanical behavior of bonded sands. The bonded specimen with higher cement content under lower confining pressure exhibits more pronounced strain softening and shear dilatancy. The peak and residual strengths, the apparent cohesion and peak/residual friction angles, and the position and slope of the critical state line increase with increase in cement content. Microscopically, bond breakage starts when the system starts to dilate and the maximum rate of bond breakage coincides with the maximum rate of dilation. Bond breakage is primarily due to tension‐shear failure and the percentage of such failures is independent of both confining pressure and cement content. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
7.
Multi‐scale investigations aided by the discrete element method (DEM) play a vital role for current state‐of‐the‐art research on the elementary behaviour of granular materials. Similar to laboratory tests, there are three important aspects to be considered carefully, which are the proper stress/strain definition and measurement, the application of target loading paths and the designed experiment setup, to be addressed in the present paper. Considering the volume sensitive characteristics of granular materials, in the proposed technique, the deformation of the tested specimen is controlled and measured by deformation gradient tensor involving both the undeformed configuration and the current configuration. Definitions of Biot strain and Cauchy stress are adopted. The expressions of them in terms of contact forces and particle displacements, respectively, are derived. The boundary of the tested specimen consists of rigid massless planar units. It is suggested that the representative element uses a convex polyhedral (polygonal) shape to minimize possible boundary arching effects. General loading paths are described by directly specifying the changes in the stress/strain invariants or directions. Loading can be applied in the strain‐controlled mode by specifying the translations and rotations of the boundary units, or in the stress‐controlled mode by using a servo‐control mechanism, or in the combination of the two methods to realize mixed boundary conditions. Taking the simulation results as the natural consequences originated from a complex system, virtual experiments provide particle‐scale information database to conduct multi‐scale investigations for better understanding in granular material behaviours and possible development of the constitutive theories provided the qualitative similarity between the simulation results from virtual experiments and observations on real material behaviour. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
8.
In this paper, numerical simulation of 3-dimensional assemblies of 1000 polydisperse sphere particles using Discrete Element
Method (DEM) is used to study the liquefaction behaviour of granular materials. Numerical simulations of cyclic triaxial shear
tests under undrained conditions are performed at different confining pressures under constant strain amplitude. Results obtained
in these numerical simulations indicate that with increase in confining pressure there is an increase in liquefaction resistance. 相似文献
9.
The paper presents Cauchy stress tensor computation over parallel grids of message passing interface (MPI) parallel three-dimensional (3D) discrete element method (DEM) simulations of granular materials, considering spherical and nonspherical particles. The stress tensor computation is studied for quasi-static and dynamic conditions, and its resulting symmetry or asymmetry is discussed within the context of classical continuum mechanics (CCM), granular materials mechanics (GMM), and micropolar continuum mechanics (MCM). The average Cauchy stress tensor computation follows Bagi's and Nicot's formulations and is verified within MPI parallel 3D DEM simulations involving dynamically adaptive compute grids. These grids allow calculation of temporal and spatial distributions of stress across granular materials under static and dynamic conditions. The vertical stress component in gravitationally deposited particle assemblies exhibits nonuniform spatial distributions under static equilibrium, and its zone of maximum value changes during the process of gravitational pluviation and collapse. These phenomena reveal a microstructural effect on stress distribution within granular materials that is attributed to their discrete particulate nature (particle size, shape, gradation, boundary conditions, etc). 相似文献
10.
为了解黄土三轴剪切过程中的微观变化,采用PFC3D建立黄土三轴试验模型,模拟围压分别为0、50、150、300 kPa的三轴剪切试验,并与室内试验进行对比分析。分析结果表明：PFC3D能够较好地模拟出不同围压下从开始到破坏到残余变形整个过程中应力-应变的变化规律,且位移场及接触应力场的变化规律与室内三轴试验宏观现象较一致;发现弹性模量、泊松比及峰值强度与数值模型中微观参数有着密切的联系,如法向刚度kn控制试样宏观弹性模量,kn /ks值控制泊松比,摩擦系数控制峰值强度;通过体应变-轴向应变曲线发现,随着围压的增加应变能增大,试样呈现出由体积膨胀到体积减缩的变化规律。其研究结果为进一步探究黄土的应力-应变性状及抗剪强度特性提供参考。 相似文献
11.
为分析砂土在复杂应力条件下的剪切力学特性,采用商业离散元软件PFC3D对单粒组中密砂的空心扭剪试验进行了仿真模拟,分析了数值试样的应力-应变关系,研究了不同剪切方向下离散介质的强度、体积应变特性以及中主应力比对它们的影响,再现了力链在加载过程中的演化,并对剪切带的倾角做了深入分析。同时,从细观上看,以颗粒接触数和纯转动率变量为中心,观察了试样内部颗粒的运动状态,对比了不同剪切方向下剪切带内外颗粒接触数与纯转动位移的变化。最后,将数值试验结果与已有的室内试验结果进行了对比。此研究实现了复杂应力条件下空心扭剪试样的三维离散元模拟,加深了对空心扭剪试验过程和结果的理解和解释。 相似文献
12.
颗粒级配对粗粒土强度与变形特性影响的研究 总被引：1，自引：0，他引：1
采用三维颗粒流软件PFC3D的内置FISH语言,进行二次开发,模拟粗粒土级配,建立一个研究高铁填方路基粗粒土变形特性的颗粒流模型。通过粗粒土三轴试验,标定了表征粗粒土细观力学性质的模型参数,并验证了模型的有效性。考虑粗粒土曲率系数和不均匀系数对颗粒级配的影响,模拟并研究了5种不同级配下粗粒土在3种不同围压下的变形特性。试验结果表明,曲率系数对粗颗粒土变形特性影响较大,当其有效粒径与中值粒径及限制粒径的差距过大时,粗粒土中颗粒填充及致密性将变差,变形增大,强度降低。以较大曲率系数与较大不均匀系数组合的粗粒土受力性能较好,但此种级配的粗粒土压缩性大。级配良好且不均匀系数较大,其应力链分布越均匀,颗粒的挤压效果越显著,力的传递和分配也越均匀。 相似文献
13.
三维离散元单轴试验模拟甲烷水合物宏观三轴强度特性 总被引：1，自引：0，他引：1
填充型水合物的砂性能源土试样可视为特殊的散粒体材料,即砂粒和水合物颗粒混合物,具有明显的非连续特征。为研究填充型水合物的能源土力学特性,初步探索了甲烷水合物在不同温度、反压条件下加荷模式的离散元模拟方法。离散元模拟中,将水合物块体视为由大量颗粒通过强胶结作用凝聚而成的整体,室内试验中的内部孔隙水压作用转化为水合物颗粒间的胶结力,故需要合理确定颗粒间胶结模型参数来实现反压的影响作用。通过参数反演建立了宏观强度、刚度参数与平行胶结模型的微观胶结参数间的宏、微观关系,基于已有室内甲烷水合物三轴试验资料,确定了给定温度和反压条件下的微观胶结参数取值,随后进行离散元单轴压缩试验。离散元单轴压缩试验模拟获得的水合物强度特性,与室内三轴试验结果符合较好;通过建立的宏、微观参数间的关系,实现了不同温度、反压下的水合物加荷模式的模拟。为进一步提出深海能源土离散元数值试验成样方法--孔隙填充水合物生成技术,形成含填充型水合物的能源土试样,研究其力学和变形特性奠定基础。 相似文献
14.
The discrete element method has been used to investigate the micro mechanics of cemented sand. High‐pressure drained triaxial tests are modelled in 3D using a flexible membrane that allows the correct deformation to develop. Simulations with up to 12 MPa confining pressure are presented, which are compared with laboratory experiments on a sand with a range of cement contents. Cementation is modelled using ‘parallel bonds’, and various parameters and strength distributions are investigated. Varying levels of cementation are successfully modelled, with the correct qualitative behaviour observed and the separate effects of cementation and confining pressures demonstrated. The triaxial behaviour is found to be highly influenced by the distribution of bond strengths. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
15.
The paper provides an in-depth exploration of the role of particle crushing on particle kinematics and shear banding in sheared granular materials. As a two-dimensional approximation, a crushable granular material may be represented by an assembly of irregularly shaped polygons to include shape diversity of realistic granular materials. Particle assemblies are subjected to biaxial shearing under flexible boundary conditions. With increasing percentage of crushed particles, mesoscale deformation becomes increasingly unstable. Fragmented deformation patterns within the granular assemblies are unable to form stable and distinct shear bands. This is confirmed by the sparsity of large fluctuating velocities in highly crushable assemblies. Without generating distinct shear bands, deformation patterns and failure modes of a highly crushable assembly are similar to those of loose particle assemblies, which are regarded as diffuse deformation. High degrees of spatial association amongst the kinematical quantities confirm the key role that non-affine deformation and particle rotation play in the generation of shear bands. Therefore, particle kinematical quantities can be used to predict the onset and subsequent development of shear zones, which are generally marked by increased particle kinematic activity, such as intense particle rotation and high granular temperature. Our results indicate that shear band thickness increases, and its speed of development slows down, with increasing percentage of crushed particles. As particles crush, spatial force correlation becomes weaker, indicating a more diffuse nature of force transmission across particle contacts. 相似文献
16.
We present a multiscale investigation on the initiation and development of compaction bands in high-porosity sandstones based on an innovative hierarchical multiscale approach. This approach couples the finite element method and the discrete element method (DEM) to offer direct, rigorous linking of the microscopic origins and mechanisms with complex macroscopic phenomena observed in granular rocks such as strain localization and failure. To simulate compaction band in granular cementitious sandstone, we adopt a bonded contact model with normal and tangential interparticle cohesions in the DEM and propose a dual-porosity structure consisting of macro-pores and interstitial voids for the representative volume element to mimic the typical meso-structure of high-porosity sandstones. In the absence of particle crushing, our multiscale analyses identify debonding and pore collapses as two major contributors to the formation of compaction bands. The critical pressures predicted by our simulations, corresponding to surges of debonding and pore collapse events, agree well with the estimations from field data. The occurrence patterns of compaction band are found closely related to specimen heterogeneity, porosity and confining pressure. Other deformation band patterns, including shear-enhanced compaction bands and compactive shear bands, were also observed under relatively low confining pressure conditions with a rough threshold at \(0.55P^{*}\) (\(P^{*}\) is the critical pressure) on the failure envelop. Key microscopic characteristics attributable to the occurrence of these various deformation patterns, including fabric anisotropy, particle rotation, debonding and pore collapse, are examined. Shear-enhanced compaction bands and pure compaction bands bear many similarities in terms of these microscopic characteristics, whereas both differ substantially from compactive shear bands. 相似文献
17.
Validation of discrete element method by simulating a 2D assembly of randomly packed elliptical rods
This paper aims at establishing the predictive capability of the discrete element method (DEM) by validating the simulated responses of granular systems against experimental observations at both the macroscale and the microscale. A previously published biaxial shearing test on a 2D assembly of randomly packed elliptical rods was chosen as the benchmark test. In carrying out the corresponding DEM simulations herein, the contact model was derived and then validated using finite element analysis; the associated parameters were calibrated experimentally. The flexible (membrane) boundary was modeled by a bonded-particle string with experimentally calibrated parameters. An iteration procedure was implemented to replicate the initial packing and also to satisfy the boundary conditions in the experiment. Overall, the DEM simulation is found effective in reproducing the stress–strain–volumetric response, the statistical observation on the fabric anisotropy and the strain localization. Furthermore, the closer the numerical packing is to the experimental one, the closer the response is reproduced, demonstrating the significance of the initial packing reconstruction. Still, there are some minor differences between the experiment and simulation, reflecting the limitations associated with the particle number and the measurement resolution used in the experiment when reproducing the initial packing. 相似文献
18.
Shaohan Zhao Ha H. Bui Vincent Lemiale Giang D. Nguyen Felix Darve 《国际地质力学数值与分析法杂志》2019,43(5):1005-1031
In this paper, a new approach to applying confining stress to flexible boundaries in the smoothed particle hydrodynamics (SPH) method is developed to facilitate its applications in geomechanics. Unlike the conventional SPH methods that impose confining boundary conditions by creating extra boundary particles, the proposed approach makes use of kernel truncation properties of SPH approximations that occur naturally at free-surface boundaries. Therefore, it does not require extra boundary particles and, as a consequence, can be utilised to apply confining stresses onto any boundary with arbitrary geometry without the need for tracking the curvature change during the computation. This enables more complicated problems that involve moving confining boundaries, such as confining triaxial tests, to be simulated in SPH without difficulties. To further enhance SPH applications in elasto-plastic computations of geomaterials, a robust numerical procedure to implement Mohr-Coulomb plasticity model in SPH is presented for the first time to avoid difficulties associated with corner singularities in Mohr-Coulomb model. The proposed approach was first validated against two-dimensional finite element (FE) solutions for confining biaxial compression tests to demonstrate its predictive capability at small deformation range when FE solutions are still valid. It is then further extended to three-dimensional conditions and utilised to simulate triaxial compression experiments. Simulation results predicted by SPH show good agreement with experiments, FE solutions, and other numerical results available in the literature. This suggests that the proposed approach of imposing confining stress boundaries is promising and can handle complex problems that involve moving confining boundary conditions. 相似文献
19.
The axial stress–strain relations of embedded granular columns encapsulated with flexible reinforcement were evaluated using an analytical procedure based on the cavity expansion method. This proposed method has firstly been verified through an experimental triaxial test on a reinforced sand specimen. A normalized relation was established between the volumetric change and the axial strain of soil to enable the analysis of granular material behavior under a continuous increase in lateral pressure. The analytical results show that the reinforced granular columns embedded in clay behave differently from granular columns subjected only to a constant confining pressure. It is found that reinforcing a column with a sleeve at the top portion will be adequate to prevent the column from bulging and also improve its load carrying capacity. The optimum skirting length that a sleeve can deter a granular column from bulging depends on the characteristics of the in situ soil and the stiffness and yield strength of the sleeve. 相似文献
20.
This paper describes a three-dimensional random network model to evaluate the thermal conductivity of particulate materials. The model is applied to numerical assemblies of poly-dispersed spheres generated using the discrete element method (DEM). The grain size distribution of Ottawa 20–30 sand is modeled using a logistic function in the DEM assemblies to closely reproduce the gradation of physical specimens. The packing density and inter-particle contact areas controlled by confining stress are explored as variables to underscore the effects of micro- and macro-scales on the effective thermal conductivity in particulate materials. It is assumed that skeletal structure of 3D granular system consists of the web of particle bodies interconnected by thermal resistor at contacts. The inter-particle contact condition (e.g., the degree of particle separation or overlap) and the particle radii determine the thermal conductance between adjacent particles. The Gauss–Seidel method allows evaluation of the evolution of temperature variation in the linear system. Laboratory measurements of thermal conductivity of Ottawa 20–30 sand corroborate the calculated results using the proposed network model. The model is extended to explore the evolution of thermal conduction depending on the nucleation habits of secondary solid phase as an anomalous material in the pore space. The proposed network model highlights that the coordination number, packing density and the inter-particle contact condition are integrated together to dominate the heat transfer characteristics in particulate materials, and allows fundamental understanding of particle-scale mechanism in macro-scale manifestation. 相似文献