A mathematical model for the auxetic response of liquid crystal elastomers.

We develop a mathematical model that builds on the surprising nonlinear mechanical response observed in recent experiments on nematic liquid crystal elastomers. Namely, under uniaxial tensile loads, the material, rather than thinning in the perpendicular directions, becomes thicker in one direction for a sufficiently large strain, while its volume remains unchanged. Motivated by this unusual large-strain auxetic behaviour, we model the material using an Ogden-type strain-energy function and calibrate its parameters to available datasets. We show that Ogden strain-energy functions are particularly suitable for modelling nematic elastomers because of their mathematical simplicity and their clear formulation in terms of the principal stretches, which have a direct kinematic interpretation. This article is part of the theme issue 'The Ogden model of rubber mechanics: Fifty years of impact on nonlinear elasticity'.

Nematic liquid crystalline elastomers are aeolotropic materials.

Continuum models describing ideal nematic solids are widely used in theoretical studies of liquid crystal elastomers. However, experiments on nematic elastomers show a type of anisotropic response that is not predicted by the ideal models. Therefore, their description requires an additional term coupling elastic and nematic responses, to account for aeolotropic effects. In order to better understand the observed elastic response of liquid crystal elastomers, we analyse theoretically and computationally different stretch and shear deformations. We then compare the elastic moduli in the infinitesimal elastic strain limit obtained from the molecular dynamics simulations with the ones derived theoretically, and show that they are better explained by including nematic order effects within the continuum framework.

Likely equilibria of the stochastic Rivlin cube.

The problem of the Rivlin cube is to determine the stability of all homogeneous equilibria of an isotropic incompressible hyperelastic body under equitriaxial dead loads. Here, we consider the stochastic version of this problem where the elastic parameters are random variables following standard probability laws. Uncertainties in these parameters may arise, for example, from inherent data variation between different batches of homogeneous samples, or from different experimental tests. As for the deterministic elastic problem, we consider the following questions: what are the likely equilibria and how does their stability depend on the material constitutive law? In addition, for the stochastic model, the problem is to derive the probability distribution of deformations, given the variability of the parameters. This article is part of the theme issue 'Rivlin's legacy in continuum mechanics and applied mathematics'.

Nonlinear scaling effects in the stiffness of soft cellular structures.

For cellular structures with uniform geometry, cell size and distribution, made from a neo-Hookean material, we demonstrate experimentally that large stretching causes nonlinear scaling effects governed by the microstructural architecture and the large strains at the cell level, which are not predicted by the linear elastic theory. For this purpose, three honeycomb-like structures with uniform square cells in stacked distribution were designed, where the number of cells varied, while the material volume and the ratio between the thickness and the length of the cell walls were fixed. These structures were manufactured from silicone rubber and tested under large uniaxial tension in a bespoke test fixture. Optical strain measurements were used to assess the deformation by capturing both the global displacements of the structure and the local deformations in the form of a strain map. The experimental results showed that, under sufficiently large strains, there was an increase in the stiffness of the structure when the same volume of material was arranged as many small cells compared to when it was organized as fewer larger cells. Finite element simulations confirmed our experimental findings. This study sheds light upon the nonlinear elastic responses of cellular structures in large-strain deformations, which cannot be captured within the linear elasticity framework.

Stochastic isotropic hyperelastic materials: constitutive calibration and model selection.

Biological and synthetic materials often exhibit intrinsic variability in their elastic responses under large strains, owing to microstructural inhomogeneity or when elastic data are extracted from viscoelastic mechanical tests. For these materials, although hyperelastic models calibrated to mean data are useful, stochastic representations accounting also for data dispersion carry extra information about the variability of material properties found in practical applications. We combine finite elasticity and information theories to construct homogeneous isotropic hyperelastic models with random field parameters calibrated to discrete mean values and standard deviations of either the stress-strain function or the nonlinear shear modulus, which is a function of the deformation, estimated from experimental tests. These quantities can take on different values, corresponding to possible outcomes of the experiments. As multiple models can be derived that adequately represent the observed phenomena, we apply Occam's razor by providing an explicit criterion for model selection based on Bayesian statistics. We then employ this criterion to select a model among competing models calibrated to experimental data for rubber and brain tissue under single or multiaxial loads.

How to characterize a nonlinear elastic material? A review on nonlinear constitutive parameters in isotropic finite elasticity.

The mechanical response of a homogeneous isotropic linearly elastic material can be fully characterized by two physical constants, the Young's modulus and the Poisson's ratio, which can be derived by simple tensile experiments. Any other linear elastic parameter can be obtained from these two constants. By contrast, the physical responses of nonlinear elastic materials are generally described by parameters which are scalar functions of the deformation, and their particular choice is not always clear. Here, we review in a unified theoretical framework several nonlinear constitutive parameters, including the stretch modulus, the shear modulus and the Poisson function, that are defined for homogeneous isotropic hyperelastic materials and are measurable under axial or shear experimental tests. These parameters represent changes in the material properties as the deformation progresses, and can be identified with their linear equivalent when the deformations are small. Universal relations between certain of these parameters are further established, and then used to quantify nonlinear elastic responses in several hyperelastic models for rubber, soft tissue and foams. The general parameters identified here can also be viewed as a flexible basis for coupling elastic responses in multi-scale processes, where an open challenge is the transfer of meaningful information between scales.

Microstructure-based hyperelastic models for closed-cell solids.

For cellular bodies involving large elastic deformations, mesoscopic continuum models that take into account the interplay between the geometry and the microstructural responses of the constituents are developed, analysed and compared with finite-element simulations of cellular structures with different architecture. For these models, constitutive restrictions for the physical plausibility of the material responses are established, and global descriptors such as nonlinear elastic and shear moduli and Poisson's ratio are obtained from the material characteristics of the constituents. Numerical results show that these models capture well the mechanical responses of finite-element simulations for three-dimensional periodic structures of neo-Hookean material with closed cells under large tension. In particular, the mesoscopic models predict the macroscopic stiffening of the structure when the stiffness of the cell-core increases.

The optimal density of cellular solids in axial tension.

For cellular bodies with uniform cell size, wall thickness, and shape, an important question is whether the same volume of material has the same effect when arranged as many small cells or as fewer large cells. To answer this question, for finite element models of periodic structures of Mooney-type material with different structural geometry and subject to large strain deformations, we identify a nonlinear elastic modulus as the ratio between the mean effective stress and the mean effective strain in the solid cell walls, and show that this modulus increases when the thickness of the walls increases, as well as when the number of cells increases while the volume of solid material remains fixed. Since, under the specified conditions, this nonlinear elastic modulus increases also as the corresponding mean stress increases, either the mean modulus or the mean stress can be employed as indicator when the optimum wall thickness or number of cells is sought.

A comparison of hyperelastic constitutive models applicable to brain and fat tissues.

In some soft biological structures such as brain and fat tissues, strong experimental evidence suggests that the shear modulus increases significantly under increasing compressive strain, but not under tensile strain, whereas the apparent Young's elastic modulus increases or remains almost constant when compressive strain increases. These tissues also exhibit a predominantly isotropic, incompressible behaviour. Our aim is to capture these seemingly contradictory mechanical behaviours, both qualitatively and quantitatively, within the framework of finite elasticity, by modelling a soft tissue as a homogeneous, isotropic, incompressible, hyperelastic material and comparing our results with available experimental data. Our analysis reveals that the Fung and Gent models, which are typically used to model soft tissues, are inadequate for the modelling of brain or fat under combined stretch and shear, and so are the classical neo-Hookean and Mooney-Rivlin models used for elastomers. However, a subclass of Ogden hyperelastic models are found to be in excellent agreement with the experiments. Our findings provide explicit models suitable for integration in large-scale finite-element computations.