Finite element modeling and static/dynamic validation of thoracolumbar-pelvic segment

Multiscale dynamics analysis of lumbar vertebral cortical bone based on the Abaqus submodel finite element method

The direct effect of macroscopic loads on the microstructure of bone tissue in a vibration environment is not yet known. Therefore, this study aims to investigate the macro- and micro-biomechanical properties of the lumbar spine system under dynamic loading in such an environment. We analyzed the dynamic characteristics of osteon by establishing a macro- and micro-scale model of the lumbar spine, using a submodel-specific boundary displacement method based on St. Venant's principle. Utilizing the results from the transient dynamic analysis of the entire lumbar spine as boundary conditions, this study simulates the dynamic behavior of osteon in each segment of the spine on a microscopic scale. The macroscopic results of the transient dynamic analysis showed that the rates of change in dynamic displacement amplitude relative to static displacement amplitude for the L1-L5 vertebrae were 212.60%, 242.11%, 314.80%, 1.17%, and 3.75%, respectively. The change in displacement amplitude under dynamic load relative to static load was highest for the L3 vertebra, as observed in the macroscopic model. The stress and strain values in the microscopic osteon of each lumbar spine segment under sinusoidal periodic loading were higher than those in the macroscopic osteon. In the microscopic bone unit, the maximum stress occurred at the cement line during the peak stress moment, while the minimum stress was observed at the innermost bone plate during the moment of minimum stress. Under dynamic loading, the microscopic bone osteon demonstrated a cyclic stress and strain response, with variations observed in different components of the osteon. These findings provide new insights into the biomechanical behavior of the lumbar spine in a vibration environment.

The impact of disc degeneration on the dynamic characteristics of the lumbar spine: a finite element analysis

The dynamics of disc degeneration was analyzed to determine the effect of disc degeneration at the L4-L5 segment on the dynamic characteristics of the total lumbar spine. A three-dimensional nonlinear finite element model of the L1-S1 normal lumbar spine was constructed and validated. This normal model was then modified to construct two degeneration models with different degrees of degeneration (mild, moderate) at the L4-L5 level. Modal analysis, harmonic response analysis, and transient dynamics analysis were performed on the total lumbar spine when experiencing following compressive loading (500 N). As the degree of disc degeneration increased, the vibration patterns corresponding to the first three orders of the model's intrinsic frequency were basically unchanged, with the first order being in the left-right lateral bending direction, the second order being in the forward-flexion and backward-extension direction, and the third order being in the axial stretching direction. The nucleus pulposus pressure peaks corresponding to the first order intrinsic frequency for the harmonic response analysis are all on the right side of the model, with sizes of 0.053 MPa, 0.061 MPa, and 0.036 MPa, respectively; the nucleus pulposus pressure peaks corresponding to the second order intrinsic frequency are all at the rear of the model, with sizes of 0.13 MPa, 0.087 MPa, and 0.11 MPa, respectively; and the nucleus pulposus pressure peaks corresponding to the third order intrinsic frequency are all at the front of the model, with sizes of 0.19 MPa, 0.22 MPa, and 0.22 MPa, respectively. The results of the transient analysis indicated that over time, the response curves of the healthy model, the mild model, and the moderate model all exhibited cyclic response characteristics. Intervertebral disc degeneration did not adversely affect the vibration characteristics of the entire lumbar spine system. Intervertebral disc degeneration significantly altered the dynamics of the degenerative segments and their neighboring normal segments. The process of disc degeneration gradually shifted the load from the nucleus pulposus to the annulus fibrosus when the entire lumbar spine was subjected to the same vibratory environment.

Finite element method-based study for spinal vibration characteristics of the scoliosis and kyphosis lumbar spine to whole-body vibration under a compressive follower preload

PURPOSE

To analyze the dynamic response of the lumbosacral vertebrae structure of a scoliosis spine and a kyphosis spine under whole-body vibration.

METHODS

Typical Lenke4 (kyphosis) and Lenke3 (scoliosis) spinal columns were used as research objects. A finite element model of the lumbosacral vertebrae segment was established and validated based on CT scanning images. Modal, harmonic response, and transient dynamic analyses were performed on the lumbar-sacral scoliosis model using the finite element software abaqus.

RESULTS

The first four resonance frequencies of kyphosis spine extracted from modal analysis were 0.86, 1.45, 8.51, and 55.71 Hz. The first four resonance frequencies of scoliosis spine extracted from modal analysis were 0.76, 1.45, 10.51, and 63.82 Hz. The scoliosis spine had the maximum resonance amplitude in the transverse direction, while the kyphosis spine had the maximum resonance amplitude in the anteroposterior direction. The dynamic response in transient analysis exhibited periodic response over time at all levels.

CONCLUSION

The scoliosis and kyphosis deformity of the spine significantly complicates the vibration response in the scoliosis and kyphosis areas at the top of the spine. Scoliosis and kyphosis patients are more likely to experience vibrational spinal diseases than healthy people. Besides, applying vertical cyclic loads on a malformed spine may cause further rotation of scoliosis and kyphosis deformities.

Biomechanical study of different bone cement distribution on osteoporotic vertebral compression Fracture-A finite element analysis

Purpose

This study aimed to compare the biomechanical effects of different bone cement distribution methods on osteoporotic vertebral compression fractures (OVCF).

Patients and methods

Raw CT data from a healthy male volunteer was used to create a finite element model of the T12-L2 vertebra using finite element software. A compression fracture was simulated in the L1 vertebra, and two forms of bone cement dispersion (integration group, IG, and separation group, SG) were also simulated. Six types of loading (flexion, extension, left/right bending, and left/right rotation) were applied to the models, and the stress distribution in the vertebra and intervertebral discs was observed. Additionally, the maximum displacement of the L1 vertebra was evaluated.

Results

Bone cement injection significantly reduced stress following L1 vertebral fractures. In the L1 vertebral body, the maximum stress of SG was lower than that of IG during flexion, left/right bending, and left/right rotation. In the T12 vertebral body, compared with IG, the maximum stress of SG decreased during flexion and right rotation. In the L2 vertebral body, the maximum stress of SG was the lowest under all loading conditions. In the T12-L1 intervertebral disc, compared with IG, the maximum stress of SG decreased during flexion, extension, and left/right bending and was basically the same during left/right rotation. However, in the L1-L2 intervertebral discs, the maximum stress of SG increased during left/right rotation compared with that of IG. Furthermore, the maximum displacement of SG was smaller than that of IG in the L1 vertebral bodies under all loading conditions.

Conclusions

SG can reduce the maximum stress in the vertebra and intervertebral discs, offering better biomechanical performance and improved stability than IG.

Biomechanical Analysis of Different Internal Fixation Combined with Different Bone Grafting for Unstable Thoracolumbar Fractures in the Elderly

This research was developed to accurately evaluate the unstable fractures of thoracolumbar before and after surgery and discuss the treatment timing and methods. Three-dimensional (3D) finite element method was adopted to construct the T₁₂-L₅ segment model of human body. The efficiency of percutaneous kyphoplasty (PKP) and percutaneous vertebroplasty (PVP), two commonly used internal fixation procedures, was retrospectively compared. A total of 150 patients with chest fracture who received PKP or PVP surgery in our hospital, and 104 patients with the same symptoms who received conservative treatment were collected and randomly rolled into PVP group (75 cases), PKP group (75 cases), and control group (104 cases). Visual analog scale (VAS) score and Oswestry disability index (ODI) of patients were collected before and after surgery and 2, 12, and 24 months after surgery. Then, the anterior and central height of the patient's cone and the kyphosis angle were calculated by X-ray. Lumbar minimally invasive fusion system and lumbar pedicle screw rod system were established by computer-aided design (CAD), and the biomechanical characteristics were analyzed. The results showed that there was no substantial difference in VAS score and ODI score between PKP and PVP (P > 0.05), but they were higher than those of the control group (P < 0.05). The anterior edge and middle height of vertebra in the two groups were higher than those in control group (P < 0.05), and the increase in PKP group was more substantial (P < 0.05). The kyphosis of the two groups was smaller than that of the control group (P < 0.05), and the decrease of the kyphosis of the PKP group was more substantial (P < 0.05). In summary, the thoracolumbar segment model established by 3D finite element method was an effective model, and it was verified on patients that both PKP and PVP could achieve relatively satisfactory efficacy. The implantation of the new internal fixation system had no obvious effect on the lumbar movement. This work provided a novel idea and method for the treatment of senile thoracolumbar unstable fracture, as well as experimental data of biomechanics for the operation of senile unstable fracture.

FINITE ELEMENT ANALYSIS OF THORACIC VERTEBRAL STABILITY SUPPORTED BY THE FOURTH SPINE

ABSTRACT Objective: In traumatic injuries of the thoracic spine, three variables are analyzed to make decisions: morphology of the injury, posterior ligamentous complex and neurological status; currently the fourth column is not evaluated; our objective was to determine the biomechanical behavior of the spine with a fracture of the fifth thoracic vertebral body when accompanied by a short oblique fracture of the sternum. Methods: An anonymous model of a healthy 25-year-old male was used, from which the thoracic spine and rib cage were obtained; in addition to the ligaments of the posterior complex and the intervertebral discs, four models were simulated. An axial section was made, a load of 400 N was applied, and the biomechanical behavior of each model was determined. Results: The area that suffered the most stress at the vertebral level was the posterior column of T4-T5 (tensile strength of 747 MPa), which exceeded the plastic limit, the load through the ribs was distributed from the first to the sixth (100 MPa), in the sternum the stress increased (200 MPa), the deformity increased to 45 mm. Conclusions: The sternum was a fundamental part of the spine’s stability; the combined injury severely increased the stress (8 MPa to 747 MPa) in the spine and exceeded the plastic limit, which generated an instability that is represented by the global deformity acquired (1 mm to 45 mm). Level of evidence II; Prospective comparative study.

Stress analysis of the thoracolumbar junction in the process of backward fall: An experimental study and finite element analysis

The aim of the present study was to evaluate the biomechanical mechanism of injuries of the thoracolumbar junction by the methods of a backward fall simulation experiment and finite element (FE) analysis (FEA). In the backward fall simulation experiment, one volunteer was selected to obtain the contact force data of the sacrococcygeal region during a fall. Utilizing the fall data, the FEA simulation of the backward fall process was given to the trunk FE model to obtain the stress status of local bone structures of the thoracolumbar junction during the fall process. In the fall simulation test, the sacrococcygeal region of the volunteer landed first; the total impact time was 1.14±0.58 sec, and the impact force was up to 4,056±263 N. The stress of thoracic (T)11 was as high as 42 MPa, that of the posterior margin and the junction of T11 was as high as 70.67 MPa, and that of the inferior articular process and the superior articular process was as high as 128 MPa. The average stress of T12 and the anterior margin of lumbar 1 was 25 MPa, and that of the endplate was as high as 21.7 MPa, which was mostly distributed in the back of the endplate and the surrounding cortex. According to the data obtained from the fall experiment as the loading condition of the FE model, the backward fall process can be simulated to improve the accuracy of FEA results. In the process of backward fall, the front edge of the vertebral body and the root of vertebral arch in the thoracolumbar junction are stress concentration areas, which have a greater risk of injury.

Effect of Model Parameters on the Biomechanical Behavior of the Finite Element Cervical Spine Model

Finite element (FE) models have frequently been used to analyze spine biomechanics. Material parameters assigned to FE spine models are generally uncertain, and their effect on the characterization of the spinal components is not clear. In this study, we aimed to analyze the effect of model parameters on the range of motion, stress, and strain responses of a FE cervical spine model. To do so, we created a computed tomography-based FE model that consisted of C2-C3 vertebrae, intervertebral disc, facet joints, and ligaments. A total of 32 FE analyses were carried out for two different elastic modulus equations and four different bone layer numbers under four different loading conditions. We evaluated the effects of elastic modulus equations and layer number on the biomechanical behavior of the FE spine model by taking the range of angular motion, stress, and strain responses into account. We found that the angular motions of the one- and two-layer models had a greater variation than those in the models with four and eight layers. The angular motions obtained for the four- and eight-layer models were almost the same, indicating that the use of a four-layer model would be sufficient to achieve a stress value converging to a certain level as the number of layers increases. We also observed that the equation proposed by Gupta and Dan (2004) agreed well with the experimental angular motion data. The outcomes of this study are expected to contribute to the determination of the model parameters used in FE spine models.

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.

Comparison of dynamic response of three TLIF techniques on the fused and adjacent segments under vibration

To explore which TLIF techniques are advantageous in reducing the risk of complications and conducive to bone fusion under the vibration. The L1-L5 finite element lumbar model was modified to simulate three different TLIF techniques (a unilateral standard cage, a crescent-shaped cage, and bilateral standard cages). The results showed that the crescent-shaped cage may reduce the risk of subsidence and provide a more stable and suitable environment for vertebral cell growth under the vibration compared to the other TLIF techniques. Unilateral cage may increase the risk of adjacent segment disease and cage failure including fatigue failure under vibration.