Comparison of control strategies for the cervical muscles of an average female head-neck finite element model

A simulation-based study for optimizing proportional-integral-derivative controller gains for different control strategies of an active upper extremity model using experimental data

This study investigates the effect of PID controller gains, reaction time, and initial muscle activation values on active human model behavior while comparing three different control strategies. The controller gains and reaction delays were optimized using published experimental data focused on the upper extremity. The data describes the reaction of five male subjects in four tests based on two muscle states (relaxed and tensed) and two states of awareness (open and closed eye). The study used a finite element model of the left arm isolated from the Global Human Body Models Consortium (GHBMC) average male simplified occupant model for simulating biomechanical simulations. Major skeletal muscles of the arm were modeled as 1D beam elements and assigned a Hill-type muscle material. Angular position control, muscle length control, and a combination of both were used as a control strategy. The optimization process was limited to 4 variables; three Proportional-Integral-Derivative (PID) controller gains and one reaction delay time. The study assumed the relaxed and tensed condition require distinct sets of controller gains and initial activation and that the closed-eye simulations can be achieved by increasing the reaction delay parameter. A post-hoc linear combination of angle and muscle length control was used to arrive at the final combined control strategy. The premise was supported by variation in the controller gains depending on muscle state and an increase in reaction delay based on awareness. The CORA scores for open-eye relaxed, closed-eye relaxed, open-eye tensed, and closed-eye tensed was 0.95, 0.90, 0.95, and 0.77, respectively using the combined control strategy.

Neck stabilization through sensory integration of vestibular and visual motion cues

BACKGROUND

To counteract gravity, trunk motion, and other perturbations, the human head-neck system requires continuous muscular stabilization. In this study, we combine a musculoskeletal neck model with models of sensory integration (SI) to unravel the role of vestibular, visual, and muscle sensory cues in head-neck stabilization and relate SI conflicts and postural instability to motion sickness.

METHOD

A 3D multisegment neck model with 258 Hill-type muscle elements was extended with postural stabilization using SI of vestibular (semicircular and otolith) and visual (rotation rate, verticality, and yaw) cues using the multisensory observer model (MSOM) and the subjective vertical conflict model (SVC). Dynamic head-neck stabilization was studied using empirical datasets, including 6D trunk perturbations and a 4 m/s2 slalom drive inducing motion sickness.

RESULTS

Recorded head translation and rotation are well matched when using all feedback loops with MSOM or SVC or assuming perfect perception. A basic version of the model, including muscle, but omitting vestibular and visual perception, shows that muscular feedback can stabilize the neck in all conditions. However, this model predicts excessive head rotations in conditions with trunk rotation and in the slalom. Adding feedback of head rotational velocity sensed by the semicircular canals effectively reduces head rotations at mid-frequencies. Realistic head rotations at low frequencies are obtained by adding vestibular and visual feedback of head rotation based on the MSOM or SVC model or assuming perfect perception. The MSOM with full vision well captures all conditions, whereas the MSOM excluding vision well captures all conditions without vision. The SVC provides two estimates of verticality, with a vestibular estimate SVCvest, which is highly effective in controlling head verticality, and an integrated vestibular/visual estimate SVCint which can complement SVCvest in conditions with vision. As expected, in the sickening drive, SI models imprecisely estimate verticality, resulting in sensory conflict and postural instability.

CONCLUSION

The results support the validity of SI models in postural stabilization, where both MSOM and SVC provide credible results. The results in the sickening drive show imprecise sensory integration to enlarge head motion. This uniquely links the sensory conflict theory and the postural instability theory in motion sickness causation.

Finite element human body models with active reflexive muscles suitable for sex based whiplash injury prediction

Previous research has not produced a satisfactory resource to study reflexive muscle activity for investigating potentially injurious whiplash motions. Various experimental and computational studies are available, but none provided a comprehensive biomechanical representation of human response during rear impacts. Three objectives were addressed in the current study to develop female and male finite element human body models with active reflexive neck muscles: 1) eliminate the buckling in the lower cervical spine of the model observed in earlier active muscle controller implementations, 2) evaluate and quantify the influence of the individual features of muscle activity, and 3) evaluate and select the best model configuration that can be used for whiplash injury predictions. The current study used an open-source finite element model of the human body for injury assessment representing an average 50th percentile female anthropometry, together with the derivative 50th percentile male morphed model. Based on the head-neck kinematics and CORelation and Analyis (CORA) tool for evaluation, models with active muscle controller and parallel damping elements showed improved head-neck kinematics agreement with the volunteers over the passive models. It was concluded that this model configuration would be the most suitable for gender-based whiplash injury prediction when different impact severities are to be studied.

Implementation and calibration of active small female and average male human body models using low-speed frontal sled tests

OBJECTIVE

The objective of this study was to implement active muscles in a computationally efficient small female finite element model (54.1 kg, 149.9 cm) suitable for predicting occupant response during precrash braking and low-speed frontal sled tests. We further calibrate and compare its results against an average male model (78.4 kg, 174.9 cm) using the same developmental approach.

METHODS

The active female model (F05-OS + Active) was developed by adding active skeletal muscle elements (n = 232) to the Global Human Body Models Consortium (GHBMC) 5th percentile female simplified occupant model (F05-OS v2.3). The muscle properties and physiological cross-sectional area (PCSA) for each muscle were taken from the M50-OS + Active v2.3 model but PCSAs were mass scaled to a 5th percentile female. A total of 8 simulations were conducted; 2 acceleration pulses (1.0 g and 2.5 g), 2 models (F05-OS + Active and M50-OS + Active), and 2 muscle states (activation and control; e.g., no activation). Each model's kinematics and reaction forces were compared with experimental data. Occupant responses of 6 5th percentile female and 6 50th percentile male volunteers (n = 12 total) were used. The data depict occupant response in precrash braking and low-speed frontal sled tests in a rigid test buck. All procedures were reviewed and approved by the Virginia Tech institutional review board. Each volunteer was in a relaxed state before the applied acceleration.

RESULTS

The occupant peak forward excursion results of both active models reasonably match the volunteer data for both pulse severities. The differences between active and control models were found to be significant by Wilcoxon signed-rank test (p < .05). The reaction loads of the active and control models lie within the experimental corridors.

CONCLUSIONS

To the authors' knowledge, this study is the first to concurrently calibrate and compare equivalently developed computational models of females and males in precrash and low-speed impacts. The modeling approach is capable of capturing the varied kinematics observed in the relaxed condition, which may be an important factor in studies focused on the effects of low-g vehicle dynamics on the occupant position. Finally, the computationally efficient modeling approach is imperative given the long duration (>500 ms) of the events simulated.

Analysis of control strategies for VIVA OpenHBM with active reflexive neck muscles

Modeling muscle activity in the neck muscles of a finite element (FE) human body model can be based on two biological reflex systems. One approach is to approximate the Vestibulocollic reflex (VCR) function, which maintains the head orientation relative to a fixed reference in space. The second system tries to maintain the head posture relative to the torso, similar to the Cervicocolic reflex (CCR). Strategies to combine these two neck muscle controller approaches in a single head-neck FE model were tested, optimized, and compared to rear-impact volunteer data. The first approach, Combined-Control, assumed that both controllers simultaneously controlled all neck muscle activations. In the second approach, Distributed-Control, one controller was used to regulate activation of the superficial muscles while a different controller acted on deep neck muscles. The results showed that any muscle controller that combined the two approaches was less effective than only using one of VCR- or CCR-based systems on its own. A passive model had the best objective rating for cervical spine kinematics, but the addition of a single active controller provided the best response for both head and cervical spine kinematics. The present study demonstrates the difficulty in completely capturing representative head and cervical spine responses to rear-impact loading and identified a controller capturing the VCR reflex as the best candidate to investigate whiplash injury mechanisms through FE modeling.

Quantitative cervical spine injury responses in whiplash loading with a numerical method of natural neural reflex consideration

BACKGROUND AND OBJECTIVE

Neural reflex is hypothesized as a regulating step in spine stabilizing system. However, neural reflex control is still in its infancy to consider in the previous finite element analysis of head-neck system for various applications. The purpose of this study is to investigate the influences of neural reflex control on neck biomechanical responses, then provide a new way to achieve an accurate biomechanical analysis for head-neck system with a finite element model.

METHODS

A new FE head-neck model with detailed active muscles and spinal cord modeling was established and globally validated at multi-levels. Then, it was coupled with our previously developed neuromuscular head-neck model to analyze the effects of vestibular and proprioceptive reflexes on biomechanical responses of head-neck system in a typical spinal injury loading condition (whiplash). The obtained effects were further analyzed by comparing a review of epidemiologic data on cervical spine injury situations.

RESULT

The results showed that the active model (AM) with neural reflex control obviously presented both rational head-neck kinematics and tissue injury risk referring to the previous experimental and epidemiologic studies, when compared with the passive model (PM) without it. Tissue load concentration locations as well as stress/strain levels were both changed due to the muscle activation forces caused by neural reflex control during the whole loading process. For the bony structures, the AM showed a peak stress level accounting for only about 25% of the PM. For the discs, the stress concentrated location was transferred from C2-C6 in the PM to C4-C6 in the AM. For the spinal cord, the strain concentrated locations were transferred from C1 segment to around C4 segment when the effects of neural reflex control were implemented, while the gray matter and white matter peak strains were reduced to 1/3 and 1/2 of the PM, respectively. All these were well correlated with epidemiological studies on clinical cervical spine injuries.

CONCLUSION

In summary, the present work demonstrated necessity of considering neural reflex in FE analysis of a head-neck system as well as our model biofidelity. Overall results also verified the previous hypothesis and further quantitatively indicated that the muscle activation caused by neural reflex is providing a protection for the neck in impact loading by decreasing the strain level and changing the possible injury to lower spinal cord level to reduce injury severity.

‘Falling heads’: investigating reflexive responses to head–neck perturbations

Background

Reflexive responses to head–neck perturbations affect the injury risk in many different situations ranging from sports-related impact to car accident scenarios. Although several experiments have been conducted to investigate these head–neck responses to various perturbations, it is still unclear why and how individuals react differently and what the implications of these different responses across subjects on the potential injuries might be. Therefore, we see a need for both experimental data and biophysically valid computational Human Body Models with bio-inspired muscle control strategies to understand individual reflex responses better.

Methods

To address this issue, we conducted perturbation experiments of the head–neck complex and used this data to examine control strategies in a simulation model. In the experiments, which we call ’falling heads’ experiments, volunteers were placed in a supine and a prone position on a table with an additional trapdoor supporting the head. This trapdoor was suddenly released, leading to a free-fall movement of the head until reflexive responses of muscles stopped the downwards movement.

Results

We analysed the kinematic, neuronal and dynamic responses for all individuals and show their differences for separate age and sex groups. We show that these results can be used to validate two simple reflex controllers which are able to predict human biophysical movement and modulate the response necessary to represent a large variability of participants.

Conclusions

We present characteristic parameters such as joint stiffness, peak accelerations and latency times. Based on this data, we show that there is a large difference in the individual reflexive responses between participants. Furthermore, we show that the perturbation direction (supine vs. prone) significantly influences the measured kinematic quantities. Finally, ’falling heads’ experiments data are provided open-source to be used as a benchmark test to compare different muscle control strategies and to validate existing active Human Body Models directly.

Importance of the cervical capsular joint cartilage geometry on head and facet joint kinematics assessed in a Finite element neck model

Finite element human neck models (NMs) aim to predict neck response and injury at the tissue level; however, contemporary models are most often assessed using global response such as head kinematics. Additionally, many NMs are developed from subject-specific imaging with limited soft tissue resolution in small structures such as the facet joints in the neck. Such details may be critical to enable NMs to predict tissue-level response. In the present study, the capsular joint cartilage (CJC) geometry in a contemporary NM was enhanced (M50-CJC) based on literature data. The M50-CJC was validated at the segment and full neck levels and assessed using relative facet joint kinematics (FJK), capsular ligament (CL) and intervertebral disc (IVD) strains, a relative vertebral rotation assessment (IV-NIC) and head kinematics in frontal and rear impact. The validation ratings at the segment level increased from 0.60 to 0.64, with improvements for modes of deformation associated with the facet joints, while no difference was noted at the head kinematic level. The improved CJC led to increased FJK rotation (188%) and IVD strain (152.2%,) attributed to the reduced facet joint gap. Further enhancements of the capsular joint representation or a link between the FJK and CL injury risk are recommended. Enhancements at the tissue level demonstrated a large effect on the IVD strain, but were not apparent in global metrics such as head kinematics. This study demonstrated that a biofidelic and detailed geometrical representation of the CJC contributes significantly to the predicted joint response, which is critical to investigate neck injury risk at the tissue level.

Estimating the Effects of Awareness on Neck-Muscle Loading in Frontal Impacts with EMG and MC Sensors

Critical traffic situations, such as vehicle collisions and emergency manoeuvres, can cause an occupant to respond with reflex and voluntary actions. These affect the occupant's position and dynamic loading during interactions with the vehicle's restraints, possibly compromising their protective function. Electromyography (EMG) is a commonly used method for measuring active muscle response and can also provide input parameters for computer simulations with models of the human body. The recently introduced muscle-contraction (MC) sensor is a wearable device with a piezo-resistive element for measuring the force of an indenting tip pressing against the surface of the body. The study aimed to compare how data collected simultaneously with EMG, video motion capture, and the novel MC sensor are related to neck-muscle loading. Sled tests with low-severity frontal impacts were conducted, assuming two different awareness conditions for seated volunteers. The activity of the upper trapezius muscle was measured using surface EMG and MC sensors. The neck-muscle load F was estimated from an inverse dynamics analysis of the head's motion captured in the sagittal plane. The volunteers' response to impact was predominantly reflexive, with significantly shorter onset latencies and more bracing observed when the volunteers were aware of the impact. Cross-correlations between the EMG and MC, EMG and F, and F and MC data were not changed significantly by the awareness conditions. The MC signal was strongly correlated (r = 0.89) with the neck-muscle loading F in the aware and unaware conditions, while the mean ΔF-MC delays were 21.0 ± 15.1 ms and 14.6 ± 12.4 ms, respectively. With the MC sensor enabling a consistent measurement-based estimation of the muscle loading, the simultaneous acquisition of EMG and MC signals improves the assessment of the reflex and voluntary responses of a vehicle's occupant subjected to low-severity loading.

Optimization of Female Head-Neck Model with Active Reflexive Cervical Muscles in Low Severity Rear Impact Collisions

ViVA Open Human Body Model (HBM) is an open-source human body model that was developed to fill the gap of currently available models that lacked the average female size. In this study, the head–neck model of ViVA OpenHBM was further developed by adding active muscle controllers for the cervical muscles to represent the human neck muscle reflex system as studies have shown that cervical muscles influence head–neck kinematics during impacts. The muscle controller was calibrated by conducting optimization-based parameter identification of published-volunteer data. The effects of different calibration objectives to head–neck kinematics were analyzed and compared. In general, a model with active neck muscles improved the head–neck kinematics agreement with volunteer responses. The current study highlights the importance of including active muscle response to mimic the volunteer’s kinematics. A simple PD controller has found to be able to represent the behavior of the neck muscle reflex system. The optimum gains that defined the muscle controllers in the present study were able to be identified using optimizations. The present study provides a basis for describing an active muscle controller that can be used in future studies to investigate whiplash injuries in rear impacts