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

Interfragmentary compression force and fixation stability of lateral tibial plateau fractures in normal and osteoporotic bones

Lateral platform collapse in fixations of lateral tibial plateau fractures (TPFs) using either double-lag screws fixation (DSF) or locking-plate fixation (LPF) is not rare. This study aimed to explore the effect of enhancing the interfragmentary compression force (IFCF) on fixation stability in lateral TPFs in normal and osteoporotic bones using finite element analysis. Finite element models of DSF in normal bone and LPF in normal and osteoporotic bones were established to simulate the fixations of lateral TPF. After model validation, axial compressive forces of 500, 1000, 1500, and 2500 N to the tibial plateau along with an IFCF of 0, 100, 200, and 300 N were applied. The maximum axial micromotion of the lateral fragment (MAM-LF), maximal translational micromotion of the lateral fragment (MTM-LF), peak von Mises stress (VMS), and peak equivalent elastic strain of the lateral fragment (EES-LF) were evaluated. The MAM-LF showed a decreasing trend as the IFCF increased in all models. For DSF models, the peak VMS of implants increased as the IFCF increased when the axial loads were 500 and 1000 N. The peak EES-LF decreased as the IFCF increased under axial loads of 1000, 1500, and 2500 N. For the normal and osteoporotic LPF models, the peak VMS of the implants decreased as the IFCF increased. Peak EES-LF decreased as IFCF increased. In conclusion, enhancing IFCF was beneficial in improving the fixation stability of lateral TPF. The optimal IFCF for DSF and LPF should be as high as reasonably feasible.

Investigation of the Role of Osteoporotic Vertebra Degeneration on the Stability of the Lumbar Spine: In Silico Modelling under Compressive Loading

An evaluation of the impact of osteoporosis on loss of spinal stability, with or without intervertebral disc degeneration, using computational analysis is presented. The research also investigates the correlation between osteoporosis and intervertebral disc degeneration. Three-dimensional finite element models of human lumbar spine segments were used to assess the influence of osteoporosis on spinal stability. Five different models of age-related degeneration were created using various material properties for trabecular bone and intervertebral discs. Calculation results indicate that in a spine with osteoporosis, the deformation of the intervertebral discs can increase by more than 30% when compared to a healthy spine. Thus, intervertebral disc deformation depends not only on the degree of degeneration of the discs themselves, but their deformation is also influenced by the degree of osteoporosis of the vertebrae. Additionally, the load-bearing capacity of the spine can decrease by up to 30% with osteoporosis, regardless of the degree of intervertebral disc deformation. In conclusion, osteoporosis can contribute to intervertebral disc degeneration.

Validation and Estimation of Obesity-Induced Intervertebral Disc Degeneration through Subject-Specific Finite Element Modelling of Functional Spinal Units

(1) Background: Intervertebral disc degeneration has been linked to obesity; its potential mechanical effects on the intervertebral disc remain unknown. This study aimed to develop and validate a patient-specific model of L3-L4 vertebrae and then use the model to estimate the impact of increasing body weight on disc degeneration. (2) Methods: A three-dimensional model of the functional spinal unit of L3-L4 vertebrae and its components were developed and validated. Validation was achieved by comparing the range of motions (RoM) and intradiscal pressures with the previous literature. Subsequently, the validated model was loaded according to the body mass index and estimated stress, deformation, and RoM to assess disc degeneration. (3) Results: During validation, L3-L4 RoM and intradiscal pressures: flexion 5.17° and 1.04 MPa, extension 1.54° and 0.22 MPa, lateral bending 3.36° and 0.54 MPa, axial rotation 1.14° and 0.52 MPa, respectively. When investigating the impact of weight on disc degeneration, escalating from normal weight to obesity reveals an increased RoM, by 3.44% during flexion, 22.7% during extension, 29.71% during lateral bending, and 33.2% during axial rotation, respectively. Also, stress and disc deformation elevated with increasing weight across all RoM. (4) Conclusions: The predicted mechanical responses of the developed model closely matched the validation dataset. The validated model predicts disc degeneration under increased weight and could lay the foundation for future recommendations aimed at identifying predictors of lower back pain due to disc degeneration.

Finite Element Analysis for Degenerative Cervical Myelopathy: Scoping Review of the Current Findings and Design Approaches, Including Recommendations on the Choice of Material Properties

BACKGROUND

Degenerative cervical myelopathy (DCM) is a slow-motion spinal cord injury caused via chronic mechanical loading by spinal degenerative changes. A range of different degenerative changes can occur. Finite element analysis (FEA) can predict the distribution of mechanical stress and strain on the spinal cord to help understand the implications of any mechanical loading. One of the critical assumptions for FEA is the behavior of each anatomical element under loading (ie, its material properties).

OBJECTIVE

This scoping review aims to undertake a structured process to select the most appropriate material properties for use in DCM FEA. In doing so, it also provides an overview of existing modeling approaches in spinal cord disease and clinical insights into DCM.

METHODS

We conducted a scoping review using qualitative synthesis. Observational studies that discussed the use of FEA models involving the spinal cord in either health or disease (including DCM) were eligible for inclusion in the review. We followed the PRISMA-ScR (Preferred Reporting Items for Systematic Reviews and Meta-Analyses extension for Scoping Reviews) guidelines. The MEDLINE and Embase databases were searched to September 1, 2021. This was supplemented with citation searching to retrieve the literature used to define material properties. Duplicate title and abstract screening and data extraction were performed. The quality of evidence was appraised using the quality assessment tool we developed, adapted from the Newcastle-Ottawa Scale, and shortlisted with respect to DCM material properties, with a final recommendation provided. A qualitative synthesis of the literature is presented according to the Synthesis Without Meta-Analysis reporting guidelines.

RESULTS

A total of 60 papers were included: 41 (68%) "FEA articles" and 19 (32%) "source articles." Most FEA articles (33/41, 80%) modeled the gray matter and white matter separately, with models typically based on tabulated data or, less frequently, a hyperelastic Ogden variant or linear elastic function. Of the 19 source articles, 14 (74%) were identified as describing the material properties of the spinal cord, of which 3 (21%) were considered most relevant to DCM. Of the 41 FEA articles, 15 (37%) focused on DCM, of which 9 (60%) focused on ossification of the posterior longitudinal ligament. Our aggregated results of DCM FEA indicate that spinal cord loading is influenced by the pattern of degenerative changes, with decompression alone (eg, laminectomy) sufficient to address this as opposed to decompression combined with other procedures (eg, laminectomy and fusion).

CONCLUSIONS

FEA is a promising technique for exploring the pathobiology of DCM and informing clinical care. This review describes a structured approach to help future investigators deploy FEA for DCM. However, there are limitations to these recommendations and wider uncertainties. It is likely that these will need to be overcome to support the clinical translation of FEA to DCM.

Understanding the etiopathogenesis of lumbar intervertebral disc herniation: From clinical evidence to basic scientific research

Lumbar intervertebral disc herniation, as a leading cause of low back pain, productivity loss, and disability, is a common musculoskeletal disorder that results in significant socioeconomic burdens. Despite extensive clinical and basic scientific research efforts, herniation etiopathogenesis, particularly its initiation and progression, is not well understood. Understanding herniation etiopathogenesis is essential for developing effective preventive measures and therapeutic interventions. Thus, this review seeks to provide a thorough overview of the advances in herniation-oriented research, with a discussion on ongoing challenges and potential future directions for clinical, translational, and basic scientific investigations to facilitate innovative interdisciplinary research aimed at understanding herniation etiopathogenesis. Specifically, risk factors for herniation are identified and summarized, including familial predisposition, obesity, diabetes mellitus, smoking tobacco, selected cardiovascular diseases, disc degeneration, and occupational risks. Basic scientific experimental and computational research that aims to understand the link between excessive mechanical load, catabolic tissue remodeling due to inflammation or insufficient nutrient supply, and herniation, are also reviewed. Potential future directions to address the current challenges in herniation-oriented research are explored by combining known progressive development in existing research techniques with ongoing technological advances. More research on the relationship between occupational risk factors and herniation, as well as the relationship between degeneration and herniation, is needed to develop preventive measures for working-age individuals. Notably, researchers should explore using or modifying existing degeneration animal models to study herniation etiopathogenesis, as such models may allow for a better understanding of how to prevent mild-to-moderately degenerated discs from herniating.

Multilevel Pedicle Subtraction Osteotomy for Correction of Thoracolumbar Kyphosis in Ankylosing Spondylitis: Clinical Effect and Biomechanical Evaluation

OBJECTIVE

To compare the clinical outcomes and biomechanical characteristics of 1-, 2-, and 3-level pedicle subtraction osteotomy (PSO), and establish selection criteria based on preoperative radiographic parameters.

METHODS

Patients undergone PSO to treat ankylosing spondylitis from February 2009 to May 2019 in Sun Yat-sen Memorial Hospital of Sun Yat-sen University were enrolled. According to the quantity of osteotomy performed, the participants were divided into group A (1-level PSO, n = 24), group B (2-level PSO, n = 19), and group C (3-level PSO, n = 11). Clinical outcomes were assessed before surgery and at the final follow-up. Comparisons of the radiographic parameters and quality-of-life indicators were performed among and within these groups, and the selection criteria were established by regression. Finite element analysis was conducted to compare the biomechanical characteristics of the spine treated with different quantity of osteotomies under different working conditions.

RESULTS

Three-level PSO improved the sagittal parameters more significantly, but resulted in longer operative time and greater blood loss (p < 0.05). Greater stress was found in the proximal screws and proximal junction area of the vertebra in the model simulating 1-level PSO. Larger stress of screws and vertebra was observed at the distal end in the model simulating 3-level PSO.

CONCLUSION

Multilevel PSO works better for larger deformity correction than single-level PSO by allowing greater sagittal parameter correction and obtaining a better distribution of stress in the hardware construct, although with longer operation time and greater blood loss. Three-level osteotomy is recommended for the patients with preoperative of global kyphosis > 85.95°, T1 pelvic angle > 62.3°, sagittal vertical alignment > 299.55 mm, and pelvic tilt+ chin-brow vertical angle > 109.6°.

Impact of osteoporosis and Cement-Augmented fusion on adjacent spinal levels Post-Fusion Surgery: Patient-Specific finite element analysis

Cement-augmentation is a technique commonly used during posterior lumbar instrumented fusion (PLIF) to reinforce compromised osteoporotic vertebral bone, minimize the risk of loosening screws, enhance stability, and improve overall surgical outcomes. In this study, we introduce a novel segmented vertebral body regional modeling approach to investigate the effects of osteoporosis and cement-augmented lumbar fusion on disc biomechanics at spinal levels adjacent to the fused vertebrae. Using our previously validated personalized-poroelastic-osteoligamentous FE model of the spine, fusion was simulated at L4-L5, and the biomechanics of adjacent levels were studied for 30 patients (non-osteoporotic patients (N = 15), osteoporotic patients (N = 15)). PLIF models, with and without cement-augmentation, were developed and compared after an 8 h-rest period (200 N), following a 16 h-cyclic compressive loading of 500-1000 N (40 and 20 min, respectively). Movement in different directions (flexion/ extension/ lateral bending/ axial rotation) was simulated using 10Nm moment before and after cyclic loading. The material mapping algorithm was validated by comparing the results of voxel-based and parametric models. The FE cement-augmented models, subject to daily activity loading, demonstrated significant differences in disc height loss and fluid loss as compared to non-cemented models. The calculated axial stress and fiber strain values were also significantly higher for these models. This work demonstrates that although osteoporosis does not significantly alter the time-dependent characteristics of adjacent IVDs post-surgery, cement-augmentation increases the risk of adjacent segment disease (ASD) incidence. A holistic understanding of the trade-offs and long-term complex interplay between structural reinforcement modalities, including cement augmentation, and altered biomechanics warrants further investigation.

An Enhanced Spine Model Validated for Simulating Dynamic Lifting Tasks in OpenSim

A fully articulated thoracolumbar spine model had been previously developed in OpenSim and had been extensively validated against experimental data during various static tasks. In the present study, we enhanced this detailed musculoskeletal model by adding the role of passive structures and adding kinematic constraints to make it suitable for dynamic tasks. We validated the spinal forces estimated by this enhanced model during nine dynamic lifting/lowering tasks. Moreover, we recently developed and evaluated five approaches in OpenSim to model the external loads applied to the hands during lifting/lowering tasks, and in the present study, we assessed which approach results in more accurate spinal forces. Regardless of the external load modeling approach, the maximum forces predicted by our enhanced spine model across all tasks, as well as the pattern of estimated spinal forces within each task, showed strong correlations (r-values and cross-correlation coefficients > 0.9) with experimental data. Given the biofidelity of our enhanced model, its accessibility via the open-source OpenSim software, and the extent to which this model has been validated, we recommend it for applications requiring estimation of spinal forces during lifting/lowering tasks using multibody-based models and inverse dynamic analyses.

Finite element mechanical analysis of ipsilateral approach and contralateral approach in unilateral bilateral endoscopic spine surgery

BACKGROUND

Unilateral bilateral endoscopic spine surgery (UBE) is often performed to treat lumbar spinal stenosis and disc herniation. It has become a prominent method in endoscopic spine surgery because of its very low learning curve and broader operative field of vision. Currently, the ipsilateral approach and contralateral approach have been established for disc herniation in the foraminal area, intervertebral foramen region, or pedicle region. The contralateral method offers many benefits over the ipsilateral approach, including less bone labour during microsurgical decompression and the preservation of facet joints. However, because it uses the interlaminar window approach, it inevitably involves osteotomy of the patient's superior and inferior articular processes, which may result in corresponding deterioration in the spine's biomechanical stability and subsequent adjacent facet joint diseases caused by facet joint degeneration postoperatively.

OBJECTIVE

As a result, the purpose of this work is to use a finite element model to evaluate how the ipsilateral approach and contralateral approach in unilateral bilateral endoscopic spine surgery affect spinal stability while treating identical intervertebral disc herniation.

STUDY DESIGN

In this study, a three-dimensional lumbar-sacral spine model was built and verified. Osteotomies were conducted for armpit-type lumbar disc herniation (LDH), periradicular-type LDH, and shoulder-type LDH. Postoperative lumbar spine models of the ipsilateral approach and contralateral approach in unilateral bilateral endoscopic spine surgery were developed. The von Mises stress on the endplate, shear force on the annulus fibrosus, pressure inside the intervertebral disc, and range of motion (ROM) of the L3 segment were all determined. The results of our well-validated model showed that osteotomy done in the ipsilateral approach deteriorated most biomechanical metrics.

RESULTS

In the majority of loading conditions, the contralateral approach caused the intervertebral disc's biomechanical properties to increase, and the ipsilateral approach caused the intervertebral disc's biomechanical properties to increase sharply more than the contralateral approach.

CONCLUSION

The contralateral approach, which is now extensively employed in unilateral bilateral endoscopic spine surgery, may be regarded as an ideal surgical alternative for treating lumbar disc herniation without producing iatrogenic instability. This approach has a low facet joint reduction rate, minimum soft tissue injury, and precisely identifies the midline of the central spinal canal during the retraction of the thecal sac and nerve roots.

Enhancing the mechanical stability of restored teeth with interfacial cracks: Finite element analysis

OBJECTIVES

This study aims to enhance the mechanical stability of restored molar teeth with class II occlusal-distal (OD) cavities. We seek to achieve this goal through a comprehensive investigation of three primary factors: (1) the choice of restoration material properties, (2) internal cavity geometries, and (3) the impact of double-layered restoration configurations.

METHODS

To achieve our objectives, we initiated by creating two-dimensional (2D) models of restored teeth featuring class II OD cavities, utilizing scanned and segmented images of maxillary molar teeth. We drew 2D profiles of dentine and enamel, which were then imported into finite element analysis (FEA) software. To explore various cavity geometries, we implemented a total of thirteen different designs, encompassing straight, oblique, grooved, curved, and double-layered configurations. We utilized a semi-circular stone to simulate the application of contact load on the restored tooth. We applied identical boundary conditions and contact loading across all models. To assign material properties, we developed a Python code, enabling the automatic assignment of seven elastic moduli ranging from 2 GPa to 26 GPa to the restoration materials. Meanwhile, constant material properties were assigned to the enamel and dentine. In total, we conducted 133 FEA simulations to comprehensively analyse the effects of the aforementioned factors on the strength and performance of restored molar teeth.

RESULTS

Our analysis revealed two key factors significantly influencing the mechanical resistance of treated teeth, particularly in the presence of a crack or debonding: (1) the marginal geometry of the OD cavity and (2) the elastic modulus of the restorative material. However, altering the internal cavity angle and implementing a double-layered restoration did not significantly influence the restored tooth's overall strength and performance in the face of crack or debonding situations.

SIGNIFICANCE

The findings of this study have substantial implications for designing and restoring class II OD cavities to enhance resistance to cracks or debonding. The use of curved marginal geometries in restoration design can significantly improve fracture resistance, with double-curved geometries reducing stress concentrations by approximately 43% compared to straight cavities. These results offer valuable guidance for strengthening the structural integrity of restored teeth, calling for further experimental investigations to explore practical applications and benefits.

Modelling for design and evaluation of industrial exoskeletons: A systematic review

Industrial exoskeletons are developed to relieve workers' physical demands in the workplace and to alleviate ergonomic issues associated with work-related musculoskeletal disorders. As a safe and economical alternative to empirical/experimental methods, modelling is considered as a powerful tool for design and evaluation of industrial exoskeletons. This systematic review aims to provide a comprehensive understanding of the current literature on the design and evaluation of industrial exoskeletons through modelling. A systematic study was conducted by general keyword searches of five electronic databases over the last two decades (2003-2022). Out of the 701 records initially retrieved, 33 eligible articles were included and analyzed in the final review, presenting a variety of model inputs, model development, and model outputs used in the modelling. This systematic review study revealed that existing modelling methods can evaluate the biomechanical and physiological effects of industrial exoskeletons and provide some design parameters. However, the modelling method is currently unable to cover some of the main evaluation metrics supported by experimental assessments, such as task performance, user experience/discomfort, change in metabolic costs etc. Standard guidelines for model construction and implementation, as well as validation of human-exoskeleton interactions, remain to be established.

Oblique lateral interbody fusion with internal fixations in the treatment for cross-segment degenerative lumbar spine disease (L2-3 and L4-5) finite element analysis

Multi-segmental lumbar degenerative disease, including intersegmental disc degeneration, is found in clinical practice. Controversy still exists regarding the treatment for cross-segment degeneration. Oblique Lateral Interbody Fusion (OLIF) with several internal fixations was used to treat cross-segment lumbar degenerative disease. A whole lumbar spine model was extracted from CT images of the whole lumbar spine of patients with lumbar degeneration. The L2-3 and L4-5 intervertebral spaces were fused with OLIF using modeling software, the Pedicle screws were performed on L2-3 and L4-5, and different internal fixations were performed on L3-4 in Finite Element (FE) software. Among the six 10 Nm moments of different directions, the L3-4 no surgery (NS) group had the relatively largest Range of Motion (ROM) in the whole lumbar spine, while the L2-5 Long segmental fixation (LSF)group had the smallest ROM and the other groups had similar ROM. The ROM in the L1-2 and L5-S1 was relatively close in the six group models, and the articular cartilage stress and disc stress on the L1-2 and L5-S1 were relatively close. In contrast, the L3-4 ROM differed relatively greatly, with the LSF ROM the smallest and the NS ROM the largest, and the L3-4 Coflex (Coflex) group more active than the L3-4 Bacfuse (Bacfuse) group and the L3-4 translaminar facet screw fixation (TFSF) group. The stress on the articular cartilage and disc at L3-4 was relatively greater in the NS disc and articular cartilage, and greater in the Coflex group than in the Bacfuse and TFSF groups, with the greatest stress on the internal fixation in the TFSF group, followed by the Coflex group, and relatively similar stress in the Bacfuse, LSF, and NS groups. In the TFSF group, the stress on the internal fixation was greater than the yield strength among different directional moments of 10 Nm, which means it is unsuitable to be an internal fixation. The LSF group had the greatest overall ROM, which may lead to postoperative low back discomfort. The NS group has the greatest overall ROM, but its increased stress on the L3-4 disc and articular cartilage may lead to accelerated degeneration of the L3-4 disc and articular cartilage. The Coflex and Bacfuse groups had a reduced L3-4 ROM but a greater stress on disc compared to the LSF group, which may lead to disc degeneration in the long term. However, their stress on the articular cartilage was relatively low. Coflex and Bacfuse can still be considered better surgical options.

Dynamic segmental kinematics of the lumbar spine during diagnostic movements

Background: In vivo measurements of segmental-level kinematics are a promising avenue for better understanding the relationship between pain and its underlying, multi-factorial basis. To date, the bulk of the reported segmental-level motion has been restricted to single plane motions. Methods: The present work implemented a novel marker set used with an optical motion capture system to non-invasively measure dynamic, 3D in vivo segmental kinematics of the lower spine in a laboratory setting. Lumbar spinal kinematics were measured for 28 subjects during 17 diagnostic movements. Results: Overall regional range of motion data and lumbar angular velocity measurement were consistent with previously published studies. Key findings from the work included measurement of differences in ascending versus descending segmental velocities during functional movements and observations of motion coupling paradigms in the lumbar spinal segments. Conclusion: The work contributes to the task of establishing a baseline of segmental lumbar movement patterns in an asymptomatic cohort, which serves as a necessary pre-requisite for identifying pathological and symptomatic deviations from the baseline.

Influence of Simulated State of Disc Degeneration and Axial Stiffness of Coupler in a Hybrid Performance Stabilisation System on the Biomechanics of a Spine Segment Model

Spinal fusion surgery leads to the restriction of mobility in the vertebral segments postoperatively, thereby causing stress to rise at the adjacent levels, resulting in early degeneration and a high risk of adjacent vertebral fractures. Thus, to address this issue, non-fusion surgery applies some pedicle screw-based dynamic stabilisation systems to provide stability and micromotion, thereby reducing stress in the fusion segments. Among these systems, the hybrid performance stabilisation system (HPSS) combines a rigid rod, transfer screw, and coupler design to offer a semi-rigid fixation method that preserves some mobility near the fusion site and reduces the adjacent segment compensatory effects. However, further research and confirmation are needed regarding the biomechanical effects of the dynamic coupler stiffness of the HPSS on the intrinsic degenerated adjacent segment. Therefore, this study utilised the finite element method to investigate the impact of the coupler stiffness of the HPSS on the mobility of the lumbar vertebral segments and the stress distribution in the intervertebral discs under flexion, extension, and lateral bending, as well as the clinical applicability of the HPSS on the discs with intrinsic moderate and severe degeneration at the adjacent level. The analytical results indicated that, regardless of the degree of disc degeneration, the use of a dynamic coupler stiffness of 57 N/mm in the HPSS may reduce the stress concentrations at the adjacent levels. However, for severely degenerated discs, the postoperative stress on the adjacent segments with the HPSS was still higher compared with that of the discs with moderate degeneration. We conclude that, when the discs had moderate degeneration, increasing the coupler stiffness led to a decrease in disc mobility. In the case of severe disc degeneration, the effect on disc mobility by coupler stiffness was less pronounced. Increasing the coupler stiffness ked to higher stress on intervertebral discs with moderate degeneration, while its effect on stress was less pronounced for discs with severe degeneration. It is recommended that patients with severe degeneration who undergo spinal dynamic stabilisation should remain mindful of the risk of accelerated adjacent segment degeneration.

Finite element study of the biomechanical effects on the rotator cuff under load

Rotator cuff injuries account for 50% of shoulder disorders that can cause shoulder pain and reduced mobility. The occurrence of rotator cuff injury is related to the variation in shoulder load, but the mechanical changes in the rotator cuff caused by load remain unclear. Therefore, the mechanical results of the rotator cuff tissue during glenohumeral abduction and adduction were analyzed based on a finite element shoulder model under non-load (0 kg) and load (7.5 kg) conditions. The results showed that the maximum von Mises stress on the supraspinatus muscle was larger than that on the subscapularis, infraspinatus, and teres minor muscles during glenohumeral abduction. Compared with the non-load condition, the maximum von Mises stress on the supraspinatus muscle increased by 75% under the load condition at 30° abduction. Under the load condition, the supraspinatus joint side exhibited an average stress that was 32% greater than that observed on the bursal side. The von Mises stress on the infraspinatus muscle was higher than that in other rotator cuff tissues during adduction. The stress on the infraspinatus muscle increased by 36% in the load condition compared to the non-load condition at 30° adduction. In summary, the increased load changed the mechanical distribution of rotator cuff tissue and increased the stress differential between the joint aspect and the bursal aspect of the supraspinatus tendon.

Biomechanical effect of proximal multifidus injury on adjacent segments during posterior lumbar interbody fusion: a finite element study

BACKGROUND

Adjacent segment degeneration (ASD) is a common complication of lumbar interbody fusion; the paraspinal muscles significantly maintain spinal biomechanical stability. This study aims to investigate the biomechanical effects of proximal multifidus injury on adjacent segments during posterior lumbar interbody fusion (PLIF).

METHODS

Data from a lumbosacral vertebral computed tomography scan of a healthy adult male volunteer were used to establish a normal lumbosacral vertebral finite element model and load the muscle force of the multifidus. A normal model, an L4/5 PLIF model (PFM) based on a preserved proximal multifidus, a total laminectomy PLIF model (TLPFM), and a hemi-laminectomy PLIF model based on a severed proximal multifidus were established, respectively. The range of motion (ROM) and maximum von Mises stress of the upper and lower adjacent segments were analyzed along with the total work of the multifidus muscle force.

RESULTS

This model verified that the ROMs of all segments with four degrees of freedom were similar to those obtained in previous research data, which validated the model. PLIF resulted in an increased ROM and maximum von Mises stress in the upper and lower adjacent segments. The ROM and maximum von Mises stress in the TLPFM were most evident in the upper adjacent segment, except for lateral bending. The ROM of the lower adjacent segment increased most significantly in the PFM in flexion and extension and increased most significantly in the TLPFM in lateral bending and axial rotation, whereas the maximum von Mises stress of the lower adjacent segment increased the most in the TLPFM, except in flexion. The muscle force and work of the multifidus were the greatest in the TLPFM.

CONCLUSIONS

PLIF increased the ROM and maximum von Mises stress in adjacent cranial segments. The preservation of the proximal multifidus muscle contributes to the maintenance of the physiological mechanical behavior of adjacent segments, thus preventing the occurrence and development of ASD.

Biomechanical effect of endplate defects on the intermediate vertebral bone in consecutive two-level anterior cervical discectomy and fusion: a finite element analysis

BACKGROUND

Intermediate vertebral collapse is a newly discovered complication of consecutive two-level anterior cervical discectomy and fusion (ACDF). There have been no analytical studies related to the effects of endplate defects on the biomechanics of the intermediate vertebral bone after ACDF. This study aimed to compare the effects of endplate defects on the intermediate vertebral bone biomechanics in the zero-profile (ZP) and cage-and-plate (CP) methods of consecutive 2-level ACDF and to determine whether collapse of the intermediate vertebra is more likely to occur using ZP.

METHODS

A three-dimensional finite element (FE) model of the intact cervical spine (C2-T1) was constructed and validated. The intact FE model was then modified to build ACDF models and imitate the situation of endplate injury, establishing two groups of models (ZP, IM-ZP and CP, IM-ZP). We simulated cervical motion, such as flexion, extension, lateral bending and axial rotation, and compared the range of motion (ROM), upper and lower endplate stress, fusion fixation device stress, C5 vertebral body stress, intervertebral disc internal pressure (intradiscal pressure, or IDP) and the ROM of adjacent segments in the models.

RESULTS

There was no significant difference between the IM-CP model and the CP model in the ROM of the surgical segment, upper and lower endplate stress, fusion fixation device stress, C5 vertebral body stress, IDP, or ROM of the adjacent segments. Compared with the CP model, the endplate stress of the ZP model is significantly higher in the flexion, extension, lateral bending and axial rotation conditions. Compared with the ZP model, endplate stress, screw stress, C5 vertebral stress and IDP in IM-ZP were significantly increased under flexion, extension, lateral bending and axial rotation conditions.

CONCLUSIONS

Compared to consecutive 2-level ACDF using CP, collapse of the intermediate vertebra is more likely to occur using ZP due to its mechanical characteristics. Intraoperative endplate defects of the anterior lower margin of the middle vertebra are a risk factor leading to collapse of the middle vertebra after consecutive 2-level ACDF using ZP.

Biomechanical evaluation of the thoracolumbar spine comparing healthy and irregular thoracic and lumbar curvatures

BACKGROUND

The geometric alignment of the spine plays an integral role in stability, biomechanical loading, and consequently, pain, and a range of healthy sagittal curvatures has been identified. Spinal biomechanics when sagittal curvature is outside the optimal range remains a debate and may provide insight into the load distribution throughout the spinal column.

METHOD

A thoracolumbar spine model (Healthy) was developed. Thoracic and lumbar curvatures were adjusted by 50% to create models with varying sagittal profiles: hypolordotic (HypoL), hyperlordotic (HyperL), hypokyphotic (HypoK), and hyperkyphotic (HyperK). In addition, lumbar spine models were constructed for the former three profiles. The models were subjected to loading conditions simulating flexion and extension. Following validation, intervertebral disc stresses, vertebral body stresses, disc heights, and intersegmental rotations were compared across all models.

RESULTS

Overall trends demonstrated that HyperL and HyperK models had a noticeable reduction in disc height and greater vertebral body stresses compared to the Healthy model. In comparison, the HypoL and HypoK models displayed opposite trends. Considering the lumbar models, the HypoL model had reduced disc stresses and flexibility, while the contrary was observed in the HyperL model. Results indicate that the models with excessive curvature may be subjected to greater stress magnitudes, while the straighter spine models may reduce these stresses.

CONCLUSIONS

Finite element modeling of spine biomechanics demonstrated that variations in sagittal profiles influence the load distribution and range of motion of the spine. Considering patient-specific sagittal profiles in finite element modeling may provide valuable insight for biomechanical analyses and targeted treatments.

Biomechanical performance of the novel assembled uncovertebral joint fusion cage in single-level anterior cervical discectomy and fusion: A finite element analysis

Introduction: Anterior cervical discectomy and fusion (ACDF) is widely accepted as the gold standard surgical procedure for treating cervical radiculopathy and myelopathy. However, there is concern about the low fusion rate in the early period after ACDF surgery using the Zero-P fusion cage. We creatively designed an assembled uncoupled joint fusion device to improve the fusion rate and solve the implantation difficulties. This study aimed to assess the biomechanical performance of the assembled uncovertebral joint fusion cage in single-level ACDF and compare it with the Zero-P device.

Methods: A three-dimensional finite element (FE) of a healthy cervical spine (C2−C7) was constructed and validated. In the one-level surgery model, either an assembled uncovertebral joint fusion cage or a zero-profile device was implanted at the C5–C6 segment of the model. A pure moment of 1.0 Nm combined with a follower load of 75 N was imposed at C2 to determine flexion, extension, lateral bending, and axial rotation. The segmental range of motion (ROM), facet contact force (FCF), maximum intradiscal pressure (IDP), and screw−bone stress were determined and compared with those of the zero-profile device.

Results: The results showed that the ROMs of the fused levels in both models were nearly zero, while the motions of the unfused segments were unevenly increased. The FCF at adjacent segments in the assembled uncovertebral joint fusion cage group was less than that that of the Zero-P group. The IDP at the adjacent segments and screw–bone stress were slightly higher in the assembled uncovertebral joint fusion cage group than in those of the Zero-P group. Stress on the cage was mainly concentrated on both sides of the wings, reaching 13.4–20.4 Mpa in the assembled uncovertebral joint fusion cage group.

Conclusion: The assembled uncovertebral joint fusion cage provided strong immobilization, similar to the Zero-P device. When compared with the Zero-P group, the assembled uncovertebral joint fusion cage achieved similar resultant values regarding FCF, IDP, and screw–bone stress. Moreover, the assembled uncovertebral joint fusion cage effectively achieved early bone formation and fusion, probably due to proper stress distributions in the wings of both sides.

Biomechanical analysis of different fixation methods for Rorabeck II supracondylar femoral fractures after total knee arthroplasty

BACKGROUND

Locking plate (LP) and retrograde intramedullary nailing (RIMN) are widely used to fix Rorabeck II supracondylar femoral fractures after total knee arthroplasty (TKA). The biomechanical properties of the implant used for treatment influence its longevity. Therefore, we aimed to evaluate the biomechanical stability of different fixations using finite element analysis.

METHODS

Seven finite element models (FEMs) were established, including LP groups (short LP, long LP, and double LP), RIMN groups (short RIMN and long RIMN), and mixed groups (long LP with short RIMN and long LP with long RIMN). The stress of the implants around the fracture area was calculated to evaluate the biomechanical stability under loads.

RESULTS

Stress was mainly distributed around the fracture area in all models. The stress-shielding phenomenon was most evident in the short LP. The trend in maximum equivalent stress values of implants around the fracture area for the seven internal fixations was: short LP (324.63 MPa) > short RIMN (306.37 MPa) > long LP (275.06 MPa) > long RIMN (262.74 MPa) > double LP (203.19 MPa) > long LP with short RIMN (124.42 MPa) > long LP with long RIMN (112.41 MPa). We found that the double LP can better disperse the stress than a single LP, and a long LP with long RIMN can prevent stress concentration and make the stress distribution more uniform.

CONCLUSION

From the perspective of biomechanics, long LP with long RIMN can stabilize fractures and avoid stress concentration in Rorabeck II supracondylar femoral fractures after TKA.

A Review of Image-Based Simulation Applications in High-Value Manufacturing

Image-Based Simulation (IBSim) is the process by which a digital representation of a real geometry is generated from image data for the purpose of performing a simulation with greater accuracy than with idealised Computer Aided Design (CAD) based simulations. Whilst IBSim originates in the biomedical field, the wider adoption of imaging for non-destructive testing and evaluation (NDT/NDE) within the High-Value Manufacturing (HVM) sector has allowed wider use of IBSim in recent years. IBSim is invaluable in scenarios where there exists a non-negligible variation between the 'as designed' and 'as manufactured' state of parts. It has also been used for characterisation of geometries too complex to accurately draw with CAD. IBSim simulations are unique to the geometry being imaged, therefore it is possible to perform part-specific virtual testing within batches of manufactured parts. This novel review presents the applications of IBSim within HVM, whereby HVM is the value provided by a manufactured part (or conversely the potential cost should the part fail) rather than the actual cost of manufacturing the part itself. Examples include fibre and aggregate composite materials, additive manufacturing, foams, and interface bonding such as welding. This review is divided into the following sections: Material Characterisation; Characterisation of Manufacturing Techniques; Impact of Deviations from Idealised Design Geometry on Product Design and Performance; Customisation and Personalisation of Products; IBSim in Biomimicry. Finally, conclusions are drawn, and observations made on future trends based on the current state of the literature.

A Critical Comparison of Comparators Used to Demonstrate Credibility of Physics-Based Numerical Spine Models

The ability of new medical devices and technology to demonstrate safety and effectiveness, and consequently acquire regulatory approval, has been dependent on benchtop, in vitro, and in vivo evidence and experimentation. Regulatory agencies have recently begun accepting computational models and simulations as credible evidence for virtual clinical trials and medical device development. However, it is crucial that any computational model undergo rigorous verification and validation activities to attain credibility for its context of use before it can be accepted for regulatory submission. Several recently published numerical models of the human spine were considered for their implementation of various comparators as a means of model validation. The comparators used in each published model were examined and classified as either an engineering or natural comparator. Further, a method of scoring the comparators was developed based on guidelines from ASME V&V40 and the draft guidance from the US FDA, and used to evaluate the pertinence of each comparator in model validation. Thus, this review article aimed to score the various comparators used to validate numerical models of the spine in order to examine the comparator's ability to lend credibility towards computational models of the spine for specific contexts of use.

Application of vertebral body compression osteotomy in pedicle subtraction osteotomy on ankylosing spondylitis kyphosis: Finite element analysis and retrospective study

BACKGROUND

Ankylosing spondylitis (AS) is a chronic inflammatory rheumatic disease, with pathological characteristics of bone erosion, inflammation of attachment point, and bone ankylosis. Due to the ossified intervertebral disc and ligament, pedicle subtraction osteotomy (PSO) is one of the mainstream surgeries of AS-related thoracolumbar kyphosis, but the large amount of blood loss and high risk of instrumental instability limit its clinical application. The purpose of our study is to propose a new transpedicular vertebral body compression osteotomy (VBCO) in PSO to reduce blood loss and improve stability.

METHODS

A retrospective analysis was performed on patients with AS-related thoracolumbar kyphosis who underwent one-level PSO in our hospital from February 2009 to May 2019. A total of 31 patients were included in this study; 6 received VBCO and 25 received eggshell vertebral body osteotomy. We collected demographic data containing gender and age at diagnosis. Surgical data contained operation time, estimated blood loss (EBL), and complications. Radiographic data contained pre-operative and follow-up sagittal parameters including chin brow-vertical angle (CBVA), global kyphosis (GK), thoracic kyphosis (TK), and lumbar lordosis (LL). A typical case with L2-PSO was used to establish a finite element model. The mechanical characteristics of the internal fixation device, vertebral body, and osteotomy plane of the two osteotomy models were analyzed under different working conditions.

RESULTS

The VBCO could provide comparable restoring of CBVA, GK, TK, and LL in the eggshell osteotomy procedure (all p > 0.05). The VBCO significantly reduced EBL compared to those with eggshell osteotomy [800.0 ml (500.0-1,439.5 ml) vs. 1,455.5 ml (1,410.5-1,497.8 ml), p = 0.033]. Compared with the eggshell osteotomy, VBCO showed better mechanical property. For the intra-pedicular screw fixation, the VBCO group had a more average distributed and lower stress condition on both nails and connecting rod. VBCO had a flattened osteotomy plane than the pitted osteotomy plane of the eggshell group, showing a lower and more average distributed maximum stress and displacement of osteotomy plane.

CONCLUSION

In our study, we introduced VBCO as an improved method in PSO, with advantages in reducing blood loss and providing greater stability. Further investigation should focus on clinical research and biomechanical analysis for the application of VBCO.

Analysis of stress and stabilization in adolescent with osteoporotic idiopathic scoliosis: finite element method

Objective: To explore the effect of osteoporosis on the stress, stability, and lumbar intervertebral disc of AIS lumbar vertebrae by finite element method. Better understand the biomechanical characteristics of osteoporotic scoliosis.Methods: Based on the CT images of normal lumbar vertebrae and lumbar vertebrae with AIS, the finite element models were established to simulate the estimated osteoporosis by changing the Young's modulus of cortical bone, cancellous bone, and endplate. Four finite element models of normal lumbar, osteoporotic lumbar, normal AIS lumbar and osteoporotic AIS lumbar were established, and the same load and boundary conditions were applied respectively. The displacement, stress, and intervertebral disc strain of the four models were compared to explore the effect of osteoporosis on the stability and injury risk of AIS.Results: After suffering from osteoporosis, under the same load, the displacement of lumbar spine increases, the stability decreases, and the stability of AIS lumbar spine decrease more obviously, especially under extension load. Suffering from osteoporosis will increase the stress of lumbar spine, AIS lumbar spine increases more obviously, and the stress is more concentrated, Osteoporotic lumbar spine only affects the strain of intervertebral disc when AIS lumbar spine bends on the concave side, resulting in greater strain behind the concave side of intervertebral disc.Conclusions: AIS patients with OP have lower lumbar stability, a higher risk of fracture of lumbar vertebrae, and spinal nerves are more likely to be compressed by intervertebral discs. OP can aggravate the scoliosis of lumbar vertebrae.

Computational lumbar spine models: A literature review

BACKGROUND

Computational spine models of various types have been employed to understand spine function, assess the risk that different activities pose to the spine, and evaluate techniques to prevent injury. The areas in which these models are applied has expanded greatly, potentially beyond the appropriate scope of each, given their capabilities. A comprehensive understanding of the components of these models provides insight into their current capabilities and limitations.

METHODS

The objective of this review was to provide a critical assessment of the different characteristics of model elements employed across the spectrum of lumbar spine modeling and in newer combined methodologies to help better evaluate existing studies and delineate areas for future research and refinement.

FINDINGS

A total of 155 studies met selection criteria and were included in this review. Most current studies use either highly detailed Finite Element models or simpler Musculoskeletal models driven with in vivo data. Many models feature significant geometric or loading simplifications that limit their realism and validity. Frequently, studies only create a single model and thus can't account for the impact of subject variability. The lack of model representation for certain subject cohorts leaves significant gaps in spine knowledge. Combining features from both types of modeling could result in more accurate and predictive models.

INTERPRETATION

Development of integrated models combining elements from different model types in a framework that enables the evaluation of larger populations of subjects could address existing voids and enable more realistic representation of the biomechanics of the lumbar spine.

Yeditepe spine mesh: Finite element modeling and validation of a parametric CAD model of lumbar spine

Finite element analysis is a powerful tool that is often used to study the biomechanical response of the spine. The primary objective of this study was to illustrate the mechanical behavior of a previously proposed parametric CAD spine model in comparison with a segmented FSU model and the literature. In this study, two finite element models of the L4-L5 spinal level were developed from the same patient's CT scan data. The first was developed using well-known segmentation methods, whereas the second was developed from the new by using a novel parametric CAD model. Both models were subjected to the same loading and boundary conditions to perform flexion, extension, lateral bending and axial rotation motions. The segmented finite element model was observed to be in good agreement with the literature. The parametric finite element model results were also observed to be in good agreement with the segmented finite element model and with the literature except under extension.

The Effect of Concave-Side Intertransverse Ligament Laxity on the Stress of AIS Lumbar Spine Based on Finite Element Method

(1) Background: Scoliosis has the mechanical characteristic of asymmetric stress distribution, which is one of the reasons for the aggravation of scoliosis. Bracing therapy is the best treatment for AIS, but it is difficult and costly to operate. Is it possible to reduce pressure in the concave side by relaxing the ITL in the concave side of scoliosis, so as to improve the abnormal stress distribution of scoliosis? In this paper, a finite element method was used to simulate the effect of the relaxation of concave-side ITL on the stress of a lumbar spine with scoliosis, which provides some guidance for the treatment of scoliosis. (2) Methods: Using CT images of a patient with scoliosis whose Cobb Angle was 43° and Lordosis Angle was 45, a scoliosis lumbar was established, and Young's modulus of the ITL of the concave-side lumbar spine was reduced by 95% to simulate ligament relaxation. By comparing the stress condition of the model vertebral body with no ligament relaxation, the effect of concave-side ITL relaxation on the mechanical characteristics of scoliosis lumbar spine was explored. (3) Results: An effective and complete model of the lumbar spine was established. The concave ITL relaxed, which only had a great impact on the bending loads. After the ligament was relaxed, the stability of the spine was reduced. Stress concentration on the concave side of vertebrae and the IVD was aggravated. Under loads on the convex side, the maximum stress on the vertebral body and the IVD increased significantly, making lumbar vertebrae more vulnerable to injury. (4) Conclusions: Laxity of the ITL on the concave side of the AIS lumbar only affects the bending load. Laxity of the concave-side ligament will reduce the stability of the lumbar, aggravate the uneven stress distribution of scoliotic lumbar vertebrae, increase the risk of IVD injury, and be unfavorable for the scoliotic lumbar spine. Relaxation of the concave ITL alone is not an effective way to treat scoliosis.

Impact of cage position on biomechanical performance of stand-alone lateral lumbar interbody fusion: a finite element analysis

Background

This study aimed to compare the biomechanical performance of various cage positions in stand-alone lateral lumbar interbody fusion(SA LLIF).

Methods

An intact finite element model of the L3-L5 was reconstructed. The model was verified and analyzed. Through changing the position of the cage, SA LLIF was established in four directions: anterior placement(AP), middle placement(MP), posterior placement(PP), oblique placement(OP). A 400 N vertical axial pre-load was imposed on the superior surface of L3 and a 10 N/m moment was applied on the L3 superior surface along the radial direction to simulate movements of flexion, extension, lateral bending, and axial rotation. Various biomechanical parameters were evaluated for intact and implanted models in all loading conditions, including the range of motion (ROM) and maximum stress.

Results

In the SA LLIF models, the ROM of L4-5 was reduced by 84.21–89.03% in flexion, 72.64–82.26% in extension, 92.5-95.85% in right and left lateral bending, and 87.22–92.77% in right and left axial rotation, respectively. Meanwhile, ROM of L3-4 was mildly increased by an average of 9.6% in all motion directions. Almost all stress peaks were increased after SA LLIF, including adjacent disc, facet joints, and endplates. MP had lower stress peaks of cage and endplates in most motion modes. In terms of the stress on facet joints and disc of the cephalad segment, MP had the smallest increment.

Conclusion

In our study, SA LLIF risked accelerating the adjacent segment degeneration. The cage position had an influence on the distribution of endplate stress and the magnitude of facet joint stress. Compared with other positions, MP had the slightest effect on the stress in the adjacent facet joints. Meanwhile, MP seems to play an important role in reducing the risk of cage subsidence.

Improve biomechanical stability using intramedullary nails with femoral neck protection in femoral shaft fractures

BACKGROUND AND OBJECTIVE

Elderly patients treated for femoral shaft fractures have a higher risk of hip fracture. We hypothesized that intramedullary nails protecting the femoral neck can improve mechanical strength and reduce the risk of subsequent hip fracture. This study aims to analyze the biomechanical stability using intramedullary nails with or without femoral neck protection through finite element analysis.

METHODS

Thirty finite element models (FEMs) were established, including five different conditions of femoral shaft fracture: Fracture healing, Proximal fractures (Transverse and oblique), Distal fractures (Transverse and oblique), and five different fixation methods. Femoral neck protection groups: cephalomedullary nail (CN), reconstruction nail (RN); No femoral neck protection groups: type-1 of antegrade intramedullary nail (AIN-1), type-2 of antegrade intramedullary nail (AIN-2), and retrograde intramedullary nail (RIN). The maximum stress of bone and internal fixation in the femoral neck region for all type of fixation were calculated to evaluate the biomechanical stability.

RESULTS

Maximum equivalent stress values of bone in the femoral neck region for five different conditions of femoral shaft fracture: AIN-2 (77.23 MPa) >RIN (77.15 MPa) > AIN-1 (76.71 MPa) > CN (60.74 MPa) > RN (57.66 MPa) for the fracture healing; RIN (80.05 MPa) > AIN-1 (79.15 MPa) > AIN-2(78.77 MPa) > RN (65.16 MPa) > CN (65.03 MPa) for the proximal transverse fracture; RIN (80.10 MPa) > AIN-2 (79.36 MPa) > AIN-1 (79.18 MPa) > RN (65.09 MPa) > CN (64.96 MPa) for the proximal oblique fracture; RIN (80.24 MPa) > AIN-2 (79.68 MPa) > AIN-1 (79.33 MPa) > CN (65.02 MPa) > RN (64.76 MPa) for the distal transverse fracture; RIN (80.23 MPa) > AIN-2 (79.61 MPa) > AIN-1 (79.35 MPa) > CN (65.06 MPa) > RN (64.76 MPa) for the distal oblique fracture. Maximum equivalent stress of internal fixation in the femoral neck region is greater than the maximum stress of bone and avoids stress concentration of bone for the femoral neck protection groups (CN and RN).

CONCLUSIONS

Intramedullary nails with femoral neck protection in the treatment of femoral shaft fractures improve mechanical strength and prevent secondary hip fractures and decrease the overall risk of reoperation postoperatively.

Development of a Computational Model of the Mechanical Behavior of the L4-L5 Lumbar Spine: Application to Disc Degeneration

Degenerative changes in the lumbar spine significantly reduce the quality of life of people. In order to fully understand the biomechanics of the affected spine, it is crucial to consider the biomechanical alterations caused by degeneration of the intervertebral disc (IVD). Therefore, this study is aimed at the development of a discrete element model of the mechanical behavior of the L4-L5 spinal motion segment, which covers all the degeneration grades from healthy IVD to its severe degeneration, and numerical study of the influence of the IVD degeneration on stress state and biomechanics of the spine. In order to analyze the effects of IVD degeneration on spine biomechanics, we simulated physiological loading conditions using compressive forces. The results of modeling showed that at the initial stages of degenerative changes, an increase in the amplitude and area of maximum compressive stresses in the disc is observed. At the late stages of disc degradation, a decrease in the value of intradiscal pressure and a shift in the maximum compressive stresses in the dorsal direction is observed. Such an influence of the degradation of the geometric and mechanical parameters of the tissues of the disc leads to the effect of bulging, which in turn leads to the formation of an intervertebral hernia.

Immuno-Modulatory Effects of Intervertebral Disc Cells

Low back pain is a highly prevalent, chronic, and costly medical condition predominantly triggered by intervertebral disc degeneration (IDD). IDD is often caused by structural and biochemical changes in intervertebral discs (IVD) that prompt a pathologic shift from an anabolic to catabolic state, affecting extracellular matrix (ECM) production, enzyme generation, cytokine and chemokine production, neurotrophic and angiogenic factor production. The IVD is an immune-privileged organ. However, during degeneration immune cells and inflammatory factors can infiltrate through defects in the cartilage endplate and annulus fibrosus fissures, further accelerating the catabolic environment. Remarkably, though, catabolic ECM disruption also occurs in the absence of immune cell infiltration, largely due to native disc cell production of catabolic enzymes and cytokines. An unbalanced metabolism could be induced by many different factors, including a harsh microenvironment, biomechanical cues, genetics, and infection. The complex, multifactorial nature of IDD brings the challenge of identifying key factors which initiate the degenerative cascade, eventually leading to back pain. These factors are often investigated through methods including animal models, 3D cell culture, bioreactors, and computational models. However, the crosstalk between the IVD, immune system, and shifted metabolism is frequently misconstrued, often with the assumption that the presence of cytokines and chemokines is synonymous to inflammation or an immune response, which is not true for the intact disc. Therefore, this review will tackle immunomodulatory and IVD cell roles in IDD, clarifying the differences between cellular involvements and implications for therapeutic development and assessing models used to explore inflammatory or catabolic IVD environments.

The Influence of Mattress Stiffness on Spinal Curvature and Intervertebral Disc Stress-An Experimental and Computational Study

Sleeping support systems can influence spinal curvature, and the misalignment of the spinal curvature can lead to musculoskeletal problems. Previous sleep studies on craniocervical support focused on pillow variants, but the mattress supporting the pillow has rarely been considered. This study used a cervical pillow and three mattresses of different stiffnesses, namely soft, medium, and hard, with an indentation load deflection of 20, 42, and 120 lbs, respectively. A novel electronic curvature measurement device was adopted to measure the spinal curvature, whereby the intervertebral disc loading was computed using the finite element method. Compared with the medium mattress, the head distance increased by 30.5 ± 15.9 mm, the cervical lordosis distance increased by 26.7 ± 14.9 mm, and intervertebral disc peak loading increased by 49% in the soft mattress environment. Considering that the pillow support may increase when using a soft mattress, a softer or thinner pillow is recommended. The head distance and cervical lordosis distance in the hard mattress environment were close to the medium mattress, but the lumbar lordosis distance reduced by 10.6 ± 6.8 mm. However, no significant increase in intervertebral disc loading was observed, but contact pressure increased significantly, which could cause discomfort and health problems.

Biomechanical Effects of a Novel Anatomic Titanium Mesh Cage for Single-Level Anterior Cervical Corpectomy and Fusion: A Finite Element Analysis

Background: The traditional titanium mesh cage (TTMC) has become common as a classical instrument for Anterior Cervical Corpectomy and Fusion (ACCF), but a series of complications such as cage subsidence, adjacent segment degeneration (ASD), and implant-related complications by using the TTMC have often been reported in the previous literature. The aim of this study was to assess whether a novel anatomic titanium mesh cage (NTMC) could improve the biomechanical condition after surgery.

Methods: The NTMC model consists of two spacers located on both sides of the TTMC which match the anatomic structure between the endplates by measuring patient preoperative cervical computed tomography (CT) data. The ranges of motion (ROMs) of the surgical segments and the stress peaks in the C6 superior endplates, titanium mesh cage (TMC), screw–bone interface, anterior titanium plate, and adjacent intervertebral disc were compared.

Results: Compared with the TTMC, the NTMC reduced the surgical segmental ROMs by 89.4% postoperatively. The C6 superior endplate stress peaks were higher in the TTMC (4.473–23.890 MPa), followed by the NTMC (1.923–5.035 MPa). The stress peaks on the TMC were higher in the TTMC (47.896–349.525 MPa), and the stress peaks on the TMC were lower in the NTMC (17.907–92.799 MPa). TTMC induced higher stress peaks in the screw–bone interface (40.0–153.2 MPa), followed by the NTMC (14.8–67.8 MPa). About the stress peaks on the anterior titanium plate, the stress of TTMC is from 16.499 to 58.432 MPa, and that of the NTMC is from 12.456 to 34.607 MPa. Moreover, the TTMC induced higher stress peaks in the C3/4 and C6/7 intervertebral disc (0.201–6.691 MPa and 0.248–4.735 MPa, respectively), followed by the NTMC (0.227–3.690 MPa and 0.174–3.521 MPa, respectively).

Conclusion: First, the application of the NTMC can effectively decrease the risks of TMC subsidence after surgery. Second, in the NTMC, the stresses at the anterior screw-plate, bone–screw, and TMC interface are much less than in the TTMC, which decreased the risks of instrument-related complications after surgery. Finally, increases in IDP at adjacent levels are associated with the internal stresses of adjacent discs which may lead to ASD; therefore, the NTMC can effectively decrease the risks of ASD.

Impact of Femoral Neck Cortical Bone Defect Induced by Core Decompression on Postoperative Stability: A Finite Element Analysis

Objective

To analyze the impact of femoral neck cortical bone defect induced by core decompression on postoperative biomechanical stability using the finite element method.

Methods

Five finite element models (FEMs) were established, including the standard operating model and four models of cortical bone defects at different portions of the femoral neck (anterior, posterior, superior, and inferior). The maximum stress of the proximal femur was evaluated during normal walking and walking downstairs.

Results

Under both weight-bearing conditions, the maximum stress values of the five models were as follows: femoral neck (inferior) > femoral neck (superior) > femoral neck (posterior) > femoral neck (anterior) > standard operation. Stress concentration occurred in the areas of femoral neck cortical bone defect. Under normal walking, the maximum stress of four bone defect models and its increased percentage comparing the standard operation were as follows: anterior (17.17%), posterior (39.02%), superior (57.48%), and inferior (76.42%). The maximum stress was less than the cortical bone yield strength under normal walking conditions. Under walking downstairs, the maximum stress of four bone defect models and its increased percentage comparing the standard operation under normal walking were as follows: anterior (36.75%), posterior (67.82%), superior (83.31%), and inferior (103.65%). Under walking downstairs conditions, the maximum stress of bone defect models (anterior, posterior, and superior) was less than the yield strength of cortical bone, while the maximum stress of bone defect model (inferior) excessed yield strength value.

Conclusions

The femoral neck cortical bone defect induced by core decompression can carry out normal walking after surgery. To avoid an increased risk of fracture after surgery, walking downstairs should be avoided when the cortical bone defect is inferior to the femoral neck except for the other three positions (anterior, posterior, and superior).

Oblique lateral interbody fusion combined with different internal fixations for the treatment of degenerative lumbar spine disease: a finite element analysis

Background

Little is known about the biomechanical performance of different internal fixations in oblique lumbar interbody fusion (OLIF). Here, finite element (FE) analysis was used to describe the biomechanics of various internal fixations and compare and explore the stability of each fixation.

Methods

CT scans of a patient with lumbar degenerative disease were performed, and the l3-S1 model was constructed using relevant software. The other five FE models were constructed by simulating the model operation and adding different related implants, including (1) an intact model, (2) a stand-alone (SA) model with no instrument, (3) a unilateral pedicle screw model (UPS), (4) a unilateral pedicle screw contralateral translaminar facet screw model (UPS-CTFS), (5) a bilateral pedicle screw (BPS) model, and (6) a cortical bone trajectory screw model (CBT). Various motion loads were set by FE software to simulate lumbar vertebral activity. The software was also used to extract the range of motion (ROM) of the surgical segment, CAGE and fixation stress in the different models.

Results

The SA group had the greatest ROM and CAGE stress. The ROM of the BPS and UPS-CTFS was not significantly different among motion loadings. Compared with the other three models, the BPS model had lower internal fixation stress among loading conditions, and the CBT screw internal fixation had the highest stress among loads.

Conclusions

The BPS model provided the best biomechanical stability for OLIF. The SA model was relatively less stable. The UPS-CTFS group had reduced ROM in the fusion segments, but the stresses on the internal fixation and CAGE were relatively higher in the than in the BPS group; the CBT group had a lower flexion and extension ROM and higher rotation and lateral flexion ROM than the BPS group. The stability of the CBT group was poorer than that of the BPS and LPS-CTFS groups. The CAGE and internal fixation stress was greater in the CBT group.

Comparative Analysis of the Biomechanical Characteristics After Different Minimally Invasive Surgeries for Cervical Spondylopathy: A Finite Element Analysis

Minimally invasive surgeries, including posterior endoscopic cervical foraminotomy (PECF), microsurgical anterior cervical foraminotomy (MACF), anterior transdiscal approach of endoscopic cervical discectomy (ATd-ECD), and anterior transcorporeal approach of endoscopic cervical discectomy (ATc-ECD), have obtained positive results for cervical spondylotic radiculopathy. Nonetheless, there is a lack of comparison among them regarding their biomechanical performance. The purpose of this study is to investigate the biomechanical changes of operated and adjacent segments after minimally invasive surgeries compared to a normal cervical spine. A three-dimensional model of normal cervical vertebrae C3–C7 was established using finite element analysis. Afterwards, four surgical models (PECF, MACF, ATd-ECD, and ATc-ECD) were constructed on the basis of the normal model. Identical load conditions were applied to simulate flexion, extension, lateral bending, and axial rotation of the cervical spine. We calculated the range of motion (ROM), intradiscal pressure (IDP), annulus fibrosus pressure (AFP), uncovertebral joints contact pressure (CPRESS), and facet joints CPRESS under different motions. For all circumstances, ATc-ECD was close to the normal cervical spine model, whereas ATd-ECD significantly increased ROM and joints CPRESS and decreased IDP in the operated segment. PECF increased more the operated segment ROM than did the MACF, but the MACF obtained maximum IDP and AFP. Except for ATc-ECD, the other models increased joints CPRESS of the operated segment. For adjacent segments, ROM, IDP, and joints CPRESS showed a downward trend in all models. All models showed good biomechanical stability. With their combination biomechanics, safety, and conditions of application, PECF and ATc-ECD could be appropriate choices for cervical spondylotic radiculopathy.

Numerical Modeling of Shockwave Treatment of Knee Joint

Arthritis is a degenerative disease that primarily affects the cartilage and meniscus of the knee joint. External acoustic stimulation is used to treat this disease. This article presents a numerical model of the knee joint aimed at the computer-aided study of the regenerative effects of shockwave treatment. The presented model was verified and validated. A numerical analysis of the conditions for the regeneration of the tissues of the knee joint under shockwave action was conducted. The results allow us to conclude that to obtain the conditions required for the regeneration of cartilage tissues and meniscus (compressive stresses above the threshold value of 0.15 MPa to start the process of chondrogenesis; distortional strains above the threshold value of 0.05% characterized by the beginning of the differentiation of the tissues in large volumes; fluid pressure corresponding to the optimal level of 68 kPa to transfer tissue cells in large volumes), the energy flux density of therapeutic shockwave loading should exceed 0.3 mJ/mm².

A Reference Database of Standardised Continuous Lumbar Intervertebral Motion Analysis for Conducting Patient-Specific Comparisons

Lumbar instability has long been thought of as the failure of lumbar vertebrae to maintain their normal patterns of displacement. However, it is unknown what these patterns consist of. Research using quantitative fluoroscopy (QF) has shown that continuous lumbar intervertebral patterns of rotational displacement can be reliably measured during standing flexion and return motion using standardised protocols and can be used to assess patients with suspected lumbar spine motion disorders. However, normative values are needed to make individualised comparisons. One hundred and thirty-one healthy asymptomatic participants were recruited and performed guided flexion and return motion by following the rotating arm of an upright motion frame. Fluoroscopic image acquisition at 15fps was performed and individual intervertebral levels from L2-3 to L5-S1 were tracked and analysed during separate outward flexion and return phases. Results were presented as proportional intervertebral motion representing these phases using continuous means and 95%CIs, followed by verification of the differences between levels using Statistical Parametric Mapping (SPM). A secondary analysis of 8 control participants matched to 8 patients with chronic, non-specific low back pain (CNSLBP) was performed for comparison. One hundred and twenty-seven asymptomatic participants’ data were analysed. Their ages ranged from 18 to 70 years (mean 38.6) with mean body mass index 23.8 kg/m² 48.8% were female. Both the flexion and return phases for each level evidenced continuous change in mean proportional motion share, with narrow confidence intervals, highly significant differences and discrete motion paths between levels as confirmed by SPM. Patients in the secondary analysis evidenced significantly less L5-S1 motion than controls (p < 0.05). A reference database of spinal displacement patterns during lumbar (L2-S1) intersegmental flexion and return motion using a standardised motion protocol using fluoroscopy is presented. Spinal displacement patterns in asymptomatic individuals were found to be distinctive and consistent for each intervertebral level, and to continuously change during bending and return. This database may be used to allow continuous intervertebral kinematics to drive dynamic models of joint and muscular forces as well as reference values against which to make patient-specific comparisons in suspected cases of lumbar spine motion disorders.

Heterogeneous material mapping methods for patient-specific finite element models of pelvic trabecular bone: A convergence study

Patient-specific finite element (FE) models of bone require the assignment of heterogeneous material properties extracted from the subject's computed tomography (CT) images. Though node-based (NB) and element-based (EB) material mapping methods (MMMs) have been proposed, the sensitivity and convergence of FE models to MMM for varying mesh sizes are not well understood. In this work, CT-derived and synthetic bone material data were used to evaluate the effect of MMM on results from FE analyses. Pelvic trabecular bone data was extracted from CT images of six subjects, while synthetic data were created to resemble trabecular bone properties. The numerical convergence of FE bone models using different MMMs were evaluated for strain energy, von-Mises stress, and strain. NB and EB MMMs both demonstrated good convergence regarding total strain energy, with the EB method having a slight edge over the NB. However, at the local level (e.g., maximum stress and strain), FE results were sensitive to the field type, MMM, and the FE mesh size. The EB method exhibited superior performance in finer meshes relative to the voxel size. The NB method converged better than did the EB method for coarser meshes. These findings may lead to higher-fidelity patient-specific FE bone models.

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.

Biomechanical modeling of spinal ligaments: finite element analysis of L4-L5 spinal segment

The complex mechanical structure of spine is usually simplified in finite element (FE) modes. In this study, different 3D models of L4-L5 spinal segment distinguished by their ligament modelling were developed (1D truss, 2D shell and 3D space truss elements). All models could be considered validated with respect to range of motion and intradiscal pressure, although their ligament stresses/forces were substantially different. The models with 2D shell and 3D space truss ligaments showed the stress distribution and identified the potential failure/injury locations in ligaments. The model with 3D space truss ligaments showed the stress/force direction (representing collagen fiber directions).

Biomechanical Evaluation of Stand-Alone Oblique Lateral Lumbar Interbody Fusion Under 3 Different Bone Mineral Density Conditions: A Finite Element Analysis

OBJECTIVE

To evaluate the biomechanical stability of stand-alone (SA) oblique lateral interbody fusion (OLIF) under different bone mineral density conditions.

METHODS

The finite element model of L2-L5 was reconstructed and verified via computed tomography scan images (M0). The L4-L5 segment of SA OLIF was created based on the validation model. By changing bone mineral density, SA OLIF was established in the normal bone mineral density group (M1), osteopenia group (M2), and osteoporosis group (M3). A 500 N vertical axial preload was imposed on the superior surface of L2, and a 10 N-m moment was applied on the L2 superior surface along the radial direction to simulate 6 different physiological motions: flexion, extension, left and right lateral bending, left and right rotation.

RESULTS

Compared with M0, the range of motion of the fusion segment was significantly reduced, and the maximum stress of the upper and lower end plates was significantly increased in all motion modes. Compared with M1, the maximum relative increases of range of motion, cephalic end-plate stress and tail end-plate stress of M2 in the L4-L5 segment were 39.1%, 9.9%, and 10.7%, and the maximum increases of the above parameters in M3 were 100%, 28.9%, and 31.6%. The maximum stress of the tail end plate of the M3 model during flexion was 54.617 MPa, which was very close to the yield stress of the lamellar bone (60 MPa).

CONCLUSIONS

With the increase of the degree of osteoporosis, the maximum stress on the upper and lower end plates of the fusion segment increased significantly, thus increasing the potential risk of implant subsidence. SA OLIF could not provide sufficient stability for patients with osteoporosis.

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.

Experimental and Computational Comparison of Intervertebral Disc Bulge for Specimen-Specific Model Evaluation Based on Imaging

Finite element modelling of the spinal unit is a promising preclinical tool to assess the biomechanical outcome of emerging interventions. Currently, most models are calibrated and validated against range of motion and rarely directly against soft-tissue deformation. The aim of this contribution was to develop an in vitro methodology to measure disc bulge and assess the ability of different specimen-specific modelling approaches to predict disc bulge. Bovine bone-disc-bone sections (N = 6) were prepared with 40 glass markers on the intervertebral disc surface. These were initially magnetic resonance (MR)-imaged and then sequentially imaged using peripheral-qCT under axial compression of 1 mm increments. Specimen-specific finite-element models were developed from the CT data, using three different methods to represent the nucleus pulposus geometry with and without complementary use of the MR images. Both calibrated specimen-specific and averaged compressive material properties for the disc tissues were investigated. A successful methodology was developed to quantify the disc bulge in vitro, enabling observation of surface displacement on qCT. From the finite element model results, no clear advantage was found in using geometrical information from the MR images in terms of the models' ability to predict stiffness or disc bulge for bovine intervertebral disc.

Biomechanical Effect of C5 /C6 Intervertebral Reconstructive Height on Adjacent Segments in Anterior Cervical Discectomy and Fusion - A Finite Element Analysis

OBJECTIVE

To investigate the biomechanical effect of different intervertebral reconstructive heights on adjacent segments following C₅ /C₆ anterior cervical discectomy and fusion (ACDF) through finite element analysis.

METHODS

A finite element model of intact C₄ -C₇ segments was developed and validated for the present study. Five additional C₄ -C₇ postoperative models were constructed with 100%, 125%, 150%, 175%, and 200% of the benchmark height of C₅ /C₆ on the basis of the intact model. The changes in intradiscal pressure (IDP) and range of motion (ROM) of adjacent segments before and after reconstruction of C₅ /C₆ were analyzed.

RESULTS

For the upper adjacent segment (C₄ /C₅ ), the IDPs under the different loading conditions all increased after reconstruction. The maximum IDPs were 0.387, 0.489, 0.491, and 0.472 MPa under flexion, extension, axial rotation, and lateral bending, respectively, observed at the reconstructive height of 200%. The minimum IDPs were observed at 150% reconstructive height under all loading conditions except extension, and were 57, 86 and 81% of the maximum IDPs under flexion, axial rotation, and lateral bending, respectively. The minimum IDP under extension occurred when the reconstructive height is 125% of the benchmark height. For the lower adjacent segment (C₆ /C₇ ), the IDPs of postoperative models under all loading conditions also increased compared to the preoperative model. The maximum IDPs after reconstruction under flexion, extension, axial rotation, and lateral bending were 0.402, 0.411, 0.461, and 0.497 MPa, respectively, when the height of the reconstruction was 200% of the benchmark. The minimum IDPs were observed after a reconstruction at 150% of the benchmark, and were 59%, 85%, 82%, and 81% of the maximum IDPs under flexion, extension, axial rotation, and lateral bending loading conditions.

CONCLUSIONS

The reconstructive height is an important factor affecting the IDP and the ROM of adjacent segments after ACDF. To delay the adjacent segment disease, an intervertebral reconstructive height of 150% is an appropriate height in C₅ /C₆ ACDF.

Calibration and validation of a novel hybrid model of the lumbosacral spine in ArtiSynth-The passive structures

In computational biomechanics, two separate types of models have been used predominantly to enhance the understanding of the mechanisms of action of the lumbosacral spine (LSS): Finite element (FE) and musculoskeletal multibody (MB) models. To combine advantages of both models, hybrid FE-MB models are an increasingly used alternative. The aim of this paper is to develop, calibrate, and validate a novel passive hybrid FE-MB open-access simulation model of a ligamentous LSS using ArtiSynth. Based on anatomical data from the Male Visible Human Project, the LSS model is constructed from the L1-S1 rigid vertebrae interconnected with hyperelastic fiber-reinforced FE intervertebral discs, ligaments, and facet joints. A mesh convergence study, sensitivity analyses, and systematic calibration were conducted with the hybrid functional spinal unit (FSU) L4/5. The predicted mechanical responses of the FSU L4/5, the lumbar spine (L1-L5), and the LSS were validated against literature data from in vivo and in vitro measurements and in silico models. Spinal mechanical responses considered when loaded with pure moments and combined loading modes were total and intervertebral range of motions, instantaneous axes and centers of rotation, facet joint contact forces, intradiscal pressures, disc bulges, and stiffnesses. Undesirable correlations with the FE mesh were minimized, the number of crisscrossed collagen fiber rings was reduced to five, and the individual influences of specific anatomical structures were adjusted to in vitro range of motions. Including intervertebral motion couplings for axial rotation and nonlinear stiffening under increasing axial compression, the predicted kinematic and structural mechanics responses were consistent with the comparative data. The results demonstrate that the hybrid simulation model is robust and efficient in reproducing valid mechanical responses to provide a starting point for upcoming optimizations and extensions, such as with active skeletal muscles.

Biomechanical effect of different plate-to-disc distance on surgical and adjacent segment in anterior cervical discectomy and fusion - a finite element analysis

Background

The plate-to-disc distance (PDD) is an important factor affecting the degeneration of adjacent segments after anterior cervical discectomy and fusion (ACDF). However, the most suitable PDD is controversial. This study examined the adjacent intervertebral disc stress, bone graft stress, titanium plate stress and screw stress to evaluate the biomechanical effect of different PDD on surgical segment and adjacent segment following C5/C6 ACDF.

Methods

We constructed 10 preoperative finite element models of intact C4–C7 segments and validated them in the present study. We simulated ACDF surgery based on the 10 intact models in software. We designed three different distance of plate-to-disc titanium plates: long PDD (10 mm), medium PDD (5 mm) and short PDD (0 mm). The changes in C4/C5 and C6/C7 intervertebral disc stress, bone graft stress, titanium plate stress and screw stress were analyzed.

Results

The von Mises stress of C4/C5 and C6/C7 intervertebral discs had no significant differences (P > 0.05) in three different PDD groups. Titanium plate stress increased as the PDD decreased. The bone graft stress and screws stress decreased as the PDD decreased. The maximum stress of each part occurred was mostly in the conditions of rotation and lateral bending.

Conclusions

The PDD has no effect on adjacent intervertebral disc stress, but it is an important factor that affecting the bone graft stress, titanium plate stress and screws stress after ACDF. Shorter PDD plate can provide better stability to reduce stress on screws and bone graft, which may be helpful to prevent cage subsidence, pseudarthrosis and instrument failure. This can serve as a reference for clinical choice of plate.

Biomechanical influence of the surgical approaches, implant length and density in stabilizing ankylosing spondylitis cervical spine fracture

Ankylosing spondylitis cervical spine fractures (ASCFs) are particularly unstable and need special consideration when selecting appropriate internal fixation technology. However, there is a lack of related biomechanical studies. This study aimed to investigate the biomechanical influence of the pattern, length, and density of instrumentation for the treatment of ASCF. Posterior, anterior, and various combined fixation approaches were constructed using the finite element model (FEM) to mimic the surgical treatment of ASCFs at C5/6. The rate of motion change (RMC) at the fractured level and the internal stress distribution (ISD) were observed. The results showed that longer segments of fixation and combined fixation approaches provided better stability and lowered the maximal stress. The RMC decreased more significantly when the length increased from 1 to 3 levels (302% decrease under flexion, 134% decrease under extension) than from 3 to 5 levels (22% decrease under flexion, 23% decrease under extension). Longer fixation seems to be more stable with the anterior/posterior approach alone, but 3-level posterior fixation may be the most cost-effective option. It is recommended to perform surgery with combined approaches, which provide the best stability. Long skipped-screwing posterior fixation is an alternative technique for use in ASCF patients.