Implementation and validation of finite element model of skull deformation and failure response during uniaxial compression

A statistical analysis of human talar shape and bone density distribution

The shape of the talus, its internal structure, and its mechanical properties are important in determining talar behavior during loading, which may be significant for the design of surgical tools and implants. Although recent studies using statistical shape modeling have described quantitative talar external shape variation, no similar quantitative study exists to describe the density distribution of internal talar structure. The goal of this study is to quantify statistical variation in talar shape and density to benefit the design of talar implants. To this end, weight-bearing computed tomography (CT) scans of the ankle were collected in neutral, bilateral standing posture, and three-dimensional models were generated for each talus. Local density derived from the Hounsfield unit of each CT voxel was extracted. A weighted spherical harmonic analysis was performed to quantify the talar external shape. One hundred and seventy-nine volumes of interest were placed in the same relative position within each talus to quantify the talar density. Additionally, a finite element analysis (FEA) was conducted on a talus with both heterogeneous and homogeneous material properties to observe the effect of these properties on the stress and strain response. Significant differences were found in the talar density in sex and age, as well as in the stress and strain response between homogeneous and heterogeneous FEA. These differences show the importance of considering heterogeneity when examining the load response of tarsal bones.

Mechanism and microstructure based concept to predict skull fracture using a hybrid-experimental-modeling-computational approach

Cellular and tissue-scale indent/impact thresholds for different mechanisms of functional impairments to the brain would be the preferred method to predict head injuries, but a comprehensive understanding of the dominant possible injury mechanisms under multiaxial stress-states and rates is currently not available. Until then, skull fracture could serve as an indication of head injury. Therefore the ability to predict the initiation of skull fracture through finite element simulation can serve as an in silico tool for assessing the effectiveness of various head protection scenarios. For this objective, skull fracture initiation was represented with a microstructurally-inspired, mechanism-based (MIMB) failure surface assuming three different dominant mechanisms of skull failure: each element, with deformation and failure properties selected based on its microstructure, was allowed to fail either in tension, compression, or shear, corresponding to clinical linear, depressed or penetrating shear-plug failure (fracture), respectively. Microstructure-inspired a priori values for the initiation threshold of each mechanism, obtained previously from uniaxial and simple-shear experiments, were iterated and optimized for the predicted load-displacement to represent that of the corresponding indentation experiment. Element-level failure enabled the visualization of the evolution of fracture by different mechanisms. The final crack pattern at the time of macroscopic (clinically-identifiable) injury was compared between the simulation and experiment obtained through 3D tomography. Even though the timing was slightly different, the simulated prediction represented remarkably well the experimental crack pattern before the appearance of the catastrophic unstable fast crack in the experiment, thus validating the implemented hybrid-experimental-modeling-computational (HEMC) concept as a tool to predict skull fracture initiation.