A parametric ribcage geometry model accounting for variations among the adult population

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 First Step Toward a Family of Morphed Human Body Models Enabling Prediction of Population Injury Outcomes

The injury risk in a vehicle crash can depend on occupant specific factors. Virtual crash testing using finite element human body models (HBMs) to represent occupant variability can enable the development of vehicles with improved safety for all occupants. In this study, it was investigated how many HBMs of different sizes that are needed to represent a population crash outcome through a metamodel. Rib fracture risk was used as an example occupant injury outcome. Morphed HBMs representing variability in sex, height, and weight within defined population ranges were used to calculate population variability in rib fracture risk in a frontal and a side crash. Two regression methods, regularized linear regression with second-order terms and Gaussian process regression (GPR), were used to metamodel rib fracture risk due to occupant variability. By studying metamodel predictive performance as a function of training data, it was found that constructing GPR metamodels using 25 individuals of each sex appears sufficient to model the population rib fracture risk outcome in a general crash scenario. Further, by utilizing the known outcomes in the two crashes, an optimization method selected individuals representative for population outcomes across both crash scenarios. The optimization results showed that 5-7 individuals of each sex were sufficient to create predictive GPR metamodels. The optimization method can be extended for more crashes and vehicles, which can be used to identify a family of HBMs that are generally representative of population injury outcomes in future work.

Assessment of a fully-parametric thoraco-lumbar spine model generator with articulated ribcage

The present paper describes a novel user-friendly fully-parametric thoraco-lumbar spine CAD model generator including the ribcage, based on 22 independent parameters (1 posterior vertebral body height per vertebra + 4 sagittal alignment parameters, namely pelvic incidence, sacral slope, L1-L5 lumbar lordosis, and T1-T12 thoracic kyphosis). Reliable third-order polynomial regression equations were implemented in Solidworks to analytically calculate 56 morphological dependent parameters and to automatically generate the spine CAD model based on primitive geometrical features. A standard spine CAD model, representing the case-study of an average healthy adult, was then created and positively assessed in terms of spinal anatomy, ribcage morphology, and sagittal profile. The immediate translation from CAD to FEM for relevant biomechanical analyses was successfully demonstrated, first, importing the CAD model into Abaqus, and then, iteratively calibrating the constitutive parameters of one lumbar and three thoracic FSUs, with particular interest on the hyperelastic material properties of the IVD, and the spinal and costo-vertebral ligaments. The credibility of the resulting lumbo-sacral and thoracic spine FEM with/without ribcage were assessed and validated throughout comparison with extensive in vitro and in vivo data both in terms of kinematics (range of motion) and dynamics (intradiscal pressure) either collected under pure bending moments and complex loading conditions (bending moments + axial compressive force).

A fast methodology for generating skeletal FEM with detailed human geometric features based on CPD and RBF algorithms

Due to the significant effects of the human anatomical characteristics on the injury mechanism of passenger in traffic accidents, it is necessary to develop human body FEM (Finite Element Model) with detailed anatomical characteristics. However, traditional development of a human body FEM is an extremely complicated process. In particular, the meshing of human body is a huge and time-consuming project. In this paper, a new fast methodology based on CPD (Coherent Point Drift) and RBF (Radial Basis Function) was proposed to achieve the rapid developing the FEM of human bone with detailed anatomical characteristics. In this methodology, the mesh morphing technology based the RBF was used to generate FEM mesh in the geometry extracted from the target CT (Computed Tomography) data. In order to further improve the accuracy and speed of mesh morphing, the target geometric feature points required in the mesh morphing process were realized via the rapid and automatic generation based on the point-cloud registration technology of the CPD algorithm. Finally, this new methodology was used to generate a 3-year-old ribcage FEM consisting of a total of 27,728 elements with mesh size 3-5 mm based on the THUMS (Total Human Model for Safety) adult model. In the entire process of generating this new ribcage model, it only took about 2.7 s. The average error between the new FEM and target geometries was only about 2.7 mm. This indicated that the new FEM well described the detailed anatomical characteristics of target geometry, thus importantly revealing that the mesh quality of the new FEM was basically similar to that of source FEM.

Approximating subject-specific brain injury models via scaling based on head-brain morphological relationships

Most human head/brain models represent a generic adult male head/brain. They may suffer in accuracy when investigating traumatic brain injury (TBI) on a subject-specific basis. Subject-specific models can be developed from neuroimages; however, neuroimages are not typically available in practice. In this study, we establish simple and elegant regression models between brain outer surface morphology and head dimensions measured from neuroimages along with age and sex information (N = 191; 141 males and 50 females with age ranging 14-25 years). The regression models are then used to approximate subject-specific brain models by scaling a generic counterpart, without using neuroimages. Model geometrical accuracy is assessed using adjusted [Formula: see text] and absolute percentage error (e.g., 0.720 and 3.09 ± 2.38%, respectively, for brain volume when incorporating tragion-to-top). For a subset of 11 subjects (from smallest to largest in brain volume), impact-induced brain strains are compared with those from "morphed models" derived from neuroimage-based mesh warping. We find that regional peak strains from the scaled subject-specific models are comparable to those of the morphed counterparts but could be considerably different from those of the generic model (e.g., linear regression slope of 1.01-1.03 for gray and white matter regions versus 1.16-1.19, or up to ~ 20% overestimation for the smallest brain studied). These results highlight the importance of incorporating brain morphological variations in impact simulation and demonstrate the feasibility of approximating subject-specific brain models without neuroimages using age, sex, and easily measurable head dimensions. The scaled models may improve subject specificity for future TBI investigations.

Influences of human thorax variability on population rib fracture risk prediction using human body models

Rib fractures remain a common injury for vehicle occupants in crashes. The risk of a human sustaining rib fractures from thorax loading is highly variable, potentially due to a variability in individual factors such as material properties and geometry of the ribs and ribcage. Human body models (HBMs) with a detailed ribcage can be used as occupant substitutes to aid in the prediction of rib injury risk at the tissue level in crash analysis. To improve this capability, model parametrization can be used to represent human variability in simulation studies. The aim of this study was to identify the variations in the physical properties of the human thorax that have the most influence on rib fracture risk for the population of vehicle occupants. A total of 15 different geometrical and material factors, sourced from published literature, were varied in a parametrized SAFER HBM. Parametric sensitivity analyses were conducted for two crash configurations, frontal and near-side impacts. The results show that variability in rib cortical bone thickness, rib cortical bone material properties, and rib cross-sectional width had the greatest influence on the risk for an occupant to sustain two or more fractured ribs in both impacts. Therefore, it is recommended that these three parameters be included in rib fracture risk analysis with HBMs for the population of vehicle occupants.

Hello, world! VIVA+: A human body model lineup to evaluate sex-differences in crash protection

Finite element Human Body Models are increasingly becoming vital tools for injury assessment and are expected to play an important role in virtual vehicle safety testing. With the aim of realizing models to study sex-differences seen in the injury- and fatality-risks from epidemiology, we developed models that represent an average female and an average male. The models were developed with an objective to allow tissue-based skeletal injury assessment, and thus non-skeletal organs and joints were defined with simplified characterizations to enhance computational efficiency and robustness. The model lineup comprises female and male representations of (seated) vehicle occupants and (standing) vulnerable road users, enabling the safety assessment of broader segments of the road user population. In addition, a new workflow utilized in the model development is presented. In this workflow, one model (the seated female) served as the base model while all the other models were generated as closely-linked derivative models, differing only in terms of node coordinates and mass distribution. This approach opens new possibilities to develop and maintain further models as part of the model lineup, representing different types of road users to reflect the ongoing transitions in mobility patterns (like bicyclists and e-scooter users). In this paper, we evaluate the kinetic and kinematic responses of the occupant and standing models to blunt impacts, mainly on the torso, in different directions (front, lateral, and back). The front and lateral impacts to the thorax showed responses comparable to the experiments, while the back impact varied with the location of impact (T1 and T8). Abdomen bar impact showed a stiffer load-deflection response at higher intrusions beyond 40 mm, because of simplified representation of internal organs. The lateral shoulder impact responses were also slightly stiffer, presumably from the simplified shoulder joint definition. This paper is the first in a series describing the development and validation of the new Human Body Model lineup, VIVA+. With the inclusion of an average-sized female model as a standard model in the lineup, we seek to foster an equitable injury evaluation in future virtual safety assessments.

A parametric head geometry model accounting for variation among adolescent and young adult populations

BACKGROUND AND OBJECTIVE

Modeling the size and shape of human skull and scalp is essential for head injury assessment, design of helmets and head-borne equipment, and many other safety applications. Finite element (FE) head models are important tools to assess injury risks and design personal protective equipment. However, current FE head models are mainly developed based on the midsize male, failing to account for the significant morphological variation that exists in the skull and brain. The objective of this study was to develop a statistical head geometry model that accounts for size and shape variations among the adolescent and young adult population.

METHODS

To represent subject-specific geometry using a homologous mesh, threshold-based segmentation of head CT scans of 101 subjects between 14 and 25 years of age was performed, followed by landmarking, mesh morphing, and projection. Skull and scalp statistical geometry models were then developed as functions of age, sex, stature, BMI, head length, head breadth, and tragion-to-top of head using generalized Procrustes analysis (GPA), principal component analysis (PCA) and multivariate regression analysis.

RESULTS

The statistical geometry models account for a high percentage of morphological variations in scalp geometry (R²=0.63), outer skull geometry (R²=0.66), inner skull geometry (R²=0.55), and skull thickness (error < 1 mm) CONCLUSIONS: Skull and scalp statistical geometry models accounts for size and shape variations among the adolescent and young adult population were developed as functions of subject covariates. These models may serve as the geometric basis to develop individualized head FE models for injury assessment and design of head-borne equipment.

Biomechanical response of human rib cage to cardiopulmonary resuscitation maneuvers: Effects of the compression location

The biomechanical response of a human rib cage to cardiopulmonary resuscitation maneuvers was investigated by means of finite element simulations. We analyzed the effect of the location where the force was applied on the achieved compression depths and stress levels experienced by the breastbone and ribs. For compression locations on the breastbone, a caudal shift of the application area toward the breastbone tip resulted in a 17% reduction of the force required to achieve a target 5 cm compression depth. We found that the use of compression regions located on the costal cartilages would involve higher risk of rib fractures.

Evaluation of a diverse population of morphed human body models for prediction of vehicle occupant crash kinematics

Morphing can be used to alter human body models (HBMs) to represent a diverse population of occupants in car crashes. The mid-sized male SAFER HBM v9 was parametrically morphed to match 22 Post Mortem Human Subjects, loaded in different configurations. Kinetics and kinematics were compared for the morphed and baseline HBMs. In frontal impacts, the morphed HBMs correlated closer with the kinematics of obese subjects, but lower to small females. In lateral impacts HBM responses were too stiff. This study outlines a necessary evaluation of all HBMs that should be morphed to represent the diverse population in vehicle safety evaluations.

A statistical lumbar spine geometry model accounting for variations by Age, Sex, Stature, and body mass index

The objective of this study was to develop a statistical lumbar spine geometry model accounting for morphological variations among the adult population. Five lumber vertebrae and lumber spine curvature were collected from CT scans of 82 adult subjects through CT segmentation, landmark identification, and template mesh mapping. Generalized Procrustes Analysis (GPA), Principal Component Analysis (PCA), and multivariate regression analysis were conducted to develop the statistical lumbar spine model. Two statistical models were established to predict the vertebrae geometry and lumbar curvature respectively. Using the statistical models, a lumbar spine finite element (FE) model could be rapidly generated with a given set of age, sex, stature, and body mass index (BMI). The results showed that the lumbar spine vertebral size was significantly affected by stature, sex and age, and the lumbar spine curvature was significantly affected by stature and age. This statistical lumbar spine model could serve as the geometric basis for quantifying effects of human characteristics on lumbar spine injury risks in impact conditions.

Predicting pelvis geometry using a morphometric model with overall anthropometric variables

Pelvic fractures have been identified as the second most common AIS2+ injury in motor vehicle crashes, with the highest early mortality rate compared to other orthopaedic injuries. Further, the risk is associated with occupant sex, age, stature and body mass index (BMI). In this study, clinical pelvic CT scans from 132 adults (75 females, 57 males) were extracted from a patient database. The population shape variance in pelvis bone geometry was studied by Sparse Principal Component Analysis (SPCA) and a morphometric model was developed by multivariate linear regression using overall anthropometric variables (sex, age, stature, BMI). In the analysis, SPCA identified 15 principal components (PCs) describing 83.6% of the shape variations. Eight of these were significantly captured (α < 0.05) by the morphometric model, which predicted 29% of the total variance in pelvis geometry. The overall anthropometric variables were significantly related to geometrical features primarily in the inferior-anterior regions while being unable to significantly capture local sacrum features, shape and position of ASIS and lateral tilt of the iliac wings. In conclusion, a new detailed morphometric model of the pelvis bone demonstrated that overall anthropometric variables account for only 29% of the variance in pelvis geometry. Furthermore, variations in the superior-anterior region of the pelvis, with which the lap belt is intended to interact, were not captured. Depending on the scenario, shape variations not captured by overall anthropometry could have important implications for injury prediction in traffic safety analysis.

Development and Application of Digital Human Models in the Field of Vehicle Collisions: A Review

In the human-vehicle-road system of collisions, the human is the most important factor, and digital human models (DHMs) are developed with the aim of preventing or at least reducing human injury. Because most of the relevant literature is focused mainly on collisions in traffic accidents (TAs), only some of the literature reviewed in this paper involves research results on other aspects of collisions. In this review, based on the background of DHMs and the application of DHMs regarding human injury biomechanics in collisions field, research results regarding the development of DHMs are described, the methods for verifying such models are introduced, and the application of the research results is discussed based on the aspect of human injury biomechanics. From the research literature, the development and validation of DHMs and their application in human injury biomechanics are summarized, and future research trends are proposed and discussed.

Morphology and growth of the pediatric lumbar vertebrae

BACKGROUND CONTEXT

The majority of existing literature describing pediatric lumbar vertebral morphology are limited to characterization of the vertebral bodies, pedicles, and spinal canal and no study has described the rates of growth for any lumbar vertebral structure. While it is known that growth of the lumbar vertebrae results in changes in vertebral shape, the dimension ratios used to quantify these shape changes do not represent the 3D morphology of the vertebral structures. Additionally, many of the previous evaluations of growth and shape are purely descriptive and do not investigate sexual dimorphism or variations across vertebral levels.

PURPOSE

This study aims to establish a database of pediatric lumbar vertebra dimension, growth, and shape data for subjects between and ages of 1 and 19 years.

STUDY DESIGN

A retrospective study of computed tomography (CT) data.

METHODS

Retrospective, abdominal, CT scans of 102 skeletally normal pediatric subjects (54 males, 48 females) between the ages of 1 and 19 years were digitally reconstructed and manually segmented. Thirty surface landmark points (LMPs), 30 vertebral measurements, the centroid size, centroid location, and the local orientation were collected for each lumbar vertebra along with the centroid size of the LMPs comprising each subject's full lumbar spine and their intervertebral disc (IVD) heights. Nonparametric statistics were used to compare dimension values across vertebral levels and between sexes. Linear models with age as the independent variable were used to characterize dimension growth for each sex and vertebral level. Age-dependent quadratic equations were fit to LMP distributions resulting from a generalized Procrustes analysis (GPA) of the vertebrae and fixed effects models were used to investigate differences in model coefficients across levels and between sexes.

RESULTS

Intervertebral level dimension differences were observed across all vertebral structures in both sexes while pedicle widths and IVDs heights were the only measurements found to be sexually dimorphic. Dimension growth rates generally varied across vertebral levels and the growth rates of males were typically larger than those of females. Differences between male and female vertebral shapes were also found for all lumbar vertebral structures.

CONCLUSIONS

To the authors' knowledge, this is the first study to report growth rates for the majority of pediatric lumbar vertebral structures and the first to describe the 3D age-dependent shapes of the pediatric lumbar spine and vertebrae. In addition to providing a quantitative database, the dimension, growth, and shape data reported here would have applications in medical device design, surgical planning, surgical training, and biomechanical modeling.

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.

Validating diverse human body models against side impact tests with post-mortem human subjects

This study aimed at evaluating the ability of morphed finite element (FE) human body models (HBMs) to reproduce the impact responses of post-mortem human subjects (PMHS) with various stature and shape. Ten side impact tests previously performed using seven PMHS under 3 m/s and 8 m/s impact velocities were selected for model evaluation. With weight, stature, sex, and age of PMHS, seven FE HBMs were developed by morphing the midsize male THUMS model into the target geometries predicted by the statistical skeleton and external body shape models. The model-predicted force histories, accelerations along the spine, and deflections in the chest and abdomen were compared to the test data. For comparison, simulations in all testing conditions were also conducted with the original midsize male THUMS, and the results from the THUMS simulations were scaled to the weight and stature from each PMHS. The CORrelation and Analysis (CORA) was used to evaluate the model accuracy, with CORA scores close to one indicating excellent agreement. Ten simulations using the morphed models exhibited 0.80 ± 0.01, 0.80 ± 0.01, 0.78 ± 0.02, and 0.78 ± 0.02 CORA scores for the impact forces to the thorax, abdomen, iliac-wings, and greater-trochanter, respectively; the corresponding CORA scores with the original THUMS were markedly lower at 0.60 ± 0.06, 0.69 ± 0.05, 0.71 ± 0.05, and 0.69 ± 0.04; while those for the scaled THUMS were 0.65 ± 0.05, 0.71 ± 0.05, 0.73 ± 0.05, and 0.72 ± 0.02, also lower than the morphed models. Across all simulations, the morphed HBMs demonstrated significantly higher accuracy than the THUMS with or without scaling. These results suggested the necessity of accounting for size and shape effects on predicting human responses in side impacts.

Defining the shape of the scapulothoracic gliding surface

PURPOSE

The aim of the study is to evaluate the difference in shape of the upper part and lower part of the Scapulothoracic Gliding Surface (STGS).

METHODS

3D-CT images of the thoracic cage of 50 patients were created in MIMICS ^®. Three anatomical landmarks (insertion m. serratus anterior on 5th rib; transverse process of 2th and 7th vertebra) were used as an anteroposterior cutting plane to define the STGS. The upper part of the STG was defined as rib 2-5 and the lower part as 5-8. Next, in MATLAB ^®, a script was used to create the sphere with best fit for upper and lower parts of STGS. The Root-Square-Mean Error (RSME) (mm) between two closest points on the fitted sphere and the STGS of both parts were calculated to determine the goodness-of-fit.

RESULTS

The RSME was found to be significantly lower for the area ribs 2-5 (mean 7.85 mm, SD 1.86) compared the area of ribs 5-8 (mean 10.08 mm, SD 1.90).

CONCLUSION

The STGS of the upper thoracic wall (2-5) is more spherical shaped than the STGS of the lower thoracic wall (rib 5-8).

Development of a simplified, reproducible, parametric 3D model of the talus

Computational foot models have significant application in surgical decision making, injury and disease diagnosis and prevention, sports performance analysis and footwear engineering. However, due to the substantial time in model building and the heavy computational costs from the complexity of the models, daily clinical application of these foot models has yet to be achieved. Much of the previous research adopted a detailed-geometry approach in modeling bones that potentially contributed to the heavy computational costs. In this research, we developed a computational talus model based on CT section image data, image reconstruction and segmentation, contact surface identification, standard shape fitting, and finite element auto meshing algorithms. Modeling the bones as rigid is common, and modeling the contact surfaces only for the rigid body saves additional computational resources. Priority, therefore, in the shape fitting with optimization is given to the contact surfaces of the talus. Thirteen sets (9 males and 4 females) of CT section data were obtained. Image reconstruction, segmentation and bone labeling were conducted on each set of CT data to identify talus and its adjacent bones. Contact surfaces of the talus were then identified based on bone spatial relationships. Apart from the talar dome surface which was fitted by a 3rd-order polynomial, standard shapes such as ellipsoids and planes were used to fit the selected contact surfaces so that the geometrical parameters maintain physical significance. Based on these parameters, we automatically recreated and meshed the least-squares fitted shapes rapidly with limited elements. Last, mean major contact surfaces of the talus were obtained and fitted by standard shapes. Although the number of samples in this study was relatively small, our method provides sufficient and accurate geometric parameters of these contact surfaces to completely describe and reproduce the talus, on both a subject specific and average basis. The method for describing the talus here helps to parametrize computational models using planes and ellipsoids, improves surgical decision making and implants with a more precise and physically significant measures, and the description provides bone geometric parameters which can later be used to relate risk analysis for bone shape specific injury rates.

Implementation of controlling strategy in a biomechanical lower limb model with active muscles for coupling multibody dynamics and finite element analysis

Computational biomechanics for human body modeling has generally been categorized into two separated domains: finite element analysis and multibody dynamics. Combining the advantages of both domains is necessary when tissue stress and physical body motion are both of interest. However, the method for this topic is still in exploration. The aim of this study is to implement unique controlling strategies in finite element model for simultaneously simulating musculoskeletal body dynamics and in vivo stress inside human tissues. A finite element lower limb model with 3D active muscles was selected for the implementation of controlling strategies, which was further validated against in-vivo human motion experiments. A unique feedback control strategy that couples together a basic Proportion-Integration-Differentiation (PID) controller and generic active signals from Computed Muscle Control (CMC) method of the musculoskeletal model or normalized EMG singles was proposed and applied in the present model. The results show that the new proposed controlling strategy show a good correlation with experimental test data of the normal gait considering joint kinematics, while stress distribution of local lower limb tissue can be also detected in real-time with lower limb motion. In summary, the present work is the first step for the application of active controlling strategy in the finite element model for concurrent simulation of both body dynamics and tissue stress. In the future, the present method can be further developed to apply it in various fields for human biomechanical analysis to monitor local stress and strain distribution by simultaneously simulating human locomotion.

Regional maps of rib cortical bone thickness and cross-sectional geometry

Here we present detailed regional bone thickness and cross-sectional measurements from full adult ribs using high resolution CT scans processed with a cortical bone mapping technique. Sixth ribs from 33 subjects ranging from 24 to 99 years of age were used to produce average cortical bone thickness maps and to provide average ± 1SD corridors for expected cross-section properties (cross-sectional areas and inertial moments) as a function of rib length. Results obtained from CT data were validated at specific rib locations using direct measurements from cut sections. Individual thickness measurements from CT had an accuracy (mean error) and precision (SD error) of -0.013 ± 0.167 mm (R² coefficient of determination of 0.84). CT-based measurement errors for rib cross-sectional geometry were -0.1 ± 13.1% (cortical bone cross-sectional area) and 4.7 ± 1.8% (total cross-sectional area). Rib cortical bone thickness maps show the expected regional variation across a typical rib's surface. The local mid-rib maxima in cortical thickness along the pleural rib aspect ranged from range 0.9 to 2.6 mm across the study population with an average map maximum of 1.4 mm. Along the cutaneous aspect, rib cortical bone thickness ranged from 0.7 to 1.9 mm with an average map thickness of 0.9 mm. Average cross-sectional properties show a steady reduction in total cortical bone area from 10% along the rib's length through to the sternal end, whereas overall cross-sectional area remains relatively constant along the majority of the rib's length before rising steeply towards the sternal end. On average, male ribs contained more cortical bone within a given cross-section than was seen for female ribs. Importantly, however, this difference was driven by male ribs having larger overall cross-sectional areas, rather than by sex differences in the bone thickness observed at specific local cortex sites. The cortical bone thickness results here can be used directly to improve the accuracy of current human body and rib models. Furthermore, the measurement corridors obtained from adult subjects across a wide age range can be used to validate future measurements from more widely available image sources such as clinical CT where gold standard reference measures (e.g. such as direct measurements obtained from cut sections) are otherwise unobtainable.

Development and validation of a semi-automatic landmark extraction method for mesh morphing

Currently, landmark-based mesh morphing technology is widely used to rapidly obtain meshes with specific geometry, which is suitable to develop parametric human finite element (FE) models. However it takes too much time for landmark extraction to obtain high geometric accuracy. The purpose of this study is to develop and validate a semi-automatic landmark extraction method to reduce the time of manual selection of landmarks without sacrificing the accuracy of identifying landmarks in the process of mesh morphing. A few contour edge landmarks were extracted manually. Mathematical landmarks and pseudo-landmarks were extracted automatically by user-defined algorithm. The radial basis function (RBF) was used to morph the baseline FE model into the target geometry based on these landmarks. The cervical vertebra (C5), rib (R7) and femur were selected as the target geometries to verify the effectiveness of the method. The maximum mean geometric error of the three types of target geometries was less than 1 mm. The mesh quality of the morphed FE model was similar to that of the baseline FE model. Compared to the traditional manual method, 2/3 to 3/4 of the time for landmark extraction was saved by the semi-automatic method.

Frontal crash simulations using parametric human models representing a diverse population

Objective: Analyses of crash data have shown that older, obese, and/or female occupants have a higher risk of injury in frontal crashes compared to the rest of the population. The objective of this study was to use parametric finite element (FE) human models to assess the increased injury risks and identify safety concerns for these vulnerable populations. Methods: We sampled 100 occupants based on age, sex, stature, and body mass index (BMI) to span a wide range of the U.S. adult population. The target anatomical geometry for each of the 100 models was predicted by the statistical geometry models for the rib cage, pelvis, femur, tibia, and external body surface developed previously. A regional landmark-based mesh morphing method was used to morph the Global Human Body Models Consortium (GHBMC) M50-OS model into the target geometries. The morphed human models were then positioned in a validated generic vehicle driver compartment model using a statistical driving posture model. Frontal crash simulations based on U.S. New Car Assessment Program (U.S. NCAP) were conducted. Body region injury risks were calculated based on the risk curves used in the US NCAP, except that scaling was used for the neck, chest, and knee-thigh-hip injury risk curves based on the sizes of the bony structures in the corresponding body regions. Age effects were also considered for predicting chest injury risk. Results: The simulations demonstrated that driver stature and body shape affect occupant interactions with the restraints and consequently affect occupant kinematics and injury risks in severe frontal crashes. U-shaped relations between occupant stature/weight and head injury risk were observed. Chest injury risk was strongly affected by age and sex, with older female occupants having the highest risk. A strong correlation was also observed between BMI and knee-thigh-hip injury risk, whereas none of the occupant parameters meaningfully affected neck injury risks. Conclusions: This study is the first to use a large set of diverse FE human models to investigate the combined effects of age, sex, stature, and BMI on injury risks in frontal crashes. The study demonstrated that parametric human models can effectively predict the injury trends for the population and may now be used to optimize restraint systems for people who are not similar in size and shape to the available anthropomorphic test devices (ATDs). New restraints that adapt to occupant age, sex, stature, and body shape may improve crash safety for all occupants.

Impact Response Comparison Between Parametric Human Models and Postmortem Human Subjects with a Wide Range of Obesity Levels

OBJECTIVE

Field data analyses have shown that obesity significantly increases the occupant injury risks in motor vehicle crashes, but the injury assessment tools for people with obesity are largely lacking. The objectives of this study were to use a mesh morphing method to rapidly generate parametric finite element models with a wide range of obesity levels and to evaluate their biofidelity against impact tests using postmortem human subjects (PMHS).

METHODS

Frontal crash tests using three PMHS seated in a vehicle rear seat compartment with body mass index (BMI) from 24 to 40 kg/m² were selected. To develop the human models matching the PMHS geometry, statistical models of external body shape, rib cage, pelvis, and femur were applied to predict the target geometry using age, sex, stature, and BMI. A mesh morphing method based on radial basis functions was used to rapidly morph a baseline human model into the target geometry. The model-predicted body excursions and injury measures were compared to the PMHS tests.

RESULTS

Comparisons of occupant kinematics and injury measures between the tests and simulations showed reasonable correlations across the wide range of BMI levels.

CONCLUSIONS

The parametric human models have the capability to account for the obesity effects on the occupant impact responses and injury risks.

The effect of age and demographics on rib shape

Elderly populations have a higher risk of rib fractures and other associated thoracic injuries than younger adults, and the changes in body morphology that occur with age are a potential cause of this increased risk. Rib centroidal path geometry for 20 627 ribs was extracted from computed tomography (CT) scans of 1042 live adult subjects, then fitted to a six-parameter mathematical model that accurately characterizes rib size and shape, and a three-parameter model of rib orientation within the body. Multivariable regression characterized the independent effect of age, height, weight, and sex on the rib shape and orientation across the adult population, and statistically significant effects were seen from all demographic factors (P < 0.0001). This study reports a novel aging effect whereby both the rib end-to-end separation and rib aspect ratio are seen to increase with age, producing elongated and flatter overall rib shapes in elderly populations, with age alone explaining up to 20% of population variability in the aspect ratio of mid-level ribs. Age was not strongly associated with overall rib arc length, indicating that age effects were related to shape change rather than overall bone length. The rib shape effect was found to be more strongly and directly associated with age than previously documented age-related changes in rib angulation. Other demographic results showed height and sex being most strongly associated with rib size, and weight most strongly associated with rib pump-handle angle. Results from the study provide a statistical model for building rib shapes typical of any given demographic by age, height, weight, and sex, and can be used to help build population-specific computational models of the thoracic rib cage. Furthermore, results also quantify normal population ranges for rib shape parameters which can be used to improve the assessment and treatment of rib skeletal deformity and disease.

Age-related changes in thoracic skeletal geometry of elderly females

OBJECTIVE

Both females and the elderly have been identified as vulnerable populations with increased injury and mortality risk in multiple crash scenarios. Particularly in frontal impacts, older females show higher risk to the chest and thorax than their younger or male counterparts. Thoracic geometry plays a role in this increase, and this study aims to quantify key parts of that geometry in a way that can directly inform human body models that incorporate the concept of person age.

METHODS

Computed tomography scans from 2 female subject groups aged 20-35 and 65-99 were selected from the International Center for Automotive Medicine scan database representing young and old female populations. A model of thoracic skeletal anatomy was built for each subject from independent parametric models of the spine, ribs, and sternum, along with further parametric models of those components' spatial relationships. Parameter values between the 2 groups are directly compared, and average parameter values within each group are used to generate statistically average skeletal geometry for young and old females. In addition to the anatomic measures explicitly used in the parameterization scheme, key measures of rib cage depth and spine curvature are taken from both the underlying subject pool and from the resultant representative geometries.

RESULTS

Statistically significant differences were seen between the young and old groups' spine and rib anatomic components, with no significant differences in local sternal geometry found. Vertebral segments in older females had higher angles relative to their inferior neighbors, providing a quantification of the kyphotic curvature known to be associated with age. Ribs in older females had greater end-to-end span, greater aspect ratio, and reduced out-of-plane deviation, producing an elongated and overall flatter curvature that leads to distal rib ends extending further anteriorly in older individuals. Combined differences in spine curvature and rib geometry led to an 18-mm difference in anterior placement of the sternum between young and old subjects.

CONCLUSIONS

This study provides new geometric data regarding the variability in anthropometry of adult females with age and has utility in advancing the veracity of current human body models. A simplified scaffold representation of underlying 3-dimensional bones within the thorax is presented, and the reported young and old female parameter sets can be used to characterize the anatomic differences expected with age and to both validate and drive morphing algorithms for aged human body models. The modular approach taken allows model parameters to hold inherent and intuitive meaning, offering advantages over more generalized methods such as principal component analysis. Geometry can be assessed on a component level or a whole thorax level, and the parametric representation of thorax shape allows direct comparisons between the current study and other individuals or human body models.

What are the differences in injury patterns of young and elderly traffic accident fatalities considering death on scene and death in hospital?

Older traffic participants have higher risks of injury than the population up to 65 years in case of comparable road traffic accidents and further, higher mortality rates at comparable injury severities. Rib fractures as risk factors are currently discussed. However, death on scene is associated with hardly survivable injuries and might not be a matter of neither rib fractures nor age. As 60% of traffic accident fatalities are estimated to die on scene, they are not captured in hospital-based trauma registries and injury patterns remain unknown. Our database comprises 309 road traffic fatalities, autopsied at the Institute of Legal Medicine Munich in 2004 and 2005. Injuries are coded according to Abbreviated Injury Scale, AIS© 2005 update 2008 [1]. Data used for this analysis are age, sex, site of death, site of accident, traffic participation mode, measures of injury severity, and rib fractures. The injury patterns of elderly, aged 65+ years, are compared to the younger ones divided by their site of death. Elderly with death on scene more often show serious thorax injuries and pelvic fractures than the younger. Some hints point towards older fatalities showing less frequently serious abdominal injuries. In hospital, elderly fatalities show lower Injury Severity Scores (ISSs) compared to the younger. The number of rib fractures is significantly higher for the elderly but is not the reason for death. Results show that young and old fatalities have different injury patterns and reveal first hints towards the need to analyze death on scene more in-depth.

The Effect of Rib Shape on Stiffness

This study investigates the isolated effect of rib shape on the mechanical characteristics of ribs subjected to multiple forms of loading. It aims to measure the variation in stiffness due to shape that is seen throughout the population and, in particular, provide a tool for researchers to better understand the influence of shape on resulting stiffness. A previously published six-parameter shape model of the central axis of human ribs was used. It has been shown to accurately model the overall rib path using intrinsic geometric properties such as size, aspect ratio, and skewness, through shapes based on logarithmic spirals with high curvature continuity. In this study the model was fitted to 19,500 ribs from 989 adult female and male CT scans having demographic distributions matching the US adult population. Mechanical loading was simulated through a simplified finite element model aimed at isolating rib shape from other factors influencing mechanical response. Four loading scenarios were used representing idealized free and constrained loading conditions in axial (body-anterior) and lateral directions. Characteristic rib stiffness and maximum stress location were tracked as simulation output measures. Regression models of rib stiffness found that all shape model parameters added information when predicting stiffness under each loading condition, with their linear combination able to account for 95% of the population stiffness variation due to shape in midlevel ribs for free axial loading, and 92%-98% in other conditions. Full regression models including interactive terms explained up to 99% of population variability. Results allow researchers to better evaluate the differences in stiffness results that are obtained from physical testing by providing a framework with which to explain variation due to rib shape.

Development, Evaluation, and Sensitivity Analysis of Parametric Finite Element Whole-Body Human Models in Side Impacts

Occupant stature and body shape may have significant effects on injury risks in motor vehicle crashes, but the current finite element (FE) human body models (HBMs) only represent occupants with a few sizes and shapes. Our recent studies have demonstrated that, by using a mesh morphing method, parametric FE HBMs can be rapidly developed for representing a diverse population. However, the biofidelity of those models across a wide range of human attributes has not been established. Therefore, the objectives of this study are 1) to evaluate the accuracy of HBMs considering subject-specific geometry information, and 2) to apply the parametric HBMs in a sensitivity analysis for identifying the specific parameters affecting body responses in side impact conditions. Four side-impact tests with two male post-mortem human subjects (PMHSs) were selected to evaluate the accuracy of the geometry and impact responses of the morphed HBMs. For each PMHS test, three HBMs were simulated to compare with the test results: the original Total Human Model for Safety (THUMS) v4.01 (O-THUMS), a parametric THUMS (P-THUMS), and a subject-specific THUMS (S-THUMS). The P-THUMS geometry was predicted from only age, sex, stature, and BMI using our statistical geometry models of skeleton and body shape, while the S-THUMS geometry was based on each PMHS's CT data. The simulation results showed a preliminary trend that the correlations between the PTHUMS- predicted impact responses and the four PMHS tests (mean-CORA: 0.84, 0.78, 0.69, 0.70) were better than those between the O-THUMS and the normalized PMHS responses (mean-CORA: 0.74, 0.72, 0.55, 0.63), while they are similar to the correlations between S-THUMS and the PMHS tests (mean-CORA: 0.85, 0.85, 0.67, 0.72). The sensitivity analysis using the PTHUMS showed that, in side impact conditions, the HBM skeleton and body shape geometries as well as the body posture were more important in modeling the occupant impact responses than the bone and soft tissue material properties and the padding stiffness with the given parameter ranges. More investigations are needed to further support these findings.