Construction and evaluation of thoracic injury risk curves for a finite element human body model in frontal car crashes

Holistic shape variation of the rib cage in an adult population

Traumatic injuries to the thorax are a common occurrence, and given the disparity in outcomes, injury risk is non-uniformly distributed within the population. Rib cage geometry, in conjunction with well-established biomechanical characteristics, is thought to influence injury tolerance, but quantifiable descriptions of adult rib cage shape as a whole are lacking. Here, we develop an automated pipeline to extract whole rib cage measurements from a large population and produce distributions of these measurements to assess variability in rib cage shape. Ten measurements of whole rib cage shape were collected from 1,719 individuals aged 25-45 years old including angular, linear, areal, and volumetric measures. The resulting pipeline produced measurements with a mean percent difference to manually collected measurements of 1.7% ± 1.6%, and the whole process takes 30 s per scan. Each measurement followed a normal distribution with a maximum absolute skew value of 0.43 and a maximum absolute excess kurtosis value of 0.6. Significant differences were found between the sexes (p < 0.001) in all except angular measures. Multivariate regression revealed that demographic predictors explain 29%-68% of the variance in the data. The angular measurements had the three lowest R2 values and were also the only three to have little correlation with subject stature. Unlike other measures, rib cage height had a negative correlation with BMI. Stature was the dominant demographic factor in predicting rib cage height, coronal area, sagittal area, and volume. Subject weight was the dominant demographic factor for rib cage width, depth, axial area, and angular measurements. Age was minimally important in this cohort of adults from a narrow age range. Individuals of similar height and weight had average rib cage measurements near the regression predictions, but the range of values across all subjects encompassed a large portion of their respective distributions. Our findings characterize the variability in adult rib cage geometry, including the variation within narrow demographic criteria. In future work, these can be integrated into computer aided engineering workflows to assess the influence of whole rib cage shape on the biomechanics of the adult human thorax.

Mechanical response of human thoracic spine ligaments under quasi-static loading: An experimental study

PURPOSE

This study aimed to investigate the geometrical and mechanical properties of human thoracic spine ligaments subjected to uniaxial quasi-static tensile test.

METHODS

Four human thoracic spines, obtained through a body donation program, were utilized for the study. The anterior longitudinal ligament (ALL), posterior longitudinal ligament (PLL), capsular ligament (CL), ligamenta flava (LF), and the interspinous ligament and supraspinous ligament complex (ISL + SSL), were investigated. The samples underwent specimen preparation, including dissection, cleaning, and reinforcement, before being immersed in epoxy resin. Uniaxial tensile tests were performed using a custom-designed mechanical testing machine equipped with an environmental chamber (T = 36.6 °C; humidity 95%). Then, the obtained tensile curves were averaged preserving the characteristic regions of typical ligaments response.

RESULTS

Geometrical and mechanical properties, such as initial length and width, failure load, and failure elongation, were measured. Analysis of variance (ANOVA) revealed significant differences among the ligaments for all investigated parameters. Pairwise comparisons using Tukey's post-hoc test indicated differences in initial length and width. ALL and PLL exhibited higher failure forces compared to CL and LF. ALL and ISL + SSL demonstrated biggest failure elongation. Comparisons with other studies showed variations in initial length, failure force, and failure elongation across different ligaments. The subsystem (Th1 - Th6 and Th7 - Th12) analysis revealed increases in initial length, width, failure force, and elongation for certain ligaments.

CONCLUSIONS

Variations of both the geometric and mechanical properties of the ligaments were noticed, highlighting their unique characteristics and response to tensile force. Presented results extend very limited experimental data base of thoracic spine ligaments existing in the literature. The obtained geometrical and mechanical properties can help in the development of more precise human body models (HBMs).

Development and validation of lumbar spine finite element model

The functional biomechanics of the lumbar spine have been better understood by finite element method (FEM) simulations. However, there are still areas where the behavior of soft tissues can be better modeled or described in a different way. The purpose of this research is to develop and validate a lumbar spine section intended for biomechanical research. A FE model of the 50th percentile adult male (AM) Total Human Model for Safety (THUMS) v6.1 was used to implement the modifications. The main modifications were to apply orthotropic material properties and nonlinear stress-strain behavior for ligaments, hyperelastic material properties for annulus fibrosus and nucleus pulposus, and the specific content of collagenous fibers in the annulus fibrosus ground substance. Additionally, a separation of the nucleus pulposus from surrounding bones and tissues was implemented. The FE model was subjected to different loading modes, in which intervertebral rotations and disc pressures were calculated. Loading modes contained different forces and moments acting on the lumbar section: axial forces (compression and tension), shear forces, pure moments, and combined loading modes of axial forces and pure moments. The obtained ranges of motion from the modified numerical model agreed with experimental data for all loading modes. Moreover, intradiscal pressure validation for the modified model presented a good agreement with the data available from the literature. This study demonstrated the modifications of the THUMS v6.1 model and validated the obtained numerical results with existing literature in the sub-injurious range. By applying the proposed changes, it is possible to better model the behavior of the human lumbar section under various loads and moments.

Evaluation of injury threshold from the number of rib fracture for predicting pulmonary injuries in blunt chest trauma

Background

Blunt chest trauma is a common presentation in emergency departments. The relationship between bone fractures and organ injuries has not been studied in detail. The purpose of this study was to examine the degree of external force represented by the number of rib fractures that causes lung injury in blunt chest trauma.

Patients and methods

This study was performed retrospectively using trauma patients who received medical examinations in a single university hospital emergency center between April 2015 and March 2020. We examined the association between the number of rib fractures and pulmonary damage using multivariable regression analysis and considered the relationship between rib fracture location and each type of lung injury.

Results

A total of 317 patients were included. The mean age was 63.1 years, 65.0% were male, and traffic accidents were the most common mechanism of injury (55.8%). The number of mean rib fractures was 4.0, and the mean Injury Severity Score was 11.3. The number of rib fractures was associated with an increased risk of pulmonary injuries: pulmonary contusion (odds ratio [OR] 1.30, 95% confidence interval [CI] 1.14-1.48, p < 0.05); hemothorax (OR 1.22, 95% CI 1.08-1.38, p < 0.05); pneumothorax (OR 1.15, 95% CI 1.02-1.30, p < 0.05); and hemopneumothorax (OR 1.14, 95% CI 1.01-1.28, p < 0.05). In addition, bilateral rib fractures were associated with fractures of the superior ribs more often and more severely, but were not associated with the occurrence of each type of lung injury.

Conclusion

The number of rib fractures was associated with an increased risk of pulmonary injuries. In addition, the type of pulmonary injury could be predicted from the number of rib fractures in blunt chest trauma.

Assessment of the sensitivity of thoracic injury criteria to subject-specific characteristics using human body models

Introduction: Chest deformation has been proposed as the best predictor of thoracic injury risk in frontal impacts. Finite Element Human Body Models (FE-HBM) can enhance the results obtained in physical crash tests with Anthropometric Test Devices (ATD) since they can be exposed to omnidirectional impacts and their geometry can be modified to reflect specific population groups. This study aims to assess the sensitivity of two thoracic injury risk criteria (PC Score and Cmax) to several personalization techniques of FE-HBMs. Methods: Three 30° nearside oblique sled tests were reproduced using the SAFER HBM v8 and three personalization techniques were applied to this model to evaluate the influence on the risk of thoracic injuries. First, the overall mass of the model was adjusted to represent the weight of the subjects. Second, the model anthropometry and mass were modified to represent the characteristics of the post-mortem human subjects (PMHS). Finally, the spine alignment of the model was adapted to the PMHS posture at t = 0 ms, to conform to the angles between spinal landmarks measured in the PMHS. The following two metrics were used to predict three or more fractured ribs (AIS3+) of the SAFER HBM v8 and the effect of personalization techniques: the maximum posterior displacement of any studied chest point (Cmax), and the sum of the upper and lower deformation of selected rib points (PC score). Results: Despite having led to statistically significant differences in the probability of AIS3+ calculations, the mass-scaled and morphed version provided, in general, lower values for injury risk than the baseline model and the postured version being the latter, which exhibited the better approximation to the PMHS tests in terms of probability of injury. Additionally, this study found that the prediction of AIS3+ chest injuries based on PC Score resulted in higher probability values than the prediction based on Cmax for the loading conditions and personalization techniques analyzed within this study. Discussion: This study could demonstrate that the personalization techniques do not lead to linear trends when they are used in combination. Furthermore, the results included here suggest that these two criteria will result in significantly different predictions if the chest is loaded more asymmetrically.

Identifying Vulnerable Impact Locations to Reduce the Occurrence of Deadly Commotio Cordis Events in Children's Baseball: A Computational Approach

Commotio cordis is the second leading cause of sudden cardiac death in young athletes. Currently available chest protectors on the market are ineffective in preventing cases of commotio cordis in young athletes who play baseball. This study focused on using contour maps to identify specific baseball impact locations to the chest that may result in instances of commotio cordis to children during baseball games. By identifying these vulnerable locations, we may design and develop chest protectors that can provide maximum protection to prevent commotio cordis in young athletes. Simulation cases were run using the validated CHARM-10 chest model, a detailed finite element model representing an average 10-year-old child's chest. A baseball model was developed in company with the chest model, and then used to impact the chest at different locations. A 7 × 8 impact location matrix was designed with 56 unique baseball impact simulations. Left ventricle strain and pressure, reaction force between the baseball and chest, and rib deformations were analyzed. Left ventricle strain was highest from baseball impacts directly over the left ventricle (0.34) as well as impacts slightly lateral and superior to the cardiac silhouette (0.34). Left ventricle pressure was highest with impacts directly over the left ventricle (82.94 kPa). We have identified the most dangerous impact locations resulting in high left ventricle strain and pressure. This novel study provided evidence of where to emphasize protective materials for establishing effective chest protectors that will minimize instances of commotio cordis in young athletes.

Biomechanics of Lumbar Spine Injury in Road Barrier Collision-Finite Element Study

Literature and field data from CIREN database have shown that lumbar spine injuries occur during car crashes. There are multiple hypotheses regarding how they occur; however, there is no biomechanical explanation for these injuries during collisions with road safety barriers (RSBs). Therefore, the objective of this study was to investigate the mechanics of vertebral fractures during car collisions with concrete RSBs. The finite element method was used for the numerical simulations. The global model of the car collision with the concrete RSB was created. The lumbar spine kinematics were extracted from the global simulation and then applied as boundary conditions to the detailed lumbar spine model. The results showed that during the collision, the occupant was elevated, and then dropped during the vehicle landing. This resulted in axial compression forces 2.6 kN with flexion bending moments 34.7 and 37.8 Nm in the L2 and L3 vertebrae. It was shown that the bending moment is the result of the longitudinal force on the eccentricity. The lumbar spine index for the L1-L5 section was 2.80, thus indicating a lumbar spine fracture. The minimum principal strain criterion of 7.4% and damage variable indicated L2 and L3 vertebrae and the inferior part of L1, as those potentially prone to fracture. This study found that lumbar spine fractures could occur as a consequence of vehicle landing during a collision with a concrete RSB mostly affecting the L1-L3 lumbar spine section. The fracture was caused by a combination of axial forces and flexion bending moments.

An Improved Method for Developing Injury Risk Curves Using the Brier Metric Score

Many injury metrics are routinely proposed from measured or derived quantities from biomechanical experiments using post mortem human subjects (PMHS). The existing literature did not provide guidance on deciding between parameters collected in an experiment that would be best to use for the development of human injury probability curves (HIPC). The objective of this study was to use the Brier Metric Score (BMS) to identify the most appropriate metric from an experiment that predicts injury outcomes. The Brier Metric Score assesses how well a metric predicts the outcome for a censored data point (a lower BMS is better). Survival analysis was then conducted with the selected metric and the best distribution was selected using Akaike information criterion (AIC). Confidence intervals (CIs) and the normalized confidence interval width (NCIS) were calculated for the injury probability curve. The testing and validation of the methods described were performed using biomechanics data in the open literature. The methods for the HIPC development procedure detailed herein have been rigorously tested and used in the generation of WIAMan HIPCs and Injury Assessment Reference Curves (IARCs) for the WIAMan ATD, but can also be used in other ATD or PMHS injury risk curve development.

On the protective capacity of a safety vest for the thoracic injury caused by falling down

Background

Aged people all over the world are prone to fall down accidentally and be injured with fracture, such as the rib fracture. To protect the elderly, the safety vest has been developed to protect them from being injured when falling down. To effectively protect the elderly, more analysis on the protective capacity of a safety vest under different situation are needed.

Results

Herein, a finite element model based on the computed tomography CT scanning data of a Chinese old female was built, and then used to simulate the process of falling down at different velocities. Analysis and comparison were done on the maximum shear stress, kinetic energy curves and internal energy curves with and without safety vest. The maximum shear stress indicated that the Abbreviated Injury Scale (AIS) 2+ injury risks of rib were 8%, 100% and 100% at the velocities of 1.5 m/s, 2.0 m/s and 2.5 m/s, respectively. The corresponding risks were lowered to 0%, 0% and 60% by the vest, respectively. Furthermore, the vest could absorb the internal energy resulted by the deformation of the thoracic osseous tissue by about 20%, thus decreasing the shear stress and the injury risk.

Conclusion

It is concluded that the safety vest decreases the injury risk when the elderly fall down, thus protects them from being injured.

A Nonlinear Viscoelastic Model for Adipose Tissue Representing Tissue Response at a Wide Range of Strain Rates and High Strain Levels

Finite element human body models (FEHBMs) are nowadays commonly used to simulate pre- and in-crash occupant response in order to develop advanced safety systems. In this study, a biofidelic model for adipose tissue is developed for this application. It is a nonlinear viscoelastic model based on the Reese et al.'s formulation. The model is formulated in a large strain framework and applied for finite element (FE) simulation of two types of experiments: rheological experiments and ramped-displacement experiments. The adipose tissue behavior in both experiments is represented well by this model. It indicates the capability of the model to be used in large deformation and wide range of strain rates for application in human body models.