Finite element modeling of potential cervical spine pain sources in neutral position low speed rear impact

Quantifying the Importance of Active Muscle Repositioning a Finite Element Neck Model in Flexion Using Kinematic, Kinetic, and Tissue-Level Responses

PURPOSE

Non-neutral neck positions are important initial conditions in impact scenarios, associated with a higher incidence of injury. Repositioning in finite element (FE) neck models is often achieved by applying external boundary conditions (BCs) to the head while constraining the first thoracic vertebra (T1). However, in vivo, neck muscles contract to achieve a desired head and neck position generating initial loads and deformations in the tissues. In the present study, a new muscle-based repositioning method was compared to traditional applied BCs using a contemporary FE neck model for forward head flexion of 30°.

METHODS

For the BC method, an external moment (2.6 Nm) was applied to the head with T1 fixed, while for the muscle-based method, the flexors and extensors were co-contracted under gravity loading to achieve the target flexion.

RESULTS

The kinematic response from muscle contraction was within 10% of the in vivo experimental data, while the BC method differed by 18%. The intervertebral disc forces from muscle contraction were agreeable with the literature (167 N compression, 12 N shear), while the BC methodology underpredicted the disc forces owing to the lack of spine compression. Correspondingly, the strains in the annulus fibrosus increased by an average of 60% across all levels due to muscle contraction compared to BC method.

CONCLUSION

The muscle repositioning method enhanced the kinetic response and subsequently led to differences in tissue-level responses compared to the conventional BC method. The improved kinematics and kinetics quantify the importance of repositioning FE neck models using active muscles to achieve non-neutral neck positions.

Assessment of Pedestrians' Head and Lower Limb Injuries in Tram-Pedestrian Collisions

Analysis of pedestrians' head and lower limb injuries at the tissue level is lacking in studies of tram-pedestrian collisions. The purpose of this paper therefore to investigate the impact response process and severity of pedestrians' injuries in tram-pedestrian collisions, using the Total Human Model for Safety (THUMS) pedestrian human body model together with the tram FE model. Two full-scale tram-pedestrian dummy crash tests were performed to validate the FE model, and the total correlation and analysis (CORA) score of head acceleration yielded values of 0.840 and 0.734, confirming a strong agreement between the FE-simulated head responses and the experimental head kinematics. The effects of different tram speeds and impact angles on pedestrians' impact response injuries and the differences were further analyzed. The results indicate that direct impact of the lower limb with the tram's obstacle deflector leads to lower limb bone shaft fractures and knee tissue damage. Neck fling contributed to worsened head injury. Coup contusions were the predominant type of brain contusion, surpassing contrecoup contusions, while diffuse axonal injury was mainly concentrated in the collision-side region of the brain. Pedestrians' injuries are influenced by tram velocity and impact angle: higher tram velocities increase the risk of lower limb and head injuries. The risk of head injury for pedestrians is higher when the impact angle is negative, while lower limb injuries are more significant when the impact angle is 0°. This study provides practical guidance for enhancing tram safety and protecting pedestrians.

Evaluation and Validation of Thorax Model Responses: A Hierarchical Approach to Achieve High Biofidelity for Thoracic Musculoskeletal System

As one of the most frequently occurring injuries, thoracic trauma is a significant public health burden occurring in road traffic crashes, sports accidents, and military events. The biomechanics of the human thorax under impact loading can be investigated by computational finite element (FE) models, which are capable of predicting complex thoracic responses and injury outcomes quantitatively. One of the key challenges for developing a biofidelic FE model involves model evaluation and validation. In this work, the biofidelity of a mid-sized male thorax model has been evaluated and enhanced by a multi-level, hierarchical strategy of validation, focusing on injury characteristics, and model improvement of the thoracic musculoskeletal system. At the component level, the biomechanical responses of several major thoracic load-bearing structures were validated against different relevant experimental cases in the literature, including the thoracic intervertebral joints, costovertebral joints, clavicle, sternum, and costal cartilages. As an example, the thoracic spine was improved by accurate representation of the components, material properties, and ligament failure features at tissue level then validated based on the quasi-static response at the segment level, flexion bending response at the functional spinal unit level, and extension angle of the whole thoracic spine. At ribcage and full thorax levels, the thorax model with validated bony components was evaluated by a series of experimental testing cases. The validation responses were rated above 0.76, as assessed by the CORA evaluation system, indicating the model exhibited overall good biofidelity. At both component and full thorax levels, the model showed good computational stability, and reasonable agreement with the experimental data both qualitatively and quantitatively. It is expected that our validated thorax model can predict thorax behavior with high biofidelity to assess injury risk and investigate injury mechanisms of the thoracic musculoskeletal system in various impact scenarios. The relevant validation cases established in this study shall be directly used for future evaluation of other thorax models, and the validation approach and process presented here may provide an insightful framework toward multi-level validating of human body models.

Load-Sharing and Kinematics of the Human Cervical Spine Under Multi-Axial Transverse Shear Loading: Combined Experimental and Computational Investigation

The cervical spine experiences shear forces during everyday activities and injurious events yet there is a paucity of biomechanical data characterizing the cervical spine under shear loading. This study aimed to (1) characterize load transmission paths and kinematics of the subaxial cervical spine under shear loading, and (2) assess a contemporary finite element cervical spine model using this data. Subaxial functional spinal units (FSUs) were subjected to anterior, posterior, and lateral shear forces (200 N) applied with and without superimposed axial compression preload (200 N) while monitoring spine kinematics. Load transmission paths were identified using strain gauges on the anterior vertebral body and lateral masses and a disc pressure sensor. Experimental conditions were simulated with cervical spine finite element model FSUs (GHBMC M50 version 5.0). The mean kinematics, vertebral strains, and disc pressures were compared to experimental results. The shear force-displacement response typically demonstrated a toe region followed by a linear response, with higher stiffness in anterior shear relative to lateral and posterior shear. Compressive axial preload decreased posterior and lateral shear stiffness and increased initial anterior shear stiffness. Load transmission patterns and kinematics suggest the facet joints play a key role in limiting anterior shear while the disc governs motion in posterior shear. The main cervical spine shear responses and trends are faithfully predicted by the GHBMC cervical spine model. These basic cervical spine biomechanics and the computational model can provide insight into mechanisms for facet dislocation in high severity impacts, and tissue distraction in low severity impacts.

A finite element-guided mathematical surrogate modeling approach for assessing occupant injury trends across variations in simplified vehicular impact conditions

A finite element (FE)-guided mathematical surrogate modeling methodology is presented for evaluating relative injury trends across varied vehicular impact conditions. The prevalence of crash-induced injuries necessitates the quantification of the human body's response to impacts. FE modeling is often used for crash analyses but requires time and computational cost. However, surrogate modeling can predict injury trends between the FE data, requiring fewer FE simulations to evaluate the complete testing range. To determine the viability of this methodology for injury assessment, crash-induced occupant head injury criterion (HIC₁₅) trends were predicted from Kriging models across varied impact velocities (10-45 mph; 16.1-72.4 km/h), locations (near side, far side, front, and rear), and angles (-45 to 45°) and compared to previously published data. These response trends were analyzed to locate high-risk target regions. Impact velocity and location were the most influential factors, with HIC₁₅ increasing alongside the velocity and proximity to the driver. The impact angle was dependent on the location and was minimally influential, often producing greater HIC₁₅ under oblique angles. These model-based head injury trends were consistent with previously published data, demonstrating great promise for the proposed methodology, which provides effective and efficient quantification of human response across a wide variety of car crash scenarios, simultaneously. This study presents a finite element-guided mathematical surrogate modeling methodology to evaluate occupant injury response trends for a wide range of impact velocities (10-45 mph), locations, and angles (-45 to 45°). Head injury response trends were predicted and compared to previously published data to assess the efficacy of the methodology for assessing occupant response to variations in impact conditions. Velocity and location were the most influential factors on the head injury response, with the risk increasing alongside greater impact velocity and locational proximity to the driver. Additionally, the angle of impact variable was dependent on the location and, thus, had minimal independent influence on the head injury risk.

Neck Muscle and Head/Neck Kinematic Responses While Bracing Against the Steering Wheel During Front and Rear Impacts

Drivers often react to an impending collision by bracing against the steering wheel. The goal of the present study was to quantify the effect of bracing on neck muscle activity and head/torso kinematics during low-speed front and rear impacts. Eleven seated subjects (3F, 8 M) experienced multiple sled impacts (Δv = 0.77 m/s; a_peak = 19.9 m/s², Δt = 65.5 ms) with their hands on the steering wheel in two conditions: relaxed and braced against the steering wheel. Electromyographic activity in eight neck muscles (sternohyoid, sternocleidomastoid, splenius capitis, semispinalis capitis, semispinalis cervicis, multifidus, levator scapulae, and trapezius) was recorded unilaterally with indwelling electrodes and normalized by maximum voluntary contraction (MVC) levels. Head and torso kinematics (linear acceleration, angular velocity, angular rotation, and retraction) were measured with sensors and motion tracking. Muscle and kinematic variables were compared between the relaxed and braced conditions using linear mixed models. We found that pre-impact bracing generated only small increases in the pre-impact muscle activity (< 5% MVC) when compared to the relaxed condition. Pre-impact bracing did not increase peak neck muscle responses during the impacts; instead it reduced peak trapezius and multifidus muscle activity by about half during front impacts. Bracing led to widespread changes in the peak amplitude and timing of the torso and head kinematics that were not consistent with a simple stiffening of the head/neck/torso system. Instead pre-impact bracing served to couple the torso more rigidly to the seat while not necessarily coupling the head more rigidly to the torso.

Material property calibration is more important than element size and number of different materials on the finite element modelling of vertebral bodies: A Taguchi study

Finite element (FE) modelling of a vertebral body (VB) is considered challenging due to the many parameters involved such as element size and type, and material properties. Previous studies have reported how these parameters affect the mechanical behaviour of a VB model; however, most studies just compared results without any specific statistical tool to quantify their influence. The Taguchi Method (TM) has been successfully used in manufacturing and biomechanics to evaluate process parameters and to determine optimum set-up conditions. This study aimed to evaluate the influence of the main finite element modelling parameters on the mechanical behaviour of a VB model using the Taguchi Method. A FE model was developed based on a C2 juvenile porcine vertebral body and three of the most commonly used modelling parameters were evaluated using TM in terms of the change in the predicted stiffness in comparison to experimental values: element size, number of different material properties for VB (based on grey-scale bins) and calibration factor for grey-scale to density to Young's Modulus equation. The influence of the combined factors was also assessed. The Taguchi analysis showed that the three factors are independent. The calibration factor is the main contributor, accounting for 97% of the predicted stiffness, with the value of 0.03 most closely aligning the numerical and experimental results. Element size accounted for 2% of the predicted stiffness, with 0.75 mm being the optimal, while the number of grey-scale bins influenced the results by less than 1%. Our findings indicate that the calibration factor is the main modelling parameter, with the element size and number of bins accounting for less than 3% of the predicted stiffness. Therefore, calibration of material properties should be done based on a large number of samples to ensure reliable results.

Application of Simulation Methods in Cervical Spine Dynamics

Neck injury is one of the most frequent spine injuries due to the complex structure of the cervical spine. The high incidence of neck injuries in collision accidents can bring a heavy economic burden to the society. Therefore, knowing the potential mechanisms of cervical spine injury and dysfunction is significant for improving its prevention and treatment. The research on cervical spine dynamics mainly concerns the fields of automobile safety, aeronautics, and astronautics. Numerical simulation methods are beneficial to better understand the stresses and strains developed in soft tissues with investigators and have been roundly used in cervical biomechanics. In this article, the simulation methods for the development and application of cervical spine dynamic problems in the recent years have been reviewed. The study focused mainly on multibody and finite element models. The structure, material properties, and application fields, especially the whiplash injury, were analyzed in detail. It has been shown that simulation methods have made remarkable progress in the research of cervical dynamic injury mechanisms, and some suggestions on the research of cervical dynamics in the future have been proposed.

Analysis of the mechanical response of damaged human cervical spine ligaments

BACKGROUND

Cervical spine ligaments that protect the spinal cord and stabilize the spine are frequently injured in motor vehicle collisions and other traumatic situations. These injuries are usually incomplete, and often difficult to notice. The focus of the presented study is placed on analysis of the effect of subfailure load on the mechanical response of the three main cervical spine ligaments: the anterior and the posterior longitudinal ligament and the ligamentum flavum.

METHODS

A total of 115 samples of human cadaveric ligaments removed within 24-48 h after death have been tested. Uniaxial tension tests along the fiber direction were performed in physiological conditions on a custom designed test equipment. The ligaments were loaded into an expected damage zone at two different subfailure values (based on previously reported reference group of 46 samples), and then reloaded to failure.

FINDINGS

The main effect of a high subfailure load has proven to be the toe elongation change. The toe elongation increase is affected by the subfailure load value. While anterior and posterior longitudinal ligament showed similar changes, the smallest subfailure effect was found in ligamentum flavum.

INTERPRETATIONS

The normal physiological region of the cervical spine ligaments mechanical response is modified by a high subfailure load. The observed ligament injury significantly compromises ligament ability to give tensile support within physiological spinal motion.

Rear-Impact Neck Whiplash: Role of Head Inertial Properties and Spine Morphological Variations on Segmental Rotations

Whiplash injuries continue to be a concern in low-speed rear impact. This study was designed to investigate the role of variations in spine morphology and head inertia properties on cervical spine segmental rotation in rear-impact whiplash loading. Vertebral morphology is rarely considered as an input parameter in spine finite element (FE) models. A methodology toward considering morphological variations as input parameters and identifying the influential variations is presented in this paper. A cervical spine FE model, with its morphology parametrized using mesh morphing, was used to study the influence of disk height, anteroposterior vertebral depth, and segmental size, as well as variations in head mass, moment of inertia, and center of mass locations. The influence of these variations on the characteristic S-curve formation in whiplash response was evaluated using the peak C2–C3 flexion marking the maximum S-curve formation and time taken for the formation of maximum S-curve. The peak C2–C3 flexion in the S-curve formation was most influenced by disk height and vertebral depth, followed by anteroposterior head center of mass location. The time to maximum S-curve was most influenced by the anteroposterior location of head center of mass. The influence of gender-dependent variations, such as the vertebral depth, suggests that they contribute to the greater segmental rotations observed in females resulting in different S-curve formation from men. These results suggest that both spine morphology and head inertia properties should be considered to describe rear-impact responses.

On the importance of retaining stresses and strains in repositioning computational biomechanical models of the cervical spine

Human body models are created in a specific posture and often repositioned and analyzed without retaining stresses that result from repositioning. For example, repositioning a human neck model within the physiological range of motion to a head-turned posture prior to an impact results in initial stresses within the tissues distracted from their neutral position. The aim of this study was to investigate the effect of repositioning on the subsequent kinetics, kinematics, and failure modes, of a lower cervical spine motion segment, to support future research at the full neck level. Repositioning was investigated for 3 modes (tension, flexion, and extension) and 3 load cases. The model was repositioned and loaded to failure in one continuous load history (case 1), or repositioned then restarted with retained stresses and loaded to failure (case 2). In case 3, the model was repositioned and then restarted in a stress-free state, representing current repositioning methods. Not retaining the repositioning stresses and strains resulted in different kinetics, kinematics, or failure modes, depending on the mode of loading. For the motion segment model, the differences were associated with the intervertebral disc fiber reorientation and load distribution, because the disc underwent the largest deformation during repositioning. This study demonstrated that repositioning led to altered response and tissue failure, which is critical for computational models intended to predict injury at the tissue level. It is recommended that stresses and strains be included and retained for subsequent analysis when repositioning a human computational neck model.

An investigation of dimensional scaling using cervical spine motion segment finite element models

The paucity of experimental data for validating computational models of different statures underscores the need for appropriate scaling methods so that models can be verified and validated using experimental data. Scaling was investigated using 50th percentile male (M50) and 5th percentile female (F05) cervical spine motion segment (C4-C5) finite element models subject to tension, flexion, and extension loading. Two approaches were undertaken: geometric scaling of the models to investigate size effects (volumetric scaling) and scaling of the force-displacement or moment-angle model results (data scaling). Three sets of scale factors were considered: global (body mass), regional (neck dimensions), and local (segment tissue dimensions). Volumetric scaling of the segment models from M50 to F05, and vice versa, produced correlations that were good or excellent in both tension and flexion (0.825-0.991); however, less agreement was found in extension (0.550-0.569). The reduced correlation in extension was attributed to variations in shape between the models leading to nonlinear effects such as different time to contact for the facet joints and posterior processes. Data scaling of the responses between the M50 and F05 models produced similar trends to volumetric scaling, with marginally greater correlations. Overall, the local tissue level and neck region level scale factors produced better correlations than the traditional global scaling. The scaling methods work well for a given subject, but are limited in applicability between subjects with different morphology, where nonlinear effects may dominate the response.

A Female Ligamentous Cervical Spine Finite Element Model Validated for Physiological Loads

Mathematical cervical spine models allow for studying of impact loading that can cause whiplash associated disorders (WAD). However, existing models only cover the male anthropometry, despite the female population being at a higher risk of sustaining WAD in automotive rear-end impacts. The aim of this study is to develop and validate a ligamentous cervical spine intended for biomechanical research on the effect of automotive impacts. A female model has the potential to aid the design of better protection systems as well as improve understanding of injury mechanisms causing WAD. A finite element (FE) mesh was created from surface data of the cervical vertebrae of a 26-year old female (stature 167 cm, weight 59 kg). Soft tissues were generated from the skeletal geometry and anatomical literature descriptions. Ligaments were modeled with nonlinear elastic orthotropic membrane elements, intervertebral disks as composites of nonlinear elastic bulk elements, and orthotropic anulus fibrosus fiber layers, while cortical and trabecular bones were modeled as isotropic plastic-elastic. The model has geometrical features representative of the female cervical spine-the largest average difference compared with published anthropometric female data was the vertebral body depth being 3.4% shorter for the model. The majority the cervical segments compare well with respect to biomechanical data at physiological loads, with the best match for flexion-extension loads and less biofidelity for axial rotation. An average female FE ligamentous cervical spine model was developed and validated with respect to physiological loading. In flexion-extension simulations with the developed female model and an existing average male cervical spine model, a greater range of motion (ROM) was found in the female model.

Injury mechanisms of the ligamentous cervical C2-C3 Functional Spinal Unit to complex loading modes: Finite Element study

The cervical spine sustains high rate complex loading modes during Motor Vehicle Crashes (MVCs) which may produce severe injuries accompanied with soft and/or hard tissue failure. Although previous numerical and experimental studies have provided insights on the cervical spine behavior under various loading scenarios, its response to complex impact loads and the resulting injury mechanisms are not fully understood. A validated Finite Element (FE) model of the ligamentous cervical C2-C3 Functional Spinal Unit (FSU) was utilized to assess the spinal response to six combined impact loading modes; flexion-extension combined with compression and distraction, and lateral bending and axial rotation combined with distraction. The FE model used time and rate-dependent material laws which permit assessing bone fracture and ligament failure. Spinal load-sharing, stresses in the spinal components, intradiscal pressure (IDP) change in the nucleus as well as contact pressure in the facet joints were predicted. Bone and ligaments failure occurrence and initiation instants were investigated. Results showed that spinal load-sharing varied with loading modes. Lateral bending combined with distraction was the most critical loading mode as it increased stresses and strains significantly and produced failure in most of the spinal components compared to other modes. The facet joints and surrounding cancellous bone as well as ligaments particularly the capsular (CL) and flavum (FL) ligaments were the most vulnerable structures to rapid flexion-extension, axial rotation and lateral bending combined with distraction or compression. The excessive stress and strain resulted from these loading modes produced rupture of the CL and FL ligaments and failure in the cancellous bone. The detection of failure initiation as well as fracture assessment demonstrated the vulnerability of ligaments to tensile combined loads and the major contribution of the bony structures in resisting compressive combined loads. Findings of this study may potentially assist in the development of injury prevention and treatment strategies.

Out-of-Position Rear Impact Tissue-Level Investigation Using Detailed Finite Element Neck Model

OBJECTIVE

Whiplash injuries can occur in automotive crashes and may cause long-term health issues such as neck pain, headache, and visual and auditory disturbance. Evidence suggests that nonneutral head posture can significantly increase the potential for injury in a given impact scenario, but epidemiological and experimental data are limited and do not provide a quantitative assessment of the increased potential for injury. Although there have been some attempts to evaluate this important issue using finite element models, none to date have successfully addressed this complex problem.

METHODS

An existing detailed finite element neck model was evaluated in nonneutral positions and limitations were identified, including musculature implementation and attachment, upper cervical spine kinematics in axial rotation, prediction of ligament failure, and the need for repositioning the model while incorporating initial tissue strains. The model was enhanced to address these issues and an iterative procedure was used to determine the upper cervical spine ligament laxities. The neck model was revalidated using neutral position impacts and compared to an out-of-position cadaver experiment in the literature. The effects of nonneutral position (axial head rotation) coupled with muscle activation were studied at varying impact levels.

RESULTS

The laxities for the ligaments of the upper cervical spine were determined using 4 load cases and resulted in improved response and predicted failure loads relative to experimental data. The predicted head response from the model was similar to an experimental head-turned bench-top rear impact experiment. The parametric study identified specific ligaments with increased distractions due to an initial head-turned posture and the effect of active musculature leading to reduced ligament distractions.

CONCLUSIONS

The incorporation of ligament laxity in the upper cervical spine was essential to predict range of motion and traumatic response, particularly for repositioning of the neck model prior to impact. The results of this study identify a higher potential for injury in out-of-position rear collisions and identified at-risk locations based on ligament distractions. The model predicted higher potential for injury by as much as 50% based on ligament distraction for the out-of-position posture and reduced potential for injury with muscle activation. Importantly, this study demonstrated that the location of injury or pain depends on the initial occupant posture, so that both the location of injury and kinematic threshold may vary when considering common head positions while driving.

Investigation of the Mass Distribution of a Detailed Seated Male Finite Element Model

Accurate mass distribution in computational human body models is essential for proper kinematic and kinetic simulations. The purpose of this study was to investigate the mass distribution of a 50th percentile male (M50) full body finite element model (FEM) in the seated position. The FEM was partitioned into 10 segments, using segment planes constructed from bony landmarks per the methods described in previous research studies. Body segment masses and centers of gravity (CGs) of the FEM were compared with values found from these studies, which unlike the present work assumed homogeneous body density. Segment masses compared well to literature while CGs showed an average deviation of 6.0% to 7.0% when normalized by regional characteristic lengths. The discrete mass distribution of the FEM appears to affect the mass and CGs of some segments, particularly those with low-density soft tissues. The locations of the segment CGs are provided in local coordinate systems, thus facilitating comparison with other full body FEMs and human surrogates. The model provides insights into the effects of inhomogeneous mass on the location of body segment CGs.