Numerical evaluation of three non‐coaxial kinematic models using the distinct element method for elliptical granular materials

This paper presents a numerical evaluation of three non‐coaxial kinematic models by performing Distinct Element Method (DEM) simple shear tests on specimens composed of elliptical particles with different aspect ratios of 1.4 and 1.7. The models evaluated are the double‐shearing model, the double‐sliding free‐rotating model and the double slip and rotation rate model (DSR² model). Two modes of monotonic and cyclic simple shear tests were simulated to evaluate the role played by the inherent anisotropy of the specimens. The main findings are supported by all the DEM simple shear tests, irrespective of particle shape, specimen density or shear mode. The evaluation demonstrates that the assumption in the double‐shearing model is inconsistent with the DEM results and that the energy dissipation requirements in the double‐sliding free‐rotating model appear to be too restrictive to describe the kinematic flow of elliptical particle systems. In contrast, the predictions made by the DSR² model agree reasonably well with the DEM data, which demonstrates that the DSR² model can effectively predict the non‐coaxial kinematic behavior of elliptical particle systems. Copyright © 2016 John Wiley & Sons, Ltd.     

Investigation of influence of particle characteristics on the non‐coaxiality of anisotropic granular materials using DEM

As a result of deposition process and particle characteristics, granular materials can be inherently anisotropic. Many researchers have strongly suggested that the inherent anisotropy is the main reason for the deformation non‐coaxiality of granular materials. However, their relationships are not unanimous because of the limited understanding of the non‐coaxial micro‐mechanism. In this study, we investigated the influence of inherent anisotropy on the non‐coaxial angle using the discrete element method. Firstly, we developed a new discrete element method approach using rough elliptic particles and proposed a novel method to produce anisotropic specimens. Secondly, the effects of initial specimen density and particle characteristics, such as particle aspect ratio A _m, rolling resistance coefficient β , and bedding plane orientation δ , were examined by a series of biaxial tests and rotational principal axes tests. Findings from the numerical simulations are summarized as follows: (1) the peak internal friction angle ? _p and the non‐coaxial angle i both increase with the initial density, A _m and β , and they both increase initially and then decrease with δ in the range of 0–90°; (2) among the particle characteristics, the influence of A _m is the most significant; and (3) for anisotropic specimens, the non‐coaxial angle can be calculated using the double slip and rotation rate model. Then, an empirical formula was proposed based on the simulation results to depict the relationship between the non‐coaxial angle and the particle characteristics. Finally, the particle‐scale mechanism of non‐coaxiality for granular materials was discussed from the perspective of energy dissipation. Copyright © 2016 John Wiley & Sons, Ltd.     

Numerical investigations of shear localization in a micro‐polar hypoplastic material

In this paper a micro‐polar continuum approach is proposed to model the essential properties of cohesionless granular materials like sand. The model takes into account the influence of particle rotations, the mean grain size, the void ratio, the stresses and couple stresses. The constitutive equations for the stresses and couple stresses are incrementally non‐linear and based on the concept of hypoplasticity. For plane strain problems the implementation of the model in a finite element program is described. Numerical studies of the evolution of micro‐polar effects within a granular strip under plane shearing are presented. It is shown that the location and evolution of shear localization is strongly influenced by the initial state and the micro‐polar boundary conditions. For large shearing the state quantities tend towards a stationary state for which a certain coupling between the norm of the stress deviator and the norm of the couple stress tensor can be derived. Copyright © 2003 John Wiley & Sons, Ltd.     

Two‐scale modeling of solute dispersion in unsaturated double‐porosity media: Homogenization and experimental validation

A two‐scale modeling of solute transport in double‐porosity (DP) media under unsaturated water flow conditions is presented. The macroscopic model was developed by applying the asymptotic homogenization method. It is based on theoretical and empirical considerations dealing with the orders of magnitude of characteristic quantities involved in the process. For this purpose a physical model that mimics the behavior of DP medium was built. The resulting two‐equation model relies on a coupling exchange term between micro‐ and macro‐porosity subdomains associated with local non‐equilibrium solute concentrations. The model was numerically implemented (Comsol Multiphysics^®) to simulate the macroscopic one‐dimensional physical process taking place into the porous medium of 3D periodic microstructure. A series of dispersion experiments of NaCl solution under unsaturated steady‐state flow conditions were performed. The experimental results were used first to calibrate the dispersion coefficient of the model, and second to validate it through two other independent experiments. The excellent agreement between the numerical simulations and the measurements of the time evolution of the non‐symmetrical breakthrough curves provides a proof of predictive capacity of the developed model. Copyright © 2010 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.     

A bounding surface elasto‐viscoplastic constitutive model for non‐isothermal cyclic analysis of asphaltic materials

An elasto‐viscoplastic constitutive model for asphaltic materials is presented within the context of bounding surface plasticity theory, taking into account the effects of the stress state, void binder degree of saturation, temperature and strain rate on the material behaviour. A stress state dependent non‐linear elasticity model is introduced to represent time‐independent recoverable portion of the deformation. The consistent visco‐plasticity framework is utilised to capture the rate‐dependent, non‐recoverable strain components. The material parameters introduced in the model are identified, and their determination from conventional laboratory tests is discussed. The capability of the model to reproduce experimentally observed response of asphaltic materials is demonstrated through numerical simulations of several laboratory test data from the literature. Copyright © 2016 John Wiley & Sons, Ltd.     

The role of non‐coaxiality in the simulation of strain localization based on classical and Cosserat continua

It is normally accepted that materials inside the shear band undergo severe rotation of the principal stress direction, which causes non‐coaxiality between the principal stress and principal plastic strain rate. However, classical plasticity flow theory implicitly assumes that the principal stress and the principal plastic strain rate are coaxial; thus, it may not correctly predict the onset of the shear band. In addition, classical continuum does not contain any internal length scales; as a result, it cannot provide a reasonable shear band thickness. In this study, the original vertex non‐coaxial plastic model based on the classical continuum is extended to the Cosserat continuum. The corresponding codes are implemented via the interface of the user defined element subroutine in ABAQUS. Through a simple shear test, the effectiveness of the user's codes is verified. Through a uniaxial compression test, the influence of non‐coaxiality on the onset, the orientation, and the thickness of the shear band is investigated. Results show that the onset of the shear localization is delayed, and the thickness of the shear band is widened when the non‐coaxial degree increases, while the orientation of the shear band is little affected by the non‐coaxial degree. In addition, it is found that the non‐coaxiality can weaken the micro‐polar effect to some extent; nonetheless, the Cosserat non‐coaxial model still has its advantage over the classical non‐coaxial model in capturing the pre‐bifurcation as well as the post‐bifurcation behaviors of strain localization. Copyright © 2016 John Wiley & Sons, Ltd.     

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.     

结构性砂土静力触探试验离散元分析

为研究结构性砂土静力触探的宏微观力学特性,在10 g重力场中生成一个净砂地基;将一个考虑胶结厚度的微观胶结接触模型引入到净砂地基中以生成结构性砂土地基;用一定速率移动探杆以模拟结构性砂土中的静力触探过程,其中,探杆由4面刚性墙组成。结果表明,随着贯入深度的增加,锥尖贯入阻力逐渐增大,增长速度逐渐减慢,在达到临界深度后贯入阻力在某一定值附近波动;锥尖部位有明显的力链集中现象,力链的集中程度和范围会随着贯入深度的增加而逐渐提高和扩大;静力触探过程中,探杆两侧的土体经历了明显的加载和卸载过程,且土体主应力方向发生偏转;离探杆越远,主应力偏转速度越慢,最终偏转角越小;不同深度处平均纯转动率（APR）的变化趋势基本相同,而APR最大值会随着土体深度而逐渐增加;探杆的贯入会使土颗粒间胶结发生破坏,胶结的破坏形式主要有拉剪破坏和压剪破坏两种,而拉剪破坏数目要比压剪破坏数目多。     

Discrete element method modeling of inherently anisotropic rocks under uniaxial compression loading

A new numerical approach is proposed in this study to model the mechanical behaviors of inherently anisotropic rocks in which the rock matrix is represented as bonded particle model, and the intrinsic anisotropy is imposed by replacing any parallel bonds dipping within a certain angle range with smooth‐joint contacts. A series of numerical models with β = 0°, 15°, 30°, 45°, 60°, 75°, and 90° are constructed and tested (β is defined as the angle between the normal of weak layers and the maximum principal stress direction). The effect of smooth‐joint parameters on the uniaxial compression strength and Young's modulus is investigated systematically. The simulation results reveal that the normal strength of smooth‐joint mainly affects the behaviors at high anisotropy angles (β > 45°), while the shear strength plays an important role at medium anisotropy angles (30°–75°). The normal stiffness controls the mechanical behaviors at low anisotropy angles. The angle range of parallel bonds being replaced plays an important role on defining the degree of anisotropy. Step‐by‐step procedures for the calibration of micro parameters are recommended. The numerical model is calibrated to reproduce the behaviors of different anisotropic rocks. Detailed analyses are conducted to investigate the brittle failure process by looking at stress‐strain behaviors, increment of micro cracks, initiation and propagation of fractures. Most of these responses agree well with previous experimental findings and can provide new insights into the micro mechanisms related to the anisotropic deformation and failure behaviors. The numerical approach is then applied to simulate the stress‐induced borehole breakouts in anisotropic rock formations at reduced scale. The effect of rock anisotropy and stress anisotropy can be captured. Copyright © 2015 John Wiley & Sons, Ltd.     

Experimental study and constitutive modelling of elasto‐plastic damage in heat‐treated mortar

This study investigates the effect of a heat‐treatment upon the thermo‐mechanical behaviour of a model cement‐based material, i.e. a normalized mortar, with a (w/c) ratio of 0.5. First, a whole set of varied experimental results is provided, in order to either identify or validate a thermo‐mechanical constitutive model, presented in the second paper part. Experimental responses of both hydraulic and mechanical behaviour are given after different heating/cooling cycling levels (105, 200, 300, 400^°C). The reference state, used for comparison purposes, is taken after mass stabilization at 60^°C. Typical uniaxial compression tests are provided, and original triaxial deviatoric compressive test responses are also given. Hydraulic behaviour is identified simultaneously to triaxial deviatoric compressive loading through gas permeability K_gas assessment. K_gas is well correlated with volumetric strain evolution: gas permeability increases hugely when ε_v testifies of a dilatant material behaviour, instead of contractile from the test start. Finally, the thermo‐mechanical model, based on a thermodynamics approach, is identified using the experimental results on uniaxial and triaxial deviatoric compression. It is also positively validated at residual state for triaxial deviatoric compression, but also by using a different stress path in lateral extension, which is at the origin of noticeable plasticity. Copyright © 2009 John Wiley & Sons, Ltd.     

Happy Jack Uraninite: A New Reference Material for High Spatial Resolution Analysis of U‐Rich Matrices

There is currently a lack of well‐characterised matrix‐matched reference materials (RMs) for forensic analysis of U‐rich materials at high spatial resolution. This study reports a detailed characterisation of uraninite (nominally UO_2+x) from the Happy Jack Mine (UT, USA). The Happy Jack uraninite can be used as a RM for the determination of rare earth element (REE) mass fractions in nuclear materials, which provide critical information for source attribution purposes. This investigation includes powder X‐ray diffraction (pXRD) data, as well as major, minor and trace element abundances determined using a variety of micro‐analytical techniques. The chemical signature of the uraninite was investigated at the macro (cm)‐scale with micro‐X‐ray fluorescence (µXRF) mapping and at high spatial resolution (tens of micrometre scale) using electron probe microanalysis (EPMA) and laser ablation‐inductively coupled plasma‐mass spectrometry (LA‐ICP‐MS) analyses. Based on EPMA results, the uraninite is characterised by homogeneous UO₂ and CaO contents of 91.57 ± 1.49% m/m (2s uncertainty) and 2.70 ± 0.38% m/m (2s), respectively. Therefore, CaO abundances were used as the internal standard when conducting LA‐ICP‐MS analyses. Overall, the major element and REE compositions are homogeneous at both the centimetre and micrometre scales, allowing this material to be used as a RM for high spatial resolution analysis of U‐rich samples.     

A non‐coaxial critical state soil model and its application to simple shear simulations

The yield vertex non‐coaxial theory is implemented into a critical state soil model, CASM (Int. J. Numer. Anal. Meth. Geomech. 1998; 22 :621–653) to investigate the non‐coaxial influences on the stress–strain simulations of real soil behaviour in the presence of principal stress rotations. The CASM is a unified clay and sand model, developed based on the soil critical state concept and the state parameter concept. Without loss of simplicity, it is capable of simulating the behaviour of sands and clays within a wide range of densities. The non‐coaxial CASM is employed to simulate the simple shear responses of Erksak sand and Weald clay under different densities and initial stress states. Dependence of the soil behaviour on the Lode angle and different plastic flow rules in the deviatoric plane are also considered in the study of non‐coaxial influences. All the predictions indicate that the use of the non‐coaxial model makes the orientations of the principal stress and the principal strain rate different during the early stage of shearing, and they approach the same ultimate values with an increase in loading. These ultimate orientations are dependent on the density of soils, and independent of their initial stress states. The use of the non‐coaxial model also softens the shear stress evolutions, compared with the coaxial model. It is also found that the ultimate shear strengths by using the coaxial and non‐coaxial models are dependent on the plastic flow rules in the deviatoric plane. Copyright © 2006 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.     

Integration and calibration of a plasticity model for granular materials

A new macroscopic constitutive model for non‐cohesive granular materials, with the focus on coarse‐sized materials (railway ballast), is presented. The model is based on the concepts of rate‐independent isotropic plasticity. The Backward Euler rule is used for integrating the pertinent evolution equations. The resulting incremental relations are solved in the strain space that is extended with the internal (hardening) variables. The model is calibrated using data from Conventional Triaxial Compression (CTC) tests, carried out at the University of Colorado at Boulder. A function evaluation method is used for the optimization, whereby a ‘multi‐vector’ strategy for choosing the appropriate start vector is proposed. Copyright © 2002 John Wiley & Sons, Ltd.     

Micromechanical approach to effective viscoelastic properties of micro‐fractured geomaterials

The aim of this paper is to formulate a micromechanics‐based approach to non‐aging viscoelastic behavior of materials with randomly distributed micro‐fractures. Unlike cracks, fractures are discontinuities that are able to transfer stresses and can therefore be regarded from a mechanical viewpoint as interfaces endowed with a specific behavior under normal and shear loading. Making use of the elastic‐viscoelastic correspondence principle together with a Mori‐Tanka homogenization scheme, the effective viscoelastic behavior is assessed from properties of the material constituents and damage parameters related to density and size of fractures. It is notably shown that the homogenized behavior thus formulated can be described in most cases by means of a generalized Maxwell rheological model. For practical implementation in structural analyses, an approximate model for the isotropic homogenized fractured medium is formulated within the class of Burger models. Although the approximation is basically developed for short‐term and long‐term behaviors, numerical applications indicate that the approximate Burger model accurately reproduce the homogenized viscoelastic behavior also in the transient conditions.     

Refinement of the Micro‐Distillation Technique for Isotopic Analysis of Geological Samples with pg‐Level Osmium Contents

In recent years, the ¹⁸⁷Re–¹⁸⁷Os isotope system has been increasingly used to study samples containing very small quantities of Os. For such samples, optimisation of measurement procedures is essential to minimise the loss of Os before mass spectrometric measurements. Micro‐distillation is a necessary purification step that is applied after the main Os chemical separation procedure, prior to Os isotope ratio measurements by negative‐thermal ionisation mass spectrometry (N‐TIMS). However, unlike the other separation steps, this procedure has not yet been optimised for small samples. In this study, we present a refined micro‐distillation method that achieved higher yields and allowed high‐precision R(¹⁸⁷Os/¹⁸⁸Os) expressed as ¹⁸⁷Os/¹⁸⁸Os measurements for small‐sized geological samples that contain only a few pg Os. The Os recovery in the micro‐distillation step was tested by changing the operating conditions including heating time and temperature, and amounts of oxidant and reductant. Recoveries were measured by the isotope dilution ICP‐MS method after the addition of ¹⁹⁰Os‐enriched spike solution. We found that the most critical factor controlling the chemical yield of Os during micro‐distillation is the extent of dilution of the reductant (HBr) by H₂O evaporated from the oxidant. A refined micro‐distillation method, in which the amount of oxidant solution is reduced from the conventional method, achieved an improved chemical yield of Os (~ 90%). This refined method was applied to the measurement of ¹⁸⁷Os/¹⁸⁸Os by N‐TIMS of varying test portions of the geological reference material BIR‐1a. The resulting ¹⁸⁷Os/¹⁸⁸Os ratios of BIR‐1a matched the literature data, with propagated uncertainties of 0.2, 1.1 and 11% digested sample quantities containing 150, 10 and 1 pg of Os, respectively.     

Role of inertia in falling head permeability test

This paper discusses the role of macro‐ and micro‐inertial forces in the falling head permeability test. The model of flow including viscous and dynamic interaction forces is formulated. Then, the model is used to perform numerical simulations of the tests for high‐permeability soils. The presented results prove the existence of a level of permeability of materials above which macro‐and micro‐inertial forces are important and that in all cases the latter forces are dominant. The results suggest limited validity of the standard interpretation of the falling head permeability test and possible usefulness of the test to determine both permeability and the parameter characterizing non‐linear interactions of fluid and solid skeleton. Copyright © 2009 John Wiley & Sons, Ltd.     

Extension of plasticity theory to debonding,grain dissolution,and chemical damage of calcarenites

The mechanical properties of calcarenites are known to be significantly affected by water saturation: both stiffness and strength decrease for wetting in the short term and for chemical dissolution in the long term. Both processes mainly affect bonds among grains: immediately after inundation depositional bonds fall in suspension, whereas diagenetic bonds dissolve more slowly. In this paper, the authors started from the micro‐structural analysis of the weathering processes to conceive a strain hardening hydro‐chemo‐mechanical coupled elastoplastic constitutive model. The concept of extended hardening rules is here enriched: weathering functions have been determined by employing a micro to macro simplified upscaling procedure. Chemical damage is incorporated into the formulation by means of a scalar damage function. Its evolution is also described by using a multiscale approach. A new term is added to the strain rate tensor in order to incorporate the dissolution induced chemical deformations developing once the soft rock is turned into a granular material. A calibration procedure for the constitutive parameters is suggested, and the model is validated by using both coupled and uncoupled chemo‐mechanical experimental test results. Copyright © 2015 John Wiley & Sons, Ltd.     

Cyclic macro‐element for soil–structure interaction: material and geometrical non‐linearities

This paper presents a non‐linear soil–structure interaction (SSI) macro‐element for shallow foundation on cohesive soil. The element describes the behaviour in the near field of the foundation under cyclic loading, reproducing the material non‐linearities of the soil under the foundation (yielding) as well as the geometrical non‐linearities (uplift) at the soil–structure interface. The overall behaviour in the soil and at the interface is reduced to its action on the foundation. The macro‐element consists of a non‐linear joint element, expressed in generalised variables, i.e. in forces applied to the foundation and in the corresponding displacements. Failure is described by the interaction diagram of the ultimate bearing capacity of the foundation under combined loads. Mechanisms of yielding and uplift are modelled through a global, coupled plasticity–uplift model. The cyclic model is dedicated to modelling the dynamic response of structures subjected to seismic action. Thus, it is especially suited to combined loading developed during this kind of motion. Comparisons of cyclic results obtained from the macro‐element and from a FE modelization are shown in order to demonstrate the relevance of the proposed model and its predictive ability. Copyright © 2001 John Wiley & Sons, Ltd.