Background. Spastic muscles of patients with cerebral palsy (CP) are considered structurally as shortened muscles, that produce high force in short muscle lengths. Yet, previous intraoperative studies in which muscles’ forces are measured directly as a function of joint angle showed consistently that spastic knee flexor muscles produce a low percentage of their maximum force in flexed knee positions. They also showed effects of epimuscular myofascial force transmission (EMFT): simultaneous activation of different muscles elevated target muscle's force. However, quantification of spastic muscle's force - muscle-tendon unit length (l. MTU. ) data during gait is lacking. Aim. Combining intraoperative experiments with participants’ musculoskeletal models developed based on their gait analyses, we aimed to test the following hypotheses: activated spastic semitendinosus (ST) muscle (1) operates at short l. MTU. 's during gait, forces are (2) low at short l. MTU. 's and (3) increase by co-activating other muscles. Methods. Ten limbs of seven children with CP (GMFCS-II) were tested. Pre-surgery, gait analyses were conducted. Intraoperatively, isometric spastic ST distal forces were measured in ten hip-knee joint angle combinations, in two conditions: (i) activation of the ST individually and (ii) simultaneously with the gracilis, biceps femoris, and rectus femoris muscles endorsing EMFT. In OpenSim, gait_2392 model was used for each limb to (a) calculate l. MTU. per each hip and knee angle combination and the gait relevant l. MTU. range, and (b) analyze gait relevant spastic muscle force - l. MTU. data. Two-way ANOVA was used to compare the patients’ l. MTU. to those of the seven age-matched typically developing (TD) children. l. MTU. values were normalized for the participants’ thigh length. (a) was used to test hypothesis (1) and (b) to test hypotheses (2) and (3): in condition (i), the percent of peak force exerted at the shortest l. MTU. calculated per limb was used as a metric for (2). In condition (ii), mean percent change in muscle force calculated within gait-relevant l. MTU. range was used as a metric for (3). Results. Modeling showed that l. MTU. of spastic ST during gait is shorter on average by 14.1% compared to TD. The ST active force at the shortest gait-relevant l. MTU. was 68.6 (20.6)% (39.9–99.2%) of the peak force. Simultaneous activation of other muscles caused substantial increases in force (minimally by 11.1%, up to several folds, with an exception for one limb). Therefore, only the first and third hypotheses were confirmed. Conclusion. The modeling showed in concert with the clinical considerations that spastic ST may be a shortened muscle that produces high force in short muscle lengths. However, this contrasts intraoperative data, which shows only low forces in flexed knee positions. Note that, the model does not distinguish the muscle-belly and tendon lengths. Therefore, it cannot isolate shorter muscle length and how this compares to the data of TD children remains unknown. Yet, the effects of co-activation of other muscles shown intraoperatively to cause an increase of the spastic ST's force are observed also in muscle force - l. MTU. data characterizing gait. Therefore, if indeed spastic ST produces high forces in short muscle-belly lengths alone, elevated forces due to co-activation of other muscles may be considered as a contributor to the patients’ pathological gait. Otherwise, such EMFT effect may be the main determinant of the pathological condition

Altered mechanical loading is a widely suggested, but poorly understood potential cause of cartilage degeneration in osteoarthritis. In rodents, osteoarthritis is induced following destabilization of the medial meniscus (DMM). This study estimates knee kinematics and contact forces in rats with DMM to gain better insight into the specific mechanisms underlying disease development in this widely-used model. Unilateral knee surgery was performed in adult male Sprague-Dawley rats (n=5 with DMM, n=5 with sham surgery). Radio-opaque beads were implanted on their femur and tibia. 8 weeks following knee surgery, rat gait was recorded using the 3D²YMOX setup (Sanctorum et al. 2019, simultaneous acquisition of biplanar XRay videos and ground reaction forces). 10 trials (1 per rat) were calibrated and processed in XMALab (Knörlein et al. 2016). Hindlimb bony landmarks were labeled on the XRay videos using transfer learning (Deeplabcut, Mathis et al. 2019; Laurence-Chasen et al. 2020). A generic OpenSim musculoskeletal model of the rat hindlimb (Johnson et al. 2008) was adapted to include a 3-degree-of-freedom knee. Inverse kinematics, inverse dynamics, static optimization of muscle forces, and joint reaction analysis were performed. In rats with DMM, knee adduction was lower compared to sham surgery. Ground reaction forces were less variable with DMM, resulting in less variability in joint external moments. The mediolateral ground reaction force was lower, resulting in lower hip adduction moment, thus less force was produced by the rectus femoris. Rats with DMM tended to break rather than propel, resulting in lower hip flexion moment, thus less force was produced by the semimembranosus. These results are consistent with lower knee contact forces in the anteroposterior and axial directions. These preliminary data indicate no overloading of the knee joint in rats with DMM, compared with sham surgery. We are currently expanding our workflow to finite element analysis, to examine mechanical cues in the cartilage of these rats (Fig1G)

Anterior cruciate ligament reconstruction (ACLR) using a semitendinosus (ST) autograft, with or without gracilis (GR), results in donor muscle atrophy and varied tendon regeneration. The effects of harvesting these muscles on muscle moment arm and torque generating capacity have not been well described. This study aimed to determine between-limb differences (ACLR vs uninjured contralateral) in muscle moment arm and torque generating capacity across a full range of hip and knee motions. A secondary analysis of magnetic resonance imaging was undertaken from 8 individuals with unilateral history of ST-GR ACLR with complete ST tendon regeneration. All hamstring muscles and ST tendons were manually segmented. Muscle length (cm), peak cross-sectional area (CSA) (cm. 2. ), and volume (cm. 3. ) were measured in ACLR and uninjured contralateral limbs. OpenSim was used to simulate and evaluate the mechanical consequences of changes in normalised moment arm (m) and torque generating capacity (N.m) between ACLR and uninjured contralateral limbs. Compared to uninjured contralateral limbs, regenerated ST tendon re-insertion varied proximal (+) (mean = 0.66cm, maximum = 3.44cm, minimum = −2.17cm, range = 5.61cm) and posterior (+) (mean = 0.38cm maximum = 0.71cm, minimum = 0.02cm, range = 0.69cm) locations relative to native anatomical positions. Compared to uninjured contralateral limbs, change in ST tendon insertion point in ACLR limbs resulted in 2.5% loss in peak moment arm and a 3.4% loss in peak torque generating capacity. Accounting for changes to both max isometric force and ST moment arm, the ST had a 14.8% loss in peak torque generating capacity. There are significant deficits in ST muscle morphology and insertion points following ST-GR ACLR. The ST atrophy and insertion point migration following ACLR may affect force transmission and distribution within the hamstrings and contribute to persistent deficits in knee flexor and internal rotator strength

This study aimed to quantify self-reported outcomes and walking gait biomechanics in patients following primary and revision THA. The specific goals of this study were to investigate: (i) if primary and revision THA patients have comparable preoperative outcomes; and (2) if revision THA patients have worse postoperative outcomes than primary THA patients. Forty-three patients undergoing primary THA for osteoarthritis and 23 patients undergoing revision THA were recruited and followed longitudinally for their first 12 postoperative months. Reasons for revision were loosening (73%), dislocation (9%), and infection (18%). Patients completed the Hip dysfunction and Osteoarthritis Outcome Score (HOOS), and underwent gait analysis preoperatively, and at 3 and 12 months postoperatively. A 10 camera motion analysis system (V5 Vantage, Vicon, UK) recorded marker trajectories (100 Hz) during walking at self- selected speeds. A generic lower-body musculoskeletal model (Gait2392) was scaled using principal component analysis [1] and the inverse kinematics tool in Opensim 3.3 was used to compute joint angles for the lower limbs in the sagittal plane. Independent samples t-test were used to compare patient reported outcomes between the primary and revision groups at each timepoint. Statistical parametric mapping was used to compare gait patterns between the two groups at each timepoint. Preoperatively, patients undergoing primary THA reported significantly worse pain (p<0.001), symptoms (p<0.001), function (p<0.001), and quality of life (p=0.004). No differences were observed at 3 and 12 months postoperatively between patients who had received a primary or revision THA. The only observed difference in gait pattern was that patients with a revision THA had reduced hip extension at 3 months, but no differences were observed preoperatively and 12 months. Despite the suggestions in the literature that revision THA is bound to have worse outcomes compared to primary THA, we found no differences in in patient-reported outcomes and gait patterns at 12 months postoperatively. This suggests that it may be possible, in some circumstances, for patients following revision THA to achieve similar outcomes to their peers undergoing primary THA

Introduction and Objective. Clinically, it is considered that spastic muscles of patients with cerebral palsy (CP) are shortened, and produce higher force in shorter muscle lengths. Yet, direct quantification of spastic muscles’ forces is rare. Remarkably, previous intraoperative tests in which muscle forces are measured directly as a function of joint angle showed for spastic gracilis (GRA) that its passive forces are low, and only a small percentage of its maximum active force is measured in flexed knee positions. However, the relationship of force characteristics of spastic GRA with its muscle-tendon unit length (l. MTU. ) is unknown. Combining intraoperative experiments with participants’ musculoskeletal models developed based on their gait analyses, we aimed to test if spastic GRA muscle (1) operates at short l. MTU. compared to that of typically developing (TD) children, and exerts higher (2) passive and (3) active forces at shorter lengths, within gait-relevant l. MTU. range. Materials and Methods. Ten limbs of seven children with CP (GMFCS-II) were tested. Pre-surgery, gait analyses were conducted. Intraoperatively, isometric spastic GRA distal forces were measured in ten hip-knee joint angle combinations, in two conditions: (i) passive state and (ii) maximal activation of the GRA exclusively. In OpenSim, gait_2392 model was used for each limb to calculate l. MTU. 's per each hip and knee angle combination and the gait-relevant l. MTU. range, and to analyze gait relevant spastic muscle force - l. MTU. data. l. MTU. values were normalized for the participants’ thigh lengths. Two-way ANOVA was used to compare the patients’ l. MTU. to those of the seven age-matched TD children to test the first hypothesis. In order to test the second and the third hypotheses, Spearman's rank correlation coefficient (ρ) was calculated to seek a correlation between the muscle's operational length (represented by mean l. MTU. within gait cycle) and muscular force characteristics (the percent force at shortest l. MTU. of peak force, either in passive or in active conditions) within gait-relevant l. MTU. range. Results. ANOVA showed that l. MTU. 's of spastic GRA are shorter (on average by 15.4%) compared to those of TD. At the shortest gait-relevant l. MTU. , the GRA passive force was 84.6 (13.7)% of the peak passive force; and the active force was 55.8 (33.9)% of the peak active force. Passive state forces show an increase at longer lengths, whereas active state force characteristics vary in a patient-specific way. Spearman's rank correlation indicated weak correlations between muscle's operational length and muscular force characteristics (ρ= −0.30 P= 0.40, and ρ= −0.27 P= 0.45, for passive and active states, respectively). Therefore, only the first hypothesis was confirmed. Conclusions. Novel muscle force - l. MTU. data for spastic GRA were obtained using intraoperative data and modelling combined. The modelling showed in concert with the clinical considerations that spastic GRA may be a shortened muscle. However, because the model does not distinguish the muscle-belly and tendon lengths, it cannot isolate shorter muscle belly length and how this compares to the data of TD children remains unknown. Moreover, the absence of a strong correlation between shorter operational muscle length and higher force production either in passive or in active conditions highlights the influence of other factors (e.g., muscle structural proteins, and muscle mechanical characteristics including intermuscular interactions etc.) on the pathology rather than ascribing it solely to the length of a spastic muscle itself

Cam-type femoroacetabular impingement (FAI) is a common cause for athletic hip injury and early hip osteoarthritis. Although corrective cam FAI surgery can improve patient reported outcome measures (PROMs), it is not clear how surgery affects muscle forces and hip joint loading. Surgery for FAI may redistribute muscle forces and contact forces at the hip joint during routine activities. The purpose of this study was to examine the muscle contributions and hip contact forces during gait in patients prior and after two years of undergoing surgery for cam FAI. Kinematics and kinetics were recorded in 11 patients with symptomatic cam FAI as they completed a gait task. Muscle and hip contact forces during the stance phase were estimated using musculoskeletal modelling and static optimization in OpenSim. All patients reported improvements in PROMs. Post-operatively, patients showed reduced forces in the long head of the biceps femoris at ipsilateral foot-strike and in the rectus femoris at the contralateral foot-strike. The reduced muscle forces decreased sagittal hip moment but did not change hip contact forces. This was the first study to evaluate hip muscle and contact forces in FAI patients post-operatively. Although hip contact forces are not altered following surgery, muscle forces are decreased even after two years. These findings can provide guidance in optimizing recovery protocols after FAI surgery to improve hip flexor and extensor muscle forces

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

Introduction. Modification in joint loading, and specifically shear stress, is found to be an important mechanical factor in the development of osteoarthritis (OA). Cartilage shear stresses can be investigated using finite element (FE) modelling, where typically in vivo joint loading as measured by an instrumented hip prosthesis is used as boundary condition. However, subject-specific gait characteristics substantially affect joint loading. The goal of this study is to investigate the effect of subject-specific joint loading as calculated using a subject-specific musculoskeletal model and integrated motion capture data on acetabular shear stress. Methods. Three healthy control subjects walked at self-selected speed while measuring marker trajectories (Vicon, Oxford Metrics, UK) and force data (two AMTI force platforms; Watertown, MA). A subject-specific MRI-based musculoskeletal model consisting of 14 segments, 19 degrees of freedom and 88 musculotendon actuators, and including wrapping surfaces around the hip joint, was used. All analyses were performed in OpenSim 3.1. The model was scaled to the dimensions of each subject using the marker positions of a static pose. A kalman smoother procedure was used to calculate joint angles. Muscle forces were calculated using static optimization, minimizing the sum of squared muscle activations, and hip contact forces (HCF) were calculated and normalized to body weight (BW). To calculate shear stress, HCFs and joint angles calculated during the stance phase of gait were imposed to a hip finite element model (hip_n10rb) using FFEbio 2.5. In the model, femoral and acetabular cartilage were represented using the Mooney-Rivlin formulation (c1=6.817, bulk modulus=1358.86) and the pelvis and femur bones as rigid bodies. Peak HCF as well as maximal acetabular shear stress, magnitude and location, and the HCF at the time of maximal shear stress were compared between subjects. Results. Maximal shear stress was lower for S3 compared to S1 and S2 (9.14, 9.48 and 7.14 MPa for S1, S2 and S3 respectively). Nevertheless, HCF at the time instance of peak stress as well as peak HCF were highest for S3 (S1: 2.40/4.54 BW, S2: 2.97/4.78 BW and S3: 3.13/6.46 BW respectively). Maximal shear stress also occurred earlier in the stance phase for S3 compared to S1 and S2 (31, 26 and 11% of the stance phase for S1, S2 and S3 respectively). In addition, the location of the peak maximal shear stress was found to be more superior for S3. Discussion. Subject-specific loading patterns clearly influence the calculated maximal shear stress in the acetabular cartilage, affecting both the magnitude and the location of the stress. In addition, higher shear stresses are not coinciding with higher HCFs. This finding highlights the need of subject-specific rather than generic loading patterns when assessing cartilage shear stresses and associated risk in OA development in individual patients

Several studies have shown that gait kinematics[1–3] and hip contact forces (HCFs)[4, 5] of patients following total hip arthroplasty (THA) do not return to normal, although improvements in kinematics are found compared to the pre-surgery. However, the evolution of HCFs after surgery has not been investigated. The goal of this study is to evaluate HCFs during gait in OA patients before and at 2 evaluation moments post-THA. Fourteen unilateral hip OA patients before and 3- and 12-months post-THA surgery walked at self-selected speed, as well as 18 healthy control subjects. 3D marker trajectories were captured using Vicon (Oxford Metrics, UK) and force data was measured using two AMTI force platforms (Watertown, MA). A musculoskeletal model consisting of 14 segments, 19 degrees of freedom and 88 musculotendon actuators and including wrapping surfaces around the hip joint was used[6]. All analyses were performed in OpenSim 3.1[7]. The model was scaled to the dimensions of each subject using the marker positions of a static pose. A kalman smoother procedure was used to calculate joint angles[8]. Muscle forces were calculated using static optimization, minimizing the sum of squared muscle activations. HCFs were calculated and normalized to body weight (BW). First and second peak HCFs were determined and used for statistical analysis. To determine differences between HCFs of OA patients at the different evaluation moments, a Friedman test was used. In case of a significant difference, post-hoc rank-based multiple comparison tests with a Bonferonni adjustment was used. To compare controls and patients at each evaluation moment separate Man-Whitney U tests were used. Differences in HCFs between the affected and non-affected legs were expressed by a symmetry index (SI), i.e. the ratio between the HCFs of the affected leg over the non-affected leg, averaged over the stance phase of the gait cycle. At the first and second HCF peaks, no significant differences were found between pre-, 3- and 12-months post-surgery (first peak average HCF: 2.68, 2.72 and 2.78BW respectively; second peak average HCF: 3.21, 3.83 and 3.77BW respectively). Compared to controls, significant differences are found for all evaluation moments at the first and second HCF peaks (average HCF controls: 3.43 and 5.15BW respectively). The SI was below 1 pre- and 3-months post-surgery (0.88 and 0.85 respectively), indicating decreased loading of the affected compared to the non-affected leg. At 12-months post-surgery SI was close to 1 (0.98). As reported before[4, 5], first or second peak HCFs do not return to normal after THA. Although HCFs increase after THA compared to pre-surgery, significant differences with controls remain. Surprisingly, no significant differences are found between the different evaluation moments of the patients, indicating no clear improvements are found after THA. Further, average HCF peaks at 3- and 12-months post-surgery are similar, indicating no further improvements are found 3-months post-surgery. However, the SI was above or close to 1 at 12-months post-surgery, indicating hip loading evolved to a more symmetrical loading 12-months post-surgery

Children with cerebral palsy (CP) often present femoral bone deformities not accounted for in generic musculoskeletal models [1,2]. MRI-based models can be used to include subject-specific muscle paths [3,4], although this is a time-demanding process. Recently, non-rigid deformation techniques have been used to transform generic bone geometry, including muscle points, onto personalized bones [5]. However, it is still unknown to what extent such an approximation of subject-specific detail affects calculated hip contact forces (HCFs) during gait in CP children. Seven children diagnosed with diplegic CP walked independently at self-selected speed. 3D marker trajectories were captured using Vicon (Oxford Metrics, UK) and force data was measured using two AMTI force platforms (Watertown, MA). MR-images were acquired (Philips Ingenia 1.5T) of all subjects lying supine. Firstly, a generic model [6] was scaled using the marker positions of a static pose. Secondly, a MRI-model containing the subject-specific bone structures and muscle paths of all hip and upper leg muscles was created [3]. Thirdly, the generic femur and pelvis geometries and muscle points were transformed onto the image-based femur and pelvis using an advanced non-rigid deformation procedure (Materialise N.V.). For all models, further analyses were performed in OpenSim 3.1 [7]. A kalman smoother procedure was used to calculate joint angles [8]. Muscle forces were calculated using a static optimization minimizing the sum of squared muscle activities. Next, HCFs were calculated and normalized to body weight (BW). First and second peak HCFs were determined and used for a Kruskal-Wallis test to determine differences between models. In case of a significant difference, a post-hoc rank-based multiple comparison test with Bonferonni adjustment was used. Further, average absolute differences in muscle points between the models was calculated, as well as average differences in moment arm lengths (MALs), reflecting muscle function. Where the scaled generic muscle points differed on average 2.49cm from the MRI points, the non-rigidly deformed points differed 1.54cm from the MRI muscle points. Specifically, the tensor fascia latae differed most between the deformed and MRI models (11.7cm). When considering MALs, the gluteii muscles present an altered function for the generic and deformed models compared to the MRI model for all degrees of freedom of the hip at the time of both HCF peaks. The differences between models resulted in a significantly increased second peak HCF for the MRI models compared to the generic models (first peak average HCF: 3.88BW, 3.95BW and 4.90BW; second peak average HCF: 3.03BW, 4.89BW and 5.32BW for the generic, MRI and non-rigidly deformed models respectively). Although not significantly different, the deformed models calculated slightly increased HCFs compare to the MRI models. The generic models underestimated HCFs compared to the MRI models, while the non-rigidly deformed models slightly overestimated HCFs. However, differences between the deformed and MRI models in terms of muscle points and MALs remain, specifically for the gluteii muscles. Therefore, further user-guided modification of the model based on MR-images will be necessary