Biomechanical Alterations in Intact Osteoporotic Spine Due to Synthetic Augmentation: Finite Element Investigation

A parametric finite element analysis of the compacted bone–cement interface following balloon kyphoplasty

Treating fractures of the spine is a major challenge for the medical community with an estimated 1.4 million fractures per annum worldwide. While a considerable volume of study exists on the biomechanical implications of balloon kyphoplasty, which is used to treat these fractures, the influence of the compacted bone–cement region properties on stress distribution within the vertebral body remains unknown. The following article describes a novel method for modelling this compacted bone–cement region using a geometry-based approach in conjunction with the knowledge of the bone volume fractions for the native and compacted bone regions. Three variables for the compacted region were examined, as follows: (1) compacted thickness, (2) compacted region Young’s modulus and (3) friction coefficient. Results from the model indicate that the properties of the compacted bone–cement region can affect stresses in the cortical bone and cement by up to +28% and −40%, respectively. These findings demonstrate the need for further investigation into the effects of the compacted bone–cement interface using computational and experimental methods on multi-segment models.

FINITE ELEMENT ANALYSIS OF A MODEL OF SIMULATED VERTEBRAL CEMENT AUGMENTATION: INFLUENCE OF THE REPRESENTATION OF THE SHAPE OF THE CEMENT DOMAIN ON BIOMECHANICAL PARAMETERS

Vertebral cement augmentation is rapidly becoming the modality of choice for treating patients who are experiencing severe and persistent pain because of osteoporosis-induced vertebral compression fracture(s). The resulting cement domain (the part of the vertebral body (VB) filled with the cement) has an irregular or complicated shape. In literature reports of finite element analysis (FEA) of models of simulated vertebral cement augmentation, a variety of representations of the shape of the cement domain have been used. In the literature, only very limited attention has been given to the issue of the influence of cement domain shape representation on biomechanical parameters for a given combination of model and loading. This issue is the subject of the present work, with the model being of the L1-L3 motion segments. Augmentation of an unfractured L2 (prophylactic augmentation) was simulated, three cement domain shapes were considered — namely, solid cylinder, with rounded edges; two prolate spheroids; and oblate spheroid — and the applied loading comprised a simultaneous application of a uniform compressive pressure of 0.53 MPa (equivalent to an 800-N compression load) and a counter-clockwise-acting axial rotation moment of 1 Nm to the superior surface of L1. It was found that (1) while the cement domain shape representation has a marked influence on the mean von Mises stress (σAVM), the maximum von Mises stress (σMVM), and the strain energy density (MSED) distribution in the cement domain, its influence on each of these parameters in each of the biological tissues in the model as well as on the total segmental range of motion is minimal and (2) for σAVM and σMVM, the lowest value of each of these parameters was obtained when the oblate spheroid model was used. From both clinical and computational perspectives, these findings are significant. For example, the latter finding suggests that there is scope for researching the combination of key process variables used, such as the cement chemistry, the cement delivery system, and the augmentation technique/approach, that would ensure that the final cement domain shape in cement-augmented VBs of patients be oblate spheroid on a consistent and predictable basis.

In silico evaluation of stress distribution after vertebral body augmentation with conventional acrylics, composites and glass polyalkenoate cements

There exists clinical evidence of fractures in adjacent vertebrae subsequent to vertebral augmentation procedures, such as vertebroplasty (VP) and kyphoplasty (KP). A potential contributory factor to such fractures may be the excessive mismatch of mechanical properties between contemporary bone cements (i.e. polymethyl methacrylate (PMMA) and bisphenol-a-glycidyl dimethacrylate (BIS-GMA)) and bone. Aluminum-free glass polyalkenoate cements (GPCs) present an interesting alternative to conventional bone cements. GPCs adhere to the philosophy that implant materials should have mechanical characteristics similar to those of the bone, and also offer chemical adhesion and intrinsic bioactivity. However, their influence on the loading patterns of augmented vertebrae (as compared with conventional bone cements) is not available in the literature. The present work investigates how the moduli of PMMA, BIS-GMA and GPC implants affect the stress distribution within a single, augmented vertebra, in both healthy and osteoporotic states. Using a finite element model of the L4 vertebra derived from computed tomography data, with simulated augmentation, it was found that, as cement stiffness increased, stress was redistributed from the cortical and trabecular bone to the cement implant. The GPC implant exhibited the least effect on stress redistribution in both the healthy and osteoporotic models compared to its acrylic counterparts. The significance of this work is that, under simulated physiological loading conditions, aluminum-free GPCs exhibit stress distribution throughout the vertebral body similar to that of the healthy bone. In comparison to conventional augmentation materials, the use of aluminum-free GPCs in VP and KP may help to ameliorate the clinical complication of adjacent vertebral body compression fractures.

Finite element analysis of the spine: towards a framework of verification, validation and sensitivity analysis

A number of papers have recently emphasised the importance of verification, validation and sensitivity testing in computational studies within the field of biomechanical engineering. This review examines the methods used in the development of spinal finite element models with a view to a standardised framework of verification, validation and sensitivity analysis. The scope of this paper is restricted to models of the vertebra, the intervertebral disc and short spinal segments. In the case of single vertebral models, specimen-specific methods have been developed, which allow direct validation against experimental tests. The focus of intervertebral disc modelling has been on representing the complex material properties and further sensitivity testing is required to fully understand the relative roles of these input parameters. In order to construct complex multi-component short segment models, many geometric and material parameters are required, some of which are yet to be fully characterised. There are also major challenges in terms of short segment model validation. Throughout the review, areas of good practise are highlighted and recommendations for future development are proposed, taking a step towards more robust spinal modelling procedures, promoting acceptance from the wider biomechanics community.