Cortical bone continuum damage mechanics constitutive model with stress triaxiality criterion to predict fracture initiation and pattern.

A primary objective of finite element human body models (HBMs) is to predict response and injury risk in impact scenarios, including cortical bone fracture initiation, fracture pattern, and the potential to simulate post-fracture injury to underlying soft tissues. Current HBMs have been challenged to predict the onset of failure and bone fracture patterns owing to the use of simplified failure criteria. In the present study, a continuum damage mechanics (CDM) model, incorporating observed mechanical response (orthotropy, asymmetry, damage), was coupled to a novel phenomenological effective strain fracture criterion based on stress triaxiality and investigated to predict cortical bone response under different modes of loading. Three loading cases were assessed: a coupon level notched shear test, whole bone femur three-point bending, and whole bone femur axial torsion. The proposed material model and fracture criterion were able to predict both the fracture initiation and location, and the fracture pattern for whole bone and specimen level tests, within the variability of the reported experiments. There was a dependence of fracture threshold on finite element mesh size, where higher mesh density produced similar but more refined fracture patterns compared to coarser meshes. Importantly, the model was functional, accurate, and numerically stable even for relatively coarse mesh sizes used in contemporary HBMs. The proposed model and novel fracture criterion enable prediction of fracture initiation and resulting fracture pattern in cortical bone such that post-fracture response can be investigated in HBMs.

Validation of a Football Helmet Finite Element Model and Quantification of Impact Energy Distribution.

Head injury in contact sports can be mitigated, in part, through the enhancement of protective helmets that may be enabled by detailed finite element models. However, many contemporary helmet FE models include simplified geometry and material properties and have limited verification and validation over a representative range of impact conditions. To address these limitations, a detailed numerical model of a modern football helmet was developed, integrated with two headforms and assessed for 60 impact conditions with excellent ratings (0.79-0.93). The strain energy of the helmet components was investigated for eight impact locations and three impact speeds. In general, the helmet shell had the highest strain energy followed by the compression shocks; however, the facemask and straps had the highest strain energy for impacts involving the facemask. The component strain energy was positively correlated with the Head Injury Criterion, while the strain energy was not strongly correlated with the Brain Injury Criterion due to the dependence on rotational kinematics. This study demonstrated the applicability of a detailed football helmet finite element model to investigate a range of impact conditions and to assess energy distribution as a function of impact location and severity as a means of future helmet optimization.

Importance of asymmetry and anisotropy in predicting cortical bone response and fracture using human body model femur in three-point bending and axial rotation.

Modeling of cortical bone response and failure is critical for the prediction of Crash Induced Injuries (CII) using advanced finite element (FE) Human Body Models (HBM). Although cortical bone is anisotropic and asymmetric in tension and compression, current HBM often utilize simple isotropic, symmetric, elastic-plastic constitutive models. In this study, a 50th percentile male femur FE model was used to quantify the effect of asymmetry and anisotropy in three-point bending and axial torsion. A complete set of cortical bone mechanical properties was identified from a literature review, and the femur model was used to investigate the importance of material asymmetry and anisotropy on the failure load/moment, failure displacement/rotation and fracture pattern. All models were able to predict failure load in bending, since this was dominated by the cortical bone material tensile response. However, only the orthotropic model was able to predict the torsional response and failure moment. Only the orthotropic model predicted the fracture initiation location and fracture pattern in bending, and the fracture initiation location in torsion; however, the anticipated spiral fracture pattern was not predicted by the models for torsional loading. The results demonstrated that asymmetry did not significantly improve the prediction capability, and that orthotropic material model with the identified material properties was able to predict the kinetics and kinematics for both three-point bending and axial torsion. This will help to provide an improved method for modeling hard tissue response and failure in full HBM.