Comparing FE human body model rib geometry to population data

Isolated Rib Response and Fracture Prediction for Young Mid-Size Male, Enabled by Population Specific Material Models and Rib Cross-Sectional Geometry

Thorax injury remains a primary contributor to mortality in car crash scenarios. Although human body models can be used to investigate thorax response to impact, isolated rib models have not been able to predict age- and sex-specific force-displacement response and fracture location simultaneously, which is a critical step towards developing human thorax models able to accurately predict injury response. Recent advancements in constitutive models and quantification of age- and sex-specific material properties, cross-sectional area, and cortical bone thickness distribution offer opportunities to improve rib computational models. In the present study, improved cortical and trabecular bone constitutive models populated with age-specific material properties, age- and sex-specific population data on rib cross-sectional area, and cortical bone thickness distribution were implemented into an isolated 6th rib from a contemporary human body model. The enhanced rib model was simulated in anterior-posterior loading for comparison to experimental age- and sex-specific (twenty-three mid-size males, age range of 22- to 57-year-old) population force-displacement response and fracture location. The improved constitutive models, populated with age-specific material properties, proved critical to predict the rib failure force and displacement, while the improved cortical bone thickness distribution and cross-sectional area improved the fracture location prediction. The enhanced young mid-size male 6th rib model was able to predict young mid-size male 6th rib experimental force-displacement response and fracture location (overpredicted the displacement at failure by 35% and underpredicted the force at failure by 8% but within ±1 SD). The results of the present study can be integrated into full body models to potentially improve thorax injury prediction capabilities.

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.

Development and preliminary validation of computationally efficient and detailed 50th percentile female human body models

OBJECTIVE

No vehicle testing standard (physical or computational) employs a mid-sized female human surrogate, despite discrepancies related to injury outcomes for female occupants amongst all vehicle users. We detail the design and preliminary validation of 50th percentile female (F50) computational human body models (HBMs) based on Global Human Body Models Consortium (GHBMC) models.

METHOD

Data for the target geometry was collected as part of the initial generation of GHBMC models. Imaging, surface data, and 15 anthropomorphic measures from a living female subject (60.8 kg and 1.61 m) served as the baseline for model development. Due to the role rib cage geometry plays in biomechanical loading, rib cage morphology from secondary retrospective data was leveraged to identify an average female rib cage based on gross anatomical features. A female rib cage was selected from an existing dataset closest to the mean depth, height, and width of the set, considering only those aged 20 - 50 years. The selected subject among this secondary set also exhibited a 7th rib angle and sternum angle within 5% of the mean measurements, and within the range of previously reported studies. The GHBMC 5th percentile, small female detailed (high biofidelity) and simplified (computationally efficient) models were morphed to match the F50 subject body surface, selected bones, and mean rib cage using established thin plate spline techniques. The models were validated vs. previously published literature studies with an emphasis on rib cage response. Model data was compared to 47 channels of experimental data across four biomechanical hub simulations, two sled test simulations (one of which included all female PMHS), and two robustness simulations to test stability. Model results were mass scaled to the average of the reported corridors. Objective evaluation was conducted using CORA. IRB approval was obtained for all prospective and retrospective data collected or used. The target rib cage was selected from retrospective image data used in prior studies (n = 339 chest CT scans).

RESULTS

The morphed HBMs closely matched the target geometry. The detailed and simplified models had masses and element counts of 61.2 kg and 61.8 kg, and 2.8 million and 0.3 million, respectively. The mass difference is due to a coarser mesh in the simplified model. The simplified model ran 23 times faster than the detailed model on the same hardware. Each model exhibited stability in robustness tests, and the average CORA scores were 0.80 and 0.72 in the detailed and simplified models, respectively. The models performed well in frontal impacts against PMHS corridors after mass scaling.

CONCLUSIONS

Numerous recent studies underscore poorer injury outcomes for female vehicle occupants compared to males. While such outcomes are multifactorial, the average female models introduced in this work offer a novel tool within a widely used family of HBMs to reduce the outcome gap in terms of injury for all drivers. HBMs can be deployed in safety studies or in future regulatory requirements faster and more economically than a resized or newly designed ATDs aimed at the same target population.

Cross-sectional properties of rib geometry from an adult population

Introduction: Human body models (HBMs) play a key role in improving modern vehicle safety systems to protect broad populations. However, their geometry is commonly derived from single individuals chosen to meet global anthropometric targets, thus their internal anatomy may not fully represent the HBM's target demographic. Past studies show sixth rib cross-sectional geometry differences between HBM ribs and population-derived ribs, and corrections to HBM ribs based on these data have improved HBM's abilities to predict rib fracture locations. Methods: We measure and report average and standard deviations (SDs) in rib cross-sectional geometric properties derived from live subject CT scans of 240 adults aged 18-90. Male and female results are given as functions of rib number and rib lengthwise position for ribs 2 through 11. Population means/SDs are reported for measures of rib total area, rib cortical bone area, and rib endosteal area, as well as inertial moment properties of these rib sections. These population corridors are compared between males and females, and against the baseline rib geometries defined in six current HBMs. Results: Total cross-sectional area results found average males ribs to be larger than those of females by between approximately 1-2 SDs depending on rib number and position, and larger in cortical bone cross-sectional area by between 0-1 SDs. Inertial moment ratios showed female ribs being between approximately 0-1 SDs more elongated than male ribs, dependent again on rib number and position. Rib cross-sectional areas from 5 of the 6 HBMs were found to be overly large along substantial portions of most ribs when compared to average population corridors. Similarly, rib aspect ratios in HBMs deviated from average population data by up to 3 SDs in regions towards sternal rib ends. Discussion: Overall, while most HBMs capture overall trends such as reductions in cross-section along shaft lengths, many also exhibit local variation that deviates from population trends. This study's results provide the first reference values for assessing the cross-sectional geometry of human ribs across a wide range of rib levels. Results also further provide clear guidelines to improve rib geometry definitions present in current HBMs in order to better represent their target demographic.

Factors affecting the numerical response and fracture location of the GHBMC M50 rib in dynamic anterior-posterior loading

Rib fractures are common traumatic injuries, with links to increased morbidity and mortality. Finite element ribs from human body models have struggled to predict the force-displacement response, force and displacement at fracture, and the fracture location for isolated rib tests. In the current study, the sensitivity of a human body model rib with updated anisotropic and asymmetric material models to changes in boundary conditions, material properties, and geometry was investigated systematically to quantify contributions to response. The updated material models using uncalibrated average material properties from literature improved the force-displacement response of the model, whereas the cross-sectional geometry was the only parameter to effect fracture location. The resulting uncalibrated model with improved material models and cross-sectional geometry closely predicted experimental average force-displacement response and fracture location.

Finite element analysis of different fixation methods of screws on absorbable plate for rib fractures

Multiple rib fractures caused by trauma are common injuries and the internal fixation methods of these injuries have been paid more and more attention by surgeons. Absorbable plates and screws are the effective way to treat rib fractures, but there are no reports on which type of screw fixation method is most effective. In this study, finite element analysis was used to study the effects of five different types of screw fixation methods on anterior rib, lateral rib and posterior rib. The finite element model of the ribs was reconstructed from CT images, and the internal pressure (40 kPa) and intercostal force (30 N) on the surfaces of the ribs were simulated accordingly. An intercostal force of 30 N was applied to the upper and lower surfaces of the ribs to simulate the effect of intercostal muscle force. The pressure of 40 kPa was applied to the inner surface of the ribs, and the normal direction was applied to the inner surface of the ribs. The positive direction was considered inspiratory pressure, and the negative direction was considered expiratory pressure. The results indicate the optimal type of screw fixation on the absorbable plate for rib fractures, and provide a basis and reference for clinical application.

Comparing FE human body model rib geometry to population data

Finite element human body models (HBMs) are used to assess injury risk in a variety of impact scenarios. The ribs are a key structural component within the chest, so their accuracy within HBMs is vitally important for modeling human biomechanics. We assessed the geometric correspondence between the ribs defined within five widely used HBMs and measures drawn from population-wide studies of rib geometry, focusing on (1) rib global shape, (2) rib cross-sectional size and shape, and (3) rib cortical bone thickness. A parametric global shape model fitted to all HBM ribs was compared to expected rib parameters calculated for each HBM's subject demographic using population reference data. The GHBMC M50 and THUMS M50 male HBMs showed 24% and 50% of their fitted rib shape parameters (6 parameters per each 12 ribs) falling outside 1SD from population expected values, respectively. For female models the GHBMC F05, THUMS F05, and VIVA F50 models had 21%, 26%, and 19% of their rib shape parameters falling outside 1SD, respectively. Cross-sectional areas and inertial moments obtained along the HBM ribs were compared to average ± 1SD corridors for male and female ribs drawn from reference population data. The GHBMC M50, THUMS M50, and VIVA F50 model ribs were all larger in overall cross-sectional area than their targeted average population values by 0.9SDs (average across the rib's full length), 1.7SDs, and 1.3SDs, respectfully. When considering cortical bone cross-sectional area, the THUMS and VIVA models-which each define a constant bone thickness value across the entire rib-overestimated bone content on average by 1.1SDs and 1.2SDs, respectively. HBMs have traditionally performed poorly when predicting rib fracture onset or fracture site, and in all HBMs in this study the rib regions with the most extreme cortical bone thickness and cross-sectional area discrepancies (compared to average reference data) corresponded to regions toward the sternal end of the ribs where rib fractures most frequently occur. Results from this study highlight geometrical components of current HBM ribs that differ from the rib geometry that would be expected from within those models' target demographics, and help researchers prioritize improvements to their biofidelity.