Correlative analysis of head kinematics and brain’s tissue response: a computational approach toward understanding the mechanisms of blast TBI

Abstract

Upon impingement of blast waves on the head, stress waves generated at the interface of the skull are transferred into the cranium and the brain tissue and may cause mild to severe blast traumatic brain injury. The intensity of the shock front, defined by the blast overpressure (BoP), that is, the blast-induced peak static overpressure, significantly affects head kinematics as well as the tissue responses of the brain. While evaluation of global linear and rotational accelerations may be feasible, an experimental determination of dynamic responses of the brain in terms of intracranial pressure (ICP), maximum shear stress (MSS), and maximum principal strain (MPS) is almost impossible. The main objective of this study is to investigate possible correlations between head accelerations and the brain’s ICP, MSS, and MPS. To this end, three different blasts were simulated by modeling the detonation of 70, 200, and 500 g of TNT at a fixed distance from the head, corresponding to peak BoPs of 0.52, 1.2, and 2 MPa, respectively. A nonlinear multi-material finite element algorithm was implemented in the LS-DYNA explicit solver. Fluid–solid interaction between the blast waves and head was modeled using a penalty-based method. Strong correlations were found between the brain’s dynamic responses and both global linear and rotational accelerations at different blast intensities ( $$R^{2 }\ge $$ 98%), implying that global kinematic parameters of the head might be strong predictors of brain tissue biomechanical parameters.