Human Body Finite Element Research Articles

Behind Amor Blunt Trauma Lung and Liver Strains from Indenter Loading via Finite Element Modeling

Abstract To determine behind armor blunt trauma injury criteria, experiments have been conducted to launch projectiles at live swine at velocities ranging up to 65 m/s using one type of indenter design. To ensure the generalizability of the developed injury criteria, tests are needed with other indenter designs. The objectives of this study were to evaluate the kinematics and injury parameters from two indenter designs using human body finite element modeling. The simulation matrix consisted of indenter designs of chord and cylindrical shape, and 150 and 230 grams (g) mass, with impacts to liver and lung regions at 30 and 60 m/s. Rib and lung strains from lung impacts, and rib strains and liver strain energy densities from liver impacts were used to evaluate the design variables of mass and shape. Both designs played a role in skeletal and organ injury parameters. Analysis revealed the increased susceptibility for skeletal and organ traumas with the high mass indenter at high velocity impacts. The cylindrical indenter may be protective for organ injuries due to the larger area of loading on the ribcage compared to the chord indenter. Results from the chord indenter may serve as a conservative estimate of injury criteria.

Translational medical bioengineering research of traumatic brain injury among Chinese and American pedestrians caused by vehicle collision based on human body finite element modeling

Based on the average human body size in China and the THUMS AM50 finite element model of the human body, the Kriging interpolation algorithm was used to model the Chinese 50th percentile human body, and the biological fidelity of the model was verified. We built three different types of passenger vehicle models, namely, sedan, sports utility vehicle (SUV), and multi-purpose vehicle (MPV), and used mechanical response analysis and finite element simulation to compare and analyze the dynamic differences and head injury differences between the Chinese 50th percentile human body and the THUMS AM50 model during passenger vehicle collisions. The results showed that there are obvious differences between the Chinese mannequin and THUMS in terms of collision time, collision position, invasion speed, and angle. When a sedan collided with the mannequins, the skull damage to the Chinese human body model was more severe, and when a sedan or SUV collided, the brain damage to the Chinese human body was more severe. The abovementioned results suggest that the existing C-NCAP pedestrian protection testing regulations may not provide the best protection for Chinese human bodies, and that the regulations need to be improved by combining collision damage mechanisms and the physical characteristics of Chinese pedestrians. This thorough investigation is positioned to shed light on the fundamental biomechanics and injury mechanisms at play. Furthermore, the amalgamation of clinically rooted translational and engineering research in the realm of traumatic brain injury has the potential to establish a solid foundation for discerning preventive methodologies. Ultimately, this endeavor holds the potential to introduce effective strategies aimed at preventing and safeguarding against traumatic brain injuries.

A comparative analysis of dimensionality reduction surrogate modeling techniques for full human body finite element impact simulations

Fast-running surrogate computational models (simpler computational models) have been successfully used to replace time-intensive finite element models. However, it is unclear how well they perform in accurately and efficiently replicating complex, full human body finite element models. Here we survey several surrogate modeling techniques and assess their accuracy in predicting full strain fields of tissues of interest during a highly dynamic behind armor blunt trauma impact to the liver. We found that coupling dimensionality reduction on the high-dimensional output space (principal component analysis or autoencoders) with machine learning techniques (Gaussian Process Regression or multi-output neural networks) provides a framework capable of accurately and efficiently replacing complex full human body models. It was found that these surrogate models can successfully predict the strain fields (<10% average strain error) of select tissues during a nonlinear impact event but careful consideration should be given to element parsing and modeling technique.

Analysis of lumbar spine injury with different back inclinations under whole-body vibration: A finite element study based on whole human body models

Whole-body vibration was found to be a cause of low back pain. Different back inclinations might change the forces on the lumbar spine, resulting in different responses to the vibration. The aim of this study was to investigate the effects of back inclinations on the lumbar spine from the perspective of analysing the internal loads and deformations on the intervertebral discs. Whole human body finite element models at 90°, 95°, 100°, 105°, 110° and 115° inclinations were used to provide a whole-body condition when predicting the behaviour of the lumbar spine. Von Mises stress on the annulus fibrosus, intradiscal pressure, and intervertebral disc height were extracted. The Risk Factors were calculated to evaluate the spinal injury risk under long-term vibration conditions. The results showed that the internal loads and deformations on the intervertebral disc decreased and then increased with the increase of the inclination, and the responses were lower at the 95° and 100° inclinations. The Risk Factors at different inclinations at the 3 Hz load were 0.61, 0.49, 0.53, 0.58, 0.58, and 0.69, respectively, at the 7 Hz load were 0.97, 0.66, 0.71, 0.86, 0.87, and 0.97, respectively, which also showed that injury risk was at a lower level at 95° and 100° inclinations. This study found a relationship between injury risk and back inclination. Occupational drivers are advised to choose a back inclination between 95° and 100° to reduce the possible adverse effects of whole-body vibration during the working condition.

Modular incorporation of deformable spine and 3D neck musculature into a simplified human body finite element model

Spinal injuries are a concern for automotive applications, requiring large parametric studies to understand spinal injury mechanisms under complex loading conditions. Finite element computational human body models (e.g. Global Human Body Models Consortium (GHBMC) models) can be used to identify spinal injury mechanisms. However, the existing GHBMC detailed models (with high computational time) or GHBMC simplified models (lacking vertebral fracture prediction capabilities) are not ideal for studying spinal injury mechanisms in large parametric studies. To overcome these limitations, a modular 50th percentile male simplified occupant model combining advantages of both the simplified and detailed models, M50-OS + DeformSpine, was developed by incorporating the deformable spine and 3D neck musculature from the detailed GHBMC model M50-O (v6.0) into the simplified GHBMC model M50-OS (v2.3). This new modular model was validated against post-mortem human subject test data in four rigid hub impactor tests and two frontal impact sled tests. The M50-OS + DeformSpine model showed good agreement with experimental test data with an average CORrelation and Analysis (CORA) score of 0.82 for the hub impact tests and 0.75 for the sled impact tests. CORA scores were statistically similar overall between the M50-OS + DeformSpine (0.79 ± 0.11), M50-OS (0.79 ± 0.11), and M50-O (0.82 ± 0.11) models (p > 0.05). This new model is computationally 6 times faster than the detailed M50-O model, with added spinal injury prediction capabilities over the simplified M50-OS model.

Effect of whole-body vibration at different frequencies on the lumbar spine: A finite element study based on a whole human body model.

Many previous studies have found that occupational drivers commonly suffered from low back pain, and low back pain and degeneration of the intervertebral disc might be associated with vibration conditions. However, the biomechanical mechanisms of whole-body vibration that caused pain and injury were not clear. In this study, a validated whole human body finite element model was used, and vibration loads at frequencies of 3, 5, 7 and 9 Hz were loaded to evaluate the frequency effects on the spine. The results showed that the responses of the spine were strong at the 5 Hz vibration load. Vibration loads would produce alternating stresses and bulges in the annulus fibrosus and change the direction of the pressure in the nucleus pulposus. The posterior region of the intervertebral disc showed greater stress fluctuations than the anterior region. The Risk Factors showed that long-term exposure to whole-body vibrations at 5 and 7 Hz might have greater adverse effects on the spine. The findings of this study confirmed that vibrations near the resonance frequency of the human body would cause more injuries to the spine than other frequencies. Alternating stress and bulge might cause fatigue and the degeneration of the intervertebral disc, which might be the mechanisms of spinal injury caused by whole-body vibration, and the posterior regions of the intervertebral disc were more susceptible to degeneration. Some appropriate measures should be taken to reduce the adverse effects of whole-body vibration on spinal health.

Biomechanical evaluation of Back injuries during typical snowboarding backward falls.

To prevent spinal and back injuries in snowboarding, back protector devices (BPDs) have been increasingly used. The biomechanical knowledge for the BPD design and evaluation remains to be explored in snowboarding accident conditions. This study aims to evaluate back-to-snow impact conditions and the associated back injury mechanisms in typical snowboarding backward falls. A previously validated snowboarder multi-body model was first used to evaluate the impact zones on the back and the corresponding impact velocities in a total of 324 snowboarding backward falls. The biomechanical responses during back-to-snow impacts were then evaluated by applying the back-to-snow impact velocity to a full human body finite element model to fall on the snow ground of three levels of stiffness (soft, hard, and icy snow). The mean values of back-to-snow normal and tangential impact velocities were 2.4m/s and 7.3m/s with maximum values up to 4.8m/s and 18.5m/s. The lower spine had the highest normal impact velocity during snowboarding backward falls. The thoracic spine was found more likely to exceed the limits of flexion-extension range of motions than the lumbar spine during back-to-snow impacts, indicating a higher injury risk. On the hard and icy snow, rib cage and vertebral fractures were predicted at the costal cartilage and the posterior elements of the vertebrae. Despite the possible back injuries, the back-to-snow impact force was always lower than the force thresholds of the current BPD testing standard. The current work provides additional biomechanical knowledge for the future design of back protections for snowboarders.

An investigation of elderly occupant injury risks based on anthropometric changes compared to young counterparts

Objective The objective of the study was to investigate the difference between elderly and young occupant injury risks using human body finite element modeling in frontal impacts. Methods Two elderly male occupant models (representative age 70–80 years) were developed using the Global Human Body Consortium (GHBMC) 50th percentile as the baseline model. In the first elderly model (EM-1), material property changes were incorporated, and in the second elderly model (EM-2), material and anthropometric changes were incorporated. Material properties were based on literature. The baseline model was morphed to elderly anthropometry for EM-2. The three models were simulated in a frontal crash vehicle environment at 56 km/h. Responses from the two elderly and baseline models were compared with cadaver experimental data in thoracic, abdominal, and frontal impacts. Correlation and analysis scores were used for correlation with experimental data. The probabilities of head, neck, and thoracic injuries were assessed. Results The elderly models showed a good correlation with experimental responses. The elderly EM-1 had higher risk of head and brain injuries compared to the elderly EM-2 and baseline GHBMC models. The elderly EM-2 demonstrated higher risk of neck, chest, and abdominal injuries than the elderly EM-1 and baseline models. Conclusions The study investigated injury risks of two elderly occupants and compared to a young occupant in frontal crashes. The change in the material properties alone (EM-1) suggested that elderly occupants may be vulnerable to a greater risk of head and thoracic injuries, whereas change in both anthropometric and material properties (EM-2) suggested that elderly occupants may be vulnerable to a greater risk of thoracic and neck injuries. The second elderly model results were in better agreement with field injury data from the literature; thus, both anthropometric and material properties should be considered when assessing the injury risks of elderly occupants. The elderly models developed in this study can be used to simulate different impact conditions and determine injury risks for this group of our population.

Development and Validation of an Active Muscle Simplified Finite Element Human Body Model in a Standing Posture.

Active muscles play an important role in postural stabilization, and muscle-induced joint stiffening can alter the kinematic response of the human body, particularly that of the lower extremities, under dynamic loading conditions. There are few full-body human body finite element models with active muscles in a standing posture. Thus, the objective of this study was to develop and validate the M50-PS+Active model, an average-male simplified human body model in a standing posture with active musculature. The M50-PS+Active model was developed by incorporating 116 skeletal muscles, as one-dimensional beam elements with a Hill-type material model and closed-loop Proportional Integral Derivative (PID) controller muscle activation strategy, into the Global Human Body Models Consortium (GHBMC) simplified pedestrian model M50-PS. The M50-PS+Active model was first validated in a gravity standing test, showing the effectiveness of the active muscles in maintaining a standing posture under gravitational loading. The knee kinematics of the model were compared against volunteer kinematics in unsuited and suited step-down tests from NASA's active response gravity offload system (ARGOS) laboratory. The M50-PS+Active model showed good biofidelity with volunteer kinematics with an overall CORA score of 0.80, as compared to 0.64 (fair) in the passive M50-PS model. The M50-PS+Active model will serve as a useful tool to study the biomechanics of the human body in vehicle-pedestrian accidents, public transportation braking, and space missions piloted in a standing posture.

Development and Validation of a Whole Human Body Finite Element Model with Detailed Lumbar Spine

Investigations showed that low back pain of occupational drivers might be closely related to the whole-body vibration. Restricted by ethical concerns, the finite element method had become a viable alternative to invasive human experiments. Many mechanical behaviors of the human spine inside of the human body were unclear; therefore, a human whole-body finite element model might be required to better understand the lumbar behavior under whole-body vibration. In this study, a human whole-body finite element model with a detailed lumbar spine segment was developed. Several validations were performed to ensure the correctness of this model. The results of anthropometry and geometry validation, static validation, and dynamic validation were presented in this study. The validation results showed that the whole human body model was reasonable and valid by comparing with published data. The model developed in this study could reflect the biomechanical response of the human lumbar spine under vibration and could be used in further vibration analysis and offer proposals for protecting human body under whole-body vibration environment.

Numerical Simulation of Tram Collision with Pedestrian

This paper is focused on the field of traffic passive safety. The main point of this research is to reduce the severity of the consequences of a collision between tramway and pedestrian. The collision scenario is based on statistical data of local research in Czech Republic. For the first assessment of influence of tram geometry, it is necessary to investigate in the numerical simulations with different shapes of tram face. For this purpose the Multibody simulations are carried. Then the detail finite element simulations are prepared. The explicit finite element code Virtual Performance solution (VPS) and human body model Virthuman where used for the simulations. The approval process and standards were prepared for trams in last years (e.g Technical Guide Tramway front-end design) and already today it is visible significant contrast with automotive industry. The main complications in design process are connected with other requirements for tram front-end construction. These requirements are from other passive and active safety field (requirement for collision tram-tram, collision tram-road vehicle, active safety eg. visibility) and from another fields (durability, design, aerodynamics). The new construction safe for pedestrian is different in more ways. Therefore it is difficult to propagate the changes in construction.

Subcutaneous adipose tissue thickness around the ASIS area for human body models in reclined positions

Objective Subcutaneous adipose tissue (SAT) thickness above the anterior superior iliac spine (ASIS) influences belt fit of a vehicle occupant. To improve finite element (FE) human body models and their application assessing future seating positions in cars, there is a need for more detailed data. Methods Anthropometric input data were used to statistically model a lower limit of the SAT thickness in the area around the ASIS (at the ASIS or in the groin) extracted from 102 postmortem computed tomography (pmCT) data sets (56 males and 46 females). Additionally, 2 pmCT scans of 1 male individual in both supine and sitting conditions were used to estimate change in SAT thickness by position. Results Distributions and locations of minimum values for SAT thickness were derived for males and females. Sex, age, and body mass index (BMI) remained in a linear regression model for the minimum SAT thickness in the ASIS area. Thirty-seven percent of the variance in the SAT distribution of the sample can be explained by these input variables. The individual with data in supine and sitting positions showed an SAT thickness value above the ASIS 6 times higher in the sitting position than in the supine position. Conclusions Individual factors influence SAT thickness around the ASIS in addition to BMI, sex, and age. The presented values need to be regarded as a lower limit of SAT thickness, because in 63% the minimum was found in the groin area and the measurements were performed in a supine position. The increase in SAT thickness in a sitting position compared to supine seen in the case example shows the need for further data acquisition to establish a transfer function interpolating between both positions. The SAT thickness minimum values in the ASIS area shown in this study can provide valuable input for soft tissue modeling in human body models with the aim to analyze seat belt fit and to computationally assess lap belt and occupant interaction sensitivity to SAT tissue thickness under load. This might be crucial in reclined sitting positions in automated driving.

Analysis of seat cushion comfort by employing a finite element buttock model as a supplement to pressure measurement

Finite element (FE) simulation provides a flexible way to evaluate seat comfort, and employing a biomechanical human body model can improve its advantage by simulating deformation behaviors of internal tissues. It is not extensively discussed whether the biomechanical responses from these internal tissues in a FE human body model have a consistent trend with commonly used body surface pressure under different loading environments. Thus, this study aims to answer this question and also explore a coupling experimental and simulation method of material selection for seat comfort. Using a validated buttock model, different biomechanical responses for ten different seat cushion foams were explored. The results showed no consistent trend between stress, strain, and peak pressure values regarding different foams. This findings supported that using FE human body models for body comfort evaluation is a useful supplement to the common experimental measurement of body surface pressure.

Applicability of Neck Injury Criteria Critical Intercepts for Human Body Finite Element Models

Applicability of Neck Injury Criteria Critical Intercepts for Human Body Finite Element Models

Pelvic Response of a Total Human Body Finite Element Model During Simulated Injurious Under Body Blast Impacts

Abstract Military operations in Iraq and Afghanistan have resulted in the increased exposure of military personnel to explosive threats. Combat-related pelvic fractures are a relatively new battlefield injury that poses a serious threat to military personnel. Injury prediction for these events continues to be a challenge due to the limited availability of blast-specific test studies and the use of established automotive-based injury criteria that do not directly translate to combat-related exposures. The objective of this study is to evaluate the pelvic response of the global human body models consortium (GHBMC) 50th percentile detailed male model (v4.3) in under body blast (UBB) loading scenarios. Nine simulations were conducted with mild or enhanced threat levels, and nominal or obtuse occupant positions. Cross-sectional force outputs from the superior pubic ramus (SPR), ilium, and sacroiliac (SI) regions were evaluated using previously developed injury risk curves (IRC). Additionally, maximum principal strain (MPS) data were extracted from the pelvic cortical bone elements. Results showed that shear force was the best predictor of fracture for the ischial and SI regions, while axial force was the best predictor for the SPR region. These outcomes were consistent with the load path of the simulated UBB events. The obtuse posture had higher peak force values for injurious and noninjurious outcomes for the SPR and SI region. The nominal posture had higher peak force values for noninjurious outcomes in the ischial region. These outcomes were supported by the MPS response present in these postures.

THE INFLUENCE OF IMPACT SPEED ON CHEST INJURY OUTCOME IN WHOLE BODY FRONTAL SLED IMPACTS

While the seatbelt restraint has significantly improved occupant safety, the protection efficiency still needs further enhance to reduce the consequence of the crash. Influence of seatbelt restraint loading on chest injury under 40 km/h has been tested and documented. However, a comprehensive profiling of the efficiency of restraint systems with various impact speeds has not yet been sufficiently reported. The purpose of this study is to analyse the effect of the seatbelt loadings on chest injuries at different impact speeds utilizing a high bio-fidelity human body Finite Element (FE) model. Based on the whole-body frontal sled test configuration, the current simulation is setup using a substitute of Post-Mortem Human Subjects (PMHS). Chest injury outcomes from simulations are analysed in terms of design variables, such as seatbelt position parameters and collision speed in a full factorial experimental design. These outcomes are specifically referred to strain-based injury probabilities and four-point chest deflections caused by the change of the parameters. The results indicate that impact speed does influence chest injury outcome. The ribcage injury risk for more than 3 fractured ribs will increase from around 40 to nearly 100% when the impact speed change from 20 to 40 km/h if the seatbelt positioned at the middle-sternum of this study. Great injuries to the chest are mainly caused by the change of inertia, which indicates that chest injuries are greatly affected by the impact speed. Furthermore, the rib fracture risk and chest deflection are nonlinearly correlated with the change of the seatbelt position parameters. The study approach can serve as a reference for seatbelt virtual design. Meanwhile, it also provides basis for the research of chest injury mechanism.

Restraint systems considering occupant diversity and pre-crash posture

Objective Use volunteer data and parametric finite element (FE) human body models to investigate how restraint systems can be designed to adapt to a diverse population and pre-crash posture changes induced by active safety features. Methods Four FE human models were generated by morphing the midsize male GHBMC simplified model into geometries representing a midsize male, midsize female, short obese female (BMI 40 kg/m2), and large obese male (BMI 40 kg/m2) based on statistical skeleton and body shape geometry models. Each human model was positioned in a generic vehicle driver environment using two occupant pre-crash postures based on volunteer test results including one resulting from 1-g abrupt braking events. Improved restraint designs were manually developed for each occupant model in a 56 km/h frontal crash condition by adding a knee airbag, adjusting the shoulder belt load limit, steering column force, and driver airbag properties (tethers, inflation, and vent size). The improved designs were then tested at both pre-crash postures. Injury risks for the head, neck, chest, and lower extremities were analyzed. Results Human size and shape dominated the occupant injury measures, while the pre-crash-braking induced posture had minimal effects. Some of the safety concerns observed for large occupants include head strike-through the airbag and a conflict between head and chest injuries, which were mitigated by a stiffer restraint system with properly-tuned driver airbag. Chest injuries were a prominent safety concern for female occupants, mitigated by a softer seatbelt and smaller airbag size near the chest. Obese occupants exhibited a higher likelihood of lower extremity injuries indicating a need for a knee airbag. A diverse set of improved restraint designs were effective in lowering injury risks, indicating that restraint adaptability is necessary for accounting for occupant diversity. Conclusions This study investigated the effects of occupant size and shape variability, posture, and restraint design on injury risk for high-speed frontal crashes. More forward initial postures due to active safety features may decrease head, neck, and lower extremity injury risk, but may also increase chest injury risk. Safety concerns observed for large occupants include head strike-through and a conflict between head and chest injuries. Obese occupants had higher knee-thigh-hip injury risk. New restraints that adapt to occupant size and body shape may improve crash safety for all occupants. Further investigation is needed to confirm and extend the findings of this study.

Inverse material characterisation of human aortic tissue for traumatic injury in motor vehicle crashes

Traumatic aortic rupture (TAR) resulting due to motor vehicle crashes (MVC) accounts for about 10%–20% fatalities. This article presents the outcome of our findings of the non-linear stress–strain response and strain rate-dependent behaviour of the aortic tissue under crash like impacts. This study uses both finite element (FE) modelling and experimental testing to enhance the understanding of injury mechanisms associated with TAR. Accurate material properties are essential for correct FE model predictions. Therefore, the objective of the current study was to experimentally characterise and identify a suitable constitutive model and model parameters for aortic tissue that can be incorporated into FE human body models for studying aortic rupture. Inverse characterisation and genetic algorithms (GA) were used to train the FE model to simulate real-life scenarios applying loading in multiple directions. A total of 32 uniaxial tensile tests were conducted on human aortic tissues along with the longitudinal and circumferential directions by loading at different strain rates ranging up to 200/s. The engineering stress–strain relationship obtained for human aorta in the longitudinal and circumferential directions via uniaxial tensile tests demonstrated an approximately bilinear behaviour and strain rate dependency in both directions. The higher stresses and modulus in the circumferential direction as compared to longitudinal direction demonstrates the anisotropic behaviour of the tissue. A constitutive model has been developed and implemented via user subroutine (UMAT) in LS-DYNA that accounts for the strain rate-dependent effects and bilinear behaviour observed in the aortic tissue. A FE model of the experimental set-up was developed, and the parameters of the model were estimated using a GA-based inverse mapping technique. The developed model enables investigation of the mechanical response of the aortic tissue under crashes and other high rate loading conditions. The work is based on the premise that a reliable constitutive model coupled with a FE model of the aorta shall help predict TARs.

KINEMATICS AND INJURY MECHANISM OF CYCLIST LOWER LIMB IN VEHICLE-TO-BICYCLE COLLISIONS

Accident data show that lower limb is one of the most frequently injured body parts for cyclists in vehicle collisions. However, studies of cyclist lower limb injuries and protection are still sparse. Therefore, the purpose of this study is to investigate the kinematics and injury mechanism of cyclist lower limb in vehicle-to-bicycle collisions considering different impact boundary conditions. To achieve this, the finite element (FE) modeling approach and an FE human body lower limb model with detailed muscles were employed, and impact boundary conditions with different vehicle front-end shapes and cycling postures were considered. Predictions of lower limb kinematics, knee ligament elongation and bending moment of upper and lower leg were used for analysis. The simulation results show that cycling posture has a significant influence on cyclist lower limb kinematics and injury risk, lateral bending toward the direction of vehicle or vehicle moving combining with lateral shearing is the main mechanism for cyclist knee ligament injuries, and injuries to long bones of cyclist leg in vehicle impacts could form lateral bending at both directions. The findings suggest that the influence of cycling posture and distinct difference in injury mechanism between cyclist and pedestrian should be considered in the assessment of vehicle safety design for cyclist lower limb protection.

Human lumbar spinal column injury criteria from vertical loading at the base: Applications to military environments

The objective of this study was to determine force-based lumbar spine injury criteria due to vertical impact using Post Mortem Human Surrogate (PMHS) experiments. Mounted personnel in military vehicles sustain loads from the pelvis in combat events such as underbody blast loadings. Forty-three post mortem human subject thoracolumbar spinal columns were obtained, screened for pre-existing trauma, bone mineral densities (BMDs) were determined, pre-test radiological images were taken, fixed at the ends in polymethylmethacrylate, load cells were attached to the ends of the fixation, positioned on custom vertical accelerator device based on a military-seating posture, and impacted at the base. Posttest images were obtained, and gross dissection was done to confirm injuries, classified into single and multilevel groups, groups A and B. Axial and resultant forces at the thoracolumbar (proximal) and lumbosacral (distal) joints were used as response variables to develop lumbar spine injury risk curves using parametric survival analysis. The Brier score metric was used to rank the variables. Age, BMD, column length, and vertebral body and intervertebral disc areas were used as covariates. The optimal metric describing the underlying response to injury was the distal resultant force for group A and proximal axial force for group B specimens. Force-BMD for group A and force-body area for group B were the best combinations. The IRCs with ±95% confidence intervals and quality of risk curves are given in the paper, and they serve as lumbar spine injury criteria. The present human cadaver Injury Risk Curves (IRCs) can be used to conduct matched pair tests to obtain dummy-based injury assessment risk curves/values to predict injury. The present IRCs can be used in human body finite element models. The relationship between covariates and primary forces presented in this study contribute to a better understanding of the role of demographic, geometric, and material factors to impact acceleration loading.