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

A benchmark of muscle models to length changes great and small

Digital human body models are used to simulate injuries that occur as a result of vehicle collisions, vibration, sports, and falls. Given enough time the body's musculature can generate force, affect the body's movements, and change the risk of some injuries. The finite-element code LS-DYNA is often used to simulate the movements and injuries sustained by the digital human body models as a result of an accident. In this work, we evaluate the accuracy of the three muscle models in LS-DYNA (MAT_156, EHTM, and the VEXAT) when simulating a range of experiments performed on isolated muscle: force-length-velocity experiments on maximally and sub-maximally stimulated muscle, active-lengthening experiments, and vibration experiments. The force-length-velocity experiments are included because these conditions are typical of the muscle activity that precedes an accident, while the active-lengthening and vibration experiments mimic conditions that can cause injury. The three models perform similarly during the maximally and sub-maximally activated force-length-velocity experiments, but noticeably differ in response to the active-lengthening and vibration experiments. The VEXAT model is able to generate the enhanced forces of biological muscle during active lengthening, while both the MAT_156 and EHTM produce too little force. In response to vibration, the stiffness and damping of the VEXAT model closely follows the experimental data while the MAT_156 and EHTM models differ substantially. The accuracy of the VEXAT model comes from two additional mechanical structures that are missing in the MAT_156 and EHTM models: viscoelastic cross-bridges, and an active titin filament. To help others build on our work we have made our simulation code publicly available.

Using muscle-tendon load limits to assess unphysiological musculoskeletal model deformation and Hill-type muscle parameter choice

Musculoskeletal simulations are a useful tool for improving our understanding of the human body. However, the physiological validity of predicted kinematics and forces is highly dependent upon the correct calibration of muscle parameters and the structural integrity of a model's internal skeletal structure. In this study, we show how ill-tuned muscle parameters and unphysiological deformations of a model's skeletal structure can be detected by using muscle elements as sensors with which modelling and parameterization inconsistencies can be identified through muscle and tendon strain injury assessment. To illustrate our approach, two modelling issues were recreated. First, a model repositioning simulation using the THUMS AM50 occupant model version 5.03 was performed to show how internal model deformations can occur during a change of model posture. Second, the muscle material parameters of the OpenSim gait2354 model were varied to illustrate how unphysiological muscle forces can arise if material parameters are inadequately calibrated. The simulations were assessed for muscle and tendon strain injuries using previously published injury criteria and a newly developed method to determine tendon strain injury threshold values. Muscle strain injuries in the left and right musculus pronator teres were detected during the model repositioning. This straining was caused by an unphysiologically large gap (12.92 mm) that had formed in the elbow joint. Similarly, muscle and tendon strain injuries were detected in the modified right-hand musculus gastrocnemius medialis of the gait2354 model where an unphysiological reduction of the tendon slack length introduced large pre-strain of the muscle-tendon unit. The results of this work show that the proposed method can quantify the internal distortion behaviour of musculoskeletal human body models and the plausibility of Hill-type muscle parameter choice via strain injury assessment. Furthermore, we highlight possible actions to avoid the presented issues and inconsistencies in literature data concerning the material characteristics of human tendons.

Predicting occupant head displacements in evasive maneuvers; tuning and comparison of a rotational based and a translational based neck muscle controller

Objective: Real-life car crashes are often preceded by an evasive maneuver, which can alter the occupant posture and muscle state. To simulate the occupant response in such maneuvers, human body models (HBMs) with active muscles have been developed. The aim of this study was to implement an omni-directional rotational head-neck muscle controller in the SAFER HBM and compare the bio-fidelity of the HBM with a rotational controller to the HBM with a translational controller, in simulations of evasive maneuvers. Methods: The rotational controller was developed using an axis-angle representation of head rotations, with x, y, and z components in the axis. Muscle load sharing was based on rotational direction in the simulation and muscle activity recorded in three volunteer experiments in these directions. The gains of the rotational and translational controller were tuned to minimize differences between translational and rotational head displacements of the HBM and volunteers in braking and lane change maneuvers using multi-objective optimizations. Bio-fidelity of the model with tuned controllers was evaluated objectively using CORrelation and Analysis (CORA). Results: The results indicated comparable performance for both controllers after tuning, with somewhat higher bio-fidelity for rotational kinematics with the translational controller. After tuning, good or excellent bio-fidelity was indicated for both controllers in the loading direction (forward in braking, and lateral in lane change), with CORA scores of 0.86-0.99 and 0.93-0.98 for the rotational and translational controllers, respectively. For rotational displacements, and translational displacements in the other directions, bio-fidelity ranged from poor to excellent, with slightly higher average CORA scores for the HBM with the translational controller in both braking and lane changing. Time-averaged muscle activity was within one standard deviation of time-averaged muscle activity from volunteers. Conclusion: Overall, the results show that when tuned, both the translational and rotational controllers can be used to predict the occupant response to an evasive maneuver, allowing for the inclusion of evasive maneuvers prior to a crash in evaluation of vehicle safety. The rotational controller shows potential in controlling omni-directional head displacements, but the translational controller outperformed the rotational controller. Thus, for now, the recommendation is to use the translational controller with tuned gains.

An investigation of tendon strains in jersey finger injury load cases using a finite element neuromuscular human body model

Introduction: A common hand injury in American football, rugby and basketball is the so-called jersey finger injury (JFI), in which an eccentric overextension of the distal interphalangeal joint leads to an avulsion of the connected musculus flexor digitorum profundus (FDP) tendon. In the field of automotive safety assessment, finite element (FE) neuromuscular human body models (NHBMs) have been validated and are employed to evaluate different injury types related to car crash scenarios. The goal of this study is to show, how such a model can be modified to assess JFIs by adapting the hand of an FE-NHBM for the computational analysis of tendon strains during a generalized JFI load case. Methods: A jersey finger injury criterion (JFIC) covering the injury mechanisms of tendon straining and avulsion was defined based on biomechanical experiments found in the literature. The hand of the Total Human Model for Safety (THUMS) version 3.0 was combined with the musculature of THUMS version 5.03 to create a model with appropriate finger mobility. Muscle routing paths of FDP and musculus flexor digitorum superficialis (FDS) as well as tendon material parameters were optimized using literature data. A simplified JFI load case was simulated as the gripping of a cylindrical rod with finger flexor activation levels between 0% and 100%, which was then retracted with the velocity of a sprinting college football player to forcefully open the closed hand. Results: The optimization of the muscle routing node positions and tendon material parameters yielded good results with minimum normalized mean absolute error values of 0.79% and 7.16% respectively. Tendon avulsion injuries were detected in the middle and little finger for muscle activation levels of 80% and above, while no tendon or muscle strain injuries of any kind occurred. Discussion: The presented work outlines the steps necessary to adapt the hand model of a FE-NHBM for the assessment of JFIs using a newly defined injury criterion called the JFIC. The injury assessment results are in good agreement with documented JFI symptoms. At the same time, the need to rethink commonly asserted paradigms concerning the choice of muscle material parameters is highlighted.

Development and verification of a physiologically motivated internal controller for the open-source extended Hill-type muscle model in LS-DYNA

Nowadays, active human body models are becoming essential tools for the development of integrated occupant safety systems. However, their broad application in industry and research is limited due to the complexity of incorporated muscle controllers, the long simulation runtime, and the non-regular use of physiological motor control approaches. The purpose of this study is to address the challenges in all indicated directions by implementing a muscle controller with several physiologically inspired control strategies into an open-source extended Hill-type muscle model formulated as LS-DYNA user-defined umat41 subroutine written in the Fortran programming language. This results in increased usability, runtime performance and physiological accuracy compared to the standard muscle material existing in LS-DYNA. The proposed controller code is verified with extensive experimental data that include findings for arm muscles, the cervical spine region, and the whole body. Selected verification experiments cover three different muscle activation situations: (1) passive state, (2) open-loop and closed-loop muscle activation, and (3) reflexive behaviour. Two whole body finite element models, the 50th percentile female VIVA OpenHBM and the 50th percentile male THUMS v5, are used for simulations, complemented by the simplified arm model extracted from the 50th percentile male THUMS v3. The obtained results are evaluated additionally with the CORrelation and Analysis methodology and the mean squared error method, showing good to excellent biofidelity and sufficient agreement with the experimental data. It was shown additionally how the integrated controller allows simplified mimicking of the movements for similar musculoskeletal models using the parameters transfer method. Furthermore, the Hill-type muscle model presented in this paper shows better kinematic behaviour even in the passive case compared to the existing one in LS-DYNA due to its improved damping and elastic properties. These findings provide a solid evidence base motivating the application of the enhanced muscle material with the internal controller in future studies with Active Human Body Models under different loading conditions.

Investigating sources for variability in volunteer kinematics in a braking maneuver, a sensitivity analysis with an active human body model

Occupant kinematics during evasive maneuvers, such as crash avoidance braking or steering, varies within the population. Studies have tried to correlate the response to occupant characteristics such as sex, stature, age, and BMI, but these characteristics explain no or very little of the variation. Therefore, hypothesis have been made that the difference in occupant response stems from voluntary behavior. The aim of this study was to investigate the effect from other sources of variability: in neural delay, in passive stiffness of fat, muscle tissues and skin, in muscle size and in spinal alignment, as a first step towards explaining the variability seen among occupants in evasive maneuvers. A sensitivity analysis with simulations of the SAFER Human Body Model in braking was performed, and the displacements from the simulations were compared to those of volunteers. The results suggest that the head and torso kinematics were most sensitive to spinal alignment, followed by muscle size. For head and torso vertical displacements, the range in model kinematics was comparable to the range in volunteer kinematics. However, for forward displacements, the included parameters only explain some of the variability seen in the volunteer experiment. To conclude, the results indicate that the variation in volunteer vertical kinematics could be partly attributed to the variability in human characteristics analyzed in this study, while these cannot alone explain the variability in forward kinematics. The results can be used in future tuning of HBMs, and in future volunteer studies, when further investigating the potential causes of the large variability seen in occupant kinematics in evasive maneuvers.

Evaluation of muscle strain injury severity in active human body models

Even though significant efforts in the field of injury detection with finite element active human body models (FE AHBMs) have been made, injuries of the muscle-tendon unit (MTU) have not yet been taken into consideration. Therefore, the goal of this study was to define a muscle strain injury criterion (MSIC) to evaluate the damage sustained by the musculature during muscle driven movement scenarios. The MSIC was derived from biomechanical tests found in the literature and the proposed threshold values were substantiated through a comparison to an estimate of the ultimate tensile strength of human skeletal muscle and the forces acting on the biceps femoris long head muscle during one sprinting gait cycle. The application of the MSIC to state-of-the-art FE AHBMs was demonstrated by evaluating the strain injury severity of selected neck muscles of a full-body AHBM during two seat rotation load cases. The results of the MSIC substantiation suggest that all three injury threshold values proposed in this work fall in a plausible corridor of forces acting on the MTU. The combined results of the AHBM simulations indicate that neither of the two examined seat rotations are likely to cause strain injury to the neck muscles and that the proposed MSIC can easily be applied to current AHBMs without further modification of the model architecture or the muscle parameters. The MSIC was also used to formulate a hypothesis on the aetiology of muscle strain injuries, through which it was demonstrated that material inhomogeneities in the MTU might be the cause for strain injuries sustained during otherwise physiological movements. This work is a first step in the direction of the definition of a wholistic injury criterion for the human skeletal muscle fibre.