The mechanical response of the ovine lumbar anulus fibrosus to uniaxial, biaxial and shear loads

Analytical and computational models of the intervertebral disc (IVD) are commonly employed to enhance understanding of the biomechanics of the human spine and spinal motion segments. The accuracy of these models in predicting physiological behaviour of the spine is intrinsically reliant on the accuracy of the material constitutive representations employed to represent the spinal tissues. There is a paucity of detailed mechanical data describing the material response of the reinforced­ground matrix in the anulus fibrosus of the IVD. In the present study, the ‘reinforced­ground matrix’ was defined as the matrix with the collagen fibres embedded but not actively bearing axial load, thus incorporating the contribution of the fibre-fibre and fibre-matrix interactions. To determine mechanical parameters for the anulus ground matrix, mechanical tests were carried out on specimens of ovine anulus, under unconfined uniaxial compression, simple shear and biaxial compression. Test specimens of ovine anulus fibrosus were obtained with an adjacent layer of vertebral bone/cartilage on the superior and inferior specimen surface. Specimen geometry was such that there were no continuous collagen fibres coupling the two endplates. Samples were subdivided according to disc region - anterior, lateral and posterior - to determine the regional inhomogeneity in the anulus mechanical response. Specimens were loaded at a strain rate sufficient to avoid fluid outflow from the tissue and typical stress-strain responses under the initial load application and under repeated loading were determined for each of the three loading types. The response of the anulus tissue to the initial and repeated load cycles was significantly different for all load types, except biaxial compression in the anterior anulus. Since the maximum applied strain exceeded the damage strain for the tissue, experimental results for repeated loading reflected the mechanical ability of the tissue to carry load, subsequent to the initiation of damage. To our knowledge, this is the first study to provide experimental data describing the response of the ‘reinforced­ground matrix’ to biaxial compression. Additionally, it is novel in defining a study objective to determine the regionally inhomogeneous response of the ‘reinforced­ground matrix’ under an extensive range of loading conditions suitable for mechanical characterisation of the tissue. The results presented facilitate the development of more detailed and comprehensive constitutive descriptions for the large strain nonlinear elastic or hyperelastic response of the anulus ground matrix.

Fit optimisation of a distal medial tibia plate

An iterative method for the fit optimisation of a pre-contoured fracture fixation plate for a given bone data set is presented. Both plate shape optimisation and plate fit quantification are conducted in a virtual environment utilising computer graphical methods and 3D bone and plate models. Two optimised shapes of the undersurface of an existing distal medial tibia plate were generated based on a dataset of 45 3D bone models reconstructed from computed tomography image data of Japanese tibiae. The existing plate shape achieved an anatomical fit on 13% of tibiae from the dataset. Modified plate 1 achieved an anatomical fit for 42% and modified plate 2 a fit for 67% of the bones. If either modified plate 1 or plate 2 is used, then the anatomical fit can be increased to 82% for the same dataset. Issues pertaining to any further improvement in plate fit/shape are discussed.

Engineering tubular bone constructs

Cell-sheet techniques have been proven effective in various soft tissue engineering applications. In this experiment, we investigated the feasibility of bone tissue engineering using a hybrid of mesenchymal stem cell (MSC) sheets and PLGA meshes. Porcine MSCs were cultured to a thin layer of cell sheets via osteogenic induction. Tube-like long bones were constructed by wrapping the cell sheet on to PLGA meshes resulting in constructs which could be cultured in spinner flasks, prior to implantation in nude rats. Our results showed that the sheets were composed of viable cells and dense matrix with a thickness of about 80–120 mm, mineral deposition was also observed in the sheet. In vitro cultures demonstrated calcified cartilage-like tissue formation and most PLGA meshes were absorbed during the 8-week culture period. In vivo experiments revealed that dense mineralized tissue was formed in subcutaneous sites and the 8- week plants shared similar micro-CT characteristics with native bone. The neo tissue demonstrated histological markers for both bone and cartilage, indicating that the bone formation pathway in constructs was akin to endochondral ossification, with the residues of PLGA having an effect on the neo tissue organization and formation. These results indicate that cell-sheet approaches in combination with custom-shaped scaffolds have potential in producing bone tissue.

Improved image contrast of the bone-muscle interface with 3T MRI compared to 1.5T MRI [Abstract]

Virtual 3D models of long bones are increasingly being used for implant design and research applications. The current gold standard for the acquisition of such data is Computed Tomography (CT) scanning. Due to radiation exposure, CT is generally limited to the imaging of clinical cases and cadaver specimens. Magnetic Resonance Imaging (MRI) does not involve ionising radiation and therefore can be used to image selected healthy human volunteers for research purposes. The feasibility of MRI as alternative to CT for the acquisition of morphological bone data of the lower extremity has been demonstrated in recent studies [1, 2]. Some of the current limitations of MRI are long scanning times and difficulties with image segmentation in certain anatomical regions due to poor contrast between bone and surrounding muscle tissues. Higher field strength scanners promise to offer faster imaging times or better image quality. In this study image quality at 1.5T is quantitatively compared to images acquired at 3T. --------- The femora of five human volunteers were scanned using 1.5T and 3T MRI scanners from the same manufacturer (Siemens) with similar imaging protocols. A 3D flash sequence was used with TE = 4.66 ms, flip angle = 15° and voxel size = 0.5 × 0.5 × 1 mm. PA-Matrix and body matrix coils were used to cover the lower limb and pelvis respectively. Signal to noise ratio (SNR) [3] and contrast to noise ratio (CNR) [3] of the axial images from the proximal, shaft and distal regions were used to assess the quality of images from the 1.5T and 3T scanners. The SNR was calculated for the muscle and bone-marrow in the axial images. The CNR was calculated for the muscle to cortex and cortex to bone marrow interfaces, respectively. --------- Preliminary results (one volunteer) show that the SNR of muscle for the shaft and distal regions was higher in 3T images (11.65 and 17.60) than 1.5T images (8.12 and 8.11). For the proximal region the SNR of muscles was higher in 1.5T images (7.52) than 3T images (6.78). The SNR of bone marrow was slightly higher in 1.5T images for both proximal and shaft regions, while it was lower in the distal region compared to 3T images. The CNR between muscle and bone of all three regions was higher in 3T images (4.14, 6.55 and 12.99) than in 1.5T images (2.49, 3.25 and 9.89). The CNR between bone-marrow and bone was slightly higher in 1.5T images (4.87, 12.89 and 10.07) compared to 3T images (3.74, 10.83 and 10.15). These results show that the 3T images generated higher contrast between bone and the muscle tissue than the 1.5T images. It is expected that this improvement of image contrast will significantly reduce the time required for the mainly manual segmentation of the MR images. Future work will focus on optimizing the 3T imaging protocol for reducing chemical shift and susceptibility artifacts.

A biomechanical study of top screw pullout in anterior scoliosis correction constructs

Study Design: Biomechanical testing of vertebral body screw pullout resistance with relevance to top screw pullout in endoscopic anterior scoliosis constructs. Objectives: To analyse the effect of screw positioning and angulation on pullout resistance of vertebral body screws, where the pullout takes place along a curved path as occurs in anterior scoliosis constructs. Summary of Background Data: Top screw pullout is a significant clinical problem in endoscopic anterior scoliosis surgery, with rates of up to 18% reported in the literature. Methods: A custom designed biomechanical test rig was used to perform pullout tests of Medtronic anterior vertebral screws where the pullout occurred along an arc of known radius. Using synthetic bone blocks, a range of pullout radii and screw angulations were tested, in order to determine an ‘optimal’ configuration. The optimal configuration was then compared with standard screw positioning using a series of tests on ovine vertebrae (n=29). Results: Screw angulation has a small but significant effect on pullout resistance, with maximum strength being achieved at 10 degree cephalad angulation. Combining 10 degree cephalad angulation with maximal spacing between the top two screws (maximum pullout radius) increased the pullout resistance by 88% compared to ‘standard’ screw positioning (screws inserted perpendicular to rod at mid-body height). Conclusions: The positioning of the top screw in anterior scoliosis constructs can significantly alter its pullout resistance.

Design methods of rhombic tensegrity structures

As a special type of novel flexible structures, tensegrity holds promise for many potential applications in such fields as materials science, biomechanics, civil and aerospace engineering. Rhombic systems are an important class of tensegrity structures, in which each bar constitutes the longest diagonal of a rhombus of four strings. In this paper, we address the design methods of rhombic structures based on the idea that many tensegrity structures can be constructed by assembling one-bar elementary cells. By analyzing the properties of rhombic cells, we first develop two novel schemes, namely, direct enumeration scheme and cell-substitution scheme. In addition, a facile and efficient method is presented to integrate several rhombic systems into a larger tensegrity structure. To illustrate the applications of these methods, some novel rhombic tensegrity structures are constructed.

The effect of friction on indenter force and pile-up in numerical simulations of bone nanoindentation

Nanoindentation is a useful technique for probing the mechanical properties of bone, and finite element (FE) modeling of the indentation allows inverse determination of elasto-plastic constitutive properties. However, FE simulations to date have assumed frictionless contact between indenter and bone. The aim of this study was to explore the effect of friction in simulations of bone nanoindentation. Two dimensional axisymmetric FE simulations were performed using a spheroconical indenter of tip radius 0.6m and angle 90°. The coefficient of friction between indenter and bone was varied between 0.0 (frictionless) and 0.3. Isotropic linear elasticity was used in all simulations, with bone elastic modulus E=13.56GPa and Poisson’s ratio =0.3. Plasticity was incorporated using both Drucker-Prager and von Mises yield surfaces. Friction had a modest effect on the predicted force-indentation curve for both von Mises and Drucker-Prager plasticity, reducing maximum indenter displacement by 10% and 20% respectively as friction coefficient was increased from zero to 0.3 (at a maximum indenter force of 5mN). However, friction has a much greater effect on predicted pile-up after indentation, reducing predicted pile-up from 0.27m to 0.11m with a von Mises model, and from 0.09m to 0.02m with Drucker-Prager plasticity. We conclude that it is important to include friction in nanoindentation simulations of bone.

Development of patient-specific computer simulation tools for endoscopic scoliosis surgery

Endoscopic approaches for anterior correction of idiopathic scoliosis are a relatively new surgical technique. This paper describes the development of patient-specific finite element modelling techniques to investigate the biomechanics of single rod anterior scoliosis correction. Spinal geometry is obtained from pre-operative CT scans and material properties for osteo-ligamentous spinal tissues are based on existing literature. The techniques being developed will allow pre-surgical prediction of stresses, forces and deformations in spinal tissues, rods and screws under post-operative physiological loads.

Deconstructing neurobiological coordination

Inherent indeterminacy of neurobiological systems has been revealed by research on coordination of multiarticular actions. We consider three important issues that these investigations raise for biomechanical measurement and performance modeling. These issues highlight the role of dynamic systems theory as a platform for integration of motor control and biomechanics in exercise and sports science.

Modelling of inflatable rescue boats (IRBs) in surf conditions to reduce injuries

The Inflatable Rescue Boat (IRB) is arguably the most effective rescue tool used by the Australian surf lifesavers. The exceptional features of high mobility and rapid response have enabled it to become an icon on Australia's popular beaches. However, the IRB's extensive use within an environment that is as rugged as it is spectacular, has led it to become a danger to those who risk their lives to save others. Epidemiological research revealed lower limb injuries to be predominant, particularly the right leg. The common types of injuries were fractures and dislocations, as well as muscle or ligament strains and tears. The concern expressed by Surf Life Saving Queensland (SLSQ) and Surf Life Saving Australia (SLSA) led to a biomechanical investigation into this unique and relatively unresearched field. The aim of the research was to identify the causes of injury and propose processes that may reduce the instances and severity of injury to surf lifesavers during IRB operation. Following a review of related research, a design analysis of the craft was undertaken as an introduction to the craft, its design and uses. The mechanical characteristics of the vessel were then evaluated and the accelerations applied to the crew in the IRB were established through field tests. The data were then combined and modelled in the 3-D mathematical modelling and simulation package, MADYMO. A tool was created to compare various scenarios of boat design and methods of operation to determine possible mechanisms to reduce injuries. The results of this study showed that under simulated wave loading the boats flex around a pivot point determined by the position of the hinge in the floorboard. It was also found that the accelerations experienced by the crew exhibited similar characteristics to road vehicle accidents. Staged simulations indicated the attributes of an optimum foam in terms of thickness and density. Likewise, modelling of the boat and crew produced simulations that predicted realistic crew response to tested variables. Unfortunately, the observed lack of adherence to the SLSA footstrap Standard has impeded successful epidemiological and modelling outcomes. If uniformity of boat setup can be assured then epidemiological studies will be able to highlight the influence of implementing changes to the boat design. In conclusion, the research provided a tool to successfully link the epidemiology and injury diagnosis to the mechanical engineering design through the use of biomechanics. This was a novel application of the mathematical modelling software MADYMO. Other craft can also be investigated in this manner to provide solutions to the problem identified and therefore reduce risk of injury for the operators.

Modelling of some biological materials using continuum mechanics

Continuum mechanics provides a mathematical framework for modelling the physical stresses experienced by a material. Recent studies show that physical stresses play an important role in a wide variety of biological processes, including dermal wound healing, soft tissue growth and morphogenesis. Thus, continuum mechanics is a useful mathematical tool for modelling a range of biological phenomena. Unfortunately, classical continuum mechanics is of limited use in biomechanical problems. As cells refashion the �bres that make up a soft tissue, they sometimes alter the tissue's fundamental mechanical structure. Advanced mathematical techniques are needed in order to accurately describe this sort of biological `plasticity'. A number of such techniques have been proposed by previous researchers. However, models that incorporate biological plasticity tend to be very complicated. Furthermore, these models are often di�cult to apply and/or interpret, making them of limited practical use. One alternative approach is to ignore biological plasticity and use classical continuum mechanics. For example, most mechanochemical models of dermal wound healing assume that the skin behaves as a linear viscoelastic solid. Our analysis indicates that this assumption leads to physically unrealistic results. In this thesis we present a novel and practical approach to modelling biological plasticity. Our principal aim is to combine the simplicity of classical linear models with the sophistication of plasticity theory. To achieve this, we perform a careful mathematical analysis of the concept of a `zero stress state'. This leads us to a formal de�nition of strain that is appropriate for materials that undergo internal remodelling. Next, we consider the evolution of the zero stress state over time. We develop a novel theory of `morphoelasticity' that can be used to describe how the zero stress state changes in response to growth and remodelling. Importantly, our work yields an intuitive and internally consistent way of modelling anisotropic growth. Furthermore, we are able to use our theory of morphoelasticity to develop evolution equations for elastic strain. We also present some applications of our theory. For example, we show that morphoelasticity can be used to obtain a constitutive law for a Maxwell viscoelastic uid that is valid at large deformation gradients. Similarly, we analyse a morphoelastic model of the stress-dependent growth of a tumour spheroid. This work leads to the prediction that a tumour spheroid will always be in a state of radial compression and circumferential tension. Finally, we conclude by presenting a novel mechanochemical model of dermal wound healing that takes into account the plasticity of the healing skin.

Development and synthesis of innovative approaches to the surgical assessment of articular cartilage degeneration

This research has established, through ultrasound, near infrared spectroscopy and biomechanics experiments, parameters and parametric relationships that can form the framework for quantifying the integrity of the articular cartilage-on-bone laminate, and objectively distinguish between normal/healthy and abnormal/degenerated joint tissue, with a focus on articular cartilage. This has been achieved by: 1. using traditional experimental methods to produce new parameters for cartilage assessment; 2. using novel methodologies to develop new parameters; and 3. investigating the interrelationships between mechanical, structural and molec- ular properties to identify and select those parameters and methodologies that can be used in a future arthroscopic probe based on points 1 and 2. By combining the molecular, micro- and macro-structural characteristics of the tissue with its mechanical properties, we arrive at a set of critical benchmarking parameters for viable and early-stage non-viable cartilage. The interrelationships between these characteristics, examined using a multivariate analysis based on principal components analysis, multiple linear regression and general linear modeling, could then to deter- mine those parameters and relationships which have the potential to be developed into a future clinical device. Specifically, this research has found that the ultrasound and near infrared techniques can subsume the mechanical parameters and combine to characterise the tissue at the molecular, structural and mechanical levels over the full depth of the cartilage matrix. It is the opinion in this thesis that by enabling the determination of the precise area of in uence of a focal defect or disease in the joint, demarcating the boundaries of articular cartilage with dierent levels of degeneration around a focal defect, better surgical decisions that will advance the processes of joint management and treatment will be achieved. Providing the basis for a surgical tool, this research will contribute to the enhancement and quanti�cation of arthroscopic procedures, extending to post- treatment monitoring and as a research tool, will enable a robust method for evaluating developing (particularly focalised) treatments.