On the Kirchhoff-Love Hypothesis (Revised and Vindicated)

The Kirchhoff-Love hypothesis expresses a kinematic constraint that is assumed to be valid for the deformations of a three-dimensional body when one of its dimensions is much smaller than the other two, as is the case for plates. This hypothesis has a long history checkered with the vicissitudes of life: even its paternity has been questioned, and recent rigorous dimension-reduction tools (based on standard Γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\varGamma $\end{document}-convergence) have proven to be incompatible with it. We find that an appropriately revised version of the Kirchhoff-Love hypothesis is a valuable means to derive a two-dimensional variational model for elastic plates from a three-dimensional nonlinear free-energy functional. The bending energies thus obtained for a number of materials also show to contain measures of stretching of the plate’s mid surface (alongside the expected measures of bending). The incompatibility with standard Γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\varGamma $\end{document}-convergence also appears to be removed in the cases where contact with that method and ours can be made.

1. Points of the plate lying initially on a normal to the middle plane of the plate remain on the normal to the middle surface of the plate after bending. 2. The distance of every point of the plate from the middle surface remains unchanged by the deformation.
For mysterious reasons, these statements have become to be known as the Kirchhoff-Love hypothesis; we shall stick to this not fully justified tradition. This hypothesis has been variously criticized in the literature, mainly for the inconsistencies that it may cause with the distribution of stresses that (depending on the specific constitutive law) should sustain the assumed deformation. Podio-Guidugli [63,64] overcame this criticism by taking the view that the Kirchhoff-Love hypothesis is a constraint on the admissible deformations, which is sustained by an appropriate reactive stress, specified in a class admitted by symmetry and determined by the equilibrium equations of three-dimensional elasticity.
More recently, the rigorous analytical tool of Γ -convergence has been employed to derive plate (and shell) theories from three-dimensional elasticity. Admittedly, the problem with this method is that in general, apart from noticeable exceptions [2,7], it only affords to derive single powers of the energy expansion in the thickness 2h. Said differently, we can obtain from a three-dimensional constitutive model either a "membrane-dominated model" or a "flexural-dominated model" (in the words of [14]), meaning that we can isolate twodimensional energies either linear or cubic in h, respectively. For example, models in the former category have been derived in [6] for linear plates and in [44] for nonlinear ones (as well as in [45] for nonlinear shells). Models in the latter category have been derived in [22,24] for nonlinear plates (as well as in [23] for nonlinear shells). These higher-order Γlimits, however, need to be evaluated on the class of deformations that minimize the lower order. In other words, we may only recover the h-cubic bending energy on the minimizers of the h-linear stretching energy. 1 Sadly, the analytic notion of Γ -limit has not yet fully evolved into that of Γ -expansion, and thus it does not yet serve the purpose of deriving blended stretching and bending energies, free to conspire together in a thin sheet of a (possibly activable) elastic material, which is our objective here.
In the engineering literature, the role of Γ -convergence as a method to validate reduced elastic theories for structure mechanics has further flourished in recent years. As lucidly reviewed in [65] and [66, § 2], Γ -convergence has been employed to justify both the Kirchhoff-Love and Reissner-Mindlin theories for linearly elastic plates, in a variety of ways, which have been classified under two general categories, standard [60,61] and improved [59]. 2 For nonlinear elastic plates, standard Γ -convergence clashed with the Kirchhoff-Love hypothesis. It was proven in [22] that the rigorous bending energy (in the flexural-dominated model) is incompatible with the deformation field assumed by the Kirchhoff-Love hypothesis. Here, we try and remedy this clash by revising (and salvaging) the classical hypothesis. In particular, we shall see that the incompatibility pointed out in [22] is resolved by our revised hypothesis.
The paper is organized as follows. In Sect. 2, we recall the basic kinematics of plates and present our revision of the classical Kirchhoff-Love hypothesis. In Sect. 3, we introduce a description for the deformation of a smooth surface that relies on a notion of Cartesian connectors, which avoid the use of coordinates and Christoffel symbols. Sections 4 and 5 are devoted to the dimension-reduction afforded by our revised kinematic hypothesis in two distinct constitutive classes, one for incompressible plates and the other for compressible ones. We consider a number of special nonlinear elastic models for the application of our method; for all we derive stretching and bending energies. The features that these results have in common is the presence in the bending energy of stretching measures of the mid surface of the plate (alongside the expected measures of bending). In Sect. 6, we summarize our conclusions and comment on some possible avenues along which this work could be extended. This paper is closed by the Appendix, where we give explicit formulae for the mean and Gaussian curvatures of the deformed mid surface in terms of the mapping that describes it.

Kinematics of Plates
Here we wish to describe the deformation of an elastic plate with a uniform width (and a planar reference configuration). The deformation will be split into two components, a planar one, which maps the reference plane mid surface onto a deformed mid surface, and an axial one, which maps vectors normal to the mid reference surface on vectors normal to the deformed mid surface. The classical Kirchhoff-Love hypothesis consists in assuming that the second mapping is an isometry (see, for example, [76, p. 551] and [12, p. 156]). Because of this isometry, assuming regularity for the planar mapping is enough to ensure an admissible deformation for sufficiently thin plates. In the following, the latter assumption will be made precise and the isometric constraint along normals will be relaxed. In this framework, an approximate right Cauchy-Green tensor will be constructed and its invariants computed.

Kinematic Preliminaries
Let S be a bounded, two-dimensional flat domain immersed in three-dimensional Euclidean space E and h > 0 a real constant. We call y : S → E an injective C 3 -immersion of S and we denote its image as S := y(S). Pursuing our aim of extending the Kirchhoff-Love hypothesis, we interpret the closed set S := S × [−h, h] ⊂ E as the reference configuration of an elastic plate whose mid surface is S. As we focus on the case 2h diam S, we set diam S = 1, for simplicity, meaning that we shall rescale all lengths to diam S.
We define the mapping f : where ν is the unit normal vector to S and φ : S → R is a C 2 -function which describes how normals to S deform into normals to S . It follows from (1) that the deformation gradient reads as where ∇ denotes the gradient in x, a prime denotes differentiation with respect to x 3 , and e 3 is the unit normal to S. 3 Furthermore, φ is assumed to obey which is justified by the requirement that f be orientation-preserving at least for x 3 = 0, as there, by (2), The classical Kirchhoff-Love hypothesis just requires that φ ≡ x 3 , and so it trivially complies with (3). In our approach, φ will rather remain free and either determined to enforce the constraint of bulk incompressibility or used to minimize the elastic energy stored across the thickness of the deformed plate. We see now how (3) can ensure that f is an orientationpreserving C 2 -diffeomorphism onto its image, for appropriately small values of h.
Proof In the following two steps, we adapt Theorem 4.1-1 of [12, p. 157] to our setting.
1. By the continuity of F, inequality (4) implies that Since S is compact, η attains it minimum in S. Moreover, the minimum of η over S must be strictly positive, otherwise y would fail to be an injective immersion. Thus, h can be chosen so that which is where we shall hereafter take it to be.
In Sects. 4 and 5, we shall use a polynomial approximation for φ in computing the invariants of the right Cauchy-Green tensor C f associated with f . Now, we justify this approximation and lay down a number of preliminary formulae for the invariants of C f .

Invariants of C f
We learned in Proposition 1 how to choose h > 0 sufficiently small so that the mapping f is a C 2 -diffeomorphism. Hypothesis (3) also implies that ∇φ(x, 0) = 0 ∀ x ∈ S; thus, it is also possible to choose h so small that This inequality will be assumed to be valid in the following, and h will be taken to comply with both (5) and (6). Within the approximation stated in (6), F in (2) will be written as Definition 1 The corresponding right Cauchy-Green tensor C f associated with f is given by where and C := (∇y) T (∇y), Here C is the right Cauchy-Green tensor associated with the deformation y, while is the left Cauchy-Green tensor associated with the same deformation. The reader should heed that all tensors C, C 1 , and C 2 act on the two-dimensional space , where V is the translation space associated with three-dimensional Euclidean space E . B(x), however, at the place y(x) ∈ S , acts on the two-dimensional space V ν := {v ∈ V : v · ν = 0}.

Remark 1
The curvature tensor ∇ s ν of S , where ∇ s denotes the surface gradient on S , is a symmetric tensor on V ν (see, for example, [27]). We easily see that both tensors C 1 and C 2 can be expressed in terms of ∇ s ν. As ∇ s ν = (∇ν)(∇y) −1 , it readily follows from (8c) and the symmetry of ∇ s ν that We now compute the principal invariants of C f .

Proposition 2
The first invariant, I 1 := tr C f , can be given the form where H := 1 2 tr(∇ s ν) and K := det(∇ s ν) (12) are the mean and Gaussian curvatures of S , respectively.
Proof We see from equations (8a)-(8c) and Remark 1 that The desired conclusion follows from the identity, which we now proceed to prove. First, we represent locally the curvature tensor ∇ s ν of S as where κ 1 , κ 2 are the principal curvatures of S and n 1 , n 2 , orthogonal unit vectors of V ν , are the corresponding principal directions of curvature, so that, at each place on S , ∇ s ν is a symmetric tensor acting on V ν . Since B is also a symmetric tensor acting on V ν , it can be represented in the frame (n 1 , n 2 ) as B = B 11 n 1 ⊗ n 1 + B 22 n 2 ⊗ n 2 + B 12 (n 1 ⊗ n 2 + n 2 ⊗ n 1 ).
Equation (11) shows that I 1 involves at most quadratic terms in φ and φ . As a consequence of the orthogonal decomposition in (8a), the second and third invariants of C f , I 2 and I 3 , will also involve higher powers of φ, but not of φ . To justify the power expansion of the stored elastic energy considered in the following, we need only retain in I 2 and I 3 the terms at most quadratic in φ; all higher powers of φ will be neglected.

Proposition 3
The third invariant, I 3 := det C f , is expressed by Proof First, we note that, by (8a), Then we consider two elementary identities valid for any second-order tensor A on V 3 : for any orthonormal basis (e 1 , e 2 ) of V 3 , det A = Ae 1 × Ae 2 · e 3 and tr A = (Ae 1 × e 2 + e 1 × Ae 2 ) · e 3 .
Making use of these identities, we readily obtain from (8b) that Basic properties of trace and determinant ensure that Moreover, by the Cayley-Hamilton theorem, which together with (19), (21), and (22a)-(22c) lead us to (18).

Proposition 4 The second invariant I 2 of C f is given by
Proof It follows from (8a) that Since C φ is tensor on V 3 , again by the Cayley-Hamilton theorem, it satisfies where I 2 is the identity on V 3 . Taking the trace of both sides of (25), we obtain from (24) that The desired conclusion then follows from (21), (22a)-(22c), and (11), since tr In Sect. 4, we shall consider materials that obey the incompressibility constraint, I 3 = 1. There, (18) will turn into a differential equation that determines φ.
In preparation for this, in the following section we refresh the preliminaries of differential geometry of surfaces in a way that avoids local charts of coordinates, but resorts instead to a number of vector fields, which describe the correspondence between local movable frames in the reference and current configurations of a material surface.

Cartesian Connectors
Here, we introduce the notion of Cartesian connectors, which in our view constitute a viable alternative to Christoffel symbols. In terms of these connectors, we reformulate the classical theorema egregium of Gauss and the Codazzi-Mainardi compatibility conditions.
We reformulate the essentials of the differential geometry of smooth surfaces embedded in three-dimensional space. For definiteness, we shall assume that the mapping y that deforms S into S is of class C 3 .
Letting (r 1 , r 2 ) be the right principal directions, that is, the (normalized) eigenvectors of C, and (l 1 , l 2 ) the left principal directions, that is, the (normalized) eigenvectors of B, with corresponding principal stretches (common to both tensors) λ 1 > 0 and λ 2 > 0, we may represent ∇y, C, and B as follows (see, for example, [28, p. 74]), We shall assume that the Cartesian frames (r 1 , r 2 , e 3 ) and (l 1 , l 2 , ν) are oriented so that e 3 = r 1 × r 2 and ν = l 1 × l 2 . It should be kept in mind that both r 1 and r 2 lie in the (x 1 , x 2 ) plane; ∇ denotes the two-dimensional gradient in this plane, whereas ∇ s denotes the surface gradient on S . The connector c is a vector field in the plane such that The existence of c and the specific form of (28a), (28b) follow from the requirement that the right principal directions (r 1 , r 2 ) be orthonormal everywhere on S. 4 Clearly, if r 1 is known then c is defined as c := (∇r 1 ) T r 2 ; on the other hand, if c is assigned, at least locally, in the class C 1 , then r 1 (and r 2 ) can be determined up to a rigid rotation by solving equations (28a), (28b). To this end, however, c must be compatible; it follows from the symmetry of both ∇ 2 r 1 and ∇ 2 r 2 that the compatibility condition reads as which for a simply connected S implies that c = ∇Φ, where Φ is an appropriate scalar potential. Since we assume that both r 1 and r 2 are determined by y, we shall here consider c as known and satisfying (29).
In complete analogy to the frame (r 1 , r 2 , e 3 ) on S, we describe the corresponding frame (l 1 , l 2 , ν) as a field of orthonormal directors on S . Equations (28a), (28b) are generalized to where the connectors c * , d * 1 , and d * 2 are planar fields defined on S. A number of consequences for these fields follow from the integrability condition that requires the second gradients of y, l 1 , l 2 , and ν to be symmetric: they are listed below.
1. For the symmetry of ∇ 2 y (in its last two legs), the following second-order tensors must be symmetric, Thus, for the symmetry of ∇ 2 y, it must be In particular, (32a) and (32b) can be combined together to yield By recalling (27b), it becomes apparent from (33) that c * is completely determined by c and C.
2. Similarly, for the symmetry of ∇ 2 l 1 , it must be which can also be written in the equivalent forms 3. For the symmetry of ∇ 2 l 2 , (34a) is supplemented by or its equivalent form 4. Finally, the symmetry of ∇ 2 ν is guaranteed by (34b) and (36).
The connectors d * 1 and d * 2 can be given a geometric interpretation by computing the curvature tensor ∇ s ν of S . It readily follows from (27a) that Letting from (38) we arrive at which is duly symmetric, as by (39) equation (32c) reduces to Both the mean curvature H and the Gaussian curvature K of S can easily be derived from (40); they are given by An important conclusion follows by combining (34a) and (43) with the aid of (39), namely Since, as shown by (33), the left-hand side of (44) is determined by C (alongside its first and second spatial derivatives), so is K. In other words, the metric on S determines the Gaussian curvature of S . This is the manifestation in our setting of the celebrated theorema egregium of Gauss. Similarly, equations (34b) and (36) are related to the Codazzi-Mainardi equations (see, for example, [72, p. 144]). In the special case where both principal stretches (but not necessarily the principal directions of stretching) are uniform in space, equation (33) reduces to where c 1 := c · r 1 and c 2 := c · r 2 . It follows from (29), (45), and the identities that where we have set c 12 := r 1 · (∇c)r 2 . A comparison with (44) readily helps us to conclude that It is not difficult to check that (48) agrees completely with equation (22) of [52], which was deduced with the more traditional use of coordinates and Christoffel symbols. As appealing as formulae (42) and (43) may be, they are not especially expedient to compute H and K, for a given deformation y, as the link between the latter and the connectors is rather intricate. In the Appendix, we shall give other formulae for H and K valid for areapreserving deformations y; they are more accessible to direct computation and also show the role played by the second gradient ∇ 2 y in determining the principal curvatures of S .

Incompressible Elastomer Plates
In this section, we consider an incompressible elastomer plate, for which we assume that the mid surface S is inextensible and the whole body S is incompressible. That is, we assume that det C = 1 and det C f = 1. ( Here we consider the former constraint as a remnant of the latter, the one that survives when, in the limit as h → 0, only stretching energy is associated with the membrane S by an appropriate dimension reduction of the elastic energy stored in the three-dimensional body S. 5

Polynomial Approximation
It is our desire to compute averages of the elastic energy stored across the (small) thickness of the plate. To this end, it will suffice to represent φ as a polynomial in x 3 . Using the expressions for the invariants of C f presented in Sect. 2, in the following proposition we shall identify this polynomial.

Proposition 5 If we let
which complies with (3) for α > 0, then (49) Proof By Proposition 3, the constraints in (49) reduce to the equation where φ is as in (50). The desired result then follows by identifying in (52) the coefficients of equal powers of x 3 up to x 3 3 , and recalling that α > 0.

Remark 2
The asymptotic expansion for φ presented in Proposition 5 is consistent with a C 3 -regularity for φ in x 3 . As this hypothesis requires more regularity than that envisaged in Sect. 2, this is clearly a particular case of the framework described there.

Remark 3
The classical Kirchhoff-Love hypothesis, which requires φ ≡ x 3 , can thus be envisaged as the lowest approximation to φ in (50). The higher order approximation represented by (51) entails a local dependence of the thickness 2h * of the deformed plate on the invariant measures of curvature for S . Explicitly, this coupling is given by

Gent's Material
The elastic energy stored in a plate made of Gent's material is given by where I 1 = tr C f is the first invariant of C f , and μ and J m are positive material constants, which can be identified with a shear modulus and a stiffening parameter, respectively. The role of the latter is illuminated by the request that to which I 1 must be subjected for W G to be meaningful. Gent's constitutive law (54) was first proposed in [25]; it represents the simplest mathematical model for rubber elasticity that accounts for the limited extensibility of the polymeric chains constituting these materials. There is a vast literature on microscopic and phenomenological theories for rubber-like materials based on limited molecular extensibility, for which Beatty [5] coined the name of limited or restricted elastic models; we refer the reader to the emphatic review [31], which focuses on Gent's material.
Taking the limit as J m → ∞ in (54), we give W G the form which is the celebrated neo-Hookean stored energy density, a special case of the Mooney-Rivlin formula, where 0 < χ 1 is a dimensionless parameter. 6 Both W nH and W MR are incapable of describing the severe stiffening that occurs even at moderates stretches for soft biological membranes [30], whereas (54) is capable.
The theory proposed in [20] extended the long-standing tradition of statistical theories for ideal molecules constituted by freely joined rigid links subject to a non-Gaussian distribution for the end-to-end distance. 7 Building on this work, Beatty [4] motivated a constitutive law for rubber elasticity depending only on I 1 and incorporating the stiffening phenomena associated with a limited extensibility of the constituting chains. It was shown in [32] that (54), which has a genuine phenomenological origin, is a very accurate approximation to Beatty's molecular-based constitutive law; it retraces all qualitative features of the latter and it reproduces its quantitative predictions, with the advantage of being mathematically simpler, even amenable to explicit, closed-form solutions. Moreover, μ and J m are related to the molecular model by where n is the number density (per unit volume) of molecular chains, N is the number of links in each chain, k is the Boltzmann constant, and T the absolute temperature.
Here, we wish to show how stretching and bending energies are blended together in a thin sheet of Gent's rubber-like material complying with (49) and (55). This result will be achieved in Proposition 6 by integrating W G across the thickness of the plate using the polynomial expression for φ found in Proposition 5.

Proposition 6 Let φ be of class C 3 in x 3 , so that it can be expressed as in Proposition 5 and let the constraints (49) be enforced. Then, the following expression is valid for all
where (59c) 6 Clearly, (57) reduces to (56) for χ = 1.
Proof By inserting (11) into (54) and making use of (50) and (51), we expand W G in powers of x 3 up to x 2 3 ; standard computations lead us to where the dependence on x has been omitted to avoid clutter. Integrating this expression for W G across the thickness of the plate, we reach our desired conclusion.
Proposition 6 is the main result of this section. The quantities hw s and h 3 w 3 introduced in (59a) are interpreted as the stretching and bending elastic energy-densities (per unit area) of Gent's plates. Following [18], we shall call w s and w b the stretching and bending contents of Gent's elastic energy, respectively.

Remark 4
There are similarities between the expression for the plate's surface energy density arrived at in Proposition 6 and the elastic energy densities posited in geometric elasticity (see, for example, [3,[17][18][19]). Geometric elasticity of plates (and shells) does blend together stretching and bending energies, which scale with different powers of h; the former, like w s in (59b), is of a pure metric nature, while in the latter, as in (59c), metric and curvature measures are combined together in an invariant way. Interesting remarks on this contamination of stretching and bending measures in w b are offered in [82].
In the vanishing thickness limit, that is, as h → 0, if both w s and w b stay bounded, the stretching energy prevails over the bending energy and provides the leading deformation mechanism; for h sufficiently small, we may consider the bending energy as a perturbation to the stretching energy.

Remark 5
Among all tensors C on V 3 such that det C = 1, the stretching content w s in (59b) attains its minimum at C = I 2 , which is its unique minimizer. Indeed, letting λ 2 = 1 λ 1 in (27b), we can write w s as which attains its unique minimum for λ 1 = 1, where w s vanishes. Thus, in the absence of obstructive boundary conditions and external forces, the surfaces S that minimize w s are isometric immersions of S in three-dimensional space. By (48), all such surfaces have K = 0. Further minimizing w b on these immersions amounts at minimizing which is the form (for K = 0) of the energy density featuring in Helfrich's functional for flexible vesicles [29].
The two-step minimization outlined in Remark 5 is clearly highly hypothetical, for at least two reasons. First, boundary conditions and external forces are always present and are responsible for shaping the equilibrium configurations of plates (especially, elastomer plates, which are more responsive to mechanical stimuli). Second, the representations of stretching and bending contents in (61) and (62) miss the main points of the full-blown representations in (59b) and (59c), that is, that the measures of stretch influence the bending content as well and that both contents are blended together in (59a) in a way that depends on h and may give rise to interesting instability scenarios driven by the plate's thickness. Limiting forms of the stretching and bending contents such as w s and w b remain however indicative and will also be used in the following section to establish contact with a branch of literature in this field, where those limits have been established differently for a number of models.

Remark 6
When stretching and bending energies compete one against the other for an equilibrium, the way the surface is stretched affects its response to bending. The coupling between the two energies is not only conveyed through tr C, the sum of the principal stretches, it also involves the relative orientation of the eigenframes of the curvature tensor ∇ s ν and left Cauchy-Green tensor B.
Letting the former be represented as in (15) and the latter as in (27c), with λ 2 = 1 λ 1 , we can write where ϕ is the angle that l 1 makes with n 1 . This is the only contribution to w b that depends on ϕ. Even for given λ 1 , κ 1 , and κ 2 , minimizing w ϕ is not trivial. While w ϕ is independent of ϕ in the special case that λ 1 = 1 or κ 1 = κ 2 , in general, it has always two stationary points at ϕ = 0 and ϕ = π 2 , which are somehow expected, as there B and ∇ s ν share the same eigenframe. However, another pair of stationary points may arise, for which provided that λ 1 , κ 1 , and κ 2 make the right-hand side of (64) positive. These extra stationary points, when they exist, make w ϕ vanish, so that it attains its infimum. This shows that at equilibrium the relative orientation of B and ∇ s ν may give rise to interesting patterns on S .

Remark 7
We have assumed at the start of this section that S is inextensible and thus C is subject to det C = 1. This constraint can be easily relaxed, while still enforcing det C f = 1. A few changes occur in our analysis, which otherwise proceeds unaltered. We record here these changes for the interested reader. The polynomial representation formula for φ in Proposition 5 becomes while the expressions for w s and w b in (59b) and (59c) are to be replaced by and respectively. It is a simple matter to check that for det C = 1 equations (65), (66), and (67) reproduce the corresponding formulae derived above.
The great advantage offered by the incompressibility constraint det C f = 1 (and amply exploited in this section) is to determine φ directly on kinematic grounds, as shown in Proposition 5. For compressible materials, this advantage is lost. We need a different criterion to determine φ. In the following section, we shall show that such a criterion can be found in minimizing the elastic energy stored in the plate, for a given deformation y of the mid surface S.

Compressible Plates
In this section, we apply the method presented in Sect. 2 to compressible materials. We shall show how the modified Kirchhoff-Love hypothesis purported in this paper actually entails non-trivial normal strains for a compressible plate. Our analysis, which again is not confined to small strains, will conduce to a blending of stretching and bending energies. To ease the comparison between these latter energies and those already proposed in the literature (mostly for small strains), we shall also consider the small-strain limit for both examples we treat in detail below. In one case, we shall derive a Koiter-like potential [37][38][39][40] for the Ciarlet-Geymonat material [13]; this potential is also shown to agree with that recently derived in [14] for the same material. In the other case, we recover the bending energy derived in [22] as a rigorous Γ -limit on isometries for a variant of the Saint-Venant-Kirchhoff material.

The Ciarlet-Geymonat Material
Ciarlet and Geymonat [13] introduced a general class of hyperelastic potentials intended to provide an extension to compressible materials of the Mooney-Rivlin stored energy (see, for example, p. 189 of [10]). Here we shall consider a special example of this general class of materials, for which the stored elastic energy is where a > 0, b > 0, c > 0, and d are material constants. Letting denote the Green-Saint-Venant strain tensor (see, for example, [28, p. 70]), we show now that, for strains of sufficiently small norm |E f |, W CG can be given the classical form of the stored elastic energy for isotropic materials, with Lamé coefficients, λ and μ, appropriately related to the material constants in (68).

Proposition 7
In the limit of small strains the energy W CG in (68) can be given the form provided we set Proof It suffices to make use in (68) of the following equations Here, we continue to represent the function φ as in (50). However, no kinematic constraint will determine the functions α(x) and β(x); we need an alternative criterion, which we identify in minimizing separately the two lowest orders in h of the elastic energy integrated across the plate's thickness, for a given deformation y of the mid surface S. Hereafter, to improve clarity, the dependence on the in-plane variable x will be omitted.

Proposition 8 Let φ be given as
, with α > 0 to ensure local orientability to f in (1). For W CG as in (68), the minimum energy density (per unit area) that can be ascribed to S is represented as where and b 1 = tr (B∇ s ν). Correspondingly, α and β are determined as Proof By (11) and (18), we can write Making use of both these equations in (68), we readily arrive at which does not depend on either β or γ and, for given tr C and det C, is minimized for positive α at the value in (74c). Choosing α as in (74c), we similarly compute which is independent of γ and is minimized for β as in (74d). Inserting (74c) and (74d) in (77) and (78), respectively, we conclude the proof.

Remark 8
Although γ features in both I 1 and I 3 as expressed in (75) and (76), it does not affect w 1 in (77) and neither it does w 3 as long as α is chosen so as to minimize w 1 . Our minimization criterion leaves γ undetermined. To determine it, we should expand further the energy density w CG , so as to include terms of order h 5 , which we renounce doing here. Both w 1 and w 3 would however remain unaffected by the value of γ .
Remark 9 Equation (74c) shows clearly how in the compressible case our method differs even more markedly from the classical Kirchhoff-Love hypothesis, as α = 1 only for det C = 1. Moreover, as for incompressible materials, the bending content w 3 also depends on the relative orientation of the eigenframes of B and ∇ s ν via b 1 .

Remark 10
It is perhaps interesting to express both w 1 and w 3 in (74a) and (74b) in terms of the Lamé coefficients, λ and μ, associated with W CG in the linearized limit (70). By use of (71), we obtain that Remark 11 In the vanishing thickness limit introduced in Sect. 4, we easily find that w 1 is minimized by C = I 2 , so that correspondingly, again by Gauss' theorema egregium (which requires K = 0), w 3 takes the form It is thus useful to consider the form acquired by w CG in (73) when C is close to I 2 and, correspondingly, ∇ s ν is close to 0. Proposition 9 Let E := 1 2 (C − I 2 ). The following asymptotic representations are valid for w 1 and w 3 in (79a) and (79b), Correspondingly, α and β in (74c) and (74d) become Proof To prove (81a) and (81b) it suffices to make use of the following (simple) estimates Similarly, (81c) and (81d) follow from inserting (71) in (74c) and (74d) and then using again (82a)-(82e).

Remark 12
Both expressions for φ and w CG provided by Proposition 9 are precisely the same as those obtained in Theorem 5.2 of [14]. 8 This ensures well-posedness to the minimum energy problem in the limit of small strains. In particular, w CG with w 1 and w 2 as in (81a) and (81b) is the form appropriate to a plate of the elastic energy density envisaged in Koiter's theory for shells [39,40], in which stretching and bending energies are blended together, but are kept in a quadratic form. Equations (79a), (79b) above provide instead the stretching and bending contents for a fully nonlinear theory of plates made of the Ciarlet-Geymonat material.

A Variant of the Saint-Venant-Kirchhoff Material
Here, to provide a further application of the method proposed in this paper, we consider a variant of the classical Saint-Venant-Kirchhoff material studied in [22]. 9 The stored energy 8 It is perhaps worth recalling that (81b) is just the same as the classical formula for the strain energy stored in a moderately bent plate comprised of a linearly isotropic elastic material, see [47, p. 133], where the Lamé coefficients, λ and μ, are replaced by Young's modulus E and Poisson's ratio σ (see, for example, [47, p. 126]). Similarly, apart from a numerical prefactor due to a difference in scaling the plate's thickness, (81b) is also the same as equation (6.4) of [22], which expresses the Γ -limit on isometries of the elastic free energy of an isotropic nonlinear material, see also footnote 10 below. 9 As shown, for example, in [10, p. 155], the classical Saint-Venant-Kirchhoff material is characterized by the following stored energy function

density (per unit volume) of this material is
where λ and μ are material constants, which can be identified with the Lamé coefficients of this material, as shown by the following small-strain approximation to W VSK .
Proposition 10 Letting E f be defined as in (69), we can give W SVK the same approximate form in (68), valid for all isotropic materials, Proof The desired conclusion follows easily from remarking that A rigorous method was devised in [22] to determine the bending content w 3 of a plate on all isometric embeddings y of S in E . There, w 3 is obtained as a Γ -limit on the class of deformations that minimize the stretching energy. It was also proved in [22] that for all isotropic materials w 3 reads as the leading term in (81b) and the normal deformation φ has a quadratic representation with coefficients 10 in accord with the leading terms in (81c) and (81d). For isometric embeddings y, we can easily relax the polynomial approximation for φ. Although this refinement makes our kinematic description more accurate, the bending content w 3 is not affected, as shown below for the material with stored energy density W SVK .

Proposition 11
Let y be such that C = I 2 . Let φ in (1) be a function of class C 2 in x 3 that obeys (3). The minimum surface energy is 3 16 3 (where K = 0, since y is an isometry), which is attained for Proof W SVK can readily be rewritten as which has the same small-strain limit as (83) (see Proposition 10).

Conclusion
We have revised the classical Kirchhoff-Love hypothesis, making it more apt to derive the blending of stretching and bending energies of a plate from the free-energy functional of three-dimensional nonlinear elasticity. In summary, we have achieved two main results: (i) we have shown that measures of stretching enter the bending energy (in addition to the expected measures of bending); (ii) we have reconciled the Kirchhoff-Love hypothesis to standard Γ -convergence results for nonlinear plates on the ground where these can be compared with ours.
We have been concerned with developing a general method to obtain two-dimensional energies from three-dimensional ones and we tested it in a number of cases, thus reviving a good practice which Truesdell [75] lamented to be forgotten: "In mathematical practice today it is, unfortunately, often forgotten that to derive basic equations is even so much a mathematician's duty as to study their properties." Of course, there is much room for improvement and further extension of the proposed method.
First, the function φ introduced in (1) was almost invariably taken to be polynomial in x 3 . One wonders whether φ could be chosen in a more general class of functions without jeopardizing our conclusions. The only exploration we did along these lines was in Proposition 11, but for isometric embeddings of S; this did not affect the bending content w 3 , but had an effect on α, which changed at the order O(h 2 ), see (95). The question is then whether we can expect that, as a rule, the bending content is not affected by letting φ vary in a wider class of functions.
Second, and more importantly, the representation of the deformation f in (1) is not the most general possible. It would be interesting to replace (1) by where the unit vector d is a director field on S, which contributes to the deformation of the whole plate S on the same footing as y, representing strains across the plate's thickness.
Were we able to retrace our entire method starting from (96) instead of (1), the surface energy density w resulting from a parent volume density W would be a function of d and ∇d, as well as of y and ∇y.
Letting d · ν > 0 throughout the deformed surface S , we find ourselves in the mist of the Cosserat director-theory for plates (and shells). This theory, which goes back to the pioneering works of the Cosserat brothers [15,16], is admirably rephrased in modern terms in the book [1] (see, in particular, Chap. XIV). A full analysis of strain and equilibrium equations were first neatly developed in [21]. In connection with this theory, the classical Kirchhoff-Love hypothesis was also used in [56], always assuming d ≡ ν. A more general thermodynamic treatment of one-director surfaces was presented in [55]. As we also learn in Sect. 1.9 of [11], this theory is intimately related to the Reissner-Mindlin theory of plates [48,68,69], which indeed allows for the normals to the mid surface in the undeformed configuration not to remain normal to the deformed mid surface (as also illustrated in Sect. 5.2 of [33]).
All this body of knowledge suggests to take (96) as a general representation of the deformation field within a plate and use it to perform a dimension-reduction of the threedimensional stored energy to derive a genuine two-dimensional energy functional; a similar pursuit was undertaken in [14] (see, in particular Sect. 6.2). 11 It remains to face the difficulties offered by assuming (96) in our entire development, a task which, if not easy, might be desirable to undertake.

Conflict of interest
The authors declare that they have no conflict of interest.
Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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/.

Appendix: Cartesian formulae for H and K
Consider a deformation y of the planar surface S such that det C = 1, which ensures that y preserves the area of any portion of S. Letting (e 1 , e 2 ) be an orthonormal frame in the plane that contains S, we represent the deformation gradient ∇y as ∇y = a 1 ⊗ e 1 + a 2 ⊗ e 2 , where a 1 := (∇y)e 1 and a 2 := (∇y)e 2 .