共查询到20条相似文献，搜索用时 500 毫秒
1.
强黏结土体或胶凝砂砾石料等颗粒材料的颗粒间存在较大黏结力,在低应力情况下主要由颗粒接触特性决定颗粒材料的宏观蠕变特性。基于颗粒细观力学方法,研究了二维各向同性颗粒材料的蠕变特性。首先引入速率相关的力与位移关系描述颗粒材料细观颗粒间特性。其次,运用Laplace变换,将时间域内的线黏弹性颗粒材料细观均匀化问题转化为拉氏空间内线弹性颗粒材料细观均匀化问题,随后基于颗粒材料线弹性问题在Reuss、Voigt和一般位移场三种假设下的解,通过Laplace逆变换得到颗粒材料相应的宏观蠕变特性解析模型,并建立了材料蠕变特性的上下限。最后,通过理论模型的解析解与商业软件（PFC2D）的数值结果的对比,验证了理论模型的合理性。 相似文献
2.
We formulate a discrete Lagrangian model for a set of interacting grains, which is purely elastic. The considered degrees of freedom for each grain include placement of barycenter and rotation. Further, we limit the study to the case of planar systems. A representative grain radius is introduced to express the deformation energy to be associated to relative displacements and rotations of interacting grains. We distinguish inter-grains elongation/compression energy from inter-grains shear and rotations energies, and we consider an exact finite kinematics in which grain rotations are independent of grain displacements. The equilibrium configurations of the grain assembly are calculated by minimization of deformation energy for selected imposed displacements and rotations at the boundaries. Behaviours of grain assemblies arranged in regular patterns, without and with defects, and similar mechanical properties are simulated. The values of shear, rotation, and compression elastic moduli are varied to investigate the shapes and thicknesses of the layers where deformation energy, relative displacement, and rotations are concentrated. It is found that these concentration bands are close to the boundaries and in correspondence of grain voids. The obtained results question the possibility of introducing a first gradient continuum models for granular media and justify the development of both numerical and theoretical methods for including frictional, plasticity, and damage phenomena in the proposed model. 相似文献
3.
This paper presents a computational homogenization approach in the framework of poromechanics. A fully coupled hydromechanics problem is formulated at the macroscopic scale. The constitutive equations are replaced by results of numerical computations on a Representative Elementary Volume in order to obtain the overall stress of the mixture as well as its transmissivity properties. At the microscale, the material is assumed to be composed of an assembly of hyperelastic grains connected by cohesive interfaces. These interfaces are also channels of a network where the fluid can percolate. The fluid acts on the boundaries of the grains and influences the behavior of the cohesive interfaces. Conversely, the opening of the interfaces induces changes in the transmissivity properties of the corresponding channels. This yields a fully coupled hydromechanical problem at the microscopic scale. The finite element method is considered for the numerical solutions at both scales, the present approach extending the purely mechanical FE2 scheme to the coupled hydromechanical framework. The local macroscopic behavior resulting from the homogenization scheme is illustrated on different numerical tests. The results clearly show the coupling between damage and fluid permeability in the overall response, as a consequence of the small-scale interaction between the action of the percolating fluid, the deformation of the solid skeleton, and the failure of the cohesive interfaces. 相似文献
4.
Jean‐Yves Delenne Moulay Saïd El Youssoufi Fabien Cherblanc Jean‐Claude Bnet 《国际地质力学数值与分析法杂志》2004,28(15):1577-1594
We investigate the stress–strain behaviour and failure of a cohesive granular material both by experiments and numerical simulations. The material is an assembly of aluminium rods glued together by means of an epoxy resin. The behaviour of cohesive bonds (force–displacement relationship, failure conditions) is characterized by performing simple loading tests (tension/compression, shear…) on a couple of rods. Then, this local behaviour is introduced in a numerical code based on a discrete element method in order to perform numerical compression tests on large samples. The validation of this approach was the main goal of the present investigation that is essentially achieved by a direct comparison between the numerical results and similar experimental tests. As a basic application, we derive the macroscopic cohesion and friction characteristics of random cohesive materials by systematic numerical simulations in a biaxial geometry. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
5.
Biaxial experiments are performed with rod assemblies to study the micro-mechanical deformation behavior of granular materials. The focus of these experiments is upon the micromechanical behavior under mixed boundary conditions, with stress-controlled lateral boundaries and displacement-controlled axial boundaries. Particle motions, i.e. displacements and rotations, are measured during the test. The particle motions are analyzed to study deformation patterns, including shear localization. Deformation patterns in a number of rod assemblies with random packing structures are investigated. 相似文献
6.
Stress–strain modeling of sand–silt mixtures is important in the analysis and design of earth structures. In this paper, we develop a stress–strain model that can predict the behavior of sand–silt mixtures with any amount of fines content. This model is based on a micromechanics approach, which involves mean‐field assumptions. For the mixtures with low amount of fines, the mechanical behavior is dominated by sand grains network. On the other hand, for the mixtures with high amount of fines, the mechanical behavior is dominated by silt grains network. Using this concept of dominant grains network, the behavior of mixtures with any amount of fines can be predicted from knowing the behavior of sand and silt, alone. We also modeled the critical state friction angle, critical state void ratio, and elastic stiffness for the mixtures as a function of fines content. The applicability of this developed stress–strain model is shown by comparing the simulated and measured results for two different types of sand–silt mixtures with full range of fines content. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
7.
Dadda Abdelali Geindreau Christian Emeriault Fabrice Rolland du Roscoat Sabine Esnault Filet Annette Garandet Aurélie 《Acta Geotechnica》2019,14(3):597-613
The mechanical efficiency of the biocementation process is directly related to the microstructural properties of the biocemented sand, such as the volume fraction of calcite, its distribution within the pore space (localized at the contact between grains, over the grain surfaces) and the contact properties: coordination number, contact surface area, contacts orientation and types of contact. In the present work, all these micromechanical properties are computed, for the first time, from 3D images obtained by X-ray tomography of intact biocemented sand samples. The evolution of all these properties with respect to the volume fraction of calcite is analyzed and compared between each other (from untreated sand to highly cemented sand). Whatever the volume fraction of calcite, it is shown that the precipitation of the calcite is localized at the contacts between grains. These results are confirmed by comparing our numerical results with analytical estimates assuming that the granular medium is made of periodic simple cubic arrangements of grains and by considering two extreme cases of precipitation: (1) The calcite is localized at the contact, and (2) the grains are covered by a uniform layer of calcite. In overall, the obtained results show that a small percentage of calcite is sufficient to get a large amount of cohesive contacts.
相似文献8.
An energy approach is proposed as a complement to the stress approach commonly considered for investigating soil desiccation cracking. The elastic strain energies before and after crack initiation are estimated by both numerical and analytical solutions. The energy released by cracking is then compared with the fracture energy to discuss crack initiation conditions. This leads to combined energy and stress conditions for crack initiation following Leguillon's theory. An approximate analytical solution is derived from a variational formulation of the porous elastic body equations. A cohesive zone model and finite element code are used to simulate crack propagation in an unsaturated porous body. This analysis shows that the energy criterion is reached before the stress criterion, and this can explain unstable crack propagation at the beginning. The approximate analytical solution allows predicting correctly the crack depth and opening in its initiation stage. 相似文献
9.
This paper is devoted to numerical analysis of strength and deformation of cohesive granular materials. The emphasis is put on the study of effects of confining pressure and loading path. To this end, the three-dimensional discrete element method is used. A nonlinear failure criterion for inter-granular interface bonding is proposed, and it is able to account for both tensile and shear failure for a large range of normal stress. This criterion is implemented in the particles flow code. The proposed failure model is calibrated from triaxial compression tests performed on representative sandstone. Numerical results are in good agreement with experimental data. In particular, the effect of confining pressure on compressive strength and failure pattern is well described by the proposed model. Furthermore, numerical predictions are studied, respectively, for compression and extension tests with a constant mean stress. It is shown that the failure strength and deformation process are clearly affected by loading path. Finally, a series of numerical simulations are performed on cubic samples with three independent principal stresses. It is found that the strength and failure mode are strongly influenced by the intermediate principal stress.
相似文献10.
In this paper, steady-state conditions for ideal monodisperse dry granular materials are both theoretically and numerically analysed. A series of discrete element (DEM) numerical simulations have been performed on a periodic cell by imposing stress paths characterized by different Lode angles, pressures, and deviatoric strain rates. The dependence of the material response on both inertial number and loading path has been discussed in terms of void ratio, fabric, and granular temperature. DEM numerical results have been finally compared with the prediction of an already conceived model based on both kinetic and critical state theories, here suitably modified to account for three-dimensional conditions. 相似文献
11.
12.
The paper develops the stiffness relationship between the movements and forces among a system of discrete interacting grains. The approach is similar to that used in structural analysis, but the stiffness matrix of granular material is inherently nonsymmetric because of the geometrics of particle interactions and of the frictional behavior of the contacts. Internal geometric constraints are imposed by the particles' shapes, in particular, by the surface curvatures of the particles at their points of contact. Moreover, the stiffness relationship is incrementally nonlinear, and even small assemblies require the analysis of multiple stiffness branches, with each branch region being a pointed convex cone in displacement space. These aspects of the particle-level stiffness relationship give rise to three types of microscale failure: neutral equilibrium, bifurcation and path instability, and instability of equilibrium. These three pathologies are defined in the context of four types of displacement constraints, which can be readily analyzed with certain generalized inverses. That is, instability and nonuniqueness are investigated in the presence of kinematic constraints. Bifurcation paths can be either stable or unstable, as determined with the Hill–Bažant–Petryk criterion. Examples of simple granular systems of three, 16, and 64 disks are analyzed. With each system, multiple contacts were assumed to be at the friction limit. Even with these small systems, microscale failure is expressed in many different forms, with some systems having hundreds of microscale failure modes. The examples suggest that microscale failure is pervasive within granular materials, with particle arrangements being in a nearly continual state of instability. 相似文献
13.
Simulation of large deformation and post‐failure of geomaterial in the framework of smoothed particle hydrodynamics (SPH) are presented in this study. The Drucker–Prager model with associated and non‐associated plastic flow rules is implemented into the SPH code to describe elastic–plastic soil behavior. In contrast to previous work on SPH for solids, where the hydrostatic pressure is often estimated from density by an equation of state, this study proposes to calculate the hydrostatic pressure of soil directly from constitutive models. Results obtained in this paper show that the original SPH method, which has been successfully applied to a vast range of problems, is unable to directly solve elastic–plastic flows of soil because of the so‐called SPH tensile instability. This numerical instability may result in unrealistic fracture and particles clustering in SPH simulation. For non‐cohesive soil, the instability is not serious and can be completely removed by using a tension cracking treatment from soil constitutive model and thereby give realistic soil behavior. However, the serious tensile instability that is found in SPH application for cohesive soil requires a special treatment to overcome this problem. In this paper, an artificial stress method is applied to remove the SPH numerical instability in cohesive soil. A number of numerical tests are carried out to check the capability of SPH in the current application. Numerical results are then compared with experimental and finite element method solutions. The good agreement obtained from these comparisons suggests that SPH can be extended to general geotechnical problems. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
14.
The paper presents some numerical results of the problem on contact of stratified elastic half‐spaces possessing unevenness of their surface layers. The approach of the Jacobi polynomials is used to reduce the singular integral equation of the problem to a system of linear algebraic equations. On the basis of the system's solution the dependence of gaps' geometric characteristics on external load and friction is investigated for two types of initial boundary disturbances, namely a local symmetric recess and a periodic set of recesses at a flat boundary. Numerical results are presented in figures for a few elastic and geometric parameters of joint components. The analysis shows that elastic properties and geometric structure of a laminated half‐space have a significant influence on the shape of interface gaps and their location. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
15.
N. Belheine J.-P. Plassiard F.-V. Donzé F. Darve A. Seridi 《Computers and Geotechnics》2009,36(1-2):320-331
A discrete element modeling of granular material was carried out using a 3D spherical discrete model with a rolling resistance, in order to take into account the roughness of grains. The numerical model of Labenne sand was generated, and the desired porosity was obtained by a radius expansion method. Using numerical triaxial tests the micro-mechanical properties of the numerical material were calibrated in order to match the macroscopic response of the real material. Numerical simulations were carried out under the same conditions as the physical experiments (porosity, boundary conditions and loading). The pre-peak, peak and post-peak behavior of the numerical material was studied. The calibration procedure revealed that the peak stress of the sand sample does not only depend on local friction parameters but also on the rolling resistance. The larger the value of the applied rolling resistance, the higher the resulting stress peak. Furthermore, the deformational response depends strongly on local friction. The numerical results are quantitatively in agreement with the laboratory test results. 相似文献
16.
Coupled theory of mixtures for clayey soils 总被引：4，自引：0，他引：4
In this work, elasto-plastic coupled equations are formulated in order to describe the time-dependent deformation of saturated cohesive soils (two-phase state). Formulation of these equations is based on the principle of virtual work and the theory of mixtures for inelastic porous media. The theory of mixtures for a linear elastic porous skeleton was first developed by Biot (Theory of elasticity and consolidation for a porous anisotropic solid, Journal of Applied Physics, 1955, 26, 188–185). An extension of Biot's theory into a nonlinear inelastic media was performed by Prevost (Mechanics of continuous porous media, International Journal of Engineering Science, 1980, 18, 787–800). The saturated soil is considered as a mixture of two deformable media, the solid grains and the water. Each medium is regarded as a continuum and follows its own motion. The flow of pore-water through the voids is assumed to follow Darcy's law. The coupled equations are developed for large deformations with finite strains in an updated Lagrangian reference frame. The coupled behavior of the two-phase materials (soil-water state) is implemented in a finite element program. A modified Cam-clay model is adopted and implemented in the finite element program in order to describe the plastic behavior of clayey soils. Penetration of a piezocone penetrometer in soil is numerically simulated and implemented into a finite element program. The piezocone penetrometer is assumed to be infinitely stiff. The continuous penetration of the cone is simulated by applying an incremental vertical movement of the cone tip boundary. Results of the finite element numerical simulation are compared with experimental measurements conducted at Louisiana State University using the calibration chamber. The numerical simulation is carried out for two cases. In the first case, the interface friction between the soil and the piezocone penetrometer is neglected. In the second case, interface friction is assumed between the soil and the piezocone. The results of the numerical simulations are compared with experimental laboratory measurements. 相似文献
17.
利用传递-反射矩阵方法和叠加原理推导了双轮胎振动荷载作用下弹性层状公路结构的动力响应解,并利用Hankel数值逆变换技术进行了数值求解。弹性层状结构的传递-反射矩阵同时具有单个矩阵尺寸小、矩阵元素只含有负指数项2个优点,计算结果不受土层单元厚度与荷载频率大小的影响,可以大幅提高计算效率和精度。通过将数值计算与室内模型试验结果进行比较,验证了方法的正确性。在此基础上,分析了面层厚度及模量对层状公路结构动应力、动位移的影响程度和规律,得出了一些有意义的结论。研究结果可为公路路面、路基结构设计提供参考。 相似文献
18.
Geometric analysis of nested Riedel structures was used to identify and quantify strain localization processes within faulted Navajo sandstone. The analysis shows systematic deviation from the basic Riedel geometry complying with the Mohr–Coulomb criterion. Using cross-cutting relations amongst deformation bands within the Riedel structures, and comparing the orientations of the deformation bands to theoretical strain calculations, we identify two coupling deformation mechanisms involved in the early stages of shear-zone evolution, namely, granular flow and discrete faulting. Both mechanisms localize during strain accumulation, and the granular flow facilitates considerable change in the initial geometry of the Riedel structures. The analysis demonstrates a systematic sequence, by which new Riedel structures form after a constant amount of shear strain takes place in the sandstone. Analysis further indicates that granular flow is the major deformation mechanism during early stages of shear-zone evolution and discrete faulting is the dominant mechanism during later deformation stages. 相似文献
19.
Recently, the shear behavior of a cohesionless granular strip that is in contact with a very rough surface of a moving bounding structure has been numerically investigated by several authors by using a micropolar hypoplastic continuum model. It was shown that the micropolar boundary conditions assumed along the interface have a strong influence on the deformations within the granular layer. In previous investigations, only interface friction angles for very rough bounding structures were assumed. In contrast, the focus of the present paper is on the influence of the interface roughness on the deformation behavior of the granular strip when the interface friction angle is lower than the peak friction angle of the granular material. In addition to the interface friction angle, particular attention is also paid to the influence of the mean grain diameter, the solid hardness, the initial void ratio, and the vertical stress on the maximum horizontal shear displacement within the granular layer before sliding is started. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
20.
Whilst faulting in the shallow crust is inevitably associated with comminution of rocks, the mechanical properties of the comminuted granular materials themselves affect the slip behavior of faults. Therefore, the mechanical behavior of any fault progresses along an evolutionary path. We analyzed granular fault rocks from four faults, and deduced an evolutionary trend of fractal size frequency. Comminution of fault rocks starts at a fractal dimension close to 1.5 (2-D measurement), at which a given grain is supported by the maximum number of grains attainable and hence is at its strongest. As comminution proceeds, the fractal dimension increases, and hence comminution itself is a slip weakening mechanism. Under the appropriate conditions, comminuted granular materials may be fluidized during seismic slip events. In this paper, we develop a new method to identify the granular fault rocks that have experienced fluidization, where the detection probability of fragmented counterparts is a key parameter. This method was applied to four fault rock samples and a successful result was obtained. Knowledge from powder technology teaches us that the volume fraction of grains normalized by maximum volume fraction attainable is the most important parameter for dynamic properties of granular materials, and once granular fault materials are fluidized, the fault plane becomes nearly frictionless. A small decrease in the normalized volume fraction of grains from 1 is a necessary condition for the phase transition to fluidization from the deformation mechanism governed by grain friction and crushing by contact stresses. This condition can be realized only when shearing proceeds under unconstrained conditions, and this demands that the gap between fault walls is widened. Normal interface vibration proposed by Brune et al. [Tectonophysics 218 (1993) 59] appears to be the most appropriate cause of this, and we presented two lines of field evidence that support this mechanism to work in nature. 相似文献