Influence prediction of injury and vibration on adjacent components of spine using finite element methods

Finite element modelling and biodynamic response prediction of the seated human body exposed to whole-body vibration

Biodynamic modelling of seat-occupant systems can assist in seat comfort design. A finite element (FE) model of the seated human body, including detailed modelling of the lumbar spine, was established to reflect the human response to vibration and biodynamic response of the lumbar spine under whole-body vibration (WBV). The lumbar spine model was established and validated against the in-vitro results and calculated data. The posture of the lumbar spine was adjusted according to the radiological research results, and the adjusted model was combined to establish a FE model of the seated human body. The present seated human model with backrest inclination angles of 10, 20, and 30°, validated by comparing the measured apparent mass and seat-to-lumbar spine transmissibility, was used to calculate the biodynamic response of the lumbar spine with three inclined backrests under WBV. The results showed that the model could characterise the apparent mass, seat-to-lumbar spine transmissibility, and the biodynamic response of the lumbar spine. Practitioner summary: Biodynamic models can represent dynamic characteristics of the human body exposed to vibration and assist in seat comfort design. The three-dimensional FE model of the human body can be used to explore the human response to vibration and the biodynamic response of the lumbar spine under WBV.

Prediction of the biomechanical behaviour of the lumbar spine under multi-axis whole-body vibration using a whole-body finite element model

Low back pain has been reported to have a high prevalence among occupational drivers. Whole-body vibration during the driving environment has been found to be a possible factor leading to low back pain. Vibration loads might lead to degeneration and herniation of the intervertebral disc, which would increase incidence of low back problems among drivers. Some previous studies have reported the effects of whole-body vibration on the human body, but studies on the internal dynamic responses of the lumbar spine under multi-axis vibration are limited. In this study, the internal biomechanical response of the intervertebral disc was extracted to investigate the biomechanical behaviour of the lumbar spine under a multi-axial vibration in a whole-body environment. A whole-body finite element model, including skin, soft tissues, the bone skeleton, internal organs and a detailed ligamentous lumbar spine, was used to provide a whole-body condition for analyses. The results showed that both vibrations close to vertical and fore-and-aft resonance frequencies would increase the transmission of vibrations in the intervertebral disc, and vertical vibration might have a greater effect on the lumbar spine than fore-and-aft vibration. The larger deformation of the posterior region of the intervertebral disc in a multi-axis vibration environment might contribute to the higher susceptibility of the posterior region of the intervertebral disc to injury. The findings of this study revealed the dynamic behaviours of the lumbar spine in multi-axis vehicle vibration conditions, and suggested that both vertical and fore-and-aft vibration should be considered for protecting the lumbar health of occupational drivers.

Effect of whole-body vibration at different frequencies on the lumbar spine: A finite element study based on a whole human body model

Many previous studies have found that occupational drivers commonly suffered from low back pain, and low back pain and degeneration of the intervertebral disc might be associated with vibration conditions. However, the biomechanical mechanisms of whole-body vibration that caused pain and injury were not clear. In this study, a validated whole human body finite element model was used, and vibration loads at frequencies of 3, 5, 7 and 9 Hz were loaded to evaluate the frequency effects on the spine. The results showed that the responses of the spine were strong at the 5 Hz vibration load. Vibration loads would produce alternating stresses and bulges in the annulus fibrosus and change the direction of the pressure in the nucleus pulposus. The posterior region of the intervertebral disc showed greater stress fluctuations than the anterior region. The Risk Factors showed that long-term exposure to whole-body vibrations at 5 and 7 Hz might have greater adverse effects on the spine. The findings of this study confirmed that vibrations near the resonance frequency of the human body would cause more injuries to the spine than other frequencies. Alternating stress and bulge might cause fatigue and the degeneration of the intervertebral disc, which might be the mechanisms of spinal injury caused by whole-body vibration, and the posterior regions of the intervertebral disc were more susceptible to degeneration. Some appropriate measures should be taken to reduce the adverse effects of whole-body vibration on spinal health.

Computational Modeling Intervertebral Disc Pathophysiology: A Review

Lower back pain is a medical condition of epidemic proportion, and the degeneration of the intervertebral disc has been identified as a major contributor. The etiology of intervertebral disc (IVD) degeneration is multifactorial, depending on age, cell-mediated molecular degradation processes and genetics, which is accelerated by traumatic or gradual mechanical factors. The complexity of such intertwined biochemical and mechanical processes leading to degeneration makes it difficult to quantitatively identify cause–effect relationships through experiments. Computational modeling of the IVD is a powerful investigative tool since it offers the opportunity to vary, observe and isolate the effects of a wide range of phenomena involved in the degenerative process of discs. This review aims at discussing the main findings of finite element models of IVD pathophysiology with a special focus on the different factors contributing to physical changes typical of degenerative phenomena. Models presented are subdivided into those addressing role of nutritional supply, progressive biochemical alterations stemming from an imbalance between anabolic and catabolic processes, aging and those considering mechanical factors as the primary source that induces morphological change within the disc. Limitations of the current models, as well as opportunities for future computational modeling work are also discussed.

Deleterious effects of whole-body vibration on the spine: A review of in vivo, ex vivo, and in vitro models

Occupational exposure to whole-body vibration is associated with the development of musculoskeletal, neurological, and other ailments. Low back pain and other spine disorders are prevalent among those exposed to whole-body vibration in occupational and military settings. Although standards for limiting exposure to whole-body vibration have been in place for decades, there is a lack of understanding of whole-body vibration-associated risks among safety and healthcare professionals. Consequently, disorders associated with whole-body vibration exposure remain prevalent in the workforce and military. The relationship between whole-body vibration and low back pain in humans has been established largely through cohort studies, for which vibration inputs that lead to symptoms are rarely, if ever, quantified. This gap in knowledge highlights the need for the development of relevant in vivo, ex vivo, and in vitro models to study such pathologies. The parameters of vibrational stimuli (eg, frequency and direction) play critical roles in such pathologies, but the specific cause-and-effect relationships between whole-body vibration and spinal pathologies remain mostly unknown. This paper provides a summary of whole-body vibration parameters; reviews in vivo, ex vivo, and in vitro models for spinal pathologies resulting from whole-body vibration; and offers suggestions to address the gaps in translating injury biomechanics data to inform clinical practice.

A comparison of the influence of three different lumbar interbody fusion approaches on stress in the pedicle screw fixation system: Finite element static and vibration analyses

This study aimed to examine breakage risk of the bilateral pedicle screw (BPS) fixation system under static and vibration loadings after three different types of lumbar interbody fusion surgery. A previously validated intact L1-sacrum finite element model was modified to simulate anterior, posterior, and transforaminal lumbar interbody fusion (ALIF, PLIF, and TLIF, respectively) with BPS fixation system (consisting of pedicle screws and rigid connecting rods) at L4-L5. As a risk parameter for breakage, the von Mises stresses in the pedicle screws and the rods for the ALIF, PLIF, and TLIF models under static loading (flexion, extension, lateral bending, and axial torsion moments) and vibration loading (sinusoidal vertical load) were calculated and compared. The calculated von Mises stresses were different in the ALIF, PLIF, and TLIF models, but these stresses for all the fusion models were found to be concentrated in neck of the pedicle screw and middle of the rod under both the static and vibration loadings. The results from static analyses showed that the maximum stress in the BPS fixation system was greater in the TLIF model than in the ALIF and PLIF models under all the applied static loadings. The results from transient dynamic analyses also showed that the TLIF generated greater dynamic responses of the stress in the BPS fixation system to the vertical vibration compared with the ALIF and PLIF. It implies that the TLIF procedure might incur a higher risk of breakage for the BPS fixation system than the ALIF and PLIF procedures.

Validation of an In Vivo Medical Image-Based Young Human Lumbar Spine Finite Element Model

Mathematical models of the human spine can be used to investigate spinal biomechanics without the difficulties, limitations, and ethical concerns associated with physical experimentation. Validation of such models is necessary to ensure that the modeled system behavior accurately represents the physics of the actual system. The goal of this work was to validate a medical image-based nonlinear lumbosacral spine finite element model of a healthy 20-yr-old female subject under physiological moments. Range of motion (ROM), facet joint forces (FJF), and intradiscal pressure (IDP) were compared with experimental values and validated finite element models from the literature. The finite element model presented in this work was in good agreement with published experimental studies and finite element models under pure moments. For applied moments of 7.5 N·m, the ROM in flexion–extension, axial rotation, and lateral bending were 39 deg, 16 deg, and 28 deg, respectively. Excellent agreement was observed between the finite element model and experimental data for IDP under pure compressive loading. The predicted FJFs were lower than those of the experimental results and validated finite element models for extension and torsion, likely due to the nondegenerate properties chosen for the intervertebral disks and morphology of the young female spine. This work is the first to validate a computational lumbar spine model of a young female subject. This model will serve as a valuable tool for predicting orthopedic spinal injuries, studying the effect of intervertebral disk replacements using advanced biomaterials, and investigating soft tissue degeneration.

Influence of different frequencies of axial cyclic loading on time-domain vibration response of the lumbar spine: A finite element study

Very few studies have quantitatively analyzed influence of the loading frequency on time-domain vibration response of the whole lumbar spine in the presence of a physiologic compressive preload. In this study, a three-dimensional non-linear finite element model of ligamentous L1-S1 segment was developed to predict time-domain dynamic response of the whole lumbar spine to axial cyclic loading with different frequencies. A compressive follower preload of 400 N was applied to the model to simulate the physiologic compressive load. Modal analysis was initially performed to extract axial resonant frequency of the model under a 40 kg upper body mass and the 400 N preload. The result showed that the axial resonant frequency was 7.77 Hz. Subsequently, transient dynamic analyses were performed on the model under a sinusoidal axial load of ±40 N at frequencies of 3, 5, 7, 9, 11, 13 and 15 Hz with the 400 N preload and 40 kg mass. The computational results (strains and stresses in the spinal components) were collected and plotted as a function of time. These predicted results were found to be frequency-dependent and consistent with the notion in engineering dynamics texts that the closer the loading frequency approaches the resonant frequency, the larger the response is. For example, the results for 5 Hz load compared to 3 Hz load showed a 68.6-111.5% increase in peak-to-bottom variations of the predicted response parameters, and the results for 13 Hz load compared to 11 Hz load showed a 26.4-37.8% decrease in these variations.

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.

Restoration of compressive loading properties of lumbar discs with a nucleus implant-a finite element analysis study

BACKGROUND CONTEXT

Discectomy is a common procedure for treating sciatica. However, both the operation and preceding herniated disc alter the biomechanical properties of the spinal segment. The disc mechanics are also altered in patients with chronic contained herniation. The biomechanical properties of the disc can potentially be restored with an elastomeric nucleus replacement implanted via minimally invasive surgery.

PURPOSE

The purpose of this study was to determine whether the compressive characteristics of the intervertebral disc after a nucleotomy can be restored with an elastomeric nucleus replacement.

STUDY DESIGN

A finite element model of the L4-L5 intervertebral disc was created to investigate the effect of the implantation of an elastomeric nucleus replacement on the biomechanical properties of the disc under axial loading.

METHOD

A L4-L5 physiologic intervertebral disc model was constructed and then modified to contain a range by volume of nucleotomies and nucleus replacements. The material properties of the nucleus replacement were based on experimental data for an elastomeric implant. The compressive stiffness, radial annular bulge, and stress distribution of the nucleotomy and nucleus replacement models were investigated under displacement-controlled loading.

RESULTS

Removal of nucleus pulposus from the physiologic disc reduced the force necessary to compress the disc 2 mm by 50%, altered the von Mises stress distribution, and reduced the outward radial annular bulge. Replacing the natural nucleus pulposus of the physiologic disc with an artificial nucleus reduced the force required to compress the disc 2 mm by 10%, indicating a restoration of disc compressive stiffness. The von Mises stress distribution and annular bulge observed in the disc with an artificial nucleus were similar to that observed in the physiologic disc.

CONCLUSION

This study demonstrates that despite having different material properties, a nucleus replacement implant can restore the axial compressive mechanical properties of a disc after a discectomy. The implant carries compressive load and transfers the load into annular hoop stress.

The immediate effect of repeated loading on the compressive strength of young porcine lumbar spine

The human spine is exposed to repeated loading during daily activities and more extremely during sports. Despite this, there remains a lack of knowledge regarding the immediate effects on the spine due to this mode of loading. Age-specific spinal injury patterns has been demonstrated and this implies differences in reaction to load mode and load history The purpose of the present study was to investigate the impact of cyclic pre-loading on the biomechanical properties and fracture patterns of the adolescent spine in an experimental model. Eight functional spinal units from four young porcine spines were harvested. The functional spinal units were cyclic loaded with 20,000 cycles and then axially compressed to failure. The compression load at failure, ultimate stress and viscoelastic parameters were calculated. The functional spinal units were examined with plain radiography, computer tomography and MRI before and after the loading, and finally macroscopically and histologically. The median compression load at failure in this study was 8.3 kN (range 5.6-8.7 kN). The median deformation for all cases was 2.24 mm (range 2.30-2.7 mm) and stiffness was 3.45 N/mm (range 3.5-4.5 N/mm). A fracture was seen on radiograph in one case, on CT and macroscopically in seven, and on MRI and histologically in all eight cases. The cyclic loaded functional spinal units in the present study were not more sensitive to axial compression than non-cyclic loaded functional spinal units from young porcine. The endplate and the growth zone were the weakest part in the cyclic loaded functional spinal units. Disc signal reduction and disc height reduction was found on MRI. The E-modulus value found in this study was of the same order of magnitude as found by others using a porcine animal model.

Vibration modes of injured spine at resonant frequencies under vertical vibration

STUDY DESIGN

A detailed three-dimensional finite element model of the spine segment T12-Pelvis was developed to investigate dynamic characteristics of whole lumbar spine with injured cases.

OBJECTIVE

This study investigates the motion mechanism of the human lumbar spine and the effect of component injuries on adjacent spinal components under whole body vibration.

SUMMARY OF BACKGROUND DATA

Several investigations have analyzed the influence of injured spines on adjacent spinal components under static loadings. However, it is not clear how the spine injury affects dynamic characteristics of whole lumbar spine and adjacent components of the injured segment under vibration.

METHODS

The T12-Pelvis model was used to obtain the modal vibration modes of the spine at resonant frequencies. Injury conditions of the spine were simulated and tested, including denucleation and/or facetectomy with removal of capsular ligaments.

RESULTS

The results indicate the first-order vertical resonant frequency of the intact model is 7.21 Hz. After the denucleation at L4-L5, it decreases by more than 4% compared with the intact condition. All the injured conditions including disc injury and ligament injury decrease the resonant frequency of the spine. Due to the denucleation at L4-L5 the anteroposterior displacements of the vertebrae from L2 to L5 decrease and the vertical displacements of the vertebrae from L1 to L4 increase under vibration. The denucleation also decreases the rotational deformations of the vertebrae from L1 to L5. The material property sensitivity analysis shows intervertebral discs have a dominating effect on variation of vertical resonant frequency of the spine.

CONCLUSION

The denucleation may decrease cushioning effects of adjacent motion segments at the injured level under vibration. The injured condition may increase the vertical displacement amplitudes of the spine above the injured level. All the injured conditions may decrease the resonant frequency of the spine system.

Early stage disc degeneration does not have an appreciable affect on stiffness and load transfer following vertebroplasty and kyphoplasty

Vertebroplasty and kyphoplasty have been reported to alter the mechanical behavior of the treated and adjacent-level segments, and have been suggested to increase the risk for adjacent-level fractures. The intervertebral disc (IVD) plays an important role in the mechanical behavior of vertebral motion segments. Comparisons between normal and degenerative IVD motion segments following cement augmentation have yet to be reported. A microstructural finite element model of a degenerative IVD motion segment was constructed from micro-CT images. Microdamage within the vertebral body trabecular structure was used to simulate a slightly (I = 83.5% of intact stiffness), moderately (II = 57.8% of intact stiffness), and severely (III = 16.0% of intact stiffness) damaged motion segment. Six variable geometry single-segment cement repair strategies (models A-F) were studied at each damage level (I-III). IVD and bone stresses, and motion segment stiffness, were compared with the intact and baseline damage models (untreated), as well as, previous findings using normal IVD models with the same repair strategies. Overall, small differences were observed in motion segment stiffness and average stresses between the degenerative and normal disc repair models. We did however observe a reduction in endplate bulge and a redistribution in the microstructural tissue level stresses across both endplates and in the treated segment following early stage IVD degeneration. The cement augmentation strategy placing bone cement along the periphery of the vertebra (model E) proved to be the most advantageous in treating the degenerative IVD models by showing larger reductions in the average bone stresses (vertebral and endplate) as compared to the normal IVD models. Furthermore, only this repair strategy, and the complete cement fill strategy (model F), were able to restore the slightly damaged (I) motion segment stiffness above pre-damaged (intact) levels. Early stage IVD degeneration does not have an appreciable effect in motion segment stiffness and average stresses in the treated and adjacent-level segments following vertebroplasty and kyphoplasty. Placing bone cement in the periphery of the damaged vertebra in a degenerative IVD motion segment, minimizes load transfer, and may reduce the likelihood of adjacent-level fractures.