Bifurcation modeling in geomaterials: From the second‐order work criterion to spectral analyses

The present paper investigates bifurcation analysis based on the second‐order work criterion, in the framework of rate‐independent constitutive models and rate‐independent boundary‐value problems. The approach applies mainly to nonassociated materials such as soils, rocks, and concretes. The bifurcation analysis usually performed at the material point level is extended to quasi‐static boundary‐value problems, by considering the stiffness matrix arising from finite element discretization. Lyapunov's definition of stability (Annales de la faculté des sciences de Toulouse 1907; 9 :203–274), as well as definitions of bifurcation criteria (Rice's localization criterion (Theoretical and Applied Mechanics. Fourteenth IUTAM Congress, Amsterdam, 1976; 207–220) and the plasticity limit criterion are revived in order to clarify the application field of the second‐order work criterion and to contrast these criteria. The first part of this paper analyses the second‐order work criterion at the material point level. The bifurcation domain is presented in the 3D stress space as well as 3D cones of unstable loading directions for an incrementally nonlinear constitutive model. The relevance of this criterion, when the nonlinear constitutive model is expressed in the classical form (dσ = Mdε) or in the dual form (dε = Ndσ), is discussed. In the second part, the analysis is extended to the boundary‐value problems in quasi‐static conditions. Nonlinear finite element computations are performed and the global tangent stiffness matrix is analyzed. For several examples, the eigenvector associated with the first vanishing eigenvalue of the symmetrical part of the stiffness matrix gives an accurate estimation of the failure mode in the homogeneous and nonhomogeneous boundary‐value problem. Copyright © 2008 John Wiley & Sons, Ltd.     

Numerical study of inter‐particle bond failure by 3D discrete element method

The cohesive‐frictional nature of cementitious geomaterials raises great interest in the discrete element method (DEM) simulation of their mechanical behavior, where a proper bond failure criterion is usually required. In this paper, the failure of bond material between two spheres was investigated numerically using DEM that can easily reproduce the failure process of brittle material. In the DEM simulations, a bonded‐grain system (composed of two particles and bond material in between) was discretized as a cylindrical assembly of very fine particles connecting two large end spheres. Then, the bonded‐grain system was subjected to compression/tension, shear, rolling and torsion loadings and their combinations until overall failure (peak state) was reached. Bonded‐grain systems with various sizes were employed to investigate bond geometry effects. The numerical results show that the compression strength is highly affected by bond geometry, with the tensile strength being dependent to a lesser degree. The shear, rolling and torsion strengths are all normal force dependent; i.e., with an increase in the normal force, these strengths first increase at a declining rate and then start to decrease upon the normal force exceeding a critical value. The combined actions of shear force, rolling moment and torque lead to a spherical failure envelope in a normalized loading space. The fitted bond geometry factors and bond failure envelopes obtained numerically in this three‐dimensional study are qualitatively consistent with those in previous two‐dimensional experiments. The obtained bond failure criterion can be incorporated into a future bond contact model. Copyright © 2015 John Wiley & Sons, Ltd.     

A stability criterion for elasto‐viscoplastic constitutive relationships

In this paper, the onset of mechanical instability in time‐sensitive elasto‐viscoplastic solids is theoretically analyzed at the constitutive level and associated with the occurrence of ‘spontaneous accelerations’ under stationary external perturbations. For this purpose, a second‐order form of Perzyna's constitutive equations is first derived by time differentiation, and a sufficient stability condition is identified for general mixed loading programs. These loading conditions are in fact the most general in both laboratory tests and real boundary value problems, where a combination of certain stress and strain components is known/prescribed. The theoretical analysis leads to find precise stability limits in terms of material hardening modulus. In the case of constitutive relationships with isotropic strain‐hardening, no instabilities are possible while the hardening modulus is larger than the so‐called ‘controllability modulus’ defined for (inviscid) elasto‐plastic materials. It is also shown that the current stress/strain rate may also directly influence the occurrence of elasto‐viscoplastic instability, which is at variance with elasto‐plastic inviscid media. Copyright © 2015 John Wiley & Sons, Ltd.     

An original semi‐discrete approach to assess gas conductivity of concrete structures

For civil engineering structures with a tightness role, structural permeability is a key issue. In this context, this paper presents a new proposition of a numerical modelling of leakage rate through a cracked concrete structure undergoing mode I cracking. The mechanical state of the material, considered in the framework of continuum mechanics based on finite element modelling, is described by means of the stress‐based nonlocal damage model which takes into account the stress state and provides realistic local mechanical fields. A semi‐discrete method based on the strong discontinuity approach to estimate crack opening is then considered in the post‐treatment phase. Using a Poiseuille's like relation, the coupling between the mechanical state of the material and its dry gas conductivity is performed. For validation purposes, an original experimental campaign is conducted on a dry concrete disc loaded in a splitting setup. During the loading, gas conductivity and digital image correlation analysis are performed. The comparison with the 3D experimental mechanical global response highlights the performance of the mechanical model. The comparison between crack openings measured by digital image correlation and estimated by the strong discontinuity method shows a good agreement. Finally, the results of the semi‐discrete approach coupled with the gas conductivity compared with experimental data show a good estimation of the structural conductivity. Consequently, if the mechanical problem is well modelled at the global scale, then the proposed approach provides good estimation of gas conductivity. Copyright © 2016 John Wiley & Sons, Ltd.     

A micromechanical finite element model for linear and damage‐coupled viscoelastic behaviour of asphalt mixture

This study presents a finite element (FE) micromechanical modelling approach for the simulation of linear and damage‐coupled viscoelastic behaviour of asphalt mixture. Asphalt mixture is a composite material of graded aggregates bound with mastic (asphalt and fine aggregates). The microstructural model of asphalt mixture incorporates an equivalent lattice network structure whereby intergranular load transfer is simulated through an effective asphalt mastic zone. The finite element model integrates the ABAQUS user material subroutine with continuum elements for the effective asphalt mastic and rigid body elements for each aggregate. A unified approach is proposed using Schapery non‐linear viscoelastic model for the rate‐independent and rate‐dependent damage behaviour. A finite element incremental algorithm with a recursive relationship for three‐dimensional (3D) linear and damage‐coupled viscoelastic behaviour is developed. This algorithm is used in a 3D user‐defined material model for the asphalt mastic to predict global linear and damage‐coupled viscoelastic behaviour of asphalt mixture. For linear viscoelastic study, the creep stiffnesses of mastic and asphalt mixture at different temperatures are measured in laboratory. A regression‐fitting method is employed to calibrate generalized Maxwell models with Prony series and generate master stiffness curves for mastic and asphalt mixture. A computational model is developed with image analysis of sectioned surface of a test specimen. The viscoelastic prediction of mixture creep stiffness with the calibrated mastic material parameters is compared with mixture master stiffness curve over a reduced time period. In regard to damage‐coupled viscoelastic behaviour, cyclic loading responses of linear and rate‐independent damage‐coupled viscoelastic materials are compared. Effects of particular microstructure parameters on the rate‐independent damage‐coupled viscoelastic behaviour are also investigated with finite element simulations of asphalt numerical samples. Further study describes loading rate effects on the asphalt viscoelastic properties and rate‐dependent damage behaviour. Copyright © 2006 John Wiley & Sons, Ltd.     

A simple three‐dimensional distinct element modeling of the mechanical behavior of bonded sands

This paper presents a simple three‐dimensional (3D) Distinct Element Method (DEM) for numerical simulation of the mechanical behavior of bonded sands. First, a series of micro‐mechanical tests on a pair of aluminum rods glued together by cement with different bond sizes were performed to obtain the contact mechanical responses of ideally bonded granular material. Second, a 3D bond contact model, which takes into account the influences of bond sizes, was established by extending the obtained 2D experimental results to 3D case. Then, a DEM incorporating the new contact model was employed to perform a set of drained triaxial compression tests on the DEM bonded specimens with different cement contents under different confining pressures. Finally, the mechanical behavior of the bonded specimens was compared with the available experimental results. The results show that the DEM incorporating the simple 3D bond contact model is able to capture the main mechanical behavior of bonded sands. The bonded specimen with higher cement content under lower confining pressure exhibits more pronounced strain softening and shear dilatancy. The peak and residual strengths, the apparent cohesion and peak/residual friction angles, and the position and slope of the critical state line increase with increase in cement content. Microscopically, bond breakage starts when the system starts to dilate and the maximum rate of bond breakage coincides with the maximum rate of dilation. Bond breakage is primarily due to tension‐shear failure and the percentage of such failures is independent of both confining pressure and cement content. Copyright © 2015 John Wiley & Sons, Ltd.     

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.     

An enhanced constitutive model for crushable granular materials

Studies in the past have tried to reproduce the mechanical behaviour of granular materials by proposing constitutive relations based on a common assumption that model parameters and parameters describing the properties, including gradation of individual grains are inevitably linked. However successful these models have proved to be, they cannot account for the changes in granular assembly behaviour if the grains start to break during mechanical loading. This paper proposes to analyse the relation between grading change and the mechanical behaviour of granular assembly. A way to model the influence of grain breakage is to use a critical state‐based model. The influence of the amount of grain breakage during loading, depending on the individual grain strength and size distribution, can be introduced into constitutive relations by means of a new parameter that controls the evolution of critical state with changes in grain size distribution. Experimental data from a calcareous sand, a quartz sand, and a rockfill material were compared with numerical results and good‐quality simulations were obtained. The main consequences of grain breakage are increased compressibility and a gradual dilatancy disappearance in the granular material. The critical state concept is also enriched by considering its overall relation to the evolution of the granular material. Copyright © 2009 John Wiley & Sons, Ltd.     

Nonlinear SEM numerical analyses of dry dense sand specimens under rapid and dynamic loading

The paper mainly concerns the mechanical response of 2D dry dense sand specimens under shock loading. The problem is numerically analysed by means of a SEM dynamic code, within which an already conceived non‐local viscoplastic constitutive model characterized by a non‐associated flow rule and by an anisotropic strain hardening has been implemented. In particular the strain localization and time dependency of the material mechanical response are taken into consideration. Both rapid/static loading and dynamic histories are numerically simulated. In the first case, the time dependency of the material mechanical response can be captured by neglecting inertial effects, while in the second one the two factors are superimposed and the propagation of the stress waves within the specimen is discussed. The interest of these analyses derives from the fact that the diffusion phenomenon takes place within a specimen already localized. Copyright © 2006 John Wiley & Sons, Ltd.     

Dynamic Model Re-construction and Corresponding Abnormal Characteristics Anatomy of Subtropical High based on Actual DataTime SeriesSamples

Aiming at the abnormal activity of West Pacific Subtropical High(WPSH) inducing the serious flood disaster in 1998, and considering the complex seasonal activity of WPSH in summer, based on the actual geopotential height field time series observational data in 1998, a non-linear dynamic model of WPSH was reconstructed by using the idea of dynamic system retrieval and genetic algorithms. Combined with the synchronous actual weather process, the corresponding dynamic characteristics and variation mechanism of WPSH, such as the equilibrium state destabilization and bifurcation of the WPSH system, were analyzed and discussed. The research showed that due to the advantages of global optimization and parallel calculation of GA, a nonlinear dynamic model could be reasonably and quickly reconstructed, and the reconstructed nonlinear dynamic model of WPSH derived from the actual abnormal case data could be used for describing objectively and explaining reasonably the basic characteristics and the abnormal activity process of WPSH. The dynamic discussion showed that the configuration and diversification of WPSH system equilibriums, such as bifurcation, catastrophe had a good matching with the actual medium short period abnormal activity of WPSH, such as northward jumping and southward droping, as well as double-ridge line pattern. The modality of WPSH geopotential height field induced by multi-equilibriums destabilization was more complex than that of by single-equilibriums bifurcation, and their exhibition patterns were more different. With the forcing parameters increaseing, closing with and exceeding critical value, the WPSH system equilibrium state destabilization or bifurcation would occur. Especially, when the equilibrium state jump/drop from low/high to the high/low value, corresponding WPSH system behaved as a leap northward /fall southward, moreover, when the stable equilibriums state from 1 to 2, corresponds with the pattern of “double-ridges” phenomenon of WPSH, contrariwise, “double-ridges” phenomenon would disappear. A meaningful technique route was provided to the mechanism research and dynamic characteristic discussion for the complex weather system of hardly modeling such as WPSH.     

Micromechanical viscoelasto‐plastic models and finite element implementation for rate‐independent and rate‐dependent permanent deformation of stone‐based materials

This paper presents parallel and serial viscoelasto‐plastic models to simulate the rate‐independent and the rate‐dependent permanent deformation of stone‐based materials, respectively. The generalized Maxwell viscoelastic and Chaboche's plastic models were employed to formulate the proposed parallel and serial viscoelasto‐plastic constitutive laws. The finite element (FE) implementation of the parallel model used a displacement‐based incremental formulation for the viscoelastic part and an elastic predictor—plastic corrector scheme for the elastoplastic component. The FE framework of the serial viscoelasto‐plastic model employed a viscoelastic predictor—plastic corrector algorithm. The stone‐based materials are consisted of irregular aggregates, matrix and air voids. This study used asphalt mixtures as an example. A digital sample was generated with imaging analysis from an optically scanned surface image of an asphalt mixture specimen. The modeling scheme employed continuum elements to mesh the effective matrix, and rigid bodies for aggregates. The ABAQUS user material subroutines defined with the proposed viscoelasto‐plastic matrix models were employed. The micromechanical FE simulations were conducted on the digital mixture sample with the viscoelasto‐plastic matrix models. The simulation results showed that the serial viscoelasto‐plastic matrix model generated more permanent deformation than the parallel one by using the identical material parameters and displacement loadings. The effect of loading rates on the material viscoelastic and viscoelasto‐plastic mixture behaviors was investigated. Permanent deformations under cyclic loadings were determined with FE simulations. The comparison studies showed that the simulation results correctly predicted the rate‐independent and rate‐dependent viscoelasto‐plastic constitutive properties of the proposed matrix models. Overall, these studies indicated that the developed micromechanical FE models have the abilities to predict the global viscoelasto‐plastic behaviors of the stone‐based materials. Copyright © 2009 John Wiley & Sons, Ltd.     

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.     

关于固体地球系统与地质作用的三点认识

文中讨论了作为极复杂巨系统的固体地球内部运动规律的 3条推论 :(1)地球系统在偏离平衡态时 ,有选择能量耗散最小方式的“惯性” ,在系统内禀熵急增时或者回到一种准稳定的定态 ,或者通过自组织迅速减少内禀熵增加率。换言之 ,地球系统可能有种尽量保存自身内能不受大量耗散的惯性 ,使其总体内禀熵产生率对时空的积分取极小值。我们把这一认识称为总体熵产生率取极小准则 ,对于孤立与渐变的地质过程 ,它等同于热力学第二定律与普里高津最小熵产生原理 ,对于极复杂地球巨系统及远离平衡态情况加了“总体”两字 ,意思是在局部或短期突变时熵产生率可能是大的 ,但它对全球与长时间的积分而言仍然取极小。换句话说 ,作为一个巨系统 ,固体地球在其局部远离平衡状态时 ,仍然能保持总体上内能消耗取极小的惯性 ,维持对全部固体地球时空系统总体熵产生率取极小准则。 (2 )关于地质作用过程演化的定态遍历准则。地球系统的复杂性不仅表现在其非线性 ,即同时存在着多种可能的定态作为其演化的趋势 ,而且表现在其时空发展过程中将经历尽可能多的定态。非线性非平衡态动力学的一般规律符合大多数包含激变事件的地质作用过程。将来在地球系统偏离平衡态足够远时 ,它可能会具有无穷多个耗散结构 ,因而使系统进入完全?     

Improved Bell's method for the stability analysis of slopes

A new procedure based on the approximation to the total normal pressure along the slip surface is developed to compute the factor of safety of slopes for slip surfaces of all shapes. By taking the whole sliding body, instead of an individual slice, as the loaded object, all the equilibrium equations are formulated according to the three‐moment equilibrium conditions rather than the two force equilibrium conditions and one‐moment equilibrium condition. The system of nonlinear equations deduced in this way is well‐scaled and enjoys excellent numerical properties such as the existence of solution with a positive factor of safety, a nearly unlimited scope of convergence and a rapid convergence rate associated with the Newton method. In the case of ?_u =0—the situation where no drainage and no consolidation are involved, furthermore, the system has a unique solution and the factor of safety has an explicit expression. Some typical examples are analyzed to demonstrate the numerical properties of the proposed procedure. Copyright © 2009 John Wiley & Sons, Ltd.     

Localized failure of fibre‐reinforced elastic–plastic materials subjected to plane strain loading

We consider discontinuous bifurcations as the indicator of a localized failure for a class of composites that are characterized by elastic fibres reinforcing an elastic–plastic matrix. A macroscopic tangent stiffness tensor for the fibre‐reinforced composite is developed by consistently homogenizing the contribution of fibres in a spherical representative volume element. Analytical solutions are derived for the critical hardening modulus and corresponding bifurcation directions for the case of plane strain loading. Properties of the solutions are further illustrated on the example of the non‐associated Drucker–Prager model at onset of yielding. Results show that presence of fibres decreases the critical hardening modulus, thus inhibiting the onset of strain localization. The rate of decrease in the critical hardening modulus is the highest for pure shear, followed by uniaxial tension, uniaxial compression, biaxial tension and biaxial compression. The main fibre parameters that control the onset of strain localization are their volumetric content and their stiffness modulus whereby very stiff fibres can produce the most significant decrease in the critical hardening modulus, especially for the state of biaxial tension. The critical hardening modulus for the non‐associated Drucker–Prager model exhibits a full range of localization modes including compaction bands, dilation bands, and transition in the form of shear bands regardless of the presence of fibres. Presence of fibres affects bifurcation directions, except in the case when Poisson's ratio of the matrix is equal to 0.25. The results demonstrate stabilizing effects of fibres by which they provide the control against the onset of strain localization. Copyright © 2006 John Wiley & Sons, Ltd.     

Drained instability analysis of sand under biaxial loading using a 3D material model

The diffused and localized instabilities in sand under drained biaxial loading have been analyzed here following a plane strain bifurcation framework, where the rate independent material is defined using a generalized 3D non-associative constitutive model. This study is focused on how various instability modes emerge with respect to initial density, confining pressure, and the applied boundary conditions. Results from large deformation framework have been compared with those from small deformation approximation and the later is noticed to fail in capturing the emergence of diffused modes and predicts delayed onset of localization. The theoretical predictions compares well with existing experimental observations.     

A constitutive model for the dynamic and high‐pressure behaviour of a propellant‐like material: Part I: Experimental background and general structure of the model

This paper is the first part of a work that aims at developing a mechanical model for the behaviour of propellant‐like materials under high confining pressure and strain rate. The behaviour of a typical material is investigated experimentally. Several microstructural deformation processes are identified and correlated with loading conditions. The resulting behaviour is complex, non‐linear, and characterized by multiple couplings. The general structure of a relevant model is sought using a thermodynamic framework. A viscoelastic‐viscoplastic‐compaction model structure is derived under suitable simplifying assumptions, in the framework of finite, though moderate, strains. Model development, identification and numerical applications are given in the companion paper. Copyright © 2001 John Wiley & Sons, Ltd.     

Diffuse instabilities with transition to localization in loose granular materials

This paper is concerned with diffuse and other ensuing failure modes in geomaterials when tested under homogeneous states of shearing in various loading programs and drainage conditions. Material instability is indeed the basic property that accounts for the instability of an initially homogeneous deformation field leading to diffuse failure and strain localization in geomaterials. The former is normally characterized by a runaway type of failure accompanied with a sudden and violent collapse of the material in the absence of any localization phenomena. Against this backdrop, we present a brief overview of material instability in elastoplastic solids where one finds a rich source of theoretical concepts including bifurcation, strain localization, diffuse failure and second‐order work, as well as a considerable body of experiments. Some compelling laboratory experimental studies of material instability with focus to diffuse failure are then presented and interpreted based on the second‐order work. Finally, various material instability analyses using an elastoplastic constitutive and a general finite element analysis of the above‐mentioned laboratory experimental tests are presented as a boundary value problem. It is shown that instability can be captured from otherwise uniform stress, density and hydraulic states, whereas uniform deviatoric loads are being applied on the external boundaries of a specimen. Although the numerical simulations reproduce well the laboratory experimental results, they also highlight the hierarchy of failure modes where localization phenomena emerge in the post‐bifurcation regime as a result of a degradation of homogeneity starting from a diffuse mode signalled by a zero second‐order work. Copyright © 2012 John Wiley & Sons, Ltd.     

公路连拱隧道土压力荷载的计算方法研究

针对目前国内公路连拱隧道尚无比较成熟的土压力荷载计算公式的现状,根据连拱隧道的受力特点,基于普氏平衡拱理论,推导了深埋连拱隧道的土压力荷载计算公式,同时结合规范中关于隧道深浅埋分界高度的规定,给出了深浅埋连拱隧道的分界埋深及荷载计算公式,形成了完整的连拱隧道深浅埋情况下的荷载计算体系。通过大量的分析计算并与其他公式对比,验证了推导出的公式的合理性,并将公式应用于湖北沪蓉西高速公路连拱隧道及分岔隧道的设计中,取得了良好的效果。     

Biaxial test simulations using a packing of polygonal particles

The mechanical response of cohesionless granular materials under monotonic loading is studied by performing molecular dynamic simulations. The diversity of shapes of soil grains is modelled by using randomly generated convex polygons as granular particles. Results of the biaxial test obtained for dense and loose media show that samples achieve the same void ratio at large strains independent of their initial density state. This limit state resembles the so‐called critical state of soil mechanics, except for some stress fluctuations, which remain for large deformations. These fluctuations are studied at the micro‐mechanical level, by following the evolution of the co‐ordination number, force chains and the fraction of the sliding contacts of the sample. Copyright © 2007 John Wiley & Sons, Ltd.