In its numerical implementation, the variational approach to brittle fracture approximates the crack evolution in an elastic solid through the use of gradient damage models. In the present paper, we first formulate the quasi-static evolution problem for a general class of such damage models. Then, we introduce a stability criterion in terms of the positivity of the second derivative of the total energy under the unilateral constraint induced by the irreversibility of damage. These concepts are applied in the particular setting of a one-dimensional traction test. That allows us to construct homogeneous as well as localized damage solutions in a closed form and to illustrate the concepts of loss of stability, of scale effects, of damage localization, and of structural failure. Considering several specific constitutive models, stress-displacement curves, stability diagrams, and energy dissipation provide identification criteria for the relevant material parameters, such as limit stress and internal length. Finally, the one-dimensional analytical results are compared with the numerical solution of the evolution problem in a two dimensional setting.

a b s t r a c tWe consider a wide class of gradient damage models which are characterized by two constitutive functions after a normalization of the scalar damage parameter. The evolution problem is formulated following a variational approach based on the principles of irreversibility, stability and energy balance. Applied to a monotonically increasing traction test of a one-dimensional bar, we consider the homogeneous response where both the strain and the damage fields are uniform in space. In the case of a softening behavior, we show that the homogeneous state of the bar at a given time is stable provided that the length of the bar is less than a state dependent critical value and unstable otherwise. However, we also show that bifurcations can appear even if the homogeneous state is stable. All these results are obtained in a closed form. Finally, we propose a practical method to identify the two constitutive functions. This method is based on the measure of the homogeneous response in a situation where this response is stable without possibility of bifurcation, and on a procedure which gives the opportunity to detect its loss of stability. All the theoretical analyses are illustrated by examples.

The paper is devoted to gradient damage models which allow us to describe all the process of degradation of a body including the nucleation of cracks and their propagation. The construction of such model follows the variational approach to fracture and proceeds into two stages: (i) definition of the energy; (ii) formulation of the damage evolution problem. The total energy of the body is defined in terms of the state variables which are the displacement field and the damage field in the case of quasi-brittle materials. That energy contains in particular gradient damage terms in order to avoid too strong damage localizations. The formulation of the damage evolution problem is then based on the concepts of irreversibility, stability and energy balance. That allows us to construct homogeneous as well as localized damage solutions in a closed form and to illustrate the concepts of loss of stability, of scale effects, of damage localization, and of structural failure. Moreover, the variational formulation leads to a natural numerical method based on an alternate minimization algorithm. Several numerical examples illustrate the ability of this approach to account for all the process of fracture including a 3D thermal shock problem where the crack evolution is very complex.