Development of a geometrically accurate and adaptable finite element head model for impact simulation: the Naval Research Laboratory–Simpleware Head Model

This study demonstrates a novel model generation methodology that addresses several limitations of conventional finite element head models (FEHM). By operating chiefly in image space, new structures can be incorporated or merged, and the mesh either decimated or refined both locally and globally. This methodology is employed in the development of a highly bio-fidelic FEHM from high-resolution scan data. The model is adaptable and presented here in a form optimised for impact and blast simulations. The accuracy and feasibility of the model are successfully demonstrated against a widely used experimental benchmark in impact loading and through the investigation of potential brain injury under blast overpressure loading.     

Age-dependent material properties of the porcine cerebrum: effect on pediatric inertial head injury criteria

During growth and development, the immature central nervous system undergoes rapid alterations in constituents and structure. We hypothesize that these alterations are accompanied by changes in the mechanical properties of brain tissue which, in turn, influence the response of the brain to traumatic inertial loads. Samples of frontal cerebrum from neonatal (2–3 days) and adult pigs were harvested and tested within 3 h post-mortem. The complex shear modulus of the samples was measured in a custom-designed oscillatory shear testing device at engineering shear strain amplitudes of 2.5% or 5% from 20–200 Hz, at 25°C and 100% humidity. In this range, the elastic and viscous components of the complex shear modulus increased significantly with the development of the cerebral region of the brain. Using an idealized model of the developing head, the age-dependent material properties of brain tissue were shown to affect the mechanical response of the brain to inertial loading. This study is a first step toward developing head injury tolerance criteria specifically for the pediatric population.     

Towards reducing impact-induced brain injury: lessons from a computational study of army and football helmet pads

We use computational simulations to compare the impact response of different football and U.S. Army helmet pad materials. We conduct experiments to characterise the material response of different helmet pads. We simulate experimental helmet impact tests performed by the U.S. Army to validate our methods. We then simulate a cylindrical impactor striking different pads. The acceleration history of the impactor is used to calculate the head injury criterion for each pad. We conduct sensitivity studies exploring the effects of pad composition, geometry and material stiffness. We find that (1) the football pad materials do not outperform the currently used military pad material in militarily relevant impact scenarios; (2) optimal material properties for a pad depend on impact energy and (3) thicker pads perform better at all velocities. Although we considered only the isolated response of pad materials, not entire helmet systems, our analysis suggests that by using larger helmet shells with correspondingly thicker pads, impact-induced traumatic brain injury may be reduced.     

Filtration, haze and foam characteristics of fermented wort mediated by yeast strain

AIMS: To investigate the influence of the choice of yeast strain on the haze, shelf life, filterability and foam quality characteristics of fermented products. METHODS AND RESULTS: Twelve strains were used to ferment a chemically defined wort and hopped ale or stout wort. Fermented products were assessed for foam using the Rudin apparatus, and filterability and haze characteristics using the European Brewing Convention methods, to reveal differences in these parameters as a consequence of the choice of yeast strain and growth medium. CONCLUSIONS: Under the conditions used, the choice of strain of Saccharomyces cerevisiae effecting the primary fermentation has an impact on all of the parameters investigated, most notably when the fermentation medium is devoid of macromolecular material. SIGNIFICANCE AND IMPACT OF THE STUDY: The filtration of fermented products has a large cost implication for many brewers and wine makers, and the haze of the resulting filtrate is a key quality criterion. Also of importance to the quality of beer and some wines is the foaming and head retention of these beverages. The foam characteristics, filterability and potential for haze formation in a fermented product have long been known to be dependant on the raw materials used, as well as other production parameters. The choice of Saccharomyces cerevisiae strain used to ferment has itself been shown here to influence these parameters.     

The role of stress waves in thoracic visceral injury from blast loading: modification of stress transmission by foams and high-density materials

Materials have been applied to the thoracic wall of anaesthetised experimental animals exposed to blast overpressure to investigate the coupling of direct stress waves into the thorax and the relative contribution of compressive stress waves and gross thoracic compression to lung injury. The ultimate purpose of the work is to develop effective personal protection from the primary effects of blast overpressure--efficient protection can only be achieved if the injury mechanism is identified and characterized. Foam materials acted as acoustic couplers and resulted in a significant augmentation of the visceral injury; decoupling and elimination of injury were achieved by application of a high acoustic impedance layer on top of the foam. In vitro experiments studying stress wave transmission from air through various layers into an anechoic water chamber showed a significant increase in power transmitted by the foams, principally at high frequencies. Material such as copper or resin bonded Kevlar incorporated as a facing upon the foam achieved substantial decoupling at high frequencies--low frequency transmission was largely unaffected. An acoustic transmission model replicated the coupling of the blast waves into the anechoic water chamber. The studies suggest that direct transmission of stress waves plays a dominant role in lung parenchymal injury from blast loading and that gross thoracic compression is not the primary injury mechanism. Acoustic decoupling principles may therefore be employed to reduce the direct stress coupled into the body and thus reduce the severity of lung injury--the most simple decoupler is a high acoustic impedance material as a facing upon a foam, but decoupling layers may be optimized using acoustic transmission models. Conventional impacts producing high body wall velocities will also lead to stress wave generation and transmission--stress wave effects may dominate the visceral response to the impact with direct compression and shear contributing little to the aetiology of the injury.     

Mechanical properties of bovine pia-arachnoid complex in shear

Traumatic brain injury (TBI) has become a major public health and socioeconomic problem that affects 1.5 million Americans annually. Finite element methods have been widely used to investigate TBI mechanisms. The pia-arachnoid complex (PAC) covering the brain plays an important role in the mechanical response of the brain during impact or inertial loading. Existing finite element brain models have tended to oversimplify the response of the PAC due to a lack of accurately defined material properties of this structure, possibly resulting in a loss of accuracy in the model predictions. The objectives of this study were to experimentally determine the material properties of the PAC under shear loading. Bovine PAC was selected in the current study in view of its availability and comparability with previous studies. Tangential shear tests were conducted at 0.8, 7.3, and 72 s(-1). The mean shear moduli were 11.73, 20.04, and 22.37 kPa at the three strain rates tested. The ultimate stress, at the three strain rates, was 9.21, 17.01, and 22.26 kPa, while the ultimate strain was 1.52, 1.58, and 1.81. Results from the current study provide essential information to properly model the PAC membrane, an important component in the skull/brain interface, in a computational model of the human/animal head. Such an improved representation of the in vivo skull/brain interface will enhance future studies investigating brain injury mechanisms under various loading conditions.     

Inertia effects on characterization of dynamic response of brain tissue

Modeling and simulation of traumatic brain injury (TBI) resulted from collision or blast loading requires characterization of mechanical response over a wide range of loading rates under valid testing conditions. In this study, mechanical response of fresh bovine brain tissue was studied using the two modified Kolsky bar techniques. Radial deformation behavior of annular specimens, which are typically used to characterize the dynamic uniaxial compressive response of biological tissues, was examined using a modified Kolsky bar and a high speed camera to collect images while the specimen deforms at an axial strain rate of 2000s(-1). The high-speed images revealed inhomogeneous specimen deformation possibly brought about by radial inertia and causing a multi-axial stress state. To acquire valid stress-strain results that can be used to produce constitutive behavior of the soft materials, a novel torsion technique was developed to obtain pure shear response at dynamic loading rates. Experimental results show clear differences in the material response using the two methods. These results indicate that the previously demonstrated annular specimen geometry aimed at reducing inertia induced stress components for high rate soft materials uniaxial-compressive testing may still possess a significant component of radial inertia induced radial stress which consequently caused the observed inhomogeneous deformation in brain tissue test samples.     

Mechanics of blast loading on the head models in the study of traumatic brain injury using experimental and computational approaches

Blast waves generated by improvised explosive devices can cause mild, moderate to severe traumatic brain injury in soldiers and civilians. To understand the interactions of blast waves on the head and brain and to identify the mechanisms of injury, compression-driven air shock tubes are extensively used in laboratory settings to simulate the field conditions. The overall goal of this effort is to understand the mechanics of blast wave–head interactions as the blast wave traverses the head/brain continuum. Toward this goal, surrogate head model is subjected to well-controlled blast wave profile in the shock tube environment, and the results are analyzed using combined experimental and numerical approaches. The validated numerical models are then used to investigate the spatiotemporal distribution of stresses and pressure in the human skull and brain. By detailing the results from a series of careful experiments and numerical simulations, this paper demonstrates that: (1) Geometry of the head governs the flow dynamics around the head which in turn determines the net mechanical load on the head. (2) Biomechanical loading of the brain is governed by direct wave transmission, structural deformations, and wave reflections from tissue–material interfaces. (3) Deformation and stress analysis of the skull and brain show that skull flexure and tissue cavitation are possible mechanisms of blast-induced traumatic brain injury.     

The influence of impact force redistribution and redirection on maximum principal strain for helmeted head impacts

Abstract

Sporting helmets with linear attenuating strategies are proficient at reducing the risk of traumatic brain injury. However, the continued high incidence of concussion in American football, has led researchers to investigate novel helmet liner strategies. These strategies typically supplement existing technologies by adding or integrating head-helmet decoupling mechanisms. Decoupling strategies aim to redirect or redistribute impact force around the head, reducing impact energy transferred to the brain. This results in decreased brain tissue strain, which is beneficial in injury risk reduction due to the link between tissue strain and concussive injury.

The purpose of this study was to mathematically demonstrate the effect of ten cases, representing theoretical redirection and redistribution helmet liner strategies, on brain tissue strain resulting from impacts to the head. The kinematic response data from twenty head impacts collected in the laboratory was mathematically modified to represent the altered response of the ten different cases and used as input parameters to determine the effect on maximum principal strain (MPS) values, calculated using finite element modeling. The results showed that a reduced dominant coordinate component (contributes the greatest to resultant) of rotational acceleration decreased maximum principal strain in American football helmets. The study theoretically demonstrates that liner strategies, if applied correctly, can influence brain motion, reduce brain tissue strain, and could decrease injury risk in an American football helmet.     

Physical model simulations of brain injury in the primate

Diffuse brain injuries resulting from non-impact rotational acceleration are investigated with the aid of physical models of the skull-brain structure. These models provide a unique insight into the relationship between the kinematics of head motion and the associated deformation of the surrogate brain material. Human and baboon skulls filled with optically transparent surrogate brain tissue are subjected to lateral rotations like those shown to produce diffuse injury to the deep white matter in the brain of the baboon. High-speed cinematography captures the deformations of the grids embedded within the surrogate brain tissue during the applied load. The overall deformation pattern is compared to the pathological portrait of diffuse brain injury as determined from animal studies and autopsy reports. Shear strain and pathology spatial distributions mirror each other. Load levels and resulting surrogate brain tissue deformations are related from one species to the other. Increased primate brain mass magnified the strain amplified without significantly altering the spatial distribution. An empirically-derived value for a critical shear strain associated with the onset of severe diffuse axonal injury in primates is determined, assuming constitutive similarity between baboon and human brain tissue. The primate skull physical model data and the critical shear strain associated with the threshold for severe diffuse axonal injury were used to scale data obtained from previous studies to man, and thus derive a diffuse axonal injury tolerance for rotational acceleration for humans.     

Brain tissue analysis of impacts to American football helmets

Concussion in American football is a prevalent concern. Research has been conducted examining frequencies, location, and thresholds for concussion from impacts. Little work has been done examining how impact location may affect risk of concussive injury. The purpose of this research was to examine how impact site on the helmet and type of impact, affects the risk of concussive injury as quantified using finite element modelling of the human head and brain. A linear impactor was used to impact a helmeted Hybrid III headform in several locations and using centric and non-centric impact vectors. The resulting dynamic response was used as input for the Wayne State Brain Injury Model to determine the risk of concussive injury by utilizing maximum principal strain as the predictive variable. The results demonstrated that impacts that occur primarily to the side of the head resulted in higher magnitudes of strain in the grey and white matter, as well as the brain stem. Finally, commonly worn American football helmets were used in this research and significant risk of injury was incurred for all impacts. These results suggest that improvements in American football helmets are warranted, in particular for impacts to the side of the helmet.     

Finite element modelling of equestrian helmet impacts exposes the need to address rotational kinematics in future helmet designs

Jockey head injuries, especially concussions, are common in horse racing. Current helmets do help to reduce the severity and incidences of head injury, but the high concussion incidence rates suggest that there may be scope to improve the performance of equestrian helmets. Finite element simulations in ABAQUS/Explicit were used to model a realistic helmet model during standard helmeted rigid headform impacts and helmeted head model University College Dublin Brain Trauma Model (UCDBTM) impacts.

Current helmet standards for impact determine helmet performance based solely on linear acceleration. Brain injury-related values (stress and strain) from the UCDBTM showed that a performance improvement based on linear acceleration does not imply the same improvement in head injury-related brain tissue loads. It is recommended that angular kinematics be considered in future equestrian helmet standards, as angular acceleration was seen to correlate with stress and strain in the brain.     

Effect of boundary conditions,impact loading and hydraulic stiffening on femoral fracture strength

Patient specific quantitative CT (QCT) imaging data together with the finite element (FE) method may provide an accurate prediction of a patient's femoral strength and fracture risk. Although numerous FE models investigating femoral fracture strength have been published, there is little consent on the effect of boundary conditions, dynamic loading and hydraulic strengthening due to intra-medullary pressure on the predicted fracture strength. We developed a QCT-derived FE model of a proximal femur that included node-specific modulus assigned based on the local bone density. The effect of three commonly used boundary conditions published in literature were investigated by comparing the resulting strain field due to an applied fracture load. The models were also augmented with viscoelastic material properties and subject to a realistic impact load profile to determine the effect of dynamic loads on the strain field. Finally, the effect of hydraulic strengthening was investigated by including node specific permeability and performing a coupled pore diffusion and stress analysis of the FE model. Results showed that all boundary conditions yield the same strain field patterns, but peak strains were 22% lower and fracture load was 18% higher when loaded at the greater trochanter than when loaded at the femoral head. Comparison of the dynamic models showed that material viscoelasticity was important, but inertial effects (vibration and shock) were not. Finally, pore pressure changes did not cause significant hydraulic strengthening of bone under fall impact loading.     

Peak linear and rotational acceleration magnitude and duration effects on maximum principal strain in the corpus callosum for sport impacts

Concussion has been linked to the presence of injurious strains in the brain tissues. Research investigating severe brain injury has reported that strains in the brain may be affected by two parameters: magnitude of the acceleration, and duration of that acceleration. However, little is known how this relationship changes in terms of creating risk for brain injury for magnitudes and durations of acceleration common in sporting environments. This has particular implications for the understanding and prevention of concussive risk of injury in sporting environments. The purpose of this research was to examine the interaction between linear and rotational acceleration and duration on maximum principal strain in the brain tissues for loading conditions incurred in sporting environments. Linear and rotational acceleration loading curves of magnitudes and durations similar to those from impact in sport were used as input to the University College Brain Trauma Model and maximum principal strain (MPS) was measured for the different curves. The results demonstrated that magnitude and duration do have an effect on the strain incurred by the brain tissue. As the duration of the acceleration increases, the magnitude required to achieve strains reflecting a high risk of concussion decreases, with rotational acceleration becoming the dominant contributor. The magnitude required to attain a magnitude of MPS representing risk of brain injury was found to be as low as 2500 rad/s² for impacts of 10–15 ms; indicating that interventions to reduce the risk of concussion in sport must consider the duration of the event while reducing the magnitude of acceleration the head incurs.     

A novel surrogate lung material for impact studies: Development and testing procedures

This work focuses on the development of a surrogate lung material (SLM) that reproduces the dynamic response of a human lung under various loading conditions and also allows for the analysis of the extent and distribution of damage. The SLM consists of polyurethane foam used to mimic the spongy lung tissue and fluid-filled gelatine microcapsules used to simulate the damage of alveoli.The bursting pressure of the microcapsules was investigated by conducting low and high rate compression tests on single microcapsules. A bursting pressure of around 5 bar was measured which is comparable to the reported lung overpressure at injury level.Low and high rate compression tests were conducted on the SLMs. From the measured mechanical properties and mass density, the stress wave speed was calculated and found to be well in the range of the reported values for human lungs (16–70 m/s).In order to study the extent and distribution of damage in the SLMs, as represented by burst microcapsules, a CT scan analysis was carried out before and after the impacts. The CT scan results clearly demonstrated the magnitude and distribution of damage within the specimen. The results are then compared to the Bowen curves, the most often used criteria for predicting blast injuries in humans. An excellent agreement was found between the observed damage in the surrogate lungs and the expected damage in real human lungs.In general, the SLM showed similar stress wave speed, bursting pressure and damage to that of the real lungs.     

A computational study of influence of helmet padding materials on the human brain under ballistic impacts

The results of a computational study of a helmeted human head are presented in this paper. The focus of the work is to study the effects of helmet pad materials on the level of acceleration, inflicted pressure and shear stress in a human brain model subjected to a ballistic impact. Four different closed cell foam materials, made of expanded polystyrene and expanded polypropylene, are examined for the padding material. It is assumed that bullets cannot penetrate the helmet shell. Finite element modelling of the helmet, padding system, head and head components is used for this dynamic nonlinear analysis. Appropriate contacts and conditions are applied between the different components of the head, as well as between the head and the pads, and the pads and the helmet. Based on the results of simulations in this work, it is concluded that the stiffness of the foam has a prominent role in reducing the level of the transferred load to the brain. A pad that is less stiff is more efficient in absorbing the impact energy and reducing the sudden acceleration of the head and consequently lowers the brain injury level. Using the pad with the least stiffness, the influence of the angle of impacts as well as the locations of the ballistic strike is studied.     

Current playground surface test standards underestimate brain injury risk for children

Playgrounds surface test standards have been introduced to reduce the number of fatal and severe injuries. However, these test standards have several simplifications to make it practical, robust and cost-effective, such as the head is represented with a hemisphere, only the linear kinematics is evaluated and the body is excluded. Little is known about how these simplifications may influence the test results. The objective of this study was to evaluate the effect of these simplifications on global head kinematics and head injury prediction for different age groups. The finite element human body model PIPER was used and scaled to seven different age groups from 1.5 up to 18 years old, and each model was impacted at three different playground surface stiffness and three head impact locations. All simulations were performed in pairs, including and excluding the body. Linear kinematics and skull bone stress showed small influence if excluding the body while head angular kinematics and brain tissue strain were underestimated by the same simplification. The predicted performance of the three different playground surface materials, in terms of head angular kinematics and brain tissue strain, was also altered when including the body. A body and biofidelic neck need to be included, together with suitable head angular kinematics based injury thresholds, in future physical or virtual playground surface test standards to better prevent brain injuries.     

Computational Study of Human Head Response to Primary Blast Waves of Five Levels from Three Directions

Human exposure to blast waves without any fragment impacts can still result in primary blast-induced traumatic brain injury (bTBI). To investigate the mechanical response of human brain to primary blast waves and to identify the injury mechanisms of bTBI, a three-dimensional finite element head model consisting of the scalp, skull, cerebrospinal fluid, nasal cavity, and brain was developed from the imaging data set of a human female. The finite element head model was partially validated and was subjected to the blast waves of five blast intensities from the anterior, right lateral, and posterior directions at a stand-off distance of one meter from the detonation center. Simulation results show that the blast wave directly transmits into the head and causes a pressure wave propagating through the brain tissue. Intracranial pressure (ICP) is predicted to have the highest magnitude from a posterior blast wave in comparison with a blast wave from any of the other two directions with same blast intensity. The brain model predicts higher positive pressure at the site proximal to blast wave than that at the distal site. The intracranial pressure wave invariably travels into the posterior fossa and vertebral column, causing high pressures in these regions. The severities of cerebral contusions at different cerebral locations are estimated using an ICP based injury criterion. Von Mises stress prevails in the cortex with a much higher magnitude than in the internal parenchyma. According to an axonal injury criterion based on von Mises stress, axonal injury is not predicted to be a cause of primary brain injury from blasts.     

Infant skull and suture properties: measurements and implications for mechanisms of pediatric brain injury

The mechanical properties of the adult human skull are well documented, but little information is available for the infant skull. To determine the age-dependent changes in skull properties, we tested human and porcine infant cranial bone in three-point bending. The measurement of elastic modulus in the human and porcine infant cranial bone agrees with and extends previous published data [McPherson, G. K., and Kriewall, T. J. (1980), J. Biomech., 13, pp. 9-16] for human infant cranial bone. After confirming that the porcine and human cranial bone properties were comparable, additional tensile and three-point bending studies were conducted on porcine cranial bone and suture. Comparisons of the porcine infant data with previously published adult human data demonstrate that the elastic modulus, ultimate stress, and energy absorbed to failure increase, and the ultimate strain decreases with age for cranial bone. Likewise, we conclude that the elastic modulus, ultimate stress, and energy absorbed to failure increase with age for sutures. We constructed two finite element models of an idealized one-month old infant head, one with pediatric and the other adult skull properties, and subjected them to impact loading to investigate the contribution of the cranial bone properties on the intracranial tissue deformation pattern. The computational simulations demonstrate that the comparatively compliant skull and membranous suture properties of the infant brain case are associated with large cranial shape changes, and a more diffuse pattern of brain distortion than when the skull takes on adult properties. These studies are a fundamental initial step in predicting the unique mechanical response of the pediatric skull to traumatic loads associated with head injury and, thus, for defining head injury thresholds for children.