共查询到20条相似文献，搜索用时 734 毫秒
1.
This study aims at providing a hybrid calibration framework to estimate Hertz-type contact parameters (particle-scale shear modulus and Poisson ratio) for both two-dimensional and three-dimensional discrete element modelling (DEM). On the basis of statistically isotropic granular packings, a set of analytical formulae between macroscopic material parameters (Young modulus and Poisson ratio) and particle-scale Hertz-type contact parameters for granular systems are derived under small-strain isotropic stress conditions. However, the derived analytical solutions are only estimated values for general models. By viewing each DEM modelling as an implicit mathematical function taking the particle-level parameters as independent variables and employing the derived analytical solutions as the initial input parameters, an automatic iterative scheme is proposed to obtain the calibrated parameters with higher accuracies. Considering highly nonlinear features and discontinuities of the macro-micro relationship in Hertz-based discrete element models, the adaptive moment estimation algorithm is adopted in this study because of its capacity of dealing with noise gradients of cost functions. The proposed method is validated with several numerical cases including randomly distributed monodisperse and polydisperse packings. Noticeable improvements in terms of calibration efficiency and accuracy have been made. 相似文献
2.
A multiscale DEM-LBM analysis on permeability evolutions inside a dilatant shear band 总被引：1，自引：1，他引：0
This paper presents a multiscale analysis of a dilatant shear band using a three-dimensional discrete element method and a lattice Boltzmann/finite element hybrid scheme. In particular, three-dimensional simple shear tests are conducted via the discrete element method. A spatial homogenization is performed to recover the macroscopic stress from the micro-mechanical force chains. The pore geometries of the shear band and host matrix are quantitatively evaluated through morphology analyses and lattice Boltzmann/finite element flow simulations. Results from the discrete element simulations imply that grain sliding and rotation occur predominately with the shear band. These granular motions lead to dilation of pore space inside the shear band and increases in local permeability. While considerable anisotropy in the contact fabric is observed with the shear band, anisotropy of the permeability is, at most, modest in the assemblies composed of spherical grains. 相似文献
3.
4.
基于离散单元法对不同密实度理想散粒体进行了双轴剪切试验的宏微观数值模拟,通过网格剖分将Voronoi多边形表征的loop单元作为散粒介质细观力学结构的基本单元,模拟了剪切过程中不同类型loop单元数量、几何形态和力学特征的演化过程,并重点分析了临界状态时散粒介质的细观力学结构特征。模拟结果显示,初始密实度不同的试样在向临界状态发展的过程中,高阶单元与低阶单元的发展规律完全不同,不同初始密实度试样中同阶loop单元的发展规律也不相同,但同阶loop单元的数量比例、几何形态、颗粒接触力及单元内滑动率最终均达到了各自的临界状态。从细观角度分析,散体介质的临界状态是高阶和低阶loop单元在荷载作用下相互转化的结果,是所有loop单元物理力学状态的综合平均与外在表现,临界状态时不同阶数的loop单元处于一个动态平衡状态,宏观上表现为常剪应力和常体积下剪切变形的不断发展。数值模拟结果也表明,loop细观结构单元包含了丰富的信息,其数量、几何形态、受力特征及接触稳定性的发展与散粒体的强度、剪胀以及临界状态的发展密切相关,可以将其作为散粒介质细观尺度的分析单元。 相似文献
5.
A three-dimensional discrete element model is used to investigate the effect of grain crushing on the tip resistance measured by cone penetration tests (CPT) in calibration chambers. To do that a discrete analogue of pumice sand, a very crushable microporous granular material, is created. The particles of the discrete model are endowed with size-dependent internal porosity and crushing resistance. A simplified Hertz–Mindlin elasto-frictional model is used for contact interaction. The model has 6 material parameters that are calibrated using one oedometer test and analogies with similar geomaterials. The calibration is validated reproducing other element tests. To fill a calibration chamber capable of containing a realistic sized CPT the discrete analogue is up-scaled by a factor of 25. CPT is then performed at two different densities and three different confinement pressures. Cone tip resistance in the crushable material is practically insensitive to initial density, as had been observed in previous physical experiments. The same CPT series is repeated but now particle crushing is disabled. The ratios of cone tip resistance between the two types of simulation are in good agreement with previous experimental comparisons of hard and crushable soils. Microscale exploration of the models indicates that crushing disrupts the buttressing effect of chamber walls on the cone. 相似文献
6.
Acta Geotechnica - Computers have been taught to clone granular soil particles for discrete element method simulations to alleviate difficulties of using three-dimensional imaging techniques for... 相似文献
7.
卵砾石地层被掘削时,其主要破坏模式并非一般的宏观整体剪切破坏,而是刀具切入卵砾石间的缝隙,将卵砾石颗粒剥离出原始地层,该过程较为契合颗粒离散元的物理力学假设。为确定掘削工况下的颗粒材料最优细观参数,对不同的颗粒接触法向刚度、接触刚度比、粒间摩擦系数、颗粒形状、颗粒尺寸等的颗粒材料开展了截割三维仿真试验,并分析研究了不同细观参数对掘削效果的影响。计算结果显示：颗粒接触法向刚度、接触刚度比的多种细观参数组合可以对应同一个宏观弹性模量值,接触刚度比大的细观参数组合的抗掘削能力更大。颗粒形状、粒间摩擦系数的多种细观参数组合可以对应同一个内摩擦角值,形状复杂的颗粒材料具有更大的抗掘削能力。在宏观参数相同的情况下,平均粒径更大的颗粒材料抗掘削能力更大。因此,利用颗粒离散元模拟掘削卵砾石地层问题时,除需保证虚拟颗粒材料的宏观整体剪切破坏响应与真实材料基本一致外,还需根据材料的抗掘削能力进一步标定细观参数,以获得更贴近实际情况的模拟效果。文章提出的方法可以为其他颗粒离散元掘削工况数值模拟中的参数标定提供参照依据。 相似文献
8.
A three-dimensional Voronoi tessellation model based on the distinct element method (DEM) is proposed to model the representative part of the microstructures of granular brittle rocks. Regularization is employed to decrease the frequency of polyhedrons with large edge ratio and contributes to a higher efficiency for element meshing. Sensitivity analyses are performed for a series of micro contact parameters in accordance with the macro responses observed in laboratory experiments (e.g. the uniaxial compression test, the Brazilian disc test and the triaxial compression test). Verifications by simulating the spalling test and plate impact test indicate that the 3D polycrystalline discrete element method (3PDEM) can be employed for efficiently simulating nonlinear mechanical behaviors, large deformation, strain softening and rock dynamics. 相似文献
9.
10.
A homogenization strategy for granular materials is presented and applied to a three-dimensional discrete element method (DEM), that uses superellipsoids as particles. Macroscopic quantities are derived from the microscopic quantities resulting from a DEM simulation by averaging over representative volume elements (RVEs). The implementation of an RVE is described in detail regarding the definition and discretization of the RVE boundary. The homogenization strategy is validated by DEM simulations of compression and shear tests of cohesionless granular assemblies. Finally, an elasto-plastic material model is fitted to the resulting stress–strain curves. 相似文献
11.
Stress-induced evolution of anisotropic thermal conductivity of dry granular materials 总被引：1，自引：1，他引：0
Jinhyun Choo Young Jin Kim Jung Hwoon Lee Tae Sup Yun Jangguen Lee Young Seok Kim 《Acta Geotechnica》2013,8(1):91-106
Various factors, such as the volumetric fraction of constituents, mineralogy, and pore fluids, affect heat flow in granular materials. Although the stress applied on granular materials controls the formation of major pathways for heat flow, few studies have focused on a detailed investigation of its significance with regard to the thermal conductivity and anisotropy of the materials. This paper presents a numerical investigation of the stress-induced evolution of anisotropic thermal conductivity of dry granular materials with supplementary experimental results. Granular materials under a variety of stress conditions in element testing are analyzed by the three-dimensional discrete element method, and quantitative variations in their anisotropic effective thermal conductivity are calculated via the network model and conductivity tensor measurements. Results show that the directional development of contact area and fabric under anisotropic stress conditions leads to the evolution of anisotropy in thermal conductivity. The anisotropy induced in thermal conductivity by shear stress is higher than that induced by compressive stress because shear stress causes more significant changes in microstructural configurations and boundary conditions. The shear-stress-induced evolution of anisotropy between principal thermal conductivities depends on dilatancy as well as shearing mode, and the shear-driven discontinuity localizes the conductivity. Factors involved in the stress-induced evolution and their implications on the thermal conductivity characterization are discussed. 相似文献
12.
This paper introduces a new generator algorithm and computer program for 3-D numerical simulation of packing configuration in a granular assemblies composed of ellipsoidal particles of different a/b aspect ratios. Each ellipsoidal particle is approximated by the revolution of an ellipse, formed by four connected arcs, about the major axis passing through its centroid. The centroid co-ordinates, major axis direction and lengths of the major and minor axes are the essential data for the packing generation and associated contact detection. The domain to be filled with particles can be a polyhedron of any shape. The packing program was coded based on a newly proposed scheme which obeys the no interpenetration kinematics of solid bodies. New contact detection algorithms for any two ellipsoids in the packing space were developed. Though simple, these algorithms effectively determine the contact condition and contact point without solving the simultaneous equations of the two ellipsoidal surfaces. Each particle's packing location, contact-point co-ordinates, and three-dimensional graphs can be created using the packing domain given boundaries, along with numbers, and geometrical information of particles to be generated. Simulation results show that this new algorithm provides an effective packing model as a required initial input for analysing the mechanics of granular material. This generation scheme potentially can explore the complex 3-D behaviours of material composed of discrete particles. Copyright © 1999 John Wiley & Sons, Ltd. 相似文献
13.
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. 相似文献
14.
The critical state is significant to the mechanical behaviors of granular materials and the foundation of the constitutive relations. Using the discrete element method (DEM), the mechanical behaviors of granular materials can be investigated on both the macroscopic and microscopic levels. A series of DEM simulations under true triaxial conditions have been performed to explore the critical state and dilatancy behavior of granular materials, which show the qualitatively similar macroscopic responses as the experimental results. The critical void ratio and stress ratio under different stress paths are presented. A unique critical state line (CSL) is shown to indicate that the intermediate principal stress ratio does not influence the CSL. Within the framework of the unique critical state, the stress–dilatancy relation of DEM simulations is found to fulfill the state-dependent dilatancy equations. As a microscopic parameter to evaluate the static determinacy of the granular system, the redundancy ratio is defined and investigated. The results show that the critical state is very close to the statically determinate state. Other particle-level indexes, including the distribution of the contact forces and the anisotropies, are carefully investigated to analyze the microstructural evolution and the underlying mechanism. The microscopic analysis shows that both the contact orientations and contact forces influence the mechanical behaviors of granular materials. 相似文献
15.
为了研究颗粒棱角对颗粒材料力学行为的影响,建立了具有不同棱角度的对称多面体颗粒,采用了一种简单并适合任意颗粒形状的接触本构模型,对三维离散元开源程序YADE进行了修改,研究了颗粒棱角度在模拟直剪试验中的影响以及接触力各向异性在剪切过程中的演化规律。研究结果表明,颗粒棱角度越小,颗粒间相互咬合自锁的作用越小,颗粒受剪更易转动,致使颗粒体系的剪切强度和剪胀性下降;竖向加载力越大,颗粒棱角度的影响越明显;法向接触力的各向异性在剪切过程中表现为先增后减最后趋向稳定的趋势;法向接触力的各向异性变化程度随颗粒棱角度的增大而增大。 相似文献
16.
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). 相似文献
17.
《Geomechanics and Geoengineering》2013,8(3):221-236
This paper presents a new point-to-face contact algorithm for contacts between two polyhedrons with planar boundaries. A new discrete numerical method called three-dimensional discontinuous deformation analysis (3-D DDA) is used and formulations of normal contact submatrices based on the proposed algorithm are derived. The presented algorithm is a simple and efficient method and it can be easily coded into a computer program. This approach does not need to use an iterative algorithm in each time step to obtain the contact plane, unlike the ‘Common-Plane’ method applied in the existing 3-D DDA. In the present 3-D DDA method, block contact constraints are enforced using the penalty method. This approach is quite simple, but may lead to inaccuracies that may be large for small values of the penalty number. The penalty method also creates block contact overlap, which violates the physical constraints of the problem. These limitations are overcome by using the augmented Lagrangian method that is used for normal contacts in this research. This point-to-face contact model has been programmed and some illustrative examples are provided to demonstrate the new contact rule between two blocks. A comparison between results obtained by using the augmented Lagrangian method and the penalty method is presented as well. 相似文献
18.
A succinct 3D discrete element model, with clumps to resemble the real shapes of granular materials, is developed. The quaternion method is introduced to transform the motion and force of a clump between local and global coordinates. The Hertz–Mindlin elastic contact force model, incorporated with the nonlinear normal viscous force and the Mohr–Coulomb friction law, is used to describe the interactions between particles. The proposed discrete element model is used to simulate direct shear tests of the irregular limestone rubbles. The simulation results of vertical displacements and shear stresses with a mixture of clumps are compared well with that of laboratory tests. The bulk friction coefficients are calculated and discussed under different contact friction coefficients and normal stresses. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
19.
非连续变形分析方法（DDA）是一种平行于有限元法的新型数值计算方法,该方法基于最小势能原理,把每个离散块体的变形、运动和块体之间的接触统一到平衡方程中进行隐式求解。然而,传统DDA方法在计算过程中需组装整体刚度矩阵并联立求解方程组,在用于大型岩土工程问题的三维数值模拟时占用内存较大、耗时较长、计算效率极低。因此,提出一种基于显式时间积分的三维球颗粒DDA方法。该方法在求解过程中不需要组装整体刚度矩阵,在求解加速度时,由于质量矩阵为对角矩阵,可存储为一维向量占用内存较少,且可分块逐自由度求解,效率较高,在接触判断上采用最大位移准则简化了接触算法,采用较小的时步,保证了计算的精确性;通过几个典型算例验证了该方法的准确性及计算效率。 相似文献
20.
Revisiting the existence of an effective stress for wet granular soils with micromechanics 
下载免费PDF全文
下载免费PDF全文 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. 相似文献