A three‐dimensional micromechanically based model

Granular materials react with complicated mechanical responses when subjected to external loading paths. This leads to sophisticated constitutive formulations requiring large numbers of parameters. A powerful and straightforward way consists in developing micro‐mechanical models embedding both micro‐scale and meso‐scale. This paper proposes a 3D micro‐mechanical model taking into account an intermediate scale (meso‐scale) that makes it possible to describe a variety of constitutive features in a natural way. The comparison between experimental tests and numerical simulations reveals the predictive capability of this model. Particularly, several simulations are carried out with different confining pressures and initial void ratios, based on the fact that the critical state is quantitatively described without requiring any critical state formulations and parameter. The model mechanism is also analyzed from a microscopic view, wherein the evolution of some key microscopic parameters is investigated. Copyright © 2017 John Wiley & Sons, Ltd.     

A micro‐mechanical study of peak strength and critical state

We present a micro‐mechanical analysis of macroscopic peak strength, critical state, and residual strength in two‐dimensional non‐cohesive granular media. Typical continuum constitutive quantities such as frictional strength and dilation angle are explicitly related to their corresponding grain‐scale counterparts (e.g., inter‐particle contact forces, fabric, particle displacements, and velocities), providing an across‐the‐scale basis for a better understanding and modeling of granular materials. These multi‐scale relations are derived in three steps. First, explicit relations between macroscopic stress and strain rate with the corresponding grain‐scale mechanics are established. Second, these relations are used in conjunction with the non‐associative Mohr–Coulomb criterion to explicitly connect internal friction and dilation angles to the micro‐mechanics. Third, the mentioned explicit connections are applied to investigate, understand, and derive micro‐mechanical conditions for peak strength, critical state, and residual strength. Copyright © 2015 John Wiley & Sons, Ltd.     

On incremental non‐linearity in granular media: phenomenological and multi‐scale views (Part I)

On the basis of fundamental constitutive laws such as elasticity, perfect plasticity, and pure viscosity, many elasto‐viscoplastic constitutive relations have been developed since the 1970s through phenomenological approaches. In addition, a few more recent micro‐mechanical models based on multi‐scale approaches are now able to describe the main macroscopic features of the mechanical behaviour of granular media. The purpose of this paper is to compare a phenomenological constitutive relation and a micro‐mechanical model with respect to a basic issue regularly raised about granular assemblies: the incrementally non‐linear character of their behaviour. It is shown that both phenomenological and micro‐mechanical models exhibit an incremental non‐linearity. In addition, the multi‐scale approach reveals that the macroscopic incremental non‐linearity could stem from the change in the regime of local contacts between particles (from plastic regime to elastic regime) in terms of the incremental macroscopic loading direction. Copyright © 2005 John Wiley & Sons, Ltd.     

基于细观试验和离散元法的盐岩力学特性

当前盐岩的宏观力学模型通常是唯象模型,不能很好地解释盐岩受力变形破坏的真正物理基础。盐岩是由于化学沉积而形成的矿物集合体,是一种主要由NaCl和少量杂质组成的多晶体,其变形机制主要由晶粒与晶界的力学特性控制。通过扫描电镜（SEM）,获得盐岩晶粒的微细观结构特征,采用分子动力学方法和纳米压痕技术,确定盐岩晶粒和晶界的微细观力学参数;将盐岩晶粒作为块体,基于Voronoi多边形技术,建立盐岩的微细观数值模型;利用离散元方法,对盐岩试件在单轴压缩和直剪条件下的宏观力学行为进行了数值模拟。数值模拟结果与宏观力学试验结果吻合度高,表明基于盐岩微细观晶粒结构特征并结合离散元数值模拟的方法能够较好地研究盐岩的宏观力学性能及其材料物理基础。     

Micromechanical analysis of the failure process of brittle rock

The failure process of brittle rock submitted to a compression state of stress with different confining pressures is investigated in this paper based on discrete element method (DEM) simulations. In the DEM model, the rock sample is represented by bonding rigid particles at their contact points. The numerical model is first calibrated by comparing the macroscopic response with the macroscopic response of Beishan granite obtained from laboratory tests. After the validation of numerical model in terms of macroscopic responses, the failure process of the DEM model under unconfined and confined compression is studied in micro‐scale in detail. The contact force network and its relation to the development of micro‐cracks and evolution of major fractures are studied. Confining pressure will prohibit the development of tensile cracks and hence alter the failure patterns. An in‐depth analysis of micro‐scale response is carried out, including the orientation distribution and probability density of stress acting on parallel bonds, the effect of particle size heterogeneity on bond breakage and the evolution of fabric tensor and coordination number of parallel bond. The proposed micromechanical analysis will allow us to extract innovative features emerged from the stresses and crack evolution in brittle rock failure process. Copyright © 2014 John Wiley & Sons, Ltd.     

Polygonal grain‐based distinct element modeling for mechanical behavior of brittle rock

The microstructure of rock was numerically reproduced by a polygonal grain‐based model, and its mechanical behavior was examined by performing the uniaxial compression test and Brazilian tests via the Universal Distinct Element Code. The numerical results of the model demonstrated good agreement with the experimental results obtained with rock specimens in terms of the stress–strain behavior, strength characteristics, and brittle fracture phenomenon. An encouraging result is that the grain‐based model‐Universal Distinct Element Code model can reproduce a low ratio of tensile to compressive strength of 1/20 to 1/10 without the need for an additional process. This finding is ascribed to the fact that the geometrical features of polygons can effectively capture the effects of angularity, finite rotation, and interlocking of grains that exist in reality. A numerical methodology to monitor the evolution of micro‐cracks was developed, which enabled us to examine the progressive process of the failure and distinguish the contribution of tensile cracking to the process from that of shear cracking. From the observations of the micro‐cracking process in reference to the stress–strain relation, crack initiation stress, and crack damage stress, it can be concluded that the failure process of the model closely resembles the microscopic observations of rock. We also carried out a parametric study to examine the relationships between the microscopic properties and the macroscopic behavior of the model. Depending on the micro‐properties, the model exhibited a variety of responses to the external load in terms of the strength and deformation characteristics, the evolution of micro‐cracks, and the post‐peak behavior. Copyright © 2016 John Wiley & Sons, Ltd.     

晶体及矿物颗粒大小对岩土材料力学性质的影响

岩石常见较大的晶体或者矿物颗粒,混凝土中是骨料,通过团簇模拟大颗粒的力学行为、团簇可以破裂。根据设计的相同数量、相同位置、不同半径的大颗粒数值单轴压缩试验,在虚拟试验条件下,考察了颗粒大小对材料力学响应的影响。通过分析颗粒材料的破裂形态、裂纹扩展过程、应力-应变曲线和破裂能量演化规律发现：大颗粒具有明显的增强特性,有阻止裂纹扩展的作用,破裂多绕大颗粒发展;增强幅度随颗粒半径的增加呈单增趋势,半径较小时,增强效果不明显。     

3‐D mesoscale simulation of crack‐permeability coupling in the Brazilian splitting test

In this paper, mesoscale hydromechanical simulations are performed to study (1) fracture features and (2) crack‐gas permeability coupling evolution in the context of the tensile splitting test. The mesostructure is based on a 2‐phase 3‐D representation of heterogeneous materials, such as concrete, where stiff aggregates are embedded into a mortar matrix. To take into account these heterogeneities without any mesh adaptation, a weak discontinuity is introduced into the strain field. In addition, a strong discontinuity is also added to take into account microcracking. This mechanical model is cast into the framework of the enhanced finite element method. Concerning the coupling with gas permeability, a double‐porosity method is used to simulate the flow through the cracks and the porosity. The apparent gas permeability is afterwards evaluated by a homogenization method. On the basis of finite element simulations, influence of aggregate size on ultimate crack opening, macroscopic ultimate tensile stress, total dissipated energy, and gas permeability evolution is numerically investigated. Furthermore, gas permeability evolution is also compared with experimental results from the literature. In addition, in the spirit of a sequential multiscale approach, macroscale gas permeability equations are identified from the hydromechanical results coming from the mesoscale computations. These equations lead to a relation between macroscale gas permeability evolution and crack opening. Besides, we show how the aggregate size influences the percolation threshold and that after this threshold, a cubic relation between macroscale gas permeability and crack opening is obtained.     

Analysis of soil‐structural interface behavior using three‐dimensional DEM simulations

The shear behavior at the interface between the soil and a structure is investigated at the macroscale and particle‐scale levels using a 3‐dimensional discrete element method (DEM). The macroscopic mechanical properties and microscopic quantities affected by the normalized interface roughness and the loading parameters are analyzed. The macro‐response shows that the shear strength of the interface increases as the normalized roughness of the interface increases, and stress softening and dilatancy of the soil material are observed in the tests that feature rough interfaces. The particle‐scale analysis illustrates that a localized band characterized by intense shear deformation emerges from the contact plane and gradually expands as shearing progresses before stabilizing at the residual stress state. The thickness of the localized band is affected by the normalized roughness of the interface and the normal stress, which ranges between 4 and 5 times that of the median grain diameter. A thicker localized band is formed when the soil has a rough shearing interface. After the localized band appears, the granular material structuralizes into 2 regions: the interface zone and the upper zone. The mechanical behavior in the interface zone is representative of the interface according to the local average stress analysis. Certain microscopic quantities in the interface zone are analyzed, including the coordination number and the material fabric. Shear at the interface creates an anisotropic material fabric and leads to the rotation of the major principal stress.     

A generalized rigid particle contact model for fracture analysis

Rigid particle models taking directly into consideration the physical mechanisms and the influence of the material meso‐structure have recently been developed for fracture studies of quasi‐brittle material such as concrete. The formulation of a generalized contact model for rigid particle simulations is presented in which the contact discretization is a model parameter. The contact model performance for different discretizations is evaluated for uniaxial tensile tests, for uniaxial compression tests and for a notched beam in mode I. Copyright © 2005 John Wiley & Sons, Ltd.     

A macroscopic damage‐plastic constitutive law for modeling quasi‐brittle fracture and ductile behavior of concrete

A new phenomenological macroscopic constitutive model for the numerical simulation of quasi‐brittle fracture and ductile concrete behavior, under general triaxial stress conditions, is presented. The model is particularly addressed to simulate a wide range of confinement stress states, as also, to capture the strong influence of the mean stress value in the concrete failure mechanisms. The model is based on a two‐surface damage‐plastic formulation. The mechanical behavior in different domains of the stress space is separately described by means of a quasi‐brittle or ductile material response:

相似文献   

Studying the effect of non‐spherical micro‐particles on Hoek–Brown strength parameter mi using numerical true triaxial compressive tests

The strength parameter m_i in the Hoek–Brown strength criterion is empirical and was developed by trial and error. To better understand the fundamental relationship between m_i and the physical characteristics of intact rock, this paper presents a systematic study of m_i by representing intact rock as a densely packed cemented particle material and simulating its mechanical behavior using particle flow modeling. Specifically, the three‐dimensional particle flow code (PFC3D) was used to conduct numerical true triaxial compression tests on intact rock and to investigate the effect of non‐spherical micro‐particle parameters on m_i. To generate numerical intact rock specimens containing non‐spherical micro‐particles, a new genesis process was proposed, and a specific loop algorithm was used based on the efficiency of the process and the acceptability of generated specimens. Four main parameters—number, aspect ratio, size, and shape—of non‐spherical micro‐particles were studied, and the results indicated that they all have great effect on m_i. The strength parameter m_i increases when the number, aspect ratio, or size is larger or the shape becomes more irregular, mainly as a result of the higher level of interlocking between particles. This confirms the observations from engineering experience and laboratory experiments. To simulate the right strength parameter m_i, it is important to use appropriate non‐spherical micro‐particles by controlling these four parameters. This is further demonstrated by the simulation of two widely studied rocks, Lac du Bonnet granite and Carrara marble. Copyright © 2014 John Wiley & Sons, Ltd.     

Numerical determination of the tensile response and the dissipated fracture energy of concrete: role of the mesostructure and influence of the loading rate

At the mesoscopic scale, concrete can be considered as a mix of coarse aggregates with a mortar paste matrix. In this paper, we investigate numerically the influence of aggregates arrangements and loading rate on the tensile response of concrete. Each coarse aggregate is assumed to be circular with six different radiuses following the aggregates size distribution of real gravel. Rate‐independent cohesive elements are used to model failure within the mesostructure. Our results show that the spatial distribution of heterogeneities does not influence the peak strength, while it changes the post‐peak macroscopic response. This implies that our specimen size is large enough for strength computation but that larger mesostructures should be considered to obtain fully reliable toughness predictions. Although the cohesive approach is able to capture the transition from one macro‐crack in quasi‐static to multiple micro‐cracks in fast dynamics, which increases the dissipated fracture energy, our results suggest that the full extent of the high‐rate strengthening of concrete observed experimentally for loading rates greater than 1/s cannot be captured with rate independent constitutive laws.Copyright © 2013 John Wiley & Sons, Ltd.     

Chemo‐mechanics of cemented granular solids subjected to precipitation and dissolution of mineral species

This paper studies the chemo‐mechanics of cemented granular solids in the context of continuum thermodynamics for fluid‐saturated porous media. For this purpose, an existing constitutive model formulated in the frame of the Breakage Mechanics theory is augmented to cope with reactive processes. Chemical state variables accounting for the reactions between the solid constituents and the solutes in the pore fluid are introduced to enrich the interactions among the microstructural units simulated by the model (i.e., grains and cement bonds). Two different reactive processes are studied (i.e., grain dissolution and cement precipitation), using the chemical variables to describe the progression of the reactions and track changes in the size of grains and bonds. Finally, a homogenization strategy is used to derive the energy potentials of the solid mixture, adopting probability density functions that depend on both mechanical and chemical indices. It is shown that the connection between the statistics of the micro‐scale attributes and the continuum properties of the solid enables the mathematical capture of numerous mechanical effects of lithification and chemical deterioration, such as changes in stiffness, expansion/contraction of the elastic domain, and development of inelastic strains during reaction. In particular, the model offers an interpretation of the plastic strains generated by aggressive environments, which are here interpreted as an outcome of chemically driven debonding and comminution. As a result, the model explains widely observed macroscopic signatures of geomaterial degradation by reconciling the energetics of the deformation/reaction processes with the evolving geometry of the microstructural attributes. Copyright © 2015 John Wiley & Sons, Ltd.     

Comparison of continuous and discontinuous constitutive models to simulate concrete behaviour under mixed‐mode failure conditions

The paper presents detailed FE simulation results of concrete elements under mixed‐mode failure conditions according to the so‐called shear‐tension test by Nooru‐Mohamed, characterized by curved cracks. A continuous and discontinuous numerical two‐dimensional approach was used. In order to describe the concrete's behaviour within continuum mechanics, two different constitutive models were used. First, an elasto‐plastic model with isotropic hardening and softening was assumed. In a compression regime, a Drucker–Prager criterion with a non‐associated flow rule was used. In turn, in a tensile regime, a Rankine criterion with an associated flow rule was adopted. Second, an isotropic damage constitutive model was applied with a single scalar damage parameter and different definitions of the equivalent strain. Both constitutive laws were enriched by a characteristic length of micro‐structure to capture properly strain localization. As an alternative approach, the extended finite element method was used. Our results were compared with the experimental ones and with results of other FE simulations reported in the literature. Copyright © 2015 John Wiley & Sons, Ltd.     

Computational mechanics of the steel–concrete interface

Concrete cracking in reinforced concrete structures is governed by two mechanisms: the activation of bond forces at the steel–concrete interface and the bridge effects of the reinforcement crossing a macro‐crack. The computational modelling of these two mechanisms, acting at different scales, is the main objective of this paper. The starting point is the analysis of the micro‐mechanisms, leading to an appropriate choice of (measurable) state variables describing the energy state in the surface systems: on the one side the relative displacement between the steel and the concrete, modelling the bond activation; on the other hand, the crack opening governing the bridge effects. These displacement jumps are implemented in the constitutive model using thermodynamics of surfaces of discontinuity. On the computational side, the constitutive model is implemented in a discrete crack approach. A truss element with slip degrees of freedom is developed. This degree of freedom represents the relative displacement due to bond activation. In turn, the bridge effect is numerically taken into account by modifying the post‐cracking behaviour of the contact elements representing discrete concrete cracks crossed by a rebar. First simulation results obtained with this model show a good agreement in crack pattern and steel stress distribution with micro‐mechanical results and experimental results. Copyright © 2001 John Wiley & Sons, Ltd.     

Parametric study of the smooth‐joint contact model on the mechanical behavior of jointed rock

The smooth‐joint contact model based on distinct element method has been widely used to represent discontinuity in the simulation of fractured rock mass, but there is rare efficient guidance for the selection of proper parameters of smooth‐joint contact model, which is the basement for using this model properly. In this paper, the effect of smooth joint parameters on the macroscopic properties and failure mechanism of jointed rock under triaxial compression test is investigated. The numerical results reveal that the friction coefficient of smooth joint plays a dominant role in controlling mechanical behaviors. The stiffness of smooth joint has a relative small influence on the mechanical behaviors. Poisson ratio decreases with the reduction of normal stiffness but increases with the reduction of shear stiffness. The reduction of smooth joint strength, which is determined by normal strength, cohesion, and friction angle of smooth joint, contributes to the breakage of bonded smooth joint and ultimately decreases the strength of the specimen. We proposed a detailed calibration process for smooth‐joint contact model according to the relationship between smooth‐joint parameters and mechanical properties. By following this process, the numerical results are validated against corresponding experimental results and good agreement between them can be found in stress‐strain curves and failure modes of different joint orientations. Further analyses from the microperspective are performed by looking at transmission of contact force, the nature and distribution of microcracks, and the particle displacement to show the failure process and failure modes.     

Coupled modelling of hydro‐mechanical behaviour of unsaturated compacted expansive soils

The Barcelona basic model cannot predict the mechanical behaviour of unsaturated expansive soils, whereas the Barcelona expansive model (BExM) can only predict the stress–strain behaviour of unsaturated expansive soils without the water‐retention behaviour being incorporated. Moreover, the micro‐parameters and the coupling function between micro‐structural and macro‐structural strains in the BExM are difficult to determine. Experimental data show that the compression curves for non‐expansive soils under constant suctions are shifted towards higher void ratios with increasing suction, whereas the opposite is true for expansive soils. According to the observed water‐retention behaviour of unsaturated expansive soils, the air‐entry value increases with density, and the relationship between the degree of saturation and void ratio is linear at constant suction. According to the above observation, an elastoplastic constitutive model is developed for predicting the hydraulic and mechanical behaviour of unsaturated expansive soils, based on the existing hydro‐mechanical model for non‐expansive unsaturated soil. The model takes into consideration the effect of degree of saturation on the mechanical behaviour and that of void ratio on the water‐retention behaviour. The concept of equivalent void ratio curve is introduced to distinguish the plastic potential curve from the yield curve. The model predictions are compared with the test results of an unsaturated expansive soil, including swelling tests under constant net stress, isotropic compression tests and triaxial shear tests under constant suction. The comparison indicates that the model offers great potential for quantitatively predicting the hydraulic and mechanical behaviour of unsaturated expansive soils. Copyright © 2011 John Wiley & Sons, Ltd.     

A continuum‐discrete hydromechanical analysis of granular deposit liquefaction

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.