Injury mechanisms of the ligamentous cervical C2-C3 Functional Spinal Unit to complex loading modes: Finite Element study

Quantifying the Importance of Active Muscle Repositioning a Finite Element Neck Model in Flexion Using Kinematic, Kinetic, and Tissue-Level Responses

PURPOSE

Non-neutral neck positions are important initial conditions in impact scenarios, associated with a higher incidence of injury. Repositioning in finite element (FE) neck models is often achieved by applying external boundary conditions (BCs) to the head while constraining the first thoracic vertebra (T1). However, in vivo, neck muscles contract to achieve a desired head and neck position generating initial loads and deformations in the tissues. In the present study, a new muscle-based repositioning method was compared to traditional applied BCs using a contemporary FE neck model for forward head flexion of 30°.

METHODS

For the BC method, an external moment (2.6 Nm) was applied to the head with T1 fixed, while for the muscle-based method, the flexors and extensors were co-contracted under gravity loading to achieve the target flexion.

RESULTS

The kinematic response from muscle contraction was within 10% of the in vivo experimental data, while the BC method differed by 18%. The intervertebral disc forces from muscle contraction were agreeable with the literature (167 N compression, 12 N shear), while the BC methodology underpredicted the disc forces owing to the lack of spine compression. Correspondingly, the strains in the annulus fibrosus increased by an average of 60% across all levels due to muscle contraction compared to BC method.

CONCLUSION

The muscle repositioning method enhanced the kinetic response and subsequently led to differences in tissue-level responses compared to the conventional BC method. The improved kinematics and kinetics quantify the importance of repositioning FE neck models using active muscles to achieve non-neutral neck positions.

Spinal Cord Stress After Anterior Cervical Diskectomy and Fusion: Results from a Patient-Specific Finite Element Model

Degenerative cervical myelopathy (DCM) is the commonest cause of cervical spinal cord dysfunction in older adults and is characterized by spinal cord compression and stress during neck motion. Although surgical decompression eliminates static spinal cord compression, cord stress resulting from flexion-extension motion of the spinal column has not been determined for single and multi-level surgical interventions. The effect of surgery on spinal cord stress is expected to change with the number of surgical levels as well as patient-specific anatomy. Using a MRI-derived patient-specific finite element model, we simulated 1-, 2- and 3-level anterior cervical diskectomy and fusion (ACDF) surgery for DCM. A substantial decrease in spinal cord stress at the level of spinal cord decompression was noted in all simulations. This was associated with a considerable increase in spinal cord stress rostral to the surgical level, and the magnitude of stress was higher in multi-level surgery. Increased spinal cord stress at the rostral adjacent segment correlated with increased segmental range of motion (r = 0.69, p = 0.002) and disk pressure (r = 0.57, p = 0.05). Together, these results indicate that ACDF for DCM is associated with adverse spinal cord stress patterns adjacent to the fusion construct, and further research is needed to determine if the altered stress is associated with clinical outcomes after surgery for DCM.

Stress Distribution on Spinal Cord According to Type of Laminectomy for Large Focal Cervical Ossification of Posterior Longitudinal Ligament Based on Finite Element Method

Most studies on the ossification of the posterior longitudinal ligament (OPLL) using the finite element method were conducted in the neutral state, and the resulting decompression was judged to be good. As these studies do not reflect the actual behavior of the cervical spine, this study conducted an analysis in the neutral state and a biomechanical analysis during flexion and extension behaviors. After validation via the construction of an intact cervical spine model, the focal OPLL model was inserted into the C4–C5 segment and a simulation was performed. The neutral state was shown by applying a fixed condition to the lower part of the T1 and Y-axis fixed condition of the spinal cord and simulating spinal cord compression with OPLL. For flexion and extension simulation, a ±30-degree displacement was additionally applied to the top of the C2 dens. Accordingly, it was confirmed that spinal cord decompression did not work well during the flexion and extension behaviors, but rather increased. Thus, if patients with focal OPLL inevitably need to undergo posterior decompression, additional surgery using an anterior approach should be considered.

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

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

The importance of intervertebral disc material model on the prediction of mechanical function of the cervical spine

BACKGROUND

Linear elastic, hyperelastic, and multiphasic material constitutive models are frequently used for spinal intervertebral disc simulations. While the characteristics of each model are known, their effect on spine mechanical response requires a careful investigation. The use of advanced material models may not be applicable when material constants are not available, model convergence is unlikely, and computational time is a concern. On the other hand, poor estimations of tissue's mechanical response are likely if the spine model is oversimplified. In this study, discrepancies in load response introduced by material models will be investigated.

METHODS

Three fiber-reinforced C2-C3 disc models were developed with linear elastic, hyperelastic, and biphasic behaviors. Three different loading modes were investigated: compression, flexion and extension in quasi-static and dynamic conditions. The deformed disc height, disc fluid pressure, range of motion, and stresses were compared.

RESULTS

Results indicated that the intervertebral disc material model has a strong effect on load-sharing and disc height change when compression and flexion were applied. The predicted mechanical response of three models under extension had less discrepancy than its counterparts under flexion and compression. The fluid-solid interaction showed more relevance in dynamic than quasi-static loading conditions. The fiber-reinforced linear elastic and hyperelastic material models underestimated the load-sharing of the intervertebral disc annular collagen fibers.

CONCLUSION

This study confirmed the central role of the disc fluid pressure in spinal load-sharing and highlighted loading conditions where linear elastic and hyperelastic models predicted energy distribution different than that of the biphasic model.

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

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

Application of Simulation Methods in Cervical Spine Dynamics

Neck injury is one of the most frequent spine injuries due to the complex structure of the cervical spine. The high incidence of neck injuries in collision accidents can bring a heavy economic burden to the society. Therefore, knowing the potential mechanisms of cervical spine injury and dysfunction is significant for improving its prevention and treatment. The research on cervical spine dynamics mainly concerns the fields of automobile safety, aeronautics, and astronautics. Numerical simulation methods are beneficial to better understand the stresses and strains developed in soft tissues with investigators and have been roundly used in cervical biomechanics. In this article, the simulation methods for the development and application of cervical spine dynamic problems in the recent years have been reviewed. The study focused mainly on multibody and finite element models. The structure, material properties, and application fields, especially the whiplash injury, were analyzed in detail. It has been shown that simulation methods have made remarkable progress in the research of cervical dynamic injury mechanisms, and some suggestions on the research of cervical dynamics in the future have been proposed.

Contribution of injured posterior ligamentous complex and intervertebral disc on post-traumatic instability at the cervical spine

Posterior ligamentous complex (PLC) and intervertebral disc (IVD) injuries are common cervical spine flexion-distraction injuries, but the residual stability following their disruption is misknown. The objective of this study was to evaluate the effect of PLC and IVD disruption on post-traumatic cervical spine stability under low flexion moment (2 Nm) using a finite element (FE) model of C2-T1. The PLC was removed first and a progressive disc rupture (one third, two thirds and complete rupture) was modeled to simulate IVD disruption at C2-C3, C4-C5 and C6-C7. At each step, a non-traumatic flexion moment was applied and the change in stability was evaluated. PLC removal had little impact at C2-C3 but increased local range of motion (ROM) at the injured level by 77.2% and 190.7% at C4-C5 and C6-C7, respectively. Complete IVD rupture had the largest impact on C2-C3, increasing C2-C3 ROM by 181% and creating a large antero-posterior displacement of the C2-C3 segment. The FE analysis showed PLC and disc injuries create spinal instability. However, the PLC played a bigger role in the stability of the middle and lower cervical spine while the IVD was more important at the upper cervical spine. Stabilization appears important when managing patients with soft tissue injuries.

Comparative analysis of the biomechanics of the adjacent segments after minimally invasive cervical surgeries versus anterior cervical discectomy and fusion: A finite element study

Purpose

Percutaneous full-endoscopic anterior cervical discectomy (PEACD) and posterior cervical foraminotomy (PCF) as alternatives to anterior cervical discectomy and fusion (ACDF) are extensively used in the treatment of patients with cervical spondylotic radiculopathy. The possibility of avoiding the risk of accelerated degeneration of the adjacent segments caused by fusion is claimed to be the theoretical advantage of these approaches; however, there is a paucity of supportive evidence from biomechanical data. Therefore, this study investigated and compared the effects of PCF, PEACD, and ACDF on the adjacent segments and operative segments of the cervical spine from a biomechanical standpoint.

Method

A normal and intact three-dimensional finite element digital model of C4-C7 was constructed and validated, and the finite element models of PEACD, PCF, and ACDF were obtained by modifying the C4-C7 model. All models were exposed to identical conditions of load during flexion, extension, axial rotation, and lateral bending. We calculated the range of motion (ROM), intervertebral disc pressure (IDP), and facet joint contact force (FJCF) of the operative segment and the adjacent segment in different motion conditions.

Result

The conventional ACDF had a remarkable influence on the ROM and IDP of the operative segment and the adjacent segments. In the PEACD model, the change of ROM was not noticeable; the IDP of the operative segment was significantly smaller, whereas the change of IDP of the adjacent segment was insignificant. In the PCF model, the ROM and IDP of all segments remained unaffected.During extension, the facet joint contact force changed significantly after ACDF, and it changed slightly after PECAD and PCF.

Conclusion

By comparatively analyzing the biomechanical changes of the cervical spine after PCF, PEACD, and ACDF using the finite element method, we suggested that PCF and PEACD were more suitable for surgical intervention of cervical spondylotic radiculopathy than ACDF from a biomechanical point of view and PCF may outperform PEACD.

Isolated Cyclic Loading During Adolescence Improves Tibial Bone Microstructure and Strength at Adulthood

Bone is a unique living tissue, which responds to the mechanical stimuli regularly imposed on it. Adolescence facilitates a favorable condition for the skeleton that enables the exercise to positively influence bone architecture and overall strength. However, it is still dubious for how long the skeletal benefits gained in adolescence is preserved at adulthood. The current study aims to use a rat model to investigate the effects of in vivo low- (LI), medium- (MI), and high- (HI) intensity cyclic loadings applied during puberty on longitudinal bone development, morphometry, and biomechanics during adolescence as well as at adulthood. Forty-two young (4-week-old) male rats were randomized into control, sham, LI, MI, and HI groups. After a 5 day/week for 8 weeks cyclic loading regime applied on the right tibia, loaded rats underwent a subsequent 41-week, normal cage activity period. Right tibias were removed at 52 weeks of age, and a comprehensive assessment was performed using μCT, mechanical testing, and finite element analysis. HI and MI groups exhibited reduced body weight and food intake at the end of the loading period compared with shams, but these effects disappeared afterward. HI cyclic loading increased BMD, bone volume fraction, trabecular thickness, trabecular number, and decreased trabecular spacing after loading. All loading-induced benefits, except BMD, persisted until the end of the normal cage activity period. Moreover, HI loading induced enhanced bone area, periosteal perimeter, and moment of inertia, which remained up to the 52nd week. After the normal cage activity at adulthood, the HI group showed increased ultimate force and stress, stiffness, postyield displacement and energy, and toughness compared with the sham group. Overall, our findings suggest that even though both trabecular and cortical bone drifted through age-related changes during aging, HI cyclic loading performed during adolescence can render lifelong benefits in bone microstructure and biomechanics. © 2020 The Authors. JBMR Plus published by Wiley Periodicals, Inc. on behalf of American Society for Bone and Mineral Research.

High Impact Exercise Improves Bone Microstructure and Strength in Growing Rats

Physical activity is beneficial for skeletal development. However, impact sports during adolescence, leading to bone growth retardation and/or bone quality improvement, remains unexplained. This study investigated the effects of in vivo low (LI), medium (MI), and high (HI) impact loadings applied during puberty on bone growth, morphometry and biomechanics using a rat model. 4-week old rats (n = 30) were divided into control, sham, LI, MI, and HI groups. The impact was applied on the right tibiae, 5 days/week for 8 weeks mimicking walking (450 µε), uphill running (850 µε) and jumping (1250 µε) conditions. Trabecular and cortical parameters were determined by micro-CT, bone growth rate by calcein labeling and toluidine blue staining followed by histomorphometry. Bio-mechanical properties were evaluated from bending tests. HI group reduced rat body weight and food consumption compared to shams. Bone growth rate also decreased in MI and HI groups despite developing thicker hypertrophic and proliferative zone heights. HI group showed significant increment in bone mineral density, trabecular thickness, cortical and total surface area. Ultimate load and stiffness were also increased in MI and HI groups. We conclude that impact loading during adolescence reduces bone growth moderately but improves bone quality and biomechanics at the end of the growing period.

Can repeated in vivo micro-CT irradiation during adolescence alter bone microstructure, histomorphometry and longitudinal growth in a rodent model?

In vivo micro-computed tomography (micro-CT) can monitor longitudinal changes in bone mass and microstructure in small rodents but imposing high doses of radiation can damage the bone tissue. However, the effect of weekly micro-CT scanning during the adolescence on bone growth and architecture is still unknown. The right proximal tibia of male Sprague-Dawley rats randomized into three dose groups of 0.83, 1.65 and 2.47 Gy (n = 11/group) were CT scanned at weekly intervals from 4th to 12th week of age. The left tibia was used as a control and scanned only at the last time point. Bone marrow cells were investigated, bone growth rates and histomorphometric analyses were performed, and bone structural parameters were determined for both left and right tibiae. Radiation doses of 1.65 and 2.47 Gy affected bone marrow cells, heights of the proliferative and hypertrophic zones, and bone growth rates in the irradiated tibiae. For the 1.65 Gy group, irradiated tibiae resulted in lower BMD, Tb.Th, Tb.N and a higher Tb.Sp compared with the control tibiae. A decrease in BMD, BV/TV, Tb.Th, Tb.N and an increase in Tb.Sp were observed between the irradiated and control tibiae for the 2.47 Gy group. For cortical bone parameters, no effects were noticed for 1.65 and 0.83 Gy groups, but a lower Ct.Th was observed for 2.47 Gy group. Tibial bone development was adversely impacted and trabecular bone, together with bone marrow cells, were negatively affected by the 1.65 and 2.47 Gy radiation doses. Cortical bone microstructure was affected for 2.47 Gy group. However, bone development and morphometry were not affected for 0.83 Gy group. These findings can be used as a proof of concept for using the reasonable high-quality image acquisition under 0.83 Gy radiation doses during the adolescent period of rats without interfering with the bone development process.

Angle-ply biomaterial scaffold for annulus fibrosus repair replicates native tissue mechanical properties, restores spinal kinematics, and supports cell viability

Annulus fibrosus (AF) damage commonly occurs due to intervertebral disc (IVD) degeneration/herniation. The dynamic mechanical role of the AF is essential for proper IVD function and thus it is imperative that biomaterials developed to repair the AF withstand the mechanical rigors of the native tissue. Furthermore, these biomaterials must resist accelerated degradation within the proteolytic environment of degenerate IVDs while supporting integration with host tissue. We have previously reported a novel approach for developing collagen-based, multi-laminate AF repair patches (AFRPs) that mimic the angle-ply architecture and basic tensile properties of the human AF. Herein, we further evaluate AFRPs for their: tensile fatigue and impact burst strength, IVD attachment strength, and contribution to functional spinal unit (FSU) kinematics following IVD repair. Additionally, AFRP resistance to collagenase degradation and cytocompatibility were assessed following chemical crosslinking. In summary, AFRPs demonstrated enhanced durability at high applied stress amplitudes compared to human AF and withstood radially-directed biaxial stresses commonly borne by the native tissue prior to failure/detachment from IVDs. Moreover, FSUs repaired with AFRPs and nucleus pulposus (NP) surrogates had their axial kinematic parameters restored to intact levels. Finally, carbodiimide crosslinked AFRPs resisted accelerated collagenase digestion without detrimentally effecting AFRP tensile properties or cytocompatibility. Taken together, AFRPs demonstrate the mechanical robustness and enzymatic stability required for implantation into the damaged/degenerate IVD while supporting AF cell infiltration and viability.

STATEMENT OF SIGNIFICANCE

The quality of life for millions of individuals globally is detrimentally impacted by IVD degeneration and herniation. These pathologies often result in the structural demise of IVD tissue, particularly the annulus fibrosus (AF). Biomaterials developed for AF repair have yet to demonstrate the mechanical strength and durability required for utilization in the spine. Herein, we demonstrate the development of an angle-ply AF repair patch (AFRP) that can resist the application of physiologically relevant stresses without failure and which contributes to the restoration of functional spinal unit axial kinematics following repair. Furthermore, we show that this biomaterial can resist accelerated degradation in a simulated degenerate environment and supports AF cell viability.