Artur L. Gower

A robust anisotropic hyperelastic formulation for the modelling of soft tissue

The Holzapfel-Gasser-Ogden (HGO) model for anisotropic hyperelastic behaviour of collagen fibre reinforced materials was initially developed to describe the elastic properties of arterial tissue, but is now used extensively for modelling a variety of soft biological tissues. Such materials can be regarded as incompressible, and when the incompressibility condition is adopted the strain energy Ψ of the HGO model is a function of one isotropic and two anisotropic deformation invariants. A compressible form (HGO-C model) is widely used in finite element simulations whereby the isotropic part of Ψ is decoupled into volumetric and isochoric parts and the anisotropic part of Ψ is expressed in terms of isochoric invariants. Here, by using three simple deformations (pure dilatation, pure shear and uniaxial stretch), we demonstrate that the compressible HGO-C formulation does not correctly model compressible anisotropic material behaviour, because the anisotropic component of the model is insensitive to volumetric deformation due to the use of isochoric anisotropic invariants. In order to correctly model compressible anisotropic behaviour we present a modified anisotropic (MA) model, whereby the full anisotropic invariants are used, so that a volumetric anisotropic contribution is represented. The MA model correctly predicts an anisotropic response to hydrostatic tensile loading, whereby a sphere deforms into an ellipsoid. It also computes the correct anisotropic stress state for pure shear and uniaxial deformations. To look at more practical applications, we developed a finite element user-defined material subroutine for the simulation of stent deployment in a slightly compressible artery. Significantly higher stress triaxiality and arterial compliance are computed when the full anisotropic invariants are used (MA model) instead of the isochoric form (HGO-C model).

Initial stress symmetry and its applications in elasticity

An initial stress within a solid can arise to support external loads or from processes such as thermal expansion in inert matter or growth and remodelling in living materials. For this reason, it is useful to develop a mechanical framework of initially stressed solids irrespective of how this stress formed. An ideal way to do this is to write the free energy density Ψ in terms of initial stress τ and the elastic deformation gradient F, so we write Ψ=Ψ(F,τ). In this paper, we present a new constitutive condition for initially stressed materials, which we call the initial stress symmetry (ISS). We focus on two consequences of this condition. First, we examine how ISS restricts the possible choices of free energy densities Ψ=Ψ(F,τ) and present two examples of Ψ that satisfy the ISS. Second, we show that the initial stress can be derived from the Cauchy stress and the elastic deformation gradient. To illustrate, we take an example from biomechanics and calculate the optimal Cauchy stress within an artery subjected to internal pressure. We then use ISS to derive the optimal target residual stress for the material to achieve after remodelling, which links nicely with the notion of homeostasis.

Shear instability in skin tissue

Partial funding by the European Community grant ERG-256605, FP7 program, and by the Hardiman Scholarship programme at the National University of Ireland Galway to the first and third authors, respectively.

On residual stresses and homeostasis: an elastic theory of functional adaptation in living matter

Living matter can functionally adapt to external physical factors by developing internal tensions, easily revealed by cutting experiments. Nonetheless, residual stresses intrinsically have a complex spatial distribution, and destructive techniques cannot be used to identify a natural stress-free configuration. This work proposes a novel elastic theory of pre-stressed materials. Imposing physical compatibility and symmetry arguments, we define a new class of free energies explicitly depending on the internal stresses. This theory is finally applied to the study of arterial remodelling, proving its potential for the non-destructive determination of the residual tensions within biological materials.

Counter-intuitive results in acousto-elasticity

Abstract We present examples of body wave and surface wave propagation in deformed solids where the slowest and the fastest waves do not travel along the directions of least and greatest stretch, respectively. These results run counter to commonly accepted theory, practice, and implementation of the principles of acousto-elasticity in initially isotropic solids. For instance, we find that in nickel and steel the fastest waves are along the direction of greatest compression, not greatest extension (and vice-versa for the slowest waves), as soon as those solids are deformed. Further, we find that when some materials are subject to a small-but-finite deformation, other extrema of wave speeds appear in non-principal directions. Examples include nickel, steel, polystyrene, and a certain hydrogel. The existence of these “oblique”, non-principal extremal waves complicates the protocols for the non-destructive determination of the directions of extreme strains.

Higher-order reverse automatic differentiation with emphasis on the third-order

It is commonly assumed that calculating third order information is too expensive for most applications. But we show that the directional derivative of the Hessian (

Comparison of finite element and analytical modeling of scattering of an acoustic wave by particles in a fluid

D^3 f(x)\cdot d

Connecting the material parameters of soft fibre-reinforced solids with the formation of surface wrinkles

D3f(x)·d) can be calculated at a cost proportional to that of a state-of-the-art method for calculating the Hessian matrix. We do this by first presenting a simple procedure for designing high order reverse methods and applying it to deduce several methods including a reverse method that calculates