A nonlocal model with strain-based damage

A thermodynamically consistent formulation of nonlocal damage in the framework of the internal variable theories of inelastic behaviours of associative type is presented. The damage behaviour is defined in the strain space and the effective stress turns out to be additively splitted in the actual stress and in the nonlocal counterpart of the relaxation stress related to damage phenomena. An important advantage of models with strain-based loading functions and explicit damage evolution laws is that the stress corresponding to a given strain can be evaluated directly without any need for solving a nonlinear system of equations. A mixed nonlocal variational formulation in the complete set of state variables is presented and is specialized to a mixed two-field variational formulation. Hence a finite element procedure for the analysis of the elastic model with nonlocal damage is established on the basis of the proposed two-field variational formulation. Two examples concerning a one-dimensional bar in simple tension and a two-dimensional notched plate are addressed. No mesh dependence or boundary effects are apparent.     

Spectral analysis of localization in nonlocal and over-nonlocal materials with softening plasticity or damage

The paper shows that spectral wave propagation analysis reveals in a simple and clear manner the effectiveness of various regularization techniques for softening materials, i.e., materials for which the yield limits soften as a function of the total strain. Both plasticity and damage models are considered. It is verified analytically in a simple way that the nonlocal integral-type model with degrading yield limit depending on the total strain works correctly if and only one adopts an unconventional nonlocal formulation introduced in 1994 by Vermeer and Brinkgreve (and in 1996 by Planas, and by Strömberg and Ristinmaa), which is here called, for the sake of brevity, ‘over-nonlocal’ because it uses a linear combination of local and nonlocal variables in which a negative weight imposed on the local variable is compensated by assigning to the nonlocal variable weight greater than 1 (this is equivalent to a nonlocal variable with a smooth positive weight function of total weight greater than 1, normalized by superposing a negative delta-function spike at the center). The spectral approach readily confirms that the nonlocal integral-type generalization of softening plasticity with an additive format gives correct localization properties only if an over-nonlocal formulation is adopted. By contrast, the nonlocal integral-type generalization of softening plasticity with a multiplicative format provides realistic localization behavior, just like the nonlocal integral-type damage model, and thus does not necessitate an over-nonlocal formulation. The localization behavior of explicit and implicit gradient-type models is also analyzed. A simple analysis shows that plasticity and damage models with gradient-type localization limiter, whether explicit or implicit, have very different localization behaviors.     

A symmetric nonlocal damage theory

The paper presents a thermodynamically consistent formulation for nonlocal damage models. Nonlocal models have been recognized as a theoretically clean and computationally efficient approach to overcome the shortcomings arising in continuum media with softening. The main features of the presented formulation are: (i) relations derived by the free energy potential fully complying with nonlocal thermodynamic principles; (ii) nonlocal integral operator which is self-adjoint at every point of the solid, including zones near to the solid’s boundary; (iii) capacity of regularizing the softening ill-posed continuum problem, restoring a meaningful nonlocal boundary value problem. In the present approach the nonlocal integral operator is applied consistently to the damage variable and to its thermodynamic conjugate force, i.e. nonlocality is restricted to internal variables only. The present model, when associative nonlocal damage flow rules are assumed, allows the derivation of the continuum tangent moduli tensor and the consistent tangent stiffness matrix which are symmetric. The formulation has been compared with other available nonlocal damage theories.Finally, the theory has been implemented in a finite element program and the numerical results obtained for 1-D and 2-D problems show its capability to reproduce in every circumstance a physical meaningful solution and fully mesh independent results.     

Nonlocal damage theory in hybrid-displacement formulations

This paper discusses three hybrid-displacement finite element formulations for the simulation of strain localization based on nonlocal damage theory. An isotropic integral nonlocal damage model is chosen. The hybrid finite element formulations adopted in this work are developed from first principles of Mechanics. The first one defines the domain approximations using the Trefftz functions derived for the linear elastic regime. When damage appears the hybrid-Trefftz displacement formulation degenerates into an hybrid-displacement formulation. The second formulation uses an enriched Trefftz basis with the consideration of local Heaviside functions. The third formulation uses orthonormal Legendre polynomials for the domain approximations. A set of benchmark tests is presented and discussed in order to compare the performance and accuracy of the different models. It is shown that the proposed hybrid-Trefftz formulation allows the reproduction of the general behavior of the structure but does not lead to a correct simulation of the strain tensor evolution. The hybrid-displacement formulation that uses orthonormal Legendre polynomials gives coherent results, so it appears to be a promising field of investigation.     

A damage model with evolving nonlocal interactions

We present a damage model for softening materials with evolving nonlocal interactions. The thermodynamic implications and the material stability issue are addressed. The proposed nonlocal averaging scheme provides the obtained constitutive models with an evolving nonlocal interaction which is activated only when damage occurs. In the analysis of structures made of quasi-brittle materials, this feature helps not only to overcome some issues with the incorrect initiation of damage but also to better control the evolving size of the active fracture process zone. This is an essential feature that is usually not considered in depth in many existing nonlocal approaches to the continuum modelling of quasi-brittle fracture. Numerical examples are given to demonstrate features of the proposed modelling approach.     

Novel variational formulations for nonlocal plasticity

A nonlocal structural model of softening plasticity is considered in the framework of the internal variable theories of inelastic behaviours of associative type. The finite-step nonlocal structural problem in a geometrically linear range is formulated according to a backward difference scheme for time integration of the flow rule. The related finite-step variational formulation in the complete set of local and nonlocal state variables is recovered. A family of mixed nonlocal variational formulations, with different combinations of state variables, is provided starting from the general variational formulation. The specialization of a mixed variational formulation to existing nonlocal models of softening plasticity, assuming both linear and nonlinear constitutive behaviour, is provided to show the effectiveness of the theory.     

Stress-based nonlocal damage model

Progressive microcracking in brittle or quasi-brittle materials, as described by damage models, presents a softening behavior that in turn requires the use of regularization methods in order to maintain objective results. Such regularization methods, which describe interactions between points, provide some general properties (including objectivity and the non-alteration of a uniform field) as well as drawbacks (damage initiation, free boundary).A modification of the nonlocal integral regularization method that takes the stress state into account is proposed in this contribution. The orientation and intensity of nonlocal interactions are modified in accordance with the stress state. The fundamental framework of the original nonlocal method has been retained, making it possible to maintain the method’s advantages. The modification is introduced through the weight function, which in this modified version depends not only on the distance between two points (as for the original model) but also on the stress state at the remote point.The efficiency of this novel approach is illustrated using several examples. The proposed modification improves the numerical solution of problems observed in numerical simulations involving regularization techniques. Damage initiation and propagation in mode I as well as shear band formation are analyzed herein.     

Micromorphic approach for finite gradient-elastoplasticity fully coupled with ductile damage: Formulation and computational aspects

It is well established that the use of inelastic constitutive equations accounting for induced softening, leads to pathological space (mesh) and time discretization dependency of the numerical solution of the associated Initial and Boundary Value Problem (IBVP). To avoid this drawback, many less or more approximate solutions have been proposed in the literature in order to regularize the IBVP and to obtain numerical solutions which are, at convergence, much less sensitive to the space and the time discretization. The basic idea behind these regularization techniques is the formulation of nonlocal constitutive equations by introducing some effects of characteristic lengths representing the materials microstructure. In this work, using the framework of generalized nonlocal continua, a thermodynamically-consistent micromorphic formulation using appropriate micromorphic state variables and their first gradients, is proposed in order to extend the classical local constitutive equations by incorporating appropriate characteristic internal lengths. The isotropic damage, the isotropic and the kinematic hardenings are supposed to carry the targeted micromorphic effects. First the theoretical aspects of this fully coupled micromorphic formulation is presented in details and the proposed generalized balance equations as well as the fully coupled micromorphic constitutive equations deduced. The associated numerical aspects in the framework of the classical Galerkin-based FE formulation are briefly discussed in the special case of micromorphic damage. Specifically, the formulation of 2D finite elements with additional degrees of freedom (d.o.f.), the dynamic explicit global resolution scheme as well as the local integration scheme, to compute the stress tensor and the state variables at each integration point of each element, are presented. Application is made to the typical uniaxial tension specimen under plane strain conditions in order to chow the predictive capabilities of the proposed micromorphic model, particularly against its ability to give (at convergence) a mesh independent solution even for high values of the ductile damage (i.e., the macroscopic cracks).     

On critical assessment of the use of local and nonlocal damage models for prediction of ductile crack growth and crack path in various loading and boundary conditions

In this work, we have assessed the results of the local and nonlocal versions of Rousselier’s damage model, which have been used here for simulation of ductile crack growth. There are several issues regarding the accuracy of the results which has been addressed in this paper, e.g., accuracy in simulation of crack path, extent and width of the damaged region, fracture resistance behaviour in situations such as symmetric vs. non-symmetric boundary-value problems, mixed-mode loading vs. mode-I loading of the crack-tip, etc. It was also observed that the shape and orientation of the elements at the crack-tip, in addition to their size, influence the results of the local damage model. In this work, it was shown that the above issues can be resolved through the use of nonlocal damage models. The predictions of the nonlocal model are also consistent with the experimental observations unlike its local counterpart. Several examples were presented, where the results as obtained by both the local and nonlocal models were compared. From this experience, it is recommended that the local damage models should not be used blindly by the analysts for all kinds of mesh design, loading, boundary conditions, etc.     

A model for high temperature creep of single crystal superalloys based on nonlocal damage and viscoplastic material behavior

A model for high temperature creep of single crystal superalloys is developed, which includes constitutive laws for nonlocal damage and viscoplasticity. It is based on a variational formulation, employing potentials for free energy, and dissipation originating from plasticity and damage. Evolution equations for plastic strain and damage variables are derived from the well-established minimum principle for the dissipation potential. The model is capable of describing the different stages of creep in a unified way. Plastic deformation in superalloys incorporates the evolution of dislocation densities of the different phases present. It results in a time dependence of the creep rate in primary and secondary creep. Tertiary creep is taken into account by introducing local and nonlocal damage. Herein, the nonlocal one is included in order to model strain localization as well as to remove mesh dependence of finite element calculations. Numerical results and comparisons with experimental data of the single crystal superalloy LEK94 are shown.     

Nonlocal implicit gradient-enhanced elasto-plasticity for the modelling of softening behaviour

An improved gradient-enhanced approach for softening elasto-plasticity is proposed, which in essence is fully nonlocal, i.e. an equivalent integral nonlocal format exists. The method utilises a nonlocal field variable in its constitutive framework, but in contrast to the integral models computes this nonlocal field with a gradient formulation. This formulation is considered ‘implicit’ in the sense that it strictly incorporates the higher-order gradients of the local field variable indirectly, unlike the common (explicit) gradient approaches. Furthermore, this implicit gradient formulation constitutes an additional partial differential equation (PDE) of the Helmholtz type, which is solved in a coupled fashion with the standard equilibrium condition. Such an approach is particularly advantageous since it combines the long-range interactions of an integral (nonlocal) model with the computational efficiency of a gradient formulation. Although these implicit gradient approaches have been successfully applied within damage mechanics, e.g. for quasi-brittle materials, the first attempts were deficient for plasticity. On the basis of a thorough comparison of the gradient-enhancements for plasticity and damage this paper rephrases the problem, which leads to a formulation that overcomes most reported problems. The two-dimensional finite element implementation for geometrically linear plain strain problems is presented. One- and two-dimensional numerical examples demonstrate the ability of this method to numerically model irreversible deformations, accompanied by the intense localisation of deformation and softening up to complete failure.     

Over-nonlocal gradient enhanced plastic-damage model for concrete

Classical continuum models exhibit strong mesh dependency during softening. One method to regularize the problem is to introduce a length scale parameter via the nonlocal formulation. However, standard nonlocal enhancement (either by integral or gradient formulation) may serve only as a partial localization limiter for many material models. The “over-nonlocal” formulation, where the weight for the nonlocal value is greater than unity and the excesses compensated by assigning a negative weight to the local value, is able to fully regularize certain material models when standard nonlocal enhancement fails to do so. A plastic-damage model for concrete is formulated with this over-nonlocal enhancement via the gradient approach and the full regularizing capabilities demonstrated.     

Variational formulations,convergence and stability properties in nonlocal elastoplasticity

A thermodynamically consistent formulation of nonlocal plasticity in the framework of the internal variable theories of inelastic behaviors of associative type is presented. A family of mixed variational formulations, with different combinations of state variables, is provided starting from the finite-step nonlocal elastoplastic structural problem. It is shown that a suitable minimum principles provides a rational basis to exploit the iterative elastic predictor-plastic corrector algorithm in terms of the dissipation functional. A sufficient condition is proved for the convergence of the iterative elastic predictor-plastic corrector algorithm based on a suitable choice of the elastic operator in the prediction phase and a necessary and sufficient condition for the existence of a unique solution (if any) of the nonlocal problem at hand is then provided. The nonlinear stability analysis of the nonlocal problem is carried out following the concept of nonexpansivity proposed in local plasticity.     

Visco-potential free-surface flows and long wave modelling

In a previous study [D. Dutykh, F. Dias, Viscous potential free-surface flows in a fluid layer of finite depth, C. R. Acad. Sci. Paris, Ser. I 345 (2007) 113–118] we presented a novel visco-potential free-surface flows formulation. The governing equations contain local and nonlocal dissipative terms. From physical point of view, local dissipation terms come from molecular viscosity but in practical computations, rather eddy viscosity should be used. On the other hand, nonlocal dissipative term represents a correction due to the presence of a bottom boundary layer. Using the standard procedure of Boussinesq equations derivation, we come to nonlocal long wave equations. In this article we analyze dispersion relation properties of proposed models. The effect of nonlocal term on solitary and linear progressive waves attenuation is investigated. Finally, we present some computations with viscous Boussinesq equations solved by a Fourier type spectral method.     

Peridynamic modelling of impact damage in three-point bending beam with offset notch

The nonlocal peridynamic theory has been proven to be a promising method for the material failure and damage analyses in solid mechanics.Based upon the integrodifferential equations,peridynamics enables predicting the complex fracture phenomena such as spontaneous crack nucleation and crack branching,curving,and arrest.In this paper,the bond-based peridynamic approach is used to study the impact damage in a beam with an offset notch,which is widely used to investigate the mixed I-II crack propagation in brittle materials.The predictions from the peridynamic analysis agree well with available experimental observations.The numerical results show that the dynamic fracture behaviors of the beam under the impact load,such as crack initiation,curving,and branching,rely on the location of the offset notch and the impact speed of the drop hammer.     

Nonlocal microplane model with strain-softening yield limits

The paper deals with the problem of nonlocal generalization of constitutive models such as microplane model M4 for concrete, in which the yield limits, called stress–strain boundaries, are softening functions of the total strain. Such constitutive models call for a different nonlocal generalization than those for continuum damage mechanics, in which the total strain is reversible, or for plasticity, in which there is no memory of the initial state. In the proposed nonlocal formulation, the softening yield limit is a function of the spatially averaged nonlocal strains rather than the local strains, while the elastic strains are local. It is demonstrated analytically as well numerically that, with the proposed nonlocal model, the tensile stress across the strain localization band at very large strain does soften to zero and the cracking band retains a finite width even at very large tensile strain across the band only if one adopts an “over-nonlocal” generalization of the type proposed by Vermeer and Brinkgreve [In: Chambon, R., Desrues, J., Vardoulakis, I. (Eds.), Localisation and Bifurcation Theory for Soils and Rocks, Balkema, Rotterdam, 1994, p. 89] (and also used by Planas et al. [Basic issue of nonlocal models: uniaxial modeling, Tecnical Report 96-jp03, Departamento de Ciencia de Materiales, Universidad Politecnica de Madrid, Madrid, Spain, 1996], and by Strömberg and Ristinmaa [Comput. Meth. Appl. Mech. Eng. 136 (1996) 127]). Numerical finite element studies document the avoidance of spurious mesh sensitivity and mesh orientation bias, and demonstrate objectivity and size effect.     

Localisation issues in local and nonlocal continuum approaches to fracture

Continuum approaches to fracture regard crack initiation and growth as the ultimate consequences of a gradual, local loss of material integrity. The material models which are traditionally used to describe the degradation process, however, may predict premature crack initiation and instantaneous, perfectly brittle crack growth. This nonphysical response is caused by localisation instabilities due to loss of ellipticity of the governing equations and—more importantly—singularity of the damage rate at the crack tip. It is argued that this singularity results in instantaneous failure in a vanishing volume, even if ellipticity is not first lost. Adding strong nonlocality to the modelling is shown to preclude localisation instabilities and remove damage rate singularities. As a result, premature crack initiation is avoided and crack growth rates remain finite. Weak nonlocality, as provided by explicit gradient models, does not suffice for this purpose. In implementing the enhanced modelling, the crack must be excluded from the equilibrium problem and the nonlocal interactions in order to avoid unrealistic damage growth.     

Nonlocal continuum theory of anisotropically damaged metals

The paper deals with a consistent and systematic general framework for the development of anisotropic continuum damage in ductile metals based on thermodynamic laws and nonlocal theories. The proposed model relies on finite strain kinematics based on the consideration of damaged as well as fictitious undamaged configurations related via metric transformation tensors which allow for the interpretation of damage tensors. The formulation is accomplished by rate-independent plasticity using a nonlocal yield condition of Drucker–Prager type, anisotropic damage based on a nonlocal damage growth criterion as well as non-associated flow and damage rules. The nonlocal theory of inelastic continua is established to be able to take into account long-range microstructural interaction. The approach incorporates macroscopic interstate variables and their higher-order gradients which properly describe the change in the internal structure and investigate the size effect of statistical inhomogeneity of the heterogeneous material. The idea of bridging length-scales is made by using higher-order gradients in the evolution equations of the equivalent inelastic strain measures which leads to a system of elliptic partial differential equations which is solved using the finite difference method at each iteration of the loading step and the displacement-based finite element procedure is governed by the standard principle of virtual work. Numerical simulations of the elastic–plastic deformation behavior of damaged solids demonstrate the efficiency of the formulation. Tension tests undergoing large strains are used to investigate the damage growth in high strength steel. The influence of various model parameters on the prediction of the deformation and localization of ductile metals is discussed.     

One-dimensional nonlocal and gradient elasticity: Closed-form solution and size effect

The equivalence between nonlocal and gradient elasticity models is investigated by making reference to one-dimensional boundary value problems equipped with two integral stress–strain laws proposed by Eringen (Nonlocal Continuum Field Theories (2002)). Corresponding closed-form solutions are derived through a procedure for the reduction of integral to differential equations. The reproduction of size effects in micro/nano rods is discussed. The differential formulation associated with the local/nonlocal model is shown to correspond to the strain-gradient formulation proposed by Aifantis (Mech. Mater. 35 (2003) 259–280).     

Generalization of cumulative damage criterion to multilevel cyclic loading

A cumulative damage theory is presented for multilevel cyclic loading. Both the crack initiation and propagation stages are included in this method. The criterion is based on the concept of the total strain energy density as a damage parameter. The effect of the mean stress/strain on the damage accumulation is also included in this formulation. It is shown that a number of earlier proposed models can be derived as a particular case of the present criterion. The predicted results are compared with the experimental data of two, three, five and six steps cyclic loading histories. A good agreement is noted between the predicted and experimental results.