Intradiscal pressure, shear strain, and fiber strain in the intervertebral disc under combined loading

The Impact of Demographic Characteristics, Obesity, Surgical Level, and Intervertebral Disc Properties on Recurrence of Lumbar Disc Herniation

BACKGROUND

This study focused on identifying factors influencing recurrent lumbar disc herniation (RLDH) by analyzing demographic data, body mass index (BMI), and radiologic disc properties in patients undergoing single-level unilateral lumbar disc herniation surgery.

METHODS

Our retrospective analysis included 2 groups from our clinic: Group 1 (n = 41) with patients experiencing RLDH requiring a second surgery, and Group 2 (n = 73) with patients having a single surgery and no recurrence over at least a 5-year follow-up. We assessed age, sex, diabetes mellitus (DM), hypertension (HT), BMI, Pfirrmann disc degeneration type, herniation types (protrusion, extrusion, and sequestration), and surgical level.

RESULTS

The mean ages of Groups 1 and 2 were 48.93 ± 13.47 and 44.4 ± 11.79 years, respectively, with no significant age difference (P = 0.064). Gender distribution was also not significantly different, with 63.41% males in Group 1 and 56.16% in Group 2 (P = 0.450). DM prevalence was similar in both groups (P = 0.727). Notably, HT was significantly lower in Group 2 (P = 0.018). The average BMI was comparable between groups (P = 0.607), and no significant difference in Pfirrmann disc degeneration scores was observed (P = 0.547). Radiologic disc type distributions did not significantly differ (P = 0.448). Most surgeries in both groups were at the L4-5 level, with no significant differences in surgical levels (P = 0.456).

CONCLUSIONS

We found that factors like gender, age, DM, obesity, surgical level, disc degeneration, and disc types do not significantly impact RLDH. However, the higher occurrence of HT in recurrent cases indicates a potential area for further research.

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.

A new method to design energy-conserving surrogate models for the coupled, nonlinear responses of intervertebral discs

The aim of this study was to design physics-preserving and precise surrogate models of the nonlinear elastic behaviour of an intervertebral disc (IVD). Based on artificial force-displacement data sets from detailed finite element (FE) disc models, we used greedy kernel and polynomial approximations of second, third and fourth order to train surrogate models for the scalar force-torque -potential. Doing so, the resulting models of the elastic IVD responses ensured the conservation of mechanical energy through their structure. At the same time, they were capable of predicting disc forces in a physiological range of motion and for the coupling of all six degrees of freedom of an intervertebral joint. The performance of all surrogate models for a subject-specific L4 | 5 disc geometry was evaluated both on training and test data obtained from uncoupled (one-dimensional), weakly coupled (two-dimensional), and random movement trajectories in the entire six-dimensional (6d) physiological displacement range, as well as on synthetic kinematic data. We observed highest precisions for the kernel surrogate followed by the fourth-order polynomial model. Both clearly outperformed the second-order polynomial model which is equivalent to the commonly used stiffness matrix in neuro-musculoskeletal simulations. Hence, the proposed model architectures have the potential to improve the accuracy and, therewith, validity of load predictions in neuro-musculoskeletal spine models.

Insights into the mechanical microenvironment within the cartilaginous endplate: An emerging role in maintaining disc homeostasis and normal function

Biomechanical factors are strongly linked with the emergence and development of intervertebral disc degeneration (IVDD). The intervertebral disc (IVD), as a unique enclosed biomechanical structure, exhibits distinct mechanical properties within its substructures. Damage to the mechanical performance of any substructure can disrupt the overall mechanical function of the IVD. Endplate degeneration serves as a significant precursor to IVDD. The endplate (EP) structure, especially the cartilaginous endplate (CEP), serves as a conduit for nutrient and metabolite transport in the IVD. It is inevitably influenced by its nutritional environment, mechanical loading, cytokines and extracellular components. Currently, reports on strategies targeting the CEP for the prevention and treatment of IVDD are scarce. This is due to two primary reasons: first, limited knowledge of the biomechanical microenvironment surrounding the degenerated CEP cells; and second, innovative biological treatment strategies, such as implanting active cells (disc or mesenchymal stem cells) or modulating natural cell activity through the addition of therapeutic factors or genes to treat IVDD often overlook a critical aspect-the restoration of the nutrient supply function and mechanical microenvironment of the endplate. Therefore, restoring the healthy structure of the CEP and maintaining a stable mechanical microenvironment within the EP are crucial for the prevention of IVDD and the repair of degenerated IVDs. We present a comprehensive literature review on the mechanical microenvironment characteristics of cartilage endplates and their associated mechanical signaling pathways. Our aim is to provide valuable insights into the development and implementation of strategies to prevent IVDD by delaying or reversing CEP degeneration.

Stepwise reduction of bony density in patients induces a higher risk of annular tears by deteriorating the local biomechanical environment

BACKGROUND CONTEXT

The relationship between osteoporosis and intervertebral disc degeneration (IDD) remains unclear. Considering that annular tear is the primary phenotype of IDD in the lumbar spine, the deteriorating local biomechanical environment may be the main trigger for annular tears.

PURPOSE

To investigate whether poor bone mineral density (BMD) in the vertebral bodies may increase the risk of annular tears via the degradation of the local biomechanical environment.

STUDY DESIGN

This study was a retrospective investigation with relevant numerical mechanical simulations.

PATIENT SAMPLE

A total of 64 patients with low back pain (LBP) and the most severe IDD in the L4-L5 motion segment were enrolled.

OUTCOME MEASURES

Annulus integration status was assessed using diffusion tensor fibre tractography (DTT). Hounsfield unit (HU) values of adjacent vertebral bodies were employed to determine BMD. Numerical simulations were conducted to compute stress values in the annulus of models with different BMDs and body positions.

METHODS

The clinical data of the 64 patients with low back pain were collected retrospectively. The BMD of the vertebral bodies was measured using the HU values, and the annulus integration status was determined according to DTT. The data of the patients with and without annular tears were compared, and regression analysis was used to identify the independent risk factors for annular tears. Furthermore, finite element models of the L4-L5 motion segment were constructed and validated, followed by estimating the maximum stress on the post and postlateral interfaces between the superior and inferior bony endplates (BEPs) and the annulus.

RESULTS

Patients with lower HU values in their vertebral bodies had significantly higher incidence rates of annular tears, with decreased HU values being an independent risk factor for annular tears. Moreover, increased stress on the BEP-annulus interfaces was associated with a stepwise reduction of bony density (ie, elastic modulus) in the numerical models.

CONCLUSIONS

The stepwise reduction of bony density in patients results in a higher risk of annular tears by deteriorating the local biomechanical environment. Thus, osteoporosis should be considered to be a potential risk factor for IDD biomechanically.

Senescent-like macrophages mediate angiogenesis for endplate sclerosis via IL-10 secretion in male mice

Endplate sclerosis is a notable aspect of spine degeneration or aging, but the mechanisms remain unclear. Here, we report that senescent macrophages accumulate in the sclerotic endplates of lumbar spine instability (LSI) or aging male mouse model. Specifically, knockout of cdkn2a (p16) in macrophages abrogates LSI or aging-induced angiogenesis and sclerosis in the endplates. Furthermore, both in vivo and in vitro studies indicate that IL-10 is the primary elevated cytokine of senescence-related secretory phenotype (SASP). Mechanistically, IL-10 increases pSTAT3 in endothelial cells, leading to pSTAT3 directly binding to the promoters of Vegfa, Mmp2, and Pdgfb to encourage their production, resulting in angiogenesis. This study provides information on understanding the link between immune senescence and endplate sclerosis, which might be useful for therapeutic approaches.

Simulation of cancer progression in bone by a virtual thermal flux with a case study on lumbar vertebrae with multiple myeloma

BACKGROUND

Two main problems examining the mechanism of cancer progression in the tissues using the computational models are lack of enough knowledge on the effective factors for such events in vivo environments and lack of specific parameters in the available computational models to simulate such complicated reactions.

METHODS

In this study, it was tried to simulate the progression of cancerous lesions in the bone tissues by an independent parameter from the anatomical and physiological characteristics of the tissues, so to degrade the orthotropic mechanical properties of the bone tissues, a virtual temperature was determined to be used by a well-known framework for simulation of damages in the composite materials. First, the reliability of the FE model to simulate hyperelastic response in the intervertebral discs (IVDs) and progressive failure in the bony components were verified by simulation of some In-Vitro tests, available in the literature. Then, the progression of the osteolytic damage was simulated in a clinical case with multiple myeloma in the lumbar vertebrae.

RESULTS

The FE model could simulate stress-shielding and diffusion of lesion in the posterior elements of the damaged vertebra which led to spinal stenosis. The load carrying shares associated with the anterior half and the posterior half of the examined vertebral body and the posterior elements were estimated equal to 41 %, 47 % and 12 %, respectively for the intact condition, that changed to 14 %, 16 % and 70 %, when lesion occupied one third of the vertebral body.

CONCLUSION

Correlation of the FE results with the deformation shapes, observed in the MRIs for the clinical case study, indicated appropriateness of the procedure, proposed for simulation of the progressive osteolytic damage in the vertebral segments. The future studies may follow simulation of tumor growth for various metastatic tissues using the method, established here.

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.

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.

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.

Bibliometric and Visualization Analysis of Biomechanical Research on Lumbar Intervertebral Disc

Background

Biomechanical research on the lumbar intervertebral disc (IVD) provides valuable information for the diagnosis, treatment, and prevention of related diseases, and has received increasing attention. Using bibliometric methods and visualization techniques, this study investigates for the first time the research status and development trends in this field, with the aim of providing guidance and support for subsequent research.

Methods

The Science Citation Index Expanded (SCI-Expanded) within the Web of Science Core Collection (WoSCC) database was used as the data source to select literature published from 2003 to 2022 related to biomechanical research on lumbar IVD. VOSviewer 1.6.19 and CiteSpace 6.2.R2 visualization software, as well as the online analysis platform of literature metrology, were utilized to generate scientific knowledge maps for visual display and data analysis.

Results

The United States is the most productive country in this field, with the Ulm University making the largest contribution. Wilke HJ is both the most prolific author and one of the highly cited authors, while Adams MA is the most cited author. Spine, J Biomech, Eur Spine J, Spine J, and Clin Biomech are not only the journals with the highest number of publications, but also highly cited journals. The main research topics in this field include constructing and validating three-dimensional (3D) finite element model (FEM) of lumbar spine, measuring intradiscal pressure, exploring the biomechanical effects and related risk factors of lumbar disc degeneration, studying the mechanical responses to different torque load combinations, and classifying lumbar disc degeneration based on magnetic resonance images (MRI), which are also the hot research themes in recent years.

Conclusion

This study systematically reviews the knowledge system and development trends in the field of biomechanics of lumbar IVD, providing valuable references for further research.

Ligamentous tethering and intradiscal pressure affecting the mechanical environment of scoliotic spines

Despite several theories have been proposed to explain the progression of Adolescent Idiopathic Scoliosis (AIS), there is no consensus on the mechanical factors that control the spinal deformities. Prominent biomechanical notions focus on the geometrical asymmetry and differential growth, however, the correlation between these phenomena remains unclear. We postulate that intradiscal pressure and its connection with the supporting ligamentous structures are the reasons behind the asymmetric growth in AIS. To investigate this hypothesis, a numerical 3D patient-specific model of a scoliotic spine is constructed to carry upper body weight. Four analyses are performed: control simulation with no ligaments followed by 3 simulations, in each, a different and stiffer set of ligaments is employed. The analyses showed that intradiscal pressure is relatively high in the spine's higher-deformity region. Moreover, the stiffness effect of the ligamentous tethering correlated directly to intradiscal pressure; the stiffer the ligaments, the higher the intradiscal pressure. Due to geometrical asymmetry, the pressure is eccentric toward the concave region of deformed vertebral units. As a result, the deformed annulus fibrosus generated uplifts in the convex side of deformed vertebral units. The eccentric pressure and the uplift are opposite in location and direction creating an imbalanced mechanical environment for the spine during growth.

The influence of over-distraction on biomechanical response of cervical spine post anterior interbody fusion: a comprehensive finite element study

Introduction: Anterior cervical discectomy and fusion (ACDF) has been considered as the gold standard surgical treatment for cervical degenerative pathologies. Some surgeons tend to use larger-sized interbody cages during ACDF to restore the index intervertebral disc height, hence, this study evaluated the effect of larger-sized interbody cages on the cervical spine with ACDF under both static and cyclic loading. Method: Twenty pre-operative personalized poro-hyperelastic finite element (FE) models were developed. ACDF post-operative models were then constructed and four clinical scenarios (i.e., 1) No-distraction; 2) 1 mm distraction; 3) 2 mm distraction; and 4) 3 mm distraction) were predicted for each patient. The biomechanical responses at adjacent spinal levels were studied subject to static and cyclic loading. Non-parametric Friedman statistical comparative tests were performed and the p values less than 0.05 were reflected as significant. Results: The calculated intersegmental range of motion (ROM) and intradiscal pressure (IDP) from 20 pre-operative FE models were within the overall ranges compared to the available data from literature. Under static loading, greater ROM, IDP, facet joint force (FJF) values were detected post ACDF, as compared with pre-op. Over-distraction induced significantly higher IDP and FJF in both upper and lower adjacent levels in extension. Higher annulus fibrosus stress and strain values, and increased disc height and fluid loss at the adjacent levels were observed in ACDF group which significantly increased for over-distraction groups. Discussion: it was concluded that using larger-sized interbody cages (the height of ≥2 mm of the index disc height) can result in remarkable variations in biomechanical responses of adjacent levels, which may indicate as risk factor for adjacent segment disease. The results of this comprehensive FE investigation using personalized modeling technique highlight the importance of selecting the appropriate height of interbody cage in ACDF surgery.

Development and validation of lumbar spine finite element model

The functional biomechanics of the lumbar spine have been better understood by finite element method (FEM) simulations. However, there are still areas where the behavior of soft tissues can be better modeled or described in a different way. The purpose of this research is to develop and validate a lumbar spine section intended for biomechanical research. A FE model of the 50th percentile adult male (AM) Total Human Model for Safety (THUMS) v6.1 was used to implement the modifications. The main modifications were to apply orthotropic material properties and nonlinear stress-strain behavior for ligaments, hyperelastic material properties for annulus fibrosus and nucleus pulposus, and the specific content of collagenous fibers in the annulus fibrosus ground substance. Additionally, a separation of the nucleus pulposus from surrounding bones and tissues was implemented. The FE model was subjected to different loading modes, in which intervertebral rotations and disc pressures were calculated. Loading modes contained different forces and moments acting on the lumbar section: axial forces (compression and tension), shear forces, pure moments, and combined loading modes of axial forces and pure moments. The obtained ranges of motion from the modified numerical model agreed with experimental data for all loading modes. Moreover, intradiscal pressure validation for the modified model presented a good agreement with the data available from the literature. This study demonstrated the modifications of the THUMS v6.1 model and validated the obtained numerical results with existing literature in the sub-injurious range. By applying the proposed changes, it is possible to better model the behavior of the human lumbar section under various loads and moments.

Long-term whole-body vibration induces degeneration of intervertebral disc and facet joint in a bipedal mouse model

Background: Whole body vibration (WBV) has been used to treat various musculoskeletal diseases in recent years. However, there is limited knowledge about its effects on the lumbar segments in upright posture mice. This study was performed to investigate the effects of axial Whole body vibration on the intervertebral disc (IVD) and facet joint (FJ) in a novel bipedal mouse model.

Methods: Six-week-old male mice were divided into control, bipedal, and bipedal + vibration groups. Taking advantage of the hydrophobia of mice, mice in the bipedal and bipedal + vibration groups were placed in a limited water container and were thus built standing posture for a long time. The standing posture was conducted twice a day for a total of 6 hours per day, 7 days per week. Whole body vibration was conducted during the first stage of bipedal building for 30 min per day (45 Hz with peak acceleration at 0.3 g). The mice of the control group were placed in a water-free container. At the 10th-week after experimentation, intervertebral disc and facet joint were examined by micro-computed tomography (micro-CT), histologic staining, and immunohistochemistry (IHC), and gene expression was quantified using real-time polymerase chain reaction. Further, a finite element (FE) model was built based on the micro-CT, and dynamic Whole body vibration was loaded on the spine model at 10, 20, and 45 Hz.

Results: Following 10 weeks of model building, intervertebral disc showed histological markers of degeneration, such as disorders of annulus fibrosus and increased cell death. Catabolism genes’ expression, such as Mmp13, and Adamts 4/5, were enhanced in the bipedal groups, and Whole body vibration promoted these catabolism genes’ expression. Examination of the facet joint after 10 weeks of bipedal with/without Whole body vibration loading revealed rough surface and hypertrophic changes at the facet joint cartilage resembling osteoarthritis. Moreover, immunohistochemistry results demonstrated that the protein level of hypertrophic markers (Mmp13 and Collagen X) were increased by long-durationstanding posture, and Whole body vibration also accelerated the degenerative changes of facet joint induced by bipedal postures. No changes in the anabolism of intervertebral disc and facet joint were observed in the present study. Furthermore, finite element analysis revealed that a larger frequency of Whole body vibration loading conditions induced higher Von Mises stresses on intervertebral disc, contact force, and displacement on facet joint.

Conclusion: The present study revealed significant damage effects of Whole body vibration on intervertebral disc and facet joint in a bipedal mouse model. These findings suggested the need for further studies of the effects of Whole body vibration on lumbar segments of humans.

Biomechanical consequences of the intervertebral disc centre of rotation kinematics during lateral bending and axial rotation

The location of the instantaneous centre of rotation (ICR) of a lumbar unit has a considerable clinical importance as a spinal health estimator. Consequently, many studies have been conducted to measure or estimate the ICR during rotations in the three anatomical planes; however the results reported are widely scattered. Even if some inter-subjects variability is to be expected, such inconsistencies are likely explained by the differences in methods and experiments. Therefore, in this paper we seek to model three behaviours of the ICR during lateral bending and axial rotation based on results published in the literature. In order to assess the metabolic and mechanical sensibility to the assumption made on the ICR kinematics, we used a previously validated three dimensional non-linear poroelastic model of a porcine intervertebral disc to simulate physiological lateral and axial rotations. The impact of the geometry was also briefly investigated by considering a 11[Formula: see text] wedge angle. From our simulations, it appears that the hypothesis made on the ICR location does not significantly affect the critical nutrients concentrations but gives disparate predictions of the intradiscal pressure at the centre of the disc (variation up to 0.7 MPa) and of the displacement fields (variation up to 0.4 mm). On the contrary, the wedge angle does not influence the estimated intradiscal pressure but leads to minimal oxygen concentration decreased up to 33% and increased maximal lactate concentration up to 13%. While we can not settle on which definition of the ICR is more accurate, this work suggests that patient-specific modeling of the ICR is required and brings new insights that can be useful for the development of new tools or the design of surgical material such as total lumbar disc prostheses.

Anatomical parameters alter the biomechanical responses of adjacent segments following lumbar fusion surgery: Personalized poroelastic finite element modelling investigations

Introduction: While the short-term post-operative outcome of lumbar fusion is satisfying for most patients, adjacent segment disease (ASD) can be prevalent in long-term clinical observations. It might be valuable to investigate if inherent geometrical differences among patients can significantly alter the biomechanics of adjacent levels post-surgery. This study aimed to utilize a validated geometrically personalized poroelastic finite element (FE) modeling technique to evaluate the alteration of biomechanical response in adjacent segments post-fusion. Methods: Thirty patients were categorized for evaluation in this study into two distinct groups [i.e., 1) non-ASD and 2) ASD patients] based on other long-term clinical follow-up investigations. To evaluate the time-dependent responses of the models subjected to cyclic loading, a daily cyclic loading scenario was applied to the FE models. Different rotational movements in different planes were superimposed using a 10 Nm moment after daily loading to compare the rotational motions with those at the beginning of cyclic loading. The biomechanical responses of the lumbosacral FE spine models in both groups were analyzed and compared before and after daily loading. Results: The achieved comparative errors between the FE results and clinical images were on average below 20% and 25% for pre-op and post-op models, respectively, which confirms the applicability of this predictive algorithm for rough pre-planning estimations. The results showed that the disc height loss and fluid loss were increased for the adjacent discs in post-op models after 16 h of cyclic loading. In addition, significant differences in disc height loss and fluid loss were observed between the patients who were in the non-ASD and ASD groups. Similarly, the increased stress and fiber strain in the annulus fibrosus (AF) was higher in the adjacent level of post-op models. However, the calculated stress and fiber strain values were significantly higher for patients with ASD. Discussion: Evaluating the biomechanical response of pre-op and post-op modeling in the non-ASD and ASD groups showed that the inherent geometric differences among patients cause significant variations in the estimated mechanical response. In conclusion, the results of the current study highlighted the effect of geometrical parameters (which may refer to the anatomical conditions or the induced modifications regarding surgical techniques) on time-dependent responses of lumbar spine biomechanics.

Mechanical stimulation promotes MSCs healing the lesion of intervertebral disc annulus fibrosus

Mesenchymal stem cells (MSCs) and scaffolds offer promising perspectives for annulus fibrosus (AF) repair. The repair effect was linked to features of the local mechanical environment related to the differentiation of MSCs. In this study, we established a Fibrinogen-Thrombin-Genipin (Fib-T-G) gel which is sticky and could transfer strain force from AF tissue to the human mesenchymal stem cells (hMSCs) embedded in the gel. After the Fib-T-G biological gel was injected into the AF fissures, the histology scores of intervertebral disc (IVD) and AF tissue showed that Fib-T-G gel could better repair the AF fissure in caudal IVD of rats, and increase the expression of AF-related proteins including Collagen 1 (COL1), Collagen 2 (COL2) as well as mechanotransduction-related proteins including RhoA and ROCK1. To clarify the mechanism that sticky Fib-T-G gel induces the healing of AF fissures and the differentiation of hMSCs, we further investigated the differentiation of hMSCs under mechanical strain in vitro. It was demonstrated that both AF-specific genes, including Mohawk and SOX-9, and ECM markers (COL1, COL2, aggrecan) of hMSCs were up-regulated in the environment of strain force. Moreover, RhoA/ROCK1 proteins were also found to be significantly up-regulated. In addition, we further -demonstrated that the fibrochondroinductive effect of the mechanical microenvironment process could be significantly blocked or up-regulated by inhibiting the RhoA/ROCK1 pathway or overexpressing RhoA in MSCs, respectively. Summarily, this study will provide a therapeutic alternative to repair AF tears and provide evidence that RhoA/ROCK1 is vital for hMSCs response to mechanical strain and AF-like differentiation.

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.

Comparative biomechanical analyses of lower cervical spine post anterior fusion versus intervertebral disc arthroplasty: A geometrically patient-specific poroelastic finite element investigation

Background/Objective

The optimal surgical technique for the treatment of cervical degenerative disc disease (CDDD) towards decreasing the risk of adjacent segment disease (ASD) remains elusive. This study aimed to comparatively investigate the biomechanics of the lower cervical spine following fusion (ACDF) and artificial disc arthroplasty (Bryan® and Prestige LP®) using a validated geometrically patient-specific poroelastic finite element modeling (FEM) approach.

Methods

Ten subject-specific pre-operative models were developed and validated based on a FEM approach. Poroelastic models were then constructed using post-operation images for three different treatment scenarios: ACDF; Prestige LP® and Bryan® artificial discs at the C5-C6 level. The biomechanical responses at both surgical and adjacent spinal levels were studied subject to static and cyclic loading.

Results

Postoperatively, greater range of motion (ROM), higher annulus fibrosus stress and strain values, and increased disc height and fluid loss at the adjacent levels were detected post ACDF, as compared with pre-op as well as artificial disc arthroplasty. The facet joint forces were larger for the Prestige LP® disc, particularly during extension. The lowest values in disc height and fluid exchange were observed in the Bryan® artificial disc arthroplasty models.

Conclusion

Biomechanical analyses revealed that ACDF poses the highest potential risk for adjacent disc degeneration. The artificial discs investigated here (Prestige LP® and Bryan®) not only preserved motion at the instrumented level, but also sustained the pre-op ROM and decreased the intradiscal pressure (IDP) and facet joint forces (FJFs) at adjacent levels, particularly during flexion/extension. The Bryan® artificial disc demonstrated the most efficacy in maintaining the natural poroelastic characteristics of adjacent discs.

The translational potential of this article

This study provided a technique for clinicians to use quantitative data towards subject-specific evaluation to comparatively evaluate the impact of ACDF and disc arthroplasty using two types of artificial discs on the biomechanics of the cervical spine. It confirms differences in the poroelastic characteristics of adjacent discs for different fixation techniques, and reveals the advantage of artificial discs with a flexible core for decreasing the risk of ASD.

The influence of the viscoelastic property of polycarbonate urethane as an artificial disc core material under various physiological motions at the L4-L5 level

Intervertebral disc (IVD) degeneration is one of the musculoskeletal disorders due to the Degenerative Disc Disease (DDD), that cause low back pain (LBP) and leads to a reduced range of motion. Spinal fusion and arthroplasty are the other surgical procedures that could replace the disc affected by DDD against artificial disc replacement (ADR). This study aims to analyse the biomechanical behaviour of proposed core material as Polycarbonate Urethane (PCU) in the L4-L5 lumbar segment for ADR with Ti-6Al-4V and Co-28Cr-6M as endplate materials and compare it to the performance of an ultra-high molecular weight polyethylene (UHMWPE) core. Finite element methods have been approached to measure the overall stress distribution along with other physiological motions like Flexion (FLEX), Extension (EXT), Axial rotation (AR) and Lateral bending (LB), respectively. Preload of 450 N compressive load, 8 N-m for Flex, 6 N-m for EXT, 6 N-m for AR and 4 N-m for LB are applied. It could be concluded that Ti-6Al-4V - PCU and Co-28Cr-6M - PCU is the best composition for the ADR for the L4-L5 level.

Psoas muscle measurement as a predictor of recurrent lumbar disc herniation: A retrospective blind study

AbstractParaspinal (erector spinae and multifidus) and psoas muscles contribute to spinal stability, but no study has yet examined the relationship between muscle mass and recurrent lumbar disc herniation (rLDH). The purpose of this study was to investigate the effect of psoas and paraspinal muscle mass on recurrent Lumbar disc herniation (LDH). This retrospective study included 49 patients with LDH (22 men, 27 women; mean age: 59.9 years; range 32-80) who underwent discectomy and partial laminectomy without fusion and underwent both pre- and postoperative magnetic resonance imaging. The presence of rLDH was determined using medical records and postoperative magnetic resonance imagings. Patients were divided into an rLDH group (26 patients) and a without-rLDH group (23 patients). Clinical characteristics, segmental motion, and paraspinal and psoas muscle mass were compared between the groups. Using ImageJ software, the cross-sectional area (CSA), lean muscle mass (LMM), and skeletal muscle index (SMI) were measured on T2 axial preoperative magnetic resonance images at L2-L3, L3-L4, and L4-L5 disc levels to represent muscle mass. Univariate and multivariate logistic regression analyses were performed. In the rLDH group, patients were younger (52.6 years vs 68.2 years; P = .001), segmental instability was more common (50.0% vs 4.3%; P = .001), and the CSA, LMM, CSASMI, and LMMSMI of psoas muscles were larger (5851.59 mm2 vs 4264.93 mm2, 5456.59 mm2 vs 4044.77 mm2, 18.77 cm2/m2 vs 13.86 cm2/m2, and 17.52 cm2/m2 vs 12.98 cm2/m2; P < .01 for all 4 variables). On multivariate logistic regression, age and segmental instability were independent risk factors for rLDH (odds ratio 0.886 and 18.527; P = .01 and P = .02, respectively). In middle-aged and elderly patients with lumbar disc herniation, relatively younger age, segmental instability, and greater psoas muscle mass may be risk factors for recurrence.

Biomechanical Properties of Paraspinal Muscles Influence Spinal Loading—A Musculoskeletal Simulation Study

Paraspinal muscles are vital to the functioning of the spine. Changes in muscle physiological cross-sectional area significantly affect spinal loading, but the importance of other muscle biomechanical properties remains unclear. This study explored the changes in spinal loading due to variation in five muscle biomechanical properties: passive stiffness, slack sarcomere length (SSL), in situ sarcomere length, specific tension, and pennation angle. An enhanced version of a musculoskeletal simulation model of the thoracolumbar spine with 210 muscle fascicles was used for this study and its predictions were validated for several tasks and multiple postures. Ranges of physiologically realistic values were selected for all five muscle parameters and their influence on L4-L5 intradiscal pressure (IDP) was investigated in standing and 36° flexion. We observed large changes in IDP due to changes in passive stiffness, SSL, in situ sarcomere length, and specific tension, often with interesting interplays between the parameters. For example, for upright standing, a change in stiffness value from one tenth to 10 times the baseline value increased the IDP only by 91% for the baseline model but by 945% when SSL was 0.4 μm shorter. Shorter SSL values and higher stiffnesses led to the largest increases in IDP. More changes were evident in flexion, as sarcomere lengths were longer in that posture and thus the passive curve is more influential. Our results highlight the importance of the muscle force-length curve and the parameters associated with it and motivate further experimental studies on in vivo measurement of those properties.

How annulus defects can act as initiation sites for herniation

PURPOSE

Both posture and loading rate are key factors in the herniation process and can determine the mechanism of disc failure. The aim of this study was to test the hypothesis that disruption visible with HR-MRI post-testing corresponds with microstructural features and further elucidate the mechanism by which this disruption weakens the disc. This will enable us to gain new insights into the herniation process.

METHODS

Thirty ovine lumbar spinal segments were subjected to combinations of four loading conditions (0-12° flexion, 0-9° lateral bending, 0-4° axial rotation, 0-1500 N axial compression) for 1000 loading cycles at 2 Hz in a dynamic disc loading simulator. The discs were scanned in an ultra-high field MRI (11.7 T) then examined using brightfield microscopy to examine their microstructure.

RESULTS

Four discs herniated and seven discs suffered nucleus displacement. These discs contained pre-existing defects in the central posterior annulus. Generally, following testing discs contained more posterior annulus disruption, Microstructural investigation revealed there was clear correspondence between HR-MRI and microstructural observations, and that the mid-outer annular-endplate junction had failed in all discs examined in this study.

CONCLUSIONS

While all discs suffered outer annulus damage, only the discs containing pre-existing defects herniated. These pre-existing defects weakened the inner and mid annulus, allowing herniation to occur once the mid and outer annular wall was compromised. We propose this can occur during the degenerative cascade.

Strength and Power-Related Measures in Assessing Core Muscle Performance in Sport and Rehabilitation

While force-velocity-power characteristics of resistance exercises, such as bench presses and squats, have been well documented, little attention has been paid to load, force, and power-velocity relationships in exercises engaging core muscles. Given that power produced during lifting tasks or trunk rotations plays an important role in most sport-specific and daily life activities, its measurement should represent an important part of the test battery in both athletes and the general population. The aim of this scoping review was 1) to map the literature related to testing methods assessing core muscle strength and stability in sport and rehabilitation, chiefly studies with particular focus on force-velocity-power characteristics of exercises involving the use of core muscles, 2) and to identify gaps in existing studies and suggest further research in this field. The literature search was conducted on Cochrane Library databases, Scopus, Web of Science, PubMed and MEDLINE, which was completed by SpringerLink, Google Scholar and Elsevier. The inclusion criteria were met in 37 articles. Results revealed that among a variety of studies investigating the core stability and core strength in sport and rehabilitation, only few of them analyzed force-velocity-power characteristics of exercises involving the use of core muscles. Most of them evaluated maximal isometric strength of the core and its endurance. However, there are some studies that assessed muscle power during lifting tasks at different loads performed either with free weights or using the Smith machine. Similarly, power and velocity were assessed during trunk rotations performed with different weights when standing or sitting. Nevertheless, there is still scant research investigating the power-velocity and force-velocity relationship during exercises engaging core muscles in able-bodied and para athletes with different demands on stability and strength of the core. Therefore, more research is needed to address this gap in the literature and aim research at assessing strength and power-related measures within cross-sectional and intervention studies. A better understanding of the power-force-velocity profiles during exercises with high demands on the core musculature has implications for designing sport training and rehabilitation programs for enhancement of athletes' performance and/or decrease their risk of back pain.

The Mismatch Between Bony Endplates and Grafted Bone Increases Screw Loosening Risk for OLIF Patients With ALSR Fixation Biomechanically

The mismatch between bony endplates (BEPs) and grafted bone (GB) triggers several complications biomechanically. However, no published study has identified whether this factor increases the risk of screw loosening by deteriorating the local stress levels. This study aimed to illustrate the biomechanical effects of the mismatch between BEP and GB and the related risk of screw loosening. In this study, radiographic and demographic data of 56 patients treated by single segment oblique lumbar interbody fusion (OLIF) with anterior lateral single rod (ALSR) fixation were collected retrospectively, and the match sufficiency between BEP and GB was measured and presented as the grafted bony occupancy rate (GBOR). Data in patients with and without screw loosening were compared; regression analyses identified independent risk factors. OLIF with different GBORs was simulated in a previously constructed and validated lumbosacral model, and biomechanical indicators related to screw loosening were computed in surgical models. The radiographic review and numerical simulations showed that the coronal plane’s GBOR was significantly lower in screw loosening patients both in the cranial and caudal vertebral bodies; the decrease in the coronal plane’s GBOR has been proven to be an independent risk factor for screw loosening. In addition, numerical mechanical simulations showed that the poor match between BEP and GB will lead to stress concentration on both screws and bone-screw interfaces. Therefore, we can conclude that the mismatch between the BEP and GB will increase the risk of screw loosening by deteriorating local stress levels, and the increase in the GBOR by modifying the OLIF cage’s design may be an effective method to optimize the patient’s prognosis.

Optimization of a lumbar interspinous fixation device for the lumbar spine with degenerative disc disease

Interspinous spacer devices used in interspinous fixation surgery remove soft tissues in the lumbar spine, such as ligaments and muscles and may cause degenerative diseases in adjacent segments its stiffness is higher than that of the lumbar spine. Therefore, this study aimed to structurally and kinematically optimize a lumbar interspinous fixation device (LIFD) using a full lumbar finite element model that allows for minimally invasive surgery, after which the normal behavior of the lumbar spine is not affected. The proposed healthy and degenerative lumbar spine models reflect the physiological characteristics of the lumbar spine in the human body. The optimum number of spring turns and spring wire diameter in the LIFD were selected as 3 mm and 2 turns, respectively—from a dynamic range of motion (ROM) perspective rather than a structural maximum stress perspective—by applying a 7.5 N∙m extension moment and 500 N follower load to the LIFD-inserted lumbar spine model. As the spring wire diameter in the LIFD increased, the maximum stress generated in the LIFD increased, and the ROM decreased. Further, as the number of spring turns decreased, both the maximum stress and ROM of the LIFD increased. When the optimized LIFD was inserted into a degenerative lumbar spine model with a degenerative disc, the facet joint force of the L3-L4 lumbar segment was reduced by 56%–98% in extension, lateral bending, and axial rotation. These results suggest that the optimized device can strengthen the stability of the lumbar spine that has undergone interspinous fixation surgery and reduce the risk of degenerative diseases at the adjacent lumbar segments.

Investigation of geometric deformations of the lumbar disc during axial body rotations

Background

Quantitative data on in vivo vertebral disc deformations are critical for enhancing our understanding of spinal pathology and improving the design of surgical materials. This study investigated in vivo lumbar intervertebral disc deformations during axial rotations under different load-bearing conditions.

Methods

Twelve healthy subjects (7 males and 5 females) between the ages of 25 and 39 were recruited. Using a combination of a dual fluoroscopic imaging system (DFIS) and CT, the images of L3–5 segments scanned by CT were transformed into three-dimensional models, which matched the instantaneous images of the lumbar spine taken by a double fluorescent X-ray system during axial rotations to reproduce motions. Then, the kinematic data of the compression and shear deformations of the lumbar disc and the coupled bending of the vertebral body were obtained.

Results

Relative to the supine position, the average compression deformation caused by rotation is between + 10% and − 40%, and the shear deformation is between 17 and 50%. Under physiological weightbearing loads, different levels of lumbar discs exhibit similar deformation patterns, and the deformation patterns of left and right rotations are approximately symmetrical. The deformation patterns change significantly under a 10 kg load, with the exception of the L3–4 disc during the right rotation.

Conclusion

The deformation of the lumbar disc was direction-specific and level-specific during axial rotations and was affected by extra weight. These data can provide new insights into the biomechanics of the lumbar spine and optimize the parameters of artificial lumbar spine devices.

Force Distribution Within Spinal Tissues During Posterior to Anterior Spinal Manipulative Therapy: A Secondary Analysis

Background

Previous studies observed that the intervertebral disc experiences the greatest forces during spinal manipulative therapy (SMT) and that the distribution of forces among spinal tissues changes as a function of the SMT parameters. However, contextualized SMT forces, relative to the ones applied to and experienced by the whole functional spinal unit, is needed to understand SMT’s underlying mechanisms.

Aim

To describe the percentage force distribution between spinal tissues relative to the applied SMT forces and total force experienced by the functional unit.

Methods

This secondary analysis combined data from 35 fresh porcine cadavers exposed to a simulated 300N SMT to the skin overlying the L3/L4 facet joint via servo-controlled linear motor actuator. Vertebral kinematics were tracked optically using indwelling bone pins. The functional spinal unit was then removed and mounted on a parallel robotic platform equipped with a 6-axis load cell. The kinematics of the spine during SMT were replayed by the robotic platform. By using serial dissection, peak and mean forces induced by the simulated SMT experienced by spinal structures in all three axes of motion were recorded. Forces experienced by spinal structures were analyzed descriptively and the resultant force magnitude was calculated.

Results

During SMT, the functional spinal unit experienced a median peak resultant force of 36.4N (IQR: 14.1N) and a mean resultant force of 25.4N (IQR: 11.9N). Peak resultant force experienced by the spinal segment corresponded to 12.1% of the total applied SMT force (300N). When the resultant force experienced by the functional spinal unit was considered to be 100%, the supra and interspinous ligaments experienced 0.3% of the peak forces and 0.5% of the mean forces. Facet joints and ligamentum flavum experienced 0.7% of the peak forces and 3% of the mean forces. Intervertebral disc and longitudinal ligaments experienced 99% of the peak and 96.5% of the mean forces.

Conclusion

In this animal model, a small percentage of the forces applied during a posterior-to-anterior SMT reached spinal structures in the lumbar spine. Most SMT forces (over 96%) are experienced by the intervertebral disc. This study provides a novel perspective on SMT force distribution within spinal tissues.

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.

A novel approach for tetrahedral-element-based finite element simulations of anisotropic hyperelastic intervertebral disc behavior

Intervertebral discs are microstructurally complex spinal tissues that add greatly to the flexibility and mechanical strength of the human spine. Attempting to provide an adjustable basis for capturing a wide range of mechanical characteristics and to better address known challenges of numerical modeling of the disc, we present a robust finite-element-based model formulation for spinal segments in a hyperelastic framework using tetrahedral elements. We evaluate the model stability and accuracy using numerical simulations, with particular attention to the degenerated intervertebral discs and their likely skewed and narrowed geometry. To this end, 1) annulus fibrosus is modeled as a fiber-reinforced Mooney-Rivlin type solid for numerical analysis. 2) An adaptive state-variable dependent explicit time step is proposed and utilized here as a computationally efficient alternative to theoretical estimates. 3) Tetrahedral-element-based FE models for spinal segments under various loading conditions are evaluated for their use in robust numerical simulations. For flexion, extension, lateral bending, and axial rotation load cases, numerical simulations reveal that a suitable framework based on tetrahedral elements can provide greater stability and flexibility concerning geometrical meshing over commonly employed hexahedral-element-based ones for representation and study of spinal segments in various stages of degeneration.

Adjuvant surgical decision-making system for lumbar intervertebral disc herniation after percutaneous endoscopic lumber discectomy: a retrospective nonlinear multiple logistic regression prediction model based on a large sample

BACKGROUND CONTEXT

Lumbar disc herniation (LDH) is a common condition that can affects an individual' quality of life. In patients for whom conservative treatment is ineffective after 3 months, surgical treatment, such as percutaneous endoscopic lumbar discectomy (PELD), is recommended. Because PELD is minimally invasive and produces thorough nerve root decompression, both surgeons and patients often prefer it to other techniques.

PURPOSE

Surgeons find it challenging to prevent postoperative recurrent LDH (rLDH) when they use PELD. We created and verified a model for evaluating patients' recurrence risk factors before surgery so that surgeons can choose other surgical techniques when necessary.

STUDY DESIGN

Retrospective study.

PATIENT SAMPLE

One thousand eight hundred seven patients who underwent PELD at our hospital between 2012 and 2015 were enrolled.

OUTCOME MEASURE

The main outcome measure was rLDH at any follow-up time point.

METHODS

Data were retrospectively analyzed for 1807 patients who underwent PELD at our hospital at some point between 2012 and 2015; all patients had been monitored for at least 5 years after surgery. They were divided into a recurrence group and a nonrecurrence group. Clinical and radiological risk factors were assessed over time to determine their correlations with recurrence and to exclude less important factors. A nonlinear multivariate logistic regression model was established to predict the recurrence rate before surgery.

RESULTS

A total of 1706 patients were monitored after PELD; data were missing for 101 additional patients. The total recurrence rate was 10.38%, and the most common time from surgery to recurrence was 1 year. Ten risk factors were assessed and included in the analysis. Regarding clinical risk factors, patients with hypertension (p < .001; correlation coefficient R [R] = 0.235; odds ratio [OR] = 4.749), diabetes (p < .001; R = 0.381; OR = 16.797), a history of smoking (p < .001; R = 0.347; OR = 9.012), and a history of performing intense physical labor (p < .001; R = 0.409; OR = 19.592) had a higher recurrence rate. Regarding radiological risk factors, patients with disc degeneration (Pfirrmann grade III) (p < .001; R = 0.228; OR = 4.919), Modic changes (level 2) (p < .001; R = 0.309; OR = 7.934), herniation in the form of extrusion (p < .001; R = 0.365; OR = 12.228), a higher disc height index (DHI) (p < .001; R = 0.336), and a larger segmental range of motion (p < .001; R = 0.243) had a higher recurrence rate. When the lumbar motion angle was negative (p < .001; R = 0.318; OR = 13.680), the recurrence rate was high. The overall accuracy of the final model was 97.6% (1665 of 1706). The recognition rate for non-rLDH cases was 99.0% (1514 of 1529), and the rate for rLDH cases was 85.3% (151 of 177); the AUC was 0.9315. A simple model was used. For those patients with postoperative trauma (p < .001; R = 0.382; OR = 13.680), a comparison model was established, and the corresponding recurrence rate was 23.0% ± 25.0% (0-76%).

CONCLUSIONS

A large cohort of patients underwent long-term monitoring, and 11 risk factors were verified for assessing each patient's risks before surgery to predict the postoperative recurrence of LDH following PELD. The risk of recurrence may be effectively reduced with the use of alternative surgical techniques in high risk cases.

Effect of different attributes of the mimic human lumbar spine biomechanics material structure change by finite element analysis

In this study, we compared stress changes and quantity effect relationships from 3D finite element models of normal and degenerative lumbar segments. We further defined the mechanisms causing alterations in mechanical stability the control of normal and degenerative lumbar segments using traditional Chinese medicine. The characteristics of the stress change and the quantity effect relationships of the three-dimensional finite element model of normal and degenerative lumbar segments were compared. The mechanism(s) leading to changes in mechanical stability and the intervention and balance between normal and degenerative lumbar segments of the traditional Chinese medicine was analyzed. The change trend of stress and strain was compared with the three dimensional finite element model under different motion states of normal lumbar vertebrae. A 3D-FEM of degenerative lumbar segments L4 ~ 5 of the human spine was established to simulate the physiological and pathological changes of the lumbar spine in response to flexion, extension, lateral bending and torsion. The stress changes in the normal and degenerative lumbar vertebrae were assessed through external force interventions and the response to TCM. Stress in the degenerative lumbar vertebrae changed according the external load. Stress and strain were compared in the FEM model under a range of motion states. Components of the human lumbar vertebrae including the cortical vertebrae, cancellous bone, endplates, fibrous rings, and facet articular processes were investigated. The elastic modulus of the nerve roots and the posterior marginal structures of the vertebral body increased with lumbar degeneration. Under stress trends in normal lumbar and different degrees of degenerative lumbar structures including cortical bone, loose bone, terminal plate, fiber ring, nucleus, small articular processes, nerve roots and posterior structures. In normal lumbar spine, 20%, 50%, 70% lumbar degeneration, 106 different lumbar anterior flexion 30 and posterior extension with different external forces showed that ANOVA F was between 3.623 and 11.381 and P changed between 0.001 and 0.05.It is clear that in the lumbar movement segments under different pressure intervention, the changes in the degree of degeneration are significantly different from each constituent structure, among which the trend of expected change between the constituent structures of the lumbar anterior flexion 30 is particularly obvious. The stress distribution in the intervertebral discs were influenced by TCM, and the space in the spinal canal enlarged so that nerve root stress decreased, vertebral body stress increased, and facet processes and pedicle stress in the posterior regions exceeded those of the anterior flexion position. The internal stress of the intervertebral disc increased in the flexion compared to the extension position, gradually increasing from top to bottom. The stress concentration point of the degenerative lumbar disc is significantly greater than the stress in the normal lumbar disc stress distribution area, and increases with the degree of degeneration. Compared with the load capacity of normal lumbar and mild (15% reduction), moderate (40% reduction) lumbar disc protrusion model in bending, extension, axial rotation, lateral bending, the results found that the load transmission of lumbar disc degeneration model to different degrees has also changed, so its compression stiffness, strain distribution and size are also different. TCM can improve and treat lumbar disc disease through its ability to regulate the mechanical environment of degenerative lumbar vertebrae. Compared to the FEM models of the lumbar vertebrae, lumbar degenerative changes could be assessed in response to alterations in the biomechanical environment. These findings provide a scientific basis for the popularization and application of TCM to prevent and treat spinal degenerative disease.

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.

Biomechanics of Lumbar Spine Injury in Road Barrier Collision-Finite Element Study

Literature and field data from CIREN database have shown that lumbar spine injuries occur during car crashes. There are multiple hypotheses regarding how they occur; however, there is no biomechanical explanation for these injuries during collisions with road safety barriers (RSBs). Therefore, the objective of this study was to investigate the mechanics of vertebral fractures during car collisions with concrete RSBs. The finite element method was used for the numerical simulations. The global model of the car collision with the concrete RSB was created. The lumbar spine kinematics were extracted from the global simulation and then applied as boundary conditions to the detailed lumbar spine model. The results showed that during the collision, the occupant was elevated, and then dropped during the vehicle landing. This resulted in axial compression forces 2.6 kN with flexion bending moments 34.7 and 37.8 Nm in the L2 and L3 vertebrae. It was shown that the bending moment is the result of the longitudinal force on the eccentricity. The lumbar spine index for the L1-L5 section was 2.80, thus indicating a lumbar spine fracture. The minimum principal strain criterion of 7.4% and damage variable indicated L2 and L3 vertebrae and the inferior part of L1, as those potentially prone to fracture. This study found that lumbar spine fractures could occur as a consequence of vehicle landing during a collision with a concrete RSB mostly affecting the L1-L3 lumbar spine section. The fracture was caused by a combination of axial forces and flexion bending moments.

Comparative biomechanical analysis of rigid vs. flexible fixation devices for the lumbar spine: A geometrically patient-specific poroelastic finite element study

BACKGROUND AND OBJECTIVE

Lumbar spinal stenosis (LSS), or the narrowing of the spinal canal, continues to be the leading preoperative diagnosis for adults older than 65 years who undergo spine surgery. Although the treatment of LSS depends on its severity, the optimal surgical technique towards decreasing the risk of adjacent segment disease (ASD) remains elusive. This study aimed to comparatively analyze spinal biomechanics with rigid and flexible fixation devices (i.e., rigid and dynamic posterolateral fusion (PLF) and interspinous process (ISP) devices) during daily activities.

METHODS

Using a validated parametric poroelastic finite element modeling approach, 8 subject-specific pre-operative models were developed, and their validity was evaluated. Parametric FE models of the lumbar spines were then regenerated based on post-operation images for (A) rigid PLF (B) dynamic PLF (C) rigid ISP device (Coflex) and (D) flexible ISP device (DIAM) at L4-L5 level. Biomechanical responses for instrumented and adjacent intervertebral discs (IVDs) were analyzed and compared subject to static and cyclic loading.

RESULTS

The preoperative models were well comparable with previous works in literature. The postoperative results for the PLF and Coflex rigid systems, demonstrated greater ROM; higher values of stress and strain in the AF region; and increased disc height and fluid loss at the adjacent levels, as compared with the pre-op models and the post-op results of the flexible systems (i.e., dynamic PLF and DIAM). The calculated forces on the facet joint were of smaller magnitude for the ISP devices as compared to the PLF, particularly during extension.

CONCLUSIONS

This study demonstrates that the dynamic PLF construct and DIAM implants could be effective to maintain the natural poroelastic characteristics of adjacent IVDs, which could be beneficial for enhancing long-term clinical outcomes. FEM provides clinicians with an invaluable patient-specific quantitative tool for informed surgical planning and discerning follow-up management.

Influence of posterior pedicle screw fixation at L4-L5 level on biomechanics of the lumbar spine with and without fusion: a finite element method

Background

Posterior pedicle screw (PS) fixation, a common treatment method for widespread low-back pain problems, has many uncertain aspects including stress concentration levels, effects on adjacent segments, and relationships with physiological motions. A better understanding of how posterior PS fixation affects the biomechanics of the lumbar spine is needed. For this purpose, a finite element (FE) model of a lumbar spine with posterior PS fixation at the L4–L5 segment level was developed by partially removing facet joints (FJs) to imitate an actual surgical procedure. This FE study aimed to investigate the influence of the posterior PS fixation system on the biomechanics of the lumbar spine before and after fusion by determining which physiological motions have the most increase in posterior instrumentation (PI) stresses and FJ loading.

Results

It was determined that posterior PS fixation increased FJ loading by approximately 35% and 23% at the L3–L4 adjacent level with extension and lateral bending motion, respectively. This increase in FJ loading at the adjacent level could point to the possibility that adjacent segment disease has developed or progressed after posterior lumbar interbody fusion. Furthermore, analyses of peak von Mises stresses on PI showed that the maximum PI stresses of 272.1 MPa and 263.7 MPa occurred in lateral bending and flexion motion before fusion, respectively.

Conclusions

The effects of a posterior PS fixation system on the biomechanics of the lumbar spine before and after fusion were investigated for all physiological motions. This model could be used as a fundamental tool for further studies, providing a better understanding of the effects of posterior PS fixation by clearing up uncertain aspects.

Stress analysis of intervertebral disc during occupational activities

BACKGROUND AND OBJECTIVE

Manual material handling activities cause large compression of the intervertebral disc of the lumbar spine. Intradiscal pressure (IDP) has generally been employed to predict the risk of low back injury. As an alternative to in vivo measurements, either motion analysis or finite element (FE) analysis has been used to estimate IDP. The purpose of this study is to propose a new biomechanical method that integrates FE analysis with motion analysis, in order to estimate the stresses and deformations of the intervertebral disc of the lumbar spine during occupational activities.

METHODS

In the proposed method, motion analysis is performed first by using motion capture data, and the results are employed as input data to FE analysis at specific times of interest during motion. In this method, an in-house interface program is used to scale an initial reference FE model to the subject of experiment, and transformed to the corresponding posture at a specific time during motion. The muscle forces and GRF obtained from motion analysis are applied to FE analysis as boundary and loading conditions. For a total of eighteen occupational activities, the IDP, shear stress, and strain of the L4-L5 segment are estimated.

RESULTS

Under each in vivo activity, the predicted IDP was in overall agreement with the available in vivo data. For lifting activities according to lift origin position, the maximum IDP occurred in the far-knee position immediately after lifting. As the lift origin position moved away from the spine, the stresses and strains in the disc increased.

CONCLUSIONS

This new proposed method is expected to allow the estimation of the stresses and deformations in the intervertebral disc during various occupational activities.

Imaging of exercise-induced spinal remodeling in elite rowers

OBJECTIVES

As in-vivo knowledge of training-induced remodeling of intervertebral discs (IVD) is scarce, this study assessed how lumbar IVDs change as a function of long-term training in elite athletes and age-matched controls using compositional Magnetic Resonance Imaging (MRI).

DESIGN

Prospective case-control study.

METHODS

Prospectively, lumbar spines of 17 elite rowers (ERs) of the German national rowing team (mean age: 23.9 ± 3.3 years) were imaged on a clinical 3.0 T MRI scanner. ERs were imaged twice during the annual training cycle, i.e., at training intensive preseason preparations (t₀) and 6 months later during post-competition recovery (t₁). Controls (n = 22, mean age: 26.3 ± 1.9 years) were imaged once at corresponding time points (t₀: n = 11; t₁: n = 11). Segment-wise, the glycosaminoglycan (GAG) content of lumbar IVDs (n = 195) was determined using glycosaminoglycan chemical exchange saturation transfer (gagCEST). Linear mixed models were set up to assess the influence of cohort and other variables on GAG content.

RESULTS

During preseason, IVD GAG values of ERs were significantly higher than those of controls (ERs(t₀): 2.58 ± 0.27% (mean ± standard deviations); controls(t₀): 1.43 ± 0.36%; p ≤ 0.001), while during post-competition recovery, such differences were not present anymore (ERs(t₁): 2.11 ± 0.18%; controls(t₁): 1.89 ± 0.24%; p = 0.362).

CONCLUSIONS

Professional elite-level rowing is transiently associated with significantly higher gagCEST values, which indicate increased lumbar IVD-GAG content and strong remodeling effects in response to training. Beyond professional rowing, core-strengthening full-body exercise may help to enhance the resilience of the lumbar spine as a potential therapeutic target in treating back pain.

Mechanical Parameters and Trajectory of Two Chinese Cervical Manipulations Compared by a Motion Capture System

Objective: To compare the mechanical parameters and trajectory while operating the oblique pulling manipulation and the cervical rotation–traction manipulation.

Methods: An experimental research measuring kinematics parameter and recording motion trajectories of two cervical manipulations were carried out. A total of 48 healthy volunteers participated in this study, who were randomly divided into two groups of 24 representing each of the two manipulations. A clinician performed two manipulations in two groups separately. A motion capture system was used to monitor and analyze kinematics parameters during the operation.

Results: The two cervical manipulations have similar thrust time, displacement, mean velocity, max velocity, and max acceleration. There were no significant differences in active and passive amplitudes between the two cervical rotation manipulations. The thrust amplitudes of the oblique pulling manipulation and the cervical rotation–traction manipulation were 5.735 ± 3.041° and 2.142 ± 1.742°, respectively. The thrust amplitudes of the oblique pulling manipulation was significantly greater than that of the cervical rotation–traction manipulation (P < 0.001).

Conclusion: Compared with the oblique pulling manipulation, the cervical rotation–traction manipulation has a less thrust amplitudes.

A microstructure-based model for a full lamellar-interlamellar displacement and shear strain mapping inside human intervertebral disc core

The determinant role of the annulus fibrosus interlamellar zones in the intervertebral disc transversal and volumetric responses and hence on their corresponding three-dimensional conducts have been only revealed and appreciated recently. Their consideration in disc modeling strategies has been proven to be essential for the reproduction of correct local strain and displacement fields inside the disc especially in the unconstrained directions of the disc. In addition, these zones are known to be the starting areas of annulus fibrosus circumferential tears and disc delamination failure mode, which is often judged as one of the most dangerous disc failure modes that could evolve with time leading to disc hernia. For this latter reason, the main goal of the current contribution is to incorporate physically for the first time, the interlamellar zones, at the scale of a complete human lumbar intervertebral disc, in order to allow a correct local vision and replication of the different lamellar-interlamellar interactions and an identification of the interlamellar critical zones. By means of a fully tridimensional chemo-viscoelastic constitutive model, which we implemented into a finite element code, the physical, mechanical and chemical contribution of the interlamellar zones is added to the disc. The chemical-induced volumetric response is accounted by the model for both the interlamellar zones and the lamellae using experimentally-based fluid kinetics. Computational simulations are performed and critically discussed upon different simple and complex physiological movements. The disc core and the interlamellar zones are numerically accessed, allowing the observation of the displacement and shear strain fields that are compared to direct MRI experiments from the literature. Important conclusions about the correct lamellar-interlamellar-nucleus interactions are provided thanks to the developed model. The critical interlamellar spots with the highest delamination potentials are defined, analyzed and related to the local kinetics and microstructure.

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.

Prediction of biomechanical responses of human lumbar discs - a stochastic finite element model analysis

BACKGROUND

Accurate biomechanical investigation of human intervertebral discs (IVDs) is difficult because of their complicated structural and material features.

AIM

To investigate probabilistic distributions of the biomechanical responses of the IVD by considering varying nonlinear structural and material properties using a stochastic finite element (FE) model.

METHODS

A FE model of a L3-4 disc was reconstructed, including the nucleus pulposus (NP), annular matrix and fibers. A Monte Carlo method was used to randomly generate 500 sets of the nonlinear material properties and fiber orientations of the disc that were implemented into the FE model. The FE model was analyzed under seven loading conditions: a 500 N compressive force, a 7.5Nm moment simulating flexion, extension, left-right lateral bending, and left-right axial rotation, respectively. The distributions of the ranges of motion (ROMs), intradiscal pressures (IDP), fiber stresses and matrix strains of the disc were analyzed.

RESULTS

Under the compressive load, the displacement varied between 0.29 mm and 0.76 mm. Under the 7.5Nm moment, the ROMs varied between 3.0° and 6.0° in primary rotations. The IDPs varied within 0.3 MPa under all the loading conditions. The maximal fiber stress (3.22 ± 0.64 MPa) and matrix strain (0.27 ± 0.12%) were observed under the flexion and extension moments, respectively.

CONCLUSION

The IVD biomechanics could be dramatically affected by the structural and material parameters used to construct the FE model. The stochastic FE model that includes the probabilistic distributions of the structural and material parameters provides a useful approach to analyze the statistical ranges of the biomechanical responses of the IVDs.

A Robust Multiscale and Multiphasic Structure-Based Modeling Framework for the Intervertebral Disc

A comprehensive understanding of multiscale and multiphasic intervertebral disc mechanics is crucial for designing advanced tissue engineered structures aiming to recapitulate native tissue behavior. The bovine caudal disc is a commonly used human disc analog due to its availability, large disc height and area, and similarities in biochemical and mechanical properties to the human disc. Because of challenges in directly measuring subtissue-level mechanics, such as in situ fiber mechanics, finite element models have been widely employed in spinal biomechanics research. However, many previous models use homogenization theory and describe each model element as a homogenized combination of fibers and the extrafibrillar matrix while ignoring the role of water content or osmotic behavior. Thus, these models are limited in their ability in investigating subtissue-level mechanics and stress-bearing mechanisms through fluid pressure. The objective of this study was to develop and validate a structure-based bovine caudal disc model, and to evaluate multiscale and multiphasic intervertebral disc mechanics under different loading conditions and with degeneration. The structure-based model was developed based on native disc structure, where fibers and matrix in the annulus fibrosus were described as distinct materials occupying separate volumes. Model parameters were directly obtained from experimental studies without calibration. Under the multiscale validation framework, the model was validated across the joint-, tissue-, and subtissue-levels. Our model accurately predicted multiscale disc responses for 15 of 16 cases, emphasizing the accuracy of the model, as well as the effectiveness and robustness of the multiscale structure-based modeling-validation framework. The model also demonstrated the rim as a weak link for disc failure, highlighting the importance of keeping the cartilage endplate intact when evaluating disc failure mechanisms in vitro. Importantly, results from this study elucidated important fluid-based load-bearing mechanisms and fiber-matrix interactions that are important for understanding disease progression and regeneration in intervertebral discs. In conclusion, the methods presented in this study can be used in conjunction with experimental work to simultaneously investigate disc joint-, tissue-, and subtissue-level mechanics with degeneration, disease, and injury.

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 Investigation Between Rigid and Semirigid Posterolateral Fixation During Daily Activities: Geometrically Parametric Poroelastic Finite Element Analyses

While spinal fusion using rigid rods remains the gold standard treatment modality for various lumbar degenerative conditions, its adverse effects, including accelerated adjacent segment disease (ASD), are well known. In order to better understand the performance of semirigid constructs using polyetheretherketone (PEEK) in fixation surgeries, the objective of this study was to analyze the biomechanical performance of PEEK versus Ti rods using a geometrically patient-specific poroelastic finite element (FE) analyses. Ten subject-specific preoperative models were developed, and the validity of the models was evaluated with previous studies. Furthermore, FE models of those lumbar spines were regenerated based on postoperation images for posterolateral fixation at the L4–L5 level. Biomechanical responses for instrumented and adjacent intervertebral discs (IVDs) were analyzed and compared subjected to static and cyclic loading. The preoperative model results were well comparable with previous FE studies. The PEEK construct demonstrated a slightly increased range of motion (ROM) at the instrumented level, but decreased ROM at adjacent levels, as compared with the Ti. However, no significant changes were detected during axial rotation. During cyclic loading, disc height loss, fluid loss, axial stress, and collagen fiber strain in the adjacent IVDs were higher for the Ti construct when compared with the intact and PEEK models. Increased ROM, experienced stress in AF, and fiber strain at adjacent levels were observed for the Ti rod group compared with the intact and PEEK rod group, which can indicate the risk of ASD for rigid fixation. Similar to the aforementioned pattern, disc height loss and fluid loss were significantly higher at adjacent levels in the Ti rod group after cycling loading which alter the fluid–solid interaction of the adjacent IVDs. This phenomenon debilitates the damping quality, which results in disc disability in absorbing stress. Such finding may suggest the advantage of using a semirigid fixation system to decrease the chance of ASD.

Biomechanical comparison of spinal cord compression types occurring in Degenerative Cervical Myelopathy

BACKGROUND

Degenerative Cervical Myelopathy results from spine degenerations narrowing the spinal canal and inducing cord compressions. Prognosis is challenging. This study aimed at simulating typical spinal cord compressions observed in patients with a realistic model to better understand pathogenesis for later prediction of patients' evolution.

METHODS

A 30% reduction in cord cross-sectional area at C5-C6 was defined as myelopathy threshold based on Degenerative Cervical Myelopathy features from literature and MRI measurements in 20 patients. Four main compression types were extracted from MRIs and simulated with a comprehensive three-dimensional finite element spine model. Median diffuse, median focal and lateral types were modelled as disk herniation while circumferential type additionally involved ligamentum flavum hypertrophy. All stresses were quantified along inferior-superior axis, compression development and across atlas-defined spinal cord regions.

FINDINGS

Anterior gray and white matter globally received the highest stress while lateral pathways were the least affected. Median diffuse compression induced the highest stresses. Circumferential type focused stresses in posterior gray matter. Along inferior-superior axis, those two types showed a peak of constraints at compression site while median focal and lateral types showed lower values but extending further.

INTERPRETATION

Median diffuse type would be the most detrimental based on stress amplitude. Anterior regions would be the most at risk, except for circumferential type where posterior regions would be equally affected. In addition to applying constraints, ischemia could be a significant component explaining the early demyelination reported in lateral pathways. Moving towards patient-specific simulations, biomechanical models could become strong predictors for degenerative changes.

Level-wise differences in in vivo lateral bending moment are associated with microstructural alterations in bovine caudal intervertebral discs

Despite its common use as a laboratory model, little is known about the in vivo forces and moments applied to the bovine caudal intervertebral disc. Such aspects are crucial, as intervertebral disc tissue is known to remodel in response to repeated loading. We hypothesized that the magnitude of loading from muscle contraction during a typical lateral bending motion varies between caudal levels and is accompanied by variations in tissue microstructure. This hypothesis was tested by estimating level-wise forces and bending moments using two independent approaches: a dynamic analytical model of the motion and analysis of muscle cross-sections obtained via computed tomography. Microstructure was assessed by measuring the collagen fiber crimp period in the annulus fibrosus, and composition was assessed via quantitative histology. Both the analytical model and muscle cross-sections indicated peak bending moments of over 3 N m and peak compressive force of over 125 N at the c1c2 level, decreasing distally. There was a significant downward trend from proximal to distal in the outer annulus fibrosus collagen crimp period in the anterior and posterior regions only, suggesting remodeling in response to the highest lateral bending moments. There were no observed trends in composition. Our results suggest that although the proximal discs in the bovine tail are subjected to forces and moments from muscle contraction that are comparable (relative to disc size) to those acting on human lumbar discs, the distal discs are not. The resulting pattern of microstructural alterations suggests that level-wise differences should be considered when using bovine discs as a research model.

Finite element investigation on the dynamic mechanical properties of low-frequency vibrations on human L2-L3 spinal motion segments with different degrees of degeneration

Exposure to low-frequency vibration is harmful to human lumbar health. However, the dynamic mechanical properties of lumbar spines with varying degrees of degeneration during time-dependent vibration remain incompletely understood. In this study, four poroelastic finite element models of human L2-L3 spinal motion segments, including the non-degeneration and the mild, moderate, and serious degeneration, were established. One-hour low-frequency vibrations with different frequencies were applied. Then, the dynamic mechanical properties of different degenerated lumbar models under the same vibration and the same lumbar model under vibrations at different frequencies were investigated. The results indicated and implied that the negative influences of 1-h vibration on the dynamic mechanical properties of the non-degenerated and mildly degenerated models were similar, but became obvious for the moderately and seriously degenerated models with time. Therefore, the damage caused by low-frequency vibration on the degenerated spinal motion segments was more serious compared with that on the healthy one. Meanwhile, the dynamic mechanical properties of the same lumbar model under vibrations at different frequencies expressed the negligible differences when the vibration frequency was not close to the lumbar natural frequency. Thus, the effects of the 1-h vibrations at different frequencies on one spinal motion segment were similar. Vibration frequency sensitivity analysis on the dynamic characteristics of human L2-L3 spinal motion segments with different degrees of degeneration.