Source author record

Jerry G. Murphy

Jerry G. Murphy appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

9works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

9 published item(s)

preprint2020arXiv

Deficiencies in numerical models of anisotropic nonlinearly elastic materials

Incompressible nonlinearly hyperelastic materials are rarely simulated in Finite Element numerical experiments as being perfectly incompressible because of the numerical difficulties associated with globally satisfying this constraint. Most commercial Finite Element packages therefore assume that the material is slightly compressible. It is then further assumed that the corresponding strain-energy function can be decomposed additively into volumetric and deviatoric parts. We show that this decomposition is not physically realistic, especially for anisotropic materials, which are of particular interest for simulating the mechanical response of biological soft tissue. The most striking illustration of the shortcoming is that with this decomposition, an anisotropic cube under hydrostatic tension deforms into another cube instead of a hexahedron with non-parallel faces. Furthermore, commercial numerical codes require the specification of a `compressibility parameter' (or `penalty factor'), which arises naturally from the flawed additive decomposition of the strain-energy function. This parameter is often linked to a `bulk modulus', although this notion makes no sense for anisotropic solids; we show that it is essentially an arbitrary parameter and that infinitesimal changes to it result in significant changes in the predicted stress response. This is illustrated with numerical simulations for biaxial tension experiments of arteries, where the magnitude of the stress response is found to change by several orders of magnitude when infinitesimal changes in `Poisson's ratio' close to the perfect incompressibility limit of 1/2 are made.

preprint2020arXiv

Dominant negative Poynting effect for soft tissues

We identify three distinct shearing modes for simple shear deformations of transversely isotropic soft tissue which allow for both positive and negative Poynting effects (that is, they require compressive and tensile lateral normal stresses, respectively, in order to maintain simple shear). The positive Poynting effect is that usually found for isotropic rubber. Here, specialisation of the general results to three strain-energy functions which are quadratic in the anisotropic invariants, linear in the isotropic strain invariants and consistent with the linear theory suggests that there are two Poynting effects which can accompany the shearing of soft tissue: a dominant negative effect in one mode of shear and a relatively small positive effect in the other two modes. We propose that the relative inextensibility of the fibres relative to the matrix is the primary mechanism behind this large negative Poynting effect.

preprint2020arXiv

Extreme softness of brain matter in simple shear

We show that porcine brain matter can be modelled accurately as a very soft rubber-like material using the Mooney-Rivlin strain energy function, up to strains as high as 60\%. This result followed from simple shear experiments performed on small rectangular fresh samples ($2.5$ cm$^3$ and $1.1$ cm$^3$) at quasi-static strain rates. They revealed a linear shear stress--shear strain relationship ($R^2> 0.97 $), characteristic of Mooney-Rivlin materials at large strains. We found that porcine brain matter is about 30 times less resistant to shear forces than a silicone gel. We also verified experimentally that brain matter exhibits the positive Poynting effect of nonlinear elasticity, and numerically that the stress and strain fields remain mostly homogeneous throughout the thickness of the samples in simple shear.

preprint2013arXiv

Bimodular rubber buckles early in bending

A block of rubber eventually buckles under severe flexure, and several axial wrinkles appear on the inner curved face of the bent block. Experimental measurements reveal that the buckling occurs earlier ---at lower compressive strains--- than expected from theoretical predictions. This paper shows that if rubber is modeled as being bimodular, and specifically, as being stiffer in compression than in tension, then flexure bifurcation happens indeed at lower levels of compressive strain than predicted by previous investigations (these included taking into account finite size effects, compressibility effects, and strain-stiffening effects.) Here the effect of bimodularity is investigated within the theory of incremental buckling, and bifurcation equations, numerical methods, dispersion curves, and field variations are presented and discussed. It is also seen that Finite Element Analysis software seems to be unable to encompass in a realistic manner the phenomenon of bending instability for rubber blocks.

preprint2013arXiv

Differences in tension and compression in the nonlinearly elastic bending of beams

The classical flexure problem of non-linear incompressible elasticity is revisited for elastic materials whose mechanical response is different in tension and compression---the so-called bimodular materials. The flexure problem is chosen to investigate this response since the two regions, one of tension and one of compression, can be identified easily using simple intuition. Two distinct problems are considered: the first is where the stress is assumed continuous across the boundary of the two regions, which assumption has a sound physical basis. The second problem considered is more speculative: it is where discontinuities of stress are allowed. It is shown that such discontinuities are necessarily small for many applications, but might nonetheless provide an explanation for the damage incurred by repeated flexure. Some experimental evidence of the possibility of bimodularity in elastomers is also presented.

preprint2013arXiv

On deforming a sector of a circular cylindrical tube into an intact tube: existence, uniqueness, and stability

Within the context of finite deformation elasticity theory the problem of deforming an open sector of a thick-walled circular cylindrical tube into a complete circular cylindrical tube is analyzed. The analysis provides a means of estimating the radial and circumferential residual stress present in an intact tube, which is a problem of particular concern in dealing with the mechanical response of arteries. The initial sector is assumed to be unstressed and the stress distribution resulting from the closure of the sector is then calculated in the absence of loads on the cylindrical surfaces. Conditions on the form of the elastic strain-energy function required for existence and uniqueness of the deformed configuration are then examined. Finally, stability of the resulting finite deformation is analyzed using the theory of incremental deformations superimposed on the finite deformation, implemented in terms of the Stroh formulation. The main results are that convexity of the strain energy as a function of a certain deformation variable ensures existence and uniqueness of the residually-stressed intact tube, and that bifurcation can occur in the closing of thick, widely opened sectors, depending on the values of geometrical and physical parameters. The results are illustrated for particular choices of these parameters, based on data available in the biomechanics literature.

preprint2013arXiv

Onset of Non-Linearity in the Elastic Bending of Blocks

The classical flexure problem of non-linear incompressible elasticity is revisited assuming that the bending angle suffered by the block is specified instead of the usual applied moment. The general moment-bending angle relationship is then obtained and is shown to be dependent on only one non-dimensional parameter: the product of the aspect ratio of the block and the bending angle. A Maclaurin series expansion in this parameter is then found. The first-order term is proportional to $μ$, the shear modulus of linear elasticity; the second-order term is identically zero, because the moment is an odd function of the angle; and the third-order term is proportional to $μ(4β-1)$, where $β$ is the non-linear shear coefficient, involving third-order and fourth-order elasticity constants. It follows that bending experiments provide an alternative way of estimating this coefficient, and the results of one such experiment are presented. In passing, the coefficients of Rivlin's expansion in exact non-linear elasticity are connected to those of Landau in weakly (fourth-order) non-linear elasticity.

preprint2013arXiv

Simple shear is not so simple

For homogeneous, isotropic, nonlinearly elastic materials, the form of the homogeneous deformation consistent with the application of a Cauchy shear stress is derived here for both compressible and incompressible materials. It is shown that this deformation is not simple shear, in contrast to the situation in linear elasticity. Instead, it consists of a triaxial stretch superposed on a classical simple shear deformation, for which the amount of shear cannot be greater than 1. In other words, the faces of a cubic block cannot be slanted by an angle greater than 45deg by the application of a pure shear stress alone. The results are illustrated for those materials for which the strain energy function does not depend on the principal second invariant of strain. For the case of a block deformed into a parallelepiped, the tractions on the inclined faces necessary to maintain the derived deformation are calculated.

preprint2013arXiv

Slight compressibility and sensitivity to changes in Poisson's ratio

Finite Element simulations of rubbers and biological soft tissue usually assume that the material being deformed is slightly compressible. It is shown here that in shearing deformations the corresponding normal stress distribution can exhibit extreme sensitivity to changes in Poisson's ratio. These changes can even lead to a reversal of the usual Poynting effect. Therefore the usual practice of arbitrarily choosing a value of Poisson's ratio when numerically modelling rubbers and soft tissue will, almost certainly, lead to a significant difference between the simulated and actual normal stresses in a sheared block because of the difference between the assumed and actual value of Poisson's ratio. The worrying conclusion is that simulations based on arbitrarily specifying Poisson's ratio close to 1/2 cannot accurately predict the normal stress distribution even for the simplest of shearing deformations. It is shown analytically that this sensitivity is due to the small volume changes which inevitably accompany all deformations of rubber-like materials. To minimise these effects, great care should be exercised to accurately determine Poisson's ratio before simulations begin.