William J. Parnell

Nonlinear pre-stress for cloaking from antiplane elastic waves

A theory is presented showing that cloaking of objects from antiplane elastic waves can be achieved by employing nonlinear elastic pre-stress in a neo-Hookean elastomeric material. This approach would appear to eliminate the requirement of metamaterials with inhomogeneous anisotropic shear moduli and density. Waves in the pre-stressed medium are bent around the cloaked (cavity) region by inducing inhomogeneous stress fields via pre-stress. The equation governing antiplane waves in the pre-stressed medium is equivalent to the antiplane equation in an unstressed medium with inhomogeneous and anisotropic shear modulus and isotropic scalar mass density. Note however that these properties are induced naturally by the pre-stress. As the magnitude of pre-stress can be altered at will, this enables objects of varying size and shape to be cloaked by placing them inside the fluid-filled deformed cavity region.

Employing pre-stress to generate finite cloaks for antiplane elastic waves

It is shown that nonlinear elastic pre-stress of neo-Hookean hyperelastic materials can be used as a mechanism to generate finite cloaks and thus render objects near-invisible to incoming antiplane elastic waves. This approach appears to negate the requirement for special cloaking metamaterials with inhomogeneous and anisotropic material properties in this case. These properties are induced naturally by virtue of the pre-stress. This appears to provide a mechanism for broadband cloaking since dispersive effects due to metamaterial microstructure will not arise.

Hyperelastic cloaking theory: transformation elasticity with pre-stressed solids

Transformation elasticity, by analogy with transformation acoustics and optics, converts material domains without altering wave properties, thereby enabling cloaking and related effects. By noting the similarity between transformation elasticity and the theory of incremental motion superimposed on finite pre-strain, it is shown that the constitutive parameters of transformation elasticity correspond to the density and moduli of small-on-large theory. The formal equivalence indicates that transformation elasticity can be achieved by selecting a particular finite (hyperelastic) strain energy function, which for isotropic elasticity is semilinear strain energy. The associated elastic transformation is restricted by the requirement of statically equilibrated pre-stress. This constraint can be cast as trF = constant, where F is the deformation gradient, subject to symmetry constraints, and its consequences are explored both analytically and through numerical examples of cloaking of anti-plane and in-plane wave motion.

A two-parameter model of the effective elastic tensor for cortical bone.

Multiscale models of cortical bone elasticity require a large number of parameters to describe the organization and composition of the tissue. We hypothesize that the macro-scale anisotropic elastic properties of different bones can be modeled retaining only two variable parameters, and setting the others to universal values identical for all bones. Cortical bone is regarded as a two-phase composite material: a dense mineralized matrix (ultrastructure) and a soft phase (pores). The ultrastructure is assumed to be a homogeneous and transversely isotropic tissue whose elastic properties in different directions are mutually dependent and can be scaled with a single parameter driving the overall rigidity. This parameter is taken to be the volume fraction of mineral f(ha). The pore network is modeled as an ensemble of water-filled cylinders and described only by the porosity p. The effective macroscopic elasticity tensor C(ij)(f(ha),p) is calculated with a multiscale micromechanics approach starting from existing models. The modeled stiffness coefficients compare favorably to four literature datasets which were chosen because they provide the full stiffness tensors of groups of human samples. Since the physical counterparts of f(ha) and p were unknown for the datasets, their values which allow the best fit of experimental tensors by the modeled ones were determined by optimization. Optimum values of f(ha) and p are found to be unique and realistic. These results suggest that a two-parameter model may be sufficient to model the elasticity of different samples of human femora and tibiae. Such a model would in particular be useful in large-scale parametric studies of bone mechanical response.

The influence of mesoscale porosity on cortical bone anisotropy. Investigations via asymptotic homogenization

Recently, the mesoscale of cortical bone has been given particular attention in association with novel experimental techniques such as nanoindentation, micro-computed X-ray tomography and quantitative scanning acoustic microscopy (SAM). A need has emerged for reliable mathematical models to interpret the related microscopic and mesoscopic data in terms of effective elastic properties. In this work, a new model of cortical bone elasticity is developed and used to assess the influence of mesoscale porosity on the induced anisotropy of the material. Only the largest pores (Haversian canals and resorption cavities), characteristic of the mesoscale, are considered. The input parameters of the model are derived from typical mesoscale experimental data (e.g. SAM data). We use the method of asymptotic homogenization to determine the local effective elastic properties by modelling the propagation of low-frequency elastic waves through an idealized material that models the local mesostructure. We use a novel solution of the cell problem developed by Parnell & Abrahams. This solution is stable for the physiological range of variation of mesoscopic porosity and elasticity found in bone. Results are computed efficiently (in seconds) and the solutions can be implemented easily by other workers. Parametric studies are performed in order to assess the influence of mesoscopic porosity, the assumptions regarding the material inside the mesoscale pores (drained or undrained bone) and the shape of pores. Results are shown to be in good qualitative agreement with existing schemes and we describe the potential of the scheme for future use in modelling more complex microstructures for cortical bone. In particular, the scheme is shown to be a useful tool with which to predict the qualitative changes in anisotropy due to variations in the structure at the mesoscale.

On nonlinear viscoelastic deformations: a reappraisal of Fung's quasi-linear viscoelastic model

This paper offers a reappraisal of Fungs model for quasi-linear viscoelasticity. It is shown that a number of negative features exhibited in other works, commonly attributed to the Fung approach, are merely a consequence of the way it has been applied. The approach outlined herein is shown to yield improved behaviour and offers a straightforward scheme for solving a wide range of models. Results from the new model are contrasted with those in the literature for the case of uniaxial elongation of a bar: for an imposed stretch of an incompressible bar and for an imposed load. In the latter case, a numerical solution to a Volterra integral equation is required to obtain the results. This is achieved by a high-order discretization scheme. Finally, the stretch of a compressible viscoelastic bar is determined for two distinct materials: Horgan–Murphy and Gent.

Estimating the dynamic effective mass density of random composites

The effective mass density of an inhomogeneous medium is discussed. Random configurations of circular cylindrical scatterers are considered, in various physical contexts: fluid cylinders in another fluid, elastic cylinders in a fluid or in another solid, and movable rigid cylinders in a fluid. In each case, time-harmonic waves are scattered, and an expression for the effective wavenumber due to Linton and Martin [J. Acoust. Soc. Am. 117, 3413-3423 (2005)] is used to derive the effective density in the low frequency limit, correct to second order in the area fraction occupied by the scatterers. Expressions are recovered that agree with either the Ament formula or the effective static mass density, depending upon the physical context.

Analytical methods to determine the effective mesoscopic and macroscopic elastic properties of cortical bone

We compare theoretical predictions of the effective elastic moduli of cortical bone at both the meso- and macroscales. We consider the efficacy of three alternative approaches: the method of asymptotic homogenization, the Mori–Tanaka scheme and the Hashin–Rosen bounds. The methods concur for specific engineering moduli such as the axial Young’s modulus but can vary for others. In a past study, the effect of porosity alone on mesoscopic properties of cortical bone was considered, taking the matrix to be isotropic. Here, we consider the additional influence of the transverse isotropy of the matrix. We make the point that micromechanical approaches can be used in two alternative ways to predict either the macroscopic (size of cortical bone sample) or mesoscopic (in between micro- and macroscales) effective moduli, depending upon the choice of representative volume element size. It is widely accepted that the mesoscale behaviour is an important aspect of the mechanical behaviour of bone but models incorporating its effect have started to appear only relatively recently. Before this only macroscopic behaviour was addressed. Comparisons are drawn with experimental data and simulations from the literature for macroscale predictions with particularly good agreement in the case of dry bone. Finally, we show how predictions of the effective mesoscopic elastic moduli can be made which retain dependence on the well-known porosity gradient across the thickness of cortical bone.

Source amplitudes for active exterior cloaking

The active cloak comprises a discrete set of multipole sources that destructively interfere with an incident time harmonic scalar wave to produce zero total field over a finite spatial region. For a given number of sources and their positions in two dimensions it is shown that the multipole amplitudes can be expressed as infinite sums of the coefficients of the incident wave decomposed into regular Bessel functions. The field generated by the active sources vanishes in the infinite region exterior to a set of circles defined by the relative positions of the sources. The results provide a direct solution to the inverse problem of determining the source amplitudes. They also define a broad class of nonradiating discrete sources. (Some figures may appear in colour only in the online journal)

Active elastodynamic cloaking

An active elastodynamic cloak destructively interferes with an incident time harmonic in-plane (coupled compressional/shear) elastic wave to produce zero total elastic field over a finite spatial region. A method is described which explicitly predicts the source amplitudes of the active field. For a given number of sources and their positions in two dimensions it is shown that the multipole amplitudes can be expressed as infinite sums of the coefficients of the incident wave decomposed into regular Bessel functions. Importantly, the active field generated by the sources vanishes in the far-field. In practice the infinite summations are clearly required to be truncated and the accuracy of cloaking is studied when the truncation parameter is modified.