
==== Front
Math Mech Solids
Math Mech Solids
MMS
spmms
Mathematics and Mechanics of Solids
1081-2865
1741-3028
SAGE Publications Sage UK: London, England

10.1177/10812865241233675
10.1177_10812865241233675
Articles
A note on the Hill–Ogden generalised strains
https://orcid.org/0000-0003-0866-1121
Federico Salvatore
Department of Mechanical & Manufacturing Engineering, University of Calgary, Calgary, AB, Canada
Salvatore Federico, Department of Mechanical & Manufacturing Engineering, University of Calgary, 2500 University Drive NW, Calgary, AB T2N 1N4, Canada. Email: salvatore.federico@ucalgary.ca
4 5 2024
10 2024
29 10 19351946
21 1 2024
31 1 2024
© The Author(s) 2024
2024
SAGE Publications
https://creativecommons.org/licenses/by-nc/4.0/ This article is distributed under the terms of the Creative Commons Attribution-NonCommercial 4.0 License (https://creativecommons.org/licenses/by-nc/4.0/) which permits non-commercial use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access pages (https://us.sagepub.com/en-us/nam/open-access-at-sage).
This brief contribution provides an overview of the Hill–Ogden generalised strain tensors, and some considerations on their representation in generalised (curvilinear) coordinates, in a fully covariant formalism that is adaptable to a more general theory on Riemannian manifolds. These strains may be naturally defined with covariant components or naturally defined with contravariant components. Each of these two macro-families is best suited to a specific geometrical context.

Metric tensor
strain tensor
generalised strain
polar decomposition
generalised coordinates
covariant formalism
Natural Sciences and Engineering Research Council of Canada https://doi.org/10.13039/501100000038 NSERC Discovery Grant RGPIN-2015-06027 typesetterts1
==== Body
pmc1. Introduction

In large-deformation Continuum Mechanics, the fundamental measure of deformation is the deformation gradient F , a two-point tensor, mapping material vectors in the body into spatial vectors in the physical space. The polar decomposition F=RU=VR decomposes the deformation gradient F into a pure deformation, described by the material right stretch tensor U , followed by a rotation, described by the two-point rotation tensor R , or into a rotation, described by the same two-point rotation tensor R , followed by a pure deformation, described by the spatial left stretch tensor V . The material Green–Lagrange strain E=12(U2−I) and the spatial Almansi–Euler strain e=12(i−V−2) could be called the fundamental strain tensors (here, we used their “mixed” form, whose components have the first index contravariant and the second covariant; I and i are the material and spatial identity tensors, respectively).

According to Truesdell [1], the Green–Lagrange strain E was first introduced by Cauchy [2] and de Saint Venant [3], and the Almansi–Euler strain e was first introduced by Almansi [4] and Hamel [5]. Another early measure of strain was the logarithmic strain, which Hencky [6] introduced in scalar form, and which Truesdell [1] defined as the strain tensor ln(U) (in this form, this strain is a “mixed” tensor) and named the Hencky strain. Yet another measure of strain, U−I (in this form, this strain is a “mixed” tensor), was introduced by Biot [7], and then called the Biot strain later by Truesdell [1]. However, Truesdell [1] did not explicitly relate the Hencky strain ln(U) and the Biot strain U−I to the standard Green–Lagrange strain E and Almansi–Euler strain e . Doyle and Ericksen [8] briefly mentioned that there is an infinity of material strain measures that can be expressed by (in their “mixed” form, and in our notation)

(1) E(m)=1m(Um−I),

where m is any non-zero real number. This family of strains includes the Green–Lagrange strain (for m=2 ), the Biot strain (for m=1 ), and the Hencky strain (for m→0 ). Hill [9] reported the strain measures (1) and noted that these can be written in the spectral decomposition

(2) E(m)=∑α=131m(λαm−1)Nα⊗Nα,

where the λα are the principal stretches, the eigenvalues of U , and Nα are the principal directions of stretch, the corresponding eigenvectors. Based on equations (1) and (2), Hill [9] proposed the generalised strain

(3) Ef=∑α=13f(λα)Nα⊗Nα,

where f:R+→R:λ↦f(λ) is a monotonically increasing function such that

(4) f(1)=0,f′(1)=1.

The properties (4) ensure that any Hill strain (3) reduces to the infinitesimal strain in the limit of small deformations [10]. This family of strain tensors, which are all material, was later extended by Ogden [10] to also include spatial strain tensors, such as the standard Almansi–Euler strain e . Ogden [10] observed that, since both functions

(5) λ↦1m(λm−1),λ↦1m(1−λ−m),

satisfy the properties (4), any strain tensor based on either U or V and on either of the forms (5) of the function f , and constructed with a spectral decomposition of the type (3), is a legitimate strain tensor. The four strain families envisioned by Ogden [10] are (in their “mixed” form and in our notation; nα=RNα are the spatial principal directions of stretch)

(6a) E(m)=1m(Um−I),E(m)=∑α=131m(λαm−1)Nα⊗Nα,

(6b) e(m)=1m(i−V−m),e(m)=∑α=131m(1−λα−m)nα⊗nα,

(6c) E(m)=1m(I−U−m),E(m)=∑α=131m(1−λα−m)Nα⊗Nα,

(6d) e(m)=1m(Vm−i),e(m)=∑α=131m(λαm−1)nα⊗nα.

In section 2.2.7 of his treatise, Ogden [10] recalls the standard material Green–Lagrange and spatial Almansi–Euler strain tensors, obtained from equation (6a) for m=2 and from equation (6b) for m=−2 , respectively, in equation (2.2.82) indicates, as examples, the families (6a) and (6d) and, in equation (2.2.88), he defines the material Almansi strain, obtained from equation (6c) for m=−2 .

This work aims at discussing the broader family (6) of Hill–Ogden generalised strains in a fully covariant formalism. We shall work in affine spaces for the sake of simplicity, but the conclusions remain valid in Riemannian manifolds. We shall see how all strain tensors in the families (6a) and (6b) are naturally defined with covariant components, and thus act ordinarily on tangent vectors, while all strain tensors in the families (6c) and (6d) are naturally defined with contravariant components, and thus act on tangent covectors parallel to the normal covectors to area elements. The latter result was already hinted in a work by Merodio and Ogden [11].

2. Theoretical background

The notation in this work follows essentially Truesdell and Noll [12] and Marsden and Hughes [13], as well as some past works [14–16]. In the modern picture of Continuum Mechanics [12,13,17–19], the physical space and deformable bodies are represented as differentiable manifolds. However, for our purposes, it suffices to consider physical space as an affine space (which is a trivial differentiable manifold; see [18]) and a body as an open subset of an affine space (again a trivial differentiable manifold).

For the full definition and properties of affine spaces, see the treatise by Epstein [18]. Briefly, an affine space is a set A , called point space, considered together with a vector space V , called modelling (or supporting) vector space and a mapping F:A×A↦V:(x,p)↦v≡p−x , called a framing and defining, for every two points x and p in A , the vector v≡p−x in V , attached at x . The difference-of-points notation can be exploited to obtain a point as the sum of a point and a vector, i.e., v≡p−x implies p=x+v .

2.1. Physical space and bodies

We represent physical space as the three-dimensional affine space S≡E3 , i.e., the familiar R3 considered as both the point space and the modelling vector space. The points in S are denoted x . The space of all vectors attached at x is the tangent space TxS and the disjoint union of all tangent spaces is the tangent bundle TS . A coordinate chart {xa} in S induces the basis {ea}a=13 in TS . Analogously, the space of all covectors attached at x (i.e., the linear maps α:TxS→R ) is called cotangent space Tx★S and the disjoint union of all cotangent spaces is the cotangent bundle T★S . In the same coordinate chart {xa} , the basis of T★S is the covector basis {ea}a=13 , associated with the vector basis {ea}a=13 via 〈ea|eb〉=δab . We recall that the action of a covector field π on a vector field y is the contraction

(7) 〈π|y〉≡〈y|π〉=πaya.

The points x in S and the vector and tensor fields on the tangent bundle TS are called spatial and, according to common practice in modern Continuum Mechanics, are denoted by lowercase symbols; similarly, lowercase indices are used for tensor components. The physical space S is equipped with the spatial metric tensor field g . By definition of metric tensor, g is symmetric (in components, gab=gba ) and positive definite. Given two vector fields y and z , the spatial metric tensor g defines their scalar product as

(8) 〈y,z〉g≡ygz=yagabzb.

The inverse metric tensor field g−1 has components gab and, given two covector fields π and φ , defines their scalar product as

(9) 〈π,φ〉g≡πg−1φ=πagabφb.

We identify a continuum body B with a conveniently (but arbitrarily) chosen placement of the body in the physical space S : this placement is assumed to be an open set in S and is called reference configuration. The points in B are denoted X and the tangent space TXB , cotangent space TX★B , tangent bundle TB , cotangent bundle T★B , vector basis {EA}A=13 in TB , and covector basis {EA}A=13 in T★B induced by a coordinate chart {XA} in B are defined as in the physical space S . The points X in B and the vector and tensor fields on the tangent bundle TB are called material and, according to common practice in modern Continuum Mechanics, are denoted by uppercase symbols; similarly, uppercase indices are used for tensor components. The body B is equipped with the material metric tensor field G , which, together with its inverse G−1 , is defined analogously to the spatial metric tensor field g and its inverse g−1 in S .

Remark 1. In a given manifold, for instance, the physical space S , there are four possible types of second-order tensor fields. The definitions of these four types of second-order tensors (seen as linear maps) and their component forms are

(10a) a:TS→TS,a=aabea⊗eb,″mixedcontravariant−covariant″,

(10b) b:T⋆S→TS,b=babea⊗eb,″contravariant″,

(10c) c:TS→T⋆S,c=cabea⊗eb,″covariant″,

(10d) d:T⋆S→T⋆S,d=dabea⊗eb,″mixedcovariant−contravariant″.

Strictly speaking, the terms “mixed”, “covariant,” and “covariant” refer to the tensor components and not to the corresponding tensors. However, as we did in section 1, we shall forgive ourselves for the venial sin of using the same terms for both components and tensors.

Remark 2. We adopt the point of view of modern tensor analysis on manifolds, and thus we strictly distinguish vectors and covectors and, in the same way, vector legs and covector legs of tensors. In contrast, in classical tensor analysis, {ea}a=13 is not considered as the covector basis associated with {ea}a=13 , but as the reciprocal vector basis, associated with {ea}a=13 by means of the metric tensor. In this sense, vectors and covectors become indistinguishable, and so do vector and covector legs of tensors. Marsden and Hughes [13] (p. 44, Box 2.1) give an exhaustive explanation on the difference between the two points of view.

2.2. Motion and deformation gradient

A motion of the body B is a map ϕ(·,t) that, at each time t , maps points X in the reference configuration B into their current placements x in the physical space S . For our purposes, we can omit writing the variable time, and we can write the motion ϕ(·,t) at time t in the lighter notation

(11) ϕ:B→S:X↦x=ϕ(X).

The motion is required to be an embedding, i.e., it must be such that its codomain restriction to the current configuration ϕ(B) is a diffeomorphism, so that ϕ is differentiable and invertible, and its inverse Φ≡ϕ−1 is differentiable.

The tangent map or differential of the motion ϕ at a point X is the two-point tensor field defined by

(12) F(X):TXB→TxS:M↦F(X)M≡(∂Mϕ)(X),

with x=ϕ(X) . The tensor field F is called the deformation gradient and its value F(X) at X , applied as a linear map on the material vector M∈TXB , gives the directional derivative of the motion ϕ with respect to M , at X . From the linearity of the differential, the components of the deformation gradient F are given by the partial derivatives of the motion ϕ :

(13) FaB(X)=ϕa,B(X).

The tangent map of the inverse motion Φ≡ϕ−1 is

(14) F−1(x):TxS→TXB:m↦F−1(x)m≡(∂mΦ)(x),

with X=Φ(x) , and it is called inverse deformation gradient. Its value F−1(x) at x , applied as a linear map on the spatial vector m∈TxS , gives the directional derivative of the inverse motion Φ≡ϕ−1 with respect to m , at x . Its components are

(15) (F−1)Ab(x)=ΦA,b(x).

Note that, while the deformation gradient F is a function of the material point X , the inverse deformation gradient F−1 is a function of the spatial point x .

The transpose (or algebraic transpose [14] or dual [13]) of the deformation gradient F is the two-point tensor field

(16) FT(x):Tx★S→TX★B,

with X=Φ(x) , defined, for every material vector field M and every spatial covector field π , by

(17a) 〈M(X)|FT(x)π(x)〉=〈π(x)|F(X)M(X)〉,

(17b) MA(X)(FT)Ab(x)πb(x)=πb(x)FbA(X)MA(X),

(note that this is the general definition of transpose of a tensor), so that, for the arbitrariness of M and π , its components are

(18) (FT)Ab(x)=FbA(X).

Note that the transpose deformation gradient FT is, like the inverse F−1 , a function of the spatial point x . Analogously, the transpose of the inverse deformation gradient F−1 is the two-point tensor field

(19) F−T(X):TX★S→Tx★S,

with x=ϕ(X) , defined, for every material covector field Π and every spatial vector field m , by

(20a) 〈m(x)|F−T(X)Π(X)〉=〈Π(X)|F−1(x)m(x)〉,

(20b) ma(x)(F−T)aB(X)ΠB(X)=ΠB(X)(F−1)Bama(x),

so that, by the arbitrariness of Π and m , its components are

(21) (F−T)aB(X)=(F−1)Ba(x).

Note that F−T is, like F , a function of the material point X .

We shall also need the definition of transpose induced by the metric tensor, which we call metric transpose (and is called simply transpose by Marsden and Hughes [13]). The most general definition of metric transpose is for a two-point tensor, so we shall introduce it for the case of the deformation gradient F . The metric transpose of the two-point tensor field F(X):TXB→TxS is the tensor field

(22) Ft(x):TxS→TXB,

with X=Φ(x) , such that, for every spatial vector field y and every material vector field M ,

(23a) 〈M(X),Ft(x)y(x)〉G(X)=〈y(x),F(X)M(X)〉g(x),

(23b) MB(X)GBL(X)(Ft)Lb(x)yb(x)=yb(x)gba(x)FaB(X)MB(X),

Enforcing the arbitrariness of y and M in equation (23) and then using the definition (18) of the transpose FT of F and the symmetry of the metric tensor, we have

(24) GBL(X)(Ft)Lb(x)=gba(x)FaB(X)=(FT)Ba(x)gab(x).

Multiplying by GAB(X) , we obtain

(25) GAB(X)GBL(X)(Ft)Lb(x)=GAB(X)(FT)Ba(x)gab(x).

Using GAB(X)GBL(X)=δAL on the left, we finally obtain the expression of the metric transpose Ft in terms of the transpose FT

(26) Ft(x)=G−1(X)FT(x)g(x),(Ft)Ab(x)=GAB(X)(FT)Ba(x)gab(x).

Note that the metric transpose Ft has the indices placed like the inverse F−1 (first material and contravariant, second spatial and covariant). Similarly, the inverse metric transpose is the tensor field

(27) F−t(X):TXB→TxS,

with x=ϕ(X) , defined as

(28) F−t(X)=g−1(x)F−T(X)G(X),(F−t)aB(X)=gab(x)(F−T)bA(X)GAB(X),

with the indices placed like F (first spatial and contravariant, second material and covariant).

2.3. Polar decomposition theorem

The Polar Decomposition Theorem applies to any tensor with non-zero determinant. The classical text by Chadwick [20] provides an elegant proof in Cartesian coordinates, generalising which to curvilinear coordinates is fairly straightforward.

Here, we present the statement and the consequences of the Polar Decomposition Theorem in the context of Continuum Mechanics, where the decomposition is applied to the deformation gradient tensor field F , whose determinant is strictly positive, i.e., detF>0 . In absolute and index notation (see, e.g., [13]), the theorem states that

(29) F(X)=R(X)U(X)=V(x)R(X),FaC(X)=RaB(X)UBC(X)=Vab(x)RbC(X).

The tensor field R is the two-point rotation, a proper two-point orthogonal tensor field with respect to the material and spatial metric tensors G and g , i.e.,

(30) Rt(x)=G−1(X)RT(x)g(x)≡R−1(x),(Rt)Ab(x)=GAB(X)(RT)Ba(x)gab(x)≡(R−1)Ab(x),

such that detR=+1 . Note that, for an orthogonal tensor field such as R , the inverse equals the metric transpose. The material tensor field U and the spatial tensor field V are the right and left stretch tensor fields and, with respect to the material and spatial metric tensors G and g , they are symmetric, i.e.,

(31a) Ut=G−1UTG=U,(Ut)AD=GAB(UT)BCGCD=UAD,

(31b) Vt=g−1VTg=V,(Vt)ad=gab(VT)bcgcd=Vad,

and positive definite, i.e.,

(32a) 〈Y,UY〉G=YAGABUBCYC>0,foreveryY≠0,

(32b) 〈y,Vy〉g=yagabVbcyc>0,foreveryy≠0.

In passing, we note that, in the general case of a two-point tensor F whose determinant is non-zero, but not necessarily positive, we have that

(33) detR=sgn(detF)={+1ifdetF>0,−1ifdetF<0,

i.e., R is proper orthogonal or improper orthogonal, depending on whether detF>0 or detF<0 .

Remark 3. In paper [14], which was written about 10 years ago, the mixed right and left Cauchy–Green deformation tensors with components UAB and Vab were denoted by the symbols U and V , in calligraphic typeface. This was because the intention was to reserve the symbols U and V to the covariant and contravariant stretch tensors with components GABUBC and Vabgbc . While there were some advantages in considering these covariant and contravariant stretch tensors as the “standard” ones, the little more maturity gained in the last 10 years suggests that it is wiser to exclusively define U and V as mixed tensors, as they arise from the Polar Decomposition Theorem written in the form (29). So, in order to avoid confusion, whenever it will be necessary to refer to results in paper [14], we shall use a sans serif typeface for the covariant and contravariant stretch tensors

(34a) U=GU,UAC=GABUBC,

(34b) V=Vg−1,Vac=Vabgbc.

The equivalence with the notation in paper [14] is illustrated in Table 1.

Table 1. Equivalence between the notation in this work and in paper [14].

	this work	paper [14]	
Mixed right stretch tensor	U , UAB	U , UAB	
Mixed left stretch tensor	V , Vab	V , Vab	
Covariant right stretch tensor	U , UAB	U , UAB	
Contravariant left stretch tensor	V , Vab	V , Vab	

2.4. Spectral and pseudo-spectral decompositions

As a corollary of the Polar Decomposition Theorem, we have that U and V are similar tensors, i.e., they are obtained from each other by rotating through R :

(35a) U(X)=Rt(x)V(x)R(X),UAQ(X)=GAB(X)(RT)Bc(x)gcd(x)Vdp(x)RpQ(X),

(35b) V(x)=R(X)U(X)Rt(x),Vaq(x)=RaB(X)UBC(X)GCD(X)(RT)Dp(x)gpq(x).

This implies that U and V share the same eigenvalues λα , called principal stretches, and their eigenvectors, the material principal directions of stretch Nα and the spatial principal directions of stretch nα , are again related through R :

(36a) Nα(X)=Rt(x)nα(x),NαA(X)=GAB(X)(RT)Bc(x)gcd(x)nαd(x),

(36b) nα(x)=R(X)Nα(X),nαa(x)=RaB(X)NαB(X).

The stretch tensors admit the spectral decompositions

(37a) U(X)=∑α=13λα(X)Nα(X)⊗Nα(X),

(37b) V(x)=∑α=13λα(X)nα(x)⊗nα(x),

where

(38a) Nα(X)=G(X)Nα(X),

(38b) nα(x)=g(x)nα(x),

are the covectors associated with Nα and nα through the material and spatial metric tensors G and g .

Through the spectral decompositions (37) and the main statement (29) of the Polar Decomposition Theorem, we also obtain the pseudo-spectral decompositions of F and R , as

(39a) F(X)=∑α=13λα(X)nα(x)⊗Nα(X),

(39b) R(X)=∑α=13nα(x)⊗Nα(X).

2.5. Fundamental deformation tensors

The fundamental deformation tensors are given by the pull-back of the spatial metric tensor g and its inverse g−1 and the push-forward of the material metric tensor G and its inverse G−1 [13]. These are the right Cauchy–Green deformation C , its inverse the Piola deformation B=C−1 , the Finger deformation c and its inverse the left Cauchy–Green deformation b=c−1 , and are given by

(40a) C(X)=FT(x)g(x)F(X),CAB(X)=(FT)Aa(x)gab(x)FbB(X),

(40b) B(X)=F−1(x)g−1(x)F−T(X),BAB(X)=(F−1)Aa(x)gab(x)(F−T)bB(X),

(40c) c(x)=F−T(X)G(X)F−1(x),cab(x)=(F−T)aA(X)GAB(X)(F−1)Bb(x),

(40d) b(x)=F(X)G−1(X)FT(x),bab(x)=FaA(X)GAB(X)(FT)Ba(x).

The relation between the deformation tensors (40) and the stretch tensors U and V is found by applying the polar decomposition (29) and exploiting the symmetry conditions (31) for U and V , and the orthogonality condition (30) for R , as

(41a) C(X)=G(X)U2(X),CAD(X)=GAB(X)UBC(X)UCD(X),

(41b) B(X)=U−2(X)G−1(X),BAD(X)=(U−1)AB(X)(U−1)BC(X)GCD(X),

(41c) c(x)=g(x)V−2(x),cad(x)=gab(x)(V−1)bc(x)(V−1)cd(x),

(41d) b(x)=V2(x)g−1(x),bad(x)=Vab(x)Vbc(x)gcd(x).

For instance, the expression (41a) of the right Cauchy–Green deformation tensor is found as

(42) C(X)=FT(x)g(x)F(X)=UT(X)RT(x)g(x)R(X)U(X)=G(X)U(X)G−1(X)RT(x)g(x)R(X)U(X)=G(X)U(X)Rt(x)R(X)U(X)=G(X)U(X)U(X),

where the orthogonality condition Rt=R−1 (equation (30)) implies Rt(x)R(X)=I , where I is the material identity tensor.

2.6. Fundamental strain tensors

The fundamental strain tensors are the material Green–Lagrange strain and the spatial Almansi–Euler strain. Both are defined by comparing a deformed metric with an undeformed metric. The Green–Lagrange strain E(X) compares the pull-back metric C(X)=FT(x)g(x)F(X) to the material metric G(X) , and the Almansi–Euler strain compares the spatial metric g(x) to the push-forward metric c(x)=F−T(X)G(X)F−1(x) , i.e.,

(43a) E(X)=12(C(X)−G(X)),EAB(X)=12(CAB(X)−GAB(X)),

(43b) e(x)=12(g(x)−c(x)),eab(x)=12(gab(x)−cab(x)).

Note that they are defined in such a way that they are obtainable from each other via pull-back and push-forward operations:

(44a) E(X)=FT(x)e(x)F−T(X),

(44b) e(x)=F−T(X)E(X)FT(x).

We also note that these two fundamental strains are both covariant and that, normally, no contravariant counterparts are considered.

The Green–Lagrange and Almansi–Euler strains can be written in terms of the right and left stretch tensors using equations (41), as

(45a) E(X)=12(G(X)U2(X)−G(X)),EAD(X)=12(GAB(X)UBC(X)UCD(X)−GAD(X)),

(45b) e(x)=12(g(x)−g(x)V−2(x)),ead(x)=12(gad(x)−gab(x)(V−1)bc(x)(V−1)cd(x)).

Alternatively, the metric tensors can be factorised, so that the terms in parentheses are always mixed tensors, i.e.,

(46a) E(X)=12G(X)(U2(X)−I),EAD(X)=12GAB(X)(UBC(X)UCD(X)−δBD),

(46b) e(x)=12g(x)(i−V−2(x)),ead(x)=12gab(x)(δbd−(V−1)bc(x)(V−1)cd(x)),

where I and i are the material identity tensor and spatial identity tensor, respectively.

3. Hill–Ogden generalised strains

Based on the formalism of section 2.6, we propose to write the Hill–Ogden generalised strain tensors (6) in the form

(47a) E(m)=1mG(Um−I),E(m)=∑α=131m(λαm−1)Nα⊗Nα,

(47b) e(m)=1mg(i−V−m),e(m)=∑α=131m(1−λα−m)nα⊗nα,

(47c) E(m)=1m(I−U−m)G−1,E(m)=∑α=131m(1−λα−m)Nα⊗Nα,

(47d) e(m)=1m(Vm−i)g−1,e(m)=∑α=131m(λαm−1)nα⊗nα.

The families (47a) and (47b) are covariant and are constructed on the model of the Green–Lagrange strain E and Almansi–Euler strain e as written in equation (46), i.e., by factorising the metric tensors on the left. The families (47c) and (47d) are instead contravariant and are constructed using the inverse metric tensors, which are factorised on the right. Let us study the covariant and contravariant families in detail.

3.1. Covariant strain families

Defining the strain families (47a) and (47b) as covariant is natural. Indeed, for m=2 , they reduce to the Green–Lagrange strain E and the Euler–Almansi strain E written in terms of U and V as in equation (46), with the metric tensors G and g , respectively, factorised on the left. The form (46) can then be transformed back into equation (43), which is written in terms of the right Cauchy–Green deformation C and the Finger deformation c . The two families (47a) and (47b) are often considered to be the only two (e.g., [14,21]).

It is worth to look at the logarithmic Hencky strains in detail. The material and spatial logarithmic strains

(48a) E(0)=Gln(U),E(0)=∑α=13ln(λα)Nα⊗Nα,

(48b) e(0)=gln(V),e(0)=∑α=13ln(λα)nα⊗nα,

are obtained from equations (47a) and (47b) for m→0 , by exploiting the fundamental limits (e.g., [14])

(49) limm→0am−1m=ln(a),limm→01−b−mm=−ln(b−1)=ln(b),

respectively. The second of the (49) is obtained from the first by changing sign, applying the substitution a=b−1 and then applying the identity −ln(b−1)=ln(b) .

In paper [14], in which we used the contravariant left stretch tensor V=Vg−1 (see Remark 3), we defined the spatial logarithmic strain as −ln(V−1)=−ln((Vg−1)−1) . With that definition in mind, in paper [14], we mentioned that applying the identity −ln(b−1)=ln(b) to −ln(V−1) would be plain wrong, because it would entail a tensor inversion transforming the covariant tensor −ln(V−1) into the contravariant tensor ln(V) , which clearly cannot be the same.

While this statement was correct, we need to do a mea culpa for Remark 4.3 in paper [14], where we mentioned that writing the spatial Hencky strain with the function ln(·) rather than with the function −ln((·)−1) was a mistake commonly found in the literature. In fact, this is not a mistake if the Hencky strain is indented as in the definition (48b). To show this, let us consider the limit of equation (47b) for m→0 . Since the metric tensor can be factorised outside of the limit, the tensor in the limit is mixed, and we have

(50) e(0)=limm→01mg(i−V−m)=glimm→01m(i−V−m)=g(−ln(V−1))=gln(V).

This is because, since V−1 is a mixed tensor, its inverse V belongs to the same space as V−1 . So, more than 10 years after we wrote [14], we realised that, in the literature, the expression ln(V) is simply either the mixed tensor form or the Cartesian form (in which the metric g can safely be omitted) of our e(0)=gln(V) of equation (48b). We felt that this was important to rectify.

3.2. Contravariant strain families

To see why defining the strain families (47c) and (47d) as contravariant is natural, let us look at the strains obtained for m=2 . The metric tensors factorised on the right result in the emergence of the Piola deformation tensor B expressed in equation (41b) in terms of U−1 and of the left Cauchy–Green deformation tensor b expressed in equation (41d) in terms of V . Indeed, we have

(51a) E(2)=12(I−U−2)G−1,E(2)=12(G−1−B),

(51b) e(2)=12(V2−i)g−1,e(2)=12(b−g−1).

Since the Piola deformation tensor B and the left Cauchy–Green deformation tensor b are the inverses of the right Cauchy–Green deformation tensor C and the Finger deformation tensor c , respectively, the strain tensors E(2) and e(2) of equation (51) are the most natural contravariant counterparts of the covariant strain tensors E(2)≡E and e(2)≡e , i.e., the standard Green–Lagrange and Almansi–Euler strain tensors (43).

Being contravariant tensors, the tensors E(2) and e(2) of equation (51) act naturally on covectors, and they take an interesting physical meaning when they act on the covectors normal to surface elements. Let us take the case of the material contravariant strain tensor E(2) as an example. To obtain the scalar strain in the direction of the material unit covector N , we write

(52) E(2,N)=〈N|E(2)N〉=12(〈N|G−1N〉−〈N|BN〉)=12(NAGABNB−NABABNB)=12(1−NABABNB).

In equation (52), 〈N|G−1N〉=NAGABNB=∥N∥=1 is the norm of N , which is equal to one by hypothesis, while

(53) 〈N|BN〉=NAGABNB≡B:(N⊗N)≡C−1:(N⊗N)

is related to the transversely isotropic invariant I5 , as shown by Merodio and Ogden [11]. We recall that there are five invariants for transverse isotropy: the first three are those for isotropy and the additional two depend on the structure tensor M⊗M , where the unit vector M is the direction of symmetry of transverse isotropy. The expressions of the invariants are (in our notation)

(54a) I1(C)=G−1:C,

(54b) I2(C)=12[(G−1:C)2−G−1:C2],

(54c) I3(C)=det(C),

(54d) I4(C)=C:(M⊗M),

(54e) I5(C)=C2:(M⊗M),

where G−1:C is the trace of C [22], C2≡CG−1C defines the square of C , and det(C)≡det(G−1C) [14]. Via the Hamilton–Cayley Theorem, which, applied to the right Cauchy–Green deformation tensor C , states that C is a root of its characteristic equation det(C−λ2G)=0 (see Appendix 1 in Maleki et al. [23] for a discussion on the eigenvalue problem for covariant tensors), Merodio and Ogden [11] showed that the invariant I5 can be written (in our notation) as

(55) I5(C)=I1(C)I4(C)−I2(C)+I3(C)[C−1:(N⊗N)],

where N=GM is the material covector associated with the material vector M through the material metric G and

(56) I3(C)C−1=(detC)C−1

is the adjoint (or adjugate) of C . Since I3(C)=J2=(detF)2 , and B=C−1=F−1g−1F−T , it is a simple exercise to show that the last term in equation (55) is the squared norm

(57) I3(C)[C−1:(N⊗N)]=∥JF−TN∥2

of the covector JF−TN , arising from the application of Nanson’s formula to transform the material area element with normal covector N . This observation justifies the definition of the strain families (47c) and (47d) as contravariant tensors.

4. Summary

This work focussed on the Hill–Ogden generalised strains, i.e., the material and spatial Hill [9] strains (3) for which the function f with the properties (4), has either form λ↦1m(λm−1) or λ↦1m(1−λ−m) proposed by Ogden [10], as shown in equation (5). Within a fully covariant formalism, we showed how the four families of strains (47) can be grouped into two macro-families.

One macro-family contains the strain tensors that are most naturally defined as being covariant: the material ones (47a) are constructed with the right stretch tensor U and by setting f:λ↦1m(λm−1) , and the spatial ones (47b) are constructed with the left stretch tensor V and by setting f:λ↦1m(1−λ−m) . The strain tensors in these two families could be considered as “standard” strain tensors, in so far as they are covariant. Moreover, they contain, as particular cases, the most common measures of strain. The covariant material family (47a) yields the Green–Lagrange strain E≡E(2)=12G(U2−I) , the Biot strain E(1)=G(U−I) , and the material Hencky strain E(0)=Gln(U) for m=2 , m=1 and m→0 , respectively. The covariant spatial family (47b) yields the Almansi–Euler strain e≡e(2)=12g(i−V−2) and the spatial Hencky strain e(0)=gln(V) for m=2 and m→0 , respectively.

The other macro-family contains the strain tensors that are most naturally defined as being contravariant, and they are constructed by switching the form of the function f with respect to the covariant macro-family: the material ones (47c) are constructed with the right stretch tensor U and by setting f:λ↦1m(1−λ−m) , and the spatial ones (47d) are constructed with the left stretch tensor V and by setting f:λ↦1m(λm−1) . The strain tensors in these two families are “non-standard” precisely in so far as they are contravariant. Contravariant strain tensors evaluate the strain along covectors (as opposed to vectors for the case of covariant strain tensors) and thus are suited to evaluate the strain relative to the normal covectors to area elements. For m=2 , the material strain tensor is E(2)=12(I−U−2)G−1 and features the Piola deformation tensor B=U−2G−1 , while the spatial strain tensor is e(2)=12(V2−i)g−1 and features the left Cauchy–Green deformation tensor b=V2g−1 .

With a didactical intent, this work has attempted to clarify some aspects of the mathematical structure and properties of the Hill–Ogden strain tensors. The covariant formalism facilitates a rigorous presentation of the relation between the stretch tensors U and V and the strain tensors of the four families (47). The hope is that this formalism may be adopted more often in the classroom, at the graduate level, and may perhaps be hinted to also at the undergraduate level.

Dedication

This work is dedicated to Professor Raymond W. Ogden to celebrate his fundamental contributions to the Theory of Elasticity and Continuum Mechanics, on occasion of his 80th birthday. Professor Ogden has had a profound influence on my formation as a theoretical mechanist, both indirectly, through his classical treatise on non-linear elasticity [10], and directly, during his visits to the University of Calgary in 2006 and 2009, and the many discussions that I was fortunate to be able to entertain with him through the years. Grazie, Ray.

Funding: The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported in part by the Natural Sciences and Engineering Research Council of Canada, through the NSERC Discovery Programme, grant no. RGPIN-2015-06027.

ORCID iD: Salvatore Federico https://orcid.org/0000-0003-0866-1121
==== Refs
References

[1] Truesdell C. The mechanical foundations of elasticity and fluid dynamics. J Ration Mech Anal 1952; 1 : 125–300.
[2] Cauchy AL. Mémoire sur les dilatations, les condensations par un changement déformé dans un système de points matériels. Math Oeuvres 1841; 12 : 343–377.
[3] de Saint Venant AJB . Sur les pressions qui se développent à l’intérieur des corps solides lorsque les déplacements de leurs points, sans alterer l’élasticité, ne peuvent cependant pas être considérés comme très-petits. Bull Soc Philomath 1844; 5 : 26–28.
[4] Almansi E. Sulle deformazioni finite dei solidi elastici isotropi, III. Rend Lincei 1911; 20 : 26–28.
[5] Hamel G. Elementare mechanik. Leipzig and Berlin: B.G. Teubner Verlag, 1912.
[6] Hencky H. Über die Form des elastizitätsgesetzes bei ideal elastischen Stoffen. Z Tech Phys 1928; 9 : 214–223.
[7] Biot MA. Elastizitätstheorie zweiter Ordnung mit Anwendungen. Z Angew Math Mech 1940; 20 : 89–99.
[8] Doyle TC Ericksen JL. Non-linear elasticity. Adv Appl Mech 1956; 4 : 53–115.
[9] Hill R. On constitutive inequalities for simple materials—I. J Mech Phys Solids 1968; 16 : 229–242.
[10] Ogden RW. Non-linear elastic deformations. New York: Dover Publications, 1997.
[11] Merodio J Ogden RW. Material instabilities in fiber-reinforced nonlinearly elastic solids under plane deformation. Arch Mech 2002; 54 : 525–552.
[12] Truesdell C Noll W . The non-linear field theories of mechanics (Encyclopedia of physics, ed. Flügge S ), vol. 3 . Berlin: Springer-Verlag, 1965.
[13] Marsden JE Hughes TJR . Mathematical foundations of elasticity. Englewood Cliff, NJ: Prentice-Hall, 1983.
[14] Federico S. Some remarks on metric and deformation. Math Mech Solids 2015; 20 : 522–539.
[15] Federico S Alhasadi MF Grillo A. Eshelby’s inclusion theory in light of Noether’s theorem. Math Mech Complex Syst 2019; 7 : 247–285.
[16] Federico S. The Truesdell rate in continuum mechanics. Z Angew Math Phys 2022; 73 : 109.
[17] Noll W. Materially uniform bodies with inhomogeneities. Arch Ration Mech An 1967; 27 : 1–32.
[18] Epstein M. The geometrical language of continuum mechanics. Cambridge: Cambridge University Press, 2010.
[19] Segev R. Foundations of geometric continuum mechanics: geometry and duality in continuum mechanics. Springer Nature, 2023.
[20] Chadwick P. Continuum mechanics, concise theory and problems. London: George Allen & Unwin, 1976.
[21] Bonet J Wood RD. Nonlinear continuum mechanics for finite element analysis. 2nd ed. Cambridge: Cambridge University Press, 2008.
[22] Federico S. Covariant formulation of the tensor algebra of non-linear elasticity. Int J Non Linear Mech 2012; 47 : 273–284.
[23] Maleki M Hashlamoun K Herzog W , et al . Effect of structural distortions on articular cartilage permeability under large deformations. Biomech Model Mechan 2020; 19 : 317–334.
