Abstract
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 the nonsmooth perturbation is established, and the topological derivative is presented in the closed form.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
Topological asymptotic analysis [1–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–11]. Concerning the theoretical development of the topological asymptotic analysis, besides the monograph [1], the reader is referred to [12–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 so-called 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–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 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.
2 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 \(\Omega \subset \mathbb {R}^3\), is assumed to be bounded and to have Lipschitz boundary \(\partial \Omega \) consisting of three mutually disjoint parts \(\Gamma _D\), \(\Gamma _N\), and \(\Gamma _C\) with positive measures, where different boundary conditions are prescribed. On \(\Gamma _D\), we prescribe Dirichlet boundary conditions (displacement), on \(\Gamma _N\) Neumann boundary conditions (traction), and on \(\Gamma _C\) a contact condition with the rigid foundation that admits interpenetration; see Fig. 1.
For the contact on \(\Gamma _C\), we consider only the normal compliance law of the type
where \(u_n:=\varvec{u} \cdot \varvec{n}\) denotes the normal component of the displacement field \(\varvec{u}\), \(\varvec{n}\) is the unit outward normal vector to the boundary \(\partial \Omega \), and \(g\) is the gap on the potential contact zone. Moreover, in (1), \(\sigma _{n}(\varvec{u})\) represents the normal component to the boundary of the stress tensor \(\varvec{\sigma }(\varvec{u})\), i.e., \(\sigma _{n}(\varvec{u})=\varvec{\sigma }(\varvec{u})\varvec{n}\cdot \varvec{n}\). The Cauchy stress tensor \(\varvec{\sigma }(\varvec{u})\) is defined as
where \(\varvec{\varepsilon }(\varvec{u})\) is the symmetric part of the gradient of the displacement field \(\varvec{u}\), i.e.,
and \(\mathbb {C}\) denotes the fourth-order elastic tensor. For an isotropic elastic body, this tensor is given by
with \(\mu \) and \(\lambda \) denoting the Lamé coefficients. In the above expression, we use \(\mathbb {I}\) and \(\mathrm {\varvec{I}}\) to denote the fourth-order and the second-order identity tensors, respectively. In terms of the engineering constant \(E\) (Young’s modulus) and \(\upsilon \) (Poisson’s ratio), the above constitutive response can be written as
The function \(p:\mathbb {R}\rightarrow \overline{\mathbb {R}}_+=[0,+\infty ]\) in (1) is used to model the interpenetration between the body and the foundation. There are two parameters, \(\beta > \alpha \), where \(\alpha \) indicates the initial contact and the value of \(\beta \) 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 \(\varvec{u}:\Omega \rightarrow \mathbb {R}^3\) such that
The last condition in (4) indicates that the contact is without friction, where \(\varvec{\sigma }_{\tau }(\varvec{u}) = \varvec{\sigma }(\varvec{u})\varvec{n} - \sigma _{n}(\varvec{u}) \varvec{n}\) denotes the tangential component of the stress tensor \(\varvec{\sigma }(\varvec{u})\).
We assume that the stress operator \(\varvec{\sigma }\) is bounded and positive definite, i.e., there exist two constants \(\underline{\sigma },\overline{\sigma }>0\) such that
and the data satisfy
The weak formulation of the problem stated in (4) is given by the following variational problem: find \(\varvec{u} \in \mathcal {U}\) with \(u_n - g \in \mathrm {dom}(p)\) such that
where the set \(\mathcal {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.
3 Topological Asymptotic Analysis
In this section, we obtain an asymptotic expansion for the energy shape functional when a small spheroidal cavity of radius \(\rho \) is introduced at an arbitrary point \(\hat{\varvec{x}}\) of \(\Omega \), far enough from the potential contact region \(\Gamma _C\); see Fig. 2.
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 \(\Omega \) and depending on the solution \(\varvec{u}\), denoted by \(\mathcal {J}_{\Omega } (\varvec{u})\). After the introduction of a singular perturbation at \(\hat{\varvec{x}}\), we have a new domain \(\Omega _\rho := \Omega \setminus \overline{\mathcal {B}_\rho }\), where \(\mathcal {B}_\rho \) is the ball of radius \(\rho \) centered at \(\hat{\varvec{x}}\), that is, \(\mathcal {B}_\rho :=\left\{ \varvec{x}\in \mathbb {R}^3: |\varvec{x}-\hat{\varvec{x}}| < \rho \right\} \); see Fig. 2.
Therefore, an asymptotic expansion of the energy shape functional defined on the perturbed domain \(\Omega _\rho \), i.e., \(\mathcal {J}_{\Omega _\rho }\), can be written as
where \(f(\rho )\) is a decreasing positive function such that \(f(\rho )\rightarrow 0\) as \(\rho \downarrow 0\), \(\mathcal {T}_\Omega (\hat{\varvec{x}})\) is the topological derivative of \(\mathcal {J}_{\Omega }\) at \(\hat{\varvec{x}}\), and \(\varvec{u}_\rho \) is the solution of the following contact problem in the perturbed domain: find a displacement field \(\varvec{u}_\rho :\Omega _\rho \rightarrow \mathbb {R}^3\) such that
Note that there is no traction applied on the boundary of the hole, i.e., the homogeneous Neumann boundary condition is assumed on \(\partial \mathcal {B}_\rho \).
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.
3.1 Domain Decomposition
We start by decomposing the domain \(\Omega _\rho \) into two parts: (i) a ball \(\mathcal {B}_R\) of radius \(R>\rho >0\) centered at \(\hat{\varvec{x}}\in \Omega \), that is, \(\mathcal {B}_R:=\left\{ \varvec{x}\in \mathbb {R}^3: |\varvec{x}-\hat{\varvec{x}}| < R \right\} \), and (ii) the domain \(\Omega _R := \Omega \setminus \overline{\mathcal {B}_R}\). Clearly, \(\mathcal {B}_R\) contains the small cavity \(\mathcal {B}_\rho \), and, for this perturbed configuration, we can define the domain as \(C(R,\rho ):=\mathcal {B}_R\setminus \overline{\mathcal {B}_\rho }\); see Fig. 3.
We use \(\Gamma _R\) to denote the exterior boundary \(\partial \mathcal {B}_R\) of the domain \(C(R,\rho )\).
First, we consider the following linear elasticity problem in \(C(R,\rho )\): given \(\overline{ \varvec{v}}\in H^{1/2}(\Gamma _R;\mathbb {R}^3)\), find a displacement field \(\varvec{\omega }_\rho :C(R,\rho ) \rightarrow \mathbb {R}^3\) such that
Using (11), we can define the Steklov–Poincaré boundary operator \(\mathcal {S}_\rho \) on \(\Gamma _R\) as
where \(\varvec{\nu }\) denotes the unit normal vector to the boundary \(\Gamma _R\) pointing outside the ball \(\mathcal B_R\). Next, we consider the following contact problem in \(\Omega _R\): find a displacement field \(\varvec{u}_\rho ^R: \Omega _R \rightarrow \mathbb {R}^3\) such that
Its variational formulation can be written as follows: find a displacement field \(\varvec{u}_\rho ^R \in \mathcal {U}_R\) with \(p(\varvec{u}^R_{\rho }\cdot \varvec{n} - g) \in \mathrm {dom}(p)\) such that
where the set \(\mathcal {U}_R\) of admissible functions is given by
From (11) and (12), it follows that the solution \(\varvec{u}_\rho \) of (9) satisfies
Consequently, the restriction of \(\varvec{u}_\rho \) to the truncated domain \(\Omega _R\) coincides with the solution \(\varvec{u}_\rho ^R\) of (13) and similarly \(\varvec{u}_\rho |_{C(R,\rho )}=\varvec{\omega }_\rho \), where \(\varvec{\omega }_\rho \) is the solution to (11) with \(\overline{\varvec{v}}=\varvec{u}_\rho |_{\Gamma _R}\). We also observe that, by the definition of the Steklov–Poincaré boundary operator in \(C(R,\rho )\), the solution \(\varvec{\omega }_\rho \) of (11) satisfies
For the unperturbed case (\(\rho =0\)), we define the Steklov–Poincaré operator as
associated with the problem
Applying the domain decomposition technique to the problem (4) on \(\Omega \), we can rewrite (5) as follows:
It is well known that \(\mathcal S_\rho \) is a positive definite operator for any \(\rho \ge 0\), and that the following asymptotic expansion holds:
with a bounded linear operator \(\mathcal S'\) [17].
3.2 Topological Derivative
For the contact model studied in this work, the energy shape functional associated to the domain \(\Omega \) is given by [16]
where \(\varvec{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 \(\mathcal {B}_\rho \), the energy shape functional associated to the perturbed domain \(\Omega _\rho \) is given by
where \(\varvec{u}_\rho \) is the solution of the problem in \(\Omega _\rho \) [(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 \(\Omega _R\):
In terms of the above functional, the contact problem in the truncated domain \(\Omega _R\), given by (13), can be written as the following optimization problem: the displacement field \(\varvec{u}_\rho ^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 \(\Omega _R\) of the minimizer \(\varvec{u}_\rho \) of the corresponding quadratic functional defined in the whole singularly perturbed domain \(\Omega _\rho \).
Proposition 3.1
Let \(\varvec{u}\) and \(\varvec{u}_\rho \) be the solutions to (5) and (14), respectively. Then
Proof
First, we show that the sequence \(\{\varvec{u}_\rho \}\), \(\rho \downarrow 0\), is bounded in \(H^1(\Omega _R;\mathbb {R}^3)\). Using \(\varvec{v}:=2 \varvec{u}_\rho -\overline{\varvec{u}}\) as a test function in (14), we obtain
The terms on the right-hand side can be estimated using the boundedness of \(\varvec{\sigma }\), the expression (16), and the properties of the data \(\overline{\varvec{u}}\) and \(\overline{\varvec{t}}\) as follows:
Using positive definiteness of \(\varvec{\sigma }\), 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 \(\underline{\sigma }\), \(\overline{\sigma }\), \(\Vert \overline{\varvec{u}}\Vert _{H^1(\Omega ;\mathbb {R}^3)}\) and \(\Vert \overline{\varvec{t}}\Vert _{(H^{1/2}(\Gamma _N;\mathbb {R}^3))^*}\) such that
Now, we use the embedding \(H^1(\Omega _R;\mathbb {R}^3)\hookrightarrow H^{1/2}(\Gamma _N;\mathbb {R}^3)\), Young’s inequality, and Korn’s inequality in \(H^1(\Omega _R;\mathbb {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 \(\varvec{v}:=\varvec{u}\) and (18) with \(\varvec{v}:=\varvec{u}_\rho \). Adding the resulting equations and multiplying by \(-1\), we obtain
Since \(\mathcal S_\rho \) is a positive definite operator which admits the asymptotic expansion (19), and the sequence \(\{\varvec{u}_\rho \}\) is bounded in \(H^1\big (\Omega _R;\mathbb {R}^3\big )\), the last term in (29) satisfies
Using this, together with the Korn inequality and the facts that \(\varvec{\sigma }\) is positive definite and \(p\) is nondecreasing, we deduce from (29) that
as desired. \(\square \)
Proposition 3.2
The functional \(\mathcal {I}_{\Omega _R}\), defined in (23), is right-differentiable at \(\rho = 0\) for any fixed \(R>\rho \) with \(\rho \ge 0\), and its derivative is
where \(\mathcal {S}^\prime \) is the main term of the asymptotic expansion of the Steklov–Poincaré boundary operator \(\mathcal {S}_\rho \) in the space of Steklov–Poincaré operators, given by
Proof
The derivative of \(\mathcal {I}_{\Omega _R}\) at \(\rho = 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 \(\varvec{u}_\rho ^R\) to \(\varvec{u}^R\) and the linearity of \(\mathcal {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 \(\rho \downarrow 0\), we obtain
Finally, the expressions (36) and (39) imply (30).\(\Box \)
Remark 3.1
Using Proposition 3.2, the asymptotic expansion of the functional \(\mathcal {I}_{\Omega _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 \(\mathcal {B}_\rho \). This is an abstract result, whose closed form for the topological derivative \(\mathcal {T}_\Omega (\hat{\varvec{x}})\) is presented in the next section.
3.3 Topological Derivative Evaluation
By the main result of the previous section, the energy shape functional has an asymptotic expansion as \(\rho \downarrow 0\); see (8) and (41). This means that the asymptotic behavior of the energy in \(C(R,\rho )\) holds in the whole domain \(\Omega \). Now, we only need to compute the topological derivative for the energy shape functional in \(C(R,\rho )\), 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–23]. Finally, an explicit and analytical formula for the topological derivative \(\mathcal {T}_\Omega (\hat{\varvec{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 \(\rho \) with homogeneous Neumann boundary condition and centered at a point \(\hat{\varvec{x}}\in \Omega \), admits for \(\rho \downarrow 0\) the following asymptotic expansion:
where \(\varvec{u}(\hat{\varvec{x}})\) is the solution of problem (4) evaluated at \(\hat{\varvec{x}}\), and \(\mathbb {H}\) is a fourth-order tensor defined as
where \(\upsilon \) is the Poisson ratio of the elastic medium, and \(\mathrm {\varvec{I}}\) and \(\mathbb {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–26]. \(\square \)
Remark 3.2
The fourth-order tensor \(\mathbb {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].
4 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.
References
Novotny, A.A., Sokołowski, J.: Topological Derivatives in Shape Optimization. Interaction of Mechanics and Mathematics. Springer, Heidelberg (2013)
Sokołowski, J., Żochowski, A.: Topological derivatives of shape functionals for elasticity systems. Mech. Struct. Mach. 29, 331–349 (2001)
Nazarov, S.A., Sokołowski, J.: Asymptotic analysis of shape functionals. J. Math. Pure Appl. 82(2), 125–196 (2003)
Sokołowski, J., Żochowski, A.: On the topological derivative in shape optimization. SIAM J. Control Optim. 37(4), 1251–1272 (1999)
Amstutz, S., Novotny, A.A., de Souza Neto, E.A.: Topological derivative-based topology optimization of structures subject to Drucker–Prager stress constraints. Comput. Meth. Appl. Mech. Eng. 233–236, 233–236 (2012)
Hintermüller, M., Laurain, A., Novotny, A.A.: Second-order topological expansion for electrical impedance tomography. Adv. Comput. Math. 36, 235–265 (2012)
Hintermüller, M., Laurain, A.: Multiphase image segmentation and modulation recovery based on shape and topological sensitivity. J. Mmath. Imaging Vis. 35, 1–22 (2009)
Giusti, S.M., Novotny, A.A., de Souza Neto, E.A., Feijóo, R.A.: Sensitivity of the macroscopic elasticity tensor to topological microstructural changes. J. Mech. Phys. Solids 57(3), 555–570 (2009)
Van Goethem, N., Novotny, A.A.: Crack nucleation sensitivity analysis. Math. Method Appl. Sci. 33(16), 1978–1994 (2010)
Novotny, A.A., Sokołowski, J., de Souza Neto, E.A.: Topological sensitivity analysis of a multi-scale constitutive model considering a cracked microstructure. Math. Method. Appl. Sci. 33(5), 676–686 (2010)
Khludnev, A.M., Sokołowski, J., Szulc, K.: Shape and topological sensitivity analysis in domains with cracks. Appl. Math. 55(6), 433–469 (2010)
Sokołowski, J., Żochowski, A.: Optimality conditions for simultaneous topology and shape optimization. SIAM J. Control Optim. 42(4), 1198–1221 (2003)
Cardone, G., Nazarov, S.A., Sokołowski, J.: Asymptotic analysis, polarization matrices, and topological derivatives for piezoelectric materials with small voids. SIAM J. Control Optim. 48(6), 3925–3961 (2010)
Iguernane, M., Nazarov, S.A., Roche, J.-R., Sokołowski, J., Szulc, K.: Topological derivatives for semilinear elliptic equations. Int. J. Appl. Math. Comput. Sci. 19(2), 191–205 (2009)
Frémiot, G., Horn, W., Laurain, A., Rao, M., Sokołowski, J.: On the analysis of boundary value problems in nonsmooth domains. Dissert. Math. 462, 149 pp. (2009).
Eck, C., Jarušek, J., Stará, J.: Normal compliance contact models with finite interpenetration. Arch. Ration. Mech. Anal. 208(1), 25–57 (2013)
Sokołowski, J., Żochowski, A.: Modelling of topological derivatives for contact problems. Numer. Math. 102(1), 145–179 (2005)
Fulmanski, P., Lauraine, A., Scheid, J.F., Sokołowski, J.: A level set method in shape and topology optimization for variational inequalities. Inter. J. Appl. Math. Comput. Sci. 17(3), 413–430 (2007)
Jarušek, J., Krbec, M., Rao, M., Sokołowski, J.: Conical differentiability for evolution variational inequalities. J. Differ. Equ. 193(1), 131–146 (2003)
Argatov, I.I., Sokołowski, J.: Asymptotics of the energy functional in the Signorini problem under small singular perturbation of the domain. Comp. Math. Math. Phys. 43(5), 710–724 (2003)
Amstutz, S.: Aspects Théoriques et Numériques en Optimisation de Forme Topologique. Phd thesis, Institut National des Sciences Appliquées de Toulouse, Toulouse (2003).
Novotny, A.A., Feijóo, R.A., Padra, C., Taroco, E.: Topological sensitivity analysis. Comput. Meth. Appl. Mech. Eng. 192(7–8), 803–829 (2003)
Laurain, A.: Singularly Perturbed Domains in Shape Optimisation. Phd thesis, Université Henri Poincaré, Nancy (2006).
Novotny, A.A., Feijóo, R.A., Taroco, E., Padra, C.: Topological sensitivity analysis for three-dimensional linear elasticity problem. Comput. Meth. Appl. Mech. Eng. 196(41–44), 4354–4364 (2007)
Hlaváček, I., Novotny, A.A., Sokołowski, J., Żochowski, A.: On topological derivatives for elastic solids with uncertain input data. J. Optimiz. Theory Appl. 141(3), 569–595 (2009)
Khludnev, A.M., Novotny, A.A., Sokołowski, J., Żochowski, A.: Shape and topology sensitivity analysis for cracks in elastic bodies on boundaries of rigid inclusions. J. Mech. Phys. Solids 57(10), 1718–1732 (2009)
Ammari, H., Kang, H.: Polarization and Moment Tensors with Applications to Inverse Problems and Effective Medium Theory. Applied mathematical sciences. Springer-Verlag, New York (2007)
Nazarov, S.A., Sokołowski, J., Specovius-Neugebauer, M.: Polarization matrices in anisotropic heterogeneous elasticity. Asymptot. Anal. 68(4), 181–221 (2010)
Acknowledgments
This research was partially supported by LabEx CARMIN–CIMPA SMV programme (France). The work of Jan Sokolowski was partially supported by the ANR-12- BS01-0007 Optiform. The work of J. Stebel was supported by the ESF grant Optimization with PDE Constraints, by the Czech Science Foundation (GAČR) Grant No. 201/09/0917 and RVO 67985840.
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by Günter Leugering.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited.
About this article
Cite this article
Giusti, S., Sokołowski, J. & Stebel, J. On Topological Derivatives for Contact Problems in Elasticity. J Optim Theory Appl 165, 279–294 (2015). https://doi.org/10.1007/s10957-014-0594-7
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10957-014-0594-7
Keywords
- Topological derivative
- Static frictionless contact problem
- Asymptotic analysis
- Domain decomposition
- Steklov–Poincaré operator