A micromechanics-based elastoplastic damage model for granular materials at low confining pressure

This paper is devoted to the formulation of a micromechanics-based constitutive model for granular materials under relatively low confining pressure. The constitutive formulation is performed within the general framework of homogenization for granular materials. However, new rigorous stress localization laws are proposed. Some local constitutive relations are established under the consideration of irreversible thermodynamics. Macroscopic plastic deformation is obtained by considering local plastic sliding in a limit number of families of contact planes. The plastic sliding at each contact plane is described by a non-associated plastic flow rule, taking into account pressure sensitivity and normal dilatancy. Nonlinear elastic deformation related to progressive compaction of contacts is also taken into account. Material softening is described by involving damage process related to degradation of microstructure fabric. The proposed model is applied to some typical granular materials (sands). The numerical predictions are compared with experimental data.     

A three-dimensional progressive damage model for fibre-composite materials

This note presents a damage model for fibre-composite materials based in the approach by Matzenmiller et al. [Matzenmiller, A., Lubliner, J., Taylor, R.L., 1995. A constitutive model for anisotropic damage in fiber-composites. Mech. Mater. 20, 125]. In this work, the model is developed in a three-dimensional context with modified formulation for the constitutive law and damage evolution. An orthotropic composite subjected to mixed failure modes is assumed in this development. Its formulation and implementation details are provided.     

A penny-shaped cohesive crack model for material damage

A novel micromechanics based damage model is proposed to address failure mechanism of defected solids with randomly distributed penny-shaped cohesive micro-cracks (Barenblatt–Dugdale type). Energy release contribution to the material damage process is estimated in a representative volume element (RVE) under macro hydrostatic stress state. Macro-constitutive relations of RVE are derived via self-consistent homogenization scheme, and they are characterized by effective nonlinear elastic properties and a class of pressure sensitive plasticity which depends on crack opening volume fraction and Poisson’s ratio. Several distinguished features of the present model are compared with Gurson model and Gurson–Tvergaard–Needleman (GTN) model, showing that the proposed model can better capture material degradation and catastrophic failure due to cohesive micro-crack growth and coalescence.     

A thermodynamic framework to model thixotropic materials

Thixotropic materials are widely used in a variety of industrial applications. The constitutive relations to describe these materials are based on one-dimensional experiments in which the material is subjected to a shear motion and there is no unique methodology to obtain proper three-dimensional models. The path towards generalization to a three-dimensional framework is invariably carried out in a ad hoc manner. Here we propose a three-dimensional model that stems from a general thermodynamic framework that has proved to be quite robust in the development of constitutive relations, namely the application of the second law of thermodynamics together with the maximization of the entropy production. This leads to a constitutive equation that has the same form of a generalized Upper Convected Maxwell equation, if we require that changes of microstructure due to the deformation of each Maxwell element that comprises the model are reversible. Changes in microstructure are governed by a potential that is a measure of the difference between the current structure and the equilibrium structure associated with it. The equilibrium structure associated with the current structure is determined by the current value of stress, considered the main break up agent. We assume that the state of equilibrium would be achieved in a Motion With Constant Stress History, starting from the current stress state, until a steady state where the kinematics is not changing.     

Investigation of the damage variable basic issues in continuum damage and healing mechanics

Consistent mathematics and mechanics are used here to properly interpret the damage variable within the confines of the concept of reduced area due to damage. In this work basic issues are investigated for the damage variable in conjunction with continuum damage and healing mechanics. First, the issue of the additive decomposition of the damage variable into damage due to voids and damage due to cracks in continuum damage mechanics is discussed. The accurate decomposition is shown to be non-additive and involves a term due to the interaction of cracks and voids. It is shown also that the additive decomposition can only be used for the special case of small damage. Furthermore, a new decomposition is derived for the evolution of the damage variable. The second issue to be discussed is the new concept of independent and dependent damage processes. For this purpose, exact expressions for the two types of damage processes are presented. The third issue addressed is the concept of healing processes occurring in series and in parallel. In this regard, systematically and consistently, the equations of healing processes occurring either consecutively or simultaneously are discussed. This is followed by introducing the new concept of small healing in damaged materials. Simplified equations that apply when healing effects are small are shown. Finally, some interesting and special damage processes using a systematic and original formulation are presented.     

A thermodynamic model for electrical current induced damage

Electromigration-induced damage, which is in principal an irreversible mass diffusion under high current density, has been a concern for VLSI design for a long time. Miniaturization of electronic device sizes down to nano-scale will make electromigration a concern for all conducting components. This paper uses thermodynamics, statistical mechanics and mass transport (diffusion) principals to propose a model for electromigration process and a damage evolution model to quantify the degradation in microelectronics (and micro electromechanical system) solder joints subjected to high current densities. Entropy production in the system is used as a damage metric. The irreversible thermodynamic damage model utilized in this work has previously been successfully applied to thermo-mechanical fatigue of microelectronic solder joints. In this paper we extend this model to electromigration-induced degradation.Electromigration process is modeled by the atomic vacancy flux (mass diffusion) process. The proposed unified model is compared with several existing analytical and empirical models. A comparison of the damage evolution model proposed in here agrees well with empirical models proposed in the literature.     

Relations between a micro-mechanical model and a damage model for ductile failure in shear

Gurson type constitutive models that account for void growth to coalescence are not able to describe ductile fracture in simple shear, where there is no hydrostatic tension in the material. But recent micro-mechanical studies have shown that in shear the voids are flattened out to micro-cracks, which rotate and elongate until interaction with neighbouring micro-cracks gives coalescence. Thus, the failure mechanism is very different from that under tensile loading. Also, the Gurson model has recently been extended to describe failure in shear, by adding a damage term to the expression for the growth of the void volume fraction, and it has been shown that this extended model can represent experimental observations. Here, numerical studies are carried out to compare predictions of the shear-extended Gurson model with the shear failures predicted by the micro-mechanical cell model. Both models show a strong dependence on the level of hydrostatic tension. Even though the reason for this pressure dependence is different in the two models, as the shear-extended Gurson model does not describe voids flattening out and the associated failure mechanism by micro-cracks interacting with neighbouring micro-cracks, it is shown that the trends of the predictions are in good agreement.     

A continuum damage model for piezoelectric materials

In this paper, a constitutive model is proposed for piezoelectric material solids containing distributed cracks. The model is formulated in a framework of continuum damage mechanics using second rank tensors as internal variables. The Helrnhotlz free energy of piezoelectric mate- rials with damage is then expressed as a polynomial including the transformed strains, the electric field vector and the tensorial damage variables by using the integrity bases restricted by the initial orthotropic symmetry of the material. By using the Talreja＇s tensor valued internal state damage variables as well as the Helrnhotlz free energy of the piezoelectric material, the constitutive relations of piezoelectric materials with damage are derived. The model is applied to a special case of piezoelectric plate with transverse matrix cracks. With the Kirchhoff hypothesis of plate, the free vibration equations of the piezoelectric rectangular plate considering damage is established. By using Galerkin method, the equations are solved. Numerical results show the effect of the damage on the free vibration of the piezoelectric plate under the close-circuit condition, and the present results are compared with those of the three-dimensional theory.     

Three-dimensional constitutive model for shape memory alloys based on microplane model

A new model for the behavior of polycrystalline shape memory alloys (SMA), based on a statically constrained microplane theory, is proposed. The new model can predict three-dimensional response by superposing the effects of inelastic deformations computed on several planes of different orientation, thus reproducing closely the actual physical behavior of the material. Due to the structure of the microplane algorithm, only a one-dimensional constitutive law is necessary on each plane. In this paper, a simple constitutive law and a robust kinetic expression are used as the local constitutive law on the microplane level. The results for SMA response on the macroscale are promising: simple one-dimensional response is easily reproduced, as are more complex features such as stress-strain subloops and tension-compression asymmetry. A key feature of the new model is its ability to accurately represent the deviation from normality exhibited by SMAs under nonproportional loading paths.     

Analytical and numerical analysis of frictional damage in quasi brittle materials

Frictional sliding and crack growth are two main dissipation processes in quasi brittle materials. The frictional sliding along closed cracks is the origin of macroscopic plastic deformation while the crack growth induces a material damage. The main difficulty of modeling is to consider the inherent coupling between these two processes. Various models and associated numerical algorithms have been proposed. But there are so far no analytical solutions even for simple loading paths for the validation of such algorithms. In this paper, we first present a micro-mechanical model taking into account the damage-friction coupling for a large class of quasi brittle materials. The model is formulated by combining a linear homogenization procedure with the Mori–Tanaka scheme and the irreversible thermodynamics framework. As an original contribution, a series of analytical solutions of stress–strain relations are developed for various loading paths. Based on the micro-mechanical model, two numerical integration algorithms are exploited. The first one involves a coupled friction/damage correction scheme, which is consistent with the coupling nature of the constitutive model. The second one contains a friction/damage decoupling scheme with two consecutive steps: the friction correction followed by the damage correction. With the analytical solutions as reference results, the two algorithms are assessed through a series of numerical tests. It is found that the decoupling correction scheme is efficient to guarantee a systematic numerical convergence.     

Coupling between martensitic phase transformations and plasticity: A microstructure-based finite element model

A microstructural finite element (MFE) model is developed to capture the interaction between martensitic transformations and plasticity in NiTi shape memory alloys (SMAs). The interaction is modeled through the grain-to-grain redistribution of stress caused by both plasticity and phase transformation, so that each mechanism affects the driving force of the other. A unique feature is that both processes are modeled at a crystallographic level and are allowed to operate simultaneously. The model is calibrated to pseudoelastic data for select single crystals of Ti–50.9at.%Ni. For polycrystals, plasticity is predicted to enhance the overall martensite volume fraction at a given applied stress. Upon unloading, residual stress can induce remnant (retained) martensite. For thermal cycling under load bias, plasticity is observed to limit the net transformation strain/cycle and increase the hysteretic width. Deformation processing, via plastic pre-straining at elevated temperature, is shown to dramatically alter subsequent pseudoelastic response, as well as induce two-way shape memory behavior during no-load thermal cycling. Overall, the model is suitable at smaller imposed strains, where martensite detwinning is not expected to dominate.     

A damage model for the blistering of polyvinylidene fluoride subjected to carbon dioxide decompression

Semicrystalline polymers (SCPs) used in complex petroleum structures are subjected to high variations of temperature and gas pressure, which induce some damage. For this reason, the Representative Volume Element (RVE) of a SCP evolving in a gaseous environment is modelled using a phenomenological approach based on the porous media theory. SCPs are in fact two-phase materials at the scale of the spherolite. One phase is crystalline (skeleton), and the other corresponds to a mixture (fluid) of gas and free amorphous medium (the latter is considered penetrable by gas). The modelling is described within the framework of the thermodynamics of irreversible processes with internal variables by considering the above RVE. An elastoviscoplastic model with cavity growth is proposed within the framework of general diffuso-mechanical continuum media. The aim is to predict the evolution of ductile damage observed during decompression. This two-phase diffuso-elastoviscoplastic-damage model is implemented in Abaqus^® software via a user subroutine. This numerical tool allows us to study the damage of polyvinylidene fluoride SCP during a rapid decompression test.The influence of the decompression rate on the evolution of initial cavities is discussed. Comparison with experiments shows the relevance of the proposed model as far as blistering is concerned. We hope that the present model can help understanding the occurrence of damage mechanisms during decompression.     

Interaction-based non-local damage model for failure in quasi-brittle materials

The purpose of this paper is to present a new macroscopic approach to describe the evolving non-local interactions which are produced at the mesoscale during damage and failure in quasi-brittle materials. A new-integral type non-local model is provided where the weight function is directly built from these interactions, and therefore takes into account their evolution during the material failure intrinsically.     

A porosity-dependent inelastic criterion for engineering materials

Many criteria have been developed to describe the yielding condition, plastic potential, and failure strength of engineering materials. In this paper, the authors first review the characteristics of a few of the more common criteria used for porous materials. It is then shown that the main features of many criteria can be represented by a unique set of recently developed equations. The ensuing multiaxial criterion becomes applicable to a variety of materials and loading states. One of the advantages of the proposed criterion, named MSDP_u, is that it is explicitly porosity-dependent. The validity of this general inelastic criterion is demonstrated using experimental results obtained from various types of materials. A brief discussion follows on the advantages and limitations of the proposed equations.     

An isotropic unilateral damage model coupled with frictional sliding for quasi-brittle materials

In this paper, we present an original extension of an isotropic damage model for quasi-brittle materials and assess its predictive capabilities. The proposed model accounts not only for unilateral behavior related to the opening and closure of microcracks but also for inelastic strains reflecting the frictional sliding along closed microcracks. More importantly, owing to its careful mathematical formulation with a particular attention paid to the continuous differentiability of the underlined thermodynamic potential, the model ensures the continuity of the inelastic stress–strain response. First applications show that it is able to predict the asymmetric behavior and hysteretic response of microcracked materials such as concrete and some rocks     

Three-dimensional phenomenological thermodynamic model of pseudoelasticity of shape memory alloys at finite strains

The familiar small strain thermodynamic 3D theory of isotropic pseudoelasticity proposed by Raniecki and Lexcellent is generalized to account for geometrical effects. The Mandel concept of mobile isoclinic, natural reference configurations is used in order to accomplish multiplicative decomposition of total deformation gradient into elastic and phase transformation (p.t.) parts, and resulting from it the additive decomposition of Eulerian strain rate tensor. The hypoelastic rate relations of elasticity involving elastic strain rate are derived consistent with hyperelastic relations resulting from free energy potential. It is shown that use of Jaumann corotational rate of stress tensor in rate constitutive equations formulation proves to be convenient. The formal equation for p.t. strain rate , describing p.t. deformation effects is proposed, based on experimental evidence. Phase transformation kinetics relations are presented in objective form. The field, coupled problem of thermomechanics is specified in rate weak form (rate principle of virtual work, and rate principle of heat transport). It is shown how information on the material behavior and motion inseparably enters the rate virtual work principle through the familiar bridging equation involving Eulerian rate of nominal stress tensor.

---

     

Models for one-variant shape memory materials based on dissipation functions

Some simple models for the macroscopic behavior of shape memory materials whose microstructure can be described as a mixture of two phases are derived on the basis of a free energy and a dissipation function. Keeping a common expression for the free energy, each model is based on a different expression for the dissipation function. Temperature-induced as well as isothermal, adiabatic and convective stress-induced transformations are studied. Attention is paid to closed form solutions, comparison among the models and parameter identification.     

A constitutive theory for creep behavior of initially isotropic materials sustaining unilateral damage

The goal of the present work is to modify structure of the creep constitutive equations existing in the literature, and simultaneously to incorporate both damage induced anisotropy and unilateral damage into the constitutive model. The proposed nonlinear-tensor constitutive equation for creep together with the damage evolution equation take into account the secondary and tertiary creep of the initially isotropic materials. The material parameters of the model are determined using basic experiments. It is shown that the creep model is capable of describing available experimental data for the lateral creep responses under uniaxial compression.     

A continuum damage model for delaminations in laminated composites

Delamination, a typical mode of interfacial damage in laminated composites, has been considered in the context of continuum damage mechanics in this paper. Interfaces where delaminations could occur are introduced between the constituent layers. A simple but appropriate continuum damage representation is proposed. A single scalar damage parameter is employed and the degradation of the interface stiffness is established. Use has been made of the concept of a damage surface to derive the damage evolution law. The damage surface is constructed so that it combines the conventional stress-based and fracture-mechanics-based failure criteria which take account of mode interaction in mixed-mode delamination problems. The damage surface shrinks as damage develops and leads to a softening interfacial constitutive law. By adjusting the shrinkage rate of the damage surface, various interfacial constitutive laws found in the literature can be reproduced. An incremental interfacial constitutive law is also derived for use in damage analysis of laminated composites, which is a non-linear problem in nature. Numerical predictions for problems involving a DCB specimen under pure mode I delamination and mixed-mode delamination in a split beam are in good agreement with available experimental data or analytical solutions. The model has also been applied to the prediction of the failure strength of overlap ply-blocking specimens. The results have been compared with available experimental and alternative theoretical ones and discussed fully.     

A micromechanical damage and fracture model for polymers based on fractional strain-gradient elasticity

We formulate a simple one-parameter macroscopic model of distributed damage and fracture of polymers that is amenable to a straightforward and efficient numerical implementation. We show that the macroscopic model can be rigorously derived, in the sense of optimal scaling, from a micromechanical model of chain elasticity and failure regularized by means of fractional strain-gradient elasticity. In particular, we derive optimal scaling laws that supply a link between the single parameter of the macroscopic model, namely, the critical energy-release rate of the material, and micromechanical parameters pertaining to the elasticity and strength of the polymer chains and to the strain-gradient elasticity regularization. We show how the critical energy-release rate of specific materials can be determined from test data. Finally, we demonstrate the scope and fidelity of the model by means of an example of application, namely, Taylor-impact experiments of polyurea 1000 rods.