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

Brain Material Properties and Integration of Arachnoid Complex for Biofidelic Impact Response for Human Head Finite Element Model

Finite element head models offer great potential to study brain-related injuries; however, at present may be limited by geometric and material property simplifications required for continuum-level human body models. Specifically, the mechanical properties of the brain tissues are often represented with simplified linear viscoelastic models, or the material properties have been optimized to specific impact cases. In addition, anatomical structures such as the arachnoid complex have been omitted or implemented in a simple lumped manner. Recent material test data for four brain regions at three strain rates in three modes of loading (tension, compression, and shear) was used to fit material parameters for a hyper-viscoelastic constitutive model. The material model was implemented in a contemporary detailed head finite element model. A detailed representation of the arachnoid trabeculae was implemented with mechanical properties based on experimental data. The enhanced head model was assessed by re-creating 11 ex vivo head impact scenarios and comparing the simulation results with experimental data. The hyper-viscoelastic model faithfully captured mechanical properties of the brain tissue in three modes of loading and multiple strain rates. The enhanced head model showed a high level of biofidelity in all re-created impacts in part due to the improved brain-skull interface associated with implementation of the arachnoid trabeculae. The enhanced head model provides an improved predictive capability with material properties based on tissue level data and is positioned to investigate head injury and tissue damage in the future.

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.

Comparison of biomechanical parameters of two Chinese cervical spine rotation manipulations based on motion capture and finite element analysis

Objective: The purpose of this study was to obtain the stress-strain of the cervical spine structure during the simulated manipulation of the oblique pulling manipulation and the cervical rotation-traction manipulation in order to compare the mechanical mechanism of the two manipulations. Methods: A motion capture system was used to record the key kinematic parameters of operating the two manipulations. At the same time, a three-dimensional finite element model of the C0-T1 full healthy cervical spine was established, and the key kinematic parameters were loaded onto the finite element model in steps to analyze and simulate the detailed process of the operation of the two manipulations. Results: A detailed finite element model of the whole cervical spine including spinal nerve roots was established, and the validity of this 3D finite element model was verified. During the stepwise simulation of the two cervical spine rotation manipulations to the right, the disc (including the annulus fibrosus and nucleus pulposus) and facet joints stresses and displacements were greater in the oblique pulling manipulation group than in the cervical rotation-traction manipulation group, while the spinal cord and nerve root stresses were greater in the cervical rotation-traction manipulation group than in the oblique pulling manipulation group. The spinal cord and nerve root stresses in the cervical rotation-traction manipulation group were mainly concentrated in the C4/5 and C5/6 segments. Conclusion: The oblique pulling manipulation may be more appropriate for the treatment of cervical spondylotic radiculopathy, while cervical rotation-traction manipulation is more appropriate for the treatment of cervical spondylosis of cervical type. Clinicians should select cervical rotation manipulations for different types of cervical spondylosis according to the patient's symptoms and needs.

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.