FEBio finite element models of the human lumbar spine

The biomechanical effect of lumbopelvic distance reduction on reconstruction after total sacrectomy: a comparative finite element analysis of four techniques

BACKGROUND CONTEXT

Following total sacrectomy, lumbopelvic reconstruction is essential to restore continuity between the lumbar spine and pelvis. However, to achieve long-term clinical stability, bony fusion between the lumbar spine and the pelvic ring is crucial. Reduction of the lumbopelvic distance can promote successful bony fusion. Although many lumbopelvic reconstruction techniques (LPRTs) have been previously analyzed, the biomechanical effect of lumbopelvic distance reduction (LPDR) has not been investigated yet.

PURPOSE

To evaluate and compare the biomechanical characteristics of four different LPRTs while considering the effect of LPDR.

STUDY DESIGN/SETTING

A comparative finite element (FE) study.

METHODS

The FE models following total sacrectomy were developed to analyze four different LPRTs, with and without LPDR. The closed-loop reconstruction (CLR), the sacral-rod reconstruction (SRR), the four-rod reconstruction (FRR), and the improved compound reconstruction (ICR) techniques were analyzed in flexion, extension, lateral bending, and axial rotation. Lumbopelvic stability was assessed through the shift-down displacement and the relative sagittal rotation of L5, while implant safety was evaluated based on the stress state at the bone-implant interface and within the rods.

RESULTS

Regardless of LPDR, both the shift-down displacement and relative sagittal rotation of L5 consistently ranked the LPRTs as ICR

CONCLUSIONS

LPDR significantly improved both lumbopelvic stability and implant safety in all reconstruction techniques after total sacrectomy. LPDR reduced the shift-down displacement of L5, the relative sagittal rotation of L5, and the stress values at the bone-implant interface. Furthermore, in the ICR and SRR techniques, LPDR decreased the peak stress values within the rods. All four investigated LPRTs demonstrated suitability for lumbopelvic reconstruction, with the ICR technique exhibiting the highest lumbopelvic stiffness.

CLINICAL SIGNIFICANCE

LPDR creates a biomechanically advantageous environment following total sacrectomy; therefore, it has the potential to impact the design of custom-made 3D-printed or traditional LPRTs. However, to confirm the findings of the current FE study, long-term clinical trials are recommended.

Collapse

Key Words

• Biomechanics
• Finite element analysis
• Lumbopelvic distance reduction
• Lumbopelvic fixation
• Lumbopelvic reconstruction
• Lumbopelvic stabilization
• Total sacrectomy

Collapse

MESH Headings

• Humans
• Finite Element Analysis
• Biomechanical Phenomena
• Sacrum/surgery
• Plastic Surgery Procedures/methods
• Lumbar Vertebrae/surgery
• Male
• Range of Motion, Articular
• Spinal Fusion/methods

Collapse

Grants Collapse

Affiliation(s)

Mate Turbucz School of PhD Studies, Semmelweis University, Üllői Str. 26, Budapest, Hungary; In Silico Biomechanics Laboratory, National Center for Spinal Disorders, Királyhágó Str. 1-3, Budapest, Hungary
Agoston Jakab Pokorni School of PhD Studies, Semmelweis University, Üllői Str. 26, Budapest, Hungary; In Silico Biomechanics Laboratory, National Center for Spinal Disorders, Királyhágó Str. 1-3, Budapest, Hungary
Benjamin Hajnal School of PhD Studies, Semmelweis University, Üllői Str. 26, Budapest, Hungary; In Silico Biomechanics Laboratory, National Center for Spinal Disorders, Királyhágó Str. 1-3, Budapest, Hungary
Kristof Koch School of PhD Studies, Semmelweis University, Üllői Str. 26, Budapest, Hungary; National Center for Spinal Disorders, Királyhágó Str. 1-3, Budapest, Hungary
Zsolt Szoverfi National Center for Spinal Disorders, Királyhágó Str. 1-3, Budapest, Hungary
Peter Pal Varga National Center for Spinal Disorders, Királyhágó Str. 1-3, Budapest, Hungary
Aron Lazary National Center for Spinal Disorders, Királyhágó Str. 1-3, Budapest, Hungary; Department of Spine Surgery, Department of Orthopaedics, Semmelweis University, Üllői Str. 78/b, Budapest, Hungary
Peter Endre Eltes In Silico Biomechanics Laboratory, National Center for Spinal Disorders, Királyhágó Str. 1-3, Budapest, Hungary; National Center for Spinal Disorders, Királyhágó Str. 1-3, Budapest, Hungary; Department of Spine Surgery, Department of Orthopaedics, Semmelweis University, Üllői Str. 78/b, Budapest, Hungary.

Collapse

The effect of polymethylmethacrylate augmentation on the primary stability of stand-alone implant construct versus posterior stabilization in oblique lumbar interbody fusion with osteoporotic bone quality- a finite element study

BACKGROUND CONTEXT

Oblique lumbar interbody fusion (OLIF) can provide an ideal minimally invasive solution for achieving spinal fusion in an older, more frail population where decreased bone quality can be a limiting factor. Stabilization can be achieved with bilateral pedicle screws (BPS), which require additional incisions and longer operative time. Alternatively, a novel self-anchoring stand-alone lateral plate system (SSA) can be used, where no additional incisions are required. Based on the relevant literature, BPS constructs provide greater primary biomechanical stability compared to lateral plate constructs, including SSA. This difference is further increased by osteoporosis. Screw augmentation in spinal fusion surgeries is commonly used; however, in the case of OLIF, it is a fairly new concept, lacking a consensus-based guideline.

PURPOSE

This comparative finite element (FE) study aimed to investigate the effect of PMMA screw augmentation on the primary stability of a stand-alone implant construct versus posterior stabilization in OLIF with osteoporotic bone quality.

STUDY DESIGN

The biomechanical effect of screw augmentation was studied inside an in-silico environment using computer-aided FE analysis.

METHODS

A previously validated and published L2-L4 FE model with normal and osteoporotic bone material properties was used. Geometries based on the OLIF implants (BPS, SSA) were created and placed inside the L3-L4 motion segment with increasing volumes (1-6 cm3) of PMMA augmentation. A follower load of 400 N and 10 Nm bending moment (in the three anatomical planes) were applied to the surgical FE models with different bone material properties. The operated L3-L4 segmental range of motion (ROM), the inserted cage's maximal caudal displacements, and L4 cranial bony endplate principal stress values were measured.

RESULTS

The nonaugmented values for the BPS construct were generally lower compared to SSA, and the difference was increased by osteoporosis. In osteoporotic bone, PMMA augmentation gradually decreased the investigated parameters and the difference between the two constructs as well. Between 3 cm3 and 4 cm3 of injected PMMA volume per screw, the difference between augmented SSA and standard BPS became comparable.

CONCLUSIONS

Based on this study, augmentation can enhance the primary stability of the constructs and decrease the difference between them. Considering leakage as a possible complication, between 3 cm3 and 4 cm3 of injected PMMA per screw can be an adequate amount for SSA augmentation. However, further in silico, and possibly in vitro and clinical testing is required to thoroughly understand the investigated biomechanical aspects.

CLINICAL SIGNIFICANCE

This study sheds light on the possible biomechanical advantage offered by augmented OLIF implants and provides a theoretical augmentation amount for the SSA construct. Based on the findings, the concept of an SSA device with PMMA augmentation capability is desirable.

Dataset of Finite Element Models of Normal and Deformed Thoracolumbar Spine

Adult spine deformity (ASD) is prevalent and leads to a sagittal misalignment in the vertebral column. Computational methods, including Finite Element (FE) Models, have emerged as valuable tools for investigating the causes and treatment of ASD through biomechanical simulations. However, the process of generating personalised FE models is often complex and time-consuming. To address this challenge, we present a dataset of FE models with diverse spine morphologies that statistically represent real geometries from a cohort of patients. These models are generated using EOS images, which are utilized to reconstruct 3D surface spine models. Subsequently, a Statistical Shape Model (SSM) is constructed, enabling the adaptation of a FE hexahedral mesh template for both the bone and soft tissues of the spine through mesh morphing. The SSM deformation fields facilitate the personalization of the mean hexahedral FE model based on sagittal balance measurements. Ultimately, this new hexahedral SSM tool offers a means to generate a virtual cohort of 16807 thoracolumbar FE spine models, which are openly shared in a public repository.

Biomechanical Analysis of a Newly Proposed Surgical Combination (MIS Screw-Rod System for Indirect Decompression+ Interspinous Fusion System for long Term Spinal Stability) in Treatment of Lumbar Degenerative Diseases

OBJECTIVE

The aim of this study is to analyze the biomechanical stability of a newly proposed surgical combination (minimally invasive surgery of screw-rod system for indirect decompression + interspinous fusion system for long term spinal stability) in treatment of lumbar degenerative diseases.

METHODS

The three-dimensional (3D) computed tomography (CT) image data of an adult healthy male volunteer were selected. An intact model of L4/5 was further established and validated by using Mimic and 3-matic, 3D slicer, abaqus, Python. Four surgical models were constructed. The biomechanical stability among these surgical modes was compared and analyzed using finite element analysis.

RESULTS

The maximum von mises on fixation system in surgical models 2 and 3 exhibited comparable values. This finding suggested that the increase in interspinous fusion did not result in a significant elevation in maximum von mises on fixation system. Compared with the third surgical model, the fourth model, which received less average von mises experienced by the screw in contact with both cancellous and cortical bone. The findings indicated that the inclusion of facet joint fusion in surgical procedures might not be necessary to increase the average von Mises stress experienced by the screw in contact with both cancellous and cortical bone.

CONCLUSIONS

The biomechanical stability of the newly proposed surgical combination (MIS screw-rod for indirect decompression + interspinous fusion for long term spinal stability technique) was not lower than that of the other surgical combination groups, and it might not be necessary to perform facet joint fusion during the surgery.

A systematic comparison between FEBio and PolyFEM for biomechanical systems

BACKGROUND AND OBJECTIVES

Finite element simulations are widely employed as a non-invasive and cost-effective approach for predicting outcomes in biomechanical simulations. However, traditional finite element software, primarily designed for engineering materials, often encountered limitations in contact detection and enforcement, leading to simulation failure when dealing with complex biomechanical configurations. Currently, a lot of model tuning is required to get physically accurate finite element simulations without failures. This adds significant human interaction to each iteration of a biomechanical model. This study addressed these issues by introducing PolyFEM, a novel finite element solver that guarantees inversion- and intersection-free solutions with completely automatic collision detection. The objective of this research is to validate PolyFEM's capabilities by comparing its results with those obtained from a well-established finite element solver, FEBio.

METHODS

To achieve this goal, five comparison scenarios were formulated to assess and validate PolyFEM's performance. The simulations were reproduced using both PolyFEM and FEBio, and the final results were compared. The five comparison scenarios included: (1) reproducing simulations from the FEBio test suite, consisting of static, dynamic, and contact-driven simulations; (2) replicating simulations from the verification paper published alongside the original release of FEBio; (3) a biomechanically based contact problem; (4) creating a custom simulation involving high-energy collisions between soft materials to highlight the difference in collision methods between the two solvers; and (5) performing biomechanical simulations of biting and quasi-stance.

RESULTS

We found that PolyFEM was capable of replicating all simulations previously conducted in FEBio. Particularly noteworthy is PolyFEM's superiority in high-energy contact simulations, where FEBio fell short, unable to complete over half of the simulations in Scenario 4. Although some of the simulations required significantly more simulation time in PolyFEM compared to FEBio, it is important to highlight that PolyFEM achieved these results without the need for any additional model tuning or contact declaration.

DISCUSSION

Despite being in the early stages of development, PolyFEM currently provides verified solutions for hyperelastic materials that are consistent with FEBio, both in previously published workflows and novel finite element scenarios. PolyFEM exhibited the ability to tackle challenging biomechanical problems where other solvers fell short, thus offering the potential to enhance the accuracy and realism of future finite element analyses.

Comparison of anterior column reconstruction techniques after en bloc spondylectomy: a finite element study

Total en bloc spondylectomy (TES) effectively treats spinal tumors. The surgery requires a vertebral body replacement (VBR), for which several solutions were developed, whereas the biomechanical differences between these devices still need to be completely understood. This study aimed to compare a femur graft, a polyetheretherketone implant (PEEK-IMP-C), a titan mesh cage (MESH-C), and a polymethylmethacrylate replacement (PMMA-C) using a finite element model of the lumbar spine after a TES of L3. Several biomechanical parameters (rotational stiffness, segmental range of motion (ROM), and von Mises stress) were assessed to compare the VBRs. All models provided adequate initial stability by increasing the rotational stiffness and decreasing the ROM between L2 and L4. The PMMA-C had the highest stiffness for flexion-extension, lateral bending, and axial rotation (215%, 216%, and 170% of intact model), and it had the lowest segmental ROM in the instrumented segment (0.2°, 0.5°, and 0.7°, respectively). Maximum endplate stress was similar for PMMA-C and PEEK-IMP-C but lower for both compared to MESH-C across all loading directions. These results suggest that PMMA-C had similar or better primary spinal stability than other VBRs, which may be related to the larger contact surface and the potential to adapt to the patient's anatomy.

Biomechanical effects of individualized artificial titanium alloy lamina implantation after laminectomy: A finite element analysis

Background and objectives: Laminectomy is a common surgical procedure in spine surgery. However, disruption of the posterior ligamentous complex of the spine may lead to a range of postoperative complications. Artificial lamina as a kind of bionic implant can well restore the posterior spinal structure. In this study, an individualized artificial titanium alloy lamina was designed to reconstruct the posterior spinal structure after laminectomy and explored its biomechanical effects, which could provide a theoretical basis for the clinical application of the artificial lamina.

Methods: Three finite element models were constructed, namely the nonlinear and non-homogeneous intact model of the whole lumbar spine, the lumbar decompression alone surgical model, and the artificial lamina implantation surgical model. The range of motion, intradiscal pressure, and annulus fibrosus peak stress were compared between the three models at the surgical and adjacent segments. The stresses of the artificial lamina and fixation screws were also analyzed for the four movement states.

Results: Compared with the intact model, the lumbar decompression alone surgical model showed an increase in range of motion, intradiscal pressure, and annulus fibrosus peak stresses at the surgical segment and adjacent segments under all conditions. The artificial lamina implantation surgical model showed an increase in these measurements only in flexion, increasing by 7.5%–22.5%, 7.6%–17.9%, and 6.4%–19.3%, respectively, over the intact model, while there was little difference under other conditions. The peak stresses in both the screw and the artificial lamina were highest in axial rotation, i. e. 46.53 MPa and 53.84 MPa, respectively. Screw stresses were concentrated on the connection between the screw and the artificial lamina, and artificial lamina stresses were concentrated on the spinous root, around the screw hole, and the contact with the vertebral body.

Conclusion: An individualized artificial titanium alloy lamina can effectively reduce the range of motion, intradiscal pressure, and annulus fibrosus stress at the surgical segment and adjacent segments. The application of artificial lamina could better preserve the biomechanical properties of the intact lumbar spine and reduce the risk of adjacent segmental disease.

Design of personalized scoliosis braces based on differentiable biomechanics—Synthetic study

This work describes a computational methodology for the design of braces for adolescent idiopathic scoliosis. The proposed methodology relies on a personalized simulation model of the patient’s trunk, and automatically searches for the brace geometry that optimizes the trade-off between clinical improvement and patient comfort. To do this, we introduce a formulation of differentiable biomechanics of the patient’s trunk, the brace, and their interaction. We design a simulation model that is differentiable with respect to both the deformation state and the brace design parameters, and we show how this differentiable model is used for the efficient update of brace design parameters within a numerical optimization algorithm. To evaluate the proposed methodology, we have obtained trunk models with personalized geometry for five patients of adolescent idiopathic scoliosis, and we have designed Boston-type braces. In a simulation setting, the designed braces improve clinical metrics by 45% on average, under acceptable comfort conditions. In the future, the methodology can be extended beyond synthetic validation, and tested with physical braces on the actual patients.

A finite element study on the effects of follower load on the continuous biomechanical responses of subaxial cervical spine

In spine biomechanics, follower loads are used to mimic the in vivo muscle forces acting on a human spine. However, the effects of the follower load on the continuous biomechanical responses of the subaxial cervical spines (C2-T1) have not been systematically clarified. This study aims at investigating the follower load effects on the continuous biomechanical responses of C2-T1. A nonlinear finite element model is reconstructed and validated for C2-T1. Six levels follower loads are considered along the follower load path that is optimized through a novel range of motion-based method. A moment up to 2 Nm is subsequently superimposed to produce motions in three anatomical planes. The continuous biomechanical responses, including the range of motion, facet joint force, intradiscal pressure and flexibility are evaluated for each motion segment. In the sagittal plane, the change of the overall range of motion arising from the follower loads is less than 6%. In the other two anatomical planes, both the magnitude and shape of the rotation-moment curves change with follower loads. At the neutral position, over 50% decrease in flexibility occurs as the follower load increases from zero to 250 N. In all three anatomical planes, over 50% and 30% decreases in flexibility occur in the first 0.5 Nm for small (≤100 N) and large (≥150 N) follower loads, respectively. Moreover, follower loads tend to increase both the facet joint forces and the intradiscal pressures. The shape of the intradiscal pressure-moment curves changes from nonlinear to roughly linear with increased follower load, especially in the coronal and transverse planes. The results obtained in this work provide a comprehensive understanding on the effects of follower load on the continuous biomechanical responses of the C2-T1.

Low-Modulus PMMA Has the Potential to Reduce Stresses on Endplates after Cement Discoplasty

Cement discoplasty has been developed to treat patients with advanced intervertebral disc degeneration. In discoplasty, poly(methylmethacrylate) (PMMA) bone cement is injected into the disc, leading to reduced pain and certain spinal alignment correction. Standard PMMA-cements have much higher elastic modulus than the surrounding vertebral bone, which may lead to a propensity for adjacent fractures. A PMMA-cement with lower modulus might be biomechanically beneficial. In this study, PMMA-cements with lower modulus were obtained using previously established methods. A commercial PMMA-cement (V-steady^®, G21 srl) was used as control, and as base cement. The low-modulus PMMA-cements were modified by 12 vol% (LA12), 16 vol% (LA16) and 20 vol% (LA20) linoleic acid (LA). After storage in 37 °C PBS from 24 h up to 8 weeks, specimens were tested in compression to obtain the material properties. A lower E-modulus was obtained with increasing amount of LA. However, with storage time, the E-modulus increased. Standard and low-modulus PMMA discoplasty were compared in a previously developed and validated computational lumbar spine model. All discoplasty models showed the same trend, namely a substantial reduction in range of motion (ROM), compared to the healthy model. The V-steady model had the largest ROM-reduction (77%), and the LA20 model had the smallest (45%). The average stress at the endplate was higher for all discoplasty models than for the healthy model, but the stresses were reduced for cements with higher amounts of LA. The study indicates that low-modulus PMMA is promising for discoplasty from a mechanical viewpoint. However, validation experiments are needed, and the clinical setting needs to be further considered.

FEBio finite element model of a pediatric cervical spine

OBJECTIVE

The underlying biomechanical differences between the pediatric and adult cervical spine are incompletely understood. Computational spine modeling can address that knowledge gap. Using a computational method known as finite element modeling, the authors describe the creation and evaluation of a complete pediatric cervical spine model.

METHODS

Using a thin-slice CT scan of the cervical spine from a 5-year-old boy, a 3D model was created for finite element analysis. The material properties and boundary and loading conditions were created and model analysis performed using open-source software. Because the precise material properties of the pediatric cervical spine are not known, a published parametric approach of scaling adult properties by 50%, 25%, and 10% was used. Each scaled finite element model (FEM) underwent two types of simulations for pediatric cadaver testing (axial tension and cardinal ranges of motion [ROMs]) to assess axial stiffness, ROM, and facet joint force (FJF). The authors evaluated the axial stiffness and flexion-extension ROM predicted by the model using previously published experimental measurements obtained from pediatric cadaveric tissues.

RESULTS

In the axial tension simulation, the model with 50% adult ligamentous and annulus material properties predicted an axial stiffness of 49 N/mm, which corresponded with previously published data from similarly aged cadavers (46.1 ± 9.6 N/mm). In the flexion-extension simulation, the same 50% model predicted an ROM that was within the range of the similarly aged cohort of cadavers. The subaxial FJFs predicted by the model in extension, lateral bending, and axial rotation were in the range of 1-4 N and, as expected, tended to increase as the ligament and disc material properties decreased.

CONCLUSIONS

A pediatric cervical spine FEM was created that accurately predicts axial tension and flexion-extension ROM when ligamentous and annulus material properties are reduced to 50% of published adult properties. This model shows promise for use in surgical simulation procedures and as a normal comparison for disease-specific FEMs.

Stability Evaluation of Different Oblique Lumbar Interbody Fusion Constructs in Normal and Osteoporotic Condition - A Finite Element Based Study

Introduction: In developed countries, the age structure of the population is currently undergoing an upward shift, resulting a decrease in general bone quality and surgical durability. Over the past decade, oblique lumbar interbody fusion (OLIF) has been globally accepted as a minimally invasive surgical technique. There are several stabilization options available for OLIF cage fixation such as self-anchored stand-alone (SSA), lateral plate-screw (LPS), and bilateral pedicle screw (BPS) systems. The constructs’ stability are crucial for the immediate and long-term success of the surgery. The aim of this study is to investigate the biomechanical effect of the aforementioned constructs, using finite element analysis with different bone qualities (osteoporotic and normal).

Method: A bi-segmental (L2–L4) finite element (FE) model was created, using a CT scan of a 24-year-old healthy male. After the FE model validation, CAD geometries of the implants were inserted into the L3–L4 motion segment during a virtual surgery. For the simulations, a 150 N follower load was applied on the models, then 10 Nm of torque was used in six general directions (flexion, extension, right/left bending, and right/left rotation), with different bone material properties.

Results: The smallest segmental (L3–L4) ROM (range of motion) was observed in the BPS system, except for right bending. Osteoporosis increased ROMs in all constructs, especially in the LPS system (right bending increase: 140.26%). Osteoporosis also increased the caudal displacement of the implanted cage in all models (healthy bone: 0.06 ± 0.03 mm, osteoporosis: 0.106 ± 0.07 mm), particularly with right bending, where the displacement doubled in SSA and LPS constructs. The displacement of the screws inside the L4 vertebra increased by 59% on average (59.33 ± 21.53%) due to osteoporosis (100% in LPS, rotation). BPS-L4 screw displacements were the least affected by osteoporosis.

Conclusions: The investigated constructs provide different levels of stability to the spine depending on the quality of the bone, which can affect the outcome of the surgery. In our model, the BPS system was found to be the most stable construct in osteoporosis. The presented model, after further development, has the potential to help the surgeon in planning a particular spinal surgery by adjusting the stabilization type to the patient’s bone quality.

A study of the sensitivity of biomechanical models of the spine for scoliosis brace design

BACKGROUND AND OBJECTIVE

The development of biomechanical models of the torso and the spine opens the door to computational solutions for the design of braces for adolescent idiopathic scoliosis. However, the design of such biomechanical models faces several unknowns, such as the correct identification of relevant mechanical elements, or the required accuracy of model parameters. The objective of this study was to design a methodology for the identification of the aforementioned elements, with the purpose of creating personalized models suited for patient-specific brace design and the definition of parameter estimation criteria.

METHODS

We have developed a comprehensive model of the torso, including spine, ribcage and soft tissue, and we have developed computational tools for the analysis of the model parameters. With these tools, we perform an analysis of the model under typical loading conditions of scoliosis braces.

RESULTS

We present a complete sensitivity analysis of the models mechanical parameters and a comparison between a reference healthy subject and a subject suffering from scoliosis. Furthermore, we make a direct connection between error bounds on the deformation and tolerances for parameter estimation, which can guide the personalization of the model.

CONCLUSIONS

Not surprisingly, the stiffness parameters that govern the lateral deformation of the spine in the frontal plane are some of the most relevant parameters, and require careful modeling. More surprisingly, their relevance is on par with the correct parameterization of the soft tissue of the torso. For scoliosis patients, but not for healthy subjects, we observe that the axial rotation of the spine also requires careful modeling.

Multiaxial validation of a finite element model of the intervertebral disc with multigenerational fibers to establish residual strain

Finite element models of the intervertebral disc are used to address research questions that cannot be tested through typical experimentation. A disc model requires complex geometry and tissue properties to be accurately defined to mimic the physiological disc. The physiological disc possesses residual strain in the annulus fibrosus (AF) due to osmotic swelling and due to inherently pre-strained fibers. We developed a disc model with residual contributions due to swelling-only, and a multigeneration model with residual contributions due to both swelling and AF fiber pre-strain and validated it against organ-scale uniaxial, quasi-static and multiaxial, dynamic mechanical tests. In addition, we demonstrated the models' ability to mimic the opening angle observed following radial incision of bovine discs. Both models were validated against organ-scale experimental data. While the swelling only model responses were within the experimental 95% confidence interval, the multigeneration model offered outcomes closer to the experimental mean and had a bovine model opening angle within one SD of the experimental mean. The better outcomes for the multigeneration model, which allowed for the inclusion of inherently pre-strained fibers in AF, is likely due to its uniform fiber contribution throughout the AF. We conclude that the residual contribution of pre-strained fibers in the AF should be included to best simulate the physiological disc and its behaviors.

Biomechanical Effects of a Cross Connector in Sacral Fractures - A Finite Element Analysis

Background: Spinopelvic fractures and approaches of operative stabilization have been a source of controversial discussion. Biomechanical data support the benefit of a spinopelvic stabilization and minimally invasive procedures help to reduce the dissatisfying complication rate. The role of a cross connector within spinopelvic devices remains inconclusive. We aimed to analyze the effect of a cross connector in a finite element model (FE model). Study Design: A FE model of the L1-L5 spine segment with pelvis and a spinopelvic stabilization was reconstructed from patient-specific CT images. The biomechanical relevance of a cross connector in a Denis zone I (AO: 61-B2) sacrum fracture was assessed in the FE model by applying bending and twisting forces with and without a cross connector. Biomechanical outcomes from the numerical model were investigated also considering uncertainties in material properties and levels of osseointegration. Results: The designed FE model showed comparable values in range-of-motion (ROM) and stresses with reference to the literature. The superiority of the spinopelvic stabilization (L5/Os ilium) ± cross connector compared to a non-operative procedure was confirmed in all analyzed loading conditions by reduced ROM and principal stresses in the disk L5/S1, vertebral body L5 and the fracture area. By considering the combination of all loading cases, the presence of a cross connector reduced the maximum stresses in the fracture area of around 10%. This difference has been statistically validated (p < 0.0001). Conclusion: The implementation of a spinopelvic stabilization (L5/Os ilium) in sacrum fractures sustained the fracture and led to enhanced biomechanical properties compared to a non-reductive procedure. While the additional cross connector did not alter the resulting ROM in L4/L5 or L5/sacrum, the reduction of the maximum stresses in the fracture area was significant.

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.

Computational Challenges in Tissue Engineering for the Spine

This paper deals with a brief review of the recent developments in computational modelling applied to innovative treatments of spine diseases. Additionally, it provides a perspective on the research directions expected for the forthcoming years. The spine is composed of distinct and complex tissues that require specific modelling approaches. With the advent of additive manufacturing and increasing computational power, patient-specific treatments have moved from being a research trend to a reality in clinical practice, but there are many issues to be addressed before such approaches become universal. Here, it is identified that the major setback resides in validation of these computational techniques prior to approval by regulatory agencies. Nevertheless, there are very promising indicators in terms of optimised scaffold modelling for both disc arthroplasty and vertebroplasty, powered by a decisive contribution from imaging methods.

Toward Patient Specific Models of Pediatric IVDs: A Parametric Study of IVD Mechanical Properties

Patient specific finite element (FE) modeling of the pediatric spine is an important challenge which offers to revolutionize the treatment of pediatric spinal pathologies, for example adolescent idiopathic scoliosis (AIS). In particular, modeling of the intervertebral disc (IVD) is a unique challenge due to its structural and mechanical complexity. This is compounded by limited ability to non-invasively interrogate key mechanical parameters of a patient's IVD. In this work, we seek to better understand the link between mechanical properties and mechanical behavior of patient specific FE models of the pediatric lumbar spine. A parametric study of IVD parameter was conducted, coupled with insights from current knowledge of the pediatric IVD. In particular, the combined effects of parameters was investigated. Recommendations are made toward areas of importance in patient specific FE modeling of the pediatric IVD. In particular, collagen fiber bundles of the IVD are found to dominate IVD mechanical behavior and are thus recommended as an area of primary focus for patient specific FE models. In addition, areas requiring further experimental research are identified. This work provides a valuable building block toward the development of patient specific models of the pediatric spine.

Automated Pipeline to Generate Anatomically Accurate Patient-Specific Biomechanical Models of Healthy and Pathological FSUs

State-of-the-art preoperative biomechanical analysis for the planning of spinal surgery not only requires the generation of three-dimensional patient-specific models but also the accurate biomechanical representation of vertebral joints. The benefits offered by computational models suitable for such purposes are still outweighed by the time and effort required for their generation, thus compromising their applicability in a clinical environment. In this work, we aim to ease the integration of computerized methods into patient-specific planning of spinal surgery. We present the first pipeline combining deep learning and finite element methods that allows a completely automated model generation of functional spine units (FSUs) of the lumbar spine for patient-specific FE simulations (FEBio). The pipeline consists of three steps: (a) multiclass segmentation of cropped 3D CT images containing lumbar vertebrae using the DenseVNet network, (b) automatic landmark-based mesh fitting of statistical shape models onto 3D semantic segmented meshes of the vertebral models, and (c) automatic generation of patient-specific FE models of lumbar segments for the simulation of flexion-extension, lateral bending, and axial rotation movements. The automatic segmentation of FSUs was evaluated against the gold standard (manual segmentation) using 10-fold cross-validation. The obtained Dice coefficient was 93.7% on average, with a mean surface distance of 0.88 mm and a mean Hausdorff distance of 11.16 mm (N = 150). Automatic generation of finite element models to simulate the range of motion (ROM) was successfully performed for five healthy and five pathological FSUs. The results of the simulations were evaluated against the literature and showed comparable ROMs in both healthy and pathological cases, including the alteration of ROM typically observed in severely degenerated FSUs. The major intent of this work is to automate the creation of anatomically accurate patient-specific models by a single pipeline allowing functional modeling of spinal motion in healthy and pathological FSUs. Our approach reduces manual efforts to a minimum and the execution of the entire pipeline including simulations takes approximately 2 h. The automation, time-efficiency and robustness level of the pipeline represents a first step toward its clinical integration.

Biomechanical effects of laminectomies in the human lumbar spine: a finite element study

BACKGROUND CONTEXT

Previous studies have analyzed the effect of laminectomy on intervertebral disc (IVD), facet-joint-forces (FJF), and range of motion (ROM), while only two have specifically reported stresses at the pars interarticularis (PI) with posterior element resection. These studies have been performed utilizing a single subject, questioning their applications to a broader population.

PURPOSE

We investigate the effect of graded PI resection in a three-dimensional manner on PI stress to provide surgical guidelines for avoidance of iatrogenic instability following lumbar laminectomy. Additionally, quantified FJF and IVD stresses can provide further insight into the development of adjacent segment disease.

STUDY DESIGN

Biomechanical finite element (FE) method investigation of the lumbar spine.

METHODS

FE models of the lumbar spine of three subjects were created using the open-source finite element software, FEBio. Single-level laminectomy, two-level laminectomy, and ventral-to-dorsal PI resection simulations were performed with varying degrees of PI resection from 0% to 75% of the native PI. These models were taken through cardinal ROM under standard loading conditions and PI stresses, FJF, IVD stresses, and ROM were analyzed.

RESULTS

The three types of laminectomy simulated in this study showed a nonlinear increase in PI stress with increased bone resection. Axial rotation generated the most stress at the PI followed by flexion, extension and lateral bending. At 75% bone resection all three types of laminectomy produced PI stresses that were near the ultimate strength of human cortical bone during axial rotation. FJF decreased with increased bone resection for the three laminectomies simulated. This was most notable in axial rotation followed by extension and lateral bending. IVD stresses varied greatly between the nonsurgical models and likewise the effect of laminectomy on IVD stresses varied between subjects. ROM was mostly unaffected by the laminectomies performed in this study.

CONCLUSIONS

Regarding the risk of iatrogenic spondylolisthesis, the combined results are sufficient evidence to suggest surgeons should be particularly cautious when PI resection exceeds 50% bone resection for all laminectomies included in this study. Lastly, the effects seen in FJF and IVD stresses indicate the degree to which the remainder of the spine must experience compensatory biomechanical changes as a result of the surgical intervention.

Febio finite element models of the human cervical spine

Finite element (FE) analysis has proven to be useful when studying the biomechanics of the cervical spine. Although many FE studies of the cervical spine have been published, they typically develop their models using commercial software, making the sharing of models between researchers difficult. They also often model only one part of the cervical spine. The goal of this study was to develop and evaluate three FE models of the adult cervical spine using open-source software and to freely provide these models to the scientific community. The models were created from computed tomography scans of 26-, 59-, and 64-year old female subjects. These models were evaluated against previously published experimental and FE data. Despite the fact that all three models were assigned identical material properties and boundary conditions, there was notable variation in their biomechanical behavior. It was therefore apparent that these differences were the result of morphological differences between the models.

Development of a finite element biomechanical whole spine model for analyzing lumbar spine loads under caudocephalad acceleration

Background:Spine injury risk due to military conflict is an ongoing concern among defense organizations throughout the world. A better understanding of spine biomechanics could assist in developing protection devices to reduce injuries caused by caudocephalad acceleration (+Gz) in under-body blasts (UBB). Although some finite element (FE) human models have demonstrated reasonable lumbar spine biofidelity, they were either partial spine models or not validated for UBB-type loading modes at the lumbar functional spinal unit (FSU) level, thus limiting their ability to analyze UBB-associated occupant kinematics.Methods:An FE functional representation of the human spine with simplified geometry was developed to study the lumbar spine responses under +Gz loading. Fifty-seven load curves obtained from post mortem human subject experiments were used to optimize the model.Results:The model was cumulatively validated for compression, flexion, extension, and anterior-, posterior-, and lateral-shears of the lumbar spine and flexion and extension of the cervical spine. The thoracic spine was optimized for flexion and compression. The cumulative CORrelation and Analysis (CORA) rating for the lumbar spine was 0.766 and the cervical spine was 0.818; both surpassed the 0.7 objective goal. The model's element size was confirmed as converged.Conclusions:An FE functional representation of the human spine was developed for +Gz lumbar load analysis. The lumbar and cervical spines were demonstrated to be quantitatively biofidelic to the FSU level for multi-directional loading and bending typically experienced in +Gz loading, filling the capability gap in current models.

A Novel Modelling Methodology Which Predicts the Structural Behaviour of Vertebral Bodies under Axial Impact Loading: A Finite Element and DIC Study

Cervical spine injuries (CSIs) arising from collisions are uncommon in contact sports, such as rugby union, but their consequences can be devastating. Several FE modelling approaches are available in the literature, but a fully calibrated and validated FE modelling framework for cervical spines under compressive dynamic-impact loading is still lacking and material properties are not adequately calibrated for such events. This study aimed to develop and validate a methodology for specimen-specific FE modelling of vertebral bodies under impact loading. Thirty-five (n = 35) individual vertebral bodies (VBs) were dissected from porcine spine segments, potted in bone cement and μCT scanned. A speckle pattern was applied to the anterior faces of the bones to allow digital image correlation (DIC), which monitored the surface displacements. Twenty-seven (n = 27) VBs were quasi-statically compressively tested to a load up to 10 kN from the cranial side. Specimen-specific FE models were developed for fourteen (n = 14) of the samples in this group. The material properties were optimised based on the experimental load-displacement data using a specimen-specific factor (kGSstatic) to calibrate a density to Young’s modulus relationship. The average calibration factor arising from this group was calculated (K¯GSstatic) and applied to a control group of thirteen (n = 13) samples. The resulting VB stiffnesses was compared to experimental findings. The final eight (n = 8) VBs were subjected to an impact load applied via a falling mass of 7.4kg at a velocity of 3.1ms⁻¹. Surface displacements and strains were acquired from the anterior VB surface via DIC, and the impact load was monitored with two load cells. Specimen-specific FE models were created for this dynamic group and material properties were assigned again based on the density–Young’s modulus relationship previously validated for static experiments, supplemented with an additional factor (KGSdynamic). The optimised conversion factor for quasi-static loading, K¯GSstatic, had an average of 0.033. Using this factor, the validation models presented an average numerical stiffness value 3.72% greater than the experimental one. From the dynamic loading experiments, the value for KGSdynamic was found to be 0.14, 4.2 times greater than K¯GSstatic. The average numerical stiffness was 2.3% greater than in the experiments. Almost all models presented similar stiffness variations and regions of maximum displacement to those observed via DIC. The developed FE modelling methodology allowed the creation of models which predicted both static and dynamic behaviour of VBs. Deformation patterns on the VB surfaces were acquired from the FE models and compared to DIC data, achieving high agreement. This methodology is now validated to be fully applied to create whole cervical spine models to simulate axial impact scenarios replicating rugby collision events.