A. RODRÍGUEZ-FERRAN, I. MORATA AND A. HUERTA quantities. In fact, a number of proposals can be found in the literature. For integral-type regularization, some examples are the use of a nonlocal damage parameter [3], nonlocal strains [7] or nonlocal strain invariants [8]. These and other existing approaches are compared in [9] by means of a simple 1D numerical test (bar under uniaxial tension). Various approaches are also possible for gradient regularization. In [10], for instance, the loading function depends on the Laplacian of the damage parameter. Note that all these strategies involve Gauss-point based quantities.A new proposal is made here: to use nonlocal displacements to regularize the problem. The two versions are proposed, discussed and compared: integral-type (nonlocal displacements obtained as the weighted average of standard, local displacements), see [11], and gradienttype (nonlocal displacements obtained as the solution of a second-order PDE). As discussed and illustrated by means of numerical examples, the regularization capabilities of this new model are very similar to that of the standard model. In addition, it is very attractive from a computational viewpoint, especially regarding (1) the computation of the consistent tangent matrix and (2) the simple and straightforward upgrade of a nonlinear FE code to account for nonlocality.An outline of this paper follows. The basic features of standard nonlocal damage models are reviewed in Section 2. The new model based on nonlocal displacements is presented in Section 3. The integral-type and gradient versions are discussed in Sections 3.1 and 3.2 respectively. The regularization capabilities are illustrated by means of a uniaxial tension test. Section 4 deals with the consistent linearization of the nonlinear equilibrium equation. It is shown how the consistent tangent matrix is much simpler to compute for the new model than for the standard models, both in the integral-type (Section 4.1) and the gradient (Section 4.2) cases. Quadratic convergence is shown for the uniaxial tension test. Section 5 shows how nonlocal displacements can be used to incorporate nonlocality into a FE code in a very simple and efficient manner, especially if the gradient regularization is chosen. The concluding remarks of Section 6 close the paper.Standard notation is used. Vector fields in the continuum are represented by slanted boldface type (u u: displacement field). Nodal vectors associated to FE discretization are denoted by upright boldface type (u: nodal displacements).
OVERVIEW OF DAMAGE MODELSFor simplicity, only elastic-scalar damage models are considered here. However, the concept of nonlocal displacements can be extended to more complex damage models exhibiting, for instance, anisotropy or plasticity [4,12].
Local damage modelsA generic local damage model consists of the following equations, summarized in table I:• A relation between Cauchy stresses σ σ and small strains ε ε -i.e. the symmetrized gradient of displacements u u, Equation (2)-, where the loss of sti...