Using a gel/plastic surrogate to study the biomechanical response of the head under air shock loading: a combined experimental and numerical investigation

Characterization of material properties and deformation in the ANGUS phantom during mild head impacts using MRI

Traumatic brain injury (TBI) is a major health concern affecting both military and civilian populations. Despite notable advances in TBI research in recent years, there remains a significant gap in linking the impulsive loadings from a blast or a blunt impact to the clinical injury patterns observed in TBI. Synthetic head models or phantoms can be used to establish this link as they can be constructed with geometry, anatomy, and material properties that match the human brain, and can be used as an alternative to animal models. This study presents one such phantom called the Anthropomorphic Neurologic Gyrencephalic Unified Standard (ANGUS) phantom, which is an idealized gyrencephalic brain phantom composed of polyacrylamide gel. Here we mechanically characterized the ANGUS phantom using tagged magnetic resonance imaging (MRI) and magnetic resonance elastography (MRE), and then compared the outcomes to data obtained in healthy volunteers. The direct comparison between the phantom's response and the data from a cohort of in vivo human subjects demonstrate that the ANGUS phantom may be an appropriate model for bulk tissue response and gyral dynamics of the human brain under small amplitude linear impulses. However, the phantom's response differs from that of the in vivo human brain under rotational impacts, suggesting avenues for future improvements to the phantom.

Revealing the Effect of Skull Deformation on Intracranial Pressure Variation During the Direct Interaction Between Blast Wave and Surrogate Head

Intracranial pressure (ICP) during the interaction between blast wave and the head is a crucial evaluation criterion for blast-induced traumatic brain injury (bTBI). ICP variation is mainly induced by the blast wave transmission and skull deformation. However, how the skull deformation influences the ICP remains unclear, which is meaningful for mitigating bTBI. In this study, both experimental and numerical models are developed to elucidate the effect of skull deformation on ICP variation. Firstly, we performed the shock tube experiment of the high-fidelity surrogate head to measure the ICP, the blast overpressure, and the skull surface strain of specific positions. The results show that the ICP profiles of all measured points show oscillations with positive and negative change, and the variation is consistent with the skull surface strain. Further numerical analysis reveals that when the blast wave reaches the measured point, the peak overpressure transmits directly through the skull to the brain, forming the local positive ICP peak, and the impulse induces the local inward deformation of the skull. As the peak overpressure passes through, the blast impulse impacts the nearby skull supported by the soft and incompressible brain tissue and extrudes the skull outward in the initial position. The inward and outward skull deformation leads to the oscillation of ICP. These numerical analyses agree with experimental results, which explain the appearance of negative and positive ICP peaks and the synchronization of negative ICP with surface strain. The study has implications for medical injury diagnosis and protective equipment design.

A Minireview on Brain Models Simulating Geometrical, Physical, and Biochemical Properties of the Human Brain

There is a growing body of evidences that brain surrogates will be of great interest for researchers and physicians in the medical field. They are currently mainly used for education and training purposes or to verify the appropriate functionality of medical devices. Depending on the purpose, a variety of materials have been used with specific and accurate mechanical and biophysical properties, More recently they have been used to assess the biocompatibility of implantable devices, but they are still not validated to study the migration of leaching components from devices. This minireview shows the large diversity of approaches and uses of brain phantoms, which converge punctually. All these phantoms are complementary to numeric models, which benefit, reciprocally, of their respective advances. It also suggests avenues of research for the analysis of leaching components from implantable devices.

Cerebrospinal Fluid Cavitation as a Mechanism of Blast-Induced Traumatic Brain Injury: A Review of Current Debates, Methods, and Findings

Cavitation has gained popularity in recent years as a potential mechanism of blast-induced traumatic brain injury (bTBI). This review presents the most prominent debates on cavitation; how bubbles can form or exist within the cerebrospinal fluid (CSF) and brain vasculature, potential mechanisms of cellular, and tissue level damage following the collapse of bubbles in response to local pressure fluctuations, and a survey of experimental and computational models used to address cavitation research questions. Due to the broad and varied nature of cavitation research, this review attempts to provide a necessary synthesis of cavitation findings relevant to bTBI, and identifies key areas where additional work is required. Fundamental questions about the viability and likelihood of CSF cavitation during blast remain, despite a variety of research regarding potential injury pathways. Much of the existing literature on bTBI evaluates cavitation based off its prima facie plausibility, while more rigorous evaluation of its likelihood becomes increasingly necessary. This review assesses the validity of some of the common assumptions in cavitation research, as well as highlighting outstanding questions that are essential in future work.

Assessment of Compression Driven Shock Tube Designs in Replicating Free-Field Blast Conditions for Traumatic Brain Injury Studies

Compression driven shock tubes are indispensable in studies of blast-induced traumatic brain injury (bTBI). The ability of shock tubes in faithfully recreating free-field blast conditions is of enormous interest and has a direct impact on injury outcomes. Toward this end, the evolution of blast wave inside and outside of the compression driven shock tube has been studied using validated, finite element based shock tube models. Several shock tube configurations (uniform cross-section, transition, conical, suddenly expanded, and end plate) have been considered. The finite element modeling approach has been used to simulate the transient, dynamic response of blast wave propagation. The response is studied for longer durations (40-100 msec) compared with the existing literature. We demonstrate that locations inside and outside of the shock tube can generate free-field blast profile in some form, but with numerous caveats. Our results indicate that the locations inside the shock tube are affected by higher underpressure and corresponding kinetic energy yield compared with free-field blast. These effects can be minimized using optimized end plate configuration at the exit of the shock tube, yet this is accompanied by secondary loading that is not representative of the free-field blast. Blast wave profile can be tailored using transition, conical, and suddenly expanded sections. We observe oscillations in the blast wave profile for suddenly expanded configuration. Locations outside the shock tube are affected by jet-wind effects because of the sudden expansion, barring a narrow region at the exit. For the desired overpressure yield inferred in bTBI, obtaining positive phase durations of <1 msec inside the shock tube, which are sought for studies in rodents, is challenging. Overall, these results underscore that replicating free-field blast conditions using a shock tube involves tradeoffs that need to be weighed carefully and their effect on injury outcomes should be evaluated during laboratory bTBI investigations.

Mathematical modeling and computer simulation of needle insertion into soft tissue

In this study we present a kinematic approach for modeling needle insertion into soft tissues. The kinematic approach allows the presentation of the problem as Dirichlet-type (i.e. driven by enforced motion of boundaries) and therefore weakly sensitive to unknown properties of the tissues and needle-tissue interaction. The parameters used in the kinematic approach are straightforward to determine from images. Our method uses Meshless Total Lagrangian Explicit Dynamics (MTLED) method to compute soft tissue deformations. The proposed scheme was validated against experiments of needle insertion into silicone gel samples. We also present a simulation of needle insertion into the brain demonstrating the method's insensitivity to assumed mechanical properties of tissue.

Role of Interfacial Conditions on Blast Overpressure Propagation Into the Brain

The complex interfacial condition between the human brain and the skull has been difficult to emulate in a surrogate system. Surrogate head models have typically been built using a homogeneous viscoelastic material to represent the brain, but the effect of different interfacial conditions between the brain and the skull on pressure transduction into the brain during blast has not been studied. In the present work, three interfacial conditions were generated in physical surrogate human head models. The first surrogate consisted of a gel brain separated from the skull by a layer of saline solution similar in thickness to the cerebrospinal fluid (CSF) layer in the human head: the fluid interface head model. The second surrogate head had the entire cranial cavity filled with the gel: the fixed interface head model. The third surrogate head contained a space-filling gel brain wrapped in a thin plastic film: the stick-slip interface head model. The human head surrogates were evaluated in a series of frontal blast tests to characterize the effect of skull-brain interfacial conditions on overpressure propagation into the gel brains. The fixed and the stick-slip interface head models showed nearly equal peak brain overpressures. In contrast, the fluid interface head model had much higher in-brain peak overpressures than the other two models, thus representing the largest transmission of forces into the gel brain. Given that the elevated peak overpressures occurred only in the fluid interface head model, the presence of the saline layer is likely responsible for this increase. This phenomenon is hypothesized to be attributed to the incompressibility of the saline and/or the impedance differences between the materials. The fixed interface head model showed pronounced high frequency energy content relative to the other two models, implying that the fluid and the stick-slip conditions provided better dampening. The cumulative impulse energy entering the three brain models were similar, suggesting that the interface conditions do not affect the total energy transmission over the positive phase duration of a blast event. This study shows that the fidelity of the surrogate human head models would improve with a CSF-emulating liquid layer.

Feasibility of using a novel instrumented human head surrogate to measure helmet, head and brain kinematics and intracranial pressure during multidirectional impact tests

OBJECTIVES

Aim of the work is to present the feasibility of using an Instrumented Human Head Surrogate (IHHS-1) during multidirectional impacts while wearing a modern ski helmet. The IHHS-1 is intended to provide reliable and repeatable data for the experimental validation of FE models and for the experimental evaluation of modern helmets designed to enhance the degree of protection against multidirectional impacts.

DESIGN

The new IHHS-1 includes 9 triaxial MEMS accelerometers embedded in a silicone rubber brain, independently molded and presenting lobes separation and cerebellum, placed into an ABS skull filled with surrogate cerebrospinal fluid. A triaxial MEMS gyroscope is placed at the brain center of mass. Intracranial pressure can be detected by eight pressure sensors applied to the skull internal surface along a transversal plane located at the brain center of mass and two at the apex. Additional MEMS sensors positioned over the skull and the helmet allow comparison between outer and inner structure kinematics and surrogate CSF pressure behavior.

METHODS

The IHHS-1 was mounted through a Hybrid III neck on a force platform and impacted with a striker connected to a pendulum tower, with the impact energies reaching 24J. Impact locations were aligned with the brain center of mass and located in the back (sagittal axis), right (90° from sagittal axis), back/right (45°), and front right (135°) locations. Following dynamic data were collected: values of the linear accelerations and angular velocities of the brain, skull and helmet; intracranial pressures inside the skull.

RESULTS

Despite the relatively low intensity of impacts (HIC at skull max value 46), the skull rotational actions reached BrIC values of 0.33 and angular accelerations of 5216rad/s², whereas brain angular acceleration resulted between 1,44 and 2,1 times lower with similar values of BrIC.

CONCLUSIONS

The IHHS-1 is a physical head surrogate that can produce repeatable data for the interpretation of inner structures behavior during multidirectional impacts with or without helmets of different characteristics.

Effective testing of personal protective equipment in blast loading conditions in shock tube: Comparison of three different testing locations

We exposed a headform instrumented with 10 pressure sensors mounted flush with the surface to a shock wave with three nominal intensities: 70, 140 and 210 kPa. The headform was mounted on a Hybrid III neck, in a rigid configuration to eliminate motion and associated pressure variations. We evaluated the effect of the test location by placing the headform inside, at the end and outside of the shock tube. The shock wave intensity gradually decreases the further it travels in the shock tube and the end effect degrades shock wave characteristics, which makes comparison of the results obtained at three locations a difficult task. To resolve these issues, we developed a simple strategy of data reduction: the respective pressure parameters recorded by headform sensors were divided by their equivalents associated with the incident shock wave. As a result, we obtained a comprehensive set of non-dimensional parameters. These non-dimensional parameters (or amplification factors) allow for direct comparison of pressure waveform characteristic parameters generated by a range of incident shock waves differing in intensity and for the headform located in different locations. Using this approach, we found a correlation function which allows prediction of the peak pressure on the headform that depends only on the peak pressure of the incident shock wave (for specific sensor location on the headform), and itis independent on the headform location. We also found a similar relationship for the rise time. However, for the duration and impulse, comparable correlation functions do not exist. These findings using a headform with simplified geometry are baseline values and address a need for the development of standardized parameters for the evaluation of personal protective equipment (PPE) under shock wave loading.

Primary Blast Brain Injury Mechanisms: Current Knowledge, Limitations, and Future Directions

Mild blast traumatic brain injury (bTBI) accounts for the majority of brain injury in United States service members and other military personnel worldwide. The mechanisms of primary blast brain injury continue to be disputed with little evidence to support one or a combination of theories. The main hypotheses addressed in this review are blast wave transmission through the skull orifices, direct cranial transmission, skull flexure dynamics, thoracic surge, acceleration, and cavitation. Each possible mechanism is discussed using available literature with the goal of focusing research efforts to address the limitations and challenges that exist in blast injury research. Multiple mechanisms may contribute to the pathology of bTBI and could be dependent on magnitudes and orientation to blast exposure. Further focused biomechanical investigation with cadaver, in vivo, and finite element models would advance our knowledge of bTBI mechanisms. In addition, this understanding could guide future research and contribute to the greater goal of developing relevant injury criteria and mandates to protect our soldiers on the battlefield.

Viscoelastic computational modeling of the human head-neck system: Eigenfrequencies and time-dependent analysis

A subject-specific 3-dimensional viscoelastic finite element model of the human head-neck system is presented and investigated based on computed tomography and magnetic resonance biomedical images. Ad hoc imaging processing tools are developed for the reconstruction of the simulation domain geometry and the internal distribution of bone and soft tissues. Material viscoelastic properties are characterized point-wise through an image-based interpolating function used then for assigning the constitutive prescriptions of a heterogenous viscoelastic continuum model. The numerical study is conducted both for modal and time-dependent analyses, compared with similar studies and validated against experimental evidences. Spatiotemporal analyses are performed upon different exponential swept-sine wave-localized stimulations. The modeling approach proposes a generalized, patient-specific investigation of sound wave transmission and attenuation within the human head-neck system comprising skull and brain tissues. Model extensions and applications are finally discussed.

Relevance of Blood Vessel Networks in Blast-Induced Traumatic Brain Injury

Cerebral vasculature is a complex network that circulates blood through the brain. However, the role of this networking effect in brain dynamics has seldom been inspected. This work is to study the effects of blood vessel networks on dynamic responses of the brain under blast loading. Voronoi tessellations were implemented to represent the network of blood vessels in the brain. The brain dynamics in terms of maximum principal strain (MPS), shear strain (SS), and intracranial pressure (ICP) were monitored and compared. Results show that blood vessel networks significantly affected brain responses. The increased MPS and SS were observed within the brain embedded with vessel networks, which did not exist in the case without blood vessel networks. It is interesting to observe that the alternation of the ICP response was minimal. Moreover, the vessel diameter and density also affected brain dynamics in both MPS and SS measures. This work sheds light on the role of cerebral vasculature in blast-induced traumatic brain injury.

A theoretical analysis of stress wave propagation in the head under primary blast loading

Traumatic brain injury due to primary blast loading has become a signature injury in recent military conflicts. Efforts have been made to study the stress wave propagation in the head. However, the relationship of incident pressure, reflected pressure and intracranial pressure is still not clear, and the experimental findings reported in the literature are contradictory. In this article, an analytical model is developed to calculate the stress wave transfer through a multiple-layered structure which is used to mimic the head. The model predicts stress at the scalp-skull and skull-brain interfaces as the functions of reflected pressure, which is further dependent on incident pressure. A numerical model is used to corroborate the theoretical predictions. It is concluded that scalp has an amplification effect on intracranial pressure. If scalp is absent, there exists a critical incident pressure, defined as P _cr at approximately 16 kPa. When peak incident pressure σ _in is higher than 16 kPa, the pressure at the skull-brain interface is greater than σ _in; otherwise, it is lower than σ _in.

The mechanics of traumatic brain injury: a review of what we know and what we need to know for reducing its societal burden

Traumatic brain injury (TBI) is a significant public health problem, on pace to become the third leading cause of death worldwide by 2020. Moreover, emerging evidence linking repeated mild traumatic brain injury to long-term neurodegenerative disorders points out that TBI can be both an acute disorder and a chronic disease. We are at an important transition point in our understanding of TBI, as past work has generated significant advances in better protecting us against some forms of moderate and severe TBI. However, we still lack a clear understanding of how to study milder forms of injury, such as concussion, or new forms of TBI that can occur from primary blast loading. In this review, we highlight the major advances made in understanding the biomechanical basis of TBI. We point out opportunities to generate significant new advances in our understanding of TBI biomechanics, especially as it appears across the molecular, cellular, and whole organ scale.

Biomechanical responses of a pig head under blast loading: a computational simulation

A series of computational studies were performed to investigate the biomechanical responses of the pig head under a specific shock tube environment. A finite element model of the head of a 50-kg Yorkshire pig was developed with sufficient details, based on the Lagrangian formulation, and a shock tube model was developed using the multimaterial arbitrary Lagrangian-Eulerian (MMALE) approach. These two models were integrated and a fluid/solid coupling algorithm was used to simulate the interaction of the shock wave with the pig's head. The finite element model-predicted incident and intracranial pressure traces were in reasonable agreement with those obtained experimentally. Using the verified numerical model of the shock tube and pig head, further investigations were carried out to study the spatial and temporal distributions of pressure, shear stress, and principal strain within the head. Pressure enhancement was found in the skull, which is believed to be caused by shock wave reflection at the interface of the materials with distinct wave impedances. Brain tissue has a shock attenuation effect and larger pressures were observed in the frontal and occipital regions, suggesting a greater possibility of coup and contrecoup contusion. Shear stresses in the brain and deflection in the skull remained at a low level. Higher principal strains were observed in the brain near the foramen magnum, suggesting that there is a greater chance of cellular or vascular injuries in the brainstem region.

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.

Physics of IED Blast Shock Tube Simulations for mTBI Research

Shock tube experiments and simulations are conducted with a spherical gelatin filled skull–brain surrogate, in order to study the mechanisms leading to blast induced mild traumatic brain injury. A shock tube including sensor system is optimized to simulate realistic improvised explosive device blast profiles obtained from full scale field tests. The response of the skull–brain surrogate is monitored using pressure and strain measurements. Fluid–structure interaction is modeled using a combination of computational fluid dynamics (CFD) simulations for the air blast, and a finite element model for the structural response. The results help to understand the physics of wave propagation, from air blast into the skull–brain. The presence of openings on the skull and its orientation does have a strong effect on the internal pressure. A parameter study reveals that when there is an opening in the skull, the skull gives little protection and the internal pressure is fairly independent on the skull stiffness; the gelatin shear stiffness has little effect on the internal pressure. Simulations show that the presence of pressure sensors in the gelatin hardly disturbs the pressure field.