Finite element study of the mechanical response in spinal cord during the thoracolumbar burst fracture

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.

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.

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.

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.

A finite element model of contusion spinal cord injury in rodents

Traumatic spinal cord injuries result from high impact forces acting on the spine and are proceeded by an extensive secondary inflammatory response resulting in motor, sensory, and autonomic dysfunction. Experimental in vivo traumatic spinal cord injuries in rodents using a contusion model have been extremely useful in elucidating the underlying pathophysiology of these injuries. However, the relationship between the pathophysiology and the biomechanical factors is still not well understood. Therefore, the aim of this research is to provide a comprehensive analysis of the biomechanics of traumatic spinal cord injury in a rat contusion model. This is achieved through the development and validation of a finite element model of the thoracic rat spinal cord and subsequently simulating controlled cortical impact-induced traumatic spinal cord injury. The effects of impactor velocity, depth, and geometry on the resulting stresses and strains within the spinal cord are investigated. Our results show that increasing impactor depth results in larger stresses and strains within the spinal cord tissue as expected. Further, for the first time ever our results show that impactor geometry (spherical versus cylindrical) plays an important role in the distribution and magnitude of stresses and strains within the cord. Therefore, finite element modelling can be a powerful tool used to predict stresses and strains that occur in spinal cord tissue during trauma.

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.

Recent Advances in Drug Design and Delivery Across Biological Barriers using Computational Models

Abstract:

The systemic delivery of pharmacological substances generally exhibits several significant limitations associated with the bio-distribution of active drugs in the body. As per consequence, human body’s defense mechanisms become impediments to drug delivery. Various technologies to overcome these limitations have been evolved including computational approaches and advanced drug delivery. As the body of human has evolved to defend itself from hostile biological as well as chemical invaders, along with that these biological barriers such as ocular barriers, blood-brain barriers, intestinal and skin barriers also limit the passage of drugs across desired sites. Therefore, efficient delivery remains an utmost challenge for researchers and scientists. The present review focuses on the techniques to deliver the drugs with efficient therapeutic efficacy at the targeted sites. This review article considered the insights into main biological barriers along with the application of computational or numerical methods dealing with different barriers by determining the drug flow, temperature and various other parameters. It also summarizes the advanced implantable drug delivery system to circumvent the inherent resistance showed by these biological barriers and in turn to improve the drug delivery.

The biomechanical effect of preexisting different types of disc herniation in cervical hyperextension injury

OBJECTIVE

Preexisting severe cervical spinal cord compression is a significant risk factor in cervical hyperextension injury, and the neurological function may deteriorate after a slight force to the forehead. There are few biomechanical studies regarding the influence of pathological factors in hyperextension loading condition. The aim of this study is to analyze the effects of preexisting different types of cervical disc herniation and different degrees of compression on the spinal cord in cervical hyperextension.

METHOD

A 3D finite element (FE) model of cervical spinal cord was modeled. Local type with median herniation, local type with lateral herniation, diffuse type with median herniation, and diffuse type with lateral herniation were simulated in neutral and extention positions. The compressions which were equivalent to 10%, 20%, 30%, and 40% of the sagittal diameter of the spinal cord were modeled.

RESULTS

The results of normal FE model were consistent with those of previous studies. The maximum von Mises stresses appeared in the pia mater for all 32 loading conditions. The maximum von Mises stresses in extension position were much higher than in neutral position. In most cases, the maximum von Mises stresses in diffuse type were higher than in local type.

CONCLUSION

Cervical spinal cord with preexisting disc herniation is more likely to be compressed in hyperextension situation than in neutral position. Diffuse type with median herniation may cause more severe compression with higher von Mises stresses concentrated at the anterior horn and the peripheral white matter, resulting in acute central cord syndrome from biomechanical point of view.

Finite element analysis of wedge and biconcave deformity in four different height restoration after augmentation of osteoporotic vertebral compression fractures

Purpose

Biomechanical comparison of wedge and biconcave deformity of different height restoration after augmentation of osteoporotic vertebral compression fractures was analyzed by three-dimensional finite element analysis (FEA).

Methods

Three-dimensional finite element model (FEM) of T11-L2 segment was constructed from CT scan of elderly osteoporosis patient. The von Mises stresses of vertebrae, intervertebral disc, facet joints, displacement, and range of motion (ROM) of wedge and biconcave deformity were compared at four different heights (Genant 0–3 grade) after T12 vertebral augmentation.

Results

In wedge deformity, the stress of T12 decreased as the vertebral height in neutral position, flexion, extension, and left axial rotation, whereas increased sharply in bending at Genant 0; L1 and L2 decreased in all positions excluding flexion of L2, and T11 increased in neutral position, flexion, extension, and right axial rotation at Genant 0. No significant changes in biconcave deformity. The stress of T11-T12, T12-L1, and L1-L2 intervertebral disc gradually increased or decreased under other positions in wedge fracture, whereas L1-L2 no significant change in biconcave fracture. The utmost overall facet joint stress is at Genant 3, whereas there is no significant change under the same position in biconcave fracture. The displacement and ROM of the wedge fracture had ups and downs, while a decline in all positions excluding extension in biconcave fracture.

Conclusions

The vertebral restoration height after augmentation to Genant 0 affects the von Mises stress, displacement, and ROM in wedge deformity, which may increase the risk of fracture, whereas restored or not in biconcave deformity.

A comprehensive finite element model of surgical treatment for cervical myelopathy

BACKGROUND

Cervical myelopathy is a common and debilitating chronic spinal cord dysfunction. Treatment includes anterior and/or posterior surgical intervention to decompress the spinal cord and stabilize the spine, but no consensus has been made as to the preferable surgical intervention. The objective of this study was to develop an finite element model of the healthy and myelopathic C2-T1 cervical spine and common anterior and posterior decompression techniques to determine how spinal cord stress and strain is altered in healthy and diseased states.

METHODS

A finite element model of the C2-T1 cervical spine, spinal cord, pia, dura, cerebral spinal fluid, and neural ligaments was developed and validated against in vivo human displacement data. To model cervical myelopathy, disc herniation and osteophytes were created at the C4-C6 levels. Three common surgical interventions were then incorporated at these levels.

FINDINGS

The finite element model accurately predicted healthy and myelopathic spinal cord displacement compared to motions observed in vivo. Spinal cord strain increased during extension in the cervical myelopathy finite element model. All surgical techniques affected spinal cord stress and strain. Specifically, adjacent levels had increased stress and strain, especially in the anterior cervical discectomy and fusion case.

INTERPRETATIONS

This model is the first biomechanically validated, finite element model of the healthy and myelopathic C2-T1 cervical spine and spinal cord which predicts spinal cord displacement, stress, and strain during physiologic motion. Our findings show surgical intervention can cause increased strain in the adjacent levels of the spinal cord which is particularly worse following anterior cervical discectomy and fusion.

Dynamics of spinal cord compression with different patterns of thoracolumbar burst fractures: Numerical simulations using finite element modelling

BACKGROUND

In thoracolumbar burst fractures, spinal cord primary injury involves a direct impact and energy transfer from bone fragments to the spinal cord. Unfortunately, imaging studies performed after the injury only depict the residual bone fragments position and pattern of spinal cord compression, with little insight on the dynamics involved during traumas. Knowledge of underlying mechanisms could be helpful in determining the severity of the primary injury, hence the extent of spinal cord damage and associated potential for recovery. Finite element models are often used to study dynamic processes, but have never been used specifically to simulate different severities of thoracolumbar burst fractures.

METHODS

Previously developed thoracolumbar spine and spinal cord finite element models were used and further validated, and representative vertebral fragments were modelled. A full factorial design was used to investigate the effects of comminution of the superior fragment, presence of an inferior fragment, fragments rotation and velocity, on maximum Von Mises stress and strain, maximum major strain, and pressure in the spinal cord.

FINDINGS

Fragment velocity clearly was the most influential factor. Fragments rotation and presence of an inferior fragment increased pressure, but rotation decreased both strains outputs. Although significant for both strains outputs, comminution of the superior fragment isn't estimated to influence outputs.

INTERPRETATION

This study is the first, to the authors' knowledge, to examine a detailed spinal cord model impacted in situ by fragments from burst fractures. This numeric model could be used in the future to comprehensively link traumatic events or imaging study characteristics to known spinal cord injuries severity and potential for recovery.

Computational simulation and modeling of the blood-brain barrier pathology

In silico methods and models in the pathology of the blood-brain barrier (BBB) or also called BBB "computational pathology", are based on using mathematical approaches together with complex, high-dimensional experimental data to evaluate and predict disease-related impacts on the CNS. These computational methods and tools have been designed to deal with BBB-linked neuropathology at the molecular, cellular, tissue, and organ levels. The molecular and cellular levels mainly include molecular docking and molecular dynamics simulations (atomistic and coarse-grain) of mutated or misfolded tight junction proteins, receptors, and various BBB transporters. The tissue and organ levels encompass the mechanistic and pharmacokinetic models as well as finite-element method and pathway analyses enriched with multiple sources of raw data (e.g., in vitro and in vivo, histopathological records, "-omics", and imaging data). Overall, this review discusses comprehensive computational techniques and strategies at different levels of complexity, providing new insights and future directions for diagnosis, treatment improvement, and a deeper understanding of BBB-related neuropathological events.

Strain rate dependent behavior of the porcine spinal cord under transverse dynamic compression

The accurate description of the mechanical properties of spinal cord tissue benefits to clinical evaluation of spinal cord injuries and is a required input for analysis tools such as finite element models. Unfortunately, available data in the literature generally relate mechanical properties of the spinal cord under quasi-static loading conditions, which is not adapted to the study of traumatic behavior, as neurological tissue adopts a viscoelastic behavior. Thus, the objective of this study is to describe mechanical properties of the spinal cord up to mechanical damage, under dynamic loading conditions. A total of 192 porcine cervical to lumbar spinal cord samples were compressed in a transverse direction. Loading conditions included ramp tests at 0.5, 5 or 50 s−1 and cyclic loading at 1, 10 or 20 Hz. Results showed that spinal cord behavior was significantly influenced by strain rate. Mechanical damage occurred at 0.64, 0.68 and 0.73 strains for 0.5, 5 or 50 s−1 loadings, respectively. Variations of behavior between the tested strain rates were explained by cyclic loading results, which revealed behavior more or less viscous depending on strain rate. Also, a parameter (stress multiplication factor) was introduced to allow transcription of a stress–strain behavior curve to different strain rates. This factor was described and was significantly different for cervical, thoracic and lumbar vertebral heights, and for the strain rates evaluated in this study.

Thoracolumbar burst fractures requiring instrumented fusion: Should reducted bone fragments be removed? A retrospective study

BACKGROUND

Thoracolumbar burst fractures are common clinical entity encountered in neurosurgical practice, accounting for 10-20% of all spinal fractures. Clinical picture could be devastating due to severe neurological deficits which lead the patients dependent both socially and emotionally.

MATERIALS AND METHODS

This study compared two groups of patients who were operated because of thoracolumbar burst fracture secondary to spinal trauma in terms of neurologic deficits, degree of improvement, and radiologic measurements at one-year follow-up. The first group (group I) included the patients who underwent posterior total laminectomy, peroperative reduction of intracanal bone fragments, and posterior spinal instrumentation and the second group (group II) included the patients who underwent total laminectomy, and spinal instrumentation without reduction of free bone fragments.

RESULTS

Neither group showed significant correlation with any measurement parameter. Radiological assessments and clinical improvements did not disclosed significant difference between the two groups at one-year follow-up.

CONCLUSION

Retropulsion of free bone fragments extend the time of surgery and causes complications. This study found that there is no need to retropulse the bone fragments in the spinal canal in patients with unstable burst fractures who underwent total laminectomy and posterior long segment stabilization.

The retrospective analysis of posterior short-segment pedicle instrumentation without fusion for thoracolumbar burst fracture with neurological deficit

This study aims to investigate the efficacy of posterior short-segment pedicle instrumentation without fusion in curing thoracolumbar burst fracture. All of the 53 patients were treated with short-segment pedicle instrumentation and laminectomy without fusion, and the restoration of retropulsed bone fragments was conducted by a novel custom-designed repositor (RRBF). The mean operation time and blood loss during surgery were analyzed; the radiological index and neurological status were compared before and after the operation. The mean operation time was 93 min (range: 62-110 min) and the mean intraoperative blood loss was 452 mL in all cases. The average canal encroachment was 50.04% and 10.92% prior to the surgery and at last followup, respectively (P < 0.01). The preoperative kyphotic angle was 17.2 degree (± 6.87 degrees), whereas it decreased to 8.42 degree (± 4.99 degrees) at last followup (P < 0.01). Besides, the mean vertebral body height increased from 40.15% (± 9.40%) before surgery to 72.34% (± 12.32%) at last followup (P < 0.01). 45 patients showed 1-2 grades improvement in Frankel's scale at last followup. This technique allows for satisfactory canal clearance and restoration of vertebral body height and kyphotic angle, and it may promote the recovery of neurological function. However, further research is still necessary to confirm the efficacy of this treatment.