Residually Stressed Fiber Reinforced Solids: A Spectral Approach

Abstract

We use a spectral approach to model residually stressed elastic solids that can be applied to carbon fiber reinforced solids with a preferred direction; since the spectral formulation is more general than the classical-invariant formulation, it facilitates the search for an adequate constitutive equation for these solids. The constitutive equation is governed by spectral invariants, where each of them has a direct meaning, and are functions of the preferred direction, the residual stress tensor and the right stretch tensor. Invariants that have a transparent interpretation are useful in assisting the construction of a stringent experiment to seek a specific form of strain energy function. A separable nonlinear (finite strain) strain energy function containing single-variable functions is postulated and the associated infinitesimal strain energy function is straightforwardly obtained from its finite strain counterpart. We prove that only 11 invariants are independent. Some illustrative boundary value calculations are given. The proposed strain energy function can be simply transformed to admit the mechanical influence of compressed fibers to be partially or fully excluded.

1. Introduction

The presence of residual stresses in solids has been the essence of numerous publications [1,2,3]. There is a considerable interest in the mechanical behaviour of residually stressed materials in recent years and attempts to comprehend the mechanical behaviour of residual stresses on solid materials can be found in the literature [2,4,5,6]. A review on the presence of residual stress, for example, in fiber reinforced composite materials can be found in [7]. In this paper, we focus on the modelling of the mechanical anisotropic response of a residually stressed fiber material with a preferred direction (RSPD) based on the spectral method (method that used the eigenvalues and eigenvectors of tensors) developed recently in the literature [8,9,10,11,12,13,14,15,16]. Applications can be found, for example, in the mathematical modelling of the mechanical behaviour of carbon fibre reinforced solids and other types of RSPD, like soft biological tissue. We note that, before the recent applications of spectral formulation in anisotropic solids, most anisotropic models used traditional classical invariants [17] (or their variants), where majority of them do not have a direct interpretation, to describe their strain energy functions (see, for example, [1,18]). However, the proposed strain energy function in this paper uses spectral invariants, each has a transparent meaning that is convenient for experimental design [19]. A discussion of the advantages of spectral invariants over classical invariants is given in [20].

In this communication, our objective is to develop a novel strain energy function using spectral invariants that contains only single-variable functions. These strain energy types are experimentally attractive [19] and have been fortunate in modelling non-elastic and elastic solids [8,9,10,11,12,13,14,15,16,20,21]. It is important to note that knowing the number of independent invariants facilitates a stringent development of a strain energy via an experiment [22], and in our spectral analysis, we derive that only 11 independent spectral invariants are required in the constitutive equation. Up to our present knowledge, we believe that, modelling RSPDs using 11 independent spectral invariants is novel. We must emphasize that, to the best of our knowledge, since we are not able to find an appropriate residual stress experiment data of materials with a preferred direction, this paper focuses on the development of a rigourous theoretical spectral constitutive equation based on a systematic and rigorous use of the restrictions imposed by thermodynamics, the derivation and use of the representation formulae, a consequence of the rigorous definitions of the different material behaviors and of the concept of material symmetry, and a priori restrictions that are required by well posed mathematical models.

The basic kinematic deforming body equations and basic properties of residual stresses are given in Section 2. We discuss spectral formulations in Section 3. In Section 4, a spectral strain energy function in the absence of residual stress is proposed and, in Appendix C, based on the wok of Shariff [14], its extension to fully or partially exclude the mechanical influence of compressed fibers is given. In Section 5, a prototype strain energy function that contains single-variable functions is proposed. This prototype function is used in Section 6 to study some boundary value problems. Conclusions are given in Section 7.

2. Main Equations

2.1. Basic Concepts

Unless stated otherwise, all subscripts i, j and k take the values 1–3 and the summation convention is not used. We only discuss quasi-static deformations of incompressible solids with negligible body forces. The right Cauchy-Green tensor is $C=F^{T}F=U^2$, where ${\displaystyle F={\displaystyle \frac{\partial x}{\partial X}}}$ is the deformation tensor, U is the right stretch tensor, and X and x denote the position vectors of a solid body particle, respectively, in the reference and current configurations. Since the material is incompressible, $det\left(F\right)=1$, where det denotes the determinant of a tensor. More details about the kinematics of deforming bodies and the equation of motion can be found, for example, in Ref. [23].

2.2. Residual Stress

Details on the definition of a residual stress can be found in Merodio et al. [6]. Briefly, we postulate the existence of an equilibrium stress field with zero traction on the surface of a body in a reference configuration ${\mathcal{B}}_r$; the term residual stress $S_R$ is often used for this equilibrium stress. Hence,

$$\begin{array}{c}\hfill \mathrm{Div}{S}_{\mathrm{R}}=0\mathrm{in}{\mathcal{B}}_r\end{array}$$

(1)

with the condition on the boundary

$$\begin{array}{c}\hfill S_{\mathrm{R}}N=0\mathrm{on}\partial {\mathcal{B}}_r,\end{array}$$

(2)

where N is unit outward normal to $\partial {\mathcal{B}}_r$, the boundary of ${\mathcal{B}}_r$ and $\mathrm{Div}$ is the divergence operator with respect to X. In account of (1) and (2), the residual stress has the mean value

$$\begin{array}{c}\hfill {\int }_{{\mathcal{B}}_r}{S}_{\mathrm{R}}\mathrm{d}V=0,\end{array}$$

(3)

and it is non-homogeneous.

3. Spectral Representation

Following the work of Shariff and Merodio [20], the strain energy $\mathsf{\Omega }$ is postulated as follows:

$$\begin{array}{c}\hfill \mathsf{\Omega }=W_{\left(a\right)}(U,S_R,A),A=a\otimes a,\end{array}$$

(4)

where ⊗ denotes the dyadic product. For an incompressible material, the Cauchy stress S for an incompressible solid is given by

$$\begin{array}{c}\hfill S=2F{\displaystyle \frac{\partial \mathsf{\Omega }}{\partial C}}{F}^T-pI,\end{array}$$

(5)

where the Lagrange multiplier p is associated with the constraint $detF=1$ and I is the identity tensor. In view of the description of $S_R$ in Section 2.2, the constitutive Equation (5) at the reference configuration ($F=I$) must satisfy the relation

$$\begin{array}{c}\hfill S_R=2{\displaystyle \frac{\partial W_{\left(a\right)}}{\partial C}}(I,S_R,A)-p_{0}I,\end{array}$$

(6)

where the Lagrange multiplier $p_0$ is the value of p at the reference configuration.

We note that the right stretch tensor

$$\begin{array}{c}\hfill U=\sum _{i=1}^3{\lambda }_i{u}_i\otimes u_i,\end{array}$$

(7)

where ${\lambda }_i$ (principal stretches) and $u_i$ are eigenvalues and eigenvectors and of U, respectively. With respect to the basis $\{u_1,u_2,u_3\}$, we express $\mathsf{\Omega }$ in terms of the 12 components

$$\begin{array}{c}\hfill {\lambda }_i,s_{ij}=u_i·S_R{u}_j,a_i=a·u_i,\end{array}$$

(8)

where · is the dot product between two vectors. The unit vector a implies

$$\begin{array}{c}\hfill a_3^2=1-a_1^2-a_2^2\end{array}$$

(9)

which proves that only 11 of the 12 components are independent. Since, for all rotation tensor Q,

$$\begin{array}{c}\hfill s_{ij}=u_i·S_R{u}_j=Q{u}_i·Q{S}_R{Q}^{T}Q{u}_j,a_i=a·u_i=Qa·Q{u}_i,\end{array}$$

(10)

it is clear that, with respect to Q, $s_{ij}$ and $a_i$ are invariants. We emphasize that although $a_i$ are invariants, they are not invariants for the tensor set $ST=\{U,S_R,A\}$. The invariants ${\lambda }_i,s_{ij}$ and $a_i^2$, for example, are invariants for the tensor set $ST$.

In Appendix A, relations between classical invariants [17] are given, which strengthen our claim that at most 11 invariants are independent. We emphasize that, unlike the spectral invariants given in (8), most of 18 classical invariants in the minimal integrity basis do not have a clear physical meaning.

According to Shariff [13], the strain energy function must satisfy the P-property [13]. To facilitate the construction of the P-property, the following six independent spectral invariants

$$\begin{array}{c}\hfill {\zeta }_i=s_{ii}=u_i·S_R{u}_i,{\overline{\zeta }}_i=u_i·S_R^2{u}_i\end{array}$$

(11)

are used rather than the invariants $s_{ij}$. Hence, the strain energy $\mathsf{\Omega }$ can be expressed as

$$\begin{array}{c}\hfill \mathsf{\Omega }=W_{\left(b\right)}({\lambda }_i,a_i,{\zeta }_i,{\overline{\zeta }}_i).\end{array}$$

(12)

The Lagrangean components [24] of ${\displaystyle {\displaystyle \frac{\partial \mathsf{\Omega }}{\partial C}}}$ with respect to the basis $\{u_1,u_2,u_3\}$ are required in our analysis and they are [19]:

$$\begin{array}{c}\hfill {\left({\displaystyle \frac{\partial \mathsf{\Omega }}{\partial C}}\right)}_{ii}={\displaystyle \frac{1}{2{\lambda }_i}}{\displaystyle \frac{\partial W_{\left(b\right)}}{\partial {\lambda }_i}}\left(i\mathrm{not}\mathrm{summed}\right)\end{array}$$

(13)

with shear components

$$\begin{array}{ccc}\hfill {\left({\displaystyle \frac{\partial \mathsf{\Omega }}{\partial C}}\right)}_{ij}& =& {\displaystyle \frac{{\displaystyle \frac{\partial W_{\left(b\right)}}{\partial u_i}}·u_j-{\displaystyle \frac{\partial W_{\left(b\right)}}{\partial u_j}}·u_i}{2({\lambda }_i^2-{\lambda }_j^2)}},i\ne j.\hfill \end{array}$$

(14)

The Eulerian components [24] of Cauchy stress S with respect to the basis $\{v_1,v_2,v_3\}$, where $v_i=F{U}^{-1}{u}_i$ are

$$\begin{array}{ccc}\hfill {\tau }_{ii}& =& {\lambda }_i{\displaystyle \frac{\partial W_{\left(b\right)}}{\partial {\lambda }_i}}-p,\hfill \\ \hfill {\tau }_{ij}& =& 2{\lambda }_i{\lambda }_j{\left({\displaystyle \frac{\partial \mathsf{\Omega }}{\partial C}}\right)}_{ij},i\ne j.\hfill \end{array}$$

(15)

4. Transversely Isotropic Elastic Solid without Residual Stress

Prior to constructing a prototype strain energy function for RSPDs, we initially discuss the spectral constitutive equation for transversely isotropic elastic solids in the absence of residual stress, see for example the work of Shariff [14]. In this section, we construct a general nonlinear (finite strain) spectral strain energy function for solids with a preferred direction. We must emphasize that a finite-strain energy function should be consistent with its infinitesimal-strain counterpart (see Ref. [20] for details). The “infinitesimal strain” approach has two advantages: (a) the nonlinear strain energy function contains separable single-variable functions [14], which are easier to analyse than multivariable functions and (b) the strain energy function can be easily amended to fully or partially exclude the mechanical influence of compressed fibers (see Appendix C).

4.1. Infinitesimal Strain

In infinitesimal deformations, the most general quadratic form strain energy function for an incompressible transversely isotropic material has the expression [14]

$$\begin{array}{c}\hfill W_{\left(T\right)}={\mu }_T{K}_1+2{\mu }_1{K}_2+{\displaystyle \frac{\beta }{2}}{K}_3^2,\end{array}$$

(16)

where

$$\begin{array}{c}\hfill K_1=\sum _{i=1}^3{\nu }_i^2=\mathrm{tr}{E}^2,K_2=\sum _{i=1}^3{\overline{a}}_i^2{\nu }_i^2=a·E^{2}a,K_3=\sum _{i=1}^3{\overline{a}}_i^2{\nu }_i=a·Ea,\end{array}$$

(17)

$$\begin{array}{c}\hfill {\mu }_1={\mu }_L-{\mu }_T,\end{array}$$

(18)

${\mu }_L$ and ${\mu }_T$ are shear moduli, $\beta$ is a ground-state constant that is related to other elastic constants, which have more direct physical interpretations, ${\overline{a}}_i=e_i·a$, $e_i$ is an eigenvector of the infinitesimal strain E and ${\nu }_i$ is an eigenvalue of E.

4.2. Finite Strain

With the aid of the infinitesimal form (16), we postulate a finite-strain energy function

$$\begin{array}{c}\hfill W_{\left(T\right)}=\sum _{i=1}^3\left[{\mu }_T{q}_1\left({\lambda }_i\right)+{\mu }_1{\alpha }_i{q}_2\left({\lambda }_i\right)\right]+{\displaystyle \frac{\beta }{2}}{\left(\sum _{i=1}^3{\alpha }_i{q}_3\left({\lambda }_i\right)\right)}^2,\end{array}$$

(19)

where ${\alpha }_i=a_i^2$. The restrictions [14]

$$\begin{array}{c}\hfill q_s\left(1\right)=0,s=1,2,3,q_1^{\prime }\left(1\right)=q_2^{\prime }\left(1\right)=0,q_1^{″}\left(1\right)=q_2^{″}\left(1\right)=2,q_3^{\prime }\left(1\right)=1\end{array}$$

(20)

are required so that the finite-strain energy function is consistent with infinitesimal strain theory. In view of $detF=1$, we could also impose $q_1^{\prime }\left(1\right)=q_1^{″}\left(1\right)=1$.

For $\mid {\lambda }_i-1\mid <<1$, the value of the nonlinear strain energy function is close to the value of the corresponding infinitesimal strain energy function, and for this range of strains the strain energy function can be approximated by letting

$$\begin{array}{c}\hfill q_1\left({\lambda }_i\right)=q_2\left({\lambda }_i\right)={({\lambda }_i-1)}^2,q_3={\lambda }_i-1.\end{array}$$

(21)

It is clear that $q_1,q_2$ and $q_3$ satisfy the properties in (20) and the property

$$\begin{array}{c}\hfill q_1^{″}>0,q_2^{″}>0.\end{array}$$

(22)

In view of (21), $q_1^{\prime },q_2^{\prime }$ and $q_3$ are monotonically increasing functions with $q_1^{\prime },q_2^{\prime }$ and $q_3$ negative and positive for ${\lambda }_i<1$ and ${\lambda }_i>1$, respectively. From a continuity point of view, we shall assume that the functions $q_1,q_2$ and $q_3$ have the above properties for all ranges of ${\lambda }_i$. The connection of $q_1,q_2$ and $q_3$ to the generalized Lagrangean strain tensor is given in Appendix B. The concepts of polyconvexity, strong ellipticity condition at the current configuration, convexity and stability can be used to put restrictions on the functions $q_1$, $q_2$ and $q_3$. However, the application of these concepts to our constitutive equation is beyond the scope of this paper. The amended strain energy function proposed in this Section to model the ramification of compressed fibers is given in Appendix C.

The general single-variable functions that depend on a principal stretch ${\lambda }_i$, appearing in $W_{\left(T\right)}$, facilitate the construction of a specific form of strain energy function from experimental data. We note that, recently, separable single-variable strain energy functions have been used to model both elastic and non-elastic solids [8,9,10,11,12,13,14]. We note in passing that, it is shown in Shariff [14], the constitutive Equation (19) has successfully predicted the mechanical behaviour of soft tissues. Below, we give an example, where our theory is compared with the experimental data of [25]. We strongly emphasize that the following specific forms given below to fit the experimental data of [25] is just a preliminary exercise; better functional forms could be obtained for $q_1,q_2$ and $q_3$ but it is not the intention of this paper to do so.

For rubberlike materials, the specific forms are used in Shariff [9]:

$$\begin{array}{c}\hfill q\left(x\right)=q_1\left(x\right)=q_2\left(x\right)=2(xln\left(x\right)-x+1)+d_0(-e^{1-x}+{\displaystyle \frac{x^2-4x+5}{2}})+d_1(e^{x-1}-{\displaystyle \frac{x^2+1}{2}})\end{array}$$

(23)

and

$$\begin{array}{c}\hfill q_3\left(x\right)=ln\left(x\right).\end{array}$$

(24)

We compare our theory with the uniaxial experiment of Ciarletta et al. [25] on fiber reinforced rubber, where the experimental characterization is performed using a uniaxial testing device with optical measures of the deformation, using two different reinforcing materials on a ground rubber matrix. In The uniaxial stretch is in the $u_1$ direction. The non-zero axial First Piola–Kirchhoff stress component is

$$\begin{array}{c}\hfill P_{11}=({\mu }_T+2{\mu }_1)q^{\prime }\left({\lambda }_1\right)+\beta q_3\left({\lambda }_1\right)q_3^{\prime }\left({\lambda }_1\right)-{\mu }_T{\displaystyle \frac{{\lambda }_3{q}^{\prime }\left({\lambda }_3\right)}{{\lambda }_1}}\end{array}$$

(25)

for the case when $a=u_1$ and in this case ${\displaystyle {\lambda }_3={\displaystyle \frac{1}{\sqrt{{\lambda }_1}}}}$. In Figure 1 we curve fit the $a=u_1$ data (visually) since we know that ${\displaystyle {\lambda }_3={\displaystyle \frac{1}{\sqrt{{\lambda }_1}}}}$.

However, we cannot curve fit for the case $a=u_2$, since the values of ${\lambda }_3$ are unknown. In view of this, we have no choice but to predict the experimental data using the relation

$$\begin{array}{c}\hfill P_{11}={\mu }_T(q^{\prime }\left({\lambda }_1\right)-{\displaystyle \frac{{\lambda }_3{q}^{\prime }\left({\lambda }_3\right)}{{\lambda }_1}}),\end{array}$$

(26)

where the dependence of ${\lambda }_3$ on ${\lambda }_1$ is obtained from solving the First Piola–Kirchhoff component equations $P_{22}=P_{33}=0$ and ${\lambda }_1{\lambda }_3{\lambda }_3=1$. It is clear from Figure 1, that we are able to predict and fit using the above specific forms.

5. Strain Energy for RSPD

Due to the lack of experimental data in the literature, we are impelled to postulate a simple prototype

$$\begin{array}{c}\hfill \mathsf{\Omega }=W_{\left(T\right)}+\sum _{i=1}^3{\zeta }_{i}r\left({\lambda }_i\right)=\tilde{\mathsf{\Omega }}({\lambda }_i,{\alpha }_i,{\zeta }_i)\end{array}$$

(27)

based on the work of Shariff et al. [26,27], where

$$\begin{array}{c}\hfill r\left(1\right)=0,r^{\prime }\left(1\right)=1.\end{array}$$

(28)

An example of $r\left({\lambda }_i\right)$ are $r\left({\lambda }_i\right)={\lambda }_i-1$ and $r\left({\lambda }_i\right)=ln{\lambda }_i$. In this paper, we use $r\left({\lambda }_i\right)={\lambda }_i-1$ for the results given in Section 6.

We could, of course, propose a more complex specific form (see, for example, Shariff et al. [26]), but in this paper, for simplicity, we only consider the specific form given in (27).

We note that for the case of $q_1^{\prime }\left(1\right)=q_1^{″}\left(1\right)=1$ Equation (6) becomes

$$\begin{array}{c}\hfill S_R=({\mu }_T-p_0)I+S_R,\end{array}$$

(29)

which implies

$$\begin{array}{c}\hfill {\mu }_T=p_0,\end{array}$$

(30)

and for $q_1^{\prime }\left(1\right)=0,q_1^{″}\left(1\right)=2$, Equation (6) takes the form

$$\begin{array}{c}\hfill S_R=-p_{0}I+S_R,\end{array}$$

(31)

and it follows that $p_0=0$.

The values of the residual stress and the ground state constants are restricted via the condition of strong ellipticity at the reference configuration (see Appendix D).

The spectral components of ${\displaystyle {\displaystyle \frac{\partial \mathsf{\Omega }}{\partial C}}}$ for the strain energy function (27) takes the form

$$\begin{array}{c}\hfill {\left({\displaystyle \frac{\partial \mathsf{\Omega }}{\partial C}}\right)}_{ii}={\displaystyle \frac{1}{2{\lambda }_i}}{\displaystyle \frac{\partial \tilde{\mathsf{\Omega }}}{\partial {\lambda }_i}}\left(i\mathrm{not}\mathrm{summed}\right)\end{array}$$

(32)

with shear components

$$\begin{array}{ccc}\hfill {\left({\displaystyle \frac{\partial \mathsf{\Omega }}{\partial C}}\right)}_{ij}& =& {\displaystyle \frac{1}{({\lambda }_i^2-{\lambda }_j^2)}}\left\{\left({\displaystyle \frac{\partial \tilde{\mathsf{\Omega }}}{\partial {\zeta }_i}}-{\displaystyle \frac{\partial \tilde{\mathsf{\Omega }}}{\partial {\zeta }_j}}\right)u_i·S_R{u}_j+\left({\displaystyle \frac{\partial \tilde{\mathsf{\Omega }}}{\partial {\alpha }_i}}-{\displaystyle \frac{\partial \tilde{\mathsf{\Omega }}}{\partial {\alpha }_j}}\right)u_i·A{u}_j\right\}.\hfill \end{array}$$

(33)

6. Boundary Value Problems

In this section we give results for a simple tension deformation of a cylinder and, expansion and contraction of a hollow sphere, which could be useful from the numerical and/or experimental point of view. The boundary value problems discussed below are for any types of RSPD and, since we are not able to find an appropriate residual stress experiment data of materials with a preferred direction, only theoretical residual stress fields are discussed in this section, which may (or may not) represent “real” residual stress fields found in RSPDs.

6.1. Residual Stress: Cylinder

The study of non-residually stressed fibre reinforced solids on a cylindrical configuration is a subject of numerous publications, see for example, Ref. [28]. The cylindrical residual stress results obtained in this paper may be used to study the mechanical influence of residual stress on cylindrical problems. To facilitate our study, we need to assume a residual stress field in the reference configuration, which is described by

$$R_a\le R\le R_b,0\le \mathsf{\Theta }\mathsf{\le }2\pi ,0\le Z\le L,$$

(34)

where R, $\mathsf{\Theta }$ and Z are reference cylindrical polar coordinates.

We consider a residual stress [6] of the form

$$S_R=s_1\left(R\right)E_R\otimes E_R+s_2\left(R\right)E_{\mathsf{\Theta }}\otimes E_{\mathsf{\Theta }},,$$

(35)

where $E_R$ and $E_{\mathsf{\Theta }}$ are cylindrical polar coordinate vectors in the reference configuration.

Satisfaction of the equilibrium equation and boundary conditions, require

$$R{\displaystyle \frac{d{s}_1}{dR}}=s_2-s_1,$$

(36)

and

$$s_1\left(R_a\right)=0,s_1\left(R_b\right)=0.$$

(37)

A simple example of $s_1$ is

$$s_1=\overline{\alpha }(R-R_a)(R-R_b),$$

(38)

where $\overline{\alpha }$ is a constant. We use (38) in the following section.

6.2. Uniform Extension of a Cylinder

In this Section, all tensor and vector components are cylindrical polar components.

We consider a uniform extension of a cylinder with $R_a=0$ described by

$$\begin{array}{c}\hfill r={\displaystyle \frac{1}{\sqrt{{\lambda }_z}}}R,\mathsf{\theta }=\Theta ,z={\lambda }_{z}Z,\end{array}$$

(39)

where $(r,\theta ,z)$ is the polar coordinate in the deformed configuration. It follows that

$$\begin{array}{c}\hfill F\equiv \mathrm{diag}({\displaystyle \frac{1}{\sqrt{{\lambda }_z}}},{\displaystyle \frac{1}{\sqrt{{\lambda }_z}}},{\lambda }_z).\end{array}$$

(40)

Therefore, ${\displaystyle {\lambda }_1={\lambda }_r={\displaystyle \frac{1}{\sqrt{{\lambda }_z}}},{\lambda }_2={\lambda }_r={\displaystyle \frac{1}{\sqrt{{\lambda }_z}}},{\lambda }_3={\lambda }_z}$, $u_1=E_R,u_2=E_{\mathsf{\Theta }}$ and $u_3=E_Z$. Here, for simplicity, we let $a=E_R$ and hence ${\alpha }_1=1,{\alpha }_2={\alpha }_3=0$.

The non-zero Cauchy stress components are ${\sigma }_{rr}$ (radial stress), ${\sigma }_{\theta \theta }$ (hoop stress) and ${\sigma }_{zz}$ (axial stress), where

$$\begin{array}{c}\hfill {\sigma }_{rr}={\lambda }_r{\displaystyle \frac{\partial \tilde{\mathsf{\Omega }}}{\partial {\lambda }_1}}-p,{\sigma }_{\theta \theta }={\lambda }_r{\displaystyle \frac{\partial \tilde{\mathsf{\Omega }}}{\partial {\lambda }_2}}-p,{\sigma }_{zz}={\lambda }_z{\displaystyle \frac{\partial \tilde{\mathsf{\Omega }}}{\partial {\lambda }_3}}-p.\end{array}$$

(41)

$\tilde{\mathsf{\Omega }}$ depends on r (or in view of (39), equivalently on R) and the Cauchy stress is inhomogeneous. It is clear that we have zero shear stresses and the equilibrium equation reduces to

$$\begin{array}{c}\hfill r\frac{\mathrm{d}{\sigma }_{rr}}{\mathrm{d}r}={\sigma }_{\theta \theta }-{\sigma }_{rr},\end{array}$$

(42)

which can be integrated to give

$$\begin{array}{c}\hfill {\sigma }_{rr}={\int }_b^r\left({\lambda }_r{\displaystyle \frac{\partial \tilde{\mathsf{\Omega }}}{\partial {\lambda }_2}}-{\lambda }_r{\displaystyle \frac{\partial \tilde{\mathsf{\Omega }}}{\partial {\lambda }_1}}\right){\displaystyle \frac{\mathrm{d}r}{r}},\end{array}$$

(43)

where ${\displaystyle b={\displaystyle \frac{R_b}{\sqrt{{\lambda }_z}}}}$. The axial stress is given by

$$\begin{array}{c}\hfill {\sigma }_{zz}={\lambda }_z{\displaystyle \frac{\partial \tilde{\mathsf{\Omega }}}{\partial {\lambda }_3}}-{\lambda }_r{\displaystyle \frac{\partial \tilde{\mathsf{\Omega }}}{\partial {\lambda }_1}}+{\sigma }_{rr}.\end{array}$$

(44)

For illustration we use the (23), (24) and (A17) for $\tilde{\mathsf{\Omega }}$ and we simply have

$$\begin{array}{c}\hfill {\sigma }_{zz}={\mu }_T[{\lambda }_z{q}^{\prime }\left({\lambda }_z\right)-{\displaystyle \frac{1}{\sqrt{{\lambda }_z}}}{q}^{\prime }\left({\displaystyle \frac{1}{\sqrt{{\lambda }_z}}}\right)],\end{array}$$

(45)

$$\begin{array}{c}\hfill {\sigma }_{rr}={\displaystyle \frac{1}{\sqrt{{\lambda }_z}}}{s}_1\left(R\right)\end{array}$$

(46)

and

$$\begin{array}{c}\hfill {\sigma }_{\theta \theta }={\displaystyle \frac{1}{\sqrt{{\lambda }_z}}}{s}_2\left(R\right),\end{array}$$

(47)

where, since $I_4\le 1$, we have considered ${\mu }_1=\beta =0$. The above indicates that ${\sigma }_{zz}$ is constant and does not depend on the residual stress. On the other hand, ${\sigma }_{rr}$ and ${\sigma }_{\theta \theta }$ are functions of the residual stress only; their absolute values decrease as ${\lambda }_z$ increases. The axial Cauchy stress ${\sigma }_{zz}$ vs ${\lambda }_z$ is shown in Figure 2 and, in Figure 3, we plot the curves of ${\sigma }_{rr}$ and ${\sigma }_{\theta \theta }$ for ${\displaystyle \overline{\alpha }=10{\displaystyle \frac{kPa}{m^2}}}$.

6.3. Spherically Symmetric Deformation of a Spherical Shell

The study of spherical shell in this section could be useful, for example, to enhance the study of cavity formation in a sphere under uniform radial tensile dead-load with the fiber in the radial direction [29]. Here, we consider a spherical shell with thick-walled having the reference geometry

$$\begin{array}{c}\hfill R_a\le R\le R_b,0\le \mathsf{\Theta }\le \pi ,0\le \mathsf{\Phi }\le 2\pi ,\end{array}$$

(48)

where $(R,\mathsf{\Theta },\mathsf{\Phi })$ is the spherical polar coordinate for the undeformed configuration. The geometry of the current deformation is described by

$$\begin{array}{c}\hfill a\le r\le b,\theta =\mathsf{\Theta },\varphi =\mathsf{\Phi },\end{array}$$

(49)

where $(r,\theta ,\varphi )$ is the spherical polar coordinate for the current configuration.

We consider a deformation defined by

$$\begin{array}{c}\hfill r=Rf\left(R\right).\end{array}$$

(50)

In the spherical polar coordinate system, the deformation gradient is

$$\begin{array}{c}\hfill F\equiv \left(\begin{array}{ccc}1+R{f}^{\prime }\left(R\right)& 0& 0\\ 0& f\left(R\right)& \\ 0& 0& f\left(R\right)\end{array}\right).\end{array}$$

(51)

Due to the incompressibility condition we have

$$\begin{array}{c}\hfill f\left(R\right)={\left(1+{\displaystyle \frac{a^3-R_a^3}{R^3}}\right)}^{\frac{1}{3}}={\displaystyle \frac{r}{R}}.\end{array}$$

(52)

The principal stretches are

$$\begin{array}{c}\hfill {\lambda }_1={\displaystyle \frac{1}{{\lambda }^2}},{\lambda }_2={\lambda }_3=\lambda ={\displaystyle \frac{r}{R}}\end{array}$$

(53)

and, we have, $u_1=E_R$, $u_2=E_{\mathsf{\Theta }}$ and $u_3=E_{\mathsf{\Phi }}$, where $\{E_R,E_{\mathsf{\Theta }},E_{\mathsf{\Phi }}\}$ is the spherical polar coordinate basis for the reference configuration.

Here, we use the residual stress

$$\begin{array}{c}\hfill S_{\mathrm{R}}=s_1\left(R\right)E_R\otimes E_R+s_2\left(R\right)E_{\mathsf{\Theta }}\otimes E_{\mathsf{\Theta }}+s_2\left(R\right)E_{\mathsf{\Phi }}\otimes E_{\mathsf{\Phi }},\end{array}$$

(54)

where

$$\begin{array}{c}\hfill s_1\left(R\right)=\kappa (R-R_a)(R-R_b)\end{array}$$

(55)

and $\kappa$ has the dimension kPa/m². The equilibrium equation requires

$$\begin{array}{c}\hfill s_2\left(R\right)={\displaystyle \frac{1}{2R}}\frac{\mathrm{d}\left(R^2{s}_1\left(R\right)\right)}{\mathrm{d}R}.\end{array}$$

(56)

The free stress surface (2) is clearly satisfied.

We only discuss, for simplicity, the case when $a=E_R$ and, we simply have, for the non-stretch invariants

$$\begin{array}{c}\hfill {\alpha }_1=1,{\alpha }_2={\alpha }_3=0,{\zeta }_1=s_1,{\zeta }_2={\zeta }_3=s_2.\end{array}$$

(57)

The non-zero Cauchy stress components (with respect to the spherical coordinate system) are ${\tau }_{rr}$, ${\tau }_{\theta \theta }$ and ${\tau }_{\varphi \varphi }$.

Since ${\zeta }_2={\zeta }_3$, from (15), we have,

$$\begin{array}{c}\hfill {\tau }_{\theta \theta }={\tau }_{\varphi \varphi }.\end{array}$$

(58)

Hence, the equilibrium equation reduces to

$$\begin{array}{c}\hfill \frac{\mathrm{d}{\tau }_{rr}}{\mathrm{d}r}+{\displaystyle \frac{2}{r}}({\tau }_{rr}-{\tau }_{\theta \theta })=0.\end{array}$$

(59)

Integrating (59) and taking account that if we assume that radial stress vanishes at $r=a$, we get

$$\begin{array}{c}\hfill {\tau }_{rr}={\int }_a^r{\displaystyle \frac{2}{r}}\left(\lambda {\displaystyle \frac{\partial \tilde{\mathsf{\Omega }}}{\partial {\lambda }_2}}-{\displaystyle \frac{1}{{\lambda }^2}}{\displaystyle \frac{\partial \tilde{\mathsf{\Omega }}}{\partial {\lambda }_1}}\right)\mathrm{d}r.\end{array}$$

(60)

The radial stress is depicted in Figure 4 for $a/R_a=0.95$ (the sphere is moved inwards). We only consider $R_b/R_a=1.182$. While the fiber is in tension with ${\displaystyle {\displaystyle \frac{b-a}{R_b-R_a}}=1.082}$, a negative radial stress is obtained; taking note that ${\displaystyle {\displaystyle \frac{b}{R_b}}=0.970}$. Figure 4 indicates that the absolute value of the radial stress ${\tau }_{rr}$ increases in the present of a residual stress. However, when the sphere is moved outwards (see Figure 5), we have, $a/R_a=1.1$, ${\displaystyle {\displaystyle \frac{b}{R_b}}=1.062787231}$, ${\displaystyle {\displaystyle \frac{b-a}{R_b-R_a}}=0.8583214670}$ and hence the fiber is in compression. We note that the radial stress is negative in some segments due to the presence of residual stress. In this section the values ${\mu }_1=\beta =0$ are used for the graphs.

7. Conclusions

In this communication, we propose a novel separable spectral strain energy function for RSPDs which contains single-variable spectral-invariant functions. The use of spectral invariants is the state of the art in modelling RSPDs and has the advantage that each of the spectral invariants has a clear physical interpretation, which is useful in aiding the design of experiments. The spectral invariants depend on the right stretch tensor, the preferred direction and the residual stress tensor. Our approach ensures that the strain energy function for infinitesimal strain can be simply obtained from its finite strain counterpart and vice-versa. The analysis of the boundary value problems given in Section 6, highlights the simplicity of the spectral approach. Modelling full or partial exclusion of compressed fibers is simply done by amending the strain energy function. In future, we require experimental data to compare our theory with different types of RSPDs. Since classical invariants can be explicitly expressed in terms of spectral invariants but not vice versa, spectral formulations are, hence, more general.