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 postbifurcation 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.We recall that a homogeneous hyperelastic model is described by a strain-energy function W (F) that depends on the deformation gradient tensor, F, with respect to a fixed reference configuration, and is characterised by a set of deterministic model parameters [14,39,57]. 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 [35,[52][53][54]. 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

In large deformations, internally pressurised elastic spherical shells and tubes may undergo a limit-point, or inflation, instability manifested by a rapid transition in which their radii suddenly increase. The possible existence of such an instability depends on the material constitutive model. Here, we revisit this problem in the context of stochastic incompressible hyperelastic materials, and ask the question: what is the probability distribution of stable radially symmetric inflation, such that the internal pressure always increases as the radial stretch increases? For the classic elastic problem, involving isotropic incompressible materials, there is a critical parameter value that strictly separates the cases where inflation instability can occur or not. By contrast, for the stochastic problem, we show that the inherent variability of the probabilistic parameters implies that there is always competition between the two cases. To illustrate this, we draw on published experimental data for rubber, and derive the probability distribution of the corresponding random shear modulus to predict the inflation responses for a spherical shell and a cylindrical tube made of a material characterised by this parameter.

Stochastic homogeneous hyperelastic solids are characterised by strain-energy densities where the parameters are random variables defined by probability density functions. These models allow for the propagation of uncertainties from input data to output quantities of interest. To investigate the effect of probabilistic parameters on predicted mechanical responses, we study radial oscillations of cylindrical and spherical shells of stochastic incompressible isotropic hyperelastic material, formulated as quasi-equilibrated motions where the system is in equilibrium at every time instant. Additionally, we study finite shear oscillations of a cuboid, which are not quasi-equilibrated. We find that, for hyperelastic bodies of stochastic neo-Hookean or Mooney-Rivlin material, the amplitude and period of the oscillations follow probability distributions that can be characterised. Further, for cylindrical tubes and spherical shells, when an impulse surface traction is applied, there is a parameter interval where the oscillatory and non-oscillatory motions compete, in the sense that both have a chance to occur with a given probability. We refer to the dynamic evolution of these elastic systems, which exhibit inherent uncertainties due to the material properties, as "likely oscillatory motions".Key words: stochastic hyperelastic models, dynamic finite strain deformation, quasi-equilibrated motion, finite amplitude oscillations, incompressibility, applied probability."Denominetur motus talis, qualis omni momento temporis t praebet configurationem capacem aequilibrii corporis iisdem viribus massalibus sollicitati, 'motus quasi aequilibratus'. Generatim motus quasi aequilibratus non congruet legibus dynamicis et proinde motus verus corporis fieri non potest, manentibus iisdem viribus masalibus." -C. Truesdell (1962) [103]

Motivated by the need to quantify uncertainties in the mechanical behaviour of solid materials, we perform simple uniaxial tensile tests on a manufactured rubber-like material that provide critical information regarding the variability in the constitutive responses between different specimens. Based on the experimental data, we construct stochastic homogeneous hyperelastic models where the parameters are described by spatially independent probability density functions at a macroscopic level. As more than one parametrised model is capable of capturing the observed material behaviour, we apply Bayes' theorem to select the model that is most likely to reproduce the data. Our analysis is fully tractable mathematically and builds directly on knowledge from deterministic finite elasticity. The proposed stochastic calibration and Bayesian model selection are generally applicable to more complex tests and materials. Keywords Stochastic elasticity • Finite strain analysis • Hyperelastic material • Bayes' factor • Experiments • Probabilities "This task is made more difficult than it otherwise would be by the fact that some of the test-pieces used have to be moulded individually, and it is difficult to make two rubber specimens having identical properties even if nominally identical procedures are followed in preparing them."-R.S. Rivlin and D.W. Saunders [55]