Peter Chadwick. 23 March 1931—12 August 2018

No abstract available

A modified Holzapfel-Ogden law for a residually stressed finite strain model of the human left ventricle in diastole

In this work, we introduce a modified Holzapfel-Ogden hyperelastic constitutive model for ventricular myocardium that accounts for residual stresses, and we investigate the effects of residual stresses in diastole using a magnetic resonance imaging–derived model of the human left ventricle (LV). We adopt an invariant-based constitutive modelling approach and treat the left ventricular myocardium as a non-homogeneous, fibre-reinforced, incompressible material. Because in vivo images provide the configuration of the LV in a loaded state even in diastole, an inverse analysis is used to determine the corresponding unloaded reference configuration. The residual stress in this unloaded state is estimated by two different methods. One is based on three-dimensional strain measurements in a local region of the canine LV, and the other uses the opening angle method for a cylindrical tube. We find that including residual stress in the model changes the stress distributions across the myocardium and that whereas both methods yield qualitatively similar changes, there are quantitative differences between the two approaches. Although the effects of residual stresses are relatively small in diastole, the model can be extended to explore the full impact of residual stress on LV mechanical behaviour for the whole cardiac cycle as more experimental data become available. In addition, although not considered here, residual stresses may also play a larger role in models that account for tissue growth and remodelling

Azimuthal shear of a transversely isotropic elastic solid

In this paper we study the problem of (plane strain) azimuthal shear of a circular cylindrical tube of incompressible transversely isotropic elastic material subject to finite deformation. The preferred direction associated with the transverse isotropy lies in the planes normal to the tube axis and is at an angle with the radial direction that depends only on the radius. For a general form of strain-energy function the considered deformation yields simple expressions for the azimuthal shear stress and the associated strong ellipticity condition in terms of the azimuthal shear strain. These apply for a sense of shear that is either “with” or “against” the preferred direction (anticlockwise and clockwise, respectively), so that material line elements locally in the preferred direction either extend or (at least initially) contract, respectively. For some specific strain-energy functions we then examine local loss of uniqueness of the shear stress–strain relationship and failure of ellipticity for the case of contraction and the dependence on the geometry of the preferred direction. In particular, for a reinforced neo-Hookean material, we obtain closed-form solutions that determine the domain of strong ellipticity in terms of the relationship between the shear strain and the angle (in general, a function of the radius) between the tangent to the preferred direction and the undeformed radial direction. It is shown, in particular, that as the magnitude of the applied shear stress increases then, after loss of ellipticity, there are two admissible values for the shear strain at certain radial locations. Absolutely stable deformations involve the lower magnitude value outside a certain radius and the higher magnitude value within this radius. The radius that separates the two values increases with increasing magnitude of the shear stress. The results are illustrated graphically for two specific forms of energy function

Two-phase piecewise homogeneous plane deformations of a fibre-reinforced neo-Hookean material with application to fibre kinking and splitting

Two-phase piecewise homogeneous plane deformations are examined in respect of a neo-Hookean matrix material reinforced with embedded aligned fibres characterized by a single stiffness parameter. The deformations are interpreted in terms of fibre kinking and fibre splitting. Previous work has shown that such a transversely isotropic material can lose ellipticity if the reinforcing stiffness is sufficiently large and the fibre direction is sufficiently compressed. In particular, it was shown that the associated failure modes are characterised by the emergence of weak surfaces of discontinuity that are normal to the fibre direction (the onset of fibre kinking) or parallel to the fibre direction (the onset of fibre splitting). Here, the analysis of strong surfaces of discontinuity, developing from weak ones, is studied. The considered model can give rise to piecewise smooth plane deformations separated by a plane stationary surface of discontinuity, interpreted as either kinking or splitting. Attention is restricted to (plane) deformations in which, on one side of the surface of discontinuity, the load axis is aligned with the fibre axis. Then the fibre stretch on this side of the discontinuity is a natural load parameter. The ellipticity status of the two-phase piecewise homogeneous plane deformations is shown to span all four possible ellipticity/non-ellipticity permutations. If both deformation states are elliptic, then a suitable intermediate deformation is shown to be non-elliptic. Moreover, it is shown that the mechanism is dissipative, and maximally dissipative quasi-static failure motion is examined in respect of both kinking and splitting. It follows that, firstly, surfaces of discontinuity perpendicular to the fibre direction, associated with fibre kinking, are nucleated followed by surfaces of discontinuity parallel to the fibre direction, associated with fibre splitting. With respect to kinking, such maximally dissipative kinks nucleate only in compression as weak surfaces of discontinuity, with the subsequent motion converting non-elliptic deformation to elliptic deformation

On fiber dispersion models: exclusion of compressed fibers and spurious model comparisons

Fiber dispersion in collagenous soft tissues has an important influence on the mechanical response, and the modeling of the collagen fiber architecture and its mechanics has developed significantly over the last few years. The purpose of this paper is twofold, first to develop a method for excluding compressed fibers within a dispersion for the generalized structure tensor (GST) model, which several times in the literature has been claimed not to be possible, and second to draw attention to several erroneous and misleading statements in the literature concerning the relative values of the GST and the angular integration (AI) models. For the GST model we develop a rather simple method involving a deformation dependent dispersion parameter that allows the mechanical influence of compressed fibers within a dispersion to be excluded. The theory is illustrated by application to simple extension and simple shear in order to highlight the effect of exclusion. By means of two examples we also show that the GST and the AI models have equivalent predictive power, contrary to some claims in the literature. We conclude that from the theoretical point of view neither of these two models is superior to the other. However, as is well known and as we now emphasize, the GST model has proved to be very successful in modeling the data from experiments on a wide range of tissues, and it is easier to analyze and simpler to implement than the AI approach, and the related computational effort is much lower

On Love-type waves in a finitely deformed magnetoelastic layered half-space

In this paper, the propagation of Love-type waves in a homogeneously and finitely deformed layered half-space of an incompressible non-conducting magnetoelastic material in the presence of an initial uniform magnetic field is analyzed. The equations and boundary conditions governing linearized incremental motions superimposed on an underlying deformation and magnetic field for a magnetoelastic material are summarized and then specialized to a form appropriate for the study of Love-type waves in a layered half-space. The wave propagation problem is then analyzed for different directions of the initial magnetic field for two different magnetoelastic energy functions, which are generalizations of the standard neo-Hookean and Mooney–Rivlin elasticity models. The resulting wave speed characteristics in general depend significantly on the initial magnetic field as well as on the initial finite deformation, and the results are illustrated graphically for different combinations of these parameters. In the absence of a layer, shear horizontal surface waves do not exist in a purely elastic material, but the presence of a magnetic field normal to the sagittal plane makes such waves possible, these being analogous to Bleustein–Gulyaev waves in piezoelectric materials. Such waves are discussed briefly at the end of the paper

Analysis and simulations for a phase‐field fracture model at finite strains based on modified invariants

Phase‐field models have already been proven to predict complex fracture patterns for brittle fracture at small strains. In this paper we discuss a model for phase‐field fracture at finite deformations in more detail. Among the identification of crack location and projection of crack growth the numerical stability is one of the main challenges in solid mechanics. Here we present a phase‐field model at finite strains, which takes into account the anisotropy of damage by applying an anisotropic split of the modified invariants of the right Cauchy‐Green strain tensor. We introduce a suitable weak notion of solution that also allows for a spatial and temporal discretization of the model. In this framework we study the existence of solutions and we show that the time‐discrete solutions converge in a weak sense to a solution of the time‐continuous formulation of the model. Numerical examples in two and three space dimensions illustrate the range of validity of the analytical results

Anisotropic behaviour of human gallbladder walls

Inverse estimation of biomechanical parameters of soft tissues from non-invasive measurements has clinical significance in patient-specific modelling and disease diagnosis. In this paper, we propose a fully nonlinear approach to estimate the mechanical properties of the human gallbladder wall muscles from in vivo ultrasound images. The iteration method consists of a forward approach, in which the constitutive equation is based on a modified Hozapfel–Gasser–Ogden law initially developed for arteries. Five constitutive parameters describing the two orthogonal families of fibres and the matrix material are determined by comparing the computed displacements with medical images. The optimisation process is carried out using the MATLAB toolbox, a Python code, and the ABAQUS solver. The proposed method is validated with published artery data and subsequently applied to ten human gallbladder samples. Results show that the human gallbladder wall is anisotropic during the passive refilling phase, and that the peak stress is 1.6 times greater than that calculated using linear mechanics. This discrepancy arises because the wall thickness reduces by 1.6 times during the deformation, which is not predicted by conventional linear elasticity. If the change of wall thickness is accounted for, then the linear model can used to predict the gallbladder stress and its correlation with pain. This work provides further understanding of the nonlinear characteristics of human gallbladder

On derivation of Euler-Lagrange Equations for incompressible energy-minimizers

We prove that any distribution $q$ satisfying the equation $\nabla q=\div{\bf f}$ for some tensor ${\bf f}=(f^i_j), f^i_j\in h^r(U)$ ($1\leq r<\infty$) -the {\it local Hardy space}, $q$ is in $h^r$, and is locally represented by the sum of singular integrals of $f^i_j$ with Calder\'on-Zygmund kernel. As a consequence, we prove the existence and the local representation of the hydrostatic pressure $p$ (modulo constant) associated with incompressible elastic energy-minimizing deformation ${\bf u}$ satisfying $|\nabla {\bf u}|^2, |{\rm cof}\nabla{\bf u}|^2\in h^1$. We also derive the system of Euler-Lagrange equations for incompressible local minimizers ${\bf u}$ that are in the space $K^{1,3}_{\rm loc}$; partially resolving a long standing problem. For H\"older continuous pressure $p$, we obtain partial regularity of area-preserving minimizers.Comment: 23 page