Feasibility of using MRIs to create subject-specific parallel-mechanism joint models

Automatic generation of knee kinematic models from medical imaging

BACKGROUND AND OBJECTIVE

Three-dimensional spatial mechanisms have been used to accurately predict passive knee kinematics, and have shown potential to be used in optimized multibody kinematic models. Such multi-body models are anatomically consistent and can accurately predict passive knee kinematics, but require extensive medical image processing and thus are not widely adopted. This study aimed to automate the generation of kinematic models of tibiofemoral (TFJ) and patellofemoral (PFJ) joints from segmented magnetic resonance imaging (MRI) and compare them against a corresponding manual pipeline.

METHODS

From segmented MRI of eight healthy participants (four females; aged 14.0 ± 2.6 years), geometric parameters (i.e., articular surfaces, ligament attachments) were determined both automatically and manually, and then assembled into TFJ and PFJ kinematic models to predict passive kinematics. The TFJ model was a six-link mechanism with deformable ligamentous constraints, whereas PFJ was a modified hinge. The ligament length changes through TFJ flexion were prescribed to literature strain profile. The geometric parameters were optimized to ensure physiological kinematic predictions through a Multiple Objective Particle Swarm Optimization.

RESULTS

Geometric parameters showed strong agreement between automatic and manual pipelines (median error of 2.8 mm for anatomical landmarks and 1.5 mm for ligament lengths). Predicted TFJ and PFJ kinematics from the two pipelines were not statistically different, except for tibial superior/inferior translation near terminal TFJ extension. The TFJ kinematics predicted from the automatic pipeline had mean errors of 3.6° and 12.4° for adduction/abduction and internal/external rotation, respectively, and <7 mm mean translational error compared to the manual pipeline. Predicted PFJ had <9° mean rotational errors and <6 mm mean translational errors.

CONCLUSIONS

The automatic pipeline developed and presented here can predict passive knee kinematics comparable to a manual pipeline, but removes laborious manual processing and provides a systematic approach to model creation.

The future of in-field sports biomechanics: wearables plus modelling compute real-time in vivo tissue loading to prevent and repair musculoskeletal injuries

This paper explores the use of biomechanics in identifying the mechanistic causes of musculoskeletal tissue injury and degeneration. It appraises how biomechanics has been used to develop training programmes aiming to maintain or recover tissue health. Tissue health depends on the functional mechanical environment experienced by tissues during daily and rehabilitation activities. These environments are the result of the interactions between tissue motion, loading, biology, and morphology. Maintaining health of and/or repairing musculoskeletal tissues requires targeting the "ideal" in vivo tissue mechanics (i.e., loading and deformation), which may be enabled by appropriate real-time biofeedback. Recent research shows that biofeedback technologies may increase their quality and effectiveness by integrating a personalised neuromusculoskeletal modelling driven by real-time motion capture and medical imaging. Model personalisation is crucial in obtaining physically and physiologically valid predictions of tissue biomechanics. Model real-time execution is crucial and achieved by code optimisation and artificial intelligence methods. Furthermore, recent work has also shown that laboratory-based motion capture biomechanical measurements and modelling can be performed outside the laboratory with wearable sensors and artificial intelligence. The next stage is to combine these technologies into well-designed easy to use products to guide training to maintain or recover tissue health in the real-world.

Quantifying changes in shoulder orientation between the prone and supine positions from magnetic resonance imaging

BACKGROUND

Predicting breast tissue motion using biomechanical models can provide navigational guidance during breast cancer treatment procedures. These models typically do not account for changes in posture between procedures. Difference in shoulder position can alter the shape of the pectoral muscles and breast. A greater understanding of the differences in the shoulder orientation between prone and supine could improve the accuracy of breast biomechanical models.

METHODS

19 landmarks were placed on the sternum, clavicle, scapula, and humerus of the shoulder girdle in prone and supine breast MRIs (N = 10). These landmarks were used in an optimization framework to fit subject-specific skeletal models and compare joint angles of the shoulder girdle between these positions.

FINDINGS

The mean Euclidean distance between joint locations from the fitted skeletal model and the manually identified joint locations was 15.7 mm ± 2.7 mm. Significant differences were observed between prone and supine. Compared to supine position, the shoulder girdle in the prone position had the lateral end of the clavicle in more anterior translation (i.e., scapula more protracted) (P < 0.05), the scapula in more protraction (P < 0.01), the scapula in more upward rotation (associated with humerus elevation) (P < 0.05); and the humerus more elevated (P < 0.05) for both the left and right sides.

INTERPRETATION

Shoulder girdle orientation was found to be different between prone and supine. These differences would affect the shape of multiple pectoral muscles, which would affect breast shape and the accuracy of biomechanical models.

Can a Musculoskeletal Model Adapted to Knee Implant Geometry Improve Prediction of 3D Contact Forces and Moments?

Tibiofemoral contact loads are crucial parameters in the onset and progression of osteoarthrosis. While contact loads are frequently estimated from musculoskeletal models, their customization is often limited to scaling musculoskeletal geometry or adapting muscle lines. Moreover, studies have usually focused on superior-inferior contact force without investigating three-dimensional contact loads. Using experimental data from six patients with instrumented total knee arthroplasty (TKA), this study customized a lower limb musculoskeletal model to consider the positioning and the geometry of the implant at knee level. Static optimization was performed to estimate tibiofemoral contact forces and contact moments as well as musculotendinous forces. Predictions from both a generic and a customized model were compared to the instrumented implant measurements. Both models accurately predict superior-inferior (SI) force and abduction-adduction (AA) moment. Notably, the customization improves prediction of medial-lateral (ML) force and flexion-extension (FE) moments. However, there is subject-dependent variability in the prediction of anterior-posterior (AP) force. The customized models presented here predict loads on all joint axes and in most cases improve prediction. Unexpectedly, this improvement was more limited for patients with more rotated implants, suggesting a need for further model adaptations such as muscle wrapping or redefinition of hip and ankle joint centers and axes.

Maintaining soldier musculoskeletal health using personalised digital humans, wearables and/or computer vision

OBJECTIVES

The physical demands of military service place soldiers at risk of musculoskeletal injuries and are major concerns for military capability. This paper outlines the development new training technologies to prevent and manage these injuries.

DESIGN

Narrative review.

METHODS

Technologies suitable for integration into next-generation training devices were examined. We considered the capability of technologies to target tissue level mechanics, provide appropriate real-time feedback, and their useability in-the-field.

RESULTS

Musculoskeletal tissues' health depends on their functional mechanical environment experienced in military activities, training and rehabilitation. These environments result from the interactions between tissue motion, loading, biology, and morphology. Maintaining health of and/or repairing joint tissues requires targeting the "ideal" in vivo tissue mechanics (i.e., loading and strain), which may be enabled by real-time biofeedback. Recent research has shown that these biofeedback technologies are possible by integrating a patient's personalised digital twin and wireless wearable devices. Personalised digital twins are personalised neuromusculoskeletal rigid body and finite element models that work in real-time by code optimisation and artificial intelligence. Model personalisation is crucial in obtaining physically and physiologically valid predictions.

CONCLUSIONS

Recent work has shown that laboratory-quality biomechanical measurements and modelling can be performed outside the laboratory with a small number of wearable sensors or computer vision methods. The next stage is to combine these technologies into well-designed easy to use products.

Multi-level personalization of neuromusculoskeletal models to estimate physiologically plausible knee joint contact forces in children

Neuromusculoskeletal models are a powerful tool to investigate the internal biomechanics of an individual. However, commonly used neuromusculoskeletal models are generated via linear scaling of generic templates derived from elderly adult anatomies and poorly represent a child, let alone children with a neuromuscular disorder whose musculoskeletal structures and muscle activation patterns are profoundly altered. Model personalization can capture abnormalities and appropriately describe the underlying (altered) biomechanics of an individual. In this work, we explored the effect of six different levels of neuromusculoskeletal model personalization on estimates of muscle forces and knee joint contact forces to tease out the importance of model personalization for normal and abnormal musculoskeletal structures and muscle activation patterns. For six children, with and without cerebral palsy, generic scaled models were developed and progressively personalized by (1) tuning and calibrating musculotendon units' parameters, (2) implementing an electromyogram-assisted approach to synthesize muscle activations, and (3) replacing generic anatomies with image-based bony geometries, and physiologically and physically plausible muscle kinematics. Biomechanical simulations of gait were performed in the OpenSim and CEINMS software on ten overground walking trials per participant. A mixed-ANOVA test, with Bonferroni corrections, was conducted to compare all models' estimates. The model with the highest level of personalization produced the most physiologically plausible estimates. Model personalization is crucial to produce physiologically plausible estimates of internal biomechanical quantities. In particular, personalization of musculoskeletal anatomy and muscle activation patterns had the largest effect overall. Increased research efforts are needed to ease the creation of personalized neuromusculoskeletal models.

Development of predictive statistical shape models for paediatric lower limb bones

BACKGROUND AND OBJECTIVE

Accurate representation of bone shape is important for subject-specific musculoskeletal models as it may influence modelling of joint kinematics, kinetics, and muscle dynamics. Statistical shape modelling is a method to estimate bone shape from minimal information, such as anatomical landmarks, and to avoid the time and cost associated with reconstructing bone shapes from comprehensive medical imaging. Statistical shape models (SSM) of lower limb bones have been developed and validated for adult populations but are not applicable to paediatric populations. This study aimed to develop SSM for paediatric lower limb bones and evaluate their reconstruction accuracy using sparse anatomical landmarks.

METHODS

We created three-dimensional models of 56 femurs, 29 pelves, 56 tibias, 56 fibulas, and 56 patellae through segmentation of magnetic resonance images taken from 29 typically developing children (15 females; 13 ± 3.5 years). The SSM for femur, pelvis, tibia, fibula, patella, haunch (i.e., combined femur and pelvis), and shank (i.e., combined tibia and fibula) were generated from manual segmentation of comprehensive magnetic resonance images to describe the shape variance of the cohort. We implemented a leave-one-out cross-validation method wherein SSM were used to reconstruct novel bones (i.e., those not included in SSM generation) using full- (i.e., full segmentation) and sparse- (i.e., anatomical landmarks) input, and then compared these reconstructions against bones segmented from magnetic resonance imaging. Reconstruction performance was evaluated using root mean squared errors (RMSE, mm), Jaccard index (0-1), Dice similarity coefficient (DSC) (0-1), and Hausdorff distance (mm). All results reported in this abstract are mean ± standard deviation.

RESULTS

Femurs, pelves, tibias, fibulas, and patellae reconstructed via SSM using full-input had RMSE between 0.89 ± 0.10 mm (patella) and 1.98 ± 0.38 mm (pelvis), Jaccard indices between 0.77 ± 0.03 (pelvis) and 0.90 ± 0.02 (tibia), DSC between 0.87 ± 0.02 (pelvis) and 0.95 ± 0.01 (tibia), and Hausdorff distances between 2.45 ± 0.57 mm (patella) and 9.01 ± 2.36 mm (pelvis). Reconstruction using sparse-input had RMSE ranging from 1.33 ± 0.61 mm (patella) to 3.60 ± 1.05 mm (pelvis), Jaccard indices ranging from 0.59 ± 0.10 (pelvis) to 0.83 ± 0.03 (tibia), DSC ranging from 0.74 ± 0.08 (pelvis) to 0.90 ± 0.02 (tibia), and Hausdorff distances ranging from 3.21 ± 1.19 mm (patella) to 12.85 ± 3.24 mm (pelvis).

CONCLUSIONS

The SSM of paediatric lower limb bones showed reconstruction accuracy consistent with previously developed SSM and outperformed adult-based SSM when used to reconstruct paediatric bones.

[Evaluation of a Bone Coordinate System Constructed Using MR Image Composing]

PURPOSE

To evaluate the accuracy of a bone coordinate system constructed using MR image composing.

METHOD

A femoral coordinate system constructed using image composing of MR images of a whole bovine femur was evaluated using CT images. The MR images were acquired by moving the table and were processed with 3D distortion correction and composing. To evaluate the reproducibility of the measurements, the same operator repeated the construction of the femoral coordinate system. In addition, distortions in the MR images were evaluated in comparison with those in the CT images.

RESULT

The center position of the femoral coordinate system constructed using the MR image composing was 1.6±0.9 mm on the X-axis, 1.5±0.8 mm on the Y-axis, and 0.2±0.3 mm on the Z-axis, and the rotation of each axis was 1° or less. The distortion of the composed MR image was about 0.3%.

CONCLUSION

The femoral coordinate system constructed using MR image composing had the same accuracy as a system constructed with CT images. The effect of MR image composing on the construction of the femoral coordinate system was small.

Biomechanical Study Concerning the Types of Resection in Arthrodesis of First Metatarsophalangeal Joint

This work concerns a biomechanical study aiming to ascertain the optimal type of joint resection when performing a joint arthrodesis. A 3-dimensional digital model of the first metatarsophalangeal joint including the entire first metatarsal bone and proximal phalanx using CT scans of the forefoot was created. Using this computer model, 4 types of resections; ball-and-socket, flat-on-flat, wedge 90°, and wedge 100° were simulated. Parameters measured using this model were the force necessary to separate the 2 fused surfaces, the surface area of the resected surfaces and the shortening of the first ray. By measuring the reactive force necessary to separate the phalanx from the first metatarsal, the 90° wedge resection was found to be the most stable, with comparable results in the case of the 100° wedge resection. Wedge resections are also more favorable when comparing the shortening of the first ray. Wedge resections, though being more technically difficult to perform prove to be the most stable for metatarsophalangeal joint-1 arthrodesis using this model.

A Novel Kinematic Model of the Tibiofemoral Joint Based on a Parallel Mechanism

This paper presents a complete kinematic model of the tibiofemoral joint (TFJ) based on a RRPP + 4-SPS parallel mechanism, where R, P, and S stand for revolute, prismatic, and spherical joints, respectively. The model accounts for the contact between tibia and femur, and the four major ligaments: anterior cruciate, posterior cruciate, medial collateral, and lateral collateral, with anatomical significance in their length variations. An experimental flexion passive motion task is performed, and the kinematic model is tested to determine its capability to reproduce the workspace of the motion task. In addition, an optimization process is performed to simulate prescribed ligament length variations during the motion task. The proposed kinematic model is capable to reproduce with high accuracy an experimental three-dimensional workspace, and at the same time, to simulate prescribed ligament length variation during the spatial flexion task. Prescribed ligament length variations are achieved through an optimization process of the ligament insertion points. This model can be used to improve the multibody kinematic optimization (MKO) process during gait analysis, and also in the design of rehabilitation devices as well as trajectories to accelerate the recovery of injured ligaments. The model shows potential to predict ligament length variations during different motion tasks, and can serve as a basis to develop complex models for kinetostatic and dynamic analyses without dealing with computationally expensive models.

Automatic generation of personalised skeletal models of the lower limb from three-dimensional bone geometries

The generation of personalised and patient-specific musculoskeletal models is currently a cumbersome and time-consuming task that normally requires several processing hours and trained operators. We believe that this aspect discourages the use of computational models even when appropriate data are available and personalised biomechanical analysis would be beneficial. In this paper we present a computational tool that enables the fully automatic generation of skeletal models of the lower limb from three-dimensional bone geometries, normally obtained by segmentation of medical images. This tool was evaluated against four manually created lower limb models finding remarkable agreement in the computed joint parameters, well within human operator repeatability. The coordinate systems origins were identified with maximum differences between 0.5 mm (hip joint) and 5.9 mm (subtalar joint), while the joint axes presented discrepancies between 1° (knee joint) to 11° (subtalar joint). To prove the robustness of the methodology, the models were built from four datasets including both genders, anatomies ranging from juvenile to elderly and bone geometries reconstructed from high-quality computed tomography as well as lower-quality magnetic resonance imaging scans. The entire workflow, implemented in MATLAB scripting language, executed in seconds and required no operator intervention, creating lower extremity models ready to use for kinematic and kinetic analysis or as baselines for more advanced musculoskeletal modelling approaches, of which we provide some practical examples. We auspicate that this technical advancement, together with upcoming progress in medical image segmentation techniques, will promote the use of personalised models in larger-scale studies than those hitherto undertaken.

Automated creation and tuning of personalised muscle paths for OpenSim musculoskeletal models of the knee joint

Computational modelling is an invaluable tool for investigating features of human locomotion and motor control which cannot be measured except through invasive techniques. Recent research has focussed on creating personalised musculoskeletal models using population-based morphing or directly from medical imaging. Although progress has been made, robust definition of two critical model parameters remains challenging: (1) complete tibiofemoral (TF) and patellofemoral (PF) joint motions, and (2) muscle tendon unit (MTU) pathways and kinematics (i.e. lengths and moment arms). The aim of this study was to develop an automated framework, using population-based morphing approaches to create personalised musculoskeletal models, consisting of personalised bone geometries, TF and PF joint mechanisms, and MTU pathways and kinematics. Informed from medical imaging, personalised rigid body TF and PF joint mechanisms were created. Using atlas- and optimisation-based methods, personalised MTU pathways and kinematics were created with the aim of preventing MTU penetration into bones and achieving smooth MTU kinematics that follow patterns from existing literature. This framework was integrated into the Musculoskeletal Atlas Project Client software package to create and optimise models for 6 participants with incrementally increasing levels of personalisation with the aim of improving MTU kinematics and pathways. Three comparisons were made: (1) non-optimised (Model 1) and optimised models (Model 3) with generic joint mechanisms; (2) non-optimised (Model 2) and optimised models (Model 4) with personalised joint mechanisms; and (3) both optimised models (Model 3 and 4). Following optimisation, improvements were consistently shown in pattern similarity to cadaveric data in comparison (1) and (2). For comparison (3), a number of comparisons showed no significant difference between the two compared models. Importantly, optimisation did not produce statistically significantly worse results in any case.

The relationship between tibiofemoral geometry and musculoskeletal function during normal activity

BACKGROUND

The effect of tibiofemoral geometry on musculoskeletal function is important to movement biomechanics.

RESEARCH QUESTION

We hypothesised that tibiofemoral geometry determines tibiofemoral motion and musculoskeletal function. We then aimed at 1) modelling tibiofemoral motion during normal activity as a function of tibiofemoral geometry in healthy adults; and 2) quantifying the effect of tibiofemoral geometry on musculoskeletal function.

METHODS

We used motion data for six activity types and CT images of the knee from 12 healthy adults. Geometrical variation of the tibia and femoral articular surfaces were measured in the CT images. The geometry-based tibiofemoral motion was calculated by fitting a parallel mechanism to geometrical variation in the cohort. Matched musculoskeletal models embedding the geometry-based tibiofemoral joint motion and a common generic tibiofemoral motion of reference were generated and used to calculate joint angles, net joint moments, muscle and joint forces for the six activities analysed. The tibiofemoral model was validated against bi-planar fluoroscopy measurements for walking for all the six planes of motion. The effect of tibiofemoral geometry on musculoskeletal function was the difference between the geometry-based model and the model of reference.

RESULTS

The geometry-based tibiofemoral motion described the pattern and the variation during walking for all six motion components, except the pattern of anterior tibial translation. Tibiofemoral geometry had moderate effect on cohort-averages of musculoskeletal function (R² = 0.60-1), although its effect was high in specific instances of the model, outputs and activities analysed, reaching 2.94 BW for the ankle reaction force during stair descent. In conclusion, tibiofemoral geometry is a major determinant of tibiofemoral motion during walking.

SIGNIFICANCE

Geometrical variations of the tibiofemoral joint are important for studying musculoskeletal function during normal activity in specific individuals but not for studying cohort averages of musculoskeletal function. This finding expands current knowledge of movement biomechanics.

An Anatomical-Based Subject-Specific Model of In-Vivo Knee Joint 3D Kinematics From Medical Imaging

Biomechanical models of the knee joint allow the development of accurate procedures as well as novel devices to restore the joint natural motion. They are also used within musculoskeletal models to perform clinical gait analysis on patients. Among relevant knee models in the literature, the anatomy-based spatial parallel mechanisms represent the joint motion using rigid links for the ligaments’ isometric fibres and point contacts for the articular surfaces. To customize analyses, therapies and devices, there is the need to define subject-specific models, but relevant procedures and their accuracy are still questioned. A procedure is here proposed and validated to define a customized knee model based on a spatial parallel mechanism. Computed tomography, magnetic resonance and 3D-video-fluoroscopy were performed on a healthy volunteer to define the personalized model geometry. The model was then validated by comparing the measured and the replicated joint motion. The model showed mean absolute difference and standard deviations in translations and rotations, respectively of 0.98 ± 0.40 mm and 0.68 ± 0.29 ° for the tibia–femur motion, and of 0.77 ± 0.15 mm and 2.09 ± 0.69 ° for the patella–femur motion. These results show that accurate personalized spatial models of knee kinematics can be obtained from in-vivo imaging.

Increasing level of neuromusculoskeletal model personalisation to investigate joint contact forces in cerebral palsy: A twin case study

BACKGROUND

Cerebral palsy is a complex neuromuscular disorder that affects the sufferers in multiple different ways. Neuromusculoskeletal models are promising tools that can be used to plan patient-specific treatments for cerebral palsy. However, current neuromusculoskeletal models are typically scaled from generic adult templates that poorly represent paediatric populations. Furthermore, muscle activations are commonly computed via optimisation methods, which may not reproduce co-contraction observed in cerebral palsy. Alternatively, calibrated EMG-informed approaches within OpenSim can capture pathology-related muscle activation abnormalities, possibly enabling more feasible estimations of muscle and joint contact forces.

METHODS

Two identical twin brothers, aged 13, one with unilateral cerebral palsy and the other typically developing, were enrolled in the study. Four neuromusculoskeletal models with increasing subject-specificity were built in OpenSim and CEINMS combining literature findings, experimental motion capture, EMG and MR data for both participants. The physiological and biomechanical validity of each model was assessed by quantifying its ability to track experimental joint moments and muscle excitations.

FINDINGS

All developed models accurately tracked external joint moments; however EMG-informed models better tracked muscle excitations compared to neural solutions generated by static optimisation. Calibrating muscle-tendon unit parameters with EMG data allowed for more physiologically plausible joint contact forces estimates. Further scaling the maximal isometric force of muscles with MR-derived muscle volumes did not affect model predictions.

INTERPRETATION

Given their ability to identify atypical joint contact forces profiles and accurately reproduce experimental data, calibrated EMG-informed models should be preferred over generic models using optimisation methods in informing the management of cerebral palsy.

Evaluation of predicted patellofemoral joint kinematics with a moving-axis joint model

The main objectives of this study were to expand the moving-axis joint model concept to the patellofemoral joint and evaluate the patellar motion against experimental patellofemoral kinematics. The experimental data was obtained through 2D-to-3D bone reconstruction of EOS images and segmented MRI data utilizing an iterative closest point optimization technique. Six knee model variations were developed using the AnyBody Modeling System and subject-specific bone geometries. These models consisted of various combinations of tibiofemoral (hinge, moving-axis, and interpolated) and patellofemoral (hinge and moving-axis) joint types. The newly introduced interpolated tibiofemoral joint is calibrated from the five EOS quasi-static lunge positions. The patellofemoral axis of the hinge model was defined by performing surface fits to the patellofemoral contact area; and the moving-axis model was defined based upon the position of the patellofemoral joint at 0° and 90° tibiofemoral-flexion. In between these angles, the patellofemoral axis moved linearly as a function of tibiofemoral-flexion, while outside these angles, the axis remained fixed. When using a moving-axis tibiofemoral joint, a hinge patellofemoral joint offers (-5.12 ± 1.23 mm, 5.81 ± 0.97 mm, 14.98 ± 2.30°, -4.35 ± 1.95°) mean differences (compared to EOS) while a moving-axis patellofemoral model provides (-2.69 ± 1.04 mm, 1.13 ± 0.80 mm, 12.63 ± 2.03°, 1.74 ± 1.46°) in terms of lateral-shift, superior translation, patellofemoral-flexion, and patellar-rotation, respectively. Furthermore, the model predictive capabilities increased as a direct result of adding more calibrated positions to the tibiofemoral model (hinge-1, moving-axis-2, and interpolated-5). Overall, a novel subject-specific moving-axis patellofemoral model has been established; that produces realistic patellar motion and is computationally fast enough for clinical applications.

A screening method to analyse the sensitivity of a lower limb multibody kinematic model

The study presents a screening method used to identify the influential parameters of a lower limb model including ligaments, at low numerical cost. Concerning multibody kinematics optimisation, the ligament parameters (isometric length) were found the most influential ones in a previous study. The screening method tested if they remain influential with minimised length variations. The most important parameters for tibiofemoral kinematics were the skin markers, segment lengths and joint parameters, including two ligaments. This was confirmed by a quantitative sensitivity analysis. The screening method has the potential to be used as a stand-alone procedure for a sensitivity analysis.

Minimal medical imaging can accurately reconstruct geometric bone models for musculoskeletal models

Accurate representation of subject-specific bone anatomy in lower-limb musculoskeletal models is important for human movement analyses and simulations. Mathematical methods can reconstruct geometric bone models using incomplete imaging of bone by morphing bone model templates, but the validity of these methods has not been fully explored. The purpose of this study was to determine the minimal imaging requirements for accurate reconstruction of geometric bone models. Complete geometric pelvis and femur models of 14 healthy adults were reconstructed from magnetic resonance imaging through segmentation. From each complete bone segmentation, three sets of incomplete segmentations (set 1 being the most incomplete) were created to test the effect of imaging incompleteness on reconstruction accuracy. Geometric bone models were reconstructed from complete sets, three incomplete sets, and two motion capture-based methods. Reconstructions from (in)complete sets were generated using statistical shape modelling, followed by host-mesh and local-mesh fitting through the Musculoskeletal Atlas Project Client. Reconstructions from motion capture-based methods used positional data from skin surface markers placed atop anatomic landmarks and estimated joint centre locations as target points for statistical shape modelling and linear scaling. Accuracy was evaluated with distance error (mm) and overlapping volume similarity (%) between complete bone segmentation and reconstructed bone models, and statistically compared using a repeated measure analysis of variance (p<0.05). Motion capture-based methods produced significantly higher distance error than reconstructions from (in)complete sets. Pelvis volume similarity reduced significantly with the level of incompleteness: complete set (92.70±1.92%), set 3 (85.41±1.99%), set 2 (81.22±3.03%), set 1 (62.30±6.17%), motion capture-based statistical shape modelling (41.18±9.54%), and motion capture-based linear scaling (26.80±7.19%). A similar trend was observed for femur volume similarity. Results indicate that imaging two relevant bone regions produces overlapping volume similarity >80% compared to complete segmented bone models, and improve analyses and simulation over current standard practice of linear scaling musculoskeletal models.

Knee medial and lateral contact forces in a musculoskeletal model with subject-specific contact point trajectories

Contact point (CP) trajectory is a crucial parameter in estimating medial/lateral tibio-femoral contact forces from the musculoskeletal (MSK) models. The objective of the present study was to develop a method to incorporate the subject-specific CP trajectories into the MSK model. Ten healthy subjects performed 45 s treadmill gait trials. The subject-specific CP trajectories were constructed on the tibia and femur as a function of extension-flexion using low-dose bi-plane X-ray images during a quasi-static squat. At each extension-flexion position, the tibia and femur CPs were superimposed in the three directions on the medial side, and in the anterior-posterior and proximal-distal directions on the lateral side to form the five kinematic constraints of the knee joint. The Lagrange multipliers associated to these constraints directly yielded the medial/lateral contact forces. The results from the personalized CP trajectory model were compared against the linear CP trajectory and sphere-on-plane CP trajectory models which were adapted from the commonly used MSK models. Changing the CP trajectory had a remarkable impact on the knee kinematics and changed the medial and lateral contact forces by 1.03 BW and 0.65 BW respectively, in certain subjects. The direction and magnitude of the medial/lateral contact force were highly variable among the subjects and the medial-lateral shift of the CPs alone could not determine the increase/decrease pattern of the contact forces. The suggested kinematic constraints are adaptable to the CP trajectories derived from a variety of joint models and those experimentally measured from the 3D imaging techniques.

Bioinspired Technologies to Connect Musculoskeletal Mechanobiology to the Person for Training and Rehabilitation

Musculoskeletal tissues respond to optimal mechanical signals (e.g., strains) through anabolic adaptations, while mechanical signals above and below optimal levels cause tissue catabolism. If an individual's physical behavior could be altered to generate optimal mechanical signaling to musculoskeletal tissues, then targeted strengthening and/or repair would be possible. We propose new bioinspired technologies to provide real-time biofeedback of relevant mechanical signals to guide training and rehabilitation. In this review we provide a description of how wearable devices may be used in conjunction with computational rigid-body and continuum models of musculoskeletal tissues to produce real-time estimates of localized tissue stresses and strains. It is proposed that these bioinspired technologies will facilitate a new approach to physical training that promotes tissue strengthening and/or repair through optimal tissue loading.

Biofeedback for Gait Retraining Based on Real-Time Estimation of Tibiofemoral Joint Contact Forces

Biofeedback assisted rehabilitation and intervention technologies have the potential to modify clinically relevant biomechanics. Gait retraining has been used to reduce the knee adduction moment, a surrogate of medial tibiofemoral joint loading often used in knee osteoarthritis research. In this paper, we present an electromyogram-driven neuromusculoskeletal model of the lower-limb to estimate, in real-time, the tibiofemoral joint loads. The model included 34 musculotendon units spanning the hip, knee, and ankle joints. Full-body inverse kinematics, inverse dynamics, and musculotendon kinematics were solved in real-time from motion capture and force plate data to estimate the knee medial tibiofemoral contact force (MTFF). We analyzed five healthy subjects while they were walking on an instrumented treadmill with visual biofeedback of their MTFF. Each subject was asked to modify their gait in order to vary the magnitude of their MTFF. All subjects were able to increase their MTFF, whereas only three subjects could decrease it, and only after receiving verbal suggestions about possible gait modification strategies. Results indicate the important role of knee muscle activation patterns in modulating the MTFF. While this paper focused on the knee, the technology can be extended to examine the musculoskeletal tissue loads at different sites of the human body.