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

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.

Implications of range of motion requirements for the laxity of ligaments in a lumbar finite element model

Finite element models of the lumbar spine often adopt ligament properties from tensile tests without accounting for possible differences between testing and in situ initial ligament length. Such differences could result in laxities or preloads at the beginning of a simulation that would affect the ligament forces, tangent stiffness, and the posture at which they fail. In vivo and in vitro human experimental data reported laxities or preloads. However, laxities or preloads, which could also result from postural differences, are often neglected in simulation studies. This study proposes a numerical methodology to identify ranges of ligament laxities or preloads compatible with the selected tensile ligament properties, the model, and the range of motion (RoM) the model aims to simulate. The approach assumes that ligaments should remain in a safe elongation range for the complete RoM, and that each ligament should play a significant mechanical role in at least one load case. The methodology was applied to the functional spinal unit (FSU) models using the RoM from healthy subjects and ligament properties from the literature. Without laxity, some ligaments reached their elongation at failure within the RoM. Laxity ranges varied considerably (from -9.2 mm preload to 10.7 mm laxity) and flexion was the most critical load case to determine them. Their effect on the mobility response was also assessed. The effect on the mobility of a FSU was also assessed. While the proposed method cannot determine an exact laxity value, it is simple and it can be applied to any model to identify a plausible range of ligament initial length.

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.