• Gill, H. S.
• Cazzola, Dario
• Preatoni, Ezio

Abstract

INTRODUCTION<br/>Musculoskeletal models of the cervical spine are valuable tools to assess the intervertebral loads and injury mechanisms of sporting head impacts and provide an intermediate step between in-vivo experimental trials and detailed finite element analysis. A previous musculoskeletal model created for estimating intervertebral joint loading during rugby contact [1] events included kinematic constraints, which neglected individual joint translation and therefore may oversimplify the behavior of the cervical spine during impacts. In order to both confidently estimate the internal loading and resulting vertebral kinematics, intervertebral viscoelastic (“bushing”) elements should be included in the model. Such viscoelastic elements are used to describe the behaviour of the joints’ passive structures [2] and allow for resulting motions to be calculated compared to kinematic constraints. A key challenge is to identify the optimal bushing element that can generate reliable kinematics of the model’s cervical spine with respect to the load applied. Therefore, the aim of this study was to assess the performance of anteroposterior shear viscoelastic elements, that were validated together with compressive elements for impulsive axial loading conditions, under loads containing shear components. <br/><br/>METHODS<br/>Axial and anteroposterior viscoelastic bushing parameters were estimated via an optimisation procedure in a combined in-vitro and in-silico study. The optimisation procedure that identified the parameters included forward dynamic simulations of five specimen specific models that were driven by the experimental axial loads in order to replicate the in-vitro experiment. We carried out a 1000 sample Monte-Carlo analysis that perturbed the viscoelastic parameters by ±50% from their optimised values to evaluate the models’ sensitivity to a range of parameter values under the same axial loads. This method classified which parameters had the greatest effect on the models’ performance. The initial parameters estimated by the optimisation procedure were then implemented in the C2-C3, C3-C4, C4-C5 and C5-C6 joints of the population specific musculoskeletal model [1]. Forward dynamic simulations were executed with the updated musculoskeletal model of the cervical spine and head positioned parallel to the horizontal. An initial simulation with a constant anterior shear load of 250 N applied to the base of the head was completed to compare the model’s anteroposterior joint displacement against available experimental data. A second simulation investigated a possible injurious loading scenario that could be experienced in a head first impact in rugby. Experimental loads were collected from an instrumented ATD headform during live scrummaging trials. The compression flexion load was scaled to included maximum values of 1100 N and 700 N for anterior shear and compression respectively. A final simulation applied impulses with higher maximal values of 2500 N and 1500 N for shear and compressive forces respectively after the experimental load was scaled further. For all simulations, peak values of the applied loads were reached in 60 ms after an initial 40 ms latent period. No muscle activation was prescribed to reduce the active muscle force effects but their passive contribution remained.<br/><br/>RESULTS AND DISCUSSION<br/>The values of the bushing parameters identified by the optimisation increased from their respective initialised values and were 23.1 MN/m for axial stiffness, 4300 Ns/m for axial damping, 75.9 kN/m for shear stiffness and 1400 Ns/m for shear damping. Perturbations in shear damping and axial stiffness had the most effect on model performance. Decreasing values of shear damping resulted in larger error in the sensitivity analysis (Figure 1). Axial stiffness also showed an effect of increasing or decreasing values however this effect was smaller than the one observed in shear damping. This is possibly caused by the considerably higher values of axial stiffness (23.1 MN/m) compared to shear damping (1400 Ns/m), which, even when perturbed by – 50 %, remain in the same order of magnitude. On the contrary, the same proportional perturbation of the optimised shear damping reduces the value by an order of magnitude.During these impulsive axial loads applied to the cervical spines in the experiments a buckling response was observed as previously observed by [3]. The anterior shear motion of the vertebrae caused by the loading of the specimens in the current study, however, did not lead to injuries because the applied load was chosen to be sub-catastrophic. The sensitivity analysis thus highlighted the importance of correct choice for anteroposterior joint damping parameters used in musculoskeletal model of the cervical as lower values chosen for axial impacts may result in misrepresentative kinematic responses of the model cervical joints in anteroposterior translation. <br/><br/>The three forward d...

Topics

• impedance spectroscopy
• experiment
• simulation
• mass spectrometry
• activation
• finite element analysis