On Topological Derivatives for Contact Problems in Elasticity

In this article, a general method for shape-topology sensitivity analysis of contact problems is proposed. The method uses domain decomposition combined with specific properties of minimizers for the energy functional. The method is applied to the static problem of an elastic body in frictionless contact with a rigid foundation. The contact model allows a small interpenetration of the bodies in the contact region. This interpenetration is modeled by means of a scalar function that depends on the normal component of the displacement field on the potential contact zone. We present the asymptotic behavior of the energy shape functional when a spheroidal void is introduced at an arbitrary point of the elastic body. For the asymptotic analysis, we use a nonoverlapping domain decomposition technique and the associated Steklov– Poincaré pseudodifferential operator. The differentiability of the energy with respect to Communicated by Günter Leugering. S. M. Giusti Departamento de Ingeniería Civil, Facultad Regional Córdoba, Universidad Tecnológica Nacional (UTN/FRC CONICET), Maestro M. López esq. Cruz Roja Argentina, X5016ZAA Córdoba, Argentina e-mail: sgiusti@civil.frc.utn.edu.ar J. Sokołowski (B) Laboratoire de Mathématiques, Institut Élie Cartan, UMR7502 (Université Lorraine, CNRS, INRIA), Université de Lorraine, B.P. 239, 54506 Vandoeuvre-lès-Nancy Cedex, France e-mail: Jan.Sokolowski@univ-lorraine.fr J. Sokołowski Systems Research Institute of the Polish Academy of Sciences, Warsaw, Poland J. Stebel Institute of Mathematics of the Academy of Sciences of the Czech Republic, Žitná 25, Praha 1, 115 67 Prague, Czech Republic e-mail: stebel@math.cas.cz


Introduction
Topological asymptotic analysis [1][2][3] allows us to obtain an asymptotic expansion of a given shape functional of linear elasticity when a geometrical domain is singularly perturbed by the insertion of holes, inclusions, source-terms, or even cracks. The main concept arising from this analysis is the topological derivative. This derivative can be seen as a first-order correction of the unperturbed shape functional to approximate the perturbed shape functional. The topological derivative was rigorously introduced in [4]. Since then, this concept has proved extremely useful in the treatment of a wide range of problems; see, for instance, [5][6][7][8][9][10][11]. Concerning the theoretical development of the topological asymptotic analysis, besides the monograph [1], the reader is referred to [12][13][14][15].
Classically, contact problems are modeled by means of a nonpenetration condition between an elastic body and a rigid obstacle or foundation. This is known as a unilateral contact condition and is modeled using variational inequalities originated by the socalled Signorini Contact Problem. A less restrictive boundary condition on the contact region is obtained by considering the normal compliance model. In this kind of models, based on the assumption of small displacement, some interpenetration between the contacting bodies is allowed, and the boundary forces are given as a function of the interpenetration. However, such models allow an arbitrarily large interpenetration of the bodies in contact, which is physically rather unrealistic. Recently, a new class of models has been presented in [16], using less restrictive boundary conditions that allow small interpenetrations of the bodies. In such a model, the small interpenetration is governed by a function that depends on the normal component of the displacement field on the boundary of the potential contact region. Clearly, this is a nonlinear boundary condition for the contact problem, leading to a new class of variational inequalities.
The shape and topological asymptotic analysis for contact problems have been studied in [17][18][19][20]. In these works, the differentiability of the energy functional with respect to a singular perturbation has been developed for the usual boundary conditions in contact problems. Due to the nonlinear condition over the contact zone, the boundary value problem becomes nonsmooth. Therefore, nonsmooth analysis is necessary since the shape differentiability of solutions to contact problems is obtained only in the framework of Hadamard differentiability of metric projections onto polyhedric sets in appropriate Sobolev spaces.
In this work, we consider the asymptotic behavior of the energy shape functional when a spheroidal void is introduced at an arbitrary point of the elastic body. We Fig. 1 Contact problem consider the energy shape functional associated to the frictionless contact problem allowing a small interpenetration between the elastic body and a rigid foundation, developed in [16]. For the asymptotic analysis, we use a domain decomposition technique and an associated Steklov-Poincaré pseudodifferential operator. The differentiability of the energy of this new class of variational problems, with respect to the nonsmooth perturbation, is established. A closed form for the topological derivative in the three-dimensional space is also presented.
The paper is organized as follows. The frictionless contact problem with small interpenetration is presented in Sect. 2. The topological asymptotic analysis, with respect to the nucleation of spherical holes (voids) in three spatial dimensions, is developed with all details in Sect. 3; here, a closed form of the topological derivatives associated with the energy shape functional is presented. The paper ends with some concluding remarks in Sect. 4.

Static Contact Model for Small Interpenetration
We consider the problem of an elastic body having contact with a rigid foundation. The domain of the body, denoted by ⊂ R 3 , is assumed to be bounded and to have Lipschitz boundary ∂ consisting of three mutually disjoint parts D , N , and C with positive measures, where different boundary conditions are prescribed. On D , we prescribe Dirichlet boundary conditions (displacement), on N Neumann boundary conditions (traction), and on C a contact condition with the rigid foundation that admits interpenetration; see Fig. 1.
For the contact on C , we consider only the normal compliance law of the type where u n := u · n denotes the normal component of the displacement field u, n is the unit outward normal vector to the boundary ∂ , and g is the gap on the potential contact zone. Moreover, in (1), σ n (u) represents the normal component to the boundary of the stress tensor σ (u), i.e., σ n (u) = σ (u)n · n. The Cauchy stress tensor σ (u) is defined as where ε(u) is the symmetric part of the gradient of the displacement field u, i.e., and C denotes the fourth-order elastic tensor. For an isotropic elastic body, this tensor is given by with μ and λ denoting the Lamé coefficients. In the above expression, we use I and I to denote the fourth-order and the second-order identity tensors, respectively. In terms of the engineering constant E (Young's modulus) and υ (Poisson's ratio), the above constitutive response can be written as The function p : R → R + = [0, +∞] in (1) is used to model the interpenetration between the body and the foundation. There are two parameters, β > α, where α indicates the initial contact and the value of β describes a limit such that no further interpenetration is possible. The function p is monotone with the following properties: The strong form of the equilibrium problem under this contact condition is to find a displacement field u : → R 3 such that The last condition in (4) indicates that the contact is without friction, where σ τ (u) = σ (u)n − σ n (u)n denotes the tangential component of the stress tensor σ (u).
We assume that the stress operator σ is bounded and positive definite, i.e., there exist two constants σ , σ > 0 such that The weak formulation of the problem stated in (4) is given by the following variational problem: where the set U of admissible functions is given by and the domain of definition of the function p is For a complete and detailed description of this model, we refer the reader to [16], where it was proved that, under the above assumptions, problem (5) has a unique solution.
The main term of this expansion is the topological derivative operator, which represents a first-order asymptotic correction term of the given shape functional with respect to the singular domain perturbation [4].
Let us consider a shape functional defined on and depending on the solution u, denoted by J (u). After the introduction of a singular perturbation atx, we have a new domain ρ := \ B ρ , where B ρ is the ball of radius ρ centered atx, that is, Therefore, an asymptotic expansion of the energy shape functional defined on the perturbed domain ρ , i.e., J ρ , can be written as is the topological derivative of J atx, and u ρ is the solution of the following contact problem in the perturbed domain: find a displacement field u ρ : Note that there is no traction applied on the boundary of the hole, i.e., the homogeneous Neumann boundary condition is assumed on ∂B ρ . From (8), in accordance with the classical definition of the topological derivative (see [4]), we have In order to perform the asymptotic expansion and evaluation of the topological derivative of problem (9), in this work, we apply a domain decomposition method and the associated Steklov-Poincaré pseudodifferential operator.

Domain Decomposition
We start by decomposing the domain ρ into two parts: (i) a ball B R of radius R > ρ > 0 centered atx ∈ , that is, We use R to denote the exterior boundary ∂B R of the domain C(R, ρ).
(a) Domains Ω R and B R . First, we consider the following linear elasticity problem in C(R, ρ): Using (11), we can define the Steklov-Poincaré boundary operator S ρ on R as where ν denotes the unit normal vector to the boundary R pointing outside the ball B R . Next, we consider the following contact problem in R : find a displacement field Its variational formulation can be written as follows: find a displacement field u R where the set U R of admissible functions is given by From (11) and (12), it follows that the solution u ρ of (9) satisfies Consequently, the restriction of u ρ to the truncated domain R coincides with the solution u R ρ of (13) and similarly u ρ | C(R,ρ) = ω ρ , where ω ρ is the solution to (11) with v = u ρ | R . We also observe that, by the definition of the Steklov-Poincaré boundary operator in C(R, ρ), the solution ω ρ of (11) satisfies For the unperturbed case (ρ = 0), we define the Steklov-Poincaré operator as Applying the domain decomposition technique to the problem (4) on , we can rewrite (5) as follows: It is well known that S ρ is a positive definite operator for any ρ ≥ 0, and that the following asymptotic expansion holds: with a bounded linear operator S [17].

Topological Derivative
For the contact model studied in this work, the energy shape functional associated to the domain is given by [16] J (u) : where u denotes the solution of the problem in the unperturbed domain (see (4)), and the function P(y) is given by Considering the singular perturbation B ρ , the energy shape functional associated to the perturbed domain ρ is given by where u ρ is the solution of the problem in ρ [(see (9)]. Now, by taking into account the domain decomposition and the Steklov-Poincaré boundary operator presented above, we can define the following functional associated to the truncated domain R : In terms of the above functional, the contact problem in the truncated domain R , given by (13), can be written as the following optimization problem: the displacement field u R ρ is the unique minimizer such that where For the optimization problem (24), we can establish the equivalence since the minimizer in (24) coincides with the restriction to R of the minimizer u ρ of the corresponding quadratic functional defined in the whole singularly perturbed domain ρ . (5) and (14), respectively. Then

Proposition 3.1 Let u and u ρ be the solutions to
Proof First, we show that the sequence {u ρ }, ρ ↓ 0, is bounded in H 1 ( R ; R 3 ). Using v := 2u ρ − u as a test function in (14), we obtain The terms on the right-hand side can be estimated using the boundedness of σ , the expression (16), and the properties of the data u and t as follows: Using positive definiteness of σ , the expression (16), and the monotonicity of p, we get a lower bound for the left-hand side of (28): Combining the above estimates with (28), we find that there is a constant C 1 > 0 depending only on σ , σ , u H 1 ( ;R 3 ) and t (H 1/2 Now, we use the embedding H 1 ( R ; R 3 ) → H 1/2 ( N ; R 3 ), Young's inequality, and Korn's inequality in H 1 ( R ; R 3 ) to obtain where C K > 0 is the constant in the Korn inequality and C 2 > 0 depends on the same quantities as C 1 . To show strong convergence, we test (14) with v := u and (18) with v := u ρ . Adding the resulting equations and multiplying by −1, we obtain Since S ρ is a positive definite operator which admits the asymptotic expansion (19), and the sequence {u ρ } is bounded in H 1 R ; R 3 , the last term in (29) satisfies Using this, together with the Korn inequality and the facts that σ is positive definite and p is nondecreasing, we deduce from (29) that as desired. (23), is right-differentiable at ρ = 0 for any fixed R > ρ with ρ ≥ 0, and its derivative is

Proposition 3.2 The functional I R , defined in
where S is the main term of the asymptotic expansion of the Steklov-Poincaré boundary operator S ρ in the space of Steklov-Poincaré operators, given by Proof The derivative of I R at ρ = 0 can be written as Let us consider the following inequalities: Now, for the left-hand side of (33), we have Considering the asymptotic expansion of the Steklov-Poincaré operator, we have Using the strong convergence of u R ρ to u R and the linearity of S , we obtain Now, the right-hand side of (33) can be written as Considering the asymptotic expansion of the Steklov-Poincaré operator, we have By taking the limit of the above expression as ρ ↓ 0, we obtain Finally, the expressions (36) and (39) imply (30).
Remark 3.1 Using Proposition 3.2, the asymptotic expansion of the functional I R can be written as and in view of the asymptotic expansion (8), we finally see that the topological derivative satisfies the following identity: Proposition 3.2 establishes the differentiability of the energy shape functional for this contact model with respect to the nonsmooth perturbation denoted by B ρ . This is an abstract result, whose closed form for the topological derivative T (x) is presented in the next section.

Topological Derivative Evaluation
By the main result of the previous section, the energy shape functional has an asymptotic expansion as ρ ↓ 0; see (8) and (41). This means that the asymptotic behavior of the energy in C(R, ρ) holds in the whole domain . Now, we only need to compute the topological derivative for the energy shape functional in C(R, ρ), with its associated elastic problem (11). In order to evaluate that topological derivative, we can use the techniques available in the literature: see for instance [4,[21][22][23]. Finally, an explicit and analytical formula for the topological derivative T (x) is given in the following result: Theorem 3.1 The energy shape functional of an elastic solid, characterized by the constitutive Eq. (2), with a spherical cavity of radius ρ with homogeneous Neumann boundary condition and centered at a pointx ∈ , admits for ρ ↓ 0 the following asymptotic expansion: where u(x) is the solution of problem (4) evaluated atx, and H is a fourth-order tensor defined as where υ is the Poisson ratio of the elastic medium, and I and I are the second-order and fourth-order identity tensors, respectively.
Proof The reader interested in the proof of this result may refer to [24][25][26].

Remark 3.2
The fourth-order tensor H in (42) can be interpreted as the polarization tensor associated to this problem. This concept is very useful, since the topological derivative formula can be written explicitly in terms of this tensor. The reader interested in this topic may refer to [13,27,28].

Conclusions
Topological derivatives are used in shape optimization [1] to determine the location of small inclusions or voids. There are applications for the optimal design in structural mechanics, or for numerical solution of inverse problems. Therefore, the knowledge of analytical form of topological derivatives is important in numerical methods of shape optimization. In general, the topological derivatives cannot be determined directly for variational inequalities. The existence of topological derivatives for contact problems with friction, dynamic contact problems, or elastic-plastic problems is not known, and it is a difficult open problem. There are no asymptotic analysis tools which can be applied in the general case of nonlinear structural models in mechanics. The only known results on topological derivatives of shape functionals are derived for the static, frictionless contact problems with the unilateral conditions [1]. An analytical expression for the topological derivative of the energy shape functional associated to a new frictionless contact model [16] has been derived. The model allows for a small interpenetration between an elastic body and a rigid foundation. The results obtained can be extended to some nonsmooth domains, including the cracks with nonpenetration conditions. This issue is a subject of current research.
The asymptotic analysis is developed for the specific case of spherical voids introduced in an elastic body. Similar results can be derived for elastic inclusions [1]. The results are obtained by an application of a nonoverlapping domain decomposition technique with the associated Steklov-Poincaré pseudodifferential operators. In this way, the differentiability of the energy can be established. The final formula is a simple analytical expression, in terms of the solution of the state equation and the constitutive parameters, evaluated at each point of the unperturbed domain. Therefore, numerical results can be obtained for this class of models and for shape-topological optimization problems.
We point out that the small interpenetration condition, prescribed in the potential contact zone, does not contribute explicitly to the first-order topological derivative of the energy functional. This means that the formula for the topological derivative of the energy functional for the contact problems takes the same form as was obtained for the classical elasticity problem for an isotropic and homogeneous medium [1,2].
The contribution of the contact model in the topological derivative is through the displacement field, the solution of the contact problem with a nonlinear boundary condition (small interpenetration). This information can be potentially used in the topological design of mechanical components, under contact conditions, to achieve a specified behavior.
The simple formula for the topological derivative of the energy functional for contact problems can be used in numerical methods of optimum design which is the subject of current research in this domain.