Assignment 2
RESEARCH ARTICLE
Cervical Spine Injuries: A Whole-Body
Musculoskeletal Model for the Analysis of
Spinal Loading
Dario Cazzola 1*, Timothy P. Holsgrove2,3, Ezio Preatoni1, Harinderjit S. Gill2,
Grant Trewartha 1
1 Department for Health, University of Bath, Bath, United Kingdom, 2 Centre for Orthopaedic Biomechanics,
Department of Mechanical Engineering, University of Bath, Bath, United Kingdom, 3 College of Engineering,
Mathematics & Physical Sciences, University of Exeter, Exeter, United Kingdom
* d.cazzola@bath.ac.uk
Abstract
Cervical spine trauma from sport or traffic collisions can have devastating consequences
for individuals and a high societal cost. The precise mechanisms of such injuries are still
unknown as investigation is hampered by the difficulty in experimentally replicating the con-
ditions under which these injuries occur. We harness the benefits of computer simulation to
report on the creation and validation of i) a generic musculoskeletal model (MASI) for the
analyses of cervical spine loading in healthy subjects, and ii) a population-specific version of
the model (Rugby Model), for investigating cervical spine injury mechanisms during rugby
activities. The musculoskeletal models were created in OpenSim, and validated against in
vivo data of a healthy subject and a rugby player performing neck and upper limb move-
ments. The novel aspects of the Rugby Model comprise i) population-specific inertial proper-
ties and muscle parameters representing rugby forward players, and ii) a custom scapula-
clavicular joint that allows the application of multiple external loads. We confirm the utility of
the developed generic and population-specific models via verification steps and validation of
kinematics, joint moments and neuromuscular activations during rugby scrummaging and
neck functional movements, which achieve results comparable with in vivo and in vitro data.
The Rugby Model was validated and used for the first time to provide insight into anatomical
loading and cervical spine injury mechanisms related to rugby, whilst the MASI introduces a
new computational tool to allow investigation of spinal injuries arising from other sporting
activities, transport, and ergonomic applications. The models used in this study are freely
available at simtk.org and allow to integrate in silico analyses with experimental approaches
in injury prevention.
1. Introduction
The population incidence of acute spinal cord injury (SCI) is reported to be in the range of
16 to 40 cases per million [1] depending upon country, with the causes of these injuries
PLOS ONE | DOI:10.1371/journal.pone.0169329 January 4, 2017 1 / 24
a1111111111
a1111111111
a1111111111
a1111111111
a1111111111
OPEN ACCESS
Citation: Cazzola D, Holsgrove TP, Preatoni E, Gill
HS, Trewartha G (2017) Cervical Spine Injuries: A
Whole-Body Musculoskeletal Model for the
Analysis of Spinal Loading. PLoS ONE 12(1):
e0169329. doi:10.1371/journal.pone.0169329
Editor: Karen Hind, Leeds Beckett University,
UNITED KINGDOM
Received: March 21, 2016
Accepted: December 13, 2016
Published: January 4, 2017
Copyright: © 2017 Cazzola et al. This is an open access article distributed under the terms of the
Creative Commons Attribution License, which
permits unrestricted use, distribution, and
reproduction in any medium, provided the original
author and source are credited.
Data Availability Statement: The musculoskeletal
models files and simulation data are available from
the FIGSHARE database [https://figshare.com/
account/home#/projects/17399]; RUGBY MODEL:
[https://dx.doi.org/10.6084/m9.figshare.4249835.
v1]; MASI model: [https://dx.doi.org/10.6084/m9.
figshare.4249808.v2]; RUGBY SCRUMMAGING
DATA: [https://dx.doi.org/10.6084/m9.figshare.
4249883.v2]. Also, all data are freely available from
Simtk.org [https://simtk.org/home/csibath].
Funding: This project is funded by the Rugby
Football Union (RFU) Injured Players Foundation.
ranging from motor vehicle collisions and community violence to recreational and work-
place-related activities [2]. Whilst traffic collisions are associated with 50% of SCIs, an
increasing proportion (currently 25%) of SCIs are due to sport and recreational activities [3].
This trend is mirrored by the demographic of SCI in the population, which is mainly repre-
sented by young men in their early thirties, who are more likely to be paraplegic, complete or
incomplete, as an outcome [4]. The societal costs for traumatic spinal cord injuries are high,
reported to be up to $9.7 billion ($US) per year in the USA [5], and up to £0.5–1 billion per
year in the UK [6].
Musculoskeletal modelling is widely used in the field of biomechanics for the analysis and
simulation of human motion, especially to investigate biomechanical variables that are not
directly measurable through in-vivo experiments. For example, musculoskeletal models have been used to: inform surgical decisions through simulation [7, 8], analyse joint load [9], and
identify motion patterns able to reduce injury risk in sporting activities [10]. Equally, musculo-
skeletal modelling can have many applications in the analysis of hazardous situations [11, 12]
such as vehicle crashes or catastrophic injuries in sport, where direct measures on internal ana-
tomical structures are not feasible or ethical.
In vitro experiments and in silico investigations have been widely performed in the field of transportation safety research [13, 14], but the primary mechanisms of acute spinal cord
injuries are yet to be fully elucidated. A key limitation in many studies is the difficulty in repli-
cating the contribution of neck and shoulder musculature to resisting external loads and,
therefore, translation of the results to real world conditions has been limited.
Rugby Union (rugby) is a full contact sport that on rare occasions can result in serious spi-
nal injuries, reported to range from 1-2/100,000 to 10/100,000 players per year [15, 16]. In the
last few years there has been a specific focus on improving the safety of the sport, for example
World Rugby (i.e. the international governing body) supported an injury prevention study
that led to international scrum law changes [17], with the ultimate aim to reduce catastrophic
injury occurrence in Rugby Union.
Similarly to transport applications, the biomechanical demands experienced during rugby
activities, such as scrummaging and tackling, have been widely analysed [18–20], but the
mechanisms of injury related to specific rugby activities are not yet fully elucidated. There is
still an open debate in the literature regarding the primary injury mechanisms for cervical
spine injury within rugby [21, 22], and both ‘buckling’ and ‘hyper flexion’ mechanisms have
been proposed. Kuster, Gibson (21) argued that buckling is the predominant mechanism,
whilst Dennison, Macri (22) stated that it is too early to draw that conclusion, and highlighted
the need for in silico analyses, together with in vivo and in vitro tests, in order to provide new insights and identify viable routes for injury prevention.
We identified rugby activities as an appropriate proof of concept for the analysis of cervical
spine loading and injury mechanisms using an in silico approach. During rugby scrummaging and tackling, the external load is applied through players’ shoulder girdle, feet and potentially
other body regions, therefore, a full body model including detailed modelling of cervical spine
and shoulder regions is required. Currently available ‘full-body models’ are limited to a com-
bined trunk-upper limb segment or include welded scapuloclavicular joints. Besides, the
majority of full body musculoskeletal models are used in gait analysis and locomotion simula-
tions [23, 24], whilst those specifically focused on upper limb [25, 26], lumbar spine [27–29],
and cervical spine [30] do not include either lower limbs or are kinematic-only models [31].
Therefore, the dynamic simulation of rugby activities necessitates a custom musculoskeletal
model, which is purposely optimised for the analysis and simulation of rugby injury scenarios.
Similar prerequisites are required also for the analysis of other injury scenarios, such as motor
vehicles collisions and human collisions in physical activity settings.
Rugby Musculoskeletal Model
PLOS ONE | DOI:10.1371/journal.pone.0169329 January 4, 2017 2 / 24
The funder had no role in study design, data
collection and analysis, decision to publish, or
preparation of the manuscript.
Competing Interests: The authors have declared
that no competing interests exist.
Musculoskeletal models are often linearly scaled for segment dimensions [32, 33], but this
procedure does not take into account subject- or population-specific musculoskeletal geome-
try and muscle parameters, which are fundamental to represent specific pathologies in clinical
applications [34] and to describe the morphology of specific sub-populations. For instance,
a rugby union Premiership forward player is on average 1.89 ± 0.07 m tall with a mass of 110.6 ± 7.7 kg [35], which is quite different from the population mean (height: 1.76 m, mass: 80.1 kg) [36] and more than 2 standard deviations greater than the anthropometrics used for
‘normal’ male subjects (height: 1.80 m, mass: 75.2 kg) [37] used in the creation of generic mus-
culoskeletal models [37]. Thus, in the case of a rugby player, the body segment inertial proper-
ties and muscle functional parameters can differ considerably from those included in generic
models.
The aim of the present study was to create, verify, and validate a generic full-body musculo-
skeletal model (‘Musculoskeletal model for the Analysis of Spinal Injuries’, or ‘MASI’) for in silico analyses of cervical spine loading during daily-living and sporting activities. A further population-specific version (the ‘Rugby Model’) of MASI that provides an in silico reproduc- tion of a rugby forward player was customised and validated for the investigation of injury
mechanisms during rugby activities.
2. Materials and Methods
The novel improvements embedded in MASI consist of i) a composite scapula-clavicular joint
that provides the linkage between cervical spine, upper limb and the remainder of the body,
and allows the application of multiple external loads, and ii) the inclusion of inertial and func-
tional parameters to permit dynamic analyses. The Rugby Model was created to test the MASI
in a specific sporting application, and differs from it only by the inclusion of inertial and mus-
cle parameters more representative of rugby forward players.
2.1 Model structure
The ‘Musculoskeletal model for the Analysis of Spinal Injury’ (MASI) and Rugby Model were
created in OpenSim (OpenSim 3.2, Simbios, Stanford, CA, USA) and Matlab software (Matlab
2013b, MathWorks, Natick, MA, USA), and share the same structure.
MASI inherited the structure of the OpenSim head and neck model [30] which we embed-
ded into a full body model (‘2354’, [23]), and was implemented to provide, for the first time,
the linkage between cervical spine, upper limb, torso and lower limbs (Fig 1). MASI comprises
35 rigid anatomical segments, 78 upper and lower cervical muscles divided into 19 muscle
groups, along with 23 torque actuators representing lower and upper limb muscles’ actions
(Fig 1). Motion between body segments was permitted via 34 joints and 30 kinematic con-
straints. To incorporate the effect of upper limb position, a new scapula-clavicular joint (SCJ)
(combining the joint motions of the acromioclavicular and sternoclavicular joints) was devel-
oped and included in the MASI, replacing the welded scapula-clavicular joint of the original
head and neck model. The model had 43 degrees of freedom, though these were reduced to 37
by locking the metatarsophalangeal and wrist joints into the neutral position.
The creation, verification and validation of the MASI and Rugby Model were split into sev-
eral phases to confirm firstly the performance of the scapula-clavicular joint, and secondly the
cervical spine joints’ moment-generating capacity due to active and passive muscle forces.
2.2 Scapula-Clavicular Joint (SCJ)
The new SCJ was created in order to reproduce shoulder girdle motion in the scapular plane
and include its effect on neck muscle behaviour while performing activities. Given the
Rugby Musculoskeletal Model
PLOS ONE | DOI:10.1371/journal.pone.0169329 January 4, 2017 3 / 24
magnitude of the external load applied in some rugby activities, the SCJ was designed to mini-
mise the number of intermediate massless bodies and kinematic constraints [38], which are
integrated in more complex models [25, 26]. This modelling choice was driven by the need to
minimise model complexity, as MASI and Rugby Model include a large number of kinematic
constraints to reproduce neck movements.
The kinematics of the shoulder girdle is a composite of the motions of three bones (humerus,
scapula, and clavicle) interacting at four articulations (glenohumeral, acromioclavicular,
scapulothoracic, and sternoclavicular joints). The novel SCJ merges acromionclavicular and
sternoclavicular joint motions, and its motion is a function of humeral elevation (Fig 2) as dem-
onstrated in cadaveric and in vivo studies [39, 40]. The sternoclavicular joint was designed fol- lowing the International Society of Biomechanics standard [41], whilst the acromionclavicular
joint used a glenoid-based reference system [42]. The glenoid plane was defined based on land-
marks on the glenoid rim, which were visually identified in OpenSim 3.2, and the glenoid-
based system reference was reconstructed using Matlab software.
The SCJ was integrated in OpenSim through two custom joints that link the torso with the
clavicle and the scapula with clavicle. Six mathematical functions (‘coordinate-coupled con-
straints’) were used to couple humeral elevation to clavicular and scapular motions. Permitted
Fig 1. The degrees of freedom (DoF), and the anatomical segments of the Rugby Model. Body
segments were divided into lower limbs (femur, tibia, talus, calcaneus and toe), pelvis, trunk, seven cervical
spine vertebrae (C1-C7), scapula, clavicle, upper limbs (humerus, radius, ulna, hand) and head (skull and
jaw). The 23 actuators of the model were associated to shoulder, elbow, spine, hip, knee, and ankle degrees
of freedom.
doi:10.1371/journal.pone.0169329.g001
Rugby Musculoskeletal Model
PLOS ONE | DOI:10.1371/journal.pone.0169329 January 4, 2017 4 / 24
scapula and clavicle motions in the coronal, sagittal and transverse planes were obtained from
in vivo data in the literature [40] and compared against in vivo kinematics measured using elec- tromagnetic tracking system and bone pins [40, 42]. The coupling functions for the clavicle
and scapula were extrapolated by correlating their motions relative to the humeral elevation,
so that the sternoclavicular joint was driven by humeral elevation and the acromioclavicular
joint by the sternoclavicular joint movements (Fig 2).
2.3 Inertial properties and neck muscles parameters
The MASI inherited body segment inertial properties mostly from the ‘2354’ model [23],
whilst the cervical spine and head segments inertial properties were derived from in vitro data [43].
Body segment inertial properties for the Rugby Model were obtained from the literature
[44–47] and combined with data from a DEXA scan (Hologic Discovery W, Bedford, MA;
QDR software version 12.4.2) of a front row rugby player (height 1.84 m, 120.4 kg). The partic-
ipant provided written informed consent prior to participation and ethical approval was
obtained from the University of Bath Institutional Ethics Committee. The DEXA imaging
allowed a more realistic distribution of the masses across the anatomical segments, while the
local centre of mass locations and moments of inertia were personalised by inputting anthro-
pometric data of the same rugby player into Yeadon’s mathematical inertia model [48]
(Table 1).
The maximum isometric force of the MASI neck muscles was initially checked by compar-
ing the neck muscles’ moment-generating capacity of the OpenSim head-neck model (e.g.
Vasavada model) [30] with measurements of healthy subjects’ neck strength from in vivo stud- ies [49, 50]. However, Vasavada’s muscle parameters were obtained from in vitro cadaveric studies [51], and the neck muscles’ moment-generating capacity was found to be lower than
healthy male subjects’ neck strength values [49, 50] (Table 2). This discrepancy was even more
evident when compared with rugby players’ neck strength values [52], and particularly for the
Fig 2. The SCJ’s reference systems on scapula (acromionclavicular) and clavicle (sternoclavicular)
joints and their coupling functions driven by humeral elevation. The sternoclavicular joint origin on the
sternal extremity of the clavicle and its reference system were designed to allow depression/elevation rotation
about the x-axis, protraction/retraction rotations about the y-axis, and axial rotation about the z-axis. The
motion of this joint is driven by the humeral elevation via 3 coordinate-coupler constraints based on spline
functions. The acromionclavicular joint was designed with a glenoid-based reference system which allows
lateral/medial rotation on x-axis, protraction/retraction rotations on y-axis, and anterior-posterior tilt on z-axis.
The glenoid based system had the z-axis perpendicular to the glenoid plane, the y-axis directed superiorly
toward the superior glenoid tubercle, and the x-axis directed anteriorly perpendicular to the other 2 axes. The
motion of the acromionclavicular joint is driven by the sternoclavicular joint motion via 3 coordinate-coupler
constraints based on linear functions.
doi:10.1371/journal.pone.0169329.g002
Rugby Musculoskeletal Model
PLOS ONE | DOI:10.1371/journal.pone.0169329 January 4, 2017 5 / 24
flexion moment generating-capacity, which was already identified as a limitation of the Vasa-
vada model due to modelling assumptions [30]. Thus, we decided to scale the maximum iso-
metric force of the Vasavada model in the MASI, aiming to match representative neck
strength values for healthy male subjects [49, 50] for extension and lateral bending (Table 2).
Two different scaling factors were used for extensors (i.e. 1.5) and flexors (i.e. 1.4); these scal-
ing factors were calculated as the ratio of healthy male subjects’ neck strength [49, 50] to the
maximum moment generated by the generic Vasavada model (Table 2). Similarly, the maxi-
mum isometric force of the Rugby Model’s muscles was scaled using different scale factors for
extensors (i.e. 1.9) and flexors (i.e. 2.7) in order to match extension and lateral bending neck
strength values of rugby players [52] from in vivo studies (Table 2). The modelling choice of specifically matching extension and lateral bending moment rather than flexion moment
avoided the use of very high scaling factors that would unrealistically increase the passive
moment of the flexors (i.e. passive stiffness).
2.4 MASI and Rugby Model verification and validation
Once the SCJ inertial properties and muscle parameters had been implemented, the models
were verified and validated following the standards detailed in the literature for musculoskele-
tal modelling research [38]. The verification and validation procedures consisted of i) a kine-
matic validation, ii) a dynamic verification, and iii) a dynamic validation. The validated MASI
and Rugby Model were then used to perform a forward simulation (i.e. computed muscle
Table 1. Mass distribution of the MASI derived from previous models and values from in vitro study [43]. Mass distribution of the Rugby Model calcu-
lated from DEXA values of a rugby forward player (1.84 m; 120.4 Kg). Masses are reported as percentage of total body mass. The principal moment of inertia
(IXX, IYY, IZZ) for the Rugby Model are shown in the las three columns of the table, and are expressed in kgm 2 .
SEGMENT MASI MASS (%) Rugby Model MASS (%) Rugby Model Ixx Rugby Model Iyy Rugby Model Izz
HEAD 5.0 4.1 0.03328 0.01250 0.03589
JAW 0.6 0.6 0.04000 0.02000 0.04000
TORSO 26.2 22.3 1.22530 0.96260 2.18700
CERV1 0.7 0.6 0.03424 0.13697 0.03424
CERV2 0.9 0.8 0.03064 0.12254 0.03064
CERV3 0.6 0.6 0.03335 0.13338 0.03335
CERV4 0.6 0.6 0.03307 0.13229 0.03307
CERV5 0.6 0.6 0.03263 0.13051 0.03263
CERV6 0.7 0.7 0.03151 0.12605 0.03151
CERV7 0.9 0.8 0.03082 0.12327 0.03082
PELVIS 14.7 18.2 0.19204 0.16271 0.10816
SCAPULA 3.0 3.2 0.00124 0.00115 0.00137
CLAV 0.6 0.7 0.00024 0.00026 0.00004
FEMUR 9.3 9.9 0.17275 0.04528 0.18217
TIBIA 4.6 4.0 0.04650 0.00820 0.04650
TALLUS 0.1 0.1 0.00100 0.00100 0.00100
CALCANEUS 1.6 1.3 0.00179 0.00499 0.00525
TOES 0.3 0.2 0.00014 0.00028 0.00139
HUMERUS 2.5 3.4 0.02410 0.00831 0.02705
ULNA 0.8 1.0 0.00585 0.00122 0.00635
RADIUS 0.8 1.0 0.00585 0.00122 0.00635
HAND 0.6 0.8 0.00175 0.00108 0.00264
doi:10.1371/journal.pone.0169329.t001
Rugby Musculoskeletal Model
PLOS ONE | DOI:10.1371/journal.pone.0169329 January 4, 2017 6 / 24
control) driven by in vivo neck functional movement kinematic data in order to simulate acti- vation of the neck muscles and compare with experimentally-obtained muscle activity data.
Finally, the Rugby Model was used to run scaling, inverse kinematics, and inverse dynamic
simulations during a rugby union scrum against an instrumented scrum machine, in order to
test its application to rugby contact events.
2.4.1 Kinematic validation. The kinematic validation was based on the comparison of
model kinematics with experimental kinematics, and included i) the calculation of the range
of motion of the SCJ, and ii) the analysis of moment arms of the muscles that originate on the
scapula and clavicle bodies throughout the humeral range of motion.
The SCJ and humeral range of motion were calculated through an inverse kinematics pro-
cedure in OpenSim, driven by experimental in vivo data of a rugby forward player (age: 22 years, height: 1.77 m, mass: 88 kg) performing a humeral elevation in the scapular plane. The
SCJ motion was compared with in vivo data of scapula and clavicle motion measured using bone pins [42], and in silico kinematics data generated using Holzbaur et al.’s previously vali- dated OpenSim ‘Upper-Limb’ model [26]. In silico kinematics data were generated through an inverse kinematics simulation driven by the same experimental in vivo kinematic data of humeral elevation.
The inverse kinematics procedures enabled the calculation of the moment arms for the ster-
nocleidomastoid and trapezius muscle groups using the Effective Torque Method [53], with an
output comparison between the MASI and Holzbaur et al.’s model [26].
2.4.2 Dynamic verification. The verification process aimed to verify that the custom SCJ
adheres to the laws of physics, and in the case of the custom joint verification process, the
main goal was to verify that the kinematic constraints did not generate any extra work, as con-
straint forces are often large but they are applied in non-movement directions so the power is
zero. The power generated by the constraint was calculated as the sum of two elements: i) the
scalar product of the external forces applied to the scapula and of the velocity of the centre of
Table 2. Neck strength values are reported from in vivo and in silico studies. Extension (Ext), flexion (Flex), axial rotation (Axial Rot), and lateral bending
(Lat Bend) scaling factors are shown for adult rugby players (Rugby Front Row) and adult healthy males (Healthy Male). The scaling factors used to scale the
MASI and Rugby Model neck muscles are calculated independently for flexors and extensors in order to match maximum extension and lateral bending neck
strength values for a healthy male (Healthy Male) and a front row rugby player (Rugby Front Row), respectively. Head-Neck Model [30] was used as scaling
reference, thus its scaling factor is 1.
Source Ext—Flex (Nm) Axial Rot (Nm) Lat Bend (Nm) Scaling Factors
Head-Neck Model
Vasavada et al. 1998 34.1–3.6 10.6 22.6 1
(in silico data)
Rugby Front Row Ext (1.9)
Olivier et al. 2008 66–43 - 61 Flex (11.8)
(in vivo data) Lat Bend (2.7)
Healthy Male Ext (1.5)
Fice et al. 2014 51–30 13 32 Flex (8.3)
(in vivo data) Axial Rot (1.2)
Lat Bend (1.4)
Rugby Front Row Player
Rugby Model 72.3–15.9 19.7 58.5 Ext (1.9)
(in silico data) Flex (2.7)
Healthy Male
MASI 50.8–10.3 12.4 31.3 Ext (1.5)
(in silico data) Flex (1.4)
doi:10.1371/journal.pone.0169329.t002
Rugby Musculoskeletal Model
PLOS ONE | DOI:10.1371/journal.pone.0169329 January 4, 2017 7 / 24
mass of the scapula; and, ii) the scalar product of the scapula generalised forces and corre-
sponding generalised speed.
2.4.3 Dynamic validation. The dynamic validation consisted of the determination of i)
the passive neck stiffness, and, ii) the maximum net joint moment generated by the muscles
throughout neck functional movements. The effect of humeral elevation on passive neck stiff-
ness was also assessed.
The maximum net joint moment was estimated to occur when agonist muscles have maxi-
mal activation and antagonists have zero activation. The MASI and Rugby Model adopted a
previously validated muscle-tendon model [54]. The passive neck stiffness was calculated
when muscles had zero activation, and the neck motion generated under specific loading con-
ditions was used to compare MASI and Rugby Model passive stiffness with in vitro and in vivo values [55]. Forward dynamics simulations generated neck movements (e.g. flexion-extension,
axial rotation and lateral bending) through torque actuators that replicated the ‘2:4:2’ loading
protocol (i.e. 2 Nm for flexion-extension and lateral bending, and 4 Nm for axial rotation
applied to C2 vertebra) used in Miura, Panjabi (55). All the other bodies, except the scapula,
the clavicle and the upper limbs, were locked in a stationary position The neck motion gener-
ated by the MASI and Rugby Model were considered acceptable if within 2 SD from in vitro values [55].
Maximal net joint moments were calculated during forward dynamics procedures that sim-
ulated a maximum isometric contraction with all the model’s bodies were locked and not free
to move. The MASI’s and Rugby Model’s maximal net joint moments were considered accept-
able if within 2 SD of in vivo reference data for maximal isometric neck contractions on dyna- mometers [49, 56].
2.4.4 Neck muscles simulated activation during functional movement. Two different
participants were recruited for the validation of muscle activation during neck functional
movements using either the MASI or Rugby Model. Experimental data of full body kinematics,
ground reaction forces, and neck muscles’ EMG of a healthy male subject (age: 64 years,
height: 1.67 m, mass: 75 kg) and a male rugby forward player (age: 22 years, height: 1.77 m,
mass: 88 kg) were collected during neck flexion, extension, lateral bending and axial rotation
movements, respectively. The participants provided written informed consent prior to partici-
pation and ethical approval was obtained from the ‘Research Ethics Approval Committee for
Health’ (REACH) of the University of Bath.
The trajectories of 68 reflective markers were measured through a 16-camera optoelectronic
system (Oqus, Qualisys, Sweden; sampling frequency: 100 Hz), ground reaction forces were
recorded via two Kistler force platforms (9287BA, Kistler Instruments Ltd, Switzerland; sam-
pling frequency: 1000 Hz), and EMG signals from bilateral sternocleidomastoid, upper trape-
zius fibers and mid trapezius fibers were collected using Delsys Trigno (DelsysInc, Boston,
Massachusetts, USA; sampling frequency: 2000 Hz). Raw EMG signals were rectified, filtered
(Butterworth order 2, bandpass 10–500 Hz), and the envelope was created applying a 50-ms
moving rectangular window. Experimental EMG signals were normalized using maximum
voluntary contraction data collected using the procedure defined in Cazzola, Stone [57].
The MASI and the Rugby Model were scaled within OpenSim procedures to match the spe-
cific participant’s anthropometry based on the 68 motion analysis markers placed on anatomi-
cal landmarks (i.e. S1 Appendix). The OpenSim inverse kinematics algorithm solved for the
joint angle time histories and then joint moments associated with participants’ motion were
calculated using an inverse dynamics algorithm. A muscle control algorithm (computed mus-
cle control or CMC) [58, 59] was then used to compute the muscle excitations required to
track the kinematics produced by the inverse kinematics algorithm, by minimizing the sum of
the square of muscle activations, while accounting for muscle activation and contraction
Rugby Musculoskeletal Model
PLOS ONE | DOI:10.1371/journal.pone.0169329 January 4, 2017 8 / 24
dynamics [60]. The output of the CMC simulation was compared against the measured EMG
data from the sternocleidomastoid and trapezius muscles during functional neck movements.
2.4.5 Rugby Model application to rugby contact events: machine scrummaging simula-
tion. A front row rugby player (age: 22 years, height 1.84 m, mass 120.4 kg) performed indi-
vidual scrummaging trials against a strain-gauge instrumented scrum machine collecting
external ‘shoulder’ forces at 500 Hz with feet positioned on force plates (9287BA, Kistler
Instruments Ltd, Switzerland; sampling frequency: 2500 Hz). During scrum trials, a total of 68
reflective markers were positioned on the player (i.e. S1 Appendix) and tracked by a motion
analysis system (Oqus, Qualisys, Sweden; sampling frequency: 250 Hz) and all data streams
(kinematics, forces and EMG) were time-synchronised. The participant provided written
informed consent prior to participation and ethical approval was obtained from the ‘Research
Ethics Approval Committee for Health’ (REACH) of the University of Bath.
Joint kinematics (i.e. joint angles) and kinetics (i.e. joint net moment and joint reaction
forces) at each individual cervical vertebra were obtained by using the Rugby Model during
‘Pre-Engagement’ and ‘Engagement’ phases of machine scrummaging. The Rugby Model was
scaled to the participant’s anthropometric data, and inverse kinematics, Residual Reduction
Algorithm (RRA), joint reaction forces analyses were run in OpenSim 3.2 (i.e. S1 Video). The
RRA [61] was run to resolve the dynamic inconsistency between the measured kinematics and
external loads applied to the system. The RRA routine was processed with force and kinematic
data filtered through a 3rd-order low-pass bidirectional Butterworth filter at 12 Hz.
The joint reaction forces analysis was run without considering neck muscles contribution,
and therefore calculating the joint reaction forces as follows:
Rj ¼ mjaj � ð X
Fext þ X
Fconst þ Rjþ1Þ
where Rj is the reaction force at joint ‘j’, mj is the mass of segment ‘j’, aj is the centre of mass
acceleration of a segment ‘j’, Fext is the sum of all external force applied to a segment ‘j’, Fconst
is the sum of all constraint forces applied to a segment ‘j’ and Rj+1 is the reaction force at an
adjacent segment ‘j+1’.
3. Results
3.1 Kinematic validation: Scapula-Clavicular Joint (SCJ) and moment
arms
The scapula and clavicle movements throughout the humeral range of motion were within two
standard deviations of experimental measurements from in vivo [40, 42] studies (Fig 3). The sternoclavicular (SC) and acromioclavicular (AC) motions modelled through the SCJ were
similar to experimental kinematics measured in the Ludewig, Phadke (40) study over the 20–
140˚ range of humeral elevation: the RMS deviation of SC motions varied from 0.1˚ to 0.3˚,
whilst the RMS deviation of AC motions varied from 2.5˚ to 3.6˚.
The AC motions generated by the SCJ were also compared with a previous in silico study [26]: the RMS deviation from Holzbaur’s model across the range of humeral elevation was
higher for scapula lateral rotation (12.4˚), than protraction-retraction (1.2˚) and anterior-pos-
terior tilt (3.3˚) motions. The higher deviation for lateral rotation was due to a different behav-
iour between the SCJ and Holzbaur’s model for humeral elevation values greater than 90˚ (Fig
3). However, the same AC motions compared with in vivo measurements [42] showed a low RMS deviation across lateral rotation (3.4˚), protraction-retraction (1.8˚), and anterior-poste-
rior tilt (3.5˚).
Rugby Musculoskeletal Model
PLOS ONE | DOI:10.1371/journal.pone.0169329 January 4, 2017 9 / 24
The moment arms of the muscles having origin or insertion points on the clavicle and scap-
ula bones during humeral elevation in the scapular plane (Fig 4) equalled the previous Vasa-
vada model estimations when the neck and the humerus were in neutral position (Fig 4a and
4b), and were comparable with Holzbaur’s model (Fig 4a and 4b) across the range of neck
motion and humeral elevation.
The MASI model had muscle moment arms which varied over humeral elevation, which is
key for correctly estimating the contribution of neck muscle forces during activities involving
upper limb motion. The sternocleidomastoid muscle group mainly showed a decrease (~0.01
m) of flexion moment arm in neutral and fully flexed head poses through the range of humeral
elevation, and an almost negligible decrease (~0.005 m) in the axial rotation moment arms
when the neck was axially rotated (Fig 4c). The flexion moment arms of the superior fibres of
trapezius (i.e. Trap Clav in Fig 4c) mainly increased (~0.020 m) during humeral elevation
when the neck was fully flexed and in neutral position (~0.017 m), whilst an extended position
minimised this change. Levator scapulae lateral bending moment arms increased during
humeral elevation when the neck was laterally bent (~ 0.01 m) and in neutral pose (< 0.01 m),
whilst its extension moment increased in fully extended pose (Fig 4c).
Fig 3. Sternoclavicular (SC) (left column) and acromioclavicular (AC) (middle and right columns) motions during humeral elevation. The
black diamonds in the left and middle columns graphs represent SC and AC kinematics during in vivo measurements [40], whereas black solid lines
are the respective SC and AC kinematics generated by the SCJ. The SC and AC joints angles (black solid line) are within 2SD from in vivo values [40]
(black dotted lines). In the right column graphs, the AC kinematics generated by SCJ is compared with another in vivo study and an in silico study. The
AC motions generated by the SCJ (black solid line) is within 2SD (black dotted lines) from in vivo studies in the literature [41] (black squares) and
comparable to the output of the Holzbaur’s model [26] (black dash-dot line). The glenoid reference system of the scapula was roto-traslated in order to
express scapula motion with respect to the acromioclavicular joint reference system, and compare it with Ludewig, Hassett (42) and Holzbaur, Murray
(26) studies.
doi:10.1371/journal.pone.0169329.g003
Rugby Musculoskeletal Model
PLOS ONE | DOI:10.1371/journal.pone.0169329 January 4, 2017 10 / 24
3.2 Dynamic verification and validation
The forward simulations confirmed that the power generated by the SCJ kinematic constraint
was nil, therefore the constraints were workless. The passive joint moments generated during
neck movement and humeral negative elevation (Fig 5a) were comparable with in vitro data [55], in the region of 2–4 Nm.
The maximal net joint moments generated by the MASI were comparable with healthy
male subjects’ neck strength for extension and lateral bending movement, whilst flexion
moment-generating capacity of the model resulted more realistic than Vasavada’s model, yet
still 3 times lower than in vivo values previously reported (Table 2). For the Rugby model, the
maximal net joint moments generated by activated neck muscles during extension and lateral
bending were comparable with rugby forward players’ neck strength in vivo values, whilst flex-
ion moment-generating capacity was 3 times lower (Table 2).
The passive neck stiffness of MASI was generally within 2SD from Miura, Panjabi (55) in vitro values and the average of in vivo values [62–67] in flexion-extension, lateral bending and axial rotation motions (Fig 5). During flexion-extension, the lower-middle cervical spine (C7/
T1-C4/C5) showed a higher level of agreement with in vivo values than in vitro values, as joint
Fig 4. Moment arms of elevator scapulae (A) and sternocleidomastoid (B) muscles throughout humeral elevation when the neck was fully extended (Ext
—black dashed and dotted line), flexed (Flex—black dotted line) and in neutral position (Neutral—black solid line). Vasavada’s (grey circles) and
Holzbaur’s moment arms (red solid line) across humeral elevation are compared with MASI’s moment arms when the neck is at neutral position. (C)
Muscles’ moment arms maximum changes of flexion (FLEX), extension (EXT), lateral bending (LAT BEND) and axial rotation (AXIAL ROT) moments
during humeral elevation when the neck is fully flexed, extended, laterally bent, axially rotated and in neutral position (NEUTRAL). Positive values
represent an increase whilst negative values a decrease of the moment arm across humeral elevation.
doi:10.1371/journal.pone.0169329.g004
Rugby Musculoskeletal Model
PLOS ONE | DOI:10.1371/journal.pone.0169329 January 4, 2017 11 / 24
motion was slightly higher (5 to 10%) than in vitro values. On the contrary C2/C3 joint showed higher stiffness and lower flexion-extension motion than both in vivo and in vitro values. Pas- sive neck stiffness in lateral bending was generally similar to in vitro values for middle-upper spine (C5-C2), except C7/T1 and C6/C7 joints that showed values greater than 2SD different.
Neck passive stiffness in axial rotation was generally lower (i.e. higher degree of motion) than
in vitro values, and more in line with in vivo values. The neck passive stiffness of the Rugby Model followed the same pattern shown by the
MASI, but the degree of motion across all the neck movements was lower than MASI (~8% in
flexion-extension and ~11% in axial rotation), in particular for lateral bending (~30%) due to
higher stiffness emanating from the doubling of maximal isometric muscle forces. Concomi-
tantly, in vivo values for rugby forward player’s cervical spine range of motion were found lower than for healthy males [68].
Fig 5. Neck joints motion under a 2 Nm (flexion-extension and lateral bending) and 4Nm (axial
rotation) load applied at C2 vertebra. MASI’s neck motion (black bars) and Rugby Model’s neck motion
(black downward diagonal bars) are compared with an average value from different in vivo studies (grey bars)
and two protocols of an in vitro study: ‘2-4-2’ protocol (white dotted bars) and ‘1:1:1’ protocol (white bars).
Average in vivo values were not available for lateral bending individual joint motion (middle graph).
doi:10.1371/journal.pone.0169329.g005
Rugby Musculoskeletal Model
PLOS ONE | DOI:10.1371/journal.pone.0169329 January 4, 2017 12 / 24
The passive neck moment in flexion and lateral bending motion, when the humerus was
maximally elevated, increased respectively to 5 Nm and 4 Nm due to the increased length of
the levator scapulae and trapezius muscles. Passive neck moment in axial rotation at maximal
humeral elevation was lower than 2 Nm.
3.3 Neck muscles activation during functional movements
Using qualitative inspection, the simulated activation of trapezius muscles showed a similar
pattern and activation level in comparison with the recorded EMGs across the neck move-
ments analysed (Figs 6 and 7).
The sternocleidomastoid muscles activations simulated by using the Rugby Model were
generally comparable with experimental EMGs, except for the left sternocleidomastoid during
extension and lateral bending (Fig 6). During neck extension the left sternocleidomastoid was
activated at the beginning of the movement, although experimental EMG showed an almost
nil activation across the movement (Fig 6). However, the sternocleidomastoid activations sim-
ulated by using the MASI showed more consistent pattern compared with the experimental
EMG (Fig 7).
During lateral bending, the activation of the left sternocleidomastoid simulated by using
the Rugby Model showed an anticipated high activation in left bending, and a late high activa-
tion during right bending with respect to experimental EMG (Fig 6). However, the same acti-
vation pattern was not present when sternocleidomastoid activation was simulated by using
MASI (Fig 7).
3.4 Biomechanical load experienced during machine scrummaging
The cervical spine was in an extended position at T1-C7 to C2-C3 joints prior to contact but
then underwent a flexion motion during impact, whereas the C2-HEAD joint remained
extended throughout the engagement phase (Fig 8b). This cervical spine motion caused a flat-
tening of cervical spine curvature during the engagement phase.
The biomechanical load experienced at vertebral joint level is representative of the joint
moments and forces that the neck muscles and other passive structures have to overcome or
resist during a machine scrummaging trial (Fig 8a). The joint net moments calculated across
all the vertebral joints were quite low and representative of the moment needed to support the
head and cervical vertebra weight (Fig 8c).
The peaks of joint net forces at cervical spine level were higher in longitudinal direction (Y
axis) than in mediolateral (Z axis) and anteroposterior (X axis) directions. Longitudinal joint
reaction forces decreased from lower to upper cervical joints (Fig 9) ranging from ~-500 N
(T1-C7) to ~-200 N (C1-HEAD). The right and left acromioclavicular joints showed the high-
est joint reaction forces in both longitudinal (~1700 N) and mediolateral (~600 N) directions
(Fig 9). The cervical spine joint reaction forces in mediolateral directions ranged from ~80 N
(T1-C7) to ~250 N (C1-HEAD). Overall, the joint reaction analysis in longitudinal direction
showed a ‘distraction’ pattern of force (i.e. tensile load) at cervical spine level. This is due to the
fact that during the engagement phase the external load is mainly acting at the shoulders with
the effect that the applied force is pushing back the trunk whilst the neck and head bodies are
accelerated forward by the inertia of the preceding movement.
4. Discussion
The MASI and the Rugby Model are novel musculoskeletal models created for injury mecha-
nism investigations, and its final aim is to provide internal measurements for cervical spine
structures that cannot be collected experimentally because they are too invasive or not ethically
Rugby Musculoskeletal Model
PLOS ONE | DOI:10.1371/journal.pone.0169329 January 4, 2017 13 / 24
Fig 6. Rugby Model simulated neck muscle activation. Simulated muscle activation from computed muscle
control (solid black line) and experimental EMG (dashed black line) of sternocleidomastoid (left column) and
upper trapezius (right column) muscles during flexion, extension, lateral bending (right and left bending) and
axial rotation (right and left rotation). Experimental EMG signal were normalized using maximum voluntary
contraction data and defined between 0% and 100%. Simulated activations are defined between 0% (no
activation) and 100% (full activation).
doi:10.1371/journal.pone.0169329.g006
Rugby Musculoskeletal Model
PLOS ONE | DOI:10.1371/journal.pone.0169329 January 4, 2017 14 / 24
Fig 7. MASI simulated muscle activation. Simulated muscle activation from computed muscle control (solid black line) and experimental average
rectified EMG (dashed black line) of sternocleidomastoid (left column) and upper trapezius (right column) muscles during flexion, extension, lateral
bending (right bending) and axial rotation (right rotation). Experimental EMG signal were normalized using maximum voluntary contraction data and
defined between 0% and 100%. Simulated activations are defined between 0% (no activation) and 100% (full activation).
doi:10.1371/journal.pone.0169329.g007
Rugby Musculoskeletal Model
PLOS ONE | DOI:10.1371/journal.pone.0169329 January 4, 2017 15 / 24
acceptable. The MASI is an in silico representation of a healthy male subject and the present validation demonstrates that it is suitable for investigating cervical spine loading through both
inverse and forward simulations. The key novel aspects of the presented models with respect
to previous head and neck models [30] are i) the integration of the head/neck model with a
comprehensive full body model, ii) the implementation of the scapula-clavicular joint (SCJ)
for the application of external loads on the shoulder girdle, and iii) the inclusion of popula-
tion-specific inertial properties (e.g. Rugby Model) and muscle parameters to permit dynamic
analyses.
The analysis of cervical spinal injuries mechanisms is one of the grand challenges of injury
prevention research, and we identified rugby union activities as a perfect proof of concept due
to impact loading being applied to the shoulder/neck region and distinctive anthropometric
properties of the participants. Thus, we created and validated the Rugby Model, which is a
population-specific version of the MASI, and includes a set of inertial properties and muscles
parameters to reproduce the features of a rugby forward player. The Rugby Model was used to
analyse for the first time the cervical spine loading in rugby-related contact events (i.e. inverse
Fig 8. (A) The three main phases of a rugby scrummaging activity: pre-engagement phase, engagement phase, and sustained-push. (B)
Average flexion-extension joint angle across cervical spine joints during ‘Pre-Eng’ and ‘Engagement’ phases. Extension motion is positive in
the graph. (C) Flexion-extension joint moment values during pre-engagement and engagement phases across the cervical spine vertebral
joints. Extensor moment is positive in the graph. The graph shows the decreasing pattern of extensor moment from C7-C6 joint to C1-Head,
which is mainly due to resist the flexion moment generate by the gravity force.
doi:10.1371/journal.pone.0169329.g008
Rugby Musculoskeletal Model
PLOS ONE | DOI:10.1371/journal.pone.0169329 January 4, 2017 16 / 24
Fig 9. Typical joint reaction force traces of scapula, clavicle, thorax (thoracic and lumbar spine plus ribcage) and cervical
spine bodies during a machine trial. The graph at the top shows the joint reaction force in anterior-posterior direction (X axis), whilst
the central and bottom graphs show respectively the joint reaction forces in the longitudinal (Y axis) and mediolateral (Z axis)
directions. The time of real engagement (tENG) is highlighted by the dashed vertical line. The red arrow highlights the increasing
pattern of tensile load at cervical spine level.
doi:10.1371/journal.pone.0169329.g009
Rugby Musculoskeletal Model
PLOS ONE | DOI:10.1371/journal.pone.0169329 January 4, 2017 17 / 24
simulation of machine scrummaging) (i.e. S1 Video). The inclusion of lower limbs and shoul-
der girdle allowed to application of the model to situations that require the inclusion of multi-
ple external loads, such as rugby scrummaging.
The MASI and Rugby Model provide validated shoulder girdle and cervical spine kinemat-
ics, and allow estimates of the moment arms during coupled neck and arm movements. From
a dynamic perspective, forward simulations confirmed their validity at reproducing in vivo maximal net joint moments and comparable in vitro [55] cervical spine passive stiffness.
The custom scapula-clavicular joint (SCJ) was created for conducting kinetic analysis of
rugby activities during which the humeral and scapula-clavicular movements and external
loading through these structures play an important role. The SCJ successfully simulates sterno-
clavicular and acromioclavicular motions in comparison with experimental values from in vivo [40, 42] and in silico [26] studies (Fig 3). This comparison also showed that SCJ and Holz- baur’s models, at humeral elevation values larger than 90˚, generates respectively lower and
higher scapular lateral rotation values with respect to in vivo measurements. However, the SCJ deviation is within 2SD from in vivo data, whilst Holzbaur’s model exceeds 2SD for humeral elevation higher than 120˚. This is due to the different regression equations and kinematic
constraints used to drive the models, which can affect trapezius and levator scapulae moment
arms at large humeral elevation.
The SCJ allows the estimation of neck muscles moment arms during combined shoulder
and neck movements (Fig 4). This improves the moment-generating capacity of the model
[69], and therefore the estimation of individual neck muscle forces for the muscles that have
origin or insertion points on the clavicle and scapula bones. During a rugby scrum or rugby
tackle a player is exposed to external loads of 2–4 kN applied on his/her shoulders [18, 19, 70],
and neck movements are typically performed while the humerus is elevated in the scapular or
coronal planes. In these postures, moment arms are dependent on humeral elevation, and
their values can affect static optimisation algorithms. A further refinement would be the imple-
mentation of an algorithm for representing the 3D interaction of muscle surfaces, allowing
moment arms to be better estimated using wrapping surfaces based on subject-specific data
[69, 71], as neck muscle geometry of rugby players might considerably differ with respect to
healthy subjects due to hypertrophy as a result of muscle conditioning. In its current form the
novel SCJ is not intended to replace shoulder models based on closed-loop kinematic chains
[25, 72–74], and optimised to provide accurate simulations for clinical applications [75].
Indeed, as in previous shoulder models [26, 76], the SCJ is driven by regression equations that
do not change according to participants’ morphological characteristics, and does not include a
constraint reproducing the gliding plane of the scapulothoracic joint. The SCJ is only intended
to reproduce shoulder kinematics during humeral elevation in the scapular plane, which is the
most representative elevation plane for the analysis of rugby activities.
The MASI and Rugby Model allowed satisfactory replication of cervical spine passive and
active behaviour in healthy subjects and rugby players. However, the passive neck stiffness cal-
culated in this study is mainly representative of passive muscle forces, as the neck kinematic
constraints can only partially replicate the stiffness from other passive structures. This is a limi-
tation of MASI and Rugby Model and future inclusion of ligaments and intervertebral discs
will be beneficial for better reproducing stiffness values in lateral bending and axial rotation.
The neck muscle activations simulated by using the Rugby Model and MASI were generally
comparable with the experimental EMGs during neck functional movements (Figs 6 and 7),
although the sternocleidomastoid showed different activation patterns when not used as pri-
mary movers of neck and head bodies (see left sternocleidomastoid during extension and right
lateral bending in Fig 6). In these movements, the sternocleidomastoid is mainly used to decel-
erate neck extension via eccentric contraction and control head motion, which is a difficult
Rugby Musculoskeletal Model
PLOS ONE | DOI:10.1371/journal.pone.0169329 January 4, 2017 18 / 24
task to simulate with current optimisation functions. In fact, muscle activations are dependent
on the objective function used to solve the redundancy problem, and an optimised function
for neck motion analysis still needs to be found. In this study, we minimised the sum of the
neck muscle forces, but different objective functions that includes the use of synergies, or
EMG driven simulations should be explored to provide more realistic results. Therefore, our
simulation results demonstrate model suitability for the calculation of muscle forces during
functional activities rather than providing benchmarking values. However, the same validation
procedure provided positive results using the MASI driven by in vivo data of a healthy male subject (Fig 7), demonstrating that the model’s output is appropriate even with different iner-
tial parameters and muscle scaling factors. Further validation procedures are still needed for
reconstructing muscle activations during specific rugby activities such as scrummaging, in
which the combination of small upper body movements and muscle co-contractions [57] may
affect the effectiveness of the current static optimization redistribution approach based on the
minimisation of the sum of muscles activation squared or sum of muscles force.
The new set of inertial parameters of the Rugby Model derives from a rugby front row
player’s DEXA scan, with segment centre of mass position and moments of inertia personal-
ised to the individual using a mathematical inertia model [48]. This combined approach is
beneficial for generating a musculoskeletal model that better replicates rugby players’
anthropometry when medical imaging is not readily available, and therefore improves lin-
ear-scaling and inverse dynamics results with respect to a generic OpenSim model (e.g.
healthy male of 1.8 m and 75.16 kg). The proposed method does not substitute for a subject-
specific scaling method based on medical images (e.g. MRI and CT scans), which is still con-
sidered the gold standard to estimate muscle origin and insertion points [24] and body seg-
ment inertial properties [77]. Also, optimisation technique methods [78] or a scaling method
based on the combination of 3D body-surface and DEXA scans would provide a better esti-
mation of body geometry [79], a specific density profile [80] and would theoretically produce
more reliable simulation outputs.
The Rugby Model includes population-specific muscle strength parameters derived from in vivo studies, which becomes fundamental to generate realistic simulations when muscle physi- ological cross-sectional area (PCSA) and maximum isometric force values differ considerably
from generic data. Individualisation of muscle parameters can be performed using simple scal-
ing methods [81] or more complex imaging-based methods [34]. Bolsterlee, Vardy [82] dem-
onstrated that a uniform scaling routine provided reliable prediction of a shoulder model’s
strength, and the scale factors can either be obtained from muscle volume data or maximum
strength measurements. Given the modelling assumptions of the Vasavada model (i.e. poten-
tially different joint centre of rotation during flexion and lack of infrahyoid muscles in the
model) [30], the maximum isometric force of the Rugby Model’s neck muscles was scaled
using two scaling factors, both obtained from maximum strength measurements of rugby for-
ward players (Table 2).
The same approach was used for the MASI and allowed a better match to in vivo neck strength values of a healthy male subject with respect to previous models [30]. However, the
neck flexion moment generating-capacity was still three times lower that in vivo values, and this was a limitation of the Vasavada model due to modelling assumptions [30] that was inher-
ited by both the Rugby Model and MASI.
The construction of subject-specific models is still needed for more accurate simulations,
but the presented models can already provide the biomechanics community with a computa-
tional tool to analyse rugby activities (i.e. the Rugby Model) or more generic applications
where cervical spine loading is of interest (i.e. the MASI). Both models are freely downloadable
from SimTK (https://simtk.org/home/csibath).
Rugby Musculoskeletal Model
PLOS ONE | DOI:10.1371/journal.pone.0169329 January 4, 2017 19 / 24
Further modelling developments will include the replacement of neck kinematic con-
straints with contact models [28, 83] and the inclusion of EMG driven simulations [84] for the
estimation of muscle forces. Future application of the models will firstly include the simulation
of injurious events, and the analysis of neck muscle contribution to cervical spine loading dur-
ing rugby activities, such as scrummaging and tackling. This first step will allow us to further
optimise the model for analysis of collisions in sports such as American Football and Ice
Hockey, as well as the analysis of cervical spine injuries during motor vehicle collisions.
Supporting Information
S1 Appendix. Experimental Set-up. Description of the experimental set up and marker set
used in the study.
(DOCX)
S1 Video. Machine Scrummaging inverse simulation. Simulation of a machine scrummaging
trial from three different views. Green arrows represent the force applied to the rugby player
during the different phases of a machine scrummaging trial.
(M4V)
Acknowledgments
The authors would like to thank the players who took part in the experimental sessions and
the laboratory technician, Mr Andreas Wallbaum, for the technical input provided. The
authors would like to thank also Dr Elena Seminati, Dr Sabina Gheduzzi, Dr Pooya Mah-
moodi, and Prof Tony Miles for their involvement in the wider research programme on the
biomechanics of rugby injuries.
Author Contributions
Conceptualization: DC GT EP TPH HSG.
Data curation: DC EP TPH.
Formal analysis: DC.
Funding acquisition: GT DC EP.
Investigation: DC TPH.
Methodology: DC GT HSG.
Project administration: GT.
Resources: DC.
Software: DC.
Supervision: GT HSG.
Validation: DC.
Visualization: DC.
Writing – original draft: DC.
Writing – review & editing: EP GT TPH HSG.
Rugby Musculoskeletal Model
PLOS ONE | DOI:10.1371/journal.pone.0169329 January 4, 2017 20 / 24
References 1. Lee BB, Cripps RA, Fitzharris M, Wing PC. The global map for traumatic spinal cord injury epidemiol-
ogy: update 2011, global incidence rate. Spinal Cord. 2014; 52(2):110–6. doi: 10.1038/sc.2012.158
PMID: 23439068
2. Tator CH. Update on the pathophysiology and pathology of acute spinal cord injury. Brain pathology.
1995; 5(4):407–13. PMID: 8974623
3. Sekhon LHS, Fehlings MG. Epidemiology, demongraphics and pathophysiology of acute spinal cord
injury. Spine. 2001; 26(24S):S2–S12.
4. Wyndaele M, Wyndaele JJ. Incidence, prevalence and epidemiology of spinal cord injury: what learns a
worldwide literature survey? Spinal Cord. 2006; 44(9):523–9. doi: 10.1038/sj.sc.3101893 PMID:
16389270
5. Ackery A, Tator C, Krassioukov A. A global perspective on spinal cord injury epidemiology. Journal of
neurotrauma. 2004; 21(10):1355–70. doi: 10.1089/neu.2004.21.1355 PMID: 15672627
6. Cripps RA, Lee BB, Wing P, Weerts E, Mackay J, Brown D. A global map for traumatic spinal cord injury
epidemiology: towards a living data repository for injury prevention. Spinal Cord. 2011; 49(4):493–501.
doi: 10.1038/sc.2010.146 PMID: 21102572
7. Steele KM, Seth A, Hicks JL, Schwartz MS, Delp SL. Muscle contributions to support and progression
during single-limb stance in crouch gait. J Biomech. 2010; 43(11):2099–105. doi: 10.1016/j.jbiomech.
2010.04.003 PMID: 20493489
8. Piazza SJ, Delp SL. Three-dimensional dynamic simulation of total knee replacement motion during a
step-up task. Journal of biomechanical engineering. 2001; 123(6):599–606. PMID: 11783731
9. McLean SG, Su A, van den Bogert AJ. Development and validation of a 3-D model to predict knee joint
loading during dynamic movement. Journal of biomechanical engineering. 2003; 125(6):864–74. PMID:
14986412
10. Reinbolt JA, Seth A, Delp SL. Simulation of human movement: applications using OpenSim. Iutam
Symposium on Human Body Dynamics. 2011; 2(0):186–98.
11. Rasmussen J, Torholm S, de Zee M. Computational analysis of the influence of seat pan inclination and
friction on muscle activity and spinal joint forces. Int J Ind Ergonom. 2009; 39(1):52–7.
12. Wu JZ, Chiou SS, Pan CS. Analysis of musculoskeletal loadings in lower limbs during stilts walking in
occupational activity. Ann Biomed Eng. 2009; 37(6):1177–89. doi: 10.1007/s10439-009-9674-5 PMID:
19296222
13. Panjabi MM, Pearson AM, Ito S, Ivancic PC, Wang JL. Cervical spine curvature during simulated whip-
lash. Clinical biomechanics (Bristol, Avon). 2004; 19(1):1–9.
14. Ivancic PC, Xiao M. Cervical spine curvature during simulated rear crashes with energy-absorbing seat.
The spine journal: official journal of the North American Spine Society. 2011; 11(3):224–33.
15. Quarrie KL, Cantu RC, Chalmers DJ. Rugby union injuries to the cervical spine and spinal cord. Sports
medicine (Auckland, NZ). 2002; 32(10):633–53.
16. Fuller CW, Brooks JH, Kemp SP. Spinal injuries in professional rugby union: a prospective cohort
study. Clin J Sport Med. 2007; 17(1):10–6. Epub 2007/02/17. doi: 10.1097/JSM.0b013e31802e9c28
PMID: 17304000
17. Preatoni E, Cazzola D, Stokes KA, England M, Trewartha G. Pre-binding prior to full engagement
improves loading conditions for front-row players in contested Rugby Union scrums. Scandinavian jour-
nal of medicine & science in sports. 2015.
18. Cazzola D, Preatoni E, Stokes KA, England ME, Trewartha G. A modified prebind engagement process
reduces biomechanical loading on front row players during scrummaging: a cross-sectional study of 11
elite teams. Br J Sports Med. 2015; 49(8):541–6. doi: 10.1136/bjsports-2013-092904 PMID: 24505041
19. Preatoni E, Stokes KA, England ME, Trewartha G. Engagement techniques and playing level impact
the biomechanical demands on rugby forwards during machine-based scrummaging. Br J Sports Med.
2015; 49(8):520–8. doi: 10.1136/bjsports-2013-092938 PMID: 24511085
20. Usman J, McIntosh AS, Frechede B. An investigation of shoulder forces in active shoulder tackles in
rugby union football. J Sci Med Sport. 2011; 14(6):547–52. doi: 10.1016/j.jsams.2011.05.006 PMID:
21689986
21. Kuster D, Gibson A, Abboud R, Drew T. Mechanisms of cervical spine injury in rugby union: a system-
atic review of the literature. Br J Sports Med. 2012; 46(8):550–4. doi: 10.1136/bjsports-2011-090360
PMID: 22378840
22. Dennison CR, Macri EM, Cripton PA. Mechanisms of cervical spine injury in rugby union: is it premature
to abandon hyperflexion as the main mechanism underpinning injury? Br J Sports Med. 2012; 46
(8):545–9. doi: 10.1136/bjsports-2012-091257 PMID: 22685172
Rugby Musculoskeletal Model
PLOS ONE | DOI:10.1371/journal.pone.0169329 January 4, 2017 21 / 24
23. Hamner SR, Seth A, Delp SL. Muscle contributions to propulsion and support during running. J Bio-
mech. 2010; 43(14):2709–16. Epub 2010/08/10. doi: 10.1016/j.jbiomech.2010.06.025 PMID: 20691972
24. Carbone V, Fluit R, Pellikaan P, van der Krogt MM, Janssen D, Damsgaard M, et al. TLEM 2.0—a com-
prehensive musculoskeletal geometry dataset for subject-specific modeling of lower extremity. J Bio-
mech. 2015; 48(5):734–41. doi: 10.1016/j.jbiomech.2014.12.034 PMID: 25627871
25. Nikooyan AA, Veeger HE, Chadwick EK, Praagman M, Helm FC. Development of a comprehensive
musculoskeletal model of the shoulder and elbow. Med Biol Eng Comput. 2011; 49(12):1425–35. Epub
2011/11/01. doi: 10.1007/s11517-011-0839-7 PMID: 22038240
26. Holzbaur KR, Murray WM, Delp SL. A model of the upper extremity for simulating musculoskeletal sur-
gery and analyzing neuromuscular control. Ann Biomed Eng. 2005; 33(6):829–40. Epub 2005/08/05.
PMID: 16078622
27. Dao TT, Pouletaut P, Charleux F, Lazary A, Eltes P, Varga PP, et al. Multimodal medical imaging (CT
and dynamic MRI) data and computer-graphics multi-physical model for the estimation of patient spe-
cific lumbar spine muscle forces. Data Knowl Eng. 2015; 96–97:3–18.
28. Senteler M, Weisse B, Rothenfluh DA, Snedeker JG. Intervertebral reaction force prediction using an
enhanced assembly of OpenSim models. Computer methods in biomechanics and biomedical engi-
neering. 2016; 19(5):538–48. doi: 10.1080/10255842.2015.1043906 PMID: 26031341
29. de Zee M, Hansen L, Wong C, Rasmussen J, Simonsen EB. A generic detailed rigid-body lumbar spine
model. J Biomech. 2007; 40(6):1219–27. doi: 10.1016/j.jbiomech.2006.05.030 PMID: 16901492
30. Vasavada AN, Li S, Delp SL. Influence of muscle morphometry and moment arms on the moment-gen-
erating capacity of human neck muscles. Spine (Phila Pa 1976). 1998; 23(4):412–22. Epub 1998/03/
28.
31. Christophy M, Faruk Senan NA, Lotz JC, O’Reilly OM. A musculoskeletal model for the lumbar spine.
Biomechanics and modeling in mechanobiology. 2012; 11(1–2):19–34. doi: 10.1007/s10237-011-0290-
6 PMID: 21318374
32. Delp SL, Loan JP, Hoy MG, Zajac FE, Topp EL, Rosen JM. An interactive graphics-based model of the
lower extremity to study orthopaedic surgical procedures. Biomedical Engineering, IEEE Transactions
on. 1990; 37(8):757–67.
33. Horsman MDK, Koopman HFJM, van der Helm FCT, Prose LP, Veeger HEJ. Morphological muscle
and joint parameters for musculoskeletal modelling of the lower extremity. Clin Biomech. 2007; 22
(2):239–47.
34. Scheys L, Spaepen A, Suetens P, Jonkers I. Calculated moment-arm and muscle-tendon lengths dur-
ing gait differ substantially using MR based versus rescaled generic lower-limb musculoskeletal models.
Gait Posture. 2008; 28(4):640–8. doi: 10.1016/j.gaitpost.2008.04.010 PMID: 18534855
35. Fuller CW, Taylor AE, Brooks JH, Kemp SP. Changes in the stature, body mass and age of English pro-
fessional rugby players: a 10-year review. J Sports Sci. 2013; 31(7):795–802. Epub 2012/12/19. doi:
10.1080/02640414.2012.753156 PMID: 23244349
36. Anthropometric Reference Data for Children and Adults: United States, 2007–2010 [Internet]. 2012.
37. Delp SL, Loan JP, Hoy MG, Zajac FE, Topp EL, Rosen JM. An interactive graphics-based model of the
lower extremity to study orthopaedic surgical procedures. IEEE transactions on bio-medical engineer-
ing. 1990; 37(8):757–67. doi: 10.1109/10.102791 PMID: 2210784
38. Hicks JL, Uchida TK, Seth A, Rajagopal A, Delp SL. Is my model good enough? Best practices for verifi-
cation and validation of musculoskeletal models and simulations of movement. Journal of biomechani-
cal engineering. 2015; 137(2):020905. doi: 10.1115/1.4029304 PMID: 25474098
39. Fung M, Kato S, Barrance PJ, Elias JJ, McFarland EG, Nobuhara K, et al. Scapular and clavicular kine-
matics during humeral elevation: a study with cadavers. Journal of shoulder and elbow surgery / Ameri-
can Shoulder and Elbow Surgeons [et al]. 2001; 10(3):278–85. Epub 2001/06/16.
40. Ludewig PM, Phadke V, Braman JP, Hassett DR, Cieminski CJ, LaPrade RF. Motion of the shoulder
complex during multiplanar humeral elevation. J Bone Joint Surg Am. 2009; 91(2):378–89. doi: 10.
2106/JBJS.G.01483 PMID: 19181982
41. Wu G, van der Helm FCT, Veeger HEJ, Makhsous M, Van Roy P, Anglin C, et al. ISB recommendation
on definitions of joint coordinate systems of various joints for the reporting of human joint motion—Part
II: shoulder, elbow, wrist and hand. Journal of Biomechanics. 2005; 38(5):981–92. PMID: 15844264
42. Ludewig PM, Hassett DR, Laprade RF, Camargo PR, Braman JP. Comparison of scapular local coordi-
nate systems. Clinical biomechanics (Bristol, Avon). 2010; 25(5):415–21.
43. Ivancic PC, Ito S, Panjabi MM. Dynamic sagittal flexibility coefficients of the human cervical spine. Acci-
dent; analysis and prevention. 2007; 39(4):688–95. doi: 10.1016/j.aap.2006.10.015 PMID: 17140545
44. de Leva P. Adjustments to Zatsiorsky-Seluyanov’s segment inertia parameters. J Biomech. 1996; 29
(9):1223–30. PMID: 8872282
Rugby Musculoskeletal Model
PLOS ONE | DOI:10.1371/journal.pone.0169329 January 4, 2017 22 / 24
45. Ivancic PC, Panjabi MM, Ito S. Cervical spine loads and intervertebral motions during whiplash. Traffic
Inj Prev. 2006; 7(4):389–99. doi: 10.1080/15389580600789127 PMID: 17114097
46. Nightingale RW, Camacho DL, Armstrong AJ, Robinette JJ, Myers BS. Inertial properties and loading
rates affect buckling modes and injury mechanisms in the cervical spine. J Biomech. 2000; 33(2):191–
7. Epub 2000/02/01. PMID: 10653032
47. Yoganandan N, Pintar FA, Zhang J, Baisden JL. Physical properties of the human head: mass, center
of gravity and moment of inertia. J Biomech. 2009; 42(9):1177–92. doi: 10.1016/j.jbiomech.2009.03.
029 PMID: 19428013
48. Yeadon MR. The simulation of aerial movement—II. A mathematical inertia model of the human body. J
Biomech. 1990; 23(1):67–74. Epub 1990/01/01. PMID: 2307693
49. Fice JB, Siegmund GP, Blouin JS. Prediction of three dimensional maximum isometric neck strength.
Ann Biomed Eng. 2014; 42(9):1846–52. doi: 10.1007/s10439-014-1046-0 PMID: 24893597
50. Vasavada AN, Li S, Delp SL. Three-dimensional isometric strength of neck muscles in humans. Spine
(Phila Pa 1976). 2001; 26(17):1904–9.
51. Kamibayashi LK, Richmond FJ. Morphometry of human neck muscles. Spine (Phila Pa 1976). 1998; 23
(12):1314–23.
52. Du Toit DE, Venter DJL, Buys FJ, Olivier PE. Kinetics of rugby union scrumming in under-19 schoolboy
rugby forwards. South African Journal for Research in Sport, Physical Education and Recreation. 2004;
26(2):33–50.
53. Sherman MA, Seth A, Delp SL. What Is a Moment Arm? Calculating Muscle Effectiveness in Bio-
mechanical Models Using Generalized Coordinates. Proceedings of the ASME Design Engineering
Technical Conferences ASME Design Engineering Technical Conferences. 2013;2013.
54. Thelen DG. Adjustment of muscle mechanics model parameters to simulate dynamic contractions in
older adults. Journal of biomechanical engineering. 2003; 125(1):70–7. Epub 2003/03/29. PMID:
12661198
55. Miura T, Panjabi MM, Cripton PA. A method to simulate in vivo cervical spine kinematics using in vitro
compressive preload. Spine (Phila Pa 1976). 2002; 27(1):43–8. Epub 2002/01/24.
56. Olivier PE, Du Toit DE. Isokinetic neck strength profile of senior elite rugby union players. J Sci Med
Sport. 2008; 11(2):96–105. doi: 10.1016/j.jsams.2007.01.009 PMID: 17560830
57. Cazzola D, Stone B, Holsgrove TP, Trewartha G, Preatoni E. Spinal muscle activity in simulated rugby
union scrummaging is affected by different engagement conditions. Scandinavian journal of medicine &
science in sports. 2015.
58. Thelen DG, Anderson FC, Delp SL. Generating dynamic simulations of movement using computed
muscle control. Journal of Biomechanics. 2003; 36(3):321–8. PMID: 12594980
59. Thelen DG, Anderson FC. Using computed muscle control to generate forward dynamic simulations of
human walking from experimental data. J Biomech. 2006; 39(6):1107–15. doi: 10.1016/j.jbiomech.
2005.02.010 PMID: 16023125
60. Zajac FE. Muscle and tendon: properties, models, scaling, and application to biomechanics and motor
control. Crit Rev Biomed Eng. 1989; 17(4):359–411. PMID: 2676342
61. Delp SL, Anderson FC, Arnold AS, Loan P, Habib A, John CT, et al. OpenSim: open-source software to
create and analyze dynamic simulations of movement. IEEE transactions on bio-medical engineering.
2007; 54(11):1940–50. Epub 2007/11/21. doi: 10.1109/TBME.2007.901024 PMID: 18018689
62. Penning L. Normal movements of the cervical spine. AJR American journal of roentgenology. 1978; 130
(2):317–26. Epub 1978/02/01. doi: 10.2214/ajr.130.2.317 PMID: 414586
63. Penning L, Wilmink JT. Rotation of the cervical spine. A CT study in normal subjects. Spine (Phila Pa
1976). 1987; 12(8):732–8. Epub 1987/10/01.
64. Dvorak J, Froehlich D, Penning L, Baumgartner H, Panjabi MM. Functional radiographic diagnosis of
the cervical spine: flexion/extension. Spine (Phila Pa 1976). 1988; 13(7):748–55. Epub 1988/07/01.
65. Dvorak J, Antinnes JA, Panjabi M, Loustalot D, Bonomo M. Age and gender related normal motion of
the cervical spine. Spine (Phila Pa 1976). 1992; 17(10 Suppl):S393–8. Epub 1992/10/01.
66. Dvorak J, Hayek J, Zehnder R. CT-functional diagnostics of the rotatory instability of the upper cervical
spine. Part 2. An evaluation on healthy adults and patients with suspected instability. Spine (Phila Pa
1976). 1987; 12(8):726–31. Epub 1987/10/01.
67. Dvorak J, Panjabi MM, Grob D, Novotny JE, Antinnes JA. Clinical validation of functional flexion/exten-
sion radiographs of the cervical spine. Spine (Phila Pa 1976). 1993; 18(1):120–7. Epub 1993/01/01.
68. Lark SD, McCarthy P. The effects of a rugby playing season on cervical range of motion. J Sports Sci.
2010; 28(6):649–55. doi: 10.1080/02640411003631968 PMID: 20397098
Rugby Musculoskeletal Model
PLOS ONE | DOI:10.1371/journal.pone.0169329 January 4, 2017 23 / 24
69. Ackland DC, Merritt JS, Pandy MG. Moment arms of the human neck muscles in flexion, bending and
rotation. J Biomech. 2011; 44(3):475–86. doi: 10.1016/j.jbiomech.2010.09.036 PMID: 21074159
70. Seminati E, Cazzola D, Preatoni E, Trewartha G. Specific tackling situations affect the biomechanical
demands experienced by rugby union players. Sports biomechanics / International Society of Biome-
chanics in Sports. 2016:1–18.
71. Scheys L, Desloovere K, Suetens P, Jonkers I. Level of subject-specific detail in musculoskeletal mod-
els affects hip moment arm length calculation during gait in pediatric subjects with increased femoral
anteversion. J Biomech. 2011; 44(7):1346–53. doi: 10.1016/j.jbiomech.2011.01.001 PMID: 21295307
72. Blana D, Hincapie JG, Chadwick EK, Kirsch RF. A musculoskeletal model of the upper extremity for use
in the development of neuroprosthetic systems. J Biomech. 2008; 41(8):1714–21. doi: 10.1016/j.
jbiomech.2008.03.001 PMID: 18420213
73. Garner BA, Pandy MG. Musculoskeletal model of the upper limb based on the visible human male data-
set. Computer methods in biomechanics and biomedical engineering. 2001; 4(2):93–126. doi: 10.1080/
10255840008908000 PMID: 11264863
74. Seth A, Matias R, Veloso AP, Delp SL. A Biomechanical Model of the Scapulothoracic Joint to Accu-
rately Capture Scapular Kinematics during Shoulder Movements. PloS one. 2016; 11(1):e0141028. doi:
10.1371/journal.pone.0141028 PMID: 26734761
75. Bolsterlee B, Veeger HE, van der Helm FC. Modelling clavicular and scapular kinematics: from mea-
surement to simulation. Med Biol Eng Comput. 2014; 52(3):283–91. doi: 10.1007/s11517-013-1065-2
PMID: 23543279
76. de Groot JH, Brand R. A three-dimensional regression model of the shoulder rhythm. Clinical biome-
chanics (Bristol, Avon). 2001; 16(9):735–43.
77. Martelli S, Valente G, Viceconti M, Taddei F. Sensitivity of a subject-specific musculoskeletal model to
the uncertainties on the joint axes location. Computer methods in biomechanics and biomedical engi-
neering. 2015; 18(14):1555–63. doi: 10.1080/10255842.2014.930134 PMID: 24963785
78. Reinbolt JA, Haftka RT, Chmielewski TL, Fregly BJ. Are patient-specific joint and inertial parameters
necessary for accurate inverse dynamics analyses of gait? IEEE transactions on bio-medical engineer-
ing. 2007; 54(5):782–93. doi: 10.1109/TBME.2006.889187 PMID: 17518274
79. Sheets AL, Corazza S, Andriacchi TP. An automated image-based method of 3D subject-specific body
segment parameter estimation for kinetic analyses of rapid movements. Journal of biomechanical engi-
neering. 2010; 132(1):011004. doi: 10.1115/1.4000155 PMID: 20524742
80. Davidson PL, Wilson SJ, Wilson BD, Chalmers DJ. Estimating subject-specific body segment parame-
ters using a 3-dimensional modeller program. J Biomech. 2008; 41(16):3506–10. doi: 10.1016/j.
jbiomech.2008.09.021 PMID: 18995857
81. Matias R, Andrade C, Veloso AP. A transformation method to estimate muscle attachments based on
three bony landmarks. J Biomech. 2009; 42(3):331–5. doi: 10.1016/j.jbiomech.2008.11.027 PMID:
19136112
82. Bolsterlee B, Vardy AN, van der Helm FC, DirkJan Veeger HE. The effect of scaling physiological
cross-sectional area on musculoskeletal model predictions. J Biomech. 2015; 48(10):1760–8. doi: 10.
1016/j.jbiomech.2015.05.005 PMID: 26050956
83. de Bruijn E, van der Helm FCT, Happee R. Analysis of isometric cervical strength with a nonlinear mus-
culoskeletal model with 48 degrees of freedom. Multibody System Dynamics. 2015; 36(4):339–62.
84. Pizzolato C, Lloyd DG, Sartori M, Ceseracciu E, Besier TF, Fregly BJ, et al. CEINMS: A toolbox to
investigate the influence of different neural control solutions on the prediction of muscle excitation and
joint moments during dynamic motor tasks. J Biomech. 2015; 48(14):3929–36. doi: 10.1016/j.jbiomech.
2015.09.021 PMID: 26522621
Rugby Musculoskeletal Model
PLOS ONE | DOI:10.1371/journal.pone.0169329 January 4, 2017 24 / 24