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.     

Consistent tangent matrices for density‐dependent finite plasticity models

The consistent tangent matrix for density‐dependent plastic models within the theory of isotropic multiplicative hyperelastoplasticity is presented here. Plastic equations expressed as general functions of the Kirchhoff stresses and density are considered. They include the Cauchy‐based plastic models as a particular case. The standard exponential return‐mapping algorithm is applied, with the density playing the role of a fixed parameter during the nonlinear plastic corrector problem. The consistent tangent matrix has the same structure as in the usual density‐independent plastic models. A simple additional term takes into account the influence of the density on the plastic corrector problem. Quadratic convergence results are shown for several representative examples involving geomaterial and powder constitutive models. Copyright © 2001 John Wiley & Sons, Ltd.     

Thermodynamic‐based model for coupling temperature‐dependent viscoelastic,viscoplastic, and viscodamage constitutive behavior of asphalt mixtures

Based on the continuum damage mechanics, a general and comprehensive thermodynamic‐based framework for coupling the temperature‐dependent viscoelastic, viscoplastic, and viscodamage behaviors of bituminous materials is presented. This general framework derives systematically Schapery‐type nonlinear viscoelasticity, Perzyna‐type viscoplasticity, and a viscodamage model analogous to the Perzyna‐type viscoplasticity. The resulting constitutive equations are implemented in the well‐known finite element code Abaqus via the user material subroutine UMAT. A systematic procedure for identifying the model parameters is discussed. Finally, the model is validated by comparing the model predictions with a comprehensive set of experimental data on hot mix asphalt that include creep‐recovery, creep, uniaxial constant strain rate, and repeated creep‐recovery tests in both tension and compression over a range of temperatures, stress levels, and strain rates. Comparisons between model predictions and experimental measurements show that the presented constitutive model is capable of predicting the nonlinear behavior of asphaltic mixes under different loading conditions. Copyright © 2011 John Wiley & Sons, Ltd.     

Continuum large cracking in a rate‐dependent plastic–damage model for cyclic‐loaded concrete structures

Formulation and algorithmic treatment of a rate‐dependent plastic–damage model modified to capture large tensile cracking in cyclic‐loaded concrete structures are presented in detail for a three‐dimensional implementation. The plastic–damage model proposed by Lee and Fenves in 1998 was founded based on isotropic damaged elasticity in combination with isotropic multi‐hardening plasticity to simulate cracking and crushing of concrete under cyclic or dynamic loadings. In order that the model can capture large crack opening displacements, which are inevitable in plain concrete structures, the excessive increase in plastic strain causing unrealistic results in cyclic behaviors is prevented when the tensile plastic–damage variable controlling the evolution of tensile damage is larger than a critical value. In such a condition, the crack opening/closing mechanism becomes similar to discrete cracking. The consistent tangent operator required to accelerate convergence rate is also formulated for the large cracking state including viscoplasticity. The validation and performance of the modified algorithm implemented in a special finite element program is exemplified through several single‐element tests as well as three structural applications. The last example examines the model in the seismic fracture analysis of Koyna dam as a benchmark problem and the resulting crack profile is compared with the available experiment. The numerical experimentations well demonstrate that the developed model whose modification is necessary to properly simulate the cyclic behavior of plain concrete subjected to large tensile strains is robust and reasonably accurate. Copyright © 2012 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.     

A visco‐elasto‐plastic model for granular materials under simple shear conditions

The numerical simulation of rapid landslides is quite complex mainly because constitutive models capable of simulating the mechanical behaviour of granular materials in the pre‐collapse and post‐collapse regimes are still missing. The goal of this paper is to introduce a constitutive model capable of capturing the response of dry granular flows from quasi‐static to dynamic conditions, in particular when the material experiences a sort of solid‐to‐fluid phase transition. An ideal assembly of identical spheres under simple shear conditions is considered. In the constitutive model, void ratio and granular temperature have been chosen as state variables, and both shear and normal stresses are computed as the sum of two contributions: the quasi‐static one and the collisional one. The former is determined by using a perfect elasto‐plastic model including the critical state concept, while the latter is derived from the kinetic theory of granular gases. The evolution of the granular temperature, fundamentally governing the material phase transition, is obtained by imposing the kinetic fluctuating energy balance. The constitutive relationship has been integrated, under both constant pressure and constant volume conditions, and the influence of shear strain rate, initial void ratio and normal pressure on the mechanical response has been investigated. Copyright © 2015 John Wiley & Sons, Ltd.     

Micromechanical model of weakly‐cemented sediments

Experiments indicate that the bulk stiffness of weakly‐cemented granular materials increases significantly with the applied load, a nonlinearity which is not captured by most micromechanical models. Here, we study the deformation of weakly‐cemented materials through grain‐scale simulations. We show that the frequently‐made assumption of linear‐elastic cement and intergranular contact area which varies according to a Hertzian‐type model cannot explain the observed nonlinearity of the bulk response. We introduce the micromechanical effect of mechanisms such as closure and opening of microcracks, granular cement and material heterogeneity through an effective contact stiffness which depends on the local deformation. We find that an exponential dependence between the effective stiffness and the local deformation, with a high exponent value, provides bulk stiffness which is in good agreement with experimental data. The inability of models with weaker intergranular stiffness‐deformation dependence to reproduce the experimental data demonstrates the highly‐nonlinear nature of the intergranular deformation. Our results highlight the importance of accounting for grain‐scale mechanisms in modeling granular materials, and provide a plausible explanation for the nonlinear behavior of weakly‐cemented sediments. Copyright © 2012 John Wiley & Sons, Ltd.     

A coupled chemo‐thermo‐hygro‐mechanical model of concrete at high temperature and failure analysis

A hierarchical mathematical model for analyses of coupled chemo‐thermo‐hygro‐mechanical behaviour in concretes at high temperature is presented. The concretes are modelled as unsaturated deforming reactive porous media filled with two immiscible pore fluids, i.e. the gas mixture and the liquid mixture, in immiscible–miscible levels. The thermo‐induced desalination process is particularly integrated into the model. The chemical effects of both the desalination and the dehydration processes on the material damage and the degradation of the material strength are taken into account. The mathematical model consists of a set of coupled, partial differential equations governing the mass balance of the dry air, the mass balance of the water species, the mass balance of the matrix components dissolved in the liquid phases, the enthalpy (energy) balance and momentum balance of the whole medium mixture. The governing equations, the state equations for the model and the constitutive laws used in the model are given. A mixed weak form for the finite element solution procedure is formulated for the numerical simulation of chemo‐thermo‐hygro‐mechanical behaviours. Special considerations are given to spatial discretization of hyperbolic equation with non‐self‐adjoint operator nature. Numerical results demonstrate the performance and the effectiveness of the proposed model and its numerical procedure in reproducing coupled chemo‐thermo‐hygro‐mechanical behaviour in concretes subjected to fire and thermal radiation. Copyright © 2005 John Wiley & Sons, Ltd.     

A rate‐dependent cohesive crack model based on anisotropic damage coupled to plasticity

In quasi‐brittle material the complex process of decohesion between particles in microcracks and localization of the displacement field into macrocracks is limited to a narrow fracture zone, and it is often modelled with cohesive crack models. Since the anisotropic nature of the decohesion process in separation and sliding is essential, it is particularly focused in this paper. Moreover, for cyclic and dynamic loading the unloading, load reversal (including crack closure) and rate dependency are essential features that are included in a new model. The modelling of degradation is based on a ‘localized’ version of anisotropic continuum damage coupled to inelasticity. The concept of strain energy equivalence between the states in the effective and nominal settings is adopted in order to define the free energy of the interface. The proposed fracture criterion is of the Mohr type, with a smooth transition of the failure and kinematics (slip and dilatation) characteristics between tension and shear. The chosen potential, of the Lemaitre‐type, for evolution of the dissipative processes is additively decomposed into plastic and damage parts, and non‐associative constitutive equations are obtained. The constitutive equations are integrated by applying the backward Euler rule and by using Newton iteration. The proposed model is assessed analytically and numerically and a typical calibration procedure for concrete is proposed. Copyright © 2006 John Wiley & Sons, Ltd.     

Template elastic‐plastic computations in geomechanics

In this paper we present a new approach to computations in elasto‐plastic geomechanics. The approach is based on the object oriented design philosophy and observations on similarity of most incremental elastic–plastic material models. This new approach to elastic–plastic computations in geomechanics allows for creation of template material models. The analysis of template material models will in turn allow for an easy implementation of other elastic–plastic material models based on the object oriented design principles. Furthermore we present some illustrative implementation details. Finally we present analysis results that emphasize features of template elastic–plastic computations in geomechanics. 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.     

Probabilistic yielding and cyclic behavior of geomaterials

In this paper, the novel concept of probabilistic yielding is used for 1‐D cyclic simulation of the constitutive behavior of geomaterials. Fokker–Planck–Kolmogorov equation‐based probabilistic elastic–plastic constitutive framework is applied for obtaining the complete probabilistic (probability density function) material response. Both perfectly plastic and hardening‐type material models are considered. It is shown that when uncertainties in material parameters are taken into consideration, even the simple, elastic‐perfectly plastic model captures some of the important features of geomaterial behavior, for example, modulus reduction with cyclic strain, which, deterministically, is only possible with more advanced constitutive models. Furthermore, it is also shown that the use of isotropic and kinematic hardening rules does not significantly improve the probabilistic material response. Copyright © 2010 John Wiley & Sons, Ltd.     

An approach for providing quasi‐convexity to yield functions and a generalized implicit integration scheme for isotropic constitutive models based on 2 unknowns

Yield and plastic potential surfaces are often affected by problems related to convexity. One such problem is encountered when the yield surface that bounds the elastic domain is itself convex; however, convexity is lost when the surface expands to pass through stress points outside the current elastic domain. In this paper, a technique is proposed, which effectively corrects this problem by providing linear homothetic expansion with respect to the centre of the yield surface. A very compact implicit integration scheme is also presented, which is of general applicability for isotropic constitutive models, provided that their yield and plastic potential functions are based on a separate mathematical definition of the meridional and deviatoric sections and that stress invariants are adopted as mechanical quantities. The elastic predictor‐plastic corrector algorithm is based on the solution of a system of 2 equations in 2 unknowns only. This further reduces to a single equation and unknown in the case of yield and plastic potential surfaces with a linear meridional section. The effectiveness of the proposed convexification technique and the efficiency and stability of the integration scheme are investigated by running numerical analyses of a notoriously demanding boundary value problem.     

Creep of saturated materials as a chemically enhanced rate‐dependent damage process

Material behaviour that exhibits characteristics of creep induced by a spontaneous mineral dissolution enhanced by material damage is studied. It is believed that the characteristic rates of the chemical processes involved determine the time‐rate dependence of the resulting strain. A basic model of a combined chemo‐plastic softening and chemically enhanced deviatoric strain hardening for saturated geomaterials is presented. Chemical softening is postulated to occur as a consequence of the net mass removal resulting from dissolution and precipitation of specific minerals occurring at the damage‐generated inter‐phase interfaces. Closed and open systems are discussed. In the former case, deformation at constant stress results entirely from a local compensation mechanism between the chemical softening and strain hardening. The classical three stages of creep are interpreted in terms of mechanisms of dissolution and precipitation, as well as the variation in the reaction surface areas involved in the mass exchange. In an open system, the above local mechanism is enhanced by the removal of mass via diffusion of species affecting the mass balance. Such a system is addressed via a boundary value problem as shown in an example. Copyright © 2007 John Wiley & Sons, Ltd.     

Stochastic micromechanical model of the deterioration of asphalt mixtures subject to moisture diffusion processes

The deleterious effect of moisture in the microstructure of asphalt mixtures, usually referred to as moisture damage, has been recognized as a main cause of early deterioration of asphalt pavements. The initiation and evolution of moisture‐related deterioration is strongly influenced by the internal air void structure of asphalt mixtures. Despite its importance, the majority of works conducted on the micromechanical modeling of asphalt mixtures overlook the role of the air void structure, mainly because of its high complexity and variability. This paper explores the influence of air void variability on the performance of asphalt mixtures subjected to moisture diffusion processes. A stochastic modeling technique based on random field theory was used to generate internal distributions of physical and mechanical properties of the asphalt matrix of the mixture that depend on probable air voids distributions. The analysis was conducted by means of a coupled numerical micromechanical model of moisture damage. The results showed that the variability and distribution of air voids are decisive in determining the moisture‐dependent performance of asphalt mixtures. Furthermore, it was also shown that a stochastic characterization of the diverse air void configurations is a feasible and useful approach to better represent and understand mechanically related deterioration processes in asphalt mixtures. Copyright © 2010 John Wiley & Sons, Ltd.     

A micromechanical modeling of ductile behavior of a porous chalk: Formulation,identification, and validation

Porosity strongly affects the overall ductile behavior of cohesive geomaterials undergoing plastic deformation. In the present paper, we proposed an original micromechanical approach that suitably couples Drucker–Prager‐type plasticity, evolving porosity under general triaxial loadings. The resulting model has the advantage to be based on a single macroscopic yield function, which also plays the role of plastic potential. It is shown that this yield function is particularly appropriate to account for the pore collapse and plastic shearing mechanisms that govern the mechanical behavior of the studied Lixhe chalk. Finally, the new model is implemented and validated by comparison to triaxial tests data, covering a wide range of confining pressures. Copyright © 2011 John Wiley & Sons, Ltd.     

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.     

Development of a Matrix‐Matched Sphalerite Reference Material (MUL‐ZnS‐1) for Calibration of In Situ Trace Element Measurements by Laser Ablation‐Inductively Coupled Plasma‐Mass Spectrometry

Sphalerite (ZnS) is an abundant ore mineral and an important carrier of elements such as Ge, Ga and In used in high‐technology applications. In situ measurements of trace elements in natural sphalerite samples using LA‐ICP‐MS are hampered by a lack of homogenous matrix‐matched sulfide reference materials available for calibration. The preparation of the MUL‐ZnS1 calibration material containing the trace elements V, Cr, Mn, Co, Ni, Cu, Ga, Ge, As, Se, Mo, Ag, Cd, In, Sn, Sb, Tl and Pb besides Zn, Fe and S is reported. Commercially available ZnS, FeS, CdS products were used as the major components, whereas the trace elements were added by doping with single‐element ICP‐MS standard solutions and natural mineral powders. The resulting powder mixture was pressed to pellets and sintered at 400 °C for 100 h using argon as an inert gas. To confirm the homogeneity of major and trace element distributions within the MUL‐ZnS1 calibration material, measurements were performed using EPMA, solution ICP‐MS, ICP‐OES and LA‐ICP‐MS. The results show that MUL‐ZnS‐1 is an appropriate material for calibrating trace element determination in sphalerite using LA‐ICP‐MS.     

Viscoelastic functionally graded finite element method with recursive time integration and applications to flexible pavements

The finite‐element (FE) method is used for modeling geotechnical and pavement structures exhibiting significant non‐homogeneity. Property gradients generated due to non‐homogeneous distributions of moisture is one such example for geotechnical materials. Aging and temperature‐induced property gradients are common sources of non‐homogeneity for asphalt pavements. Investigation of time‐dependent behavior combined with functionally graded property gradation can be accomplished by means of the non‐homogeneous viscoelastic analysis procedure. This paper describes the development of a generalized isoparametric FE formulation to capture property gradients within elements, and a recursive formulation for solution of hereditary integral equations. The formulation is verified by comparison with analytical and numerical solutions. Two application examples are presented: the first describes stationary crack‐tip fields for viscoelastic functionally graded materials, and the second example demonstrates the application of the proposed procedures for efficient and accurate simulations of interfaces between layers of flexible pavement. Copyright © 2011 John Wiley & Sons, Ltd.