The importance of fluid-structure interaction in spinal trauma models

Cervical Column and Cord and Column Responses in Whiplash With Stenosis: A Finite Element Modeling Study

Spine degeneration is a normal aging process. It may lead to stenotic spines that may have implications for pain and quality of life. The diagnosis is based on clinical symptomatology and imaging. Magnetic resonance images often reveal the nature and degree of stenosis of the spine. Stenosis is concerning to clinicians and patients because of the decreased space in the spinal canal and potential for elevated risk of cord and/or osteoligamentous spinal column injuries. Numerous finite element models of the cervical spine have been developed to study the biomechanics of the osteoligamentous column such as range of motion and vertebral stress; however, spinal cord modeling is often ignored. The objective of this study was to determine the external column and internal cord and disc responses of stenotic spines using finite element modeling. A validated model of the subaxial spinal column was used. The osteoligamentous column was modified to include the spinal cord. Mild, moderate, and severe degrees of stenosis commonly identified in civilian populations were simulated at C5-C6. The column-cord model was subjected to postero-anterior acceleration at T1. The range of motion, disc pressure, and cord stress-strain were obtained at the index and superior and inferior adjacent levels of the stenosis. The external metric representing the segmental motion was insensitive while the intrinsic disc and cord variables were more sensitive, and the index level was more affected by stenosis. These findings may influence surgical planning and patient education in personalized medicine.

Finite element modeling and analysis of effect of preexisting cervical degenerative disease on the spinal cord during flexion and extension

Recent studies have emphasized the importance of dynamic activity in the development of myelopathy. However, current knowledge of how degenerative factors affect the spinal cord during motion is still limited. This study aimed to investigate the effect of various types of preexisting herniated cervical disc and the ligamentum flavum ossification on the spinal cord during cervical flexion and extension. A detailed dynamic fluid-structure interaction finite element model of the cervical spine with the spinal cord was developed and validated. The changes of von Mises stress and maximum principal strain within the spinal cord in the period of normal, hyperflexion, and hyperextension were investigated, considering various types and grades of disc herniation and ossification of the ligamentum flavum. The flexion and extension of the cervical spine with spinal canal encroachment induced high stress and strain inside the spinal cord, and this effect was also amplified by increased canal encroachments and cervical hypermobility. The spinal cord might evade lateral encroachment, leading to a reduction in the maximum stress and principal strain within the spinal cord in local-type herniation. Although the impact was limited in the case of diffuse type, the maximum stress tended to appear in the white matter near the encroachment site while compression from both ventral and dorsal was essential to make maximum stress appear in the grey matter. The existence of canal encroachment can reduce the safe range for spinal cord activities, and hypermobility activities may induce spinal cord injury. Besides, the ligamentum flavum plays an important role in the development of central canal syndrome.Significance. This model will enable researchers to have a better understanding of the influence of cervical degenerative diseases on the spinal cord during extension and flexion.

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

BACKGROUND

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Effect of Different Types of Ossification of the Posterior Longitudinal Ligament on the Dynamic Biomechanical Response of the Spinal Cord: A Finite Element Analysis

Ossification of the posterior longitudinal ligament (OPLL) has been identified as an important cause of cervical myelopathy. However, the biomechanical mechanism between the OPLL type and the clinical characteristics of myelopathy remains unclear. The aim of this study was to evaluate the effect of different types of OPLL on the dynamic biomechanical response of the spinal cord. A three-dimensional finite element model of the fluid-structure interaction of the cervical spine with spinal cord was established and validated. The spinal cord stress and strain, cervical range of motion (ROM) in different types of OPLL models were predicted during dynamic flexion and extension activity. Different types of OPLL models showed varying degrees of increase in stress and strain under the process of flexion and extension, and there was a surge toward the end of extension. Larger spinal cord stress was observed in segmental OPLL. For continuous and mixed types of OPLL, the adjacent segments of OPLL showed a dramatic increase in ROM, while the ROM of affected segments was limited. As a dynamic factor, flexion and extension of the cervical spine play an amplifying role in OPLL-related myelopathy, while appropriate spine motion is safe and permitted. Segmental OPLL patients are more concerned about the spinal cord injury induced by large stress, and patients with continuous OPLL should be noted to progressive injuries of adjacent level.

The biomechanical effect of different types of ossification of the ligamentum flavum on the spinal cord during cervical dynamic activities

Ossification of the ligamentum flavum (OLF) is thought to be an influential etiology of myelopathy, as thickened ligamentum flavum causes the stenosis of the vertebral canal, which could subsequently compress the spinal cord. Unfortunately, there was little information available on the effects of cervical OLF on spinal cord compression, such as the relationship between the progression of cervical OLF and nervous system symptoms during dynamic cervical spine activities. In this research, a finite element model of C1-C7 including the spinal cord featured by dynamic fluid-structure interaction was reconstructed and utilized to analyze how different types of cervical OLF affect principal strain and stress distribution in spinal cord during spinal activities towards six directions. For patients with cervical OLF, cervical extension induces higher stress within the spinal cord among all directions. From the perspective of biomechanics, extension leads to stress concentration in the lateral corticospinal tracts or the posterior of gray matter. Low energy damage to the spinal cord would be caused by the high and fluctuating stresses during cervical movements to the affected side for patients with unilateral OLF at lower grades.

Finite element modeling of the human cervical spinal cord and its applications: A systematic review

Background Context

Finite element modeling (FEM) is an established tool to analyze the biomechanics of complex systems. Advances in computational techniques have led to the increasing use of spinal cord FEMs to study cervical spinal cord pathology. There is considerable variability in the creation of cervical spinal cord FEMs and to date there has been no systematic review of the technique. The aim of this study was to review the uses, techniques, limitations, and applications of FEMs of the human cervical spinal cord.

Methods

A literature search was performed through PubMed and Scopus using the words finite element analysis, spinal cord, and biomechanics. Studies were selected based on the following inclusion criteria: (1) use of human spinal cord modeling at the cervical level; (2) model the cervical spinal cord with or without the osteoligamentous spine; and (3) the study should describe an application of the spinal cord FEM.

Results

Our search resulted in 369 total publications, 49 underwent reviews of the abstract and full text, and 23 were included in the study. Spinal cord FEMs are used to study spinal cord injury and trauma, pathologic processes, and spine surgery. Considerable variation exists in the derivation of spinal cord geometries, mathematical models, and material properties. Less than 50% of the FEMs incorporate the dura mater, cerebrospinal fluid, nerve roots, and denticulate ligaments. Von Mises stress, and strain of the spinal cord are the most common outputs studied. FEM offers the opportunity for dynamic simulation, but this has been used in only four studies.

Conclusions

Spinal cord FEM provides unique insight into the stress and strain of the cervical spinal cord in various pathological conditions and allows for the simulation of surgical procedures. Standardization of modeling parameters, anatomical structures and inclusion of patient-specific data are necessary to improve the clinical translation.

Tensile properties of human spinal dura mater and pericranium

Autologous pericranium is a promising dural graft material. An optimal graft should exhibit similar mechanical properties to the native dura, but the mechanical properties of human pericranium have not been characterized, and studies of the biomechanical performance of human spinal dura are limited. The primary aim of this study was to measure the tensile structural and material properties of the pericranium, in the longitudinal and circumferential directions, and of the dura in each spinal region (cervical, thoracic and lumbar) and in three directions (longitudinal anterior and posterior, and circumferential). The secondary aim was to determine corresponding constitutive stress-strain equations using a one-term Ogden model. A total of 146 specimens were tested from 7 cadavers. Linear regression models assessed the effect of tissue type, region, and orientation on the structural and material properties. Pericranium was isotropic, while spinal dura was anisotropic with higher stiffness and strength in the longitudinal than the circumferential direction. Pericranium had lower strength and modulus than spinal dura across all regions in the longitudinal direction but was stronger and stiffer than dura in the circumferential direction. Spinal dura and pericranium had similar strain at peak force, toe, and yield, across all regions and directions. Human pericranium exhibits isotropic mechanical behavior that lies between that of the longitudinal and circumferential spinal dura. Further studies are required to determine if pericranium grafts behave like native dura under in vivo loading conditions. The Ogden parameters reported may be used for computational modeling of the central nervous system. Graphical abstract.

Effect of degenerative factors on cervical spinal cord during flexion and extension: a dynamic finite element analysis

Spinal cord injury (SCI) is a global problem that brings a heavy burden to both patients and society. Recent investigations indicated degenerative disease is taking an increasing part in SCI with the growth of the aging population. However, little insight has been gained about the effect of cervical degenerative disease on the spinal cord during dynamic activities. In this work, a dynamic fluid-structure interaction model was developed and validated to investigate the effect of anterior and posterior encroachment caused by degenerative disease on the spinal cord during normal extension and flexion. Maximum von-Mises stress and maximum principal strain were observed at the end of extension and flexion. The abnormal stress distribution caused by degenerative factors was concentrated in the descending tracts of the spinal cord. Our finding indicates that the excessive motion of the cervical spine could potentially exacerbate spinal cord injury and enlarge injury areas. Stress and strain remained low compared to extension during moderate flexion. This suggests that patients with cervical degenerative disease should avoid frequent or excessive flexion and extension which could result in motor function impairment, whereas moderate flexion is safe. Besides, encroachment caused by degenerative factors that are not significant in static imaging could also cause cord compression during normal activities.

Effect of Velocity and Contact Stress Area on the Dynamic Behavior of the Spinal Cord Under Different Testing Conditions

Knowledge of the dynamic behavior of the spinal cord under different testing conditions is critical for our understanding of biomechanical mechanisms of spinal cord injury. Although velocity and contact stress area are known to affect external mechanical stress or energy upon sudden traumatic injury, quantitative investigation of the two clinically relevant biomechanical variables is limited. Here, freshly excised rat spinal-cord–pia-arachnoid constructs were tested through indentation using indenters of different sizes (radii: 0.25, 0.50, and 1.00 mm) at various loading rates ranging from 0.04 to 0.20 mm/s. This analysis found that the ex vivo specimen displayed significant nonlinear viscoelasticity at <10% of specimen thickness depth magnitudes. At higher velocity and larger contact stress area, the cord withstood a higher peak load and exhibited more sensitive mechanical relaxation responses (i.e., increasing amplitude and speed of the drop in peak load). Additionally, the cord became stiffer (i.e., increasing elastic modulus) and softer (i.e., decreasing elastic modulus) at a higher velocity and larger contact stress area, respectively. These findings will improve our understanding of the real-time complex biomechanics involved in traumatic spinal cord injury.

Comparison of numerical methods for cerebrospinal fluid representation and fluid-structure interaction during transverse impact of a finite element spinal cord model

Spinal cord impacts can have devastating consequences. Computational models can investigate such impacts but require biofidelic numerical representations of the neural tissues and fluid-structure interaction with cerebrospinal fluid. Achieving this biofidelity is challenging, particularly for efficient implementation of the cerebrospinal fluid in full computational human body models. The goal of this study was to assess the biofidelity and computational efficiency of fluid-structure interaction methods representing the cerebrospinal fluid interacting with the spinal cord, dura, and pia mater using experimental pellet impact test data from bovine spinal cords. Building on an existing finite element model of the spinal cord and pia mater, an orthotropic hyperelastic constitutive model was proposed for the dura mater and fit to literature data. The dura mater and cerebrospinal fluid were integrated with the existing finite element model to assess four fluid-structure interaction methods under transverse impact: Lagrange, pressurized volume, smoothed particle hydrodynamics, and arbitrary Lagrangian-Eulerian. The Lagrange method resulted in an overly stiff mechanical response, whereas the pressurized volume method over-predicted compression of the neural tissues. Both the smoothed particle hydrodynamics and arbitrary Lagrangian-Eulerian methods were able to effectively model the impact response of the pellet on the dura mater, outflow of the cerebrospinal fluid, and compression of the spinal cord; however, the arbitrary Lagrangian-Eulerian compute time was approximately five times higher than smoothed particle hydrodynamics. Crucial to implementation in human body models, the smoothed particle hydrodynamics method provided a computationally efficient and representative approach to model spinal cord fluid-structure interaction during transverse impact.

A Hyper-Viscoelastic Continuum-Level Finite Element Model of the Spinal Cord Assessed for Transverse Indentation and Impact Loading

Finite Element (FE) modelling of spinal cord response to impact can provide unique insights into the neural tissue response and injury risk potential. Yet, contemporary human body models (HBMs) used to examine injury risk and prevention across a wide range of impact scenarios often lack detailed integration of the spinal cord and surrounding tissues. The integration of a spinal cord in contemporary HBMs has been limited by the need for a continuum-level model owing to the relatively large element size required to be compatible with HBM, and the requirement for model development based on published material properties and validation using relevant non-linear material data. The goals of this study were to develop and assess non-linear material model parameters for the spinal cord parenchyma and pia mater, and incorporate these models into a continuum-level model of the spinal cord with a mesh size conducive to integration in HBM. First, hyper-viscoelastic material properties based on tissue-level mechanical test data for the spinal cord and hyperelastic material properties for the pia mater were determined. Secondly, the constitutive models were integrated in a spinal cord segment FE model validated against independent experimental data representing transverse compression of the spinal cord-pia mater complex (SCP) under quasi-static indentation and dynamic impact loading. The constitutive model parameters were fit to a quasi-linear viscoelastic model with an Ogden hyperelastic function, and then verified using single element test cases corresponding to the experimental strain rates for the spinal cord (0.32–77.22 s⁻¹) and pia mater (0.05 s⁻¹). Validation of the spinal cord model was then performed by re-creating, in an explicit FE code, two independent ex-vivo experimental setups: 1) transverse indentation of a porcine spinal cord-pia mater complex and 2) dynamic transverse impact of a bovine SCP. The indentation model accurately matched the experimental results up to 60% compression of the SCP, while the impact model predicted the loading phase and the maximum deformation (within 7%) of the SCP experimental data. This study quantified the important biomechanical contribution of the pia mater tissue during spinal cord deformation. The validated material models established in this study can be implemented in computational HBM.

Numerical Models of Spinal Cord Trauma: The Effect of Cerebrospinal Fluid Pressure and Epidural Fat on the Results

Spinal cord injury (SCI) is commonly caused by traumatic mechanical damage. Although numerical models can help predict the mechanics of SCI without putting the subjects in danger, previous studies did not focus on alternations in cerebrospinal fluid (CSF) pressure and did not account for the presence of epidural fat. This study aims to numerically compare the mechanical behavior of the human spine when subjected to contusion and burst fracture with varying CSF pressure, either normal or elevated pressure that represents intracranial hypertension. An additional aim is to find out how the presence of the fat in the model affects the SCI calculations. CSF and epidural fat were modeled as smoothed-particle hydrodynamics (SPH) and the soft tissues were modeled as hyperelastic. This approach made it possible to account for CSF pressure alteration and its effect on the cord. Validation models resulted in good correlation with previous numerical and experimental studies. The results were able to capture the fluid dynamics of the CSF while demonstrating a considerable change in the stresses of the spinal cord. The comparison of the CSF pressures demonstrated that SCI in patients with elevated pressure and in regions where insufficient epidural fat exists might lead to higher spinal cord stresses. Yet, in regions with enough fat, the fat can absorb energy and counteract the effect of the elevated pressure. These results indicate important aspects that need to be accounted for in future numerical models of SCI while also demonstrating how the injury might be aggravated by preexisting conditions.

Effects of a contusion load on spinal cord with different curvatures

The cervical spine injury is a complicated procedure in the combination of different injury loads and postures. The aim of this study is to investigate the injury mechanism considering different types of cervical curvatures subjected to contusion loads. A finite element model of a cervical spinal cord was constructed. Gray matter, white matter and pia matter were modeled and hyperelastic material properties were assigned. Convergence analysis and validation analysis were carried out. The model was simulated in 3 different spinal curvatures and loaded by 2 directions with 4 compression degrees. The maximum von Mises stress in the whole model was concentrated in the pia matter in all loading cases. When investigating spinal cord injury, the pia matter must be considered. For all three curvatures, the stress in the gray matter and white matter was higher in front-to-back loading condition than that in back-to-front loading condition. The front-to-back impact may cause a larger damage. A back-to-front load damaged the structure around the central canal and a front-to-back contusion load damaged the anterior horn of the spinal cord at most time. From the view of the maximum stress, the lordotic curvature did not show significant buffering effect. However, the pathological curvature had large areas affected and the lordotic curvature showed some benefits to some degree from the view of stress distribution.

Accurate simulation of the herniated cervical intervertebral disc using controllable expansion: a finite element study

Expansions were carried out in finite element (FE) models of disc hernia including symmetric (median, lateral, paramedian) and asymmetric types. In all models, lubricous disk bulging that applied a linear compression to the anterior part of the cord was observed at the posterior surfaces of the expansion zone, respectively. The shape and position of protrusions varyed with the temperature, magnitude, and location of expanding elements. The geometric deformation and stress distribution of the spinal cord increased as the extent of compression grew. This method is in possession of enormous potential in promoting further individualized research of cervical spondylotic myelopathy.

Correlating Tissue Mechanics and Spinal Cord Injury: Patient-Specific Finite Element Models of Unilateral Cervical Contusion Spinal Cord Injury in Non-Human Primates

Non-human primate (NHP) models are the closest approximation of human spinal cord injury (SCI) available for pre-clinical trials. The NHP models, however, include broader morphological variability that can confound experimental outcomes. We developed subject-specific finite element (FE) models to quantify the relationship between impact mechanics and SCI, including the correlations between FE outcomes and tissue damage. Subject-specific models of cervical unilateral contusion SCI were generated from pre-injury MRIs of six NHPs. Stress and strain outcomes were compared with lesion histology using logit analysis. A parallel generic model was constructed to compare the outcomes of subject-specific and generic models. The FE outcomes were correlated more strongly with gray matter damage (0.29 < R² < 0.76) than white matter (0.18 < R² < 0.58). Maximum/minimum principal strain, Von-Mises and Tresca stresses showed the strongest correlations (0.31 < R² < 0.76) with tissue damage in the gray matter while minimum principal strain, Von-Mises stress, and Tresca stress best predicted white matter damage (0.23 < R² < 0.58). Tissue damage thresholds varied for each subject. The generic FE model captured the impact biomechanics in two of the four models; however, the correlations between FE outcomes and tissue damage were weaker than the subject-specific models (gray matter [0.25 < R² < 0.69] and white matter [R² < 0.06] except for one subject [0.26 < R² < 0.48]). The FE mechanical outputs correlated with tissue damage in spinal cord white and gray matters, and the subject-specific models accurately mimicked the biomechanics of NHP cervical contusion impacts.

Biomechanical Effects on Cervical Spinal Cord and Nerve Root Following Laminoplasty for Ossification of the Posterior Longitudinal Ligament in the Cervical Spine: A Comparison Between Open-Door and Double-Door Laminoplasty Using Finite Element Analysis

Many clinical case series have reported the predisposing factors for C5 palsy and have presented comparisons of the two types of laminoplasty. However, there have been no biomechanical studies focusing on cervical spinal cord and nerve root following laminoplasty. The purpose of this study is to investigate biomechanical changes in the spinal cord and nerve roots following the two most common types of laminoplasty, open-door and double-door laminoplasty, for cervical ossification of the posterior longitudinal ligament (OPLL). A finite element (FE) model of the cervical spine and spinal cord with nerve root complex structures was developed. Stress changes in the spinal cord and nerve roots, posterior shift of the spinal cord, and displacement of the cervical nerve roots were analyzed with two types of cervical laminoplasty models for variations in the degree of canal occupying ratio and shape of the OPLL. The shape and degree of spinal cord compression caused by the OPLL had more influence on the changes in stress, posterior shift of the spinal cord, and displacement of the nerve root than the type of laminoplasty. The lateral-type OPLL resulted in imbalanced stress on the nerve roots and the highest nerve root displacement. Type of laminoplasty and shape and degree of spinal cord compression caused by OPLL were found to influence the changes in stress and posterior displacement of the cervical spinal cord and nerve roots. Lateral-type OPLL might contribute to the development of C5 palsy due to the imbalanced stress and tension on the nerve roots after laminoplasty.

Engineering approaches to understanding mechanisms of spinal column injury leading to spinal cord injury

BACKGROUND

The mechanical interactions occurring between the spinal column and spinal cord during an injury event are complex and variable, and likely have implications for the clinical presentation and prognosis of the individual.

METHODS

The engineering approaches that have been developed to better understand spinal column and cord interactions during an injury event are discussed. These include injury models utilising human and animal cadaveric specimens, in vivo anaesthetised animals, finite element models, inanimate physical systems and combinations thereof.

FINDINGS

The paper describes the development of these modelling approaches, discusses the advantages and disadvantages of the various models, and the major outcomes that have had implications for spinal cord injury research and clinical practice.

INTERPRETATION

The contribution of these four engineering approaches to understanding the interaction between the biomechanics and biology of spinal cord injury is substantial; they have improved our understanding of the factors contributing to the spinal column disruption, the degree of spinal cord deformation or motion, and the resultant neurological deficit and imaging features. Models of the injury event are challenging to produce, but technological advances are likely to improve these models and, consequently, our understanding of the mechanical context in which the biological injury occurs.

Viscoelasticity of spinal cord and meningeal tissues

Compared to the outer dura mater, the mechanical behavior of spinal pia and arachnoid meningeal layers has received very little attention in the literature. This is despite experimental evidence of their importance with respect to the overall spinal cord stiffness and recovery following compression. Accordingly, inclusion of the mechanical contribution of the pia and arachnoid maters would improve the predictive accuracy of finite element models of the spine, especially in the distribution of stresses and strain through the cord's cross-section. However, to-date, only linearly elastic moduli for what has been previously identified as spinal pia mater is available in the literature. This study is the first to quantitatively compare the viscoelastic behavior of isolated spinal pia-arachnoid-complex, neural tissue of the spinal cord parenchyma, and intact construct of the two. The results show that while it only makes up 5.5% of the overall cross-sectional area, the thin membranes of the innermost meninges significantly affect both the elastic and viscous response of the intact construct. Without the contribution of the pia and arachnoid maters, the spinal cord has very little inherent stiffness and experiences significant relaxation when strained. The ability of the fitted non-linear viscoelastic material models of each condition to predict independent data within experimental variability supports their implementation into future finite element computational studies of the spine.

STATEMENT OF SIGNIFICANCE

The neural tissue of the spinal cord is surrounded by three fibrous layers called meninges which are important in the behavior of the overall spinal-cord-meningeal construct. While the mechanical properties of the outermost layer have been reported, the pia mater and arachnoid mater have received considerably less attention. This study is the first to directly compare the behavior of the isolated neural tissue of the cord, the isolated pia-arachnoid complex, and the construct of these individual components. The results show that, despite being very thin, the inner meninges significantly affect the elastic and time-dependent response of the spinal cord, which may have important implications for studies of spinal cord injury.

Comparison of in vivo and ex vivo viscoelastic behavior of the spinal cord

Despite efforts to simulate the in vivo environment, post-mortem degradation and lack of blood perfusion complicate the use of ex vivo derived material models in computational studies of spinal cord injury. In order to quantify the mechanical changes that manifest ex vivo, the viscoelastic behavior of in vivo and ex vivo porcine spinal cord samples were compared. Stress-relaxation data from each condition were fit to a non-linear viscoelastic model using a novel characterization technique called the direct fit method. To validate the presented material models, the parameters obtained for each condition were used to predict the respective dynamic cyclic response. Both ex vivo and in vivo samples displayed non-linear viscoelastic behavior with a significant increase in relaxation with applied strain. However, at all three strain magnitudes compared, ex vivo samples experienced a higher stress and greater relaxation than in vivo samples. Significant differences between model parameters also showed distinct relaxation behaviors, especially in non-linear relaxation modulus components associated with the short-term response (0.1-1 s). The results of this study underscore the necessity of utilizing material models developed from in vivo experimental data for studies of spinal cord injury, where the time-dependent properties are critical. The ability of each material model to accurately predict the dynamic cyclic response validates the presented methodology and supports the use of the in vivo model in future high-resolution finite element modeling efforts.

STATEMENT OF SIGNIFICANCE

Neural tissues (such as the brain and spinal cord) display time-dependent, or viscoelastic, mechanical behavior making it difficult to model how they respond to various loading conditions, including injury. Methods that aim to characterize the behavior of the spinal cord almost exclusively use ex vivo cadaveric or animal samples, despite evidence that time after death affects the behavior compared to that in a living animal (in vivo response). Therefore, this study directly compared the mechanical response of ex vivo and in vivo samples to quantify these differences for the first time. This will allow researchers to draw more accurate conclusions about spinal cord injuries based on ex vivo data (which are easier to obtain) and emphasizes the importance of future in vivo experimental animal work.

Multimodal sensor-based weight drop spinal cord impact system for large animals

BACKGROUND CONTEXT

A conventional weight drop spinal cord (SC) impact system for large animals is composed of a high-speed video camera, a vision system, and other things. However, a camera with high speed at over 5,000 frames per second (FPS) is very expensive. In addition, the use of the vision system involves complex pattern recognition algorithms and accurate arrangement of the camera and the target.

PURPOSE

The purpose of this study was to develop a large animal spinal cord injury (SCI) modeling system using a multimodal sensor instead of a high-speed video camera and vision system. Another objective of this study was to demonstrate the possibility of the developed system to measure the impact parameters in the experiments using different stiffness materials and an in vivo porcine SC.

STUDY DESIGN

A multimodal sensor-based SCI impact system was developed for large animals. The experiments to measure SC impact parameters were then performed using three different stiffness materials and a Yucatan miniature pig to verify the performance of system developed.

METHODS

A comparative experiment was performed using three different stiffness materials such as high-density (HD) sponge, rubber, and clay to demonstrate the system and perform measurement for impact parameters such as impact velocity, impulsive force, and maximally compressed displacement reflecting physical properties of materials. In the animal experiment, a female Yucatan miniature pig of 60-kg weight was used. Impact conditions for all experiments were fixed at freefalling object mass of 50 g and height of 20 cm.

RESULTS

In the impact test, measured impact velocities were almost the same for the three different stiffness materials at 1.84±0.0153 m/s. Impulsive forces for the three materials of rubber, HD sponge, and clay were 50.88 N, 32.35 N, and 6.68 N, respectively. Maximally compressed displacements for rubber, HD sponge, and clay were 1.93 mm, 3.35 mm, and 15.01 mm, respectively. In the pig experiment, impact velocity, impulsive force, and maximally compressed dural displacement were measured at 1.84 m/s, 13.35 N, and 3.04 mm, respectively. After 3 days from the experiment, paralysis was confirmed for the lower half body of the experimental pig.

CONCLUSIONS

Through experiments, it was verified that our proposed system could be used to measure the SC impact parameters and induce SCI for large animals.

Spinal Meninges and Their Role in Spinal Cord Injury: A Neuroanatomical Review

Current recommendations support early surgical decompression and blood pressure augmentation after traumatic spinal cord injury (SCI). Elevated intraspinal pressure (ISP), however, has probably been underestimated in the pathophysiology of SCI. Recent studies provide some evidence that ISP measurements and durotomy may be beneficial for individuals suffering from SCI. Compression of the spinal cord against the meninges in SCI patients causes a "compartment-like" syndrome. In such cases, intentional durotomy with augmentative duroplasty to reduce ISP and improve spinal cord perfusion pressure (SCPP) may be indicated. Prior to performing these procedures routinely, profound knowledge of the spinal meninges is essential. Here, we provide an in-depth review of relevant literature along with neuroanatomical illustrations and imaging correlates.

Increased stress and strain on the spinal cord due to ossification of the posterior longitudinal ligament in the cervical spine under flexion after laminectomy

Myelopathy in the cervical spine due to cervical ossification of the posterior longitudinal ligament could be induced by static compression and/or dynamic factors. It has been suggested that dynamic factors need to be considered when planning and performing the decompression surgery on patients with the ossification of the posterior longitudinal ligament. A finite element model of the C2-C7 cervical spine in the neutral position was developed and used to generate flexion and extension of the cervical spine. The segmental ossification of the posterior longitudinal ligament on the C5 was assumed, and laminectomy was performed on C4-C6 according to a conventional surgical technique. For various occupying ratios of the ossified ligament between 20% and 60%, von-Mises stresses, maximum principal strains in the spinal cord, and cross-sectional area of the cord were investigated in the pre-operative and laminectomy models under flexion, neutral position, and extension. The results were consistent with previous experimental and computational studies in terms of stress, strain, and cross-sectional area. Flexion leads to higher stresses and strains in the cord than the neutral position and extension, even after decompression surgery. These higher stresses and strains might be generated by residual compression occurring at the segment with the ossification of the posterior longitudinal ligament. This study provides fundamental information under different neck positions regarding biomechanical characteristics of the spinal cord in cervical ossification of the posterior longitudinal ligament.

Cervical ossification of the posterior longitudinal ligament: factors affecting the effect of posterior decompression

OBJECTIVE

Decompression procedures for cervical myelopathy of ossification of the posterior longitudinal ligament (OPLL) are anterior decompression with fusion, laminoplasty, and posterior decompression with fusion. Preoperative and postoperative stress analyses were performed for compression from hill-shaped cervical OPLL using 3-dimensional finite element method (FEM) spinal cord models.

METHODS

Three FEM models of vertebral arch, OPLL, and spinal cord were used to develop preoperative compression models of the spinal cord to which 10%, 20%, and 30% compression was applied; a posterior compression with fusion model of the posteriorly shifted vertebral arch; an advanced kyphosis model following posterior decompression with the spinal cord stretched in the kyphotic direction; and a combined model of advanced kyphosis following posterior decompression and intervertebral mobility. The combined model had discontinuity in the middle of OPLL, assuming the presence of residual intervertebral mobility at the level of maximum cord compression, and the spinal cord was mobile according to flexion of vertebral bodies by 5°, 10°, and 15°.

RESULTS

In the preoperative compression model, intraspinal stress increased as compression increased. In the posterior decompression with fusion model, intraspinal stress decreased, but partially persisted under 30% compression. In the advanced kyphosis model, intraspinal stress increased again. As anterior compression was higher, the stress increased more. In the advanced kyphosis  + intervertebral mobility model, intraspinal stress increased more than in the only advanced kyphosis model following decompression. Intraspinal stress increased more as intervertebral mobility increased.

CONCLUSION

In high residual compression or instability after posterior decompression, anterior decompression with fusion or posterior decompression with instrumented fusion should be considered.

Biomechanical Behaviors in Three Types of Spinal Cord Injury Mechanisms

Clinically, spinal cord injuries (SCIs) are radiographically evaluated and diagnosed from plain radiographs, computed tomography (CT), and magnetic resonance imaging. However, it is difficult to conclude that radiographic evaluation of SCI can directly explain the fundamental mechanism of spinal cord damage. The von-Mises stress and maximum principal strain are directly associated with neurological damage in the spinal cord from a biomechanical viewpoint. In this study, the von-Mises stress and maximum principal strain in the spinal cord as well as the cord cross-sectional area (CSA) were analyzed under various magnitudes for contusion, dislocation, and distraction SCI mechanisms, using a finite-element (FE) model of the cervical spine with spinal cord including white matter, gray matter, dura mater with nerve roots, and cerebrospinal fluid (CSF). A regression analysis was performed to find correlation between peak von-Mises stress/peak maximum principal strain at the cross section of the highest reduction in CSA and corresponding reduction in CSA of the cord. Dislocation and contusion showed greater peak stress and strain values in the cord than distraction. The substantial increases in von-Mises stress as well as CSA reduction similar to or more than 30% were produced at a 60% contusion and a 60% dislocation, while the maximum principal strain was gradually increased as injury severity elevated. In addition, the CSA reduction had a strong correlation with peak von-Mises stress/peak maximum principal strain for the three injury mechanisms, which might be fundamental information in elucidating the relationship between radiographic and mechanical parameters related to SCI.

Stress analysis of the cervical spinal cord: Impact of the morphology of spinal cord segments on stress

OBJECTIVE

Although there are several classifications for cervical myelopathy, these do not take differences between spinal cord segments into account. Moreover, there has been no report of stress analyses for individual segments to date.

METHODS

By using the finite element method, we constructed 3-dimensional spinal cord models comprised of gray matter, white matter, and pia mater of the second to eighth cervical vertebrae (C2-C8). We placed compression components (disc and yellow ligament) at the front and back of these models, and applied compression to the posterior section covering 10%, 20%, 30%, or 40% of the anteroposterior diameter of each cervical spinal cord segment.

RESULTS

Our results revealed that, under compression applied to an area covering 10%, 20%, or 30% of the anteroposterior diameter of the cervical spinal cord segment, sites of increased stress varied depending on the morphology of each cervical spinal cord segment. Under 40% compression, stress was increased in the gray matter, lateral funiculus, and posterior funiculus of all spinal cord segments, and stress differences between the segments were smaller.

CONCLUSION

These results indicate that, under moderate compression, sites of increased stress vary depending on the morphology of each spinal cord segment or the shape of compression components, and also that the variability of symptoms may depend on the direction of compression. However, under severe compression, the differences among the cervical spinal segments are smaller, which may facilitate diagnosis.

Mechanical Design and Analysis of a Unilateral Cervical Spinal Cord Contusion Injury Model in Non-Human Primates

Non-human primate (NHP) models of spinal cord injury better reflect human injury and provide a better foundation to evaluate potential treatments and functional outcomes. We combined finite element (FE) and surrogate models with impact data derived from in vivo experiments to define the impact mechanics needed to generate a moderate severity unilateral cervical contusion injury in NHPs (Macaca mulatta). Three independent variables (impactor displacement, alignment, and pre-load) were examined to determine their effects on tissue level stresses and strains. Mechanical measures of peak force, peak displacement, peak energy, and tissue stiffness were analyzed as potential determinants of injury severity. Data generated from FE simulations predicted a lateral shift of the spinal cord at high levels of compression (>64%) during impact. Submillimeter changes in mediolateral impactor position over the midline increased peak impact forces (>50%). Surrogate cords established a 0.5 N pre-load protocol for positioning the impactor tip onto the dural surface to define a consistent dorsoventral baseline position before impact, which corresponded with cerebrospinal fluid displacement and entrapment of the spinal cord against the vertebral canal. Based on our simulations, impactor alignment and pre-load were strong contributors to the variable mechanical and functional outcomes observed in in vivo experiments. Peak displacement of 4 mm after a 0.5N pre-load aligned 0.5-1.0 mm over the midline should result in a moderate severity injury; however, the observed peak force and calculated peak energy and tissue stiffness are required to properly characterize the severity and variability of in vivo NHP contusion injuries.

Effect of posterior decompression extent on biomechanical parameters of the spinal cord in cervical ossification of the posterior longitudinal ligament

Ossification of the posterior longitudinal ligament is a common cause of the cervical myelopathy due to compression of the spinal cord. Patients with ossification of the posterior longitudinal ligament usually require the decompression surgery, and there is a need to better understand the optimal surgical extent with which sufficient decompression without excessive posterior shifting can be achieved. However, few quantitative studies have clarified this optimal extent for decompression of cervical ossification of the posterior longitudinal ligament. We used finite element modeling of the cervical spine and spinal cord to investigate the effect of posterior decompression extent for continuous-type cervical ossification of the posterior longitudinal ligament on changes in stress, strain, and posterior shifting that occur with three different surgical methods (laminectomy, laminoplasty, and hemilaminectomy). As posterior decompression extended, stress and strain in the spinal cord decreased and posterior shifting of the cord increased. The location of the decompression extent also influenced shifting. Laminectomy and laminoplasty were very similar in terms of decompression results, and both were superior to hemilaminectomy in all parameters tested. Decompression to the extents of C3-C6 and C3-C7 of laminectomy and laminoplasty could be considered sufficient with respect to decompression itself. Our findings provide fundamental information regarding the treatment of cervical ossification of the posterior longitudinal ligament and can be applied to patient-specific surgical planning.

Influence of sagittal and axial types of ossification of posterior longitudinal ligament on mechanical stress in cervical spinal cord: A finite element analysis

BACKGROUND

There are few studies focusing on the prediction of stress distribution according to the types of ossification of the posterior longitudinal ligament, which can be fundamental information associated with clinical aspects such as the relationship between stress level and neurological symptom severity. In this study, the influence of sagittal and axial types of ossification of the posterior longitudinal ligament on mechanical stress in the cervical spinal cord was investigated.

METHODS

A three-dimensional finite element model of the cervical spine with spinal cord was developed and validated. The von Mises stresses in the cord and the reduction in cross-sectional areas and volume of the cord were investigated for various axial and sagittal types according to the occupying ratio of ossification of the posterior longitudinal ligament in the spinal canal.

FINDINGS

The influence of axial type was less than that of the sagittal type, even though the central type showed higher maximum stresses in the cord, especially for the continuous type. With a 60% occupying ratio of ossification of the posterior longitudinal ligament, the maximum stress was significantly high and the cross-sectional area of the spinal cord was reduced by more than 30% of the intact area regardless of sagittal or axial types. Finally, a higher level of sagittal extension would increase the peak cord tissue stress, which would be related to the neurological dysfunction and tissue damage.

INTERPRETATION

Quantitative investigation of biomechanical characteristics such as mechanical stress may provide fundamental information for pre-operative planning of treatment for ossification of the posterior longitudinal ligament.

Cervical ossification of the posterior longitudinal ligament: Biomechanical analysis of the influence of static and dynamic factors

OBJECTIVE

Cervical myelopathy due to ossification of the posterior longitudinal ligament (OPLL) is induced by static factors, dynamic factors, or a combination of both. We used a three-dimensional finite element method (3D-FEM) to analyze the stress distributions in the cervical spinal cord under static compression, dynamic compression, or a combination of both in the context of OPLL.

METHODS

Experimental conditions were established for the 3D-FEM spinal cord, lamina, and hill-shaped OPLL. To simulate static compression of the spinal cord, anterior compression at 10, 20, and 30% of the anterior-posterior diameter of the spinal cord was applied by the OPLL. To simulate dynamic compression, the OPLL was rotated 5°, 10°, and 15° in the flexion direction. To simulate combined static and dynamic compression under 10 and 20% anterior static compression, the OPLL was rotated 5°, 10°, and 15° in the flexion direction.

RESULTS

The stress distribution in the spinal cord increased following static and dynamic compression by cervical OPLL. However, the stress distribution did not increase throughout the entire spinal cord. For combined static and dynamic compression, the stress distribution increased as the static compression increased, even for a mild range of motion (ROM).

CONCLUSION

Symptoms may appear under static or dynamic compression only. However, under static compression, the stress distribution increases with the ROM of the responsible level and this makes it very likely that symptoms will worsen. We conclude that cervical OPLL myelopathy is induced by static factors, dynamic factors, and a combination of both.

Effect of spinal cord compression on local vascular blood flow and perfusion capacity

Spinal cord injury (SCI) can induce prolonged spinal cord compression that may result in a reduction of local tissue perfusion, progressive ischemia, and potentially irreversible tissue necrosis. Due to the combination of risk factors and the varied presentation of symptoms, the appropriate method and time course for clinical intervention following SCI are not always evident. In this study, a three-dimensional finite element fluid-structure interaction model of the cervical spinal cord was developed to examine how traditionally sub-clinical compressive mechanical loads impact spinal arterial blood flow. The spinal cord and surrounding dura mater were modeled as linear elastic, isotropic, and incompressible solids, while blood was modeled as a single-phased, incompressible Newtonian fluid. Simulation results indicate that anterior, posterior, and anteroposterior compressions of the cervical spinal cord have significantly different ischemic potentials, with prediction that the posterior component of loading elevates patient risk due to the concomitant reduction of blood flow in the arterial branches. Conversely, anterior loading compromises flow through the anterior spinal artery but minimally impacts branch flow rates. The findings of this study provide novel insight into how sub-clinical spinal cord compression could give rise to certain disease states, and suggest a need to monitor spinal artery perfusion following even mild compressive loading.

Biomechanical analysis of cervical myelopathy due to ossification of the posterior longitudinal ligament: Effects of posterior decompression and kyphosis following decompression

Cervical ossification of the posterior longitudinal ligament (OPLL) results in myelopathy. Conservative treatment is usually ineffective, thus, surgical treatment is required. One of the reasons for the poor surgical outcome following laminoplasty for cervical OPLL is kyphosis. In the present study, a 3-dimensional finite element method (3D-FEM) was used to analyze the stress distribution in preoperative, posterior decompression and kyphosis models of OPLL. The 3D-FEM spinal cord model established in this study consisted of gray and white matter, as well as pia mater. For the preoperative model, 30% anterior static compression was applied to OPLL. For the posterior decompression model, the lamina was shifted backwards and for the kyphosis model, the spinal cord was studied at 10, 20, 30, 40 and 50° kyphosis. In the preoperative model, high stress distributions were observed in the spinal cord. In the posterior decompression model, stresses were lower than those observed in the preoperative model. In the kyphosis model, an increase in the angle of kyphosis resulted in augmented stress on the spinal cord. Therefore, the results of the present study indicated that posterior decompression was effective, but stress distribution increased with the progression of kyphosis. In cases where kyphosis progresses following surgery, detailed follow-ups are required in case the symptoms worsen.

Nonlinear viscoelastic characterization of the porcine spinal cord

Although quasi-static and quasi-linear viscoelastic properties of the spinal cord have been reported previously, there are no published studies that have investigated the fully (strain-dependent) nonlinear viscoelastic properties of the spinal cord. In this study, stress relaxation experiments and dynamic cycling were performed on six fresh porcine lumbar cord specimens to examine their viscoelastic mechanical properties. The stress relaxation data were fitted to a modified superposition formulation and a novel finite ramp time correction technique was applied. The parameters obtained from this fitting methodology were used to predict the average dynamic cyclic viscoelastic behavior of the porcine cord. The data indicate that the porcine spinal cord exhibited fully nonlinear viscoelastic behavior. The average weighted root mean squared error for a Heaviside ramp fit was 2.8 kPa, which was significantly greater (p<0.001) than that of the nonlinear (comprehensive viscoelastic characterization method) fit (0.365 kPa). Further, the nonlinear mechanical parameters obtained were able to accurately predict the dynamic behavior, thus exemplifying the reliability of the obtained nonlinear parameters. These parameters will be important for future studies investigating various damage mechanisms of the spinal cord and studies developing high-resolution finite elements models of the spine.

Biomechanical analysis of cervical spondylotic myelopathy: the influence of dynamic factors and morphometry of the spinal cord

OBJECTIVE

Patients with cervical spondylotic myelopathy (CSM) have the same clinical symptoms that vary according to the degree of spinal cord compression and the cross-sectional cord shape. We used a three-dimensional finite element method (3D-FEM) to analyze the stress distributions of the spinal cord with neck extension under three cross-sectional cord shapes.

METHODS

Experimental condition for the 3D-FEM spinal cord, ligamentum flavum, and anterior compression shape (central, lateral, and diffuse types) was established. To simulate neck extension, the spinal cord was extended by 20° and the ligamentum flavum was shifted distally according to movement of the cephalad lamina.

RESULTS

The stress distribution in the spinal cord increased due to invagination of the ligamentum flavum into the neck extension. The range of stress distribution observed for the diffuse type was wider than for the central and lateral types. In addition, the stress distribution in the spinal cord was increased by the pincer movement of the ligamentum flavum and by the anterior compression of the spinal cord. The range of stress distribution observed for the diffuse type under antero-posterior compression was also wider than for the central and lateral types.

CONCLUSION

This simulation model showed that the clinical symptoms of CSM due to compression of the diffuse type may be stronger than for the central and lateral types. Therefore, careful follow-up is recommended for anterior compression of the spinal cord of diffuse type.

Maximum principal strain correlates with spinal cord tissue damage in contusion and dislocation injuries in the rat cervical spine

The heterogeneity of the primary mechanical mechanism of spinal cord injury (SCI) is not currently used to tailor treatment strategies because the effects of these distinct patterns of acute mechanical damage on long-term neuropathology have not been fully investigated. A computational model of SCI enables the dynamic analysis of mechanical forces and deformations within the spinal cord tissue that would otherwise not be visible from histological tissue sections. We created a dynamic, three-dimensional finite element (FE) model of the rat cervical spine and simulated contusion and dislocation SCI mechanisms. We investigated the relationship between maximum principal strain and tissue damage, and compared primary injury patterns between mechanisms. The model incorporated the spinal cord white and gray matter, the dura mater, cerebrospinal fluid, spinal ligaments, intervertebral discs, a rigid indenter and vertebrae, and failure criteria for ligaments and vertebral endplates. High-speed (∼ 1 m/sec) contusion and dislocation injuries were simulated between vertebral levels C3 and C6 to match previous animal experiments, and average peak maximum principal strains were calculated for several regions at the injury epicenter and at 1-mm intervals from +5 mm rostral to -5 mm caudal to the lesion. Average peak principal strains were compared to tissue damage measured previously in the same regions via axonal permeability to 10-kD fluorescein-dextran. Linear regression of tissue damage against peak maximum principal strain for pooled data within all white matter regions yielded similar and significant (p<0.0001) correlations for both contusion (R(2)=0.86) and dislocation (R(2)=0.52). The model enhances our understanding of the differences in injury patterns between SCI mechanisms, and provides further evidence for the link between principal strain and tissue damage.

The Effect of Cerebrospinal Fluid Thickness on Traumatic Spinal Cord Deformation

A spinal cord injury may lead to loss of motor and sensory function and even death. The biomechanics of the injury process have been found to be important to the neurological damage pattern, and some studies have found a protective effect of the cerebrospinal fluid (CSF). However, the effect of the CSF thickness on the cord deformation and, hence, the resulting injury has not been previously investigated. In this study, the effects of natural variability (in bovine) as well as the difference between bovine and human spinal canal dimensions on spinal cord deformation were studied using a previously validated computational model. Owing to the pronounced effect that the CSF thickness was found to have on the biomechanics of the cord deformation, it can be concluded that results from animal models may be affected by the disparities in the CSF layer thickness as well as by any difference in the biological responses they may have compared with those of humans.

Biomechanical study of the spinal cord in thoracic ossification of the posterior longitudinal ligament

BACKGROUND

Ossification of the posterior longitudinal ligament (OPLL) in the thoracic spine produces myelopathy. This is often progressive and is not affected by conservative treatment. Therefore, decompressive surgery is usually chosen.

OBJECTIVE

To conduct a stress analysis of the thoracic OPLL.

METHODS

The three-dimensional finite element spinal cord model was established. We used local ossification angle (LOA) for the degree of compression of spinal cord. LOA was the medial angle at the intersection between a line from the superior posterior margin at the cranial vertebral body of maximum OPLL to the top of OPLL with beak type, and a line from the lower posterior margin at the caudal vertebral body of the maximum OPLL to the top of OPLL with beak type. LOA 20°, LOA 25°, and LOA 30° compression was applied to the spinal cord in a preoperative model, the posterior decompressive model, and a model for the development of kyphosis.

RESULTS

In a preoperative model, at more than LOA 20° compression, high stress distributions in the spinal cord were observed. In a posterior decompressive model, the stresses were lower than in the preoperative model. In the model for development of kyphosis, high-stress distributions were observed in the spinal cord at more than LOA 20° compression.

CONCLUSIONS

Posterior decompression was an effective operative method. However, when the preoperative LOA is more than 20°, it is very likely that symptoms will worsen. If operation is performed at greater than LOA 20°, then correction of kyphosis by fixation of instruments or by forward decompression should be considered.