Uncertainty analysis of material properties and morphology parameters in numerical models regarding the motion of lumbar vertebral segments

Validation of a patient-specific finite element analysis framework for identification of growing rod-failure regions in early onset scoliosis patients

PURPOSE

Growing rods are the gold-standard for treatment of early onset scoliosis (EOS). However, these implanted rods experience frequent fractures, requiring additional surgery. A recent study by the U.S. Food and Drug Administration (FDA) identified four common rod fracture locations. Leveraging this data, Agarwal et al. were able to correlate these fractures to high-stress regions using a novel finite element analysis (FEA) framework for one patient. The current study aims to further validate this framework through FEA modeling extended to multiple patients.

METHODS

Three patient-specific FEA models were developed to match the pre-operative patient data taken from both registry and biplanar radiographs. The surgical procedure was then simulated to match the post-operative deformity. Body weight and flexion bending (1 Nm) loads were then applied and the output stress data on the rods were analyzed.

RESULTS

Radiographic data showed fracture locations at the mid-construct, adjacent to the distal and tandem connector across the patients. Stress analysis from the FEA showed these failure locations matched local high-stress regions for all fractures observed. These results qualitatively validate the efficacy of the FEA framework by showing a decent correlation between localized high-stress regions and the actual fracture sites in the patients.

CONCLUSIONS

This patient-specific, in-silico framework has huge potential to be used as a surgical tool to predict sites prone to fracture in growing rod implants. This prospective information would therefore be vital for surgical planning, besides helping optimize implant design for reducing rod failures.

Comparative FEM study on intervertebral disc modeling: Holzapfel-Gasser-Ogden vs. structural rebars

Introduction: Numerical modeling of the intervertebral disc (IVD) is challenging due to its complex and heterogeneous structure, requiring careful selection of constitutive models and material properties. A critical aspect of such modeling is the representation of annulus fibers, which significantly impact IVD biomechanics. This study presents a comparative analysis of different methods for fiber reinforcement in the annulus fibrosus of a finite element (FE) model of the human IVD. Methods: We utilized a reconstructed L4-L5 IVD geometry to compare three fiber modeling approaches: the anisotropic Holzapfel-Gasser-Ogden (HGO) model (HGO fiber model) and two sets of structural rebar elements with linear-elastic (linear rebar model) and hyperelastic (nonlinear rebar model) material definitions, respectively. Prior to calibration, we conducted a sensitivity analysis to identify the most important model parameters to be calibrated and improve the efficiency of the calibration. Calibration was performed using a genetic algorithm and in vitro range of motion (RoM) data from a published study with eight specimens tested under four loading scenarios. For validation, intradiscal pressure (IDP) measurements from the same study were used, along with additional RoM data from a separate publication involving five specimens subjected to four different loading conditions. Results: The sensitivity analysis revealed that most parameters, except for the Poisson ratio of the annulus fibers and C01 from the nucleus, significantly affected the RoM and IDP outcomes. Upon calibration, the HGO fiber model demonstrated the highest accuracy (R2 = 0.95), followed by the linear (R2 = 0.89) and nonlinear rebar models (R2 = 0.87). During the validation phase, the HGO fiber model maintained its high accuracy (RoM R2 = 0.85; IDP R2 = 0.87), while the linear and nonlinear rebar models had lower validation scores (RoM R2 = 0.71 and 0.69; IDP R2 = 0.86 and 0.8, respectively). Discussion: The results of the study demonstrate a successful calibration process that established good agreement with experimental data. Based on our findings, the HGO fiber model appears to be a more suitable option for accurate IVD FE modeling considering its higher fidelity in simulation results and computational efficiency.

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.

Torque- and Muscle-Driven Flexion Induce Disparate Risks of In Vitro Herniation: A Multiscale and Multiphasic Structure-Based Finite Element Study

The intervertebral disc is a complex structure that experiences multiaxial stresses regularly. Disc failure through herniation is a common cause of lower back pain, which causes reduced mobility and debilitating pain, resulting in heavy socioeconomic burdens. Unfortunately, herniation etiology is not well understood, partially due to challenges in replicating herniation in vitro. Previous studies suggest that flexion elevated risks of herniation. Thus, the objective of this study was to use a multiscale and multiphasic finite element model to evaluate the risk of failure under torque- or muscle-driven flexion. Models were developed to represent torque-driven flexion with the instantaneous center of rotation (ICR) located on the disc, and the more physiologically representative muscle-driven flexion with the ICR located anterior of the disc. Model predictions highlighted disparate disc mechanics regarding bulk deformation, stress-bearing mechanisms, and intradiscal stress-strain distributions. Specifically, failure was predicted to initiate at the bone-disc boundary under torque-driven flexion, which may explain why endplate junction failure, instead of herniation, has been the more common failure mode observed in vitro. By contrast, failure was predicted to initiate in the posterolateral annulus fibrosus under muscle-driven flexion, resulting in consistent herniation. Our findings also suggested that muscle-driven flexion combined with axial compression could be sufficient for provoking herniation in vitro and in silico. In conclusion, this study provided a computational framework for designing in vitro testing protocols that can advance the assessment of disc failure behavior and the performance of engineered disc implants.

Influence of structural and material property uncertainties on biomechanics of intervertebral discs - Implications for disc tissue engineering

This study investigated how variations of structural and material properties of human intervertebral discs (IVDs) affect the biomechanical responses of the IVDs under simulated physiological loading conditions using a stochastic finite element (SFE) model. An SFE method, which combined an anatomic FE model of human lumbar L3-4 segment and probabilistic analysis of its structural and material properties, was used to generate a dataset of 500 random disc samples with varying structural and material properties. The sensitivity of the biomechanical responses, including intervertebral displacements/rotations, intradiscal pressures (IDP), fiber stresses and matrix strains of annulus fibrosus (AF), were systematically quantified under various physiological loading conditions, including a 500N compression and 7.5Nm moments in the 3 primary rotations. Significant variations of the IDPs, IVD displacements/rotations, and stress/strain distributions were found using the dataset of 500 ramdom disc samples. Under all the loading conditions, the IDPs were positively correlated with the Poisson's ratio of the NP (r = 0.46 to 0.75, p = 0.004-0.001) and negatively with the Young's modulus of the annulus matrix (r = -0.48 to -0.65, p = 0.003-0.001). The primary intervertebral rotations were significantly affected by the Young's modulus of the annulus matrix (r = -0.44 to -0.71, p = 0.001-0.032) and the orientations of the annular fibers (r = -0.45 to -0.69, p = 0.001-0.029). The heterogeneity of structures and material properties of the IVD had distinct effects on the biomechanical performances of the IVD. These data could help improve our understanding of the intrinsic biomechanics of the IVD and provide references for optimal design of tissue engineered discs by controlling structural and material properties of the disc components.

Towards a validated patient-specific computational modeling framework to identify failure regions in traditional growing rods in patients with early onset scoliosis

Background

While growing rods are an important contribution to early-onset scoliosis treatment, rod fractures are a common complication that require reoperations. A recent retrieval analysis study performed on failed traditional growing rods revealed that there are commonalities among patient characteristics based on the location of rod fracture. However, it remains unknown if these locations correspond to high stress regions in the implanted construct.

Methods

A patient-specific finite element scoliotic model was developed to match the pre-operative (pre-op) scoliotic curve of a patient as described in previously published articles, and by using the patient registry information along with biplanar radiographs. A dual stainless-steel traditional growing rod construct was implanted into this scoliotic model and the surgical procedure was simulated to match the post-operative (post-op) scoliotic curve parameters. Muscle stabilization and gravity was simulated through follower load application. Rod distraction magnitudes were chosen based on pre-op to post-op cobb angle correction, and flexion bending load was simulated to identify the high stress regions on the rods.

Results

The patient-specific finite element model identified two high stress regions on the posterior surface of the rods, one at mid construct and the other adjacent to the distal anchors. This correlated well with the data obtained from the retrieval analysis performed by researchers at U.S. Food and Drug Administration (FDA) which showed the posterior surface of the rod as the fracture initiation site, and the three locations of failure as mid-construct, adjacent to distal anchors, and adjacent to tandem connector.

Conclusions

The result of this study confirms that the high stress regions on the growing rods, as identified by the FEA, match the fracture prone sites identified in the retrieval analysis performed at the FDA. This proof-of-concept patient-specific approach can be used to predict sites prone to fracture in growing rods.

The importance of curve severity, type and instrumentation strategy in the surgical correction of adolescent idiopathic scoliosis: an in silico clinical trial on 64 cases

Adolescent idiopathic scoliosis is a three-dimensional deformity of the spine which is frequently corrected with the implantation of instrumentation with generally good or excellent clinical results; mechanical post-operative complications such as implant loosening and breakage are however relatively frequent. The rate of complications is associated with a lack of consensus about the surgical decision-making process; choices about the instrumentation length, the anchoring implants and the degree of correction are indeed mostly based on personal views and previous experience of the surgeon. In this work, we performed an in silico clinical trial on a large number of subjects in order to clarify which factors have the highest importance in determining the risk of complications by quantitatively analysing the mechanical stresses and loads in the instrumentation after the correction maneuvers. The results of the simulations highlighted the fundamental role of the curve severity, also in its three-dimensional aspect, and of the instrumentation strategy, whereas the length of the fixation had a lower importance.

Neck mobility in the Jurassic plesiosaur Cryptoclidus eurymerus: finite element analysis as a new approach to understanding the cervical skeleton in fossil vertebrates

The sauropterygian clade Plesiosauria arose in the Late Triassic and survived to the very end of the Cretaceous. Plesiosauria evolved the greatest species diversity of any marine reptile clade, attaining a global distribution. Plesiosauria consist of two clades, Rhomaleosauridae and Neoplesiosauria. Basal Neoplesiosauria have long necks with at least 30 cervicals, but show qualitative osteological evidence for a stiff neck. Here we quantify neck mobility in lateral, ventral, and dorsal directions based on finite element modeling of neck vertebrae from the Middle Jurassic plesiosaur Cryptoclidus eurymerus. We model the mobility in a single motion segment, consisting of two adjacent cervical vertebrae and the joints connecting them. Based on the model with a maximum intervertebral spacing of 3 mm, we find that in Cryptoclidus, the maximum angle of lateral deflection in the motion segment was 2°. The maximum angle of ventral deflection was 5° and of dorsal deflection was 5°. When these values are multiplied by the number of cervical vertebrae, it becomes apparent that neck mobility was limited in all directions. The maximum angle of total lateral deflection in the neck was 67°. The maximum angle of total ventral deflection was 148° and of total dorsal deflection was 157°. This raises the question of the function of such a long, multi-segment but immobile neck. We posit that the long neck served in hydrodynamic and visual camouflage, hiding the bulk of the body from the small but abundant prey, such as schooling fish and squid. Neck immobility may have been advantageous in withstanding strong hydrodynamic forces acting on the neck during predatory strikes.