A comparison of Cartesian-only vs. Cartesian-spherical hybrid coordinates for statistical shape modeling in the lumbar spine.

OBJECTIVE: The purpose of this study was to compare two methods for quantifying differences in geometric shapes of human lumbar vertebra using statistical shape modeling (SSM). METHODS: A novel 3D implementation of a previously published 2D, nonlinear SSM was implemented and compared to a commonly used, Cartesian method of SSM. The nonlinear method, or Hybrid SSM, and Cartesian SSM were applied to lumbar vertebra shapes from a cohort of 18 full lumbar triangle meshes derived from CT scans. The comparison included traditional metrics for cumulative variance, generality, and specificity and results from application-based biomechanics using finite element simulation. RESULTS: The Hybrid SSM has less compactness - likely due to the increased number of mathematical constraints in the SSM formulation. Similar results were found between methods for specificity and generality. Compared to the previously validated, manually-segmented FE model, both SSM methods produced similar and agreeable results. CONCLUSION: Visual, statistical, and biomechanical findings did not convincingly support the superiority of the Hybrid SSM over the simpler Cartesian SSM. SIGNIFICANCE: This work suggests that, of the two methods compared, the Cartesian SSM is adequate to capture the variations in shape of the posterior spinal structures for biomechanical modeling applications.

Investigation of traumatic brain injuries using the next generation of simulated injury monitor (SIMon) finite element head model.

The objective of this study was to investigate potential for traumatic brain injuries (TBI) using a newly developed, geometrically detailed, finite element head model (FEHM) within the concept of a simulated injury monitor (SIMon). The new FEHM is comprised of several parts: cerebrum, cerebellum, falx, tentorium, combined pia-arachnoid complex (PAC) with cerebro-spinal fluid (CSF), ventricles, brainstem, and parasagittal blood vessels. The model's topology was derived from human computer tomography (CT) scans and then uniformly scaled such that the mass of the brain represents the mass of a 50th percentile male's brain (1.5 kg) with the total head mass of 4.5 kg. The topology of the model was then compared to the preliminary data on the average topology derived from Procrustes shape analysis of 59 individuals. Material properties of the various parts were assigned based on the latest experimental data. After rigorous validation of the model using neutral density targets (NDT) and pressure data, the stability of FEHM was tested by loading it simultaneously with translational (up to 400 g) combined with rotational (up to 24,000 rad/s2) acceleration pulses in both sagittal and coronal planes. Injury criteria were established in the manner shown in Takhounts et al. (2003a). After thorough validation and injury criteria establishment (cumulative strain damage measure--CSDM for diffuse axonal injuries (DAI), relative motion damage measure--RMDM for acute subdural hematoma (ASDH), and dilatational damage measure--DDM for contusions and focal lesions), the model was used in investigation of mild TBI cases in living humans based on a set of head impact data taken from American football players at the collegiate level. It was found that CSDM and especially RMDM correlated well with angular acceleration and angular velocity. DDM was close to zero for most impacts due to their mild severity implying that cavitational pressure anywhere in the brain was not reached. Maximum principal strain was found to correlate well with RMDM and angular head kinematic measures. Maximum principal stress didn't correlate with any kinematic measure or injury metric. The model was then used in the investigation of brain injury potential in NHTSA conducted side impact tests. It was also used in parametric investigations of various "what if" scenarios, such as side versus frontal impact, to establish a potential link between head kinematics and injury outcomes. The new SIMon FEHM offers an advantage over the previous version because it is geometrically more representative of the human head. This advantage, however, is made possible at the expense of additional computational time.

On the Development of the SIMon Finite Element Head Model.

The SIMon (Simulated Injury Monitor) software package is being developed to advance the interpretation of injury mechanisms based on kinematic and kinetic data measured in the advanced anthropomorphic test dummy (AATD) and applying the measured dummy response to the human mathematical models imbedded in SIMon. The human finite element head model (FEHM) within the SIMon environment is presented in this paper. Three-dimensional head kinematic data in the form of either a nine accelerometer array or three linear CG head accelerations combined with three angular velocities serves as an input to the model. Three injury metrics are calculated: Cumulative strain damage measure (CSDM) - a correlate for diffuse axonal injury (DAI); Dilatational damage measure (DDM) - to estimate the potential for contusions; and Relative motion damage measure (RMDM) - a correlate for acute subdural hematoma (ASDH). During the development, the SIMon FEHM was tuned using cadaveric neutral density targets (NDT) data and further validated against the other available cadaveric NDT data and animal brain injury experiments. The hourglass control methods, integration schemes, mesh density, and contact stiffness penalty coefficient were parametrically altered to investigate their effect on the model's response. A set of numerical and physical parameters was established that allowed a satisfactory prediction of the motion of the brain with respect to the skull, when compared with the NDT data, and a proper separation of injury/no injury cases, when compared with the brain injury data. Critical limits for each brain injury metric were also established. Finally, the SIMon FEHM performance was compared against HIC15 through the use of NHTSA frontal and side impact crash test data. It was found that the injury metrics in the current SIMon model predicted injury in all cases where HIC15 was greater than 700 and several cases from the side impact test data where HIC15 was relatively small. Side impact was found to be potentially more injurious to the human brain than frontal impact due to the more severe rotational kinematics.