Uniqueness theorem in coupled strain gradient elasticity with mixed boundary conditions

The equilibrium equations and the traction boundary conditions are evaluated on the basis of the condition of the stationarity of the Lagrangian for coupled strain gradient elasticity. The quadratic form of strain energy can be written as a function of the strain and the second gradient of displacement and contains a fourth-, a fifth- and a sixth-order stiffness tensor C4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbb {C}}_4$$\end{document}, C5\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbb {C}}_5$$\end{document} and C6\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbb {C}}_6$$\end{document}, respectively. Assuming invariance under rigid body motions the balance of linear and angular momentum is obtained. The uniqueness theorem (Kirchhoff) for the mixed boundary value problem is proved for the case of the coupled linear strain gradient elasticity (novel). To this end, the total potential energy is altered to be presented as an uncoupled quadratic form of the strain and the modified second gradient of displacement vector. Such a transformation leads to a decoupling of the equation of the potential energy density. The uniqueness of the solution is proved in the standard manner by considering the difference between two solutions.

. In these theories it was assumed that the potential energy density depends not only on the strain, but also on higher derivatives of the displacement vector. More recently, the generalized non-classical theories have been also applied to modeling of materials at the micro-and nano-scale [10,23] to describing of phenomena like dislocations [21], to analyzing of composites with a high difference of the material properties at a lower scale [2,20,42,45,51] to describing some phenomena in regions with stress concentrations [5], to accounting for boundary and surface energies [14,28] or to removing singularities caused by discontinues of boundary conditions (e.g., [6,24,46,49]). It has been shown in numerous papers (see, for example, [22,34,35,41]) that some restrictions of the classical theory of elasticity can be overcome with such gradient expansion.
The uniqueness of solution of equilibrium equations is a crux of any theory. Within the framework of the classical linear elasticity theory the first uniqueness theorem is given for isotropic materials in Kirchhoff [29], which has been completed by Cosserat and Cosserat [11] for boundary value problems with displacement boundary conditions. The generalization to anisotropic materials can be found in Knops and Payne [30] or Bertram and Glüge [9]. A uniqueness theorem for linearly elastic body with surface stresses is proved in Gurtin and Murdoch [26]. Sufficient condition for the uniqueness of the solution of the mixed boundary value problem is the definiteness of the stiffness tetrad C 4 . Requiring additionally stability of the solution yields the condition that C 4 needs to be positive definite. The necessity of this condition has been shown by Cosserat and Cosserat [11].
Modern existence and uniqueness theorems employ different versions of Korn's inequality [31,32], which is an adaption of Poincaré's inequality [43] to problems that involve a projection of the differential operator into the symmetric part. With its help, the requirements for existence and uniqueness of a solution to a variational form of the boundary value problem of classical linear elasticity are reduced to requirements according to the theorem of Babuška-Lax-Milgram [7,33,48].
Conditions sufficient for uniqueness solution are considered for uncoupled isotropic gradient elasticity in Mindlin and Eshel [38]. Usually, positive definiteness of the stored elastic energy is assumed, see Sect. 2, Eq. (1), which reduces to positive definiteness of the stiffness tensors C 4 , and C 6 , while C 5 = O was presumed as well. These are sufficient conditions, but they are surely not necessary. The number of possibilities of boundary value problems in gradient elasticity is much larger compared to classical elasticity. Nevertheless, Mindlin's uniqueness proof is of reasonable generality: it holds for mixed boundary conditions, and no restrictions regarding the shape of the domain are necessary. Neglecting of the constitutive coupling tensor of fifthrank C 5 leads to requirement of positive definiteness of both first and second gradient constitutive stiffness tensors of fourth-and sixth-ranks C 4 and C 6 . Inequality constraints on first and on second gradient constitutive parameters in absence of coupling terms were also presented in Mindlin [36], dell'Isola et al. [12]. In Nazarenko et al. [39] results of Mindlin [36] and dell'Isola et al. [12] have been extended for the case of coupled strain gradient elasticity. In this regard, it has been introduced a block diagonalization of the composite stiffness in strain gradient elasticity. By such a formal transformation, which contains a fourth-, a fifth-and a sixth-order stiffness tensor C 4 , C 5 and C 6 , necessary conditions for positive definiteness and convexity of the isotropic strain and strain gradient energy, accounting for the coupling stiffness C 5 have been obtained. Closed-form relations for the compliance tensors in the frame of the linear theory of coupled gradient elasticity have been given for arbitrary material symmetry classes in Nazarenko et al. [40].
The uniqueness of the solution for an isotropic material featured by four constants (simplified Aifantis model, see, e.g., [3,47]) has been proved for the dynamic evolution problem in Polizzotto [44]. An existence and uniqueness theorem for weak solutions of the equilibrium problem for the simplified case of the linear isotropic strain gradient elasticity (dilatational strain gradient elasticity) is given in Eremeyev et al. [14] and was extended to nonlinear case in Eremeyev et al. [18] and to non-smooth domains with edges in Eremeyev and dell'Isola [15]. In particular, the presence of edges changes the regularity of solutions. In the case of degeneration of the elastic moduli, uniqueness of solutions is proved for a specific boundary conditions, see, e.g., Eremeyev et al. [13,17].
The aim of the paper is to study the uniqueness of solution of the mixed boundary problem in the frame of the linear coupled strain gradient elasticity, which was not done in the previous papers (novelty). This paper has the following structure. Notations used in the paper are introduced in the next section. The equilibrium equation and the traction boundary conditions for the coupled strain gradient elasticity are re-derived on the basis of the principle of virtual power in Sect. 3. The difference between the equations presented here and in Mindlin and Eshel [38] is in index associations of the scalar products. This is the result of the presence of a coupled term in the strain and strain gradient energy Eq. (1) and of the symmetry of the stiffness tensor of fifth-rank. The potential energy density Eq. (1) of Sect. 2 is assumed to be a function of the strain and the second gradient of the displacement vector and contains a fourth-, a fifth-and a sixth-order stiffness tensors C 4 , C 5 and C 6 , respectively. In Sect. 4, the balances of the linear and the angular momentum are employed in the derivation of the equations of strain gradient elasticity. Starting with the balance of the linear and the angular momentum we recover the principle of virtual power. The uniqueness theorem (Kirchhoff) for mixed boundary value problem is extended to the coupled linear strain gradient elasticity in Sect. 5. Conclusions and discussion are presented in the last section.

Notation
Scalars, vectors, second-and higher-rank tensors are denoted by italic letters (like a or A), bold minuscules (like a), bold majuscules (like A), and blackboard bold majuscules (like A), respectively. The basic mathematical operations are given, for example, in Altenbach [4], Eremeyev et al. [16].
The elastic energy density taking into account strain and strain gradient tensors can be written as where C 4 , C 5 , C 6 are the stiffness tensors of 4th, 5th and 6th order, respectively. The strain tensor and the second gradient of displacement vector are is the displacement field, with the position vector of a material point x. For the sake of simplicity, the independent variable x is dropped. ∇ is the three-dimensional nabla operator defined as where e i is an orthonormal base vector. It is implied a summation by the repeating indices. ⊗ denotes the dyadic product. A gradient of the displacement field u is defined by action of the nabla operator on u: On a surface, the spatial nabla operator can be split into a normal and a tangential (surface) part, with and Here I is the three-dimensional unit dyadic. The stresses and the double stresses are defined as The subscript indices correspond the order of tensors E 2 , E 3 , C 4 , C 5 , C 6 , which have the following symmetries (see, e.g., [27]) The dots denote scalar contractions The double and triple scalar contractions in Eqs. (1), (7), (8) act in accordance with the following rules where δ i j is the Kronecker symbol. For a dyadic D 2 , the symmetric part is sym(D 2 ), the skew part is skw(D 2 ), and the axial vector of the skew part of D 2 is noted as axi(D 2 ). axi(D 2 ) is defined by its action on arbitrary vector a according to the rule Here × denotes the vector product.

Variational formulation of the equilibrium equation in the linear strain gradient elasticity
The equilibrium equations and the corresponding natural boundary conditions can be derived on the basis of the Lagrangian variational principle (see, e.g., [1,14]) applied to the coupled linear strain gradient continua.

Variation of the strain energy for coupled linear strain gradient continua
The total potential energy in a volume V with a variation of u is where the stresses T 2 and the double stresses T 3 are determined in Eqs. (7) and (8).
Applying the chain rule and the divergence theorem, the right-hand of Eq. (18) is reduced to the sum of volume and surface integrals where S is the boundary of V and n is the unit outward normal to S. It should be noted, that gradient of the displacement variation δu ⊗ ∇ and variation of the displacement δu on the surface S are dependent. Indeed, knowing δu on S it is possible to determine the surface gradient of δu. Following Toupin [50] and Mindlin [37] we decompose the gradient of δu into its normal and tangential parts in accordance with Eqs. (4)-(6) Using identity Eq. (20) the last term, in the integrand in Eq. (19), can be written down as a sum After applying the chain rule, the last term in Eq. (21) containing the surface gradient can be presented in the following form: We use the surface divergence theorem in the form (see, [37], p. 435) where C is the union of all edges of the domain V , t is the unit tangent to the edge, and m = t × n is the unit outward normal to C tangent to t. Taking v = δu · T 3 · n in (23) gives Using identity (24) and Eq. (21) the last term, in Eq. (19) for the variation of the total potential energy can be written down in the following form: Substituting this into Eq. (19) we obtain Remark 1 There is a qualitative difference between converting volume integrals into surface integrals and surface integrals into edge integrals. While the surface of the volume divides space into inside and outside of the body, this does not hold for edges on the surface. The surface is, so to say, unbounded, unlike the volume. The surface divergence theorem produces two integrals (Eq. 24), the first contains (n · ∇ s )dS and the second contains δu · [T 3 · ·n ⊗ m] dC. Both can be zero, but not at the same time: -for bodies bounded only by planar facets, n · ∇ s = 0, hence the first integral vanishes, and -for bodies without edges, the edge integral disappears because the surface is unbounded.
One might argue that the second case is the more general case, as one can regularize sharp corners by a tiny but smooth edge radius, and get the first case as the limit to a zero edge radius. Therefore, for simplicity the edge integrals can be summarized into the surface integrals, in the sense of a regularization, which allows to drop the edge integrals. In the remainder we keep the edge integrals for the sake of completeness, but it may be possible to neglect the edges and all the quantities that are defined on them if one is willing to regularize sharp corners, which simplifies the derivation considerably.

Variational equation of equilibrium and boundary conditions
The Lagrangian l can be written down in the following form: where a is a work of external forces and double forces and w is a strain and strain gradient energy density Eq. (1). Then the principle of the Lagrangian stationary (the principle of virtual power) is The variation of the strain and strain gradient energy requires an admissible form of the work: Here f is a body force per unit volume, p is a surface traction, R is a surface normal double force on S (e.g., [14]), and c is a line force on edge C.
Since the normal gradient of displacement is described as or the work of external forces can be written as where r n is a double traction in normal direction on S (e.g., [8]) For all admissible functions δu the variation of the potential energy has a form which leads with Eq. (26) to the equation of equilibrium: and to the dynamic boundary conditions for the surface tractions and the edge forces, which can be prescribed: -the vector field of the tractions on the part of the surface of the body S d -the double tractions in normal direction on the S d -the line forces on edge on the part of edge C d c pr = T 3 · ·n ⊗ m .
Here, the subscript pr denotes prescribed. The displacement or kinematic boundary conditions in terms of the displacement fields u and its normal gradient Du on the part of the surface of the body S g (S d ∪ S g = S), and the displacement u on the part of the edge are apparent from Eq. (33)

Equilibrium equation based on the balance of the virtual power
Let us consider the Lagrangian variational principle (principle of virtual power) Eq. (33) for rigid body motion where u 0 and ω are two constant vectors. For such δu from Eq. (33) it is possible to obtain the equilibrium equation for a free solid body (e.g., [9,14]) Thus we have two balance laws for the total force and the total torque in the form: -the balance of linear momentum -the balance of angular momentum where the tractions p on S, and the line forces on edge c on C are described as and m is the torque induced by double forces on the surface of the body by the rotational field ω × x with For an arbitrary rotational field ω × x, where x is position vector and ω -any constant vector is antisymmetric. Thus and, consequently Here the axial part of the vector axi(ω×n⊗n) is defined according to Eq. (17) and its normal part axi n (ω×n⊗n) analogically as for normal gradient Eqs. (4)-(6).
With the above identity we achieve the balance of angular momentum in the form [8] V It is possible to obtain from the balance of linear and angular momentum the local balance of linear momentum or first Eulerian law of motion Since p is determined from Eq. (45) and using the following integral transformation one can show that axi(T 2 ) = 0 and as a consequence the symmetry of the stress tensor (second Eulerian law of motion) follows.
We are now able to reformulate the principle of virtual power.

Theorem 1 A motion of the body is dynamically admissible if and only if the balance of virtual power in the form
holds for all differentiable vector fields δu.
Proof Let us assume that a motion of the body is dynamically admissible so that Eqs. (53)-(55) hold. Then the balance of linear momentum is fulfilled if holds in every point, and also for any differentiable vector fields δu. The integral over the body also vanishes Accounting for Eqs. (18), (26) we obtain In addition, the balance of moment of momentum Eq. (55) is fulfilled if or V T 2 · ·skw(δu ⊗ ∇) dV = 0 ( 6 2 ) holds for any differentiable vector fields δu. We subtract this from equation (60) and obtain Eq. (56).

Extension of Kirchhoff's uniqueness theorem to coupled strain gradient elasticity
In this subsection it is demonstrated that the solution of the mixed boundary value problem, if exists, is unique. We consider a mixed boundary value problem: -the stress equation of equilibrium Both parts of the surface S d and S g and the edge C d and C g should be almost everywhere disjoint and should form the entire surface and its edges Since all equations of the gradient elasticity theory are linear, we can use the principle of superposition: where u 0 and x 0 are two constant vectors and is a constant anti-symmetric tensor. If the displacements are prescribed for at least three not collinear points then the solution is unique Proof The theorem assumptions and the principle of superposition lead to the following stating: is a solution of the mixed boundary value problem with zero-boundary conditions: f = 0 in V , p = 0 and r n = 0 on S d , c = 0 on C d , δu = 0 and Dδu = 0 on S g . Following the principle of virtual displacements Eqs. (33) and (56) we have for the zero-boundary value problem considering that at least one factor under each integral on the left-hand side of the above equation is zero in the integration domain. Here the constitutive relations are defined as where With Eqs. (69) and (70) the integrand in equation (68) can be presented as The presence on the right-hand side of the coupling term C 5 complicates essential obtaining the inequality constraints needed for positive definiteness of w.
In order to separate the contributions of the strain and of the strain gradient, this equation can be transformed by using a formal modification of the second gradient of displacement and of the stiffness tensor of fourth-rank as shown in Nazarenko et al. [39]: The modified second gradient of displacement and the modified stiffness tensor are marked by superscript m, and are specified as and Because of the assumed positive definiteness of the quadratic form of the potential energy density, δE 2 · ·C m 4 · ·δE 2 + δE m 3 · · · C 6 · · · δE m 3 (76) must be nonnegative everywhere. Then integrating the last expression gives V δE 2 · ·C m 4 · ·δE 2 + δE m 3 · · · C 6 · · · δE m 3 dV and after Korn's inequality (see, e.g., [19,30]) can be zero, if and only if the integrand δu ⊗ ∇ is zero everywhere. Here the magnitude of a second-rank tensor M is defined as and k is a constant depending only on the region of integration V . Therefore, the displacement δu can be strain-free only, and with Eq. (67) It is the proof of the first part of the theorem. The second part of the theorem proof is rather evident, indeed, if S g has at least three not collinear points, for which the displacements are prescribed, then δu must vanish in these points, and that δu will be zero field everywhere Remark 2 We would like to point out that Eq. (72) can also be modified in order to present the variation of the strain energy density as block matrices by introducing the variation of the modified strain and the modified stiffness tensor of sixth-rank (see [40]) with and Such a transformation is mathematically equivalent to Eq. (73) and leads to the same proving scheme. Thus, the both modifications can be used for proving and can be considered as an additional validation of the presented proof.

Conclusions
The equilibrium equation and the corresponding natural boundary conditions are re-derived on the basis of the Lagrange variational principle modified for the coupled linear strain gradient continua. The potential energy density Eq. (1) is assumed to be a function of the strain and the second gradient of displacement and contains a fourth-, a fifth-and a sixth-order stiffness tensor C 4 , C 5 and C 6 . Assuming rigid body motion the balance lows for the total forces and total torque are obtained from the principle of virtual power. These balance lows are employed in the derivation of the equilibrium equations of the coupled strain gradient elasticity. It is shown that the balance of angular momentum leads to the symmetry of the stress tensor T 2 as a consequence. Then starting with the balance of the linear and the angular momentum the principle of the virtual power is recovered.
Kirchhoff's uniqueness theorem for the mixed boundary value problem is proven for the equilibrium problem for the coupled linear strain gradient elasticity. The existence of the coupling term C 5 makes the problem more complicated. By a transformation which may be considered as a block diagonalization, the equation of the potential energy density Eq. (1) was modified to present w as an uncoupled quadratic form of the strain and the modified second gradient of displacement (see details in [39,40]). Then, one may apply Mindlin's uniqueness proof for isotropic materials. Here we proved the uniqueness of the solution in the usual way by contradiction by considering two solutions and their difference, similar to Kirchhoff's proof in classical elasticity. The proof requires positive definiteness of the potential energy density. The positive definiteness conditions for constitutive parameters are presented in Nazarenko et al. [39] in a tensorial form without restricting the symmetry class, but the inequalities are therein given for the 8 elasticity parameters in hemitropy (invariance under the action of the special orthogonal group).