On the existence of elastic minimizers for initially stressed materials

TL;DR

This paper explores the mathematical conditions under which elastic minimizers exist for initially stressed materials in nonlinear elasticity, extending existing theorems to include initial stress effects.

Contribution

It introduces a framework linking elastic minimizers with constitutive restrictions for initially stressed materials, extending Ball's theorem to these cases.

Findings

Proves local existence of relaxed states as implicit functions of initial stress.

Shows initial stress tensor belongs to divergence-free symmetric tensors group.

Extends existence theorem for elastic minimizers to initially stressed materials.

Abstract

A soft solid is said to be initially stressed if it is subjected to a state of internal stress in its unloaded reference configuration. Developing a sound mathematical framework to model initially stressed solids in nonlinear elasticity is key for many applications in engineering and biology. This work investigates the links between the existence of elastic minimizers and the constitutive restrictions for initially stressed materials subjected to finite deformations. In particular, we consider a subclass of constitutive responses in which the strain energy density is taken as a scalar valued function of both the deformation gradient and the initial stress tensor. The main advantage of this approach is that the initial stress tensor belongs to the group of the divergence-free symmetric tensors satisfying the boundary condition in any given reference configuration. However, it is still unclear which physical restrictions must be imposed for the well-posedness of this elastic problem. Assuming that the constitutive response depends on the choice of the reference configuration only through the initial stress tensor, under given conditions we prove the local existence of a relaxed state given by an implicit tensor function of the initial stress distribution. This tensor function is generally not unique, and can be transformed accordingly to the symmetry group of the material at fixed initial stresses. These results allow to extend Ball's existence theorem of elastic minimizers for the proposed constitutive choice of initially stressed materials.

Figures

Tables

Table 1. Table 1. List of functional spaces

Symbol	Definition
$ℒ\left({ℝ}^n\right)$	Set of all the linear applications from ${ℝ}^n$ to ${ℝ}^n$.
${ℒ}^{+}\left({ℝ}^n\right)$	Set of all the $𝖫\in ℒ\left({ℝ}^n\right)$ such that $det𝖫>0$.
$𝒪\left({ℝ}^n\right)$	Set of all the orthogonal tensors $𝖰\in ℒ\left({ℝ}^n\right)$, namely all the $𝖰$ such that ${𝖰}^T𝖰=𝖨$.
${𝒪}^{+}\left({ℝ}^n\right)$	Set of all $𝖰\in 𝒪\left({ℝ}^n\right)$ such that $det𝖰>0$.
$𝒮\left({ℝ}^n\right)$	Set of all the symmetric linear applications from ${ℝ}^n$ to ${ℝ}^n$, namely all the $𝖫\in ℒ\left({ℝ}^n\right)$ such that ${𝖫}^T=𝖫$
${ℒ}_1^{+}\left({ℝ}^n\right)$	Special unitary group, namely the subset of ${ℒ}^{+}\left({ℝ}^n\right)$ with determinant $1$.
$𝒟$	Denotes ${ℒ}^{+}\left({ℝ}^n\right)$ for a compressible material, or ${ℒ}_1^{+}\left({ℝ}^n\right)$ for an incompressible material.
$C^0(U,V)$	Set of all the continuous function from the set $U\subseteq {ℝ}^n$ to the set $V\subseteq {ℝ}^N$.
$C^k(U,V)$	Set of all the function from the set $U\subseteq {ℝ}^n$ to the set $V\subseteq {ℝ}^N$ admitting continuous derivatives of order $k$.
$L^p(U,V)$	Set of all the function from the set $U\subseteq {ℝ}^n$ to the set $V\subseteq {ℝ}^N$ with finite $L^p$ norm.
$W^{1,p}(U,V)$	Sobolev space of all the functions from the set $U\subseteq {ℝ}^n$ to the set $V\subseteq {ℝ}^N$, where both the functions and their weak partial derivatives belong to $L^p(U,V)$.

Equations

$${\mathsf{Q}}^T\mathsf{Q}=\mathsf{I}$$

$${\mathsf{T}}_0(\mathsf{I})=0,$$

$${\mathsf{P}}_0(\mathsf{F})=\frac{\partial {\psi }_0}{\partial \mathsf{F}}{\mathsf{T}}_0(\mathsf{F})=\frac{1}{\det \mathsf{F}}{\mathsf{P}}_0(\mathsf{F}){\mathsf{F}}^T,$$

$${\mathcal{L}}_{\delta }^{+}({\mathbb{R}}^3)=\left\{\mathsf{F}\in {\mathcal{L}}^{+}({\mathbb{R}}^n)|\det \mathsf{F}=\delta \right\}.$$

$${\tilde{\psi }}_0(\mathsf{F})={\psi }_0(\mathsf{F})\forall \mathsf{F}\in {\mathcal{L}}_1^{+}({\mathbb{R}}^3),$$

$${\tilde{\psi }}_0(\mathsf{F})={\psi }_0((\det \mathsf{F}{)}^{-1/3}\mathsf{F}).$$

$${\mathsf{P}}_0(\mathsf{F},p)=\frac{\partial {\tilde{\psi }}_0}{\partial \mathsf{F}}-p{\mathsf{F}}^{-T}{\mathsf{T}}_0(\mathsf{F},p)={\mathsf{P}}_0(\mathsf{F},p){\mathsf{F}}^T,$$

$$I_1(\mathsf{C})=\mathrm{tr}(\mathsf{C}),I_2(\mathsf{C})=\frac{(\mathrm{tr}\mathsf{C}{)}^2-\mathrm{tr}({\mathsf{C}}^2)}{2},I_3(\mathsf{C})=\det (\mathsf{C}),$$

$$\{\begin{array}{rlrl}& {\psi }_0(\mathsf{F})\to +\infty & & \text{when }\det \mathsf{F}\to 0^{+}\\ & {\psi }_0(\mathsf{F})\to +\infty & & \text{when }\left|\mathsf{F}\right|+\left|{\mathsf{F}}^{-1}\right|\to +\infty \end{array}$$

$$\mathsf{F}={\mathsf{F}}_{\text{e}}\mathsf{G}.$$

$$\det {\mathsf{F}}_{\mathrm{e}}=1\Rightarrow \det \mathsf{F}=\det \mathsf{G}=\delta .$$

$${\psi }_{\mathsf{G}}(\mathsf{F})=(\det \mathsf{G}){\psi }_0(\mathsf{F}{\mathsf{G}}^{-1}).$$

$${\mathsf{P}}_{\mathsf{G}}(\mathsf{F})=(\det \mathsf{G})\frac{\partial {\psi }_0(\mathsf{F}{\mathsf{G}}^{-1})}{\partial \mathsf{F}},{\mathsf{T}}_{\mathsf{G}}(\mathsf{F})=\frac{1}{\det \mathsf{F}}{\mathsf{P}}_{\mathsf{G}}(\mathsf{F}){\mathsf{F}}^T.$$

$$\mathsf{T}(\mathsf{F})=\mathsf{T}(\mathsf{FQ}),\forall \mathsf{F}\in \mathcal{D};$$

$$\psi (\mathsf{F})=\psi (\mathsf{FQ}),\forall \mathsf{F}\in \mathcal{D}.$$

$$\begin{array}{rl}{\psi }_{\mathsf{G}}(\mathsf{F})& =(\det \mathsf{G}){\psi }_0(\mathsf{F}{\mathsf{G}}^{-1})=\\ & =(\det \mathsf{G}){\psi }_0(\mathsf{F}{\mathsf{G}}^{-1}\mathsf{Q})=\\ & =(\det \mathsf{G}){\psi }_0(\mathsf{F}{\mathsf{G}}^{-1}\mathsf{QG}{\mathsf{G}}^{-1})=\\ & ={\psi }_{\mathsf{G}}(\mathsf{F}{\mathsf{G}}^{-1}\mathsf{QG})\end{array}\forall \mathsf{Q}\in {\mathcal{G}}_0$$

$${\mathcal{G}}_{\mathsf{G}}={\mathsf{G}}^{-1}{\mathcal{G}}_0\mathsf{G}.$$

$$\mathsf{T}=\mathsf{T}(\mathsf{F};\mathsf{\Sigma }(\mathbf{X})),$$

$$\mathsf{\Sigma }:{\Omega }_0\to \mathcal{S}({\mathbb{R}}^3)$$

$$\mathsf{T}(\mathsf{I};\mathsf{S})=\mathsf{S}\forall \mathsf{S}\in \mathcal{S}({\mathbb{R}}^3).$$

$$\mathsf{T}({\mathsf{F}}_2{\mathsf{F}}_1;\mathsf{S})=\mathsf{T}\left({\mathsf{F}}_2;\mathsf{T}({\mathsf{F}}_1;\mathsf{S})\right),\forall \mathsf{S}\in \mathcal{S}({\mathbb{R}}^3).$$

$$\psi :{\mathcal{L}}^{+}({\mathbb{R}}^3)\times \mathcal{S}({\mathbb{R}}^3)\to \mathbb{R}.$$

$$\mathsf{P}(\mathsf{F};\mathsf{S})=\frac{\partial \psi }{\partial \mathsf{F}}(\mathsf{F};\mathsf{S}),\mathsf{T}(\mathsf{F};\mathsf{S})=\frac{1}{\det \mathsf{F}}\mathsf{P}(\mathsf{F};\mathsf{S}){\mathsf{F}}^T,\mathsf{F}\in {\mathcal{L}}^{+}({\mathbb{R}}^3),\mathsf{S}\in \mathcal{S}({\mathbb{R}}^3).$$

$$\psi :{\mathcal{L}}_1^{+}({\mathbb{R}}^3)\times \mathcal{S}({\mathbb{R}}^3)\to \mathbb{R},$$

$$\tilde{\psi }(\mathsf{F};\mathsf{S})=\psi (\mathsf{F};\mathsf{S})\forall \mathsf{F}\in {\mathcal{L}}_1^{+}({\mathbb{R}}^3),\forall \mathsf{S}\in \mathcal{S}({\mathbb{R}}^3);$$

$$\tilde{\psi }(\mathsf{F};\mathsf{S})=\psi ((\det \mathsf{F}{)}^{-1/3}\mathsf{F};\mathsf{S}).$$

$$\{\begin{array}{rl}& \mathsf{P}(\mathsf{F},p;\mathsf{S})=\frac{\partial \tilde{\psi }}{\partial \mathsf{F}}(\mathsf{F};\mathsf{S})-p{\mathsf{F}}^{-T},\\ & \mathsf{T}(\mathsf{F},p;\mathsf{S})=\mathsf{P}(\mathsf{F},p;\mathsf{S}){\mathsf{F}}^T,\end{array}\mathsf{F}\in {\mathcal{L}}_1^{+}({\mathbb{R}}^3),\mathsf{S}\in \mathcal{S}({\mathbb{R}}^3).$$

$$\psi ({\mathsf{F}}_2{\mathsf{F}}_1;\mathsf{S})=(\det {\mathsf{F}}_1)\psi ({\mathsf{F}}_2;\mathsf{T}({\mathsf{F}}_1;\mathsf{S})).$$

$$\frac{1}{\det {\mathsf{F}}_1\det {\mathsf{F}}_2}\frac{\partial \psi }{\partial \mathsf{F}}({\mathsf{F}}_2{\mathsf{F}}_1;\mathsf{S}){\mathsf{F}}_1^T{\mathsf{F}}_2^T=\frac{1}{\det {\mathsf{F}}_2}\frac{\partial \psi }{\partial {\mathsf{F}}_2}({\mathsf{F}}_2;\mathsf{T}({\mathsf{F}}_1;\mathsf{S})){\mathsf{F}}_2^T.$$

$$\frac{\partial \psi }{\partial {\mathsf{F}}_2}=\frac{\partial \psi }{\partial \mathsf{F}}{\mathsf{F}}_1^T$$

Full text

Developing a sound mathematical framework to model initially stressed solids in nonlinear elasticity is key for many applications in engineering and biology. This work investigates the links between the existence of elastic minimizers and the constitutive restrictions for initially stressed materials subjected to finite deformations. In particular, we consider a subclass of constitutive responses in which the strain energy density is taken as a scalar valued function of both the deformation gradient and the initial stress tensor. The main advantage of this approach is that the initial stress tensor belongs to the group of the divergence-free symmetric tensors satisfying the boundary condition in any given reference configuration. However, it is still unclear which physical restrictions must be imposed for the well-posedness of this elastic problem. Assuming that the constitutive response depends on the choice of the reference configuration only through the initial stress tensor, under given conditions we prove the local existence of a relaxed state given by an implicit tensor function of the initial stress distribution. This tensor function is generally not unique, and can be transformed accordingly to the symmetry group of the material at fixed initial stresses. These results allow to extend Ball’s existence theorem of elastic minimizers for the proposed constitutive choice of initially stressed materials.

1. Rivlin’s legacy on constitutive equations in nonlinear physics

Ronald Rivlin made seminal contribution for the development of constitutive theories in continuum physics.

As collected by Barenblatt and Joseph [53], Rivlin wrote 31 papers on isotropic finite elasticity, 8 on the anisotropic theory of elasticity and more than 50 papers on the theory of constitutive equations.

Pioneering contributions include the rigorous mathematical theory of isotropic [50] and anisotropic [21, 62] nonlinear elasticity, that shaped our modern approach to the formulation of constitutive laws in continuum mechanics [54, 52].

Moreover, Pipkin and Rivlin [42] studied the constitutive restrictions for enforcing the invariance of the physical law by rotating simultaneously both the actual and the reference system, with particular interest in nonlinear elasticity [51]. Rivlin’s work also focused on the physical implication of the material symmetries [48] and the well-posedness of the non-linear elastic problems [47, 49].

Inspired by his works, in this paper we aim at studying the role of physical restrictions for the development of a constitutive theory of initially stressed bodies, namely elastic media that are described exploiting a reference configuration that does not coincide with the relaxed one.

2. Introduction to initially stressed materials

Existence theorems for nonlinear elastic materials have been developed over the last decades on a different mathematical ground than the ones for linear elasticity, that are classically based on the Korn’s inequality [22, 28, 11]. The required objectivity of the strain energy function of a soft materials indeed conflicts with its convexity [15]. Thus, less restrictive conditions for the existence of the finite elastic solution are imposed using a direct approach in the calculus of variations, such as quasi-convexity [37], quasi-convexity at the border [5] and polyconvexity [4]. Without attempting an exhaustive mathematical characterization of such seminal results, we emphasize that proving the existence of elastic minimizers is interwoven with the search for physical restrictions on the elastic strain energy function, that classically contains a nonlinear dependence only on the deformation gradient $\mathsf{F}$.

This work investigates the link between the existence of elastic minimizers and the constitutive assumptions for initially stressed materials subjected to finite deformations. If the unloaded reference configuration is not unstressed, we say that the elastic body is subjected to an initial stress $\mathsf{\Sigma}$, that must satisfy the equilibrium equations. Initial stresses are commonly observed in soft materials. In inert matter they can be actively controlled by an external stimulus, e.g. in hydrogels [39] and in dielectric elastomers [8]. Such stimuli, e.g. an electric field in dielectric elastomers, generate a distortion of the microstructure of the material that is made physically compatible by the emergence of an internal state of stress. In living matter, initial stresses are also known as residual stresses [26, 29, 30, 31], and they result from incompatible growth processes both in healthy and pathological conditions [63, 14]. Such residual stresses not only may enhance the functionality and the efficiency of biological structures, e.g. in arteries [10], but they may also be used to trigger a programmed shape transition through a mechanical instability, forming complex patterns such as the intestinal villi [12] or the brain sulci [6].

From a constitutive viewpoint, a well established approach to account for initial stresses is based on the multiplicative decomposition of the deformation gradient into an elastic deformation tensor and an incompatible tensor [7, 33, 34]. This method was initially applied to provide a kinematic description of crystal plasticity [45], and later adapted to describe the volumetric growth in nonlinear elastic media [55]. Assuming a material isomorphism for the strain energy function, the initial stress is constitutively related to the elastic deformation tensor from the virtual incompatible state. Thus, it has been shown that the resulting elasticity tensor constitutively depends on the initial stress [29, 30, 31], that in turn may affect the symmetry group of the material [30]. Since the incompatible tensor is not necessarily the gradient of a deformation, it maps the unloaded configuration into a virtual state that may not possess a Euclidean metric [31]. Accordingly, the main drawback of this approach is that such a virtual state may not be achieved in physical practice, not even by cutting procedures, and the incompatible tensor must be assumed a priori in order to generate a self-equilibrated state of initial stresses.

A less restrictive mathematical framework accounts for initial stresses by formulating implicit constitutive equations linking the Helmoltz free energy, the initial stress and the kinematic quantities possibly mapping the evolving natural states of the materials [44]. For soft solids, this constitutive approach has shown that there exists a far richer class of non-dissipative materials than the class of bodies that is usually understood as being elastic [43].

In this work, we consider a subclass of constitutive responses in which the strain energy function is taken as a scalar valued function of both the deformation gradient and of the initial stress. As first discussed in [57], objectivity is enforced for an initially stressed material made by an isotropic material by considering a dependence on the ten invariants of the two above mentioned tensors. Under the incompressibility constraint, it has been shown that only eight invariants are independent [59]. This method has been widely used to model initially stressed materials; applications of this theory include wave propagation in soft media [58, 41], the modeling of residual stress in living tissues [66] and the stability of residually stressed materials [13, 56, 46].

The main advantage of this approach is that the initial stress tensor $\mathsf{\Sigma}$ belongs to the group of the divergence-free symmetric tensors satisfying the boundary condition in the given reference configuration, whilst it is still unclear which physical restrictions must be imposed for the well-posedness of the elastic problem. A basic constitutive restriction known as the initial stress compatibility condition (ISCC) imposes that the Cauchy stress reduces to the initial stress when the deformation tensor is equal to the matrix identity [57, 23]. By imposing ISCC and the polyconvexity of the resulting strain energy function in the absence of initial stresses few constitutive relations have been proposed. A simple functional expression has been proposed in [36], containing material parameters that also depend on the particular choice of the reference configuration, as generally prescribed by [65]. A more restrictive constitutive class has been proposed in [23], assuming that the material parameters do not change under a change of reference configuration. This assumption has lead to define a new condition. i.e. the initial stress reference independence (ISRI) [24], that is inspired by the multiplicative decomposition approach.

This work aims at clarifying some constitutive aspects of the mathematical theory of initially stressed materials, unraveling the main implications of imposing the ISRI condition on the existence of elastic minimizers.

The article is organized as follows. In Section 3, we provide some basic kinematic and constitutive notions for nonlinear elastic materials. In Section 4, we introduce the main differences of the proposed mathematical framework for initially stressed materials with respect to the theory of elastic distortions, discussing the mechanical signification of the ISCC and the ISRI conditions. In Section 5, we prove the local existence of a relaxed state for each material point. In Section 6, we prove that the residual stresses provoke an elastic distortion on the transformation of the symmetry group. We also give en existence theorem for the elastic minimizers for the proposed constitutive choice of the initially stressed material. In Section 7, we use the proposed framework to solve the physical problem of an elastic disc subjected to an anisotropic initial stress. Finally, the results are summarized and critically discussed in the last section.

3. Background and notation

We denote by $\mathcal{L}(\mathbb{R}^{n})$ the set of all the automorphisms of $\mathbb{R}^{n}$, and with $\mathcal{L}^{+}(\mathbb{R}^{n})$ the group (with respect to the operation of function composition) of all the linear applications belonging to $\mathcal{L}(\mathbb{R}^{n})$ with positive determinant.

Let $\mathcal{O}(\mathbb{R}^{n})$ be the group such that

[TABLE]

where $\mathsf{Q}\in\mathcal{L}(\mathbb{R}^{n})$ and $\mathsf{I}$ is the identity.

We indicate with $\mathcal{O}^{+}(\mathbb{R}^{n})\subset\mathcal{O}(\mathbb{R}^{n})$ the group of all the elements of $\mathcal{O}(\mathbb{R}^{n})$ with positive determinant; if $n=3$, this group coincides with the set of the rigid rotations. We also introduce the set $\mathcal{S}(\mathbb{R}^{n})$ of all the symmetric linear applications that belong to $\mathcal{L}(\mathbb{R}^{n})$.

Let the open set $\Omega_{0}\subset\mathbb{R}^{3}$ be the reference configuration of a body and $\boldsymbol{X}\in\Omega_{0}$ the material point. We denote the deformation field by $\boldsymbol{\varphi}\in C^{2}(\Omega_{0},\,\mathbb{R}^{3})$that maps the reference domain $\Omega_{0}$ to the actual configuration $\Omega$.

Accordingly, the deformation gradient reads $\mathsf{F}=\operatorname{Grad}\boldsymbol{\varphi}$. If the body is made of a homogeneous elastic material, we assume a purely elastic constitutive behavior such that the Cauchy stress tensor $\mathsf{T}_{0}$ depends on the deformation gradient $\mathsf{F}$.

We say that the body has a relaxed reference configuration if

[TABLE]

where $\mathsf{T}_{0}$ is the Cauchy stress and $\mathsf{I}$ is the identity tensor.

If the body is composed of a hyperelastic material, we denote its strain energy density in a point $\boldsymbol{X}$ with $\psi_{0}(\mathsf{F}(\boldsymbol{X})):\mathcal{L}^{+}(\mathbb{R}^{3})\rightarrow\mathbb{R}$. Whenever appropriate, we omit the explicit dependence of the physical quantities on the material position $\boldsymbol{X}$. The first Piola–Kirchhoff and the Cauchy stress tensor are given by

[TABLE]

respectively.

In order to account for an incompressibility constraint, we introduce the following group:

[TABLE]

Accordingly, the domain of the strain energy density $\psi_{0}$ is given by $\mathcal{L}^{+}_{1}(\mathbb{R}^{3})$, where the argument is the special unitary group. However, it is convenient to introduce an extension of $\psi_{0}$ to all $\mathcal{L}^{+}(\mathbb{R}^{3})$ and then to use the method of Lagrangian multiplier to enforce the incompressibility constraint. Let $\tilde{\psi}_{0}:\mathcal{L}^{+}(\mathbb{R}^{3})\rightarrow\mathbb{R}$ such that

[TABLE]

a possible extension is given by

[TABLE]

So, the first Piola–Kirchhoff and the Cauchy stress tensors are given by

[TABLE]

where $p$ is the Lagrangian multiplier. For the sake of simplicity, in the following we will omit the distinction between $\tilde{\psi}_{0}$ and $\psi_{0}$ wherever appropriate and we denote by $\mathcal{D}$ either the group $\mathcal{L}^{+}(\mathbb{R}^{3})$, if the material is unconstrained, or the group $\mathcal{L}^{+}_{1}(\mathbb{R}^{3})$, if the material is incompressible.

We denote by

[TABLE]

the principal invariants of $\mathsf{C}$, where $\mathsf{C}=\mathsf{F}^{T}\mathsf{F}$ is the right Cauchy–Green strain tensor.

We finally introduce a fundamental notion for the existence of minimizers in nonlinear elastic materials, known as the non-degeneracy axiom [4]:

Axiom 3.1** (Non-degeneracy for a hyperelastic body).**

Let $\psi_{0}$ be a strain energy density, we say that $\psi_{0}$ is non-degenerate if

[TABLE]

where $|\mathsf{F}|=\sqrt{\operatorname{tr}(\mathsf{F}^{T}\mathsf{F})}$.

The last condition of (4) indeed ensures that the hyperelastic energy goes to infinity as soon as one of the principal invariants (3) goes to $+\infty$. If the material is incompressible, only the second equation of (4) applies. For the ease of the readers, we collect all the symbols used to denote the functional spaces in Table 1.

4. Mathematical frameworks for initially stressed materials

In this section, we summarize the basic features of two mathematical frameworks used to model nonlinear elastic materials whose unloaded reference configuration is not stress-free, namely the theory of elastic distortions and the theory of initially stressed bodies.

4.1. The theory of elastic distortions

If the relation (1) does not hold, the material is subjected to a state of stress in the reference configuration. A classical constitutive approach consists in assuming a multiplicative decomposition of the deformation gradient [55], such that:

[TABLE]

where $\mathsf{G}$ is the tensor field that describes the elastic distortion from the reference configuration to the relaxed one, whilst $\mathsf{F}_{\mathrm{e}}$ represents the elastic distortion that restores the geometrical compatibility under the action of external tractions (as depicted in Fig. 1). Since the underlying metric is not Euclidean whence $\mathsf{G}$ is not a gradient of a deformation field, it may be impossible to attain a stress-free configuration in the physical world. In the last decades, the distortion tensor $\mathsf{G}$ has been advocated to model different biological processes, such as volumetric growth [16], remodelling [20] and active strains [32, 64].

In physical practice, it is assumed the initial stress in the body is generated by a distortion of the reference configuration. Consequently, the strain energy function depends on the distorted metric, given by $\mathsf{F}\mathsf{G}^{-1}$.

If the material is incompressible, this constraint is imposed on the elastic tensor, whilst the distortion tensor also describes the local change of volume, such that:

[TABLE]

Accordingly, the strain energy density of the material is given by :

[TABLE]

From standard application of the second law of the thermodynamics in the Clausius-Duhem form, the first Piola–Kirchhoff and Cauchy stress read

[TABLE]

The theory of distortions provides a transparent explanation for the transformation law of the material properties. Let $\mathcal{G}$ be the material symmetry group of a hyperelastic material, it is defined as the set of all the tensors $\mathsf{Q}\in\mathcal{L}^{+}_{1}(\mathbb{R}^{3})$ such that

[TABLE]

where the response function $\mathsf{T}$ may eventually depend on the local distortion $\mathsf{G}$.

An equivalent definition for a hyperelastic material can be given as the set of all $\mathsf{Q}\in\mathcal{L}^{+}_{1}(\mathbb{R}^{3})$ such that

[TABLE]

If we exploit the theory of elastic distortions, let $\mathcal{G}$ be the material symmetry group of the corresponding strain energy density $\psi_{0}$. It has been shown that [19]:

[TABLE]

Thus, the material symmetry group of the initially stressed material is given by

[TABLE]

Notably, $\mathcal{G}_{\mathsf{G}}$ is the conjugate group of $\mathcal{G}_{0}$ through $\mathsf{G}$.

The main drawback of the theory of elastic distortions is that $\mathsf{G}$ has to be provided by means of a constitutive assumption. Nonetheless, since the underlying metric may not be Euclidean, the values of its components cannot be directly inferred in many physical problems. An experimental attempt to search for a stress–free configuration consists in performing several (ideally infinite) cuts in the body to release the local stresses stored inside the material [2, 63, 14, 10, 3]. Although successful in simple system models [1, 38, 12], this approach is unsuitable when interested in investigating the effect of a generic state of initial stress on the material response. In the following, we describe how this difficulty can be circumvented by building a constitutive theory that explicitly depends on the underlying spatial distribution of internal stresses.

4.2. The theory of initially stressed bodies

Alternatively, it can be assumed that the material response depends on both the deformation applied on the body and on the initial stress, intended as the existing stress field $\mathsf{\Sigma}$ in the undeformed reference configuration, i.e. the Cauchy stress when the body is undeformed. This assumption has been discussed in [29, 30, 31], such that the material response in a point $\boldsymbol{X}$ of the body reads

[TABLE]

where $\mathsf{T}$ is the Cauchy stress. We remark that the initial stress tensor field

[TABLE]

generally depends on the material position vector $\boldsymbol{X}$; we omit such an explicit notation in the following for the sake of brevity wherever appropriate. We denote by $\mathsf{S}\in\mathcal{S}(\mathbb{R}^{3})$ the specific expression of the initial stress in a point $\boldsymbol{X}$, namely $\mathsf{S}=\mathsf{\Sigma}(\boldsymbol{X})$ for a given $\boldsymbol{X}\in\Omega_{0}$.

The function $\mathsf{T}:\mathcal{D}\times\mathcal{S}(\mathbb{R}^{3})\rightarrow\mathcal{S}(\mathbb{R}^{3})$ must satisfy certain restrictions. First, in the absence of initial stresses, the strain energy function must obey the standard requirements ensuring the existence of elastic minimizers in nonlinear elasticity. Second, the constitutive response should be such that the Cauchy stress is equal to the initial stress in the absence of any elastic deformations. This is referred to as ISCC, i.e. initial stress compatibility condition, [57] and reads:

[TABLE]

A subclass of material responses in which the strain energy function depends only on the elastic deformation and the initial stress, but not explicitly on the choice of the reference configuration, has been proposed in [23, 24]. Under this constitutive assumption, it is possible to introduce another restriction called Initial Stress Reference Independence (ISRI), stating that

[TABLE]

In this work we give a new mechanical interpretation of such a restrictive condition and we discuss its mathematical implications for the existence of elastic minimizers. The condition (12) imposes that there is no energy dissipation resulting from the elastic deformation and represents a frame invariance requirement: the deformation field solution of the elastic problem must not depend on the choice of the reference configuration [23, 24].

If the material is hyperelastic, we can assume that the strain energy function reads [57]:

[TABLE]

Recalling the values of the material parameters are assumed to be independent on the choice of the initially stressed configuration, the first Piola–Kirchhoff tensor $\mathsf{P}$ and the Cauchy stress tensor $\mathsf{T}$ read:

[TABLE]

Under the incompressibility constraint, the strain energy density is a function such that

[TABLE]

as done in (2), we introduce an extension $\tilde{\psi}$ of $\psi$ to all $\mathcal{L}^{+}(\mathbb{R}^{3})$ to define the stress tensors, namely

[TABLE]

a possible extension is given by:

[TABLE]

The Piola–Kirchhoff and the Cauchy stress tensors are given by

[TABLE]

For the sake of simplicity, we omit the difference between $\tilde{\psi}$ and $\psi$ wherever appropriate. For hyperelastic materials, we remind that (12) can be reformulated as an equivalent condition to be imposed on the functional dependence of the strain energy function [24]. We give further mathematical details of this important result in the following, proving that the restriction imposed on the strain energy density is a consequence of (12).

Proposition 4.1**.**

Let $\psi:\mathcal{D}\times\mathcal{S}(\mathbb{R}^{3})\rightarrow\mathbb{R}$ be a strain energy density, and assume that the ISCC (11) and the ISRI conditions (12) hold. Then, for all $\mathsf{F}_{1},\,\mathsf{F}_{2}\in\mathcal{D}$ and for all $\mathsf{S}\in\mathcal{S}(\mathbb{R}^{3})$, the following relation must hold:

[TABLE]

Proof.

For the sake of brevity, let $\psi$ be the strain energy of a compressible material (the incompressible case is analogous). The ISRI condition (12) reads

[TABLE]

Since $\mathsf{F}=\mathsf{F}_{2}\mathsf{F}_{1}$, then, by using the chain rule, we obtain

[TABLE]

and the equation (17) becomes

[TABLE]

We find that $\psi(\mathsf{F}_{2}\mathsf{F}_{1};\,\mathsf{S})=(\det{\mathsf{F}_{1}})\psi(\mathsf{F}_{2};\,\mathsf{T}(\mathsf{F}_{1};\,\mathsf{S}))+C$. Setting $\mathsf{F}_{1}=\mathsf{I}$ and making use of (11) we find that $C=0$ and we get the claim. ∎

In the next sections, we prove the local existence of a relaxed state around each material point and a theorem on the existence of elastic minimizers for a strain energy of the form given by (13).

5. Existence of a relaxed state

In the theory of elastic distortions, we must provide a constitutive form for the tensor field that we should apply locally to each point in the reference configuration to obtain the (virtual) relaxed one [55]. This theoretical framework has strong mathematical properties. Indeed, if $\psi_{0}$ is polyconvex, then also $\psi_{\mathsf{G}}$ (defined in (6)) inherits such a property [40]. As discussed earlier, this approach is straightforward but only suitable in simple system models, since it requires the a priori knowledge of the virtual relaxed state.

In this section we prove a theorem on the existence of relaxed configuration using the constitutive framework of initially stressed bodies. Moreover, we prove that, if a strain energy satisfy the ISRI (16) and $\psi(\cdot,\mathsf{0})$ is polyconvex, then $\psi(\cdot,\,\mathsf{S})$ is polyconvex for all $\mathsf{S}\in\mathcal{S}(\mathbb{R}^{3})$.

First, we give the following statement of the non-degeneracy axiom for this class of materials.

Axiom 5.1** (Non-degeneracy for an initially stressed body).**

Let $\psi:\mathcal{L}^{+}(\mathbb{R}^{3})\times\mathcal{S}(\mathbb{R}^{3})\rightarrow\mathbb{R}$ be the strain energy density of an initially stressed body. We say that $\psi$ is non-degenerate if

[TABLE]

We now derive the existence of a point-like relaxed state associated to each point of the initially stressed configuration. In [30], a stress-free virtual state is defined for each material point in the initially stressed configuration by considering the limiting behavior as the radius of the spherical neighborhood tends to zero. Its existence required the following hypotheses: $\mathsf{\Sigma}\in C^{1}(\Omega_{0},\mathcal{S}(\mathbb{R}^{3}))$, $\psi=\psi(\mathsf{F};\,\mathsf{\Sigma}(\boldsymbol{X}))$ being twice differentiable with respect to both arguments, and the distortion from the neighborhood of each point to the free state to be once differentiable in space. Here we are going to obtain a proof of existence of a virtual state using weaker hypotheses.

Theorem 5.1** (Existence of a relaxed state).**

Let $\psi$ be a non-degenerate strain energy density in the sense of (18). We also assume that $\psi(\cdot,\mathsf{S})$ is at least $C^{1}$ and proper, i.e. it is not identically equal to $+\infty$, for all $\mathsf{S}\in\mathcal{S}(\mathbb{R}^{3})$.

Then, for each $\boldsymbol{X}\in\Omega_{0}$, given $\mathsf{\Sigma}(\boldsymbol{X})\in\mathcal{S}(\mathbb{R}^{3})$, there exists a local distortion $\mathsf{G}_{\mathsf{\Sigma}(\boldsymbol{X})}$ such that

[TABLE]

Proof.

Let $\psi$ be the strain energy of a compressible material.

We denote by $f_{\mathsf{\Sigma}(\boldsymbol{X})}=\psi(\cdot;\,\mathsf{\Sigma}(\boldsymbol{X}))$. The domain of the function $f_{\mathsf{\Sigma}(\boldsymbol{X})}$ is given by $\mathcal{D}=\mathcal{L}^{+}(\mathbb{R}^{3})$. Thus, from the non-degenericity axiom (5.1), we get

[TABLE]

Since the function $f_{\mathsf{\Sigma}(\boldsymbol{X})}$ is continuous and proper, it must be bounded from below, hence there exists a value $m\in\mathbb{R}$ such that $f_{\mathsf{\Sigma}(\boldsymbol{X})}(\mathsf{F})>m$ for all $\mathsf{F}$. Moreover, there exists a value $M>m$ such that

[TABLE]

where with $f_{\mathsf{\Sigma}(\boldsymbol{X})}^{-1}(A)$ we denote the pre–image of the subset $A\subseteq\operatorname{Im}(f_{\mathsf{\Sigma}(\boldsymbol{X})})$ through the function $f_{\mathsf{\Sigma}(\boldsymbol{X})}$.

The non-empty set $U$ is bounded as a direct consequence of the coercivity property expressed in (18). Thus, there exists a minimum of $f_{\mathsf{\Sigma}(\boldsymbol{X})}$ in $\bar{U}$, where $\bar{U}$ is the closure of $U$.

The tensor that realizes such a minimum may be not unique as exposed in the following Remark 5.1. Let us denote with $\mathsf{G}_{\mathsf{\Sigma}(\boldsymbol{X})}$ one of them. Since (19) holds, the tensors $\mathsf{G}_{\mathsf{\Sigma}(\boldsymbol{X})}$ cannot belong to the boundary of the set $\mathcal{L}^{+}(\mathbb{R}^{3})$ and it is a critical point for $\psi(\cdot;\,\mathsf{\Sigma}(\boldsymbol{X}))$.

From (14), we get

[TABLE]

If the material is incompressible, following the same argument, there exists a tensor $\mathsf{G}_{\mathsf{\Sigma}(\boldsymbol{X})}$ such that

[TABLE]

where in this case $\mathcal{D}=\mathcal{L}^{+}_{1}(\mathbb{R}^{3})$. We define a new function $\hat{\psi}$ such that

[TABLE]

where $\tilde{\psi}$ is an extension of $\psi$ as defined in (15).

Since $\mathsf{G}_{\mathsf{\Sigma}(\boldsymbol{X})}$ is a minimum for $\psi$ in $\mathcal{L}^{+}_{1}(\mathbb{R}^{3})$, there exists a $p_{\mathsf{\Sigma}(\boldsymbol{X})}\in\mathbb{R}$ such that $(\mathsf{G}_{\mathsf{\Sigma}(\boldsymbol{X})},p_{\mathsf{\Sigma}(\boldsymbol{X})})$ is a critical point for $\hat{\psi}$ [18], so that

[TABLE]

and thus $\mathsf{T}(\mathsf{G}_{\mathsf{\Sigma}(\boldsymbol{X})},\,p_{\mathsf{\Sigma}(\boldsymbol{X})};\,\mathsf{\Sigma}(\boldsymbol{X}))=\mathsf{0}$. This concludes the proof. ∎

Remark 5.1**.**

Given an initial stress tensor, this Theorem implies that there exists a tensor $\mathsf{G}_{\mathsf{\Sigma}(\boldsymbol{X})}$ that locally maps the body to an unstressed state. Such a distortion is not unique in general: if $\mathsf{Q}$ belongs to the material symmetry group of $\psi$, then also $\mathsf{T}(\mathsf{G}_{\mathsf{\Sigma}(\boldsymbol{X})}\mathsf{Q};\,\mathsf{\Sigma}(\boldsymbol{X}))=\mathsf{0}$.

Remark 5.2**.**

The collection of local maps

[TABLE]

which transform each point of the reference configuration into a point in the local unstressed virtual state, satisfies $\widehat{\mathsf{G}}[\mathsf{\Sigma}]\in\mathcal{B}(\Omega_{0},\,\mathcal{D})$, where $\mathcal{B}(\Omega_{0},\,\mathcal{D})$ denotes the set of all the bounded function $f:\Omega_{0}\rightarrow\mathcal{D}$. Moreover, the tensor map $\widehat{\mathsf{G}}[\mathsf{\Sigma}]$ may not be geometrically compatible, i.e. there could not exist any differentiable vector field $\boldsymbol{\varphi}_{\widehat{\mathsf{G}}}$ such that $\operatorname{Grad}\boldsymbol{\varphi}_{\widehat{\mathsf{G}}}=\widehat{\mathsf{G}}[\mathsf{\Sigma}]$. In this case, there does not exist a deformation that maps the reference configuration of the residually stressed material into a relaxed one. In fact, assuming that the reference configuration is simply connected, such a deformation exists if and only if

[TABLE]

*In the following, we call $\widehat{\mathsf{G}}[\mathsf{\Sigma}]$ the relaxing map. *

Remark 5.3**.**

By simple application of the mean stress theorem [25], in the absence of surface tractions and body forces we obtain

[TABLE]

so that the mean value of the initial stress tensor is zero. Thus, the Cartesian components of the residual stress tensor are necessarily spatially inhomogeneous whenever $\mathsf{\Sigma}\neq\mathsf{0}$ [27]. Accordingly, the functional form of the map $\widehat{\mathsf{G}}[\mathsf{\Sigma}]$ is also inhomogeneous. A homogeneous initial stress $\mathsf{\Sigma}$ can only exist if surface tractions or body forces are applied.

Remark 5.4**.**

Note that in the case in which $\mathsf{\Sigma}$ has singular values over a set $S_{\infty}\subset\Omega_{0}$ we are requiring the hypothesis that $\psi(\cdot;\,\mathsf{\Sigma}(\boldsymbol{X}))$ remains a proper function, i.e. it is not identically equal to $+\infty$, when $\boldsymbol{X}\in S_{\infty}$. Due to the continuity of $\psi(\cdot,\mathsf{S})$, this ensures that $\widehat{\mathsf{G}}[\mathsf{\Sigma}](\boldsymbol{X})$ is bounded when $\boldsymbol{X}\in S_{\infty}$, as it will be shown in Section 7.

6. Existence of elastic minimizers for initially stressed bodies

In this section, we prove that if the strain energy density $\psi$ satisfies the assumption of Theorem 5.1, then $\psi$ satisfies the ISRI (16) if and only if it is expressible using the theory of elastic distortion (6).

Theorem 6.1**.**

Let $\psi$ satisfy the hypotheses of Theorem 5.1 and the ISCC condition. We denote by $\psi_{0}$ the strain energy of the material in the absence of initial stresses, being

[TABLE]

Then, the function $\psi$ satisfy the ISRI (16) if and only if we can express it as

[TABLE]

where $\widehat{\mathsf{G}}[\mathsf{\Sigma}](\boldsymbol{X})$ is a function such that

[TABLE]

Proof.

It is proved in [24] that the strain energy $\psi_{\widehat{\mathsf{G}}[\mathsf{\Sigma}](\boldsymbol{X})}$ satisfies the ISRI.

Let $\psi$ be a strain energy which satisfies the ISRI and such that $\psi(\mathsf{F};\,\mathsf{0})=\psi_{0}(\mathsf{F})$. The existence of the function $\widehat{\mathsf{G}}[\mathsf{\Sigma}](\boldsymbol{X})$ is guaranteed by the Theorem 5.1.

Omitting the explicit dependence on $\boldsymbol{X}$ for the sake of compactness, we obtain:

[TABLE]

that concludes the proof. ∎

Let us now introduce the following Definiton:

Definition 6.1**.**

A strain energy density $\psi_{0}(\mathsf{F})$ is said polyconvex if there exists a convex function $h:\mathbb{R}^{19}\rightarrow\mathbb{R}$ such that

[TABLE]

for all $\mathsf{F}\in\mathcal{L}^{+}(\mathbb{R}^{3}).$

We also introduce the following useful Lemma:

Lemma 6.1**.**

Let $\psi(\mathsf{F};\,\mathsf{\Sigma}(\boldsymbol{X}))$ be a strain energy density satisfying the hypotheses of Theorem 5.1, the ISRI and such that $\psi(\mathsf{F};\,\mathsf{0})$ is polyconvex. Then $\psi(\mathsf{F};\,\mathsf{\Sigma}(\boldsymbol{X}))$ is polyconvex for all $\boldsymbol{X}\in\Omega_{0}$.

Proof.

From the Theorem 6.1, we get

[TABLE]

Following the Remark 5.2, we have that $\widehat{\mathsf{G}}[\mathsf{\Sigma}]\in\mathcal{B}(\Omega_{0},\mathcal{L}^{+}(\mathbb{R}^{3}))$ and the Lemma is a direct consequence of the Lemma $6.5$ in [40]. ∎

Indeed, under some regularity assumptions, if the strain energy density $\psi$ is polyconvex in the relaxed case, then it is polyconvex for all $\mathsf{\Sigma}:\Omega_{0}\to\mathcal{S}(\mathbb{R}^{3})$. It is now possible to prove a theorem of existence of elastic minimizers for initially stressed bodies.

Theorem 6.2** (Existence of elastic minimizers).**

Let $\Omega_{0}\subset\mathbb{R}^{3}$ be a connected, bounded and open subset with a regular boundary and let $\psi(\mathsf{F};\,\mathsf{\Sigma}(\boldsymbol{X}))$ be a strain energy density for an initially stressed material, with $\psi(\cdot;\,\mathsf{\Sigma}(\boldsymbol{X}))\in C^{1}(\mathcal{L}^{+}(\mathbb{R}^{3}))$ and $\psi(\mathsf{F};\cdot)\in C^{0}(\mathcal{S}(\mathbb{R}^{3}))$. Let $\mathsf{\Sigma}:\Omega_{0}\to\mathcal{S}(\mathbb{R}^{3})$ be a measurable function.

We assume that:

• (i)

(initial stress independence and non–degeneracy) $\psi$ fulfills the hypotheses of Theorem 5.1 and the ISRI condition (16);

• (ii)

(polyconvexity of the relaxed energy) in the absence of initial stresses, the strain energy density $\psi(\mathsf{F};\,\mathsf{0})$ is polyconvex with respect to $\mathsf{F}$, namely there exists a convex function $h:\mathcal{L}(\mathbb{R}^{3})\times\mathcal{L}(\mathbb{R}^{3})\times(0,\,+\infty)\rightarrow\mathbb{R}$ such that

[TABLE]

• (iii)

(coercivity of the relaxed energy) there exist $\alpha>0$, $\beta\in\mathbb{R}$, $p\geq 2$, $q\geq p/(p-1)$, $r>1$ such that:

[TABLE]

We assume that there exist two disjointed subset $\Gamma_{0},\,\Gamma_{1}$ such that $\partial\Omega_{0}=\Gamma_{0}\cup\Gamma_{1}$ and such that $|\Gamma_{0}|>0$. Let $\boldsymbol{f}:\Omega_{0}\rightarrow\mathbb{R}^{3}$ and $\boldsymbol{t}:\Gamma_{1}\rightarrow\mathbb{R}^{3}$ measurable such that the application

[TABLE]

is continuous on $W^{1,p}(\Omega_{0},\,\mathbb{R}^{3})$. Finally let $\boldsymbol{\varphi}_{0}:\Gamma_{0}\rightarrow\mathbb{R}^{3}$ be a measurable function and such that the set

[TABLE]

is non-empty.

Then, defining the functional $\mathcal{F}:U\rightarrow\mathbb{R}\cup\{+\infty\}$ as

[TABLE]

and assuming that $\inf\mathcal{F}[\boldsymbol{\varphi}]<+\infty$, there exists an elastic minimizer

[TABLE]

Proof.

Using the Theorem 6.1, from (i) we have that

[TABLE]

We prove the claim as a direct application of the Theorem 7.3 in [4]. Here we only sketch the proof, pointing to [4] for the details.

Since $\mathsf{\Sigma}$ is measurable and $\psi(\mathsf{F},\cdot)$ is continuous for all $\mathsf{F}$, then $\Psi(\boldsymbol{X},\mathsf{F})$ is a Carathéodory function, i.e. it is continuous with respect to $\mathsf{F}$ a.e. in $\Omega_{0}$ and measurable in $\Omega_{0}$ for all $\mathsf{F}\in\mathcal{L}^{+}(\mathbb{R}^{3})$. Hence, the functional $\mathcal{F}$ is well defined.

By simple application of Lemma 6.1 and (ii), $\Psi(\boldsymbol{X},\,\mathsf{F})$ is polyconvex a.e. in $\Omega_{0}$. The coercivity of $\mathcal{F}$ is enforced a.e. in $\Omega_{0}$ by the hypothesis (iii), the boundedness of $\widehat{\mathsf{G}}[\mathsf{\Sigma}]$ in Remark 5.2 and the continuity of $L$ in $W^{1,p}(\Omega_{0},\,\mathbb{R}^{3})$. The non-degeneracy of $\Psi(\boldsymbol{X},\,\mathsf{F})$ for $\det\mathsf{F}\rightarrow 0^{+}$ is given by (i). Hence, by applying the standard methods of the calculus of variations, we can show the existence of infimizing sequences $\boldsymbol{\varphi}_{k}\in U$ which admit weakly converging subsequences to a limit point $\boldsymbol{\varphi}\in U$. Since the functional $\mathcal{F}$ is lower semicontinuous as a consequence of its policonvexity, the weak limit $\boldsymbol{\varphi}\in U$ minimizes $\mathcal{F}$.

∎

Such a Theorem is a standard application of Ball’s theorem on the existence of solutions in nonlinear elasticity [4]. The main result obtained in this section is that the ISRI automatically guarantees that the polyconvexity is preserved for all the initial stress fields if it holds for $\mathsf{\Sigma}=\mathsf{0}$. Conversely, if we do not assume the ISRI, the polyconvexity of the strain energy density should be imposed by a suitable constitutive restriction on the dependence with respect to the initial stress field.

According to the Theorem 6.1, imposing the ISRI condition is equivalent to require that the initial stress tensor $\mathsf{\Sigma}$ is generated by an elastic distortion given by $\widehat{\mathsf{G}}[\mathsf{\Sigma}]$.

Conversely, if the ISRI does not hold, the dependence of the stored elastic energy on the choice of the reference configuration is not solely related to the the corresponding variation of the initial stress. Thus, the material properties may depend on the specific initial stress field.

For the sake of clarity, let us investigate how the material symmetry group depends on the presence of an initial stress within the body. We denote by $\mathcal{G}_{0}$ the material symmetry group of the relaxed state around a material point $\boldsymbol{X}$. In view of Theorem 6.1 and following the same computation of (9), if the strain energy fulfills the ISRI condition for a generic initial stress field $\mathsf{S}$ and for all $\mathsf{Q}\in\mathcal{G}_{0}$, we get that the material symmetry group $\mathcal{G}_{\mathsf{S}}$ of the initially stressed body is given by

[TABLE]

where the tensor $\mathsf{G}_{\mathsf{S}}$ is defined in Theorem 5.1. Hence, the group $\mathcal{G}_{\mathsf{S}}$ is conjugated to the group $\mathcal{G}_{0}$ through the tensor $\mathsf{G}_{\mathsf{S}}$, exactly as in the theory of elastic distortions (9).

Conversely, we now consider a strain energy of the form

[TABLE]

where $J_{1}=\operatorname{tr}(\mathsf{S}\mathsf{C})$, and the function $f$ and $g$ must be such that the energy density (24) satisfies the ISCC (11) and $g$ is non-constant. If the material is initially unstressed (i.e. $\mathsf{S}=\mathsf{0}$), the strain energy density (24) defines a general isotropic nonlinear elastic response and the material symmetry group is given by

[TABLE]

However, if we consider an initial stress $\mathsf{S}=\alpha\boldsymbol{M}\otimes\boldsymbol{M}$, where $\boldsymbol{M}$ is a unit vector, we observe a change in the nature of the material symmetry group. In fact, considering that

[TABLE]

The material symmetry group $\mathcal{G}_{\mathsf{S}}$ is given by

[TABLE]

so that the material is not anymore isotropic but transversely isotropic.

Thus, if the material does not satisfy the ISRI, the material symmetry group $\mathcal{G}_{\mathsf{S}}$ is not conjugated with $\mathcal{G}_{0}$ and it is not possible to obtain the material symmetry group $\mathcal{G}_{0}$ by an elastic distortion of the material. In other words, if the ISRI condition is not fulfilled, the body may change its material symmetry group depending on the imposed initial stress field, leading to a modification of the material response.

7. An illustrating example: the relaxed state of a soft disc with anisotropic initial stress

As an example, we consider a disc of radius $R_{0}$ composed of an incompressible nonlinear elastic material subjected to planar strains and initial stresses. Let $(\boldsymbol{E}_{R},\,\boldsymbol{E}_{\Theta})$ and $(\boldsymbol{e}_{r},\,\boldsymbol{e}_{\theta})$ be the cylindrical vector basis in Lagrangian and Eulerian coordinates respectively. We assume that the initial stress is axis-symmetric, having the following general form

[TABLE]

The body in the reference configuration must obey the linear momentum balance, that in the absence of bulk forces reads

[TABLE]

Since the residual stress tensor $\mathsf{\Sigma}$ depends only on the radial coordinate $R$, (26) reduces to the following scalar equation

[TABLE]

that is fulfilled if and only if

[TABLE]

If the disc is not subjected to any external traction, then $\Sigma_{RR}(R_{0})=0$, so that $\alpha=0$. We now aim at calculating the elastic minimizer corresponding to this particular choice of the initial stresses. Let $\psi(\mathsf{F};\,\mathsf{\Sigma})$ be the strain energy density of the initially stressed disc. We assume that in the absence of residual stresses, the material behaves as a general isotropic material, such that

[TABLE]

where $f:[0,+\infty[\rightarrow\mathbb{R}$ is a convex function of its scalar argument. In view of Theorem 6.1, $\psi(\mathsf{F};\,\mathsf{\Sigma})=\psi(\mathsf{F}\widehat{\mathsf{G}}^{-1}[\mathsf{\Sigma}];\,\mathsf{0})$. Using the polar decomposition $\widehat{\mathsf{G}}[\mathsf{\Sigma}]=\mathsf{R}\mathsf{U}_{\widehat{\mathsf{G}}}$ where the tensor $\mathsf{R}$ is a proper orthogonal tensor and $\mathsf{U}_{\widehat{\mathsf{G}}}$ is the corresponding right stretch tensor, we denote the metric tensor of the initial elastic distortion by

[TABLE]

where $\lambda_{\widehat{\mathsf{G}}}$ is the principal eigenvalue of $\mathsf{U}_{\widehat{\mathsf{G}}}$. Accordingly, the ISCC condition (11) imposes:

[TABLE]

where $p_{\mathsf{\Sigma}}$ acts as the Lagrange multiplier enforcing the incompressibility of the metric tensor. From the expression of the initial stress (25) and the equations (27)-(28), we get that $\mathsf{U}_{\widehat{\mathsf{G}}}$ is diagonal with respect to the cylindrical vector basis, thus

[TABLE]

Considering that

[TABLE]

by enforcing the ISRI condition we can write the strain energy density as

[TABLE]

By applying the trace and the determinant operator on both sides of (28), we obtain respectively

[TABLE]

where $I_{1}(\mathsf{\Sigma})=\operatorname{tr}\mathsf{\Sigma}$ and $I_{3}(\mathsf{\Sigma})=\det\mathsf{\Sigma}$. After substituting in (29) the first equation into the second one, we get

[TABLE]

The term $\frac{I_{1}(\mathsf{\Sigma})^{2}}{4}-I_{3}(\mathsf{\Sigma})=\frac{(\Sigma_{RR}-\Sigma_{\Theta\Theta})^{2}}{4}$ is always positive. Since $f(x)$ is strictly convex with a minimum in $x=0$, the rhs of (30) is a positive–definite, strictly monotone function of $I_{1}(\tilde{\mathsf{B}})$ Thus, (30) is invertible and the principal eigenvalue $\lambda_{\widehat{\mathsf{G}}}=\lambda_{\widehat{\mathsf{G}}}(\Sigma_{RR},\Sigma_{\Theta\Theta})$ is given by:

[TABLE]

We multiply each side of (28) by $\mathsf{C}$ on the right, by applying the trace operator we get

[TABLE]

Accordingly, the strain energy function $\psi(\mathsf{F};\,\mathsf{\Sigma})$ for an initially stressed isotropic material reads

[TABLE]

Note that the relaxing map corresponding to (32) is the map $\widehat{\mathsf{G}}[\mathsf{\Sigma}]$ as defined by (27) and (28), since we have written $\psi(\mathsf{F};\,\mathsf{\Sigma})$ as $\psi(\mathsf{F}\widehat{\mathsf{G}}^{-1}[\mathsf{\Sigma}];\,\mathsf{0})$.

A mapping whose deformation gradient corresponds to $\mathsf{U}_{\widehat{\mathsf{G}}}$ is given by

[TABLE]

This relaxing map corresponds to a controllable deformation for isotropic materials, meaning that it can be supported by surface tractions alone at equilibrium. It describes the opening of the initial disc into a circular sector, corresponding to non-homogeneous displacements and homogeneous strains [61, 60]. In fact, we remark that (33) does not globally map a physically compatible configuration even if the Riemann curvature of the underlying metric tensor is zero. This can be easily checked since the curl operator of the deformation tensor corresponding to (33) is not zero if $\lambda_{\widehat{\mathsf{G}}}\neq 1$. From (31), this condition implies $\Sigma_{RR}\neq\Sigma_{\Theta\Theta}$, or equivalently $\gamma\neq 0$. Therefore, the relaxing map given by (33) is a non-uniform controllable stress state with uniform deviatoric invariants. The latter is the necessary condition for stress controllability given in [9].

8. Concluding remarks

This work proved novel insights on the link between the existence of elastic minimizers and the constitutive assumptions for initially stressed materials subjected to finite deformations.

Assuming a strain energy density in the form $\psi(\mathsf{F};\,\mathsf{\Sigma})$ and a non-degeneracy axiom, we clarified the mathematical implications of assuming the ISRI condition as a constitutive restriction. Theorem 5.1 proves the existence of a relaxed state given by the tensor function $\widehat{\mathsf{G}}[\mathsf{\Sigma}]$ as an implicit function of the initial stress distribution. The tensor $\widehat{\mathsf{G}}[\mathsf{\Sigma}]$ is generally not unique, and can be transformed accordingly to the symmetry group of $\psi$. Moreover, Theorem 6.1 proves that each strain energy density function $\psi(\mathsf{F};\,\mathsf{\Sigma})$ that satisfies the ISRI condition can be written as $\psi(\mathsf{F};\,\mathsf{\Sigma})=(\det\widehat{\mathsf{G}}[\mathsf{\Sigma}])\,\psi(\mathsf{F}\widehat{\mathsf{G}}[\mathsf{\Sigma}]^{-1};\,\mathsf{0})$. Thus, we prove that the material symmetry group of the initially stressed material satisfying the ISRI condition locally changes as we vary $\mathsf{\Sigma}$ according to the theory of elastic distortions.

Furthermore, we have used the previous results of Ball to prove the existence Theorem 6.2 of the elastic minimizers for a strain energy density in the form $\psi(\mathsf{F};\,\mathsf{\Sigma})$, that satisfies the ISRI condition under suitable constitutive restrictions. Such a result is based on the proof that the polyconvexity of the strain energy density of an initially stressed material is automatically inherited for all $\mathsf{\Sigma}$ if it holds in the case $\mathsf{\Sigma}=\mathsf{0}$, given some necessary conditions on the non-degeneracy and the regularity of $\psi(\cdot,\,\mathsf{\Sigma})$.

Whilst the theory of elastic distortion requires an a priori choice of the virtual incompatible state, the constitutive restrictions on $\psi(\mathsf{F};\,\mathsf{\Sigma})$ ensure the existence of the elastic minimizers corresponding to the physically observable distribution of the initial stresses. In an illustrative example, we have shown how to calculate the relaxed state of an incompressible isotropic disc as a function of the axis-symmetric distribution of initial stresses.

We finally remark that the ISRI condition should be assumed for the materials that do not undergo a change in the underlying material structure, so that the initial stresses arise only in response to an elastic distortion. This happens, for example, for the residual stresses generated by a differential growth. By using the classification proposed by Epstein [19], the ISRI condition is indeed well suited for modeling the growth or the remodelling of a soft material, namely a change of shape that does not affect the material properties and the microstructure of the material. On the contrary, when there is a modification of the microstructure that involves a change in the material properties, the ISRI condition would be physically flawed and other constitutive choices should be done.

Our results prove useful guidelines for the constitutive restrictions on the strain energy densities of initially stresses materials, having important applications for the study of the morphological stability and wave propagation analysis in soft tissues [14], and the non-destructive evaluation of residual stresses generated by a differential growth in biological materials [35, 17].

Funding

This work has been partially supported by Progetto Giovani GNFM 2017 funded by the National Group of Mathematical Physics (GNFM – INdAM), and by the AIRC MFAG grant 17412.

Acknowledgements

We are thankful to Davide Ambrosi, Giulia Bevilacqua and Alessandro Musesti for useful discussions about the contents of this paper. The authors are members of the National Group of Mathematical Physics of the Istituto Nazionale di Alta Matematics, GNFM–INdAM.