Likely cavitation in stochastic elasticity
L. Angela Mihai Danielle Fitt Thomas E. Woolley Alain Goriely


	
Likely cavitation in stochastic elasticity
    L. Angela Mihai[School of Mathematics, Cardiff University, Senghennydd Road, Cardiff, CF24 4AG, UK, Email: ]
	    Danielle Fitt[School of Mathematics, Cardiff University, Senghennydd Road, Cardiff, CF24 4AG, UK, Email: ]
	    Thomas E. Woolley[School of Mathematics, Cardiff University, Senghennydd Road, Cardiff, CF24 4AG, UK, Email: ]
	    Alain Goriely[Mathematical Institute, University of Oxford, Woodstock Road, Oxford, OX2 6GG, UK, Email: ]

    October 13, 2018
====================================================================================================================================================================================================================================================================================================================================================================================================================================================================




We revisit the classic problem of elastic cavitation within the framework of stochastic elasticity. For the deterministic elastic problem, involving homogeneous isotropic incompressible hyperelastic spheres under radially symmetric tension, there is a critical dead-load traction at which cavitation can occur for some materials. In addition to the well-known case of stable cavitation post-bifurcation at the critical dead load, we show the existence of unstable snap cavitation for some isotropic materials satisfying Baker-Ericksen inequalities. For the stochastic problem, we derive the probability distribution of the deformations after bifurcation. In this case, we find that, due to the probabilistic nature of the material parameters, there is always a competition between the stable and unstable states. Therefore, at a critical load, stable or unstable cavitation occurs with a given probability, and there is also a probability that the cavity may form under smaller or greater loads than the expected critical value. We refer to these phenomena as `likely cavitation'. Moreover, we provide examples of homogeneous isotropic incompressible materials exhibiting stable or unstable cavitation together with their stochastic equivalent.  

	
Key words: stochastic hyperelastic models, stable or unstable cavitation, isotropic incompressible spheres, Baker-Ericksen inequalities, dead-load traction, probability.







§ INTRODUCTION


Experiments carried out by Gent and Lindley in 1958 <cit.>, on rubber cylinders, revealed that some materials can rupture under relatively small tensile dead loads by opening an internal cavity. Following this work, the theoretical analysis of Ball (1982) <cit.> provided an explanation for the formation of a spherical cavity at the centre of a sphere of isotropic hyperelastic incompressible material in radially symmetric tension under prescribed surface displacements or dead loads. There, the word `cavitation' was used to describe such void-formation within a solid by analogy to the similar phenomenon observed in fluids. As cavitation in solids is an inherently nonlinear mechanical effect, not captured by the linear elasticity theory, many different studies have been devoted to the modelling of this effect within the finite elasticity framework. For instance: spheres of particular homogeneous isotropic incompressible materials were discussed in <cit.>; homogeneous anisotropic spheres with transverse isotropy about the radial direction were examined in <cit.>; concentric homogeneous spheres of different hyperelastic material were analysed in <cit.>; non-spherical cavities were investigated in <cit.>; cavities with non-zero pressure were presented in <cit.>; cavitation in an elastic membrane was studied in <cit.>; the homogenisation problem of nonlinear elastic materials was treated in <cit.>; growth-induced cavitation in nonlinearly elastic solids was explored in <cit.>. Recent experimental results on the onset, healing and growth of cavities in elastomers were reported in <cit.>. For many other results on cavitation in solids, we refer to the review articles <cit.>, focusing on rubberlike materials, and the references therein.

The present work focuses on the phenomenon of “cavitation” contained within the theoretical context of finite elastostatics. Finite elasticity theory covers the simplest case where internal forces only depend on the current deformation of the material and not on its history, and is based on average data values. Within this framework, hyperelastic materials are the class of material models described by a strain-energy function characterised by a set of deterministic model parameters. In addition, for solid elastic materials, uncertainties in the observational data generally arise from the inherent variation in material properties and testing protocols <cit.>. In view of these uncertainties, recently, stochastic representations of isotropic incompressible hyperelastic materials characterised by a stochastic strain-energy function, for which the model parameters are random variables following standard probability laws, were proposed in <cit.>, while compressible versions of these models were constructed in <cit.>. Ogden-type stochastic strain-energy functions were calibrated to experimental data for rubber and soft tissue materials in <cit.>, and anisotropic stochastic models with the model parameters as spatially-dependent random field variables were calibrated to vascular tissue data in <cit.>. These models employ the maximum entropy principle for a discrete probability distribution introduced by Jaynes (1957) <cit.> and based on the notion of  entropy (or uncertainty) defined by Shannon (1948) <cit.>. Such models can be useful for stochastic finite element implementations <cit.>.

For stochastic hyperelastic models, the immediate question is: what is the influence of the random model parameters on the predicted nonlinear elastic responses? This question was previously considered by us in <cit.>, for the stochastic Rivlin cube, and in <cit.>, for the symmetric inflation of internally pressurised stochastic spherical shells and tubes. These idealised problems illustrate some important effects on the likely elastic responses of stochastic hyperelastic materials under large strains.

Here, we address this question by employing a similar approach as in <cit.> to revisit, in the context of stochastic elasticity, the cavitation problem of incompressible spheres of stochastic homogeneous isotropic hyperelastic materials under uniform radial tensile dead loads. Moreover, for all homogeneous isotropic hyperelastic models considered so far in the literature, cavitation appears as a supercritical bifurcation, where typically, after bifurcation, the cavity radius monotonically increases as the applied load increases (see, e.g., <cit.>). However, as we demonstrate here, the usual restriction that a material satisfies the Baker-Ericksen (BE) inequalities <cit.> is not sufficient to exclude the possibility of a subcritical bifurcation. In this case, one expects a snap cavitation for which there is a jump in the radius of the cavity immediately after bifurcation. Indeed, we obtain the general conditions under which a cavitation can appear through a supercritical or subcritical bifurcation and construct explicitly, for the first time, examples of isotropic incompressible hyperelastic models that exhibit snap cavitation. The stochastic version of these models are then explored. In this case, we find that, due to the probabilistic nature of the model parameters, supercritical or subcritical bifurcation occurs with a given probability, and there is also a probability that the cavity may form under smaller or greater loads that the expected critical value. We refer to these phenomena as `likely cavitation'.

We begin, in Section <ref>, with a detailed presentation of the stochastic elastic framework. Then, in Section <ref>, for the stochastic sphere, after we review the elastic solution to the cavitation problem under uniformly applied tensile dead load, we recast the problem in the stochastic setting, and find the probabilistic solution. Concluding remarks are provided in Section <ref>.






§ STOCHASTIC ISOTROPIC HYPERELASTIC MODELS


We recall that a homogeneous hyperelastic model is described by a strain-energy function W(𝐅) that depends on the deformation gradient tensor, 𝐅, with respect to a fixed reference configuration, and is characterised by a set of deterministic model parameters <cit.>. In contrast, a stochastic homogeneous hyperelastic model is defined by a stochastic strain-energy function, for which the model parameters are random variables that satisfy standard probability laws <cit.>. In this case, each model parameter is described in terms of its mean value and its variance, which contains information about the range of values about the mean value. While it is rarely possible if ever to obtain complete information about a random quantity in an elastic sample of material, the partial information provided by the mean value and the variance is the most commonly used in many practical applications <cit.>. Here, we combine finite elasticity  and information theory, and rely on the following general hypotheses <cit.>:


  (A1) Material objectivity: The principle of material objectivity (frame indifference) states that constitutive equations must be invariant under changes of frame of reference. It requires that the scalar strain-energy function, W=W(F), depending only on the deformation gradient F, with respect to the reference configuration, is unaffected by a superimposed rigid-body transformation (which involves a change of position) after deformation, i.e., W(R^TF)=W(F), where R∈ SO(3) is a proper orthogonal tensor (rotation). Material objectivity is guaranteed by considering strain-energy functions defined in terms of invariants.
	

  (A2) Material isotropy: The principle of isotropy requires that the strain-energy function is unaffected by a superimposed rigid-body transformation prior to deformation, i.e., W(FQ)=W(F), where Q∈ SO(3). For isotropic materials, the  strain-energy  function is a symmetric function of the  principal stretches {λ_i}_i=1,2,3 of F, i.e., W(F)=𝒲(λ_1,λ_2,λ_3).
	

  (A3) Baker-Ericksen inequalities: In addition to the fundamental principles of objectivity and material symmetry, in order for the behaviour of a hyperelastic material to be physically realistic, there are some universally accepted constraints on the constitutive equations. Specifically, for a hyperelastic body, the Baker-Ericksen (BE) inequalities, which state that the greater principal (Cauchy) stress occurs in the direction of the greater principal stretch, are <cit.>:

    (T_i-T_j)(λ_i-λ_j)>0    λ_i≠λ_j,   i,j=1,2,3,

where {λ_i}_i=1,2,3 and {T_i}_i=1,2,3 denote the principal stretches and the principal Cauchy stresses, respectively, and “≥” replaces the strict inequality “>” if any two principal stretches are equal. The BE inequalities (<ref>) take the equivalent form <cit.>

    (λ_i∂𝒲/∂λ_i-λ_j∂𝒲/∂λ_j)(λ_i-λ_j)>0    λ_i≠λ_j,   i,j=1,2,3,

where the strict inequality “>” is replaced by “≥” if any two principal stretches are equal.


  (A4) Finite mean and variance for the random shear modulus: We assume that, for any given finite deformation, the random shear modulus, μ, and its inverse, 1/μ, are second-order random variables, i.e., they have finite mean value and finite variance <cit.>.

While (A4) contains physically realistic expectations on the random shear modulus, which will be drawn from a probability distribution, assumptions (A1)-(A3) are well-known principles in isotropic finite elasticity <cit.>.

Specifically, we focus our attention on homogeneous incompressible hyperelastic materials characterised by the following stochastic strain-energy function <cit.>,

    𝒲(λ_1,λ_2,λ_3)=μ_1/2m^2(λ_1^2m+λ_2^2m+λ_3^2m-3)
    	+μ_2/2n^2(λ_1^2n+λ_2^2n+λ_3^2n-3),

where m and n are deterministic constants, and μ_1 and μ_2 are random variables following given probability distributions. In the deterministic elastic case,  μ_1, μ_2, m and n are constants, and the model contains, as special cases, the neo-Hookean model, the Mooney-Rivlin model, and the one- and two-term Ogden models. In both the deterministic elastic and stochastic cases, the shear modulus for infinitesimal deformations of these models is defined as μ=μ_1+μ_2 <cit.>. Note that we could easily extend our description to include m and n as stochastic variables as well. However, increasing the complexity in this way is not relevant for the present discussion. Including additional sources of randomness is an avenue of future research.

As it is well known, the deformation of an homogeneous isotropic hyperelastic material under uniaxial tension is a simple extension in the direction of the tensile force if and only if the BE inequalities hold <cit.>. Under these conditions, the shear modulus is positive, but the individual coefficients may be either positive or negative, allowing for some interesting nonlinear elastic effects to be captured (see <cit.> and the references therein). In particular, in the present paper, the initiation of either stable or unstable snap cavitation in a homogeneous isotropic sphere will be presented.

For the stochastic materials described by (<ref>), condition (A4) is guaranteed by setting the following mathematical expectations <cit.>:

    E[μ]=μ>0,   
    
    		E[log μ]=ν,   .

Then, under the constraints (<ref>), the random shear modulus, μ, with mean value μ and standard deviation μ=√(Var[μ]), defined as the square root of the variance, Var[μ], follows a Gamma probability distribution <cit.>, with hyperparameters ρ_1>0 and ρ_2>0 satisfying

    μ=ρ_1ρ_2,    μ=√(ρ_1)ρ_2.

The corresponding probability density function takes the form <cit.>

    g(μ;ρ_1,ρ_2)=μ^ρ_1-1e^-μ/ρ_2/ρ_2^ρ_1Γ(ρ_1),     μ>0  ρ_1, ρ_2>0,

where Γ:ℝ^*_+→ℝ is the complete Gamma function

    Γ(z)=∫_0^+∞t^z-1e^-tdt.

For technical convenience, we set a finite constant value b>-∞, such that μ_i>b, i=1,2 (e.g., b=0 if μ_1>0 and μ_2>0, but b is not unique in general), and introduce the auxiliary random variable <cit.>

    R_1=μ_1-b/μ-2b,

such that 0<R_1<1. Consequently, we can equivalently express the random model parameters μ_1 and μ_2 as follows,

    μ_1=R_1(μ-2b)+b,    μ_2=μ-μ_1=(1-R_1)(μ-2b)+b.

It is reasonable to assume <cit.>

    E[log R_1]=ν_1,   ,
    
    E[log(1-R_1)]=ν_2,   ,

in which case, the random variable R_1, with mean value R_1 and variance Var[R_1], follows a standard Beta distribution <cit.>, with hyperparameters ξ_1>0 and ξ_2>0 satisfying

    R_1=ξ_1/ξ_1+ξ_2,    Var[R_1]=ξ_1ξ_2/(ξ_1+ξ_2)^2(ξ_1+ξ_2+1).

The associated probability density function is

    β(r;ξ_1,ξ_2)=r^ξ_1-1(1-r)^ξ_2-1/B(ξ_1,ξ_2),         r∈(0,1)  ξ_1, ξ_2>0,

where B:ℝ^*_+×ℝ^*_+→ℝ is the Beta function

    B(x,y)=∫_0^1t^x-1(1-t)^y-1dt.

Thus, for the random coefficients given by (<ref>), the corresponding mean values take the form,

    μ_1=R_1(μ-2b)+b,
        μ_2=μ-μ_1=(1-R_1)(μ-2b)+b,

and the variances and covariance are, respectively,

    Var[μ_1]
    =(μ-2b)^2Var[R_1]+(R_1)^2Var[μ]+Var[μ]Var[R_1],
          Var[μ_2]
    =(μ-2b)^2Var[R_1]+(1-R_1)^2Var[μ]+Var[μ]Var[R_1],
          Cov[μ_1,μ_2]=1/2(Var[μ]-Var[μ_1]-Var[μ_2]).

It should be noted that the random variables μ and R_1 are independent, depending on parameters (ρ_1,ρ_2) and (ζ_1,ζ_2), respectively, which are derived by fitting distributions to given data. However, μ_1 and μ_2 are dependent variables as they both require (μ,R_1) to be defined. Explicit derivations of the probability distributions for the random parameters when stochastic isotropic hyperelastic models are calibrated to experimental data are presented in <cit.>.

Our aim here is to analyse the radially symmetric finite deformations of a sphere of stochastic hyperelastic material defined by (<ref>), under tension, when subject to prescribed surface dead loads applied uniformly in the radial direction. One can view the stochastic sphere as an ensemble (or population) of  spheres, where each sphere has the same initial radius and is made from a homogeneous isotropic incompressible hyperelastic material, with the elastic parameters not known with certainty, but drawn from known probability distributions. Then, for every hyperelastic sphere, the finite elasticity theory applies. For the stochastic hyperelastic body, the question is: what is the probability distribution of stable radially symmetric deformation under a given surface dead load?






§ INCOMPRESSIBLE SPHERES


In this section, we consider a sphere of stochastic incompressible hyperelastic material described by (<ref>), subject to a radially symmetric deformation, caused by the sole action of a given radial tensile dead load. As for the deterministic elastic sphere <cit.>, we obtain conditions on the constitutive law, such that, setting the internal pressure equal to zero, where the radius tends to zero, the required external dead load is finite, and therefore cavitation occurs. We further analyse the stability of the cavitated solution, and distinguish between supercritical cavitation, where the cavity radius monotonically increases as the dead load increases, and subcritical (snap) cavitation, with a sudden jump to a finite internal radius immediately after initiation. To the best of our knowledge, in the deterministic elastic case, the onset of snap cavitation  in a homogeneous isotropic sphere has not been discussed before. Therefore, we start our analysis in the deterministic elastic context before  extending it to the stochastic case.

For the stochastic sphere, the radially symmetric deformation takes the form

    r=g(R),    θ=Θ,    ϕ=Φ,

where (R,Θ,Φ) and (r,θ,ϕ) are the spherical polar coordinates in the reference and current configuration, respectively, such that 0≤ R≤ B, and g(R)≥ 0 is to be determined. The corresponding deformation gradient is equal to F=diag(λ_1,λ_2,λ_3), with

    λ_1=dg/dR=λ^-2,    λ_2=λ_3=g(R)/R=λ,

where λ_1 and  λ_2=λ_3 are the radial and hoop stretches, respectively, and dg/dR denotes the derivative of g with respect to R. By (<ref>),

    g^2dg/dR=R^2,

hence,

    g(R)=(R^3+c^3)^1/3,

where c≥ 0 is a constant to be calculated. If c>0, then g(R)→ c>0 as R→ 0_+, and a spherical cavity of radius c forms at the centre of the sphere, from zero initial radius (see Figure <ref>), otherwise the sphere remains undeformed.





Assuming that the deformation (<ref>) is due to a prescribed radial tensile dead load, applied uniformly on the sphere surface in the reference configuration, in the absence of body forces, the radial equation of equilibrium is

    dP_11/dR+2/R(P_11-P_22)=0,

or equivalently,

    dP_11/dλλ^-2+2P_11-P_22/1-λ^3=0,

where P=(P_ij)_i,j=1,2,3 is the first Piola-Kirchhoff stress tensor. For an incompressible material,

    P_11=∂𝒲/∂λ_1-p/λ_1,    
    P_22=∂𝒲/∂λ_2-p/λ_2.

Denoting

    W(λ)=𝒲(λ^-2,λ,λ),

where λ=r/R=g(R)/R=(1+c^3/R^3)^1/3>1, we obtain

    dW/dλ=-2/λ^3∂𝒲/∂λ_1+2∂𝒲/∂λ_2
    =-2P_11/λ^3+2P_22.

Then, setting the internal pressure (at R→ 0_+) equal to zero, by (<ref>) and (<ref>), the external tension (at R=B) is equal to

    T=P_11/λ^2|_λ=λ_b.=∫_λ_b^λ_cdW/dλdλ/λ^3-1,

and the applied dead load, in the reference configuration, is

    P=Tλ_b^2=λ_b^2∫_λ_b^λ_cdW/dλdλ/λ^3-1,

where λ_c and λ_b represent the stretches at the centre and outer surface, respectively. The value of the required dead load, P_0, for the onset of cavitation (bifurcation from the reference state) is obtained by taking λ_c→∞ and λ_b=(1+c^3/B^3)^1/3→ 1 as c→ 0_+ in (<ref>), i.e.,

    P_0=∫_1^∞dW/dλdλ/λ^3-1.

The BE inequalities (<ref>) imply

    dW/dλ1/λ^3-1>0,

hence, P_0>0. Then, if the critical dead load given by (<ref>) is finite, cavitation takes place, else, the sphere remains undeformed.

Before considering the stochastic setting, we briefly revisit the deterministic elastic case.




 §.§ Deterministic elastic spheres


For a sphere made of a hyperelastic material with the strain-energy function

    𝒲(λ_1,λ_2,λ_3)=μ_1/2m^2(λ_1^2m+λ_2^2m+λ_3^2m-3)
    +μ_2/2n^2(λ_1^2n+λ_2^2n+λ_3^2n-3),

where μ_1 and μ_2 are positive constants,  (<ref>) takes the form

    W(λ)=μ_1/2m^2(λ^-4m+2λ^2m-3)+μ_2/2n^2(λ^-4n+2λ^2n-3).

For the onset of cavitation, the critical dead load traction, defined by (<ref>), is equal to

    P_0=2μ_1/m∫_1^∞λ^2m-1-λ^-4m-1/λ^3-1dλ+2μ_2/n∫_1^∞λ^2n-1-λ^-4n-1/λ^3-1dλ,

or equivalently, by the change of variable x=λ^3-1,

    P_0   =2μ_1/3m∫_0^∞(x+1)^(2m-3)/3-(x+1)^-(4m+3)/3/xdx
       +2μ_2/3n∫_0^∞(x+1)^(2n-3)/3-(x+1)^-(4n+3)/3/xdx.

By (<ref>), P_0 is finite, and hence, a spherical cavity forms, if and only if the following conditions are simultaneously satisfied:  2m-3<0, -4m-3<0, 2n-3<0, -4n-3<0, or equivalently  <cit.> (see also Example 5.1 of <cit.>), if and only if

    -3/4< m, n<3/2.

In particular, cavitation is found in a neo-Hookean sphere (with m=1 and n=0), but not in a Mooney-Rivlin sphere (with m=1 and n=-1). The special cases when m∈{-1/2, 1} and n=0, are given as examples in <cit.>, and when m∈{1/2, 3/4, 1, 5/4} and n=0, the explicit critical loads are provided in <cit.>. When these bounds and the BE inequalities are satisfied, the critical pressure P_0 is finite and the problem is to find the behavior of the cavity in a neighborhood of this critical value. In each of those previously studied cases (see, e.g., Figure 2 of <cit.>), cavitation forms from zero radius and then presents itself  as a supercritical bifurcation with stable cavitation (i.e. the new bifurcated solution exists locally for values of P>P_0, and the radius of the cavity monotonically increases with the applied load post-bifurcation).

Another theoretical possibility is that the bifurcation could be subcritical (i.e., the cavitated solution exists locally for values less than P_0 and is unstable). One would then expect an unstable snap cavitation with a sudden jump to a cavitated solution with a finite internal radius. This subcritical behaviour of the homogeneous isotropic elastic sphere has not been explicitly demonstrated in the literature before. Here, we show that, depending on the model parameters, the family of materials (<ref>) can exhibit both behaviours. General conditions for a given material to exhibit either a subcritical or supercritical bifurcation are provided in Appendix <ref>. 





As an example, we illustrate the variety of behaviours when m=1 and n=-1/2, such that (<ref>) takes the form

    W(λ)=μ_1/2(λ^-4+2λ^2-3)+2μ_2(λ^2+2λ^-1-3).

In this case, under the deformation (<ref>), the BE inequalities (<ref>) are reduced to

    μ_1+2μ_2λ^3/1+λ^3>0.

The inequality (<ref>) implies that, when λ→ 1, the  shear modulus must be positive, i.e., μ=μ_1+μ_2>0, while if λ→∞, then μ_1+2μ_2>0. Noting that the function of λ on the left-hand side is monotonically increasing when μ_2 is positive, and decreasing if μ_2 is negative, and taking μ_1>0, the two limits imply that the BE inequalities are satisfied for all values of λ if

    0<μ_1/μ<2.

For sufficiently small c/B, the corresponding dead-load traction, defined by (<ref>), is equal to

    P   =2μ_1[(1+c^3/B^3)^1/3+1/4(1+c^3/B^3)^-2/3]+4μ_2(1+c^3/B^3)^1/3
       =4μ(1+c^3/B^3)^1/3-2μ_1[(1+c^3/B^3)^1/3-1/4(1+c^3/B^3)^-2/3].

Then, the critical tensile dead load given by (<ref>) takes the form

    P_0=4μ-3μ_1/2.

As P_0 is positive in (<ref>), we have 0<μ_1/μ<8/3, which is guaranteed by (<ref>).

The question is now to find the possible behaviour of the cavity opening c as a function of P in a neighbourhood of P_0. On differentiating (<ref>) with respect to c/B, we obtain

    dP/d(c/B)=2c^2/B^2{2μ(1+c^3/B^3)^-2/3-μ_1[(1+c^3/B^3)^-2/3+1/2(1+c^3/B^3)^-5/3]}.

Hence, by Proposition <ref> given in Appendix <ref> (with n=3), when

    0<μ_1/μ<4/3=inf_0<c/B<1[2(1+c^3/B^3)(3/2+c^3/B^3)^-1],

where “inf” denotes infimum, the bifurcation is supercritical and the radius of the cavity monotonically increases as the tensile dead load increases. However, if there exists c_0>0, such that

    2(1+c_0^3/B^3)(3/2+c_0^3/B^3)^-1<μ_1/μ<2,

then the bifurcation is subcritical and the required applied load starts to decrease at c=c_0, where there is a sudden jump in the opening of cavity. In particular, if (<ref>) holds for c_0=0, i.e.,

    4/3<μ_1/μ<2,

then (<ref>) is valid while the cavitation becomes unstable.

Thus, dP/d(c/B)→ 0 as c→ 0_+, and by Proposition <ref>, the bifurcation at the critical load, P_0, is supercritical (respectively, subcritical) if dP/d(c/B)>0 (respectively, dP/d(c/B)<0) for arbitrarily small c/B. Examples of both these behaviours are illustrated in Figures <ref> and <ref>.





 §.§ Stochastic elastic spheres


We now turn our attention to the stochastic model described by (<ref>), with m=1 and n=-1/2, and the other parameters drawn from probability distributions. In this case, recalling that μ follows a Gamma distribution g(u;ρ_1,ρ_2), defined by (<ref>), the probability distribution of stable cavitation is equal to

    P_1(μ_1)=1-∫_0^3μ_1/4g(u;ρ_1,ρ_2)du,

and that of unstable cavitation is

    P_2(μ_1)=1-P_1(μ_1).















For example, taking ρ_1=405 and  ρ_2= 0.01 (see Figure <ref>), the mean value of the shear modulus is μ=ρ_1ρ_2=4.05, and the probability distributions given by equations (<ref>)-(<ref>) are illustrated numerically in Figure <ref> (with blue lines for P_1 and red lines for P_2). In this case, if μ_1=5<5.4=4μ/3 say, then stable cavitation is expected, but there is also about 10% chance that unstable snap cavitation occurs. Similarly, when 4μ/3=5.4<μ_1=5.8<8.1=2μ, unstable cavitation is expected, but there is also about 10% chance that the cavitation is stable. Stable and unstable cavitation of a stochastic sphere are illustrated numerically in Figure <ref>. Specifically:


  (a) In Figure <ref>(a), b=0 in (<ref>), and the random variable R_1=μ_1/μ is drawn from a Beta distribution with ξ_1=287 and ξ_2=36. In this case, μ_1=3.6<5.4=4μ/3, and stable cavitation, with supercritical bifurcation after the spherical cavity opens, is expected.


  (b) In Figure <ref>(b), b=-3 in (<ref>), and the random variable R_1=(μ_1+3)/(μ+6) draws its values from a Beta distribution with ξ_1=325 and ξ_2=10. Thus, 4μ/3=5.4<μ_1=6.75<8.1=2μ, and unstable cavitation, with subcritical bifurcation after the spherical cavity forms, is expected.


For the numerical examples shown in Figure <ref> also, the critical dead load is P_0=4μ-3μ_1/2, as given by (<ref>), with μ and μ_1 following probability distributions. In each case, the expectation is that the onset of cavitation occurs at the mean value P_0=4μ-3μ_1/2, found at the intersection of the dashed black line with the horizontal axis. However, there is a chance that cavity can form under smaller or greater critical loads that the expected load value, as shown by the coloured interval about the mean value along the horizontal axis.

To summarise, for a stochastic elastic sphere under uniform tensile dead load, we obtain the probabilities of stable or unstable cavitation, given that the material parameters are generated from known probability density functions. In the deterministic elastic case, there is a single critical parameter value that strictly separates the cases where the initiation of either stable or unstable cavitation occurs. By contrast, in the stochastic case, there is a probabilistic interval, containing the deterministic critical value, where there is always a competition between the stable and unstable states in the sense that both have a quantifiable chance to be found. For the onset of cavitation, there is also a probabilistic interval where a cavity may form, with a given probability, under smaller or greater loads that the expected critical value.






§ CONCLUSION


This work is motivated by the fact that a crucial part in assessing the physical properties of many solid materials is to quantify the uncertainties in their mechanical responses, which cannot be ignored. In particular, the idealised problem of the formation of a spherical cavity at the centre of a solid sphere illustrates some important effects on the likely elastic responses of stochastic hyperelastic materials under large strains.

For homogeneous isotropic incompressible spheres of stochastic hyperelastic material, subject to radial tensile dead loads applied uniformly on the sphere surface, we examined the possible radially symmetric deformations and determined which of these deformations are stable. Homogeneous stochastic hyperelastic material models satisfying certain theoretical assumptions were recently introduced to capture the dispersion in experimental data in addition to the traditional mean-data values <cit.>.

For the deterministic elastic problem, where the model parameters are single-valued constants, non-trivial deformations, whereby a spherical cavity forms at the centre of the sphere, are possible for some classes of materials when the applied tensile dead loads are sufficiently large <cit.>. In some materials, cavitation is stable, in the sense that the cavity radius monotonically increases as the applied dead load increases <cit.>. Here, we showed that a sudden jump in the cavity opening, causing unstable snap cavitation, at the critical dead load can also occur in a homogeneous isotropic incompressible sphere, provided that the material satisfies Baker-Ericksen inequalities. If such a material could be found, a sphere made of this material would suddenly increase its volume at a critical load and show some form of hysteresis as the load is removed.

In the stochastic case, the probabilistic nature of the solution reflects the probability in the constitutive law, and bifurcation and stability can be quantified in terms of probabilities. By contrast to the deterministic elastic problem, where deterministic critical parameter values strictly separate the cases where either the stable or unstable cavitation occurs, for the stochastic problem, we obtained probabilistic intervals where both states have a quantifiable chance to exist. For the onset of cavitation, there is a probabilistic interval where the cavity may form, with a given probability, under smaller or greater loads that the expected critical value.

As a direct application of our approach, one could consider the cavitation of an inhomogeneous sphere made of concentric homogeneous spheres of different stochastic material, similar to the concentric homogeneous spheres of deterministic elastic material treated in <cit.> and <cit.>. Such composite spheres would require comparing both ensemble and spatial averages.







§ STABILITY ANALYSIS


We provide a corrected version of Proposition 5.2 of <cit.> and its proof. In particular, we show that, in the deterministic elastic case, both subcritical and supercritical behaviours close to the bifurcation are possible, depending on the material.


Let W(λ) be twice differentiable at λ=1, and

    P(c)=(1+c^n)^(n-1)/n∫_(1+c^n)^1/n^∞dW/dλdλ/λ^n-1,

where n>1. Then lim_c→0_+(dP/dc)=0, and if

    lim_c→0_+(P(c)-lim_λ→11/n(n-1)d^2W/dλ^2)>0,

then dP/dc>0 for sufficiently small c>0 (i.e., the bifurcation is supercritical), while if

    lim_c→0_+(P(c)-lim_λ→11/n(n-1)d^2W/dλ^2)<0,

then dP/dc<0 for sufficiently small c>0 (i.e., the bifurcation is subcritical).

These cases are illustrated, for the particular example of material presented in this paper, in Figures <ref> and <ref>.



  
Proof. We denote θ=(1+c^n)^(n-1)/n and define P(θ)=P(c). Then

    P(θ)=θ∫_θ^1/(n-1)^∞dW/dλdλ/λ^n-1

and

    dP/dc=dP/dθdθ/dc,

where

    dP/dθ=∫_θ^1/(n-1)^∞dW/dλdλ/λ^n-1-θ^1/(n-1)/n-1(dW/dλ1/λ^n-1)|_λ=θ^1/(n-1)..

It follows that

    lim_θ→1dP/dθ   =lim_θ→1∫_θ^1/(n-1)^∞dW/dλdλ/λ^n-1-lim_θ→1θ^1/(n-1)/n-1dW/dλ1/λ^n-1|_λ=θ^1/(n-1).
       =lim_c→0_+P(c)-lim_θ→1θ^1/(n-1)/n-1(dW/dλ1/λ^n-1)|_λ=θ^1/(n-1).
       =lim_c→0_+P(c)-lim_λ→11/n(n-1)dW/dλ1/λ-1
       =lim_c→0_+P(c)-lim_λ→11/n(n-1)d^2W/dλ^2.

Hence,

    lim_c→0_+dP/dc=lim_θ→1dP/dθ=lim_c→0_+P(c)-lim_λ→11/n(n-1)d^2W/dλ^2

and dP/dc>0 (respectively, dP/dc<0) for sufficiently small c>0 if and only if (<ref>) (respectively, (<ref>)) holds. This concludes the proof.


Note that the difference between this result and Proposition 5.2 of <cit.> comes from the (correct) minus sign between the two terms on the right-hand side of (<ref>) (whereas a plus sign is found in the corresponding unlabelled expression appearing between Equations (5.25) and (5.26) of <cit.>).











  
Acknowledgement. We thank  John Ball for a discussion on the corrected version of Proposition 5.2 of <cit.>, as presented here in Appendix <ref>. The support for Alain Goriely by the Engineering and Physical Sciences Research Council of Great Britain under research grant EP/R020205/1 is gratefully acknowledged.




9
	
Abramowitz:1964 Abramowitz M, Stegun IA. 1964. Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, National Bureau of Standards, Applied Mathematics Series, vol. 55, Washington.

Antman:1987:AN Antman SS, Negro n-Marrero PV. 1987. The remarkable nature of radially symmetric equilibrium states of aleotropic nonlinearly elastic bodies, Journal of Elasticity 18, 131-164.

ArreguiMena:2016:AMMM Arregui-Mena JD, Margetts L, Mummery PM. 2016. Practical application of the stochastic finite element method, Archives of Computational Methods in Engineering 23, 171-190.

Babuska:2005:BTZ Babuška I, Tempone R, Zouraris GE. 2005. Solving elliptic boundary value problems with uncertain coefficients by the finite element method: the stochastic formulation, Computer Methods in Applied Mechanics and Engineering 194, 1251-1294.

BakerEricksen:1954 Baker M, Ericksen JL. 1954. Inequalities restricting the form of stress-deformation relations for isotropic elastic solids and Reiner-Rivlin fluids, Journal of the Washington Academy of Sciences 44, 24-27.

Ball:1982 Ball JM. 1982. Discontinuous equilibrium solutions and cavitation in nonlinear elasticity, Philosophical Transactions of the Royal Society A 306, 557-611.

Bayes:1763 Bayes T. 1763. An essay towards solving a problem in the doctrine of chances, Philosophical Transactions 53, 370-418.

ChouWang:1989:CWH Chou-Wang M-S, Horgan CO. 1989. Void nucleation and growth for a class of incompressible nonfinearly elastic materials, International Journal of Solids and Structures 25, 1239-1254.

Caylak:2018:etal Caylak I, Penner E, Dridger A, Mahnken R. 2018. Stochastic hyperelastic modeling considering dependency of material parameters, Computational Mechanics (doi: 10.1007/s00466-018-1563-z).

Farmer:2017 Farmer CL. 2017. Uncertainty quantification and optimal decisions, Proceedings of the Royal Society A 473, 20170115.

Fond:2001 Fond, C. 2001. Cavitation criterion for rubber materials: a review of void-growth models, Journal of Polymer Science: Part B 39, 2081-2096.

Gent:1991 Gent AN. 1991. Cavitation in rubber: a cautionary tale, Rubber Chemistry and Technology 63, G49-G53.

Gent:1959:GL Gent AN, Lindley PB. 1959. Internal rupture of bonded rubber cylinders in tension, Proceedings of the Royal Society of London A 249, 195-205.

goriely17 Goriely A. 2017. The Mathematics and Mechanics of Biological Growth, Springer-Verlag, New York.

Goriely:2010:GMV Goriely A, Moulton DE, Vandiver R. 2010. Elastic cavitation, tube hollowing, and differential growth in plants and biological tissues, Europhysics Letters 91, 18001.

Hauseux:2017:HHB Hauseux P, Hale JS, Bordas SPS. 2017. Accelerating Monte Carlo estimation with derivatives of high-level finite element models, Computer Methods in Applied Mechanics and Engineering 318, 917-936.

Hauseux:2018:HHCB Hauseux P, Hale JS, Cotin S, Bordas SPS. 2018. Quantifying the uncertainty in a hyperelastic soft tissue model with stochastic parameters, Applied Mathematical Modelling (doi: 10.1016/j.apm.2018.04.021)

Horgan:1989:HP Horgan CO, Pence TJ. 1989. Cavity formation at the center of a composite incompressible nonlinearly elastic sphere, Journal of Applied Mechanics 56, 302-308.

Horgan:1995:HP Horgan CO, Polignone DA. 1995. Cavitation in nonlinearly elastic solids: a review, Applied Mechanics Reviews 48, 471-485.

Hughes:2010:HH Hughes I, Hase TPA. 2010. Measurements and Their Uncertainties : A Practical Guide to Modern Error Analysis, Oxford University Press, Oxford.

James:1992:JS James RD, Spector SJ. 1991. The formation of filamentary voids in solids, Journal of Mechanics and Physics of Solids 39, 783-813.

Jaynes:1957a Jaynes ET. 1957. Information theory and statistical mechanics i, Physical Review 108, 171-190.

Jaynes:1957b Jaynes ET. 1957. Information theory and statistical mechanics ii, Physical Review 106, 620-630.

Jaynes:2003 Jaynes ET. 2003. Probability Theory: The Logic of Science, Cambridge University Press, Cambridge, UK.
	
Johnson:1994:JKB Johnson NL, Kotz S, Balakrishnan N. 1994. Continuous Univariate Distributions, Vol 1, 2nd edition, John Wiley & Sons, New York.

LopezPamies:2009 Lopez-Pamies O. 2009. Onset of cavitation in compressible, isotropic, hyperelastic solids, Journal of Elasticity 94, 115-145.

LopezPamies:2011:etal Lopez-Pamies O, Idiart MI, Nakamura T. 2011. Cavitation in elastomeric solids: I - A defect-growth theory, Journal of the Mechanics and Physics of Solids 59, 1464-1487.

Marzano:1983 Marzano M. 1983. An interpretation of Baker-Ericksen inequalities in uniaxial deformation and stress, Meccanica 18, 233-235.

McCoy:1973 McCoy JJ. 1973. A statistical theory for predicting response of materials that possess a disordered structure, Technical report ARPA 2181, AMCMS Code 5911.21.66022, Army Materials and Mechanics Research Center, Watertown, Massachusetts.

McMahon:2010:McMG McMahon J, Goriely A. 2010. Spontaneous cavitation in growing elastic membranes, Mathematics and Mechanics of Solids 15, 57-77.

Merodio:2006:MS Merodio J, Saccomandi G. 2006. Remarks on cavity formation in fiber-reinforced incompressible non-linearly elastic solids, European Journal of Mechanics A/Solids 25, 778–792.

Mihai:2011:MG Mihai LA, Goriely A. 2011. Positive or negative Poynting effect? The role of adscititious inequalities in hyperelastic materials, Proceedings of the Royal Society A 467, 3633-3646.

Mihai:2013:MG Mihai LA, Goriely A. 2013. Numerical simulation of shear and the Poynting effects by the finite element method: An application of the generalised empirical inequalities in non-linear elasticity, International Journal of Non-Linear Mechanics 49, 1-14.

Mihai:2017:MG Mihai LA, Goriely A. 2017. How to characterize a nonlinear elastic material? A review on nonlinear constitutive parameters in isotropic finite elasticity, Proceedings of the Royal Society A 473, 20170607 (doi: 10.1098/rspa.2017.0607).

Mihai:2018a:MWG Mihai LA, Woolley TE, Goriely A. 2018. Stochastic isotropic hyperelastic materials: constitutive calibration and model selection, Proceedings of the Royal Society A 474, 20170858.

Mihai:2018b:MWG Mihai LA, Woolley TE, Goriely A. 2018. Likely equilibria of the stochastic Rivlin cube, Philosophical Transactions of the Royal Society A, 20180068 (doi: 10.1098/rsta.2018.0068).

Mihai:2018:MDWG Mihai LA, Fitt D, Woolley TE, Goriely A. 2018. Likely equilibria of stochastic hyperelastic spherical shells and tubes, arxiv: 1808.02110.

Oden:2018 Oden JT. 2018. Adaptive multiscale predictive modelling, Acta Numerica 27, 353-450.

Ogden:1997 Ogden RW. 1997. Non-Linear Elastic Deformations, 2nd ed, Dover, New York.

Pence:2007:PT Pence TJ, Tsai SJ. 2007. Bulk cavitation and the possibility of localized deformation due to surface layer swelling, Journal of Elasticity 87, 161-185.

Polignone:1993a:PH Polignone DA, Horgan CO. 1993. Cavitation for incompressible anisotropic nonlinearly elastic spheres, Journal of Elasticity 33, 27-65.

Polignone:1993b:PH Polignone DA, Horgan CO. 1993. Effects of material unisotropy and inhomogeneity on cavitation for composite incompressible anisotropic nonlinearly elastic spheres, IInternational Journal of Solids and Structures 30, 3381-3416.

Poulain:2017:PLLPRC Poulain X, Lefèvre V, Lopez-Pamies O, Ravi-Chandar K. 2017. Damage in elastomers: nucleation and growth of cavities, micro-cracks, and macro-cracks, International Journal of Fracture 205, 1-21.

Poulain:2018:PLPRC Poulain X, Lopez-Pamies O, Ravi-Chandar K. 2018. Damage in elastomers: Healing of internally nucleated cavities and micro-cracks, Soft Matter 14, 4633-4640.

Rivlin:1953 Rivlin RS. 1953. The solution of problems in second order elasticity theory, Journal of Rational Mechanics and Analysis 2, 53-81.

Shannon:1948 Shannon CE. 1948. A mathematical theory of communication, Bell System Technical Journal 27, 379-423, 623-659.

Soni:2017:SG Soni J, Goodman R. 2017. A Mind at Play: How Claude Shannon Invented the Information Age, Simon & Schuster, New York.

Sivaloganathan:1992 Sivaloganathan, I. 1992. Cavitation, the incompressible limit, and material inhomogeneity, Quarterly of Applied Mathematics 49, 521-541.

Sivaloganathan:1999 Sivaloganathan J. 1999. On cavitation and degenerate cavitation under internal hydrostatic pressure, Proceedings of the Royal Society A 455, 3645-3664.

Soize:2000 Soize C. 2000. A nonparametric model of random uncertainties for reduced matrix models in structural dynamics, Probabilistic Engineering Mechanics 15, 277-294.

Soize:2001 Soize C. 2001. Maximum entropy approach for modeling random uncertainties in transient elastodynamics, Journal of the Acoustical Society of America 109, 1979-1996.

Staber:2015:SG Staber B, Guilleminot J. 2015. Stochastic modeling of a class of stored energy functions for incompressible hyperelastic materials with uncertainties, Comptes Rendus Mécanique 343, 503-514.

Staber:2016:SG Staber B, Guilleminot J. 2016. Stochastic modeling of the Ogden class of stored energy functions for hyperelastic materials: the compressible case, Journal of Applied Mathematics and Mechanics/Zeitschrift für Angewandte Mathematik und Mechanik 97, 273-295.

Staber:2017:SG Staber B, Guilleminot J. 2017. Stochastic hyperelastic constitutive laws and identification procedure for soft biological tissues with intrinsic variability, Journal of the Mechanical Behavior of Biomedical Materials 65, 743-752.

Staber:2018:SG Staber B, Guilleminot J. 2018. A random field model for anisotropic strain energy functions and its application for uncertainty quantification in vascular mechanics, Computer Methods in Applied Mechanics and Engineering 333, 94-113.

Steigmann:1992 Steigmann DJ. 1992. Cavitation in elastic membranes - an example, Journal of Elasticity 28, 277-287.

TruesdellNoll:2004 Truesdell C, Noll W. 2004. The Non-Linear Field Theories of Mechanics, 3rd ed, Springer-Verlag, New York.





