Numerical investigations of rib fracture failure models in different dynamic loading conditions

Rib fracture is one of the most common thoracic injuries in vehicle traffic accidents that can result in fatalities associated with seriously injured internal organs. A failure model is critical when modelling rib fracture to predict such injuries. Different rib failure models have been proposed in prediction of thorax injuries. However, the biofidelity of the fracture failure models when varying the loading conditions and the effects of a rib fracture failure model on prediction of thoracic injuries have been studied only to a limited extent. Therefore, this study aimed to investigate the effects of three rib failure models on prediction of thoracic injuries using a previously validated finite element model of the human thorax. The performance and biofidelity of each rib failure model were first evaluated by modelling rib responses to different loading conditions in two experimental configurations: (1) the three-point bending on the specimen taken from rib and (2) the anterior–posterior dynamic loading to an entire bony part of the rib. Furthermore, the simulation of the rib failure behaviour in the frontal impact to an entire thorax was conducted at varying velocities and the effects of the failure models were analysed with respect to the severity of rib cage damages. Simulation results demonstrated that the responses of the thorax model are similar to the general trends of the rib fracture responses reported in the experimental literature. However, they also indicated that the accuracy of the rib fracture prediction using a given failure model varies for different loading conditions.     

Effect of intercostal muscle and costovertebral joint material properties on human ribcage stiffness and kinematics

Current finite element (FE) models of the human thorax are limited by the lack of local-level validation, especially in the ribcage. This study exercised an existing FE ribcage model for a 50th percentile male under quasi-static point loading and dynamic sternal loading. Both force-displacement and kinematic responses of the ribcage were compared against experimental data. The sensitivity of the model response to changes in the material properties of the costovertebral (CV) joints and intercostal muscles was assessed. The simulations found that adjustments to the CV joints tended to change the amount of rib rotation in the sagittal plane, while changes to the elastic modulus and thickness of the intercostal muscles tended to alter both the stiffness and the direction and magnitude of rib motions. This study can lend insight into the role that the material properties of these two thoracic structures play in the dynamics of the ribcage during a frontal loading condition.     

The rib cage reduces intervertebral disc pressures in cadaveric thoracic spines by sharing loading under applied dynamic moments

The effects of the rib cage on thoracic spine loading are not well studied, but the rib cage may provide stability or share loads with the spine. Intervertebral disc pressure provides insight into spinal loading, but such measurements are lacking in the thoracic spine. Thus, our objective was to examine thoracic intradiscal pressures under applied pure moments, and to determine the effect of the rib cage on these pressures. Human cadaveric thoracic spine specimens were positioned upright in a testing machine, and Dynamic pure moments (0 to ±5 N·m) with a compressive follower load of 400 N were applied in axial rotation, flexion - extension, and lateral bending. Disc pressures were measured at T4-T5 and T8-T9 using needle-mounted pressure transducers, first with the rib cage intact, and again after the rib cage was removed. Changes in pressure vs. moment slopes with rib cage removal were examined. Pressure generally increased with applied moments, and pressure-moment slope increased with rib cage removal at T4-T5 for axial rotation, extension, and lateral bending, and at T8-T9 for axial rotation. The results suggest the intact rib cage carried about 62% and 56% of axial rotation moments about T4-T5 and T8-T9, respectively, as well as 42% of extension moment and 36–43% of lateral bending moment about T4-T5 only. The rib cage likely plays a larger role in supporting moments than compressive loads, and may also play a larger role in the upper thorax than the lower thorax.     

Effects of the rib cage on thoracic spine flexibility.

Besides protecting the internal organs of the thorax, the rib cage is the site of numerous muscle attachments. It also decreases the overall flexibility of the thoracic spine. This study developed finite element (FE) models of the thoracic spine with and without the rib cage, and the effects of the rib cage on thoracic spine flexibility were determined. The numerical models were validated by comparing the maximum rotation of the models for several loading cases with experimental data in the literature. After adapting the material properties for the discs and ligaments, the calculated maximum rotations differed from the measured median values by less than 1 degrees without the rib cage and by less than 2.5 degrees with it. The rib cage decreased the mean flexibility of the thoracic spine by 23% to 47%, depending on the loading plane. Assuming the ribs to be rigid beams required a corresponding reduction of ligament stiffnesses in order to achieve the same agreement of the maximum rotations with the measured median values. Interconnecting the FE thoracic spine model plus rib cage with the existing detailed FE lumbar spine model improves the simulation of force directions of muscles attached to the rib cage or thoracolumbar spine. In addition, such a model is suitable for determining the effects of lumbar spine implants on spinal balance.     

Thoracic force-deflection studies in living and embalmed primates

In order to evaluate the validity of the use of the cadaver in making estimates of the living human response to cardio-thoracic impact a series of tests were performed under similar conditions comparing the force-deflection response of the anaesthetized and the embalmed Rhesus monkey. Twelve primates (100 tests) were subjected to controlled thoracic impact with a Remington Humane Stunner while seated in a sled held firmly in place. Tests on 8 of the monkeys were performed while they were anaesthetized and again 30 days after being embalmed. Four monkeys were tested only after embalming. X-rays were used to evaluate rib fractures.

Results showed clear significant differences between anaesthetized primates and those tested 30 days after being embalmed. Forces developed under the same test conditions were 132 per cent higher for the embalmed primates compared to anaesthetized monkeys under static test conditions at maximum displacement and 25–50 per cent higher under dynamic testing conditions for all but the initial part of the tests. The results of these studies indicate that the mechanical characteristics of the primate thorax are dramatically affected by death and/or subsequent embalming procedures.     

Robust human body model injury prediction in simulated side impact crashes

This study developed a parametric methodology to robustly predict occupant injuries sustained in real-world crashes using a finite element (FE) human body model (HBM). One hundred and twenty near-side impact motor vehicle crashes were simulated over a range of parameters using a Toyota RAV4 (bullet vehicle), Ford Taurus (struck vehicle) FE models and a validated human body model (HBM) Total HUman Model for Safety (THUMS). Three bullet vehicle crash parameters (speed, location and angle) and two occupant parameters (seat position and age) were varied using a Latin hypercube design of Experiments. Four injury metrics (head injury criterion, half deflection, thoracic trauma index and pelvic force) were used to calculate injury risk. Rib fracture prediction and lung strain metrics were also analysed. As hypothesized, bullet speed had the greatest effect on each injury measure. Injury risk was reduced when bullet location was further from the B-pillar or when the bullet angle was more oblique. Age had strong correlation to rib fractures frequency and lung strain severity. The injuries from a real-world crash were predicted using two different methods by (1) subsampling the injury predictors from the 12 best crush profile matching simulations and (2) using regression models. Both injury prediction methods successfully predicted the case occupant's low risk for pelvic injury, high risk for thoracic injury, rib fractures and high lung strains with tight confidence intervals. This parametric methodology was successfully used to explore crash parameter interactions and to robustly predict real-world injuries.     

Three-dimensional osseo-ligamentous model of the thorax representing initiation of scoliosis by asymmetric growth

A biomechanical model of the human thorax was constructed to investigate how asymmetric growth of the thorax might initiate spinal lateral curvature and axial rotation as seen in scoliosis deformities. Geometric data specifying nodal points of the model were taken from stereo-radiographs of an adolescent subject. An initially symmetrical geometry was created by 'mirroring' measurements of a hemi-thorax and spine. Published data provided cross-sectional measurements of the ribs, material properties of tissues and global flexibilities of the intervertebral motion segments. The ribs, sternum, intervertebral motion segments and intercostal ligaments were represented by elastic elements. Model deformations were calculated by the direct stiffness finite element method, with growth represented by an initial strain term in the constitutive law. Non-linear behavior was accommodated by running the model recursively, with updated node locations at each step. Both stress relaxation and stress modulation of growth in the component tissues were simulated. Thoracic growth of 20% with asymmetric growth of the ribs was simulated to give rib length asymmetries of 11%. similar to that observed in a previous study of patients with idiopathic scoliosis. This resulted in the model having a small thoracic scoliosis curvature convex toward the side of the longer ribs. Variations of the model which permitted free motion at the costo-vertebral joints or produced changes in the curvature of the posterior parts of the ribs resulted in axial rotation of the vertebrae similar to that observed clinically. The model supports the idea that growth asymmetry could initiate a small scoliosis during adolescence.     

Chest wall motion during spontaneous breathing and mechanical ventilation in dogs

We measured the volume change of the thoracic cavity (delta Vth) and the volumes displaced by the diaphragm (delta Vdi) and rib cage (delta Vrc) in six pentobarbital-anesthetized dogs lying supine. A high-speed X-ray scanner (dynamic spatial reconstructor) provided three-dimensional images of the thorax during spontaneous breathing and during mechanical ventilation with paralysis. Tidal volume (VT) was measured by integrating gas flow. Changes in thoracic liquid volume (delta Vliq, presumably caused by changes in thoracic blood volume) were calculated as delta Vth - VT. Absolute volume displaced by the rib cage was not significantly different during the two modes of ventilation. During spontaneous breathing, thoracic blood volume increased during inspiration; delta Vliq was 12.3 +/- 4.1% of delta Vth. During mechanical ventilation, delta Vliq was nearly zero. Configuration of the relaxed chest wall was similar during muscular relaxation induced by either pharmacological paralysis or hyperventilation. Expiratory muscle activity produced 50 +/- 11% of the delta Vth during spontaneous breathing. We conclude that at constant VT the volume displaced by the rib cage is remarkably similar during the transition from spontaneous breathing to mechanical ventilation, while both diaphragmatic volume displacement and changes in intrathoracic blood volume decrease by a similar amount.     

In vitro analysis of kinematics and elastostatics of the human rib cage during thoracic spinal movement for the validation of numerical models

Neither kinematic nor stiffness properties of the rib cage during thoracic spinal motion were investigated in previous studies, while being essential for the accurate validation of numerical models of the whole thorax. The aim of this in vitro study therefore was to quantify the kinematics and elastostatics of the human rib cage under defined boundary conditions. Eight fresh frozen human thoracic spine specimens (C7-L1, median age 55 years, ranging from 40 to 60 years) including entire rib cages were loaded quasi-statically in flexion/extension, lateral bending, and axial rotation using pure moments of 5 Nm. Relative motions of ribs, thoracic vertebrae, and sternal structures as well as strains on the ribs were measured using optical motion tracking of 150 reflective markers per specimen, while specimens were loaded displacement-controlled with a constant rate of 1°/s for 3.5 cycles. The third full cycle was used to determine relative angles and strains at full loading of the spine for all motion directions. Largest relative angles were found in the main loading directions with only small motions at the mid-thoracic levels. Highest strains of the intercostal spaces were detected in the anterior section of the lowest fourth of the rib cage, showing compressions and elongations of more than 10% in all spinal motion planes. Elastostatic rib deformation was generally less than 1%. Rib-sternum relative motions exhibited complex motion patterns, overall showing relative angles below 2°. The results indicate that rib cage structures are not macroscopically deformed during spinal motion, but exhibit characteristic reproducible kinematics patterns.     

Blunt Liver Injury with Intact Ribs under Impacts on the Abdomen: A Biomechanical Investigation

Abdominal trauma accounts for nearly 20% of all severe traffic injuries and can often result from intentional physical violence, from which blunt liver injury is regarded as the most common result and is associated with a high mortality rate. Liver injury may be caused by a direct impact with a certain velocity and energy on the abdomen, which may result in a lacerated liver by penetration of fractured ribs. However, liver ruptures without rib cage fractures were found in autopsies in a series of cases. All the victims sustained punches on the abdomen by fist. Many studies have been dedicated to determining the mechanism underlying hepatic injury following abdominal trauma, but most have been empirical. The actual process and biomechanism of liver injury induced by blunt impact on the abdomen, especially with intact ribs remained, are still inexhaustive. In order to investigate this, finite element methods and numerical simulation technology were used. A finite element human torso model was developed from high resolution CT data. The model consists of geometrically-detailed liver and rib cage models and simplified models of soft tissues, thoracic and abdominal organs. Then, the torso model was used in simulations in which the right hypochondrium was punched by a fist from the frontal, lateral, and rear directions, and in each direction with several impact velocities. Overall, the results showed that liver rupture was primarily caused by a direct strike of the ribs induced by blunt impact to the abdomen. Among three impact directions, a lateral impact was most likely to cause liver injury with a minimum punch speed of 5 m/s (the momentum was about 2.447 kg.m/s). Liver injuries could occur in isolation and were not accompanied by rib fractures due to different material characteristics and injury tolerance.     

Development and validation of subject-specific finite element models for blunt trauma study

This study developed and validated finite element (FE) models of swine and human thoraxes and abdomens that had subject-specific anatomies and could accurately and efficiently predict body responses to blunt impacts. Anatomies of the rib cage, torso walls, thoracic, and abdominal organs were reconstructed from X-ray computed tomography (CT) images and extracted into geometries to build FE meshes. The rib cage was modeled as an inhomogeneous beam structure with geometry and bone material parameters determined directly from CT images. Meshes of soft components were generated by mapping structured mesh templates representative of organ topologies onto the geometries. The swine models were developed from and validated by 30 animal tests in which blunt insults were applied to swine subjects and CT images, chest wall motions, lung pressures, and pathological data were acquired. A comparison of the FE calculations of animal responses and experimental measurements showed a good agreement. The errors in calculated response time traces were within 10% for most tests. Calculated peak responses showed strong correlations with the experimental values. The stress concentration inside the ribs, lungs, and livers produced by FE simulations also compared favorably to the injury locations. A human FE model was developed from CT images from the Visible Human project and was scaled to simulate historical frontal and side post mortem human subject (PMHS) impact tests. The calculated chest deformation also showed a good agreement with the measurements. The models developed in this study can be of great value for studying blunt thoracic and abdominal trauma and for designing injury prevention techniques, equipments, and devices.     

Abdominal muscle activation changes if the purpose is to control pelvis motion or thorax motion

The aim of this study was to compare trunk muscular recruitment and lumbar spine kinematics when motion was constrained to either the thorax or the pelvis. Nine healthy women performed four upright standing planar movements (rotations, anterior–posterior translations, medial–lateral translations, and horizontal circles) while constraining pelvis motion and moving the thorax or moving the pelvis while minimizing thorax motion, and four isometric trunk exercises (conventional curl-up, reverse curl-up, cross curl-up, and reverse cross curl-up). Surface EMG (upper and lower rectus abdominis, lateral and medial aspects of external oblique, internal oblique, and latissimus dorsi) and 3D lumbar displacements were recorded. Pelvis movements produced higher EMG amplitudes of the oblique abdominals than thorax motions in most trials, and larger lumbar displacements in the medial–lateral translations and horizontal circles. Conversely, thorax movements produced larger rotational lumbar displacement than pelvis motions during rotations and higher EMG amplitudes for latissimus dorsi during rotations and anterior–posterior translations and for lower rectus abdominis during the crossed curl-ups. Thus, different neuromuscular compartments appear when the objective changes from pelvis to thorax motion. This would suggest that both movement patterns should be considered when planning spine stabilization programs, to optimize exercises for the movement and muscle activations desired.     

Effect of assumed stiffness and mass density on the impact response of the human chest using a three-dimensional FE model of the human body

The mass density, Young's modulus (E), tangent modulus (Et), and yield stress (sigma y) of the human ribs, sternum, internal organs, and muscles play important roles when determining impact responses of the chest associated with pendulum impact. A series of parametric studies was conducted using a commercially available three-dimensional finite element (FE) model, Total HUman Model for Safety (THUMS) of the whole human body, to determine the effect of changing these material properties on the predicted impact force, chest deflection, and the number of rib fractures and fractured ribs. Results from this parametric study indicate that the initial chest apparent stiffness was mainly influenced by the stiffness and mass density of the superficial muscles covering the torso. The number of rib fractures and fractured ribs was primarily determined by the stiffness of the ribcage. Similarly, the stiffness of the ribcage and internal organs contributed to the maximum chest deflection in frontal impact, while the maximum chest deflection for lateral impact was mainly affected by the stiffness of the ribcage. Additionally, the total mass of the whole chest had a moderately effect on the number of rib fractures.     

Finite element model of the human neck during omni-directional impacts. Part II: relation between cervical curvature and risk of injury

A detailed 3D FE model of the human neck was used to assess a possible relationship between risk of injury and cervical spine curvature for various impacts. A FE model was previously developed, representing the head and neck of a 50th percentile human with a normal lordotic curvature. The model behaviour was omni-directionally validated for various impacts using published results. For the present study, the model was deformed in order to obtain a straight and a kyphotic curvature, and for each geometry, rear-end, frontal, lateral and oblique impact were simulated. Although results showed similar kinematic patterns, significant differences were found in the distribution and peak values of ligament elongations, forces and moments along the cervical spine for the three configurations. It was concluded that the variability observed on the curvature of the human cervical spine may have a significant influence both on the behaviour and on the risk of injury of the neck during impact.     

Finite element model of the human neck during omni-directional impacts. Part II: relation between cervical curvature and risk of injury

A detailed 3D FE model of the human neck was used to assess a possible relationship between risk of injury and cervical spine curvature for various impacts. A FE model was previously developed, representing the head and neck of a 50th percentile human with a normal lordotic curvature. The model behaviour was omni-directionally validated for various impacts using published results. For the present study, the model was deformed in order to obtain a straight and a kyphotic curvature, and for each geometry, rear-end, frontal, lateral and oblique impact were simulated. Although results showed similar kinematic patterns, significant differences were found in the distribution and peak values of ligament elongations, forces and moments along the cervical spine for the three configurations. It was concluded that the variability observed on the curvature of the human cervical spine may have a significant influence both on the behaviour and on the risk of injury of the neck during impact.     

Chest wall motion during epidural anesthesia in dogs

To determine the relative contribution of rib cage and abdominal muscles to expiratory muscle activity during quiet breathing, we used lumbar epidural anesthesia in six pentobarbital sodium-anesthetized dogs lying supine to paralyze the abdominal muscles while leaving rib cage muscle motor function substantially intact. A high-speed X-ray scanner (Dynamic Spatial Reconstructor) provided three-dimensional images of the thorax. The contribution of expiratory muscle activity to tidal breathing was assessed by a comparison of chest wall configuration during relaxed apnea with that at end expiration. We found that expiratory muscle activity was responsible for approximately half of the changes in thoracic volume during inspiration. Paralysis of the abdominal muscles had little effect on the pattern of breathing, including the contribution of expiratory muscle activity to tidal breathing, in most dogs. We conclude that, although there is consistent phasic expiratory electrical activity in both the rib cage and the abdominal muscles of pentobarbital-anesthetized dogs lying supine, the muscles of the rib cage are mechanically the most important expiratory muscles during quiet breathing.     

Virtual traumatology of pregnant women: The PRegnant car Occupant Model for Impact Simulations (PROMIS)

This study report documents the development of a finite element (FE) model for analyzing trauma in pregnant women involved in road accidents and help the design of a specific safety device. The model is representative of a 50th percentile pregnant woman at 26 weeks of pregnancy in sitting position. To achieve this, the HUMOS 2 model, which has been validated in a wide range of dynamic tests, was scaled to the morphology of a woman in the 50th percentile and coupled with a model of gravid uterus. During scaling, special attention was paid to the pelvic region which is known to differ considerably in morphological terms between men and women. The gravid uterus model includes a placenta, a fetus, uterosacral ligaments and the amniotic fluid by means of fluid structure interaction formulation. The uterus and the female model were coupled using an original method whereby the growth of an uterus was simulated to compress the abdominal organs in a realistic manner. The model was validated based on experimental tests described in the literature. Additional tests based on abdominal loadings with a seatbelt on Post Mortem Human Surrogates (PMHS) coupled to silicone uterus were also performed.     

A model for studies of mechanical interactions between the human spine and rib cage

A three-dimensional mathematical model useful for studies of the mechanics of the human skeletal thorax is described. To construct this model, rib cage elements are incorporated into a previously reported model of the thoracolumbar spine. The vertebrae and bony portions of the ribs and sternum are idealized as rigid bodies. The behavior of the discs, ligaments and costal cartilages are modelled by deformable elements. Appropriate geometric and stiffness property data are assigned to the elements of the model. In constructing the model, it was found that the mechanical response of the costo-vertebral joint is strongly influenced by articulation geometry. Although rigid bodies were used to model calcified portions of the ribs, the model predicted rib cage deformations in close agreement with those measured experimentally. These studies indicate that the rigid body motion of calcified portions of the rib makes a major contribution to the deformation of the rib cage in response to certain types of loadings. Quantitative results are also reported on the roles the rib cage plays in bending responses of the spine, the lateral stability of the spine, and the production and correction of several scoliotic deformities.     

Structural change of the thorax in chronic obstructive pulmonary disease.

This study examines structural changes of the thorax in hyperinflated subjects with chronic obstructive pulmonary disease (COPD). Age-matched normal subjects were used for comparison. Thoracic dimensions were determined using anteroposterior and lateral chest radiographs performed at total lung capacity, functional residual capacity, and residual volume. Rib cage dimensions (lateral diameter, rib angle, anteroposterior diameter) and diaphragm position were determined at each lung volume. There were no significant differences in rib cage dimension between the COPD and normal subjects for all lung volumes. In contrast, the diaphragm was significantly lower in the COPD subjects. The change of rib cage dimensions in the COPD subjects (for a similar volume change) was not different from that in normal subjects, whereas the change of diaphragm position in the COPD subjects (for a similar volume change) was reduced. In conclusion, the primary structural change of the thorax in COPD with chronic hyperinflation is confined to the diaphragm, with no appreciable structural change in the rib cage.     

Implementation and validation of thoracic side impact injury prediction metrics in a human body model

This study's purpose was to implement injury metrics into the Total Human Model for Safety (THUMS) mirroring the spinal accelerometers, rib accelerometers and chest band instrumentation from two lateral post-mortem human subject sled test configurations. In both sled configurations, THUMS contacted a flat rigid surface (either a wall or beam) at 6.7 m/s. Sled A maximum simulated wall forces for the thorax, abdomen and pelvis were 7.1, 5.0 and 10.0 kN versus 5.7 ± 0.8, 3.4 ± 1.2 and 6.2 ± 2.7 kN experimentally. Sled B maximum simulated beam forces for the torso and pelvis were 8.0 and 7.6 kN versus 8.5 ± 0.2 and 7.9 ± 2.5 kN experimentally. Quantitatively, force magnitude contributed more to variation between simulated and experimental forces than phase shift. Acceleration-based injury metrics were within one standard deviation of experimental means except for the lower spine in the rigid wall sled test. These validated metrics will be useful for quantifying occupant loading conditions and calculating injury risks in various loading configurations.