Abstract

Digital image correlation (DIC) is rapidly increasing in popularity in biomechanical studies of the musculoskeletal system. DIC allows the re-construction of full field displacement and/or strain maps of the surface of an object. DIC systems typically consist of two cameras focussing on the same region of interest. This constrains the angle between the cameras to be relatively narrow when studying specimens characterised by complex geometrical features, giving rise to concerns on the accuracy of the out of plane estimates of movement.

The aim of this research was to compare the movement profiles of bony segments measured by DIC and by an optoelectronic motion capture system.

Five porcine cervical spine segments (C2-C6) were obtained from the local butcher. These were stripped of all anterior soft tissues while the posterior structures were left intact. A speckle pattern was applied to the anterior aspect of the specimens, while custom made infrared clusters were rigidly attached to the 3 middle vertebral bodies (C3-C5). The specimens were mounted in a custom made impact rig which fully constrained C6 but allowed C2 to translate in the axial direction of the segment. Images were acquired at 4kHz, both for the DIC (Photron Europe Ltd, UK) and motion capture cameras (Qualisys Oqus 400, Sweden). The in-plane and out of plane displacements of each of the VBs were plotted as a function of time and the similarity between the curves thus obtained was analysed using the SPM1D technique which allowed a comparison to be made in terms of t-statistics. No statistical differences were found between the two techniques in all axis of movement, however the out of plane movements were characterised by higher variance which is attributed to the uncertainty arising from the near parallel positioning of the cameras in the experimental set-up.

Head collisions in sport can result in catastrophic cervical spine injuries. Musculo-skeletal (MSK) modelling can help analyse the relationship between players' motion, external loading and internal stresses that lead to injury. However, the literature lacks sport specific MSK models. In automotive research the intervertebral disc behaviour has been represented as viscoelastic elements (“bushing”), whose stiffness and damping parameters are often estimated under quasi-static conditions and may lack validity in dynamic impacts. The aim of this study was to develop a validated cervical spine model for axial impacts for future use in the analysis of head-first rugby collisions.

A drop test rig was used to replicate a sub-catastrophic axial head impact. A load of 80 N from 0.5 m was applied to the cranial aspect of a C2-C6 porcine spinal specimen mounted in the neutral position. The 3D motion of C3-C5 vertebras (4 kHz) and the cranial axial load (1 MHz) were measured via motion capture (Qualysis, Sweden) and a uniaxial load cell (RDP Electronics Ltd, UK). Specimen specific models were created in NMSBuilder and OpenSim after the vertebrae geometries were obtained from the segmentation of micro-CT images of the specimens. The compressive viscoelastic properties of four vertebral joints (C2-C3 through to C5-C6) were optimised via a Genetic Algorithm (MATLAB v2016b, The Mathworks Inc) to minimise tracking errors.

The optimisation converged to a solution of 140–49000 kN/m and 2000–8000 Ns/m for stiffness and damping respectively (RMSE=5.1 mm). Simulated joint displacements ranged between 0.09 – 1.75 mm compared to experimental 0.1 – 0.8 mm.

Optimal bushing parameters were higher than previously reported values measured through quasi-static testing. Higher stiffness and damping values could be explained by the higher-dynamics nature of the event analysed related to a different part of the non-linear intervertebral disc load-displacement curve.