Multilevel Validation of a Male Neck Finite Element Model With Active Musculature

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.

Subaxial Cervical Spine Motion With Different Sizes of Head-supported Mass Under Accelerative Forces

INTRODUCTION

The evolution of military helmet devices has increased the amount of head-supported mass (HSM) worn by warfighters. HSM has important implications for spine biomechanics, and yet, there is a paucity of studies that investigated the effects of differing HSM and accelerative profiles on spine biomechanics. The aim of this study is to investigate the segmental motions in the subaxial cervical spine with different sizes of HSM under Gx accelerative loading.

METHODS

A three-dimensional finite element model of the male head-neck spinal column was used. Three different size military helmets were modeled and incorporated into head-neck model. The models were exercised under Gx accelerative loading by inputting low and high pulses to the cervical vertebra used in the experimental studies. Segmental motions were obtained and normalized with respect to the non-HSM case to quantify the effect of HSM.

RESULTS

Segmental motions increased with an increase in velocity at all segments of the spine. Increasing helmet size resulted in larger motion increases. Angulations ranged from 0.9° to 9.3° at 1.8 m/s and from 1.3° to 10.3° at 2.6 m/s without a helmet. Helmet increased motion between 5% to 74% at 1.8 m/s. At 2.6 m/s, the helmet increased segmental motion anywhere from 10% to 105% in the subaxial cervical spine. The greatest motion was seen at the C5-C6 level, followed by the C6-C7 level.

CONCLUSIONS

The subaxial cervical spine experiences motion increases at all levels at both velocity profiles with increasing HSM. Larger helmet and greater impact velocity increased motion at all levels, with C5-C6 demonstrating the largest range of motion. HSM should be minimized to reduce the risk of cervical spine injury to the warfighter.

Effect of muscle pre-tension and pre-impact neck posture on the kinematic response of the cervical spine in simulated low-speed rear impacts

Computational human body models (HBMs) can identify potential injury pathways not easily accessible through experimental studies, such as whiplash induced injuries. However, previous computational studies investigating neck response to simulated impact conditions have neglected the effect of pre-impact neck posture and muscle pre-tension on the intervertebral kinematics and tissue-level response. The purpose of the present study was addressing this knowledge gap using a detailed neck model subjected to simulated low-acceleration rear impact conditions, towards improved intervertebral kinematics and soft tissue response for injury assessment. An improved muscle path implementation in the model enabled the modeling of muscle pre-tension using experimental muscle pre-stretch data determined from previous cadaver studies. Cadaveric neck impact tests and human volunteer tests with the corresponding cervical spine posture were simulated using a detailed neck model with the reported boundary conditions and no muscle activation. Computed intervertebral kinematics of the model with pre-tension achieved, for the first time, the S-shape behavior of the neck observed in low severity rear impacts of both cadaver and volunteer studies. The maximum first principal strain in the muscles for the model with pre-tension was 27% higher than that without pre-tension. Although, the pre-impact neck posture was updated to match the average posture reported in the experimental tests, the change in posture was generally small with only small changes in vertebral kinematics and muscle strain. This study provides a method to incorporate muscle pre-tension in HBM and quantifies the importance of pre-tension in calculating tissue-level distractions.

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

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

The Influence of Mattress Stiffness on Spinal Curvature and Intervertebral Disc Stress-An Experimental and Computational Study

Sleeping support systems can influence spinal curvature, and the misalignment of the spinal curvature can lead to musculoskeletal problems. Previous sleep studies on craniocervical support focused on pillow variants, but the mattress supporting the pillow has rarely been considered. This study used a cervical pillow and three mattresses of different stiffnesses, namely soft, medium, and hard, with an indentation load deflection of 20, 42, and 120 lbs, respectively. A novel electronic curvature measurement device was adopted to measure the spinal curvature, whereby the intervertebral disc loading was computed using the finite element method. Compared with the medium mattress, the head distance increased by 30.5 ± 15.9 mm, the cervical lordosis distance increased by 26.7 ± 14.9 mm, and intervertebral disc peak loading increased by 49% in the soft mattress environment. Considering that the pillow support may increase when using a soft mattress, a softer or thinner pillow is recommended. The head distance and cervical lordosis distance in the hard mattress environment were close to the medium mattress, but the lumbar lordosis distance reduced by 10.6 ± 6.8 mm. However, no significant increase in intervertebral disc loading was observed, but contact pressure increased significantly, which could cause discomfort and health problems.

Anteroposterior shear stiffness of the upper thoracic spine at quasi-static and dynamic loading rates-An in vitro biomechanical study

To evaluate the biomechanical properties of the upper thoracic spine in anterior-posterior shear loading at various displacement rates. These data broaden our understanding of thoracic spine biomechanics and inform efforts to model the spine and spinal cord injuries. Seven T1-T2 thoracic functional spinal units were loaded non-destructively by a pure shear force up to 200 N, starting from a neutral posture. Tests were run in both posterior and anterior directions, at displacement rates of 1, 10, and 100 mm/s. The three-dimensional motion of the specimen was recorded at 1000 Hz. Individual and averaged load-displacement curves were generated and specimen stiffnesses were calculated. Due to a nonlinear response of the specimens, stiffness was defined separately for both the lower half and the upper half of the specimen range of motion. Specimens were significantly stiffer in the anterior direction than in the posterior direction, across all rates. At low displacements, the anterior stiffness averaged 230 N/mm, 76% higher than the low displacement posterior stiffness of 131 N/mm. At high displacements, anterior stiffness averaged 258 N/mm, 51% stiffer than the high displacement posterior stiffness of 171 N/mm. Shear displacement rate had a small effect on the load response, with the 100 mm/s rate causing a mildly stiffer response at low displacements in the anterior direction. Overall, the load-displacement response exhibited pseudo-quadratic behavior at 1 and 10 mm/s but became more linear at 100 mm/s. The shear stiffness in the upper thoracic spine is greatest in the anterior loading direction, being 51%-76% greater than posterior, most likely due to facet interactions. The effect of the shear displacement rate is low.

Sex, Age and Stature Affects Neck Biomechanical Responses in Frontal and Rear Impacts Assessed Using Finite Element Head and Neck Models

The increased incidence of injury demonstrated in epidemiological data for the elderly population, and females compared to males, has not been fully understood in the context of the biomechanical response to impact. A contributing factor to these differences in injury risk could be the variation in geometry between young and aged persons and between males and females. In this study, a new methodology, coupling a CAD and a repositioning software, was developed to reposture an existing Finite element neck while retaining a high level of mesh quality. A 5th percentile female aged neck model (F05_75YO) and a 50th percentile male aged neck model (M50_75YO) were developed from existing young (F05_26YO and M50_26YO) neck models (Global Human Body Models Consortium v5.1). The aged neck models included an increased cervical lordosis and an increase in the facet joint angles, as reported in the literature. The young and the aged models were simulated in frontal (2, 8, and 15 g) and rear (3, 7, and 10 g) impacts. The responses were compared using head and relative facet joint kinematics, and nominal intervertebral disc shear strain. In general, the aged models predicted higher tissue deformations, although the head kinematics were similar for all models. In the frontal impact, only the M50_75YO model predicted hard tissue failure, attributed to the combined effect of the more anteriorly located head with age, when compared to the M50_26YO, and greater neck length relative to the female models. In the rear impacts, the F05_75YO model predicted higher relative facet joint shear compared to the F05_26YO, and higher relative facet joint rotation and nominal intervertebral disc strain compared to the M50_75YO. When comparing the male models, the relative facet joint kinematics predicted by the M50_26YO and M50_75YO were similar. The contrast in response between the male and female models in the rear impacts was attributed to the higher lordosis and facet angle in females compared to males. Epidemiological data reported that females were more likely to sustain Whiplash Associated Disorders in rear impacts compared to males, and that injury risk increases with age, in agreement with the findings in the present study. This study demonstrated that, although the increased lordosis and facet angle did not affect the head kinematics, changes at the tissue level were considerable (e.g., 26% higher relative facet shear in the female neck compared to the male, for rear impact) and relatable to the epidemiology. Future work will investigate tissue damage and failure through the incorporation of aged material properties and muscle activation.

A Hyper-Viscoelastic Continuum-Level Finite Element Model of the Spinal Cord Assessed for Transverse Indentation and Impact Loading

Finite Element (FE) modelling of spinal cord response to impact can provide unique insights into the neural tissue response and injury risk potential. Yet, contemporary human body models (HBMs) used to examine injury risk and prevention across a wide range of impact scenarios often lack detailed integration of the spinal cord and surrounding tissues. The integration of a spinal cord in contemporary HBMs has been limited by the need for a continuum-level model owing to the relatively large element size required to be compatible with HBM, and the requirement for model development based on published material properties and validation using relevant non-linear material data. The goals of this study were to develop and assess non-linear material model parameters for the spinal cord parenchyma and pia mater, and incorporate these models into a continuum-level model of the spinal cord with a mesh size conducive to integration in HBM. First, hyper-viscoelastic material properties based on tissue-level mechanical test data for the spinal cord and hyperelastic material properties for the pia mater were determined. Secondly, the constitutive models were integrated in a spinal cord segment FE model validated against independent experimental data representing transverse compression of the spinal cord-pia mater complex (SCP) under quasi-static indentation and dynamic impact loading. The constitutive model parameters were fit to a quasi-linear viscoelastic model with an Ogden hyperelastic function, and then verified using single element test cases corresponding to the experimental strain rates for the spinal cord (0.32–77.22 s⁻¹) and pia mater (0.05 s⁻¹). Validation of the spinal cord model was then performed by re-creating, in an explicit FE code, two independent ex-vivo experimental setups: 1) transverse indentation of a porcine spinal cord-pia mater complex and 2) dynamic transverse impact of a bovine SCP. The indentation model accurately matched the experimental results up to 60% compression of the SCP, while the impact model predicted the loading phase and the maximum deformation (within 7%) of the SCP experimental data. This study quantified the important biomechanical contribution of the pia mater tissue during spinal cord deformation. The validated material models established in this study can be implemented in computational HBM.

Importance of passive muscle, skin, and adipose tissue mechanical properties on head and neck response in rear impacts assessed with a finite element model

OBJECTIVE

The objective of this study was to improve head-neck kinematic predictions of a contemporary finite element (FE) head-neck model, assessed in rear impact scenarios (3-10 g), by including an accurate representation of the skin, adipose tissue, and passive muscle mechanical properties. The soft tissues of the neck have a substantial contribution to kinematic response, with the contribution being inversely proportional to the impact severity. Thus accurate representation of these passive tissues is critical for the assessment of kinematic response and the potential for crash induced injuries. Contemporary Human Body Models (HBMs) often incorporate overly stiff mechanical properties of passive tissues for numerical stability, which can affect the predicted kinematic response of the head and neck.

METHODS

Soft tissue material properties including non-linearity, compression-tension asymmetry, and viscoelasticity were implemented in constitutive models for the skin, adipose, and passive muscle tissues, based on experimental data in the literature. A quasi-linear viscoelastic formulation was proposed for the skin, while a phenomenological hyper-viscoelastic model was used for the passive muscle and adipose tissues. A head-neck model extracted from a contemporary FE HBM was updated to include the new tissue models and assessed using head rotation angle for rear impact scenarios (3 g, 7 g, and 10 g peak accelerations), and compared to postmortem human surrogate (PMHS) data for 7 g impacts.

RESULTS

The head rotation angle increased with the new material models for all three rear impact cases: (3 g: +43%, 7 g: +52%, 10 g: +71%), relative to the original model. The increase in head rotation was primarily attributed to the improved skin model, with the passive muscle being a secondary contributor to the increase in response. A 52% increase in head rotation for the 7 g impact improved the model response with respect to PMHS data, placing it closer to the experimental average, compared to the original model.

CONCLUSIONS

The improved skin, adipose tissue, and passive muscle material model properties, based on published experimental data, increased the neck compliance in rear impact, with improved correspondence to published PMHS test data for medium severity impacts. Future studies will investigate the coupled effect of passive and active muscle tissue for low severity impacts.

Vestibulocollic and Cervicocollic Muscle Reflexes in a Finite Element Neck Model During Multidirectional Impacts

Active neck musculature plays an important role in the response of the head and neck during impact and can affect the risk of injury. Finite element Human Body Models (HBM) have been proposed with open and closed-loop controllers for activation of muscle forces; however, controllers are often calibrated to specific experimental loading cases, without considering the intrinsic role of physiologic muscle reflex mechanisms under different loading conditions. This study aimed to develop a single closed-loop controller for neck muscle activation in a contemporary male HBM based on known reflex mechanisms and assess how this approach compared to current open-loop controllers across a range of impact directions and severities. Controller parameters were optimized using volunteer data and independently assessed across twelve impact conditions. The kinematics from the closed-loop controller simulations showed good average CORA rating to the experimental data (0.699) for the impacts following the ISO/TR9790 standard. Compared to previously optimized open-loop activation strategy, the average difference was less than 9%. The incorporation of the reflex mechanisms using a closed-loop controller can provide robust performance for a range of impact directions and severities, which is critical to improving HBM response under a larger spectrum of automotive impact simulations.