Finite element–based injury metrics for pulmonary contusion via concurrent model optimization

Strain-rate-dependent material properties of human lung parenchymal tissue using inverse finite element approach

Automobile crashes and blunt trauma often lead to life-threatening thoracic injuries, especially to the lung tissues. These injuries can be simulated using finite element-based human body models that need dynamic material properties of lung tissue. The strain-rate-dependent material parameters of human parenchymal tissues were determined in this study using uniaxial quasi-static (1 mm/s) and dynamic (1.6, 3, and 5 m/s) compression tests. A bilinear material model was used to capture the nonlinear behavior of the lung tissue, which was implemented using a user-defined material in LS-DYNA. Inverse mapping using genetic algorithm-based optimization of all experimental data with the corresponding FE models yielded a set of strain-rate-dependent material parameters. The bilinear material parameters are obtained for the strain rates of 0.1, 100, 300, and 500 s-1. The estimated elastic modulus increased from 43 to 153 kPa, while the toe strain reduced from 0.39 to 0.29 when the strain rate was increased from 0.1 to 500 s-1. The optimized bilinear material properties of parenchymal tissue exhibit a piecewise linear relationship with the strain rate.

Evaluating the Limits in the Biomechanics of Blunt Lung Injury

Thoracic blunt trauma is evident in up to one fifth of all hospital admissions, and is second only to head trauma in motor vehicle crashes. One of the most problematic injury mechanisms associated with blunt thoracic trauma is pulmonary contusion, occurring in up to 75% of blunt thoracic trauma cases. The source and effects of pulmonary contusion caused by blunt lung injury are not well defined, especially within the field of continuum biomechanics. This, paired with unreliable diagnostics for pulmonary contusion, leads to uncertainty in both the clinical entity and mechanics of how to predict presence of injury. There is a distinct need to combine the clinical aspects with mechanical insights through the identification and mitigation of blunt lung trauma and material testing and modeling. This is achieved through using the mechanical insights of lung tissue behavior in order to better understand the injurious mechanisms and courses of treatment of blunt-caused pulmonary contusion. This paper hopes to act as a step forward in connecting two perspectives of blunt lung injury, the clinical entity and mechanical testing and modeling, by reviewing the known literature and identifying the unknowns within the two related fields. Through a review of related literature, clinical evidence is correlated to mechanical data to gain a better understanding of what is being missed in identification and response to blunt lung injury as a whole.

Development and implementation of a time- and computationally-efficient methodology for reconstructing real-world crashes using finite element modeling to improve crash injury research investigations

Eleven Crash Injury Research and Engineering Network (CIREN) frontal crashes were reconstructed using a novel, time-efficient methodology involving a simplified vehicle model. Kinematic accuracy was assessed using novel kinematic scores between 0-1 and chest injury was assessed using literature-defined injury metric time histories. The average kinematic score across all simulations was 0.87, indicating good kinematic accuracy. Time histories for chest compression, rib strain, shoulder belt force, and steering column force discerned the most causative components of chest injury in all cases. Abbreviated Injury Scale (AIS) 2+ and AIS 3+ chest injury risk functions using belt force identified chest injury with 81.8% success.

Biomechanically Based Correlate for Localized Lung Contusion From Nonlethal Blunt Impact Projectiles

INTRODUCTION

Injury mechanics of blunt impact projectiles differ from those experienced in whole body motor vehicle collisions because the effects are localized around the point of impact, and thus, injury thresholds based upon gross chest kinematics (e.g., force, velocity) may not be applicable across impact types. Therefore, knowledge of biomechanically based tissue injury correlates for blunt impact projectiles are needed to better guide design and development of protective systems as well as assess injury risks from blunt impact projectile weapons.

MATERIALS AND METHODS

In this study, subject-specific swine finite element models were used to quantify the tissue-level stresses and strains resulting from high speed projectile impact. These tissue-level injury doses were correlated to pathology injury outcomes to produce injury risk curves for lung contusion. Details of the pathology data and finite element results are provided in Appendix 1. Survival analysis regression methods were applied to develop lung injury regression curves and a number of statistical methods were used to evaluate several biomechanical metrics as correlates to lung contusion. Uncertainty and sensitivity analyses were used to further confirm the selection of the correlate.

RESULTS

Statistical analysis revealed that normalized strain-energy density was the best correlate for prediction of lung tissue damage. Going further, normalized strain-energy density also proved to be suitable for prediction of the percentage of contused lung volume, a more meaningful medical diagnosis. As expected, peak strain-energy density is most sensitive to muscle-skin properties, as quantified through a comprehensive uncertainty and sensitivity analysis over three sets of projectile weights and speeds.

CONCLUSIONS

Normalized strain-energy density was found to be the best correlate for prediction of lung tissue damage and correlate well to extent of contused lung volume.

Modeling lung tissue dynamics and injury under pressure and impact loading

A nonlinear viscoelastic model for the lung is implemented and evaluated for high-rate loading. Principal features of the model include a closed-cell approximation of the bulk compressibility accounting for air inside the lung and a damage-injury component by which local trauma is induced by cumulative normalized internal energy and amplified by gradients of energy density. The latter feature is adapted for use in standard numerical (i.e., explicit finite element) simulations in terms of the local rate of strain energy density and the longitudinal wave speed. Injury predictions for direct loading of a block of extracted lung material, rather than the entire thorax, via pressure pulses are in reasonably close agreement with experimental observations for an extracted rabbit lung: a threshold applied pressure exists above which edema is observed experimentally, correlating with low but non-negligible damage in the numerical results. Responses to impact by cylindrical and spherical projectiles are also interrogated. Penetration depths are comparable to those observed experimentally, as is drastically increasing damage with increasing impact velocity. Damage initiates and propagates from the impact surface, with local severity of injury decreasing with distance from the impact zone, in agreement with some empirical evidence. The model predicts more severe local injury, relative to the aforementioned surface pressure loading, than what is observed experimentally. Possible reasons for the discrepancy are analyzed, and adjustments to the model, with caveats, are suggested accordingly.

UNCERTAINTIES OF IMPACT CONFIGURATION FOR NUMERICAL REPLICATIONS OF REAL-WORLD TRAUMA: A FE ANALYSIS

This study deals with free fall accident analysis involving adults, and their numerical replications using a finite element model of the human thorax. The main purpose is to determine the role of body position at impact in the thorax injury risk appearance. For this study, cases of real-world free-fall provided by an emergency department were selected and investigated. These cases involved both male and female with an age range of 20 to 63 years, who sustained accidental free-fall with both injured and uninjured cases. The examination of the patients' medical record provided helpful information to accurately perform numerical replications with the finite element model HUByx (Hermaphrodite Universal Biomechanical yx model) which was already validated for various experimental tests in the field of automobile, ballistic impacts and blast. The results of simulations at different impact location allowed highlighting the crucial influence of the body orientation in the risk of thoracic injury occurrence.

Novel Method to Track Soft Tissue Deformation by Micro-Computed Tomography: Application to the Mitral Valve

Increasing availability of micro-computed tomography (µCT) as a structural imaging gold-standard is bringing unprecedented geometric detail to soft tissue modeling. However, the utility of these advances is severely hindered without analogous enhancement to the associated kinematic detail. To this end, labeling and following discrete points on a tissue across various deformation states is a well-established approach. Still, existing techniques suffer limitations when applied to complex geometries and large deformations and strains. Therefore, we herein developed a non-destructive system for applying fiducial markers (minimum diameter: 500 µm) to soft tissue and tracking them through multiple loading conditions by µCT. Using a novel applicator to minimize adhesive usage, four distinct marker materials were resolvable from both tissue and one another, without image artifacts. No impact on tissue stiffness was observed. µCT addressed accuracy limitations of stereophotogrammetry (inter-method positional error 1.2 ± 0.3 mm, given marker diameter 1.9 ± 0.1 mm). Marker application to ovine mitral valves revealed leaflet Almansi areal strains (45 ± 4%) closely matching literature values, and provided radiographic access to previously inaccessible regions, such as the leaflet coaptation zone. This system may meaningfully support mechanical characterization of numerous tissues or biomaterials, as well as tissue-device interaction studies for regulatory standards purposes.

SIMULATION STUDY TO COMPARE THE RELATIONSHIP BETWEEN THE BLAST-INDUCED OBJECTIVE INDICES AND THE SURVIVAL RATE OF SUBJECTS WITH PRIMARY BLAST LUNG INJURY

When an injury due to blast overpressure (BOP) is generated, it is important to estimate the severity of the injury using information about the blast conditions and to supply proper treatments according to the degree of the damage. However, there have been no investigations that have tried to verify the relationship between the blast-related objective indices and the degree of blast-induced injury. In this study, the correlations between the survival rate of the subjects with BOP-induced lung damage and each of four blast-induced indices, first principal strain, first principal strain rate, first principal stress and pulmonary inner pressure, were investigated using a simplified thorax model by introducing the concept of the V ACC –V LUNG ratio graph which represents the volume ratio between the seriously-damaged meshes and the overall meshes of the thorax model in respect to each index. Experimental results demonstrated that the decay parameters of the sigmoidal curve-fitted graphs of the first principal stress are the most effective of the analyzed indices for the estimation of the survival rate in patients with blast-induced lung damage. The results have a potential clinical application to improve the efficacy of treatment for blast injury patients.

Experimental Evidence of Mechanical Isotropy in Porcine Lung Parenchyma

Pulmonary injuries are a major source of morbidity and mortality associated with trauma. Trauma includes injuries associated with accidents and falls as well as blast injuries caused by explosives. The prevalence and mortality of these injuries has made research of pulmonary injury a major priority. Lungs have a complex structure, with multiple types of tissues necessary to allow successful respiration. The soft, porous parenchyma is the component of the lung which contains the alveoli responsible for gas exchange. Parenchyma is also the portion which is most susceptible to traumatic injury. Finite element simulations are an important tool for studying traumatic injury to the human body. These simulations rely on material properties to accurately recreate real world mechanical behaviors. Previous studies have explored the mechanical properties of lung tissues, specifically parenchyma. These studies have assumed material isotropy but, to our knowledge, no study has thoroughly tested and quantified this assumption. This study presents a novel methodology for assessing isotropy in a tissue, and applies these methods to porcine lung parenchyma. Briefly, lung parenchyma samples were dissected so as to be aligned with one of the three anatomical planes, sagittal, frontal, and transverse, and then subjected to compressive mechanical testing. Stress-strain curves from these tests were statistically compared by a novel method for differences in stresses and strains at percentages of the curve. Histological samples aligned with the anatomical planes were also examined by qualitative and quantitative methods to determine any differences in the microstructural morphology. Our study showed significant evidence to support the hypothesis that lung parenchyma behaves isotropically.

Finite element modeling of blast lung injury in sheep

A detailed 3D finite element model (FEM) of the sheep thorax was developed to predict heterogeneous and volumetric lung injury due to blast. A shared node mesh of the sheep thorax was constructed from a computed tomography (CT) scan of a sheep cadaver, and while most material properties were taken from literature, an elastic-plastic material model was used for the ribs based on three-point bending experiments performed on sheep rib specimens. Anesthetized sheep were blasted in an enclosure, and blast overpressure data were collected using the blast test device (BTD), while surface lung injury was quantified during necropsy. Matching blasts were simulated using the sheep thorax FEM. Surface lung injury in the FEM was matched to pathology reports by setting a threshold value of the scalar output termed the strain product (maximum value of the dot product of strain and strain-rate vectors over all simulation time) in the surface elements. Volumetric lung injury was quantified by applying the threshold value to all elements in the model lungs, and a correlation was found between predicted volumetric injury and measured postblast lung weights. All predictions are made for the left and right lungs separately. This work represents a significant step toward the prediction of localized and heterogeneous blast lung injury, as well as volumetric injury, which was not recorded during field testing for sheep.

Finite element model prediction of pulmonary contusion in vehicle-to-vehicle simulations of real-world crashes

OBJECTIVE

Pulmonary contusion (PC) is a common chest injury following motor vehicle crash (MVC). Because this injury has an inflammatory component, studying PC in living subjects is essential. Medical and vehicle data from the Crash Injury Research and Engineering Network (CIREN) database were utilized to examine pulmonary contusion in case occupants with known crash parameters.

METHOD

The selected CIREN cases were simulated with vehicle finite element models (FEMs) with the Total HUman Model for Safety (THUMS) version 4 as the occupant. To match the CIREN crash parameters, vehicle simulations were iteratively improved to optimize maximum crush location and depth. Fifteen cases were successfully modeled with the simulated maximum crush matching the CIREN crush to within 10%. Following the simulations, stress and strain metrics for the elements within the lungs were calculated. These injury metrics were compared to patient imaging data to determine the best finite element predictor of pulmonary contusion.

RESULTS

When the thresholds were evaluated using volumetric criteria, first principal strain was the metric with the least variation in the FEM prediction of PC.

CONCLUSIONS

A preliminary threshold for maximum crush was calculated to predict a clinically significant volume of pulmonary contusion.

Driver Injury Risk Variability in Finite Element Reconstructions of Crash Injury Research and Engineering Network (CIREN) Frontal Motor Vehicle Crashes

OBJECTIVE

A 3-phase real-world motor vehicle crash (MVC) reconstruction method was developed to analyze injury variability as a function of precrash occupant position for 2 full-frontal Crash Injury Research and Engineering Network (CIREN) cases.

METHOD

Phase I: A finite element (FE) simplified vehicle model (SVM) was developed and tuned to mimic the frontal crash characteristics of the CIREN case vehicle (Camry or Cobalt) using frontal New Car Assessment Program (NCAP) crash test data. Phase II: The Toyota HUman Model for Safety (THUMS) v4.01 was positioned in 120 precrash configurations per case within the SVM. Five occupant positioning variables were varied using a Latin hypercube design of experiments: seat track position, seat back angle, D-ring height, steering column angle, and steering column telescoping position. An additional baseline simulation was performed that aimed to match the precrash occupant position documented in CIREN for each case. Phase III: FE simulations were then performed using kinematic boundary conditions from each vehicle's event data recorder (EDR). HIC15, combined thoracic index (CTI), femur forces, and strain-based injury metrics in the lung and lumbar vertebrae were evaluated to predict injury.

RESULTS

Tuning the SVM to specific vehicle models resulted in close matches between simulated and test injury metric data, allowing the tuned SVM to be used in each case reconstruction with EDR-derived boundary conditions. Simulations with the most rearward seats and reclined seat backs had the greatest HIC15, head injury risk, CTI, and chest injury risk. Calculated injury risks for the head, chest, and femur closely correlated to the CIREN occupant injury patterns. CTI in the Camry case yielded a 54% probability of Abbreviated Injury Scale (AIS) 2+ chest injury in the baseline case simulation and ranged from 34 to 88% (mean = 61%) risk in the least and most dangerous occupant positions. The greater than 50% probability was consistent with the case occupant's AIS 2 hemomediastinum. Stress-based metrics were used to predict injury to the lower leg of the Camry case occupant. The regional-level injury metrics evaluated for the Cobalt case occupant indicated a low risk of injury; however, strain-based injury metrics better predicted pulmonary contusion. Approximately 49% of the Cobalt occupant's left lung was contused, though the baseline simulation predicted 40.5% of the lung to be injured.

CONCLUSIONS

A method to compute injury metrics and risks as functions of precrash occupant position was developed and applied to 2 CIREN MVC FE reconstructions. The reconstruction process allows for quantification of the sensitivity and uncertainty of the injury risk predictions based on occupant position to further understand important factors that lead to more severe MVC injuries.

Penetrating thorax injury leads to mild systemic activation of neutrophils without inflammatory complications

INTRODUCTION

Trauma is one of the major causes of morbidity and mortality. Thoracic injuries are associated with inflammatory complications such as ARDS. The pathogenesis of this complication after pulmonary injury is incompletely understood, but neutrophils are thought to play a pivotal role. The aim of this project was to gain more insight in the role of thoracic injuries in the pathophysiological processes that link systemic neutrophil activation with inflammatory complications after trauma.

METHODS

In this prospective cohort study fifty-five patients with isolated penetrating thoracic injury were included at a level one Trauma Unit. Blood samples were analysed for neutrophil phenotype with the use of flowcytometry within 3 h of trauma and repeated six and 24 h after injury. The presence of inflammatory complications (e.g. ARDS or sepsis/septic shock) was assessed during admission, and this was related to the neutrophil phenotpe.

RESULTS

The clinical follow-up of fifty-three patients was uneventful. Only two patients developed an inflammatory complication. Within 3 h after trauma, neutrophils showed a decreased expression of FcγRII (p=0.007) and FcγRIII (p=0.001) compared to healthy individuals. After 6 h, expression of active FcγRII (p=0.017), C5aR (p=0.004) and CAECAM8 (p=0.043) increased, whereas L-selectin (p=0.002) decreased. After 24 h also CXCR-2 (CD182) expression increased compared to healthy individuals (p=0.001).

CONCLUSIONS

Penetrating thoracic trauma leads to a distinct primed activation status of circulating neutrophils within hours. In addition to activation of cells, both young and reverse migrated neutrophils are released into the circulation. This degree of systemic inflammation does not exceed a threshold of inflammation that is needed for the development of inflammatory complications like ARDS.

Investigation of pulmonary contusion extent and its correlation to crash, occupant, and injury characteristics in motor vehicle crashes

BACKGROUND

Pulmonary contusion (PC) is a leading injury in blunt chest trauma and is most commonly caused by motor vehicle crashes (MVC). To improve understanding of the relationship between insult and outcome, this study relates PC severity to crash, occupant, and injury parameters in MVCs.

METHODS

Twenty-nine subjects with PC were selected from the Crash Injury Research and Engineering Network (CIREN) database, which contains detailed crash and medical information on MVC occupants. Computed tomography scans of these subjects were segmented using a semi-automated protocol to quantify the volumetric percentage of injured tissue in each lung. Techniques were used to quantify the geometry and location of PC, as well as the location of rib fractures. Injury extent including percent PC volume and the number of rib fractures was analyzed and its relation to crash and occupant characteristics was explored.

RESULTS

Frontal and near-side crashes composed 72% of the dataset and the near-side door was the component most often associated with PC causation. The number of rib fractures increased with age and fracture patterns varied with crash type. In near-side crashes, occupant weight and BMI were positively correlated with percent PC volume and the number of rib fractures, and the impact severity was positively correlated with percent PC volume in the lung nearest the impact.

CONCLUSIONS

This study quantified PC morphology in 29 MVC occupants and examined the relationship between injury severity and crash and occupant parameters to better characterize the mechanism of injury. The results of this study may contribute to the prevention, mitigation, and treatment of PC.

Response corridors for the medial-lateral compressive stiffness of the human arm: Implications for side impact protection

The biofidelity of side impact anthropomorphic test devices (ATDs) is crucial in order to accurately predict injury risk of human occupants. Although the arm serves as a load path to the thorax, there are currently no biofidelity response requirements for the arm. The purpose of this study was to characterize the compressive stiffness of male and female arms in medial-lateral loading and develop corresponding stiffness response corridors. This was accomplished by performing a series of pendulum tests on 18 isolated post-mortem human surrogate (PMHS) arms, obtained from four male and five female surrogates, at impact velocities of 2m/s and 4m/s. Matched tests were performed on the arm of the SID-IIs ATD for comparison. The arms were oriented vertically with the medial side placed against a rigid wall to simulate loading during a side impact automotive collision. The force versus deflection response data were normalized to that of a 50th percentile male and a 5th percentile female using a new normalizing technique based on initial arm width, and response corridors were developed for each impact velocity. A correlation analysis showed that all non-normalized dependent variables (initial stiffness, secondary stiffness, peak force, and peak deflection) were highly correlated with the initial arm width and initial arm circumference. For both impact velocities the PMHS arms exhibited a considerable amount of deflection under very low force before any substantial increase in force occurred. The compression at which the force began to increase considerably was consistent with the average tolerable medial-lateral arm compression experienced by volunteers. The initial stiffness (K1), secondary stiffness (K2), peak force, and peak deflection were found to significantly increase (p<0.05) with respect to impact velocity for both the non-normalized and normalized PMHS data. Although the response of the SID-IIs arm was similar to that of the female PMHS arms for both impact velocities, the SID-IIs arm did not exhibit a considerable toe region and therefore did not fall within the response corridors for the 5th percentile female. Overall, the results of the current study could lead to improved biofidelity of side impact ATDs by providing valuable data necessary to validate the compressive response of the ATD arm independent of the global ATD thoracic response.