Effect of impact velocity and ligament mechanical properties on lumbar spine injuries in posterior-anterior impact loading conditions: a finite element study

Biomechanical Study of Symmetric Bending and Lifting Behavior in Weightlifter with Lumbar L4-L5 Disc Herniation and Physiological Straightening Using Finite Element Simulation

BACKGROUND

Physiological curvature changes of the lumbar spine and disc herniation can cause abnormal biomechanical responses of the lumbar spine. Finite element (FE) studies on special weightlifter models are limited, yet understanding stress in damaged lumbar spines is crucial for preventing and rehabilitating lumbar diseases. This study analyzes the biomechanical responses of a weightlifter with lumbar straightening and L4-L5 disc herniation during symmetric bending and lifting to optimize training and rehabilitation.

METHODS

Based on the weightlifter's computed tomography (CT) data, an FE lumbar spine model (L1-L5) was established. The model included normal intervertebral discs (IVDs), vertebral endplates, ligaments, and a degenerated L4-L5 disc. The bending angle was set to 45°, and weights of 15 kg, 20 kg, and 25 kg were used. The flexion moment for lifting these weights was theoretically calculated. The model was tilted at 45° in Abaqus 2021 (Dassault Systèmes Simulia Corp., Johnston, RI, USA), with L5 constrained in all six degrees of freedom. A vertical load equivalent to the weightlifter's body mass and the calculated flexion moments were applied to L1 to simulate the weightlifter's bending and lifting behavior. Biomechanical responses within the lumbar spine were then analyzed.

RESULTS

The displacement and range of motion (ROM) of the lumbar spine were similar under all three loading conditions. The flexion degree increased with the load, while extension remained unchanged. Right-side movement and bending showed minimal change, with slightly more right rotation. Stress distribution trends were similar across loads, primarily concentrated in the vertebral body, increasing with load. Maximum stress occurred at the anterior inferior margin of L5, with significant stress at the posterior joints, ligaments, and spinous processes. The posterior L5 and margins of L1 and L5 experienced high stress. The degenerated L4-L5 IVD showed stress concentration on its edges, with significant stress also on L3-L4 IVD. Stress distribution in the lumbar spine was uneven.

CONCLUSIONS

Our findings highlight the impact on spinal biomechanics and suggest reducing anisotropic loading and being cautious of loaded flexion positions affecting posterior joints, IVDs, and vertebrae. This study offers valuable insights for the rehabilitation and treatment of similar patients.

Subject-Specific Geometry of FE Lumbar Spine Models for the Replication of Fracture Locations Using Dynamic Drop Tests

For traumatic lumbar spine injuries, the mechanisms and influence of anthropometrical variation are not yet fully understood under dynamic loading. Our objective was to evaluate whether geometrically subject-specific explicit finite element (FE) lumbar spine models based on state-of-the-art clinical CT data combined with general material properties from the literature could replicate the experimental responses and the fracture locations via a dynamic drop tower-test setup. The experimental CT datasets from a dynamic drop tower-test setup were used to create anatomical details of four lumbar spine models (T12 to L5). The soft tissues from THUMS v4.1 were integrated by morphing. Each model was simulated with the corresponding loading and boundary conditions from the dynamic lumbar spine tests that produced differing injuries and injury locations. The simulations resulted in force, moment, and kinematic responses that effectively matched the experimental data. The pressure distribution within the models was used to compare the fracture occurrence and location. The spinal levels that sustained vertebral body fracture in the experiment showed higher simulation pressure values in the anterior elements than those in the levels that did not fracture in the reference experiments. Similarly, the spinal levels that sustained posterior element fracture in the experiments showed higher simulation pressure values in the vertebral posterior structures compared to those in the levels that did not sustain fracture. Our study showed that the incorporation of the spinal geometry and orientation could be used to replicate the fracture type and location under dynamic loading. Our results provided an understanding of the lumbar injury mechanisms and knowledge on the load thresholds that could be used for injury prediction with explicit FE lumbar spine models.

Biomechanical evaluation of Back injuries during typical snowboarding backward falls

To prevent spinal and back injuries in snowboarding, back protector devices (BPDs) have been increasingly used. The biomechanical knowledge for the BPD design and evaluation remains to be explored in snowboarding accident conditions. This study aims to evaluate back-to-snow impact conditions and the associated back injury mechanisms in typical snowboarding backward falls. A previously validated snowboarder multi-body model was first used to evaluate the impact zones on the back and the corresponding impact velocities in a total of 324 snowboarding backward falls. The biomechanical responses during back-to-snow impacts were then evaluated by applying the back-to-snow impact velocity to a full human body finite element model to fall on the snow ground of three levels of stiffness (soft, hard, and icy snow). The mean values of back-to-snow normal and tangential impact velocities were 2.4 m/s and 7.3 m/s with maximum values up to 4.8 m/s and 18.5 m/s. The lower spine had the highest normal impact velocity during snowboarding backward falls. The thoracic spine was found more likely to exceed the limits of flexion-extension range of motions than the lumbar spine during back-to-snow impacts, indicating a higher injury risk. On the hard and icy snow, rib cage and vertebral fractures were predicted at the costal cartilage and the posterior elements of the vertebrae. Despite the possible back injuries, the back-to-snow impact force was always lower than the force thresholds of the current BPD testing standard. The current work provides additional biomechanical knowledge for the future design of back protections for snowboarders.

The mechanical behavior of bovine spinal cord white matter under various strain rate conditions: tensile testing and visco-hyperelastic constitutive modeling

The mechanical behavior of the white matter is important for estimating the damage of the spinal cord during accidents. In this study, we conducted uniaxial tension testing in vitro of bovine spinal cord white matter under extremely high strain rate conditions (up to 100 s^-1). A visco-hyperelastic constitutive law for modeling the strain rate-dependent behavior of the bovine spinal cord white matter was developed. A set of material constants was obtained using a Levenberg-Marquardt fitting algorithm to match the uniaxial tension experimental data with various strain rates. Our experimental data confirmed that the modulus and tensile strength increased when the strain rate is higher. For the extremely high strain rate condition (100 s^-1), we found that both the modulus and failure stress significantly increased compared with the low strain rate case. These new data in terms of mechanical response at high strain rate provide insight into the spine injury mechanism caused by high-speed impact. Moreover, the developed constitutive model will allow researchers to perform more realistic finite element modeling and simulation of spinal cord injury damage under various complicated conditions.

Finite Element Method for the Evaluation of the Human Spine: A Literature Overview

The finite element method (FEM) represents a computer simulation method, originally used in civil engineering, which dates back to the early 1940s. Applications of FEM have also been used in numerous medical areas and in orthopedic surgery. Computing technology has improved over the years and as a result, more complex problems, such as those involving the spine, can be analyzed. The spine is a complex anatomical structure that maintains the erect posture and supports considerable loads. Applications of FEM in the spine have contributed to the understanding of bone biomechanics, both in healthy and abnormal conditions, such as scoliosis, fractures (trauma), degenerative disc disease and osteoporosis. However, since FEM is only a digital simulation of the real condition, it will never exactly simulate in vivo results. In particular, when it concerns biomechanics, there are many features that are difficult to represent in a FEM. More FEM studies and spine research are required in order to examine interpersonal spine stiffness, young spine biomechanics and model accuracy. In the future, patient-specific models will be used for better patient evaluations as well as for better pre- and inter-operative planning.

The importance of intervertebral disc material model on the prediction of mechanical function of the cervical spine

BACKGROUND

Linear elastic, hyperelastic, and multiphasic material constitutive models are frequently used for spinal intervertebral disc simulations. While the characteristics of each model are known, their effect on spine mechanical response requires a careful investigation. The use of advanced material models may not be applicable when material constants are not available, model convergence is unlikely, and computational time is a concern. On the other hand, poor estimations of tissue's mechanical response are likely if the spine model is oversimplified. In this study, discrepancies in load response introduced by material models will be investigated.

METHODS

Three fiber-reinforced C2-C3 disc models were developed with linear elastic, hyperelastic, and biphasic behaviors. Three different loading modes were investigated: compression, flexion and extension in quasi-static and dynamic conditions. The deformed disc height, disc fluid pressure, range of motion, and stresses were compared.

RESULTS

Results indicated that the intervertebral disc material model has a strong effect on load-sharing and disc height change when compression and flexion were applied. The predicted mechanical response of three models under extension had less discrepancy than its counterparts under flexion and compression. The fluid-solid interaction showed more relevance in dynamic than quasi-static loading conditions. The fiber-reinforced linear elastic and hyperelastic material models underestimated the load-sharing of the intervertebral disc annular collagen fibers.

CONCLUSION

This study confirmed the central role of the disc fluid pressure in spinal load-sharing and highlighted loading conditions where linear elastic and hyperelastic models predicted energy distribution different than that of the biphasic model.

The Effect of Thoracolumbar Injury Classification in the Clinical Outcome of Operative and Non-Operative Treatments

This review assesses the validity of a biomechanical approach using finite element analysis in the Thoracolumbar Injury Classification and Severity Score System (TLICS) by addressing the “gray zone” decision discrepancy of thoracolumbar spinal injuries. A systematic review was performed using the keywords “Thoracolumbar Injury Classification” AND “finite element analysis of the spinal column” to evaluate the validity of the TLICS and finite element analysis of the thoracolumbar spinal column. Results were classified according to the main conclusions and level of evidence. Thirteen articles are included. Four of the articles evaluated the TLICS in comparison to other classification systems of thoracolumbar spinal injuries. A notable finding is that the TLICS had inconsistencies with other classification systems in the treatment of burst fractures without neurological deficits. One article evaluated the TLICS with the inclusion of magnetic resonance imaging (MRI) in the evaluation, which decreased the agreement between the suggested and actual treatment. Among the three finite element analysis studies, limited data have been published on the posterior ligamentous complex (PLC) status when an injury is suspected or indeterminate. The TLICS has been a reliable classification system in the management of single-column fractures and three-column injuries treated with surgical stabilization. Special attention to enhancing the TLICS classification system by eliminating the “gray zone” of a TLICS score of 4 is essential. Biomedical computational modeling evaluating the PLC status of indeterminate or injury suspected is needed to enhance the current TLICS system and to clarify the decision discrepancy in the “gray zone.”

Fiber splay precludes the direct identification of ligament material properties: Implications for ACL graft selection

Anterior cruciate ligament (ACL) injuries typically require surgical reconstruction to restore adequate knee stability. The middle third of an injured patient's patellar tendon (PT) is a commonly used graft for ACL reconstruction. However, many clinicians and researchers question whether it is the best option, as several studies have suggested that it is a stiffer material than the ACL. Still, there is little to no consensus on even the most basic material property of ligaments/tendons: the tangent modulus in the fiber direction, or slope of the linear portion of the uniaxial stress-strain curve. In this study, we investigate the effect of fiber splay (the tendency of collagen fibers to spread out near the enthesis) on the apparent tangent modulus. Using a simplified theoretical model, we establish a quantity we call the splay ratio, which describes the relationship between splay geometry and the apparent tangent modulus. We then more rigorously investigate the effect of the splay ratio on the apparent tangent modulus of the ovine PT and anteromedial and posterolateral regions of the ACL using experimental and computational methods. Both approaches confirmed that splay geometry significantly affects the apparent material behavior. Because true material properties are independent of geometry, we conclude that the macroscopic response of ligaments and tendons is not sufficient for the characterization of their material properties, but rather is reflective of both material and structural properties. We further conclude that the PT is probably not a stiffer material than ACL, but that the PT graft is likely a stiffer structure than either ACL region.

Biomechanical modelling of the facet joints: a review of methods and validation processes in finite element analysis

There is an increased interest in studying the biomechanics of the facet joints. For in silico studies, it is therefore important to understand the level of reliability of models for outputs of interest related to the facet joints. In this work, a systematic review of finite element models of multi-level spinal section with facet joints output of interest was performed. The review focused on the methodology used to model the facet joints and its associated validation. From the 110 papers analysed, 18 presented some validation of the facet joints outputs. Validation was done by comparing outputs to literature data, either computational or experimental values; with the major drawback that, when comparing to computational values, the baseline data was rarely validated. Analysis of the modelling methodology showed that there seems to be a compromise made between accuracy of the geometry and nonlinearity of the cartilage behaviour in compression. Most models either used a soft contact representation of the cartilage layer at the joint or included a cartilage layer which was linear elastic. Most concerning, soft contact models usually did not contain much information on the pressure-overclosure law. This review shows that to increase the reliability of in silico model of the spine for facet joints outputs, more needs to be done regarding the description of the methods used to model the facet joints, and the validation for specific outputs of interest needs to be more thorough, with recommendation to systematically share input and output data of validation studies.

Simulation analysis of impact damage to the bone tissue surrounding a dental implant

Dental implant may suffer transient external impacts. To simulate the effect of impact forces on bone damage is very important for evaluation of damage and guiding treatment in clinics. In this study, an animal model was established by inserting an implant into the femoral condyle of New Zealand rabbit. Implant with good osseointegration was loaded with impact force. A three-dimensional finite element model was established based on the data of the animal model. Damage process to bone tissue was simulated with Abaqus 6.13 software combining dynamic mechanical properties of the femur. The characteristics of bone damage were analyzed by comparing the results of animal testing with numerical simulation data. After impact, cortical bone around the implant and trabecular at the bottom of the implant were prone to damage. The degree of damage correlated with the direction of loading and the magnitude of the impact. Lateral loading was most likely performed to damage cancellous bone. The stress wave formed by the impact force can damage the implant–bone interface and peri-implant trabeculae. The data from numerical simulations were consistent with data from animal experiments, highlighting the importance of a thorough examination and evaluation based on the patient’s medical history.

Subject-specific finite element analysis of a lumbar cage produced by electron beam melting

The aim of this study was the analysis of the mechanical behaviour of a partially porous lumbar custom-made cage by means of a subject-specific finite element analysis (FEA). The cage, made of Ti6Al4V ELI alloy, was produced via electron beam melting (EBM) process and surgically implanted in a female subject, 50 years old. The novelty of this study was the customized design of the cage and of its internal structure, which is impossible to obtain with the traditional production techniques. The 3D model of the spine was obtained from the computed tomography (CT) of the patient. Moreover, high-resolution industrial CT was also used to reconstruct a 3D model of the cage, with its real (as-produced) features, such as superficial roughness, morphology of the bulk and of the porous structure. The workflow was divided in several steps: the main finite element analyses were non-linear and quasi-static regarding: the rhombic dodecahedron (RD) unit cell of the porous structure; the device; the whole L4-L5 motion segment with the implanted cage. Stress distribution was calculated under compression load for all models. For the RD unit cell, the maximum stress appeared at the connected cross nodes, where notch effect was present. For the cage subjected to a load of 1 kN, the porous structure did not present any functional failure. For the whole biomechanical system subjected to a physiological load of 360 N, the calculated stress in the bone was smaller than its yield strength value. On the axial view, a zone with higher compressive stresses was present on the L5 vertebral body. This was due to the contact stress between the cage and the vertebra. From the comparison between FE results and the CT images of the spine, bone remodelling was supposed, with the formation of new bone. Graphical abstract Workflow showing the phases of the research.