Source author record

Michael D. Gilchrist

Michael D. Gilchrist 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

12works
6topics
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

12 published item(s)

preprint2020arXiv

A combined experimental and numerical study of stab-penetration forces

The magnitude of force used in a stabbing incident can be difficult to quantify, although the estimate given by forensic pathologists is often seen as `critical' evidence in medico-legal situations. The main objective of this study is to develop a quantitative measure of the force associated with a knife stabbing biological tissue, using a combined experimental and numerical technique. A series of stab-penetration tests were performed to quantify the force required for a blade to penetrate skin at various speeds and using different `sharp' instruments. A computational model of blade penetration was developed using ABAQUS/EXPLICIT, a non-linear finite element analysis (FEA) commercial package. This model, which incorporated element deletion along with a suitable failure criterion, is capable of systematically quantifying the effect of the many variables affecting a stab event. This quantitative data could, in time, lead to the development of a predictive model that could help indicate the level of force used in a particular stabbing incident.

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

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.

preprint2020arXiv

Mechanical characterization of brain tissue in simple shear at dynamic strain rates

During severe impact conditions, brain tissue experiences a rapid and complex deformation, which can be seen as a mixture of compression, tension and shear. Moreover, diffuse axonal injury (DAI) occurs in animals and humans when both the strains and strain rates exceed 10% and 10/s, respectively. Knowing the mechanical properties of brain tissue in shear at these strains and strain rates is thus of particular importance, as they can be used in finite element simulations to predict the occurrence of brain injuries under different impact conditions. In this research, an experimental setup was developed to perform simple shear tests on porcine brain tissue at strain rates < 120/s. The maximum measured shear stress at strain rates of 30, 60, 90 and 120/s was 1.15 +/- 0.25 kPa, 1.34 +/- 0.19 kPa, 2.19 +/- 0.225 kPa and 2.52 +/- 0.27 kPa, (mean +/- SD), respectively, at the maximum amount of shear, K = 1. Good agreement of experimental, theoretical (Ogden and MooneyRivlin models) and numerical shear stresses was achieved (p = 0.7866 - 0.9935). Specimen thickness effects (2.0 - 10.0 mm thick specimens) were also analyzed numerically and we found that there is no significant difference (p = 0.9954) in the shear stress magnitudes, indicating a homogeneous deformation of the specimens during simple shear tests. Stress relaxation tests in simple shear were also conducted at different strain magnitudes (10% - 60% strain) with the average rise time of 14 ms. This allowed us to estimate elastic and viscoelastic parameters (mu = 4942.0 Pa and Prony parameters: g1 = 0.520, g2 = 0.3057, tau1 = 0.0264 s, tau2 = 0.011 s) that can be used in FE software to analyze the hyperviscoelastic behavior of brain tissue.

preprint2020arXiv

Mechanical characterization of brain tissue in tension at dynamic strain rates

Mechanical characterization of brain tissue at high loading velocities is crucial for modeling Traumatic Brain Injury (TBI). During severe impact conditions, brain tissue experiences compression, tension and shear. Limited experimental data is available for brain tissue in extension at dynamic strain rates. In this research, a High Rate Tension Device (HRTD) was developed to obtain dynamic properties of brain tissue in extension at strain rates of < 90/s. In vitro tensile tests were performed to obtain properties of brain tissue at strain rates of 30, 60 and 90/s up to 30% strain. The brain tissue showed a stiffer response with increasing strain rates, showing that hyperelastic models are not adequate. Specifically, the tensile engineering stress at 30% strain was 3.1 +/- 0.49 kPa, 4.3 +/- 0.86 kPa, 6.5 +/- 0.76 kPa (mean +/- SD) at strain rates of 30, 60 and 90/s, respectively. Force relaxation tests in tension were also conducted at different strain magnitudes (10-60% strain) with the average rise time of 24 ms, which were used to derive time dependent parameters. One-term Ogden, Fung and Gent models were used to obtain material parameters from the experimental data. Numerical simulations were performed using a one-term Ogden model to analyze hyperelastic behavior of brain tissue up to 30% strain. The material parameters obtained in this study will help to develop biofidelic human brain finite element models, which can subsequently be used to predict brain injuries under impact conditions and as a reconstruction and simulation tool for forensic investigations.

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

Characterising the Anisotropic Mechanical Properties of Excised Human Skin

The mechanical properties of skin are important for a number of applications including surgery, dermatology, impact biomechanics and forensic science. In this study we have investigated the influence of location and orientation on the deformation characteristics of 56 samples of excised human skin. Uniaxial tensile tests were carried out at a strain rate of 0.012s$^{-1}$ on skin from the back. Digital Image Correlation was used for 2D strain measurement and a histological examination of the dermis was also performed. The mean ultimate tensile strength (UTS) was 21.6$\pm$8.4MPa, the mean failure strain 54$\pm$17%, the mean initial slope 1.18$\pm$0.88MPa, the mean elastic modulus 83.3$\pm$34.9MPa and the mean strain energy was 3.6$\pm$1.6MJ/m$^3$. A multivariate analysis of variance has shown that these mechanical properties of skin are dependent upon the orientation of Langer lines (P$<$0.0001-P=0.046). The location of specimens on the back was also found to have a significant effect on the UTS (P =0.0002), the elastic modulus (P=0.001) and the strain energy (P=0.0052). The histological investigation concluded that there is a definite correlation between the orientation of Langer Lines and the preferred orientation of collagen fibres in the dermis (P$<$0.001). The data obtained in this study will provide essential information for those wishing to model the skin using a structural constitutive model.

preprint2013arXiv

Inhomogeneous Deformation of Brain Tissue During Tension Tests

Mechanical characterization of brain tissue has been investigated extensively by various research groups over the past fifty years. These properties are particularly important for modelling Traumatic Brain Injury (TBI). In this research, we present the design and calibration of a High Rate Tension Device (HRTD) capable of performing tests up to a maximum strain rate of 90/s. We use experimental and numerical methods to investigate the effects of inhomogeneous deformation of porcine brain tissue during tension at different specimen thicknesses (4.0-14.0 mm), by performing tension tests at a strain rate of 30/s. One-term Ogden material parameters (mu = 4395.0 Pa, alpha = -2.8) were derived by performing an inverse finite element analysis to model all experimental data. A similar procedure was adopted to determine Young's modulus (E= 11200 Pa) of the linear elastic regime. Based on this analysis, brain specimens of aspect ratio (diameter/thickness) S < 1.0 are required to minimise the effects of inhomogeneous deformation during tension tests.

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

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.

preprint2013arXiv

Third- and fourth-order elasticity of biological soft tissues

In the theory of weakly non-linear elasticity, Hamilton et al. [J. Acoust. Soc. Am. \textbf{116} (2004) 41] identified $W = μI_2 + (A/3)I_3 + D I_2^2$ as the fourth-order expansion of the strain-energy density for incompressible isotropic solids. Subsequently, much effort focused on theoretical and experimental developments linked to this expression in order to inform the modeling of gels and soft biological tissues. However, while many soft tissues can be treated as incompressible, they are not in general isotropic, and their anisotropy is associated with the presence of oriented collagen fiber bundles. Here the expansion of $W$ is carried up to fourth-order in the case where there exists one family of parallel fibers in the tissue. The results are then applied to acoustoelasticity, with a view to determining the second- and third-order nonlinear constants by employing small-amplitude transverse waves propagating in a deformed soft tissue.

preprint2012arXiv

Automated Estimation of Collagen Fibre Dispersion in the Dermis and its Contribution to the Anisotropic Behaviour of Skin

Collagen fibres play an important role in the mechanical behaviour of many soft tissues. Modelling of such tissues now often incorporates a collagen fibre distribution. However, the availability of accurate structural data has so far lagged behind the progress of anisotropic constitutive modelling. Here, an automated process is developed to identify the orientation of collagen fibres using inexpensive and relatively simple techniques. The method uses established histological techniques and an algorithm implemented in the MATLAB image processing toolbox. It takes an average of 15 s to evaluate one image, compared to several hours if assessed visually. The technique was applied to histological sections of human skin with different Langer line orientations and a definite correlation between the orientation of Langer lines and the preferred orientation of collagen fibres in the dermis was observed. The structural parameters of the Gasser-Ogden-Holzapfel (GOH) model were all successfully evaluated. It is expected that the results of this study will assist those wishing to model skin, and that the algorithm described will be of benefit to those who wish to evaluate the collagen dispersion of other soft tissues.