Soft tissue structure modelling for use in orthopaedic applications and musculoskeletal biomechanics

We present our methodology for the three-dimensional anatomical and geometrical description of soft tissues, relevant for orthopaedic surgical applications and musculoskeletal biomechanics. The technique involves the segmentation and geometrical description of muscles and neurovascular structures from high-resolution computer tomography scanning for the reconstruction of generic anatomical models. These models can be used for quantitative interpretation of anatomical and biomechanical aspects of different soft tissue structures. This approach should allow the use of these data in other application fields, such as musculoskeletal modelling, simulations for radiation therapy, and databases for use in minimally invasive, navigated and robotic surgery

Effect of a rigid ankle-foot orthosis on hamstring length in children with hemiplagia

Eighteen children with hemiplegia, mean age 8 years 5 months, underwent gait analysis and musculoskeletal modelling using specially designed software. The maximum lengths of the hamstrings were determined for each child walking in and out of an ankle–foot orthosis (AFO). The muscles were deemed to be short if shorter than the normal average – 1SD. In bare feet 8 participants had short medial hamstrings with a higher proportion of these in the less involved individuals. All participants showed an increase in maximum hamstring length when wearing an AFO. In all but one child this was sufficient to restore hamstring length to within normal limits. These finding suggest that hamstring pathology in hemiplegic gait is usually secondary to more distal lower limb pathology

The influence of barefoot and shod running on Triceps-surae muscle strain characteristics

The aim of the current investigation was to determine the effects of barefoot and shod running on the kinematics of the Triceps-Surae muscle group. Twelve male participants ran at 4.0 m.s-1 (± 5%) in both barefoot and shod conditions. Kinematics were measured using an eight-camera motion analysis system. Muscle kinematics from the lateral Gastrocnemius, medial Gastrocnemius and Soleus were obtained using musculoskeletal modelling software (Opensim v3.2). The results showed that muscle strain for the lateral Gastrocnemius (barefoot = 1.10 & shod = 0.33 %), medial Gastrocnemius (barefoot = 1.07 & shod = 0.32 %) and Soleus (barefoot = 3.43 & shod = 2.18 %) were significantly larger for the barefoot condition. Given the proposed association between the extent of muscle strain and the etiology of chronic muscle strain pathologies, the current investigation shows that running barefoot may place runners at greater risk from Triceps-Surae strain injuries

Comments and corrections on 3D modeling studies of locomotor muscle moment arms in archosaurs

© 2015 Bates et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, reproduction and adaptation in any medium and for any purpose provided that it is properly attributed. For attribution, the original author(s), title, publication source (PeerJ) and either DOI or URL of the article must be cited. The attached file is the published version of the article

Musculoskeletal modelling of an ostrich (Struthio camelus) pelvic limb: influence of limb orientation on muscular capacity during locomotion

We developed a three-dimensional, biomechanical computer model of the 36 major pelvic limb muscle groups in an ostrich (Struthio camelus) to investigate muscle function in this, the largest of extant birds and model organism for many studies of locomotor mechanics, body size, anatomy and evolution. Combined with experimental data, we use this model to test two main hypotheses. We first query whether ostriches use limb orientations (joint angles) that optimize the moment-generating capacities of their muscles during walking or running. Next, we test whether ostriches use limb orientations at mid-stance that keep their extensor muscles near maximal, and flexor muscles near minimal, moment arms. Our two hypotheses relate to the control priorities that a large bipedal animal might evolve under biomechanical constraints to achieve more effective static weight support. We find that ostriches do not use limb orientations to optimize the moment-generating capacities or moment arms of their muscles. We infer that dynamic properties of muscles or tendons might be better candidates for locomotor optimization. Regardless, general principles explaining why species choose particular joint orientations during locomotion are lacking, raising the question of whether such general principles exist or if clades evolve different patterns (e.g., weighting of muscle force–length or force–velocity properties in selecting postures). This leaves theoretical studies of muscle moment arms estimated for extinct animals at an impasse until studies of extant taxa answer these questions. Finally, we compare our model’s results against those of two prior studies of ostrich limb muscle moment arms, finding general agreement for many muscles. Some flexor and extensor muscles exhibit self-stabilization patterns (posture-dependent switches between flexor/extensor action) that ostriches may use to coordinate their locomotion. However, some conspicuous areas of disagreement in our results illustrate some cautionary principles. Importantly, tendon-travel empirical measurements of muscle moment arms must be carefully designed to preserve 3D muscle geometry lest their accuracy suffer relative to that of anatomically realistic models. The dearth of accurate experimental measurements of 3D moment arms of muscles in birds leaves uncertainty regarding the relative accuracy of different modelling or experimental datasets such as in ostriches. Our model, however, provides a comprehensive set of 3D estimates of muscle actions in ostriches for the first time, emphasizing that avian limb mechanics are highly three-dimensional and complex, and how no muscles act purely in the sagittal plane. A comparative synthesis of experiments and models such as ours could provide powerful synthesis into how anatomy, mechanics and control interact during locomotion and how these interactions evolve. Such a framework could remove obstacles impeding the analysis of muscle function in extinct taxa

Boosting Personalised Musculoskeletal Modelling with Physics-informed Knowledge Transfer

Data-driven methods have become increasingly more prominent for musculoskeletal modelling due to their conceptually intuitive simple and fast implementation. However, the performance of a pre-trained data-driven model using the data from specific subject(s) may be seriously degraded when validated using the data from a new subject, hindering the utility of the personalised musculoskeletal model in clinical applications. This paper develops an active physics-informed deep transfer learning framework to enhance the dynamic tracking capability of the musculoskeletal model on the unseen data. The salient advantages of the proposed framework are twofold: 1) For the generic model, physics-based domain knowledge is embedded into the loss function of the data-driven model as soft constraints to penalise/regularise the data-driven model. 2) For the personalised model, the parameters relating to the feature extraction will be directly inherited from the generic model, and only the parameters relating to the subject-specific inference will be finetuned by jointly minimising the conventional data prediction loss and the modified physics-based loss. In this paper, we use the synchronous muscle forces and joint kinematics prediction from surface electromyogram (sEMG) as the exemplar to illustrate the proposed framework. Moreover, convolutional neural network (CNN) is employed as the deep neural network to implement the proposed framework, and the physics law between muscle forces and joint kinematics is utilised as the soft constraints. Results of comprehensive experiments on a self-collected dataset from eight healthy subjects indicate the effectiveness and great generalization of the proposed framework.Comment: arXiv admin note: text overlap with arXiv:2207.0143

Articular human joint modelling

Copyright @ Cambridge University Press 2009.The work reported in this paper encapsulates the theories and algorithms developed to drive the core analysis modules of the software which has been developed to model a musculoskeletal structure of anatomic joints. Due to local bone surface and contact geometry based joint kinematics, newly developed algorithms make the proposed modeller different from currently available modellers. There are many modellers that are capable of modelling gross human body motion. Nevertheless, none of the available modellers offer complete elements of joint modelling. It appears that joint modelling is an extension of their core analysis capability, which, in every case, appears to be musculoskeletal motion dynamics. It is felt that an analysis framework that is focused on human joints would have significant benefit and potential to be used in many orthopaedic applications. The local mobility of joints has a significant influence in human motion analysis, in understanding of joint loading, tissue behaviour and contact forces. However, in order to develop a bone surface based joint modeller, there are a number of major problems, from tissue idealizations to surface geometry discretization and non-linear motion analysis. This paper presents the following: (a) The physical deformation of biological tissues as linear or non-linear viscoelastic deformation, based on spring-dashpot elements. (b) The linear dynamic multibody modelling, where the linear formulation is established for small motions and is particularly useful for calculating the equilibrium position of the joint. This model can also be used for finding small motion behaviour or loading under static conditions. It also has the potential of quantifying the joint laxity. (c) The non-linear dynamic multibody modelling, where a non-matrix and algorithmic formulation is presented. The approach allows handling complex material and geometrical nonlinearity easily. (d) Shortest path algorithms for calculating soft tissue line of action geometries. The developed algorithms are based on calculating minimum ‘surface mass’ and ‘surface covariance’. An improved version of the ‘surface covariance’ algorithm is described as ‘residual covariance’. The resulting path is used to establish the direction of forces and moments acting on joints. This information is needed for linear or non-linear treatment of the joint motion. (e) The final contribution of the paper is the treatment of the collision. In the virtual world, the difficulty in analysing bodies in motion arises due to body interpenetrations. The collision algorithm proposed in the paper involves finding the shortest projected ray from one body to the other. The projection of the body is determined by the resultant forces acting on it due to soft tissue connections under tension. This enables the calculation of collision condition of non-convex objects accurately. After the initial collision detection, the analysis involves attaching special springs (stiffness only normal to the surfaces) at the ‘potentially colliding points’ and motion of bodies is recalculated. The collision algorithm incorporates the rotation as well as translation. The algorithm continues until the joint equilibrium is achieved. Finally, the results obtained based on the software are compared with experimental results obtained using cadaveric joints

Personalized hip joint kinetics during deep squatting in young, athletic adults

The goal of this study was to report deep squat hip kinetics in young, athletic adults using a personalized numerical model solution based on inverse dynamics. Thirty-five healthy subjects underwent deep squat motion capture acquisitions and MRI scans of the lower extremities. Musculoskeletal models were personalized using each subject's lower limb anatomy. The average peak hip joint reaction force was 274 percent bodyweight. Average peak hip and knee flexion angles were 107 degrees and 112 degrees respectively. These new findings show that deep squatting kinetics in the younger population differ substantially from the previously reported in vivo data in older subjects