Detailed subject-specific FE rib modeling for fracture prediction

Probabilistic Finite Element Analysis of Human Rib Biomechanics: A Framework for Improved Generalizability

In dynamic impact events, thoracic injuries often involve rib fractures, which are closely related to injury severity. Previous studies have investigated the behavior of isolated ribs under impact loading conditions, but often neglected the variability in anatomical shape and tissue material properties. In this study, we used probabilistic finite element analysis and statistical shape modeling to investigate the effect of population-wide variability in rib cortical bone tissue mechanical properties and rib shape on the biomechanical response of the rib to impact loading. Using the probabilistic finite element analysis results, a response surface model was generated to rapidly investigate the biomechanical response of an isolated rib under dynamic anterior-posterior load given the variability in rib morphometry and tissue material properties. The response surface was used to generate pre-fracture force-displacement computational corridors for the overall population and a population sub-group of older mid-sized males. When compared to the experimental data, the computational mean response had a RMSE of 4.28N (peak force 94N) and 6.11N (peak force 116N) for the overall population and sub-group respectively, whereas the normalized area metric when comparing the experimental and computational corridors ranged from 3.32% to 22.65% for the population and 10.90% to 32.81% for the sub-group. Furthermore, probabilistic sensitivities were computed in which the contribution of uncertainty and variability of the parameters of interest was quantified. The study found that rib cortical bone elastic modulus, rib morphometry and cortical thickness are the random variables that produce the largest variability in the predicted force-displacement response. The proposed framework offers a novel approach for accounting biological variability in a representative population and has the potential to improve the generalizability of findings in biomechanical studies.

A deep learning-based pipeline for developing multi-rib shape generative model with populational percentiles or anthropometrics as predictors

Rib cross-sectional shapes (characterized by the outer contour and cortical bone thickness) affect the rib mechanical response under impact loading, thereby influence the rib injury pattern and risk. A statistical description of the rib shapes or their correlations to anthropometrics is a prerequisite to the development of numerical human body models representing target demographics. Variational autoencoders (VAE) as anatomical shape generators remain to be explored in terms of utilizing the latent vectors to control or interpret the representativeness of the generated results. In this paper, we propose a pipeline for developing a multi-rib cross-sectional shape generative model from CT images, which consists of the achievement of rib cross-sectional shape data from CT images using an anatomical indexing system and regular grids, and a unified framework to fit shape distributions and associate shapes to anthropometrics for different rib categories. Specifically, we collected CT images including 3193 ribs, surface regular grid is generated for each rib based on anatomical coordinates, the rib cross-sectional shapes are characterized by nodal coordinates and cortical bone thickness. The tensor structure of shape data based on regular grids enable the implementation of CNNs in the conditional variational autoencoder (CVAE). The CVAE is trained against an auxiliary classifier to decouple the low-dimensional representations of the inter- and intra- variations and fit each intra-variation by a Gaussian distribution simultaneously. Random tree regressors are further leveraged to associate each continuous intra-class space with the corresponding anthropometrics of the subjects, i.e., age, height and weight. As a result, with the rib class labels and the latent vectors sampled from Gaussian distributions or predicted from anthropometrics as the inputs, the decoder can generate valid rib cross-sectional shapes of given class labels (male/female, 2nd to 11th ribs) for arbitrary populational percentiles or specific age, height and weight, which paves the road for future biomedical and biomechanical studies considering the diversity of rib shapes across the population.

Explaining and predicting the increased thorax injury in aged females: age and subject-specific thorax geometry coupled with improved bone constitutive models and age-specific material properties evaluated in side impact conditions

Predicting and understanding thorax injury is fundamental for the assessment and development of safety systems to mitigate injury risk to the increasing and vulnerable aged population. While computational human models have contributed to the understanding of injury biomechanics, contemporary human body models have struggled to predict rib fractures and explain the increased incidence of injury in the aged population. The present study enhanced young and aged human body models (HBMs) by integrating a biofidelic cortical bone constitutive model and population-based bone material properties. The HBMs were evaluated using side impact sled tests assessed using chest compression and number of rib fractures. The increase in thoracic kyphosis and the associated change in rib angle with increasing age, led to increased rib torsional moment increasing the rib shear stress. Coupled with and improved cortical bone constitutive model and aged material properties, the higher resulting shear stress led to an increased number of rib fractures in the aged model. The importance of shear stress resulting from torsional load was further investigated using an isolated rib model. In contrast, HBM chest compression, a common thorax injury-associated metric, was insensitive to the aging factors studied. This study proposes an explanation for the increased incidence of thorax injury with increasing age reported in epidemiological data, and provides an enhanced understanding of human rib mechanics that will benefit assessment and design of future safety systems.

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.

Analysis for Predictors of Failure of Orthodontic Mini-implant Using Patient-Specific Finite Element Models

In this study, we analyzed the clinical factors and mechanical parameters for predicting orthodontic mini-implant (OMI) failure in the mandible, which has different properties from the maxilla. A patient-specific finite element analysis was applied to 32 OMIs (6 failures and 26 successes) implanted between the mandibular second premolars and first molars used for anchorage. The peak stress and strain parameters were calculated for each sample. A logistic regression of the failure (vs. success) of OMIs on the mechanical parameters in the models was conducted. In addition, the influence of clinical factors on the mechanical parameters considered to be related to OMI failure was examined by a regression analysis. The mechanical parameter which best predicts OMI failure in the mandible was found to be a minimum principal strain of between 0.5 to 1.0 mm from the OMI surface (R² = 0.8033). The results indicate the patient's bone density, distance between the OMIs and adjacent root, and vertical implantation angle of the OMIs are potential clinical predictors of OMI failure in the mandible.

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.

Revisable and high-strength wheel-spun alginate/graphene oxide based fibrous rods towards a flexible and biodegradable rib internal fixation system

The intramedullary splint insertion fixation system is the mainstream clinical solution to severe rib fractures. However, the titanium alloy scaffolds have limitations in biocompatibility, flexibility and complexity of surgeries. Here we present a revisable wheel-spun alginate (Alg)/graphene oxide (GO)-based fibrous rod as a potential alternative for a rib internal fixation system. The reversible fusion and fission capability was obtained by optimized Alg/GO blended spinning and GO coating post-treatment. The mechanical performance of the demonstrated rod samples matches the properties of the human rib. A self-designed cubic matrix was used to conduct in situ cell culture. In vitro evaluation not only confirms the cell viability and migration on the fibers' surface, but also investigated the degradation and fission performance of fibrous rods. With a simple, minimally invasive implantation and controlled swelling, Alg/GO fibrous rods are able to tightly fix the rib fracture wound while maintaining sufficient flexibility.

Rib cortical bone thickness variation in adults by age and sex

Rib fractures are a common and serious outcome of blunt thoracic trauma and their likelihood is greater in older individuals. Osteoporotic bone loss is a well-documented aging phenomenon with sex-specific characteristics, but within rib bones, neither baseline maps of regional thickness nor the rates of bone thinning with age have been quantified across whole ribs. This study presents such data from 4014 ribs of 240 adult subjects aged 20-90. A validated cortical bone mapping technique was applied to clinical computed tomography scans to obtain local rib cortical bone thickness measurements over the surfaces of ribs 2 through 11. Regression models to age and sex gave rates of cortex thinning in local zones and aggregated across whole ribs. The statistical parametric mapping provided these relationships regionally as a function of rib surface location. All models showed significant reductions in bone thickness with age (p < 0.01). Average whole-rib thinning occurred at between 0.011 to 0.032 mm/decade (males) and 0.035 to 0.043 mm/decade (females), with sex and age accounting for up to 37% of population variability (R² ). Rates of thinning differed regionally and by rib, with the highest bone loss of up to 0.074 mm/decade occurring in mid-rib cutaneous and superior regions of ribs 2-6. Rates were consistently higher in females than males (significantly so across whole ribs but not all local regions) and were more pronounced in cutaneous, superior, and inferior rib aspects (average 0.025 mm/decade difference in ribs 4-8) compared to pleural aspects which had the thickest cortices but saw only minor differences in thinning rates by sex (0.045 mm/decade for females and 0.040 mm/decade for males). Regional analysis showed male and female bone thickness differences that were not statistically significant at 20 years of age (p > 0.05 across practically all regions) but subsequent cortex thinning meant that substantial pleural and cutaneous regions were thinner (p < 0.05) in females than males by 55 years of age. The techniques and results from this study can be applied to assess rib bone content loss in clinical settings across wide populations. Additionally, average cortex thickness results can be mapped directly to finite element models of the thorax, and regression results are used to modify such models to represent the ribs of men and women across their full adult lifespan.

Methodology to geometrically age human body models to average and subject-specific anthropometrics, demonstrated using a small stature female model assessed in a side impact

The aged population has been associated with an increased risk of injury in car-crash, creating a critical need for improved assessment of safety systems. Finite element human body models (HBMs) have been proposed, but require representative geometry of the aged population and high mesh quality. A new hybrid Morphing-CAD methodology was applied to a 26-year-old (YO) 5^th percentile female model to create average 75YO and subject-specific 86YO HBMs. The method achieved accurate morphing targets while retaining high mesh quality. The three HBMs were integrated into a side sled impact test demonstrating similar kinematic response but differing rib fracture patterns.

Global and local characterization explains the different mechanisms of failure of the human ribs

Knowledge of the mechanics and mechanistic reasons inducing rib fracture is fundamental for forensic investigations and for the design of implants and cardiopulmonary resuscitation devices. A mechanical rationale to explain the different rib mechanisms of failure is still a challenge. The aim of this work was to experimentally characterize human ribs to test the hypothesis that a correlation exists between the ribs properties and the mechanism of failure. 89 ribs were tested in antero-posterior compression. The full-field strain distribution was measured through Digital Image Correlation. The fracture load ranged 7-132 N. Two main different mechanisms of failure were observed: brittle and buckling. The strain analysis showed that the direction of principal strains was either aligned with the ribs, or oblique, around 45°, with a rather uniform direction in the most strained area. The maximum principal strains were in the range between 1000 and 30000 microstrain and the minimum principal strain between -30000 and -800 microstrain. The ribs undergoing brittle fracture had significantly thicker cortical bone than those undergoing buckling. Also, larger tensile strains were observed in the specimens with brittle fracture than in the buckling ones. These findings support the focus of cortical thickness modelling which could help in sharpening computational models for the aforesaid purposes.

Assessment of Thorax Finite Element Model Response for Behind Armor Blunt Trauma Impact Loading Using an Epidemiological Database

Nonperforating ballistic impacts on thoracic armor can cause blunt injuries, known as behind-armor blunt trauma (BABT). To evaluate the potential for this injury, the back face deformation (BFD) imprinted into a clay backing is measured; however, the link between BFD and potential for injury is uncertain. Computational human body models (HBMs) have the potential to provide an improved understanding of BABT injury risk to inform armor design but require assessment with relevant loading scenarios. In this study, a methodology was developed to apply BABT loading to a computational thorax model, enhanced with refined finite element mesh and high-deformation rate mechanical properties. The model was assessed using an epidemiological BABT survivor database. BABT impact boundary conditions for 10 cases from the database were recreated using experimentally measured deformation for specific armor/projectile combinations, and applied to the thorax model using a novel prescribed displacement methodology. The computational thorax model demonstrated numerical stability under BABT impact conditions. The predicted number of rib fractures, the magnitude of pulmonary contusion, and injury rank, increased with armor BFD, back face velocity, and input energy to the thorax. In three of the 10 cases, the model overpredicted the number of rib fractures, attributed to impact location positional sensitivity and limited details from the database. The integration of an HBM with the BABT loading method predicted rib fractures and injury ranks that were in good agreement with available medical records, providing a potential tool for future armor evaluation and injury assessment.

Subject-specific rib finite element models with material data derived from coupon tests under bending loading

Rib fractures are common thoracic injuries in motor vehicle crashes. Several human finite element (FE) human models have been created to numerically assess thoracic injury risks. However, the accurate prediction of rib biomechanical response has shown to be challenging due to human variation and modeling approaches. The main objective of this study was to better understand the role of modeling approaches on the biomechanical response of human ribs in anterior-posterior bending. Since the development of subject specific rib models is a time-consuming process, the second objective of this study was to develop an accurate morphing approach to quickly generate high quality subject specific rib meshes. The exterior geometries and cortical-trabecular boundaries of five human 6th-level ribs were extracted from CT-images. One rib mesh was developed in a parametric fashion and the other four ribs were developed with an in-house morphing algorithm. The morphing algorithm automatically defined landmarks on both the periosteal and endosteal boundaries of the cortical layer, which were used to morph the template nodes to target geometries. Three different cortical bone material models were defined based on the stress-strain data obtained from subject-specific tensile coupon tests for each rib. Full rib anterior-posterior bending tests were simulated based on data recorded in testing. The results showed similar trends to test data with some sensitivity relative to the material modeling approach. Additionally, the FE models were substantially more resistant to failure, highlighting the need for better techniques to model rib fracture. Overall, the results of this work can be used to improve the biofidelity of human rib finite element models.

Morphology of direct and indirect rib fractures

Rib fractures are a common finding in legal medicine and information on the impact mechanism is relevant for trauma reconstruction. This study focuses on morphological characteristics of rib fractures resulting from direct or indirect force. Fresh human ribs (n = 312) were divided into two groups and broken through local force (direct) and bending (indirect) in anterolateral areas. The ribs were macerated, visually investigated and the results statistically analysed. The indirect fractures showed a significant larger lateral offset of the internal and external fracture ends while the fracture ends of the direct fractures were more often straight, in line. Also, the morphology of the inner and outer fracture edges was significantly related to fracture type. Direct fractures mostly had rough and jagged inner edges (tension side) and straight, smooth outer edges (compression side), whereas indirect fractures more often showed the characteristics vice versa. The results were more convincing in combination and in ribs from persons aged ≤ 75 years at death. In summary, the direct and indirect rib fractures showed significantly different characteristics regarding orientation and offset of the fracture lines and roughness of the inner and outer fracture edges, which can be helpful to distinguish the traumatizing impact mechanisms in forensic autopsy routine.

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.