Inverse problem for cracked inhomogeneous Kirchhoff–Love plate with two hinged rigid inclusions

We consider a family of variational problems on the equilibrium of a composite Kirchhoff–Love plate containing two flat rectilinear rigid inclusions, which are connected in a hinged manner. It is assumed that both inclusions are delaminated from an elastic matrix, thus forming an interfacial crack between the inclusions and the surrounding elastic media. Displacement boundary conditions of an inequality type are set on the crack faces that ensure a mutual nonpenetration of opposite crack faces. The problems of the family depend on a parameter specifying the coordinate of a connection point of the inclusions. For the considered family of problems, we formulate a new inverse problem of finding unknown coordinates of a hinge joint point. The continuity of solutions of the problems on this parameter is proved. The solvability of this inverse problem has been established. Using a passage to the limit, a qualitative connection between the problems for plates with flat and bulk hinged inclusions is shown.


Introduction
The success of the modern industry is largely based on the widespread use of composite materials. Mathematical modeling of the properties of composites, their behavior under certain conditions provide opportunities for optimizing the parameters and characteristics of composites. As a result, qualitatively new models as well as various complex mathematical approaches leading to new formulations of problems are being developed. Along with the advantages of composite materials, the difference in the physical characteristics of the components can lead to cracks and, as a result, to partial damages to structural elements or destructions of entire composite structures. In some cases, in practice, it is required to assess the possible consequences due to partial damage to some structural elements. Since it is not always possible to find out damage characteristics of hidden structural elements directly, it is necessary to look for methods to determine them indirectly.
For a loaded body, the presence of an inclusion or a crack means that there can be significant stress values in the vicinity of these inhomogeneities. For this reason, the choice of one or another method for specifying the boundary conditions characterizing mechanical processes near a crack can affect the physical adequacy of an entire mathematical model as a whole.
In this work, we follow the nonlinear modeling that uses Signorini-type boundary conditions on the crack faces which was developed in [1][2][3][4][5][6][7][8][9][10][11][12][13] and other works. These nonlinear conditions need to use variational formulations. Within the variational approach, equilibrium problems for composite bodies with elastic or rigid inclusions have been successfully formulated and investigated, see for example [14][15][16][17][18][19][20][21][22][23][24][25][26][27][28]. The other challenge of our work concerns the variation of inclusion shapes. The shape and topological sensitivity analysis of related variational problems was developed in [29][30][31][32][33][34][35] and other papers. In particular, for a crack problem subject to nonpenetration, an optimal control problem concerning minimization of the Griffith functional with respect to the shape and the location of a small elastic inclusion was investigated in [36,37]. Explicit formulae of derivatives of the energy functionals with respect to the shape variation of rigid inclusions were obtained in [19,24]. Inverse problems for the nondestructive determination of elastic and rigid inclusions embedded in Kirchhoff-Love bending elastic plates were investigated in [38][39][40]. We refer the reader to convergence results for elliptic variational-hemivariational inequalities [12] and for double phase obstacle problems with multivalued convection term [13].
Currently we are dealing with a nonlinear model describing a plate with two flat rigid inclusions which are connected in a hinged manner. Namely, we suppose that both inclusions have a common hinge fiber. In addition to this, we assume that both inclusions are delaminated from an elastic media, thus forming an interfacial crack. A new inverse problem predicting the structural integrity of a composite plate is formulated. It is assumed that we know that the plate has two hinged inclusions, but simultaneously we do not know the location of a hinge fiber. Under the assumption that the sum of lengths of the inclusions is constant, a solvability of this inverse problem is proved.

Family of equilibrium problems
Let ⊂ R 2 be a bounded domain with a boundary of class C 1,1 . Suppose that a rectilinear curve γ = (-1, 1) × {0} lies strictly inside , i.e.,γ ⊂ . We assume that consists of two disjoint curves = 0 ∪ 1 , where meas( 0 ) > 0 and meas( 1 ) > 0. We suppose that the domain can be split into two subdomains 1 and 2 with Lipschitz boundaries ∂ 1 and ∂ 2 such that γ ⊂ ∂ 1 ∩ ∂ 2 , and meas(∂ i ∩ 0 ) > 0, i = 1, 2. The latter condition is important for Korn's inequality to hold in the non-Lipschitz domain γ = \γ . For simplicity, the plate is assumed to have a uniform thickness 2h = 2. Let us assign a threedimensional Cartesian space {x 1 , x 2 , z} with the set { γ } × {0} ⊂ R 3 corresponding to the middle plane of the plate. The curve γ defines a crack (a cut) in the plate. This means that the cylindrical surface of the through crack may be defined by the relations x = (x 1 , x 2 ) ∈ γ , -1 ≤ z ≤ 1, where |z| is the distance to the middle plane. Depending on the direction of the normal ν = (0, 1) to γ , we will keep in mind the positive face γ + and the negative face γof the curve γ . The jump [q] of an arbitrary function q across the curve γ is found by the formula [q] = q| γ + -q| γ -. We assume that the plate has two hinged flat inclusions at the positive face of the curve γ . For a fixed parameter δ ∈ [- 1 2 , 1 2 ], the curve γ is split into two rectilinear curves γ δ 1 = (-1, δ) × {0} and γ δ 2 = (δ, 1) × {0} corresponding to two flat rigid inclusions (see Fig. 1). Denote by χ = χ(x) = (W , w) the displacement vector of the mid-surface points (x ∈ γ ), by W = (w 1 , w 2 ) the displacements in the plane {x 1 , x 2 }, and by w the displacements along the axis z. The strain and integrated stress tensors are denoted by ε ij = ε ij (W ), σ ij = σ ij (W ), respectively [1]: where {a ijkl } is the given elasticity tensor, assumed to be symmetric and positive definite A summation convention over repeated indices is used in the sequel. Next we denote the bending moments by the formulas [1] where the tensor {d ijkl } has the same properties as the tensor {a ijkl }. Let B(O, ·, ·) be a bilinear form defined by the equality where χ = (W , w),χ = (W ,w), and O is some measurable subset of . The potential energy functional of the plate has the following representation [1]: 3 describes given body forces [1]. Introduce the Sobolev spaces where n = (n 1 , n 2 ) is the normal vector to . Note that the following inequality with a constant c > 0 independent of χ holds for the bilinear form B( γ , ·, ·) [1].
Remark 1 Inequality (1) yields the equivalence of the standard norm and the semi-norm determined by the bilinear form B( γ , ·, ·) in the space H( γ ). The condition of mutual nonpenetration of opposite faces of the crack is given by [1,16] [ Due to presence of flat rigid inclusions in the plate, restrictions of the functions describing displacements χ to the corresponding curves γ δ 1 , γ δ 2 satisfy a special kind of relations. We introduce the following space which allows us to characterize the properties of the flat rigid inclusions: , and S is some subset of (see [15,16]). Suppose that displacements on two flat rigid inclusions have the following properties [11]: where χ = (W , w), ζ i ∈ R(γ δ i ), i = 1, 2. Since there is the hinge joint of inclusions, we assume that the displacements and angles of rotation of the normal fibers of both inclusions are equal at the corresponding point x δ = (δ, 0) of the middle plane of the plate The variational formulation describing the equilibrium state for the elastic plate with the flat delaminated rigid inclusions can be formulated as follows: For fixed δ ∈ [-1/2, 1/2], let K δ denote the set of admissible displacements .
Remark 3 We should note that the notions of flat inclusion and thin inclusions (see [15,16]) in the same way are specified with the help of cylindrical surfaces, but for flat inclusions there are additional relations describing a connection of plate deflections w and angles of rotation of the normal fibers ∂w ∂x i , i = 1, 2, of an inclusion, see (4). Since bulk (volume) rigid inclusions in Kirchhoff-Love plates are defined with the help of sets of the following type O × [-1, 1], where O is some subdomain lying in the middle plane of the plate [14]; then these relationships (4) for partial derivatives are naturally fulfilled for bulk (volume) rigid inclusions in O.
One can check by simple reasonings that the set K δ is convex and closed. The coercivity and weak semicontinuity of the energy functional together with (1) ensure that problem (7) has a unique solution. Moreover, by the Gateaux differentiability of , problem (7) is equivalent to the following variational inequality (see [1]):

Inverse problem for a hinge fiber
In this section we formulate a new inverse problem for the composite plate containing two flat rigid inclusions. This problem is motivated by prediction of the structural integrity of composite bodies. In order to formulate the inverse problem, we suppose that the parameter δ is assumed to be unknown. Let us introduce the cost functional describing the deviation from some given observation χ g ∈ L 2 ( 1 ) 3 More generally, we consider an arbitrary weakly continuous functional of cost J(χ) : H( γ ) → R. As commonly adopted in optimal control theory, we should either minimize: find p 1 ∈ R such that or maximize the cost: find p 2 ∈ R such that Then, for given intermediate p ∈ [p 1 , p 2 ], we can formulate the inverse problem as follows. Find a parameter δ ∈ [-1/2, 1/2] and a displacement field ξ δ of the body defined in γ such that We prove a solution existence of problem (11), (12) fitting suitable values of p. The following statement takes place.
As the next step we reveal a suitable relationship between the sets of admissible displacements. Namely, we show that for the agreed value δ ∈ [-1/2, 1/2], the sequence {δ k } ⊂ [-1/2, 1/2], and an arbitrary test function χ ∈ K δ , there exists a sequence {χ k } of functions χ k ∈ K δ k converging strongly in H( γ ) to χ as k → ∞. Since x 2 = 0 on γ δ i , i = 1, 2, the following relations hold for χ = (W , w): where the former equation in (15) describes rigid properties of flat inclusions, and (16), (17) reflect the hinged joint at the point (δ, 0). Next, we construct a sequence {χ k } of appropriate vector functions χ k = (W k , w k ) with W k = (w 1k , w 2k ), k = 1, 2, . . . . We define its two components in the following simplest way: while a construction of the second components w 2k of W k needs more detailed reasonings. We start with constructing auxiliary functions μ 1 (x), μ 2 (x), μ 1 k (x), μ 2 k (x), k = 1, 2, . . . , defined by the relations (see Fig. 2): For functionsμ,μ k defined on γ by relations the difference θ k (x) =μ k -μ is well defined for all x ∈ γ . We write out the following representation for this function: for δ k < δ; and for cases when δ k > δ. Note that in both expressions for θ k (x), the following constant value can be estimated with the help of (17): It is evident that θ k (x) ∈ H 1 (γ ) in both cases when δ k < δ and δ k > δ. Since θ k (x) is a piecewise linear function, taking into account (19), we can easily obtain the following estimate: where C > 0 does not depend on k ∈ N.
Using the above construction, the following functions are well defined in the corresponding Sobolev spaces: Applying the linear continuous lifting operator (the corresponding result for lifting operators in the framework of domains with cuts can be found in [1], Lemma 1.12) we can construct functionsw 2 ,w 2k ∈ H 1 ( γ ), k = 1, 2, . . . , with the properties Due to continuity of the lifting operator and the uniform estimate (20), we have Then, for the functions w 2k =w 2k -w 2 + w 2 , it follows By construction, we have the equalities which guarantee the inclusion χ k ∈ K δ k , k ∈ N. Relations (21) and (18) provide the convergence Keeping in mind the convergences (14), (22), we can pass to the limit as k → ∞ in inequalities which gives the limit Therefore, by the arbitrariness of χ ∈ K δ , we reveal thatξ = ξ δ .  (11), (12). The theorem is proved.
Remark 4 As can be seen from the previous reasonings, the segment [-1 2 , 1 2 ] can be replaced with an arbitrary closed segment lying in (-1, 1).