## Abstract

In this paper, it is proven an existence and uniqueness theorem for weak solutions of the equilibrium problem for linear isotropic dilatational strain gradient elasticity. Considered elastic bodies have as deformation energy the classical one due to Lamé but augmented with an additive term that depends on the norm of the gradient of dilatation: only one extra second gradient elastic coefficient is introduced. The studied class of solids is therefore related to Korteweg or Cahn–Hilliard fluids. The postulated energy naturally induces the space in which the aforementioned well-posedness result can be formulated. In this energy space, the introduced norm does involve the linear combination of some specific higher-order derivatives only: it is, in fact, a particular example of anisotropic Sobolev space. It is also proven that aforementioned weak solutions belongs to the space \(H^1(div,V)\), i.e. the space of \(H^1\) functions whose divergence belongs to \(H^1\). The proposed mathematical frame is essential to conceptually base, on solid grounds, the numerical integration schemes required to investigate the properties of dilatational strain gradient elastic bodies. Their energy, as studied in the present paper, has manifold interests. Mathematically speaking, its singularity causes interesting mathematical difficulties whose overcoming leads to an increased understanding of the theory of second gradient continua. On the other hand, from the mechanical point of view, it gives an example of energy for a second gradient continuum which can sustain externally applied surface forces and double forces but cannot sustain externally applied surface couples. In this way, it is proven that couple stress continua, introduced by Toupin, represent only a particular case of the more general class of second gradient continua. Moreover, it is easily checked that for dilatational strain gradient continua, balance of force and balance of torques (or couples) are not enough to characterise equilibrium: to this aim, externally applied surface double forces must also be specified. As a consequence, the postulation scheme based on variational principles seems more suitable to study second gradient continua. It has to be remarked finally that dilatational strain gradient seems suitable to model the experimentally observed behaviour of some material used in 3D printing process.

## Introduction

Strain gradient elasticity deals with those models of continuum media where a strain energy density depends on the first and second gradients of placements. The mechanical motivation which requires the introduction of such an energy is related to the existence, at lower scales with respect to the scale where continuum behaviour is observed, of longer distance interactions among the basic constituents of the considered deformable body, see, e.g. [58, 59, 68].

Strain gradient elasticity has been originated not so recently. It is a particular case of the higher gradient continua already considered in the works by Gabrio Piola, who also recognised the mechanical mechanism, at lower scales, which may require the introduction of this class of generalised continua. Later, Le Roux and other scholars see, e.g. [21, 51, 52] and the reference therein, continued the line of thought started by Piola and, more recently, also by using the by using a more modern formalism and the powerful tools given by functional analysis, Paul Germain gave a further impulse to generalised continuum mechanics, see [28, 42, 43].

In the case of fluids, the strain gradient model, which was also at first formulated by Piola, was studied by Korteweg [49] and by Cahn and Hilliard [15, 16] see also [8, 29, 67]. It has to be remarked that Piola deduced the strong form of the evolution equations for continua (and in particular for fluids) using a variational principle, by calculating the corresponding Euler-Lagrange stationarity conditions in the referential (Lagrangian) description and finally transporting (using so-called Piola transport) these conditions in the current (Eulerian) description [8, 29].

Within rational mechanics approach, the corresponding and needed boundary value problems were fully formulated by Toupin [73, 74] and Mindlin [53, 54] by using some techniques from the differential geometry of surfaces, see also [25] for *N*th gradient media. More recently, the model was also applied for modelling materials at the micro- and nanometer scale [20, 40], or composites with high contrast (at a lower scale) in material properties [2, 22, 24, 26, 37, 61, 62, 69, 77].

In statics, the boundary value problems, which are deduced from the stationarity condition for a second gradient energy functional, constitute a system of fourth-order partial differential equations complemented by consistent boundary conditions. Consistent boundary conditions, as already recognised by Lagrange and Piola, can be easily deduced by formulating a variational principle (either a minimisation principle or a stationarity condition) and then integrating by parts to get the so-called strong formulation of equilibrium condition.

For infinitesimal deformations one assumes that a strain energy density *W* depends on the first and second gradients of displacements \({\mathbf {u}}\), so that we have \(W=W(\nabla {\mathbf {u}},\nabla \nabla {\mathbf {u}})\), where *W* is a quadratic form of its arguments. We call a model *gradient complete* if the following inequality is valid

where \(C_1\) and \(C_2\) are positive constants which are independent on \({\mathbf {u}}\) and \(\Vert \cdot \Vert \) denotes the Euclidean norm. In this case, the analysis of the well posedness is based on the general theory of elliptic PDEs of higher order [4, 38, 48]. For the strain gradient elasticity such analysis was provided in [47, 50]. In this case, we have a unique weak solution which belongs to \(H^2(V)\).

It is however remarkable that some interesting mechanical problems require the introduction of non-complete second gradient continuum models [9, 13, 24, 26, 32]. Therefore, we must conclude that in addition to gradient complete models such as those included in the so-called Toupin–Mindlin strain gradient elasticity, one must consider also models where a strain energy density depends on a subset of second derivatives of displacements or on particular combinations of them. We call such models *gradient incomplete*. It is worth to mention explicitly now the so-called couple stress theory [56, 78] with the strain energy in the form \(W=W(\nabla {\mathbf {u}},\nabla \times {\mathbf {u}})\), where “\(\times \)” stands for the cross product. This special case of second gradient continuum models have a peculiarity which made them more acceptable by those who want to postulate continuum mechanics on the primitive concept of “balance laws” for forces and couples (or torques). In facts, [21] in the framework of couple-stress theories, the admissible boundary conditions (naturally deduced from the adapted variational principle) assume a particular for which double-forces are vanishing, see [28, 43] for a detailed discussion of the concept of double force and Alibert et al. [5] for its interpretation at micro-scales. Let us also note that the concept of double forces and other hypermoments are known in structural mechanics, see comments in [31, pp. 282–283]. As it is rather difficult to include the boundary conditions involving double forces in the postulation scheme based on balance laws, then its supporters have the tendency to deny the physical importance of said double forces. Instead, we believe in the importance of the role of double forces in continuum mechanics and, following the point of view expressed by Toupin [73, 74], we believe that the postulation scheme based on variational principles is more efficient in producing generalised continuum models. Moreover, if one considers the need of exploiting numerical integration schemes for getting predictions by using the introduced mathematical models, it is clear that a great economy of thought is gained by starting the modelling procedure basing it on variational principles. It is remarkable that in [73, 74], the more general theory of second gradient elasticity is considered together with couple stress one, albeit in the title only the particular case of couple stress models is evoked (and the reader is explicitly warned about this circumstance).

Another example of incomplete second gradient models, which already was mentioned above, is the Korteweg fluid with \(W=W(\rho ,\nabla \rho )\), where \(\rho \) is a mass density [8, 29]. The cited examples are not exhausting the possibilities explored up to now. Indeed, as already mentioned, it has been proven that also some materials having a beam-lattice micro-structure can be modelled using gradient incomplete models [2, 7, 22, 24, 26, 37, 62], for pantographic beam lattice microstructured materials see [30, 32]. In all aforementioned cases of incomplete second gradient models, the analysis of the well posedness of boundary value problems requires a proper modification of the standard techniques [18, 33, 39]. We believe that the correct framework for studying, in an unitary way, many (if not all) incomplete second gradient models is given by the theory of anisotropic Sobolev’s spaces as formulated by Nikol’skii [55] and other nonclassical Sobolev’s spaces. Most likely, Nikol’skii fundamental contributions supply a general and powerful tool for studying the mathematical problems related to the more general class of *N*-th gradient micro-structured continua [25].

The aim of this paper is to prove the uniqueness and existence of weak solutions when considering the particular case of incomplete second gradient model given by the dilatational strain gradient continuum. Within this specific model, the strain energy density depends on strains and the gradient of dilatation, i.e. on the gradient of the trace of the strain tensor. Clearly, the model is closely related to the model of Korteweg or Cahn–Hilliard fluids [8, 29, 67] and could be treated as a gradient extension of generalised models with one-dimensional microstructure [17] or of the damage and poroelasticity models presented in [6, 19, 41, 44, 45, 60, 60, 63]. It is also a particular case of the model [80] where the dependence on the gradient of deviatoric part of strains is neglected. Moreover, the model studied in the present paper could be useful also when modelling the behaviour of porous media undergoing small deformations including ones used for 3D printing [79].

The paper is organised as follows. After the short mathematical preliminaries needed to set the used functional spaces, which are given in Sect. 2, we consider the constitutive equation for deformation energy and external interactions, and consequently, we formulate the principle of virtual work used to determine equilibrium configurations (Sect. 3). The class of external interactions which can be applied to dilatational strain gradient continua were formally determined in [23]. In the present paper, we give the rigorous proof of the fact that dilatational strain gradient continua can support forces per unit area but also double forces per unit area at their (Lipschitz continuous) boundaries.

In Sect. 4, based on the formulated principle of virtual work, we introduce the definition of weak solution that we use, and based on our listed assumptions, we prove that such a solution exists and is unique. Moreover, we prove that it is a minimiser of the total energy functional. Also in Sect. 4, we remark that the two Euler laws of dynamics are necessary but not sufficient conditions for the equilibrium of dilatational strain gradient continua. This remark naturally implies that variational principles seem to be more suitable to formulate novel mathematical models in applied science, see, e.g. [27, 36, 70, 71].

Some conclusions are presented in the final section, together with some future research perspectives.

## Mathematical preliminaries

Let *V* be a bounded connected domain in \({\mathbb {R}}^n\), \(n=2, 3\), with \(C^1\)-regular boundary \(S\equiv \partial V\). For other possible requirements to *S*, we refer to [3]. In what follows, we use Sobolev’s spaces \(H^k(V)\), \(k=1, 2\), defined as follows

where \(L^2(V)\) is the Lebesgue space of square integrable functions, and \(\nabla \) is the spatial nabla operator.

We denote vectors and tensors of any order using semibold font. So for simplicity, we use the notation for vector- or tensor-valued functions \({\mathbf {u}}\in {\mathbf {L}}^2(V)\), etc., as an equivalent to \((L^2(V))^n\), that means that each Cartesian component of \({\mathbf {u}}\equiv (u_1,\ldots ,u_n)\) belongs to \(L^2(V)\).

In addition following Jacques Louis Lions [72, p. 100], we introduce the functional space

where \(\nabla \cdot {\mathbf {u}}\) is the divergence of \({\mathbf {u}}\). Obviously, \({\mathbf {H}}^1(V)\subset {\mathbf {H}}(\mathrm {div},V)\). Considering vector-functions \({\mathbf {u}}\in {\mathbf {H}}^2(V)\), we also introduce the space

Obviously, \({\mathbf {H}}^1(\mathrm {div},V)\) is a separable Hilbert space with the inner product

and the norm

where “\(\cdot \)” and “ : ” stands for the dot and double dot products, respectively. Using Tartar’s example [72, p. 101], we can see that \({\mathbf {u}}=(f_1(x_1)g_1(x_2,x_3), f_2(x_2)g_2(x_1,x_3),f_3(x_3)g_3(x_1,x_3)\) belongs to \({\mathbf {H}}^1(\mathrm {div},V)\) if \(f_i\in H^2(V)\) and \(g_i\in H^1(V)\). As for Tartar’s example given for \({\mathbf {H}}(\mathrm {div},V)\), this example gives the possibility to understand which kinds of discontinuity are possible for functions belonging to \({\mathbf {H}}^1(\mathrm {div},V)\). For further properties of these spaces, we refer to [3, 18, 46, 72].

## Strain energy density and the principle of virtual work

For infinitesimal deformations, we introduce a vector of displacements

where \({\mathbf {x}}\) is a position vector. The corresponding strain tensor and its trace, i.e. dilatation, are given by

where *T* stands for a transpose tensor.

Within the dilatational strain gradient elasticity, there exists a strain energy density *W* which depends on \({\mathbf {e}}\) and \(\nabla \theta \):

Considering quadratic form of *W*, we came to

where \({\mathbf {C}}\), \({\mathbf {B}}\), and \({\mathbf {D}}\) are fourth-order tensors of elastic moduli. In what follows, we restrict ourselves to isotropic materials, so *W* takes the form

where \(\lambda \) and \(\mu \) are the Lamé moduli, \(\lambda >-2/3\mu \), \(\mu >0\), and \(\alpha >0\) is an additional elastic modulus responsible for strain gradient effects.

The equilibrium equations and the corresponding natural boundary conditions within the considered model can be derived using the Lagrange variational principle modified for the strain gradient media as in [1, 8, 29]. Let \(F({\mathbf {u}})\) be the functional of the strain energy

The first variation of *F* takes the form

where \({\mathbf {v}}=\delta {\mathbf {u}}\) and \(\delta \) is the variation symbol. After integration by parts, \(\delta F\) can be transformed into

where \({\mathbf {n}}\) is the vector of unit outward normal to *S*, \(\frac{\partial }{\partial n}\) means the normal derivative, \(\nabla _s=({\mathbf {I}}-{\mathbf {n}}\otimes {\mathbf {n}})\cdot \nabla \) is the surface nabla-operator, \({\mathbf {I}}\) is the unit tensor, “\(\otimes \)” means dyadic product, \(H=-1/2\nabla _s\cdot {\mathbf {n}}\) is the mean curvature of *S*, and the surface divergence theorem was used as in [31], see details in Appendix. In addition, here \({\mathbf {T}}\), \(\varvec{\sigma }\), and \({\mathbf {m}}\) are the total stress, stress tensors, and double force vector given by the relations

In what follows, we assume kinematic boundary conditions on a part \(S_0\) of *S*

which describe the case when \(S_0\) is fixed. So the external loads can be applied on the rest of *S* that is on \(S_1=S/S_0\).

The form of the first variation (3.4) dictates an admissible form of the work of external forces and double forces

where \({\mathbf {f}}\) is a volume force vector, \({\mathbf {t}}\) is a traction vector, *c* is a surface normal double force.

As a result, we formulate the principle of virtual work

for all admissible functions \({\mathbf {v}}\).

Let us consider (3.7) for rigid body motion, i.e. \({\mathbf {v}}\) given by

where \({\mathbf {a}}\) and \({\mathbf {b}}\) are constant vectors and \(\times \) is the cross product. Let us recall that for such \({\mathbf {v}}\) from (3.7) can obtain the conditions of equilibrium of a free solid body (\(S_1=S\)) in classic linear elasticity [33]. These conditions say that the total force and the total torque must be zero. Substituting (3.8) into (3.7), we get

So we have again the total force and the total torque as zero

Note that double forces are not included in these balance equations. For derivation, we used the identities

Just calculated Eqs. (3.10) and (3.11) express the two Euler laws for dilatational strain gradient continua. They are necessary conditions for equilibrium. However, one has to remark that they cannot be sufficient to equilibrium. In facts, in the following section, we prove the existence and uniqueness of weak solutions for the equilibrium elastic problem when surface forces, volume forces and surface double forces are given. As the field of double forces *c* does not appear in found Euler laws, for a given system of balanced external forces, there exists a different equilibrium solution for every different choice of the field *c*. An important consequence of this circumstance seems now evident: the balance laws of forces and torques (couples) are not sufficient to characterise the equilibrium configurations of considered continua. Therefore, postulation schemes based *only* on these two basic laws cannot give the right conceptual schemes for studying them.

## Weak solutions and their properties

From (3.7), one can derive the equilibrium equations and the natural boundary conditions as given in Appendix, see also, e.g. [1, 8, 29] for more details. Here, instead, we use (3.7) as the principal equation for determination of weak solutions of the considered problem without using its strong formulation. It has to be remarked that numerical integration schemes calculate approximations of weak solutions and that strong formulation does not play any role in the context of numerical computations. The importance of the integration by part process leading to strong equilibrium conditions resides in the need of determining which externally applied interactions (i.e. forces, moments, double forces, etc.) can be sustained by a continuum whose internal structure is characterised by a certain deformation energy.

Following the classic approach to the analysis of weak solutions, we introduce the natural functional space for equation (3.7) called energy space. For \(C^2(V)\) functions, we introduce the inner product as follows

Let us note that we treat all quantities as dimensionless ones. Equation (4.1) produces the energy norm

Now, we introduce

### Definition 4.1

The energy space *E* is the closure of \(C^2(V)\) functions satisfying (3.5) in the norm \(\Vert \cdot \Vert _E\).

Obviously, *E* is a separable Hilbert space. Moreover, using the Korn inequality, we can prove that *E* is identical to \({\mathbf {H}}^1(\mathrm {div},V)\). Indeed, we have the Korn inequality [18]

with a positive constant *C*. So the norm \(\Vert \cdot \Vert _E\) is equivalent to \(\Vert \cdot \Vert _{{\mathbf {H}}^1(\mathrm {div},V)}\) as these norms coincide up to positive factors. Thus, we can use properties of \({\mathbf {H}}^1(\mathrm {div},V)\) for characterisation of weak solutions.

### Definition 4.2

We call \({\mathbf {u}}_0\in {\mathbf {H}}^1(\mathrm {div},V)\) a weak solution of the boundary value problem under consideration if it satisfies boundary conditions (3.5) and the relation

where

for all \({\mathbf {v}}\in {\mathbf {H}}^1(\mathrm {div},V)\) also satisfying (3.5).

With this definition, we formulate the main theorem.

### Theorem 4.3

Let \(V\subset {\mathbb {R}}^3\) be a bounded connected domain with \(C^1\)-regular boundary \(S=S_0\cup S_1\), \({\mathbf {f}}\in {\mathbf {L}}^{6/5}(V)\), \({\mathbf {t}}\in {\mathbf {L}}^{4/3}(S_1)\), and \(c\in H^{1}(S_1)\). Then, there exists a weak solution \({\mathbf {u}}_0\in {\mathbf {H}}^1(\mathrm {div},V)\). It is unique. Furthermore,

where

### Proof

The key point of the proof is to show that \(L({\mathbf {v}})\) is a linear bounded functional in *E*. Then, using the Riesz representation theorem, it follows that there is an unique element \(\varvec{\ell }\in E\) such that

As a result, Eq. (4.3) takes the form

so we get the unique solution \({\mathbf {u}}_0=\varvec{\ell }\).

As \(C^2(V)\) is dense in *E* in what follows, we can use functions in \(E\cap C^2(V)\). The necessary properties of \(L({\mathbf {v}})\) follow from the Sobolev imbedding theorems. First, let us consider

As \({\mathbf {H}}^1(V)\subset E\), and the imbedding operators from \(H^1(V)\) to \(L^6(V)\) and \(L^4(S_1)\) are continuous, we find that

where Hölder’s inequalities were used, and *C* stands for a positive constant depending on \({\mathbf {f}}\) and \({\mathbf {t}}\).

Now, let us consider

By definition of \({\mathbf {H}}^1(\mathrm {div},V)\), we have that if \({\mathbf {v}}\in {\mathbf {H}}^1(\mathrm {div},V)\), then \(\theta ({\mathbf {v}})=\nabla \cdot {\mathbf {v}}\) belongs to \(H^1(V)\). So \(\theta \) has a trace on \(S_1\) and \(\theta \big |_{S_1}\in L^4(S_1)\) and \(\theta \big |_{S_1}\in L^2(S_1)\). Obviously, \({\mathbf {v}}|_{S_1}\in {\mathbf {L}}^2(S_1)\). As \(\theta = \frac{\partial {\mathbf {v}}}{\partial n}\cdot {\mathbf {n}}+\nabla _s\cdot {\mathbf {v}}\), we can transform \(L_2({\mathbf {v}})\) using integration by parts as follows

Here, we take into account that \({\mathbf {v}}|_{\partial S_1}={\mathbf {0}}\) along \(\partial S_1\) as this contour is an interface between \(S_0\) and \(S_1\), whereas \({\mathbf {v}}|_{ S_0}={\mathbf {0}}\). As a result, there is no a contour integral in \(L_2({\mathbf {v}})\).

As a result, with Hölder’s inequality, we get

Thus, \(L({\mathbf {v}})=L_1({\mathbf {v}})+L_2({\mathbf {v}})\) is a linear and bounded functional.

The minimisation property of \({\mathbf {u}}_0\) follows from the fact that \(J({\mathbf {u}})\) is a quadratic functional and Eq. (4.3) states that the first variation of \( J({\mathbf {u}})\) vanishes at \({\mathbf {u}}={\mathbf {u}}_0\), see, e.g. [33] for more details. \(\square \)

From the proof, it follows that the weak solution is bounded

with a positive constant *C* which depends on *V* and \(S_1\). As a solution is unique, \({\mathbf {u}}_0={\mathbf {0}}\) if and only if \({\mathbf {f}}={\mathbf {0}}\), \({\mathbf {t}}={\mathbf {0}}\), and \(c=0\). In particular, if \({\mathbf {f}}={\mathbf {0}}\), \({\mathbf {t}}={\mathbf {0}}\), but \(c\ne 0\) we have \({\mathbf {u}}_0\ne {\mathbf {0}}\).

### Remark 1

If \(V\subset {\mathbb {R}}^2\), then the theorem is also valid for \({\mathbf {f}}\in {\mathbf {L}}^p(V)\) and \({\mathbf {t}}\in {\mathbf {L}}^q(V)\) for \(1<p,q<\infty \).

### Remark 2

In the proof, we use the same assumptions on \({\mathbf {f}}\) and \({\mathbf {t}}\) as for the linear elasticity, see [18, 33], and even more strong assumption was applied to *c*. In fact, using dual spaces, one can consider more weak assumptions, such as \({\mathbf {t}}\in {\mathbf {H}}^{-1/2}(S_1)\), \(c\in H^{1/2}(S_1)\), see [14, 46] for trace properties in \(H(\mathrm {div},V)\). We leave this analysis to forthcoming papers.

### Remark 3

Instead of \(C^1\)-regular boundary, one can consider less regular surface for which we have required Sobolev’s imbedding theorems [3]. We leave this again to forthcoming papers.

## Conclusions

We proved the existence and uniqueness of a weak solution of a specific class of equilibrium problems within the framework of newly proposed strain gradient elasticity model specified with the adjective “*dilatational*”. The external interactions applied to the considered dilatational strain gradient continuum are dead loads constituted by (i) forces per unit area and (ii) purely normal double forces per unit area applied at the boundary of the considered continuum, which is occupying a regular domain, and (iii) volume forces applied at every material point belonging to the same continuum. This result has a conceptual impact as it proves that a consistent mathematical problem can be formulated in which applied external dead loads cannot be reduced to forces per unit area and couples per unit area, reduction that is assumed, instead, in [76]. This means that the postulation scheme described there does not include all conceivable and logically consistent continuum models and that there are continuum models, based on variational principles, which cannot be obtained by postulating the balance of force and couples for every subbody of considered continuum. The statement of the theorems presented in Sect. 4 gives a solid mathematical ground to the discussion presented by Sedov [64,65,66] and should be sufficient to conclude the debate between the variational and “balancist” schools in continuum mechanics by establishing the superiority of variational principle postulation schemes. One has to remark explicitly here that the concept of stress, as developed by “balancists”, is specific for first gradient continua and cannot be easily generalised to higher gradient continua, see, e.g. [23]. In fact, assuming only balance of forces and couples for any sub-body, we necessarily come to the Cosserat continuum [34, 56] or to the Cauchy continuum if we neglect couples as independent on forces. For media with microstructure, one has to introduce additional balance equations, see [17, 35]. Instead, for hyperelastic, the concept of deformation energy can be more easily generalised and leads easily, via a process based on stationarity conditions and integration by parts, to the strong formulations of mathematically consistent equilibrium problems which naturally include well-posed natural boundary conditions. Moreover, the concept of deformation energy is based on the fundamental concept of deformation, which is pure kinematical and is strongly linked to experimental evidence. Instead, the concept of stress has a complex mathematical nature: in fact, it can be defined as the linear and continuous functional mapping displacement fields and their gradients into the expended work, see Eq. (3.3). It is rather difficult to measure stress without having developed a dynamical theory, and these measurements are always based on the indirect determination based on the direct measure of deformations and the use of ad hoc postulated constitutive equations. It seems to us rather difficult to attribute a “more intuitive” physical nature to stress than to deformation.

In the specific considered modelling instance, we consider a very specific case of the energy space adapted to the postulated deformation energy. We prove that weak solutions to the equilibrium problem belong to \({\mathbf {H}}^1(\mathrm {div},V)\). This is an intermediate functional space between \({\mathbf {H}}^1(V)\) and \({\mathbf {H}}^2(V)\) and belongs to the wider class of anisotropic Sobolev spaces introduced by Nikol’skii [55], see also [10,11,12, 57, 75] for a study about traces of functions in these spaces.

This characterisation of the solutions is essential if one wants to solve with numerical methods the problem of the deformation of dilatational strain gradient continua under the specified class of externally applied dead loads. In facts, in order to get a reasonable accurate approximation and a convergent integration scheme, one has to suppose that the applied loads have a regularity compatible with formula (3.7) and that the used mixed finite elements method exploits a discretisation based on a set of test functions dense in \({\mathbf {H}}^1(\mathrm {div},V)\).

## References

- 1.
Abali, B.E., Müller, W.H., dell’Isola, F.: Theory and computation of higher gradient elasticity theories based on action principles. Arch. Appl. Mech.

**87**(9), 1495–1510 (2017) - 2.
Abdoul-Anziz, H., Seppecher, P.: Strain gradient and generalized continua obtained by homogenizing frame lattices. Math. Mech. Complex Syst.

**6**(3), 213–250 (2018) - 3.
Adams, R.A., Fournier, J.J.F.: Sobolev Spaces, Pure and Applied Mathematics, vol. 140, 2nd edn. Academic Press, Amsterdam (2003)

- 4.
Agranovich, M.: Elliptic boundary problems. In: Agranovich, M., Egorov, Y., Shubin, M. (eds.) Partial Differential Equations IX: Elliptic Boundary Problems. Encyclopaedia of Mathematical Sciences, vol. 79, pp. 1–144. Springer, Berlin (1997)

- 5.
Alibert, J.J., Seppecher, P., dell’Isola, F.: Truss modular beams with deformation energy depending on higher displacement gradients. Math. Mech. Solids

**8**(1), 51–73 (2003) - 6.
Andreaus, U., Giorgio, I., Madeo, A.: Modeling of the interaction between bone tissue and resorbable biomaterial as linear elastic materials with voids. Z. Angew. Math. Phys.

**66**(1), 209–237 (2015) - 7.
Andreaus, U., Spagnuolo, M., Lekszycki, T., Eugster, S.R.: A Ritz approach for the static analysis of planar pantographic structures modeled with nonlinear Euler–Bernoulli beams. Cont. Mech. Thermodyn.

**30**(5), 1103–1123 (2018) - 8.
Auffray, N., dell’Isola, F., Eremeyev, V.A., Madeo, A., Rosi, G.: Analytical continuum mechanics à la Hamilton–Piola least action principle for second gradient continua and capillary fluids. Math. Mech. Solids

**20**(4), 375–417 (2015) - 9.
Barchiesi, E., Spagnuolo, M., Placidi, L.: Mechanical metamaterials: a state of the art. Math. Mech. Solids

**24**(1), 212–234 (2019) - 10.
Besov, O.V., II’in, V.P., Nikol’skii, S.M.: Integral Representations of Functions and Imbedding Theorems, vol. 1. Wiley, New York (1978)

- 11.
Besov, O.V., II’in, V.P., Nikol’skii, S.M.: Integral Representations of Functions and Imbedding Theorems, vol. 2. Wiley, New York (1979)

- 12.
Besov, O.V., II’in, V.P., Nikol’skii, S.M.: Integral Representations of Functions and Imbedding Theorems (in Russian). Nauka, Moscow (1996)

- 13.
Boutin, C., Giorgio, I., dell’Isola, F., Placidi, L.: Linear pantographic sheets: asymptotic micro-macro models identification. Math. Mech. Complex Syst.

**5**(2), 127–162 (2017) - 14.
Buffa, A., Ciarlet Jr., P.: On traces for functional spaces related to Maxwell’s equations Part II: Hodge decompositions on the boundary of Lipschitz polyhedra and applications. Math. Methods Appl. Sci.

**24**(1), 31–48 (2001) - 15.
Cahn, J.W., Hilliard, J.E.: Free energy of a nonuniform system. I. Interfacial free energy. J. Chem. Phys.

**28**(2), 258–267 (1958) - 16.
Cahn, J.W., Hilliard, J.E.: Free energy of a nonuniform system. III. Nucleation in a two-component incompressible fluid. J. Chem. Phys.

**31**(3), 688–699 (1959) - 17.
Capriz, G.: Continua with Microstructure. Springer, New York (1989)

- 18.
Ciarlet, P.G.: Mathematical Elasticity. Vol. I: Three-Dimensional Elasticity. North-Holland, Amsterdam (1988)

- 19.
Contrafatto, L., Cuomo, M., Greco, L.: Meso-scale simulation of concrete multiaxial behaviour. Eur. J. Environ. Civ. Eng.

**21**(7–8), 896–911 (2017) - 20.
Cordero, N.M., Forest, S., Busso, E.P.: Second strain gradient elasticity of nano-objects. J. Mech. Phys. Solids

**97**, 92–124 (2016) - 21.
dell’Isola, F., Della Corte, A., Giorgio, I.: Higher-gradient continua: the legacy of Piola, Mindlin, Sedov and Toupin and some future research perspectives. Math. Mech. Solids

**22**(4), 852–872 (2017) - 22.
dell’Isola, F., Giorgio, I., Pawlikowski, M., Rizzi, N.: Large deformations of planar extensible beams and pantographic lattices: heuristic homogenisation, experimental and numerical examples of equilibrium. Proc. R. Soc. Lond. Ser. A

**472**(2185), 20150,790 (2016) - 23.
dell’Isola, F., Seppecher, P.: Edge contact forces and quasi-balanced power. Meccanica

**32**(1), 33–52 (1997) - 24.
dell’Isola, F., Seppecher, P., Alibert, J.J., Lekszycki, T., Grygoruk, R., Pawlikowski, M., Steigmann, D., Giorgio, I., Andreaus, U., Turco, E., Gołaszewski, M., Rizzi, N., Boutin, C., Eremeyev, V.A., Misra, A., Placidi, L., Barchiesi, E., Greco, L., Cuomo, M., Cazzani, A., Corte, A.D., Battista, A., Scerrato, D., Eremeeva, I.Z., Rahali, Y., Ganghoffer, J.F., Müller, W., Ganzosch, G., Spagnuolo, M., Pfaff, A., Barcz, K., Hoschke, K., Neggers, J., Hild, F.: Pantographic metamaterials: an example of mathematically driven design and of its technological challenges. Cont. Mech. Thermodyn.

**31**(4), 851–884 (2019) - 25.
dell’Isola, F., Seppecher, P., Madeo, A.: How contact interactions may depend on the shape of Cauchy cuts in \(N\)th gradient continua: approach “à la D’Alembert”. Z. Angew. Math. Phys.

**63**(6), 1119–1141 (2012) - 26.
dell’Isola, F., Seppecher, P., Spagnuolo, M., Barchiesi, E., Hild, F., Lekszycki, T., Giorgio, I., Placidi, L., Andreaus, U., Cuomo, M., Eugster, S.R., Pfaff, A., Hoschke, K., Langkemper, R., Turco, E., Sarikaya, R., Misra, A., De Angelo, M., D’Annibale, F., Bouterf, A., Pinelli, X., Misra, A., Desmorat, B., Pawlikowski, M., Dupuy, C., Scerrato, D., Peyre, P., Laudato, M., Manzari, L., Göransson, P., Hesch, C., Hesch, S., Franciosi, P., Dirrenberger, J., Maurin, F., Vangelatos, Z., Grigoropoulos, C., Melissinaki, V., Farsari, M., Muller, W., Abali, B.E., Liebold, C., Ganzosch, G., Harrison, P., Drobnicki, R., Igumnov, L., Alzahrani, F., Hayat, T.: Advances in pantographic structures: design, manufacturing, models, experiments and image analyses. Cont. Mech. Thermodyn.

**31**(4), 1231–1282 (2019) - 27.
dell’Isola, F., Steigmann, D.J.: Discrete and Continuum Models for Complex Metamaterials. Cambridge University Press, Cambridge (2020)

- 28.
Epstein, M., Smelser, R.: An appreciation and discussion of Paul Germain’s “The method of virtual power in the mechanics of continuous media, I: second-gradient theory”. Math. Mech. Complex Syst.

**8**(2), 191–199 (2020) - 29.
Eremeyev, V.A., Altenbach, H.: Equilibrium of a second-gradient fluid and an elastic solid with surface stresses. Meccanica

**49**(11), 2635–2643 (2014) - 30.
Eremeyev, V.A., Alzahrani, F.S., Cazzani, A., dell’Isola, F., Hayat, T., Turco, E., Konopińska-Zmysłowska, V.: On existence and uniqueness of weak solutions for linear pantographic beam lattices models. Cont. Mech. Thermodyn.

**31**(6), 1843–1861 (2019) - 31.
Eremeyev, V.A., Cloud, M.J., Lebedev, L.P.: Applications of Tensor Analysis in Continuum Mechanics. World Scientific, NJ (2018)

- 32.
Eremeyev, V.A., dell’Isola, F., Boutin, C., Steigmann, D.: Linear pantographic sheets: existence and uniqueness of weak solutions. J. Elast.

**132**(2), 175–196 (2018) - 33.
Eremeyev, V.A., Lebedev, L.P.: Existence of weak solutions in elasticity. Math. Mech. Solids

**18**(2), 204–217 (2013) - 34.
Eremeyev, V.A., Lebedev, L.P., Altenbach, H.: Foundations of Micropolar Mechanics. Springer, Heidelberg (2013)

- 35.
Eringen, A.C.: Microcontinuum Field Theory. I. Foundations and Solids. Springer, New York (1999)

- 36.
Eugster, S.R., Steigmann, D.J.: Variational methods in the theory of beams and lattices. In: Altenbach H., Öchsner A. (eds) Encyclopedia of Continuum Mechanics, pp. 1–9. Springer, Berlin. https://doi.org/10.1007/978-3-662-55771-6_176 (2020)

- 37.
Eugster, S., Steigmann, D., et al.: Continuum theory for mechanical metamaterials with a cubic lattice substructure. Math. Mech. Complex Syst.

**7**(1), 75–98 (2019) - 38.
Fichera, G.: Linear Elliptic Differential Systems and Eigenvalue Problems. Lecture Notes in Mathematics, vol. 8. Springer, Berlin (1965)

- 39.
Fichera, G.: Existence theorems in elasticity. In: Flügge, S. (ed.) Handbuch der Physik, vol. VIa/2, pp. 347–389. Springer, Berlin (1972)

- 40.
Forest, S., Cordero, N.M., Busso, E.P.: First vs. second gradient of strain theory for capillarity effects in an elastic fluid at small length scales. Comput. Mater. Sci.

**50**(4), 1299–1304 (2011) - 41.
Gagneux, G., Millet, O.: Modeling capillary hysteresis in unsatured porous media. Math. Mech. Complex Syst.

**4**(1), 67–77 (2016) - 42.
Germain, P.: Functional concepts in continuum mechanics. Meccanica

**33**(5), 433–444 (1998) - 43.
Germain, P.: The method of virtual power in the mechanics of continuous media, I: Second-gradient theory. Math. Mech. Complex Syst.

**8**(2), 153–190 (2020) - 44.
Giorgio, I., Andreaus, U., Scerrato, D., Braidotti, P.: Modeling of a non-local stimulus for bone remodeling process under cyclic load: application to a dental implant using a bioresorbable porous material. Math. Mech. Solids

**22**(9), 1790–1805 (2017) - 45.
Giorgio, I., Scerrato, D.: Multi-scale concrete model with rate-dependent internal friction. Eur. J. Environ. Civ. Eng.

**21**(7–8), 821–839 (2017) - 46.
Girault, V., Raviart, P.A.: Finite Element Methods for Navier–Stokes Equations: Theory and Algorithms. Springer, Berlin (1986)

- 47.
Healey, T.J., Krömer, S.: Injective weak solutions in second-gradient nonlinear elasticity. ESAIM: Control, Optim. Calc. Var.

**15**(4), 863–871 (2009) - 48.
Hörmander, L.: The Analysis of Linear Partial Differential Operators. II. Differential Operators with Constant Coefficients, A Series of Comprehensive Studies in Mathematics, vol. 257. Springer, Berlin (1983)

- 49.
Korteweg, D.J.: Sur la forme que prennent les équations des mouvements des fluides si l’on tient compte des forces capillaires par des variations de densité. Arch. Néerl. Sci. Exactes Nat. Sér.

**II**(6), 1–24 (1901) - 50.
Mareno, A., Healey, T.J.: Global continuation in second-gradient nonlinear elasticity. SIAM J. Math. Anal.

**38**(1), 103–115 (2006) - 51.
Maugin, G.A.: Generalized continuum mechanics: various paths. In: Continuum Mechanics Through the Twentieth Century: A Concise Historical Perspective, pp. 223–241. Springer, Dordrecht (2013)

- 52.
Maugin, G.A.: Non-Classical Continuum Mechanics: A Dictionary. Springer, Singapore (2017)

- 53.
Mindlin, R.D.: Micro-structure in linear elasticity. Arch. Ration. Mech. Anal.

**16**(1), 51–78 (1964) - 54.
Mindlin, R.D., Eshel, N.N.: On first strain-gradient theories in linear elasticity. Int. J. Solids Struct.

**4**(1), 109–124 (1968) - 55.
Nikol’skii, S.M.: On imbedding, continuation and approximation theorems for differentiable functions of several variables. Rus. Math. Surv.

**16**(5), 55 (1961) - 56.
Nowacki, W.: Theory of Asymmetric Elasticity. Pergamon-Press, Oxford (1986)

- 57.
Ohno, M., Shizuta, Y., Yanagisawa, T.: The trace theorem on anisotropic Sobolev spaces. Tohoku Math. J., Second Ser.

**46**(3), 393–401 (1994) - 58.
Pideri, C., Seppecher, P.: A homogenization result for elastic material reinforced periodically with high rigidity elastic fibres. Comptes Rendus l’Acad. Sci. Ser. IIB Mech. Phys. Chem. Astron.

**8**(324), 475–481 (1997) - 59.
Pideri, C., Seppecher, P.: A second gradient material resulting from the homogenization of an heterogeneous linear elastic medium. Cont. Mech. Thermodyn.

**9**(5), 241–257 (1997) - 60.
Placidi, L., Barchiesi, E., Misra, A.: A strain gradient variational approach to damage: a comparison with damage gradient models and numerical results. Math. Mech. Complex Syst.

**6**(2), 77–100 (2018) - 61.
Placidi, L., Barchiesi, E., Turco, E., Rizzi, N.L.: A review on 2D models for the description of pantographic fabrics. Z. Angew. Math. Phys.

**67**(5), 121 (2016) - 62.
Rahali, Y., Giorgio, I., Ganghoffer, J.F., dell’Isola, F.: Homogenization à la Piola produces second gradient continuum models for linear pantographic lattices. Int. J. Eng. Sci.

**97**, 148–172 (2015) - 63.
Sciarra, G., dell’Isola, F., Coussy, O.: Second gradient poromechanics. Int. J. Solids Struct.

**44**(20), 6607–6629 (2007) - 64.
Sedov, L.I.: Mathematical methods for constructing new models of continuous media. Rus. Math. Surv.

**20**(5), 123 (1965) - 65.
Sedov, L.I.: Models of continuous media with internal degrees of freedom: PMM vol. 32, n 5, 1968, pp. 771–785. J. Appl. Math. Mech.

**32**(5), 803–819 (1968) - 66.
Sedov, L.I.: Variational methods of constructing models of continuous media. In: Parkus, H., Sedov, L.I. (eds.) Irreversible Aspects of Continuum Mechanics and Transfer of Physical Characteristics in Moving Fluids, pp. 346–358. Springer, Vienna (1968)

- 67.
Seppecher, P.: Second-gradient theory: application to Cahn–Hilliard fluids. In: Maugin G.A., Drouot R., Sidoroff F. (eds) Continuum Thermomechanics. Solid Mechanics and Its Applications, vol 76, pp. 379–388. Springer, Dordrecht (2000)

- 68.
Soubestre, J., Boutin, C.: Non-local dynamic behavior of linear fiber reinforced materials. Mech. Mater.

**55**, 16–32 (2012) - 69.
Spagnuolo, M., Yildizdag, M.E., Andreaus, U., Cazzani, A.M.: Are higher-gradient models also capable of predicting mechanical behavior in the case of wide-knit pantographic structures? Math. Mech. Solids. https://doi.org/10.1177/1081286520937339 (2020)

- 70.
Steigmann, D.J.: The variational structure of a nonlinear theory for spatial lattices. Meccanica

**31**(4), 441–455 (1996) - 71.
Steigmann, D., Baesu, E., Rudd, R.E., Belak, J., McElfresh, M.: On the variational theory of cell-membrane equilibria. Interfaces Free Bound.

**5**(4), 357–366 (2003) - 72.
Tartar, L.: An Introduction to Sobolev Spaces and Interpolation Spaces. Springer, Berlin (2007)

- 73.
Toupin, R.A.: Elastic materials with couple-stresses. Arch. Ration. Mech. Anal.

**11**(1), 385–414 (1962) - 74.
Toupin, R.A.: Theories of elasticity with couple-stress. Arch. Ration. Mech. Anal.

**17**(2), 85–112 (1964) - 75.
Triebel, H.: Theory of Function Spaces III, Monographs in Mathematics, vol. 100. Birkhäuser, Basel (2006)

- 76.
Truesdell, C., Noll, W.: The Non-linear Field Theories of Mechanics, 3rd edn. Springer, Berlin (2004)

- 77.
Volkov-Bogorodsky, D.B., Evtushenko, Y.G., Zubov, V.I., Lurie, S.A.: Calculation of deformations in nanocomposites using the block multipole method with the analytical-numerical account of the scale effects. Comput. Math. Math. Phys.

**46**(7), 1234–1253 (2006) - 78.
Yang, F.A.C.M., Chong, A.C.M., Lam, D.C.C., Tong, P.: Couple stress based strain gradient theory for elasticity. Int. J. Solids Struct.

**39**(10), 2731–2743 (2002) - 79.
Yang, H., Ganzosch, G., Giorgio, I., Abali, B.E.: Material characterization and computations of a polymeric metamaterial with a pantographic substructure. Z. Angew. Math. Phys.

**69**(4), 105 (2018) - 80.
Zhou, S., Li, A., Wang, B.: A reformulation of constitutive relations in the strain gradient elasticity theory for isotropic materials. Int. J. Solids Struct.

**80**, 28–37 (2016)

## Acknowledgements

This work was supported by the Russian Science Foundation under Grant 20-41-04404 issued to the Institute of Applied Mechanics of Russian Academy of Sciences.

## Author information

### Affiliations

### Corresponding author

## Additional information

### Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Appendix: First variation

### Appendix: First variation

Let us discuss the derivation of the first variation given by (3.4). In what follows, for simplicity, we consider \(C^4(V)\) functions. Integrating by parts in (3.3), we get

where \(S=\partial V\) and \({\mathbf {n}}\) is a unit vector of outward normal to *S*. Here, we also used that \({\mathbf {e}}\) is a symmetric tensor, so \({\mathbf {e}}({\mathbf {u}}):{\mathbf {e}}({\mathbf {v}})={\mathbf {e}}({\mathbf {u}}):\nabla {\mathbf {v}}\). Introducing stress tensor \(\varvec{\sigma }\) and double force vector \({\mathbf {m}}\) by

so the strain energy density takes the form

we re-write (5.1) in a compact form

Integrating by part again, we get

Substituting (5.3) into (5.2), we transform \(\delta F({\mathbf {u}};{\mathbf {v}})\) as follows

Using the identity \(\nabla (\nabla \cdot {\mathbf {m}})=\nabla \cdot (\nabla {\mathbf {m}})^T\) and introducing the total stress tensor \({\mathbf {T}}\) as

where *T* stands for a transpose tensor, we re-write (5.4) in more compact form

Let us consider the last surface integral in (5.5). We recall that \(\theta ({\mathbf {v}})=\nabla \cdot {\mathbf {v}}\). We use the following representation for \(\nabla \cdot {\mathbf {v}}\big |_S\)

where \(\nabla _s=({\mathbf {I}}-{\mathbf {n}}\otimes {\mathbf {n}})\cdot {\mathbf {n}}\) and \(\frac{\partial }{\partial n}\) is the normal derivative. So we get

We consider *S* as a smooth enough surface with a contour \(C=\partial S\). Using the surface divergence theorem as in [31], we can integrate by parts in the last integral in (5.6) as follows

where *H* is the mean curvature of *S*, \(\varvec{\nu }\) is the unit vector of the outward normal to *C* such that \(\varvec{\nu }\cdot {\mathbf {n}}=0\).

Summarising, we get

Note that for bounded domain *V*, \(S=\partial V\), so \(C=\partial S=\partial \partial V\) is empty, but we keep it here for the case of mixed boundary conditions when a part of *S* is clamped. Then, *C* is an interface between clamped part of the boundary and the free one. As on a clamped surface, we have \({\mathbf {v}}={\mathbf {0}}\), the integral over *C* also vanishes.

Using standard techniques of calculus of variations from variational equation \(\delta F({\mathbf {u}};{\mathbf {v}})=0\) for all admissible \({\mathbf {v}}\) it follows the homogeneous boundary-value problem for the determination of \({\mathbf {u}}({\mathbf {x}})\in C^4(V)\)

Obviously, using (5.11), we can simplify (5.10). Indeed, substituting (5.11) into (5.10), we get

Boundary-value problem (5.9)–(5.11) can be easily extended to the case of external loadings. Indeed, in this case, we have the equilibrium equation and the natural boundary conditions in the form

where \({\mathbf {f}}\) is a vector of volume forces, \({\mathbf {t}}\) is a traction, and *c* is a scalar surface double force.

Let us note that in this case can exclude \({\mathbf {n}}\cdot {\mathbf {m}}({\mathbf {u}})\) from (5.14). We have

At first look, we can replace boundary conditions (5.14) and (5.15) by more simple equations

with two given functions \({\mathbf {t}}\) and *c*. But here, we cannot treat \({\mathbf {t}}\) and *c* independently. For example, we cannot apply \({\mathbf {t}}={\mathbf {0}}\) with \(c\ne 0\). For example, for a curved surface (\(H\ne 0 \)), a constant double force *c* produces normal pressure \({\tilde{{\mathbf {t}}}}=2H c {\mathbf {n}}\). So one should aware of application of such “simplified” boundary conditions in the case of strain-gradient continua. This situation is similar to the case of Kirchhoff plates when at the boundary a combination of transverse forces and moments, see, e.g. [31].

As kinematical counterparts of (5.10) and (5.11), we consider the following essential boundary conditions

Considering mixed boundary conditions on the base of (5.8) and (5.13)–(5.15), we come to the formulation of the virtual work principle in the form (3.7). In the case of \(C^4\)-functions, (3.7) is equivalent to (5.13)–(5.15).

## Rights and permissions

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

## About this article

### Cite this article

Eremeyev, V.A., Lurie, S.A., Solyaev, Y.O. *et al.* On the well posedness of static boundary value problem within the linear dilatational strain gradient elasticity.
*Z. Angew. Math. Phys.* **71, **182 (2020). https://doi.org/10.1007/s00033-020-01395-5

Received:

Accepted:

Published:

### Keywords

- Strain gradient elasticity
- Weak solution
- Dilatational elasticity
- Sobolev spaces
- \(H(\mathrm{div},V)\)

### Mathematics Subject Classification

- Primary 74A30
- Secondary 74G30
- 74G65
- 46E35