Experimental Ultrasound Approach for Studying Knee Intra-Articular Femur-Tibia Movements under Different Loads

The purpose of the present study was to develop an experimental model for the study of intra-articular knee movements depending on the function of the knee joint and involved muscle groups under isometric stretching conditions with different loads. The experimental procedure included an ultrasound examination of a knee joint after isometric stretching in healthy men (n = 32). The changes (in millimeters) in the distances between the femur and tibia were measured using an ultrasound sonographer at three stages. The first stage was performed on ten (n = 10) healthy men in five different sitting and upright positions. In the second and third experimental model stages, lower limbs loading was applied to 22 participants. Our hypothesis, which was confirmed, was that as a result of increased loads on the participant's back, an intra-articular decrease in the femur-tibia cartilage surface distance would be observed. The accuracy of the created experimental model was improved over its three stages from 30% to 9%. Quantitative model data can help to create a mathematical model of the mechanical effects during the deformation of knee joint bone cartilage and it can also help outline some future tasks: increasing loading weights, enlarging participant groups, performing comparisons of men and women, and performing comparisons of healthy and pathological individuals.

Manual wheelchair biomechanics while overcoming various environmental barriers: A systematic review

During manual wheelchair (MWC) locomotion, the user's upper limbs are subject to heavy stresses and fatigue because the upper body is permanently engaged to propel the MWC. These stresses and fatigue vary according to the environmental barriers encountered outdoors along a given path. This study aimed at conducting a systematic review of the literature assessing the biomechanics of MWC users crossing various situations, which represent physical environmental barriers. Through a systematic search on PubMed, 34 articles were selected and classified according to the investigated environmental barriers: slope; cross-slope; curb; and ground type. For each barrier, biomechanical parameters were divided into four categories: spatiotemporal parameters; kinematics; kinetics; and muscle activity. All results from the different studies were gathered, including numerical data, and assessed with respect to the methodology used in each study. This review sheds light on the fact that certain situations (cross-slopes and curbs) or parameters (kinematics) have scarcely been studied, and that a wider set of situations should be studied. Five recommendations were made at the end of this review process to standardize the procedure when reporting materials, methods, and results for the study of biomechanics of any environmental barrier encountered in MWC locomotion: (i) effectively reporting barriers' lengths, grades, or heights; (ii) striving for standardization or a report of the approach conditions of the barrier, such as velocity, especially on curbs; (iii) reporting the configuration of the used MWC, and if it was fitted to the subject's morphology; (iv) reporting rotation sequences for the expression of moments and kinematics, and when used, the definition of the musculoskeletal model; lastly (v) when possible, reporting measurement uncertainties and model reconstruction errors.

Contribution of passive moments to inter-segmental moments during gait: A systematic review

Inter-segmental moments computed by inverse dynamic during gait come from active moments, due to muscle contraction, but also from passive moments, resulting from the resistance of the periarticular structures to their deformation. The evaluation of the proportion of the inter-segmental moments that can be attributed to passive moments has led to divergent results. Thus, the purpose of this study was to systematically search and synthesize the evidence of the contribution of passive moments to inter-segmental moments during healthy and pathological gait. A broad systematic search was performed including four databases. Thirteen studies met all inclusion criteria. Results showed that passive moments participate to inter-segmental moments during gait in a non-negligible way. For the ankle, the evaluation of the proportion of inter-segmental moment attributed to passive structures is 5-20% around the push-off. For the knee, this proportion is 40-98% during late swing and 10-80% during the single support phase. For the hip, it is 20-50% at push-off. For pathological population, it has been shown that this contribution may sometimes be more important, either due to a smaller inter-segmental moment or a larger passive moment. These results suggest that passive mechanisms can contribute substantially to normal human gait, facilitating the propulsion or the braking of the joint. Passive structures, acting as elastic springs, thus help to reduce the energy cost of gait. For pathological gait, studying the contribution of passive moments to inter-segmental moments can help to better understand the aetiology of the pathology.

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.

Multibody Kinematics Optimization for the Estimation of Upper and Lower Limb Human Joint Kinematics: A Systematized Methodological Review

Multibody kinematics optimization (MKO) aims to reduce soft tissue artefact (STA) and is a key step in musculoskeletal modeling. The objective of this review was to identify the numerical methods, their validation and performance for the estimation of the human joint kinematics using MKO. Seventy-four papers were extracted from a systematized search in five databases and cross-referencing. Model-derived kinematics were obtained using either constrained optimization or Kalman filtering to minimize the difference between measured (i.e., by skin markers, electromagnetic or inertial sensors) and model-derived positions and/or orientations. While hinge, universal, and spherical joints prevail, advanced models (e.g., parallel and four-bar mechanisms, elastic joint) have been introduced, mainly for the knee and shoulder joints. Models and methods were evaluated using: (i) simulated data based, however, on oversimplified STA and joint models; (ii) reconstruction residual errors, ranging from 4 mm to 40 mm; (iii) sensitivity analyses which highlighted the effect (up to 36 deg and 12 mm) of model geometrical parameters, joint models, and computational methods; (iv) comparison with other approaches (i.e., single body kinematics optimization and nonoptimized kinematics); (v) repeatability studies that showed low intra- and inter-observer variability; and (vi) validation against ground-truth bone kinematics (with errors between 1 deg and 22 deg for tibiofemoral rotations and between 3 deg and 10 deg for glenohumeral rotations). Moreover, MKO was applied to various movements (e.g., walking, running, arm elevation). Additional validations, especially for the upper limb, should be undertaken and we recommend a more systematic approach for the evaluation of MKO. In addition, further model development, scaling, and personalization methods are required to better estimate the secondary degrees-of-freedom (DoF).

A subject-specific integrative biomechanical framework of the pelvis for gait analysis

Analysis of the human locomotor system using rigid-body musculoskeletal models has increased in the biomechanical community with the objective of studying muscle activations of different movements. Simultaneously, the finite element method has emerged as a complementary approach for analyzing the mechanical behavior of tissues. This study presents an integrative biomechanical framework for gait analysis by linking a musculoskeletal model and a subject-specific finite element model of the pelvis. To investigate its performance, a convergence study was performed and its sensitivity to the use of non-subject-specific material properties was studied. The total hip joint force estimated by the rigid musculoskeletal model and by the finite element model showed good agreement, suggesting that the integrative approach estimates adequately (in shape and magnitude) the hip total contact force. Previous studies found movements of up to 1.4 mm in the anterior-posterior direction, for single leg stance. These results are comparable with the displacement values found in this study: 0-0.5 mm in the sagittal axis. Maximum von Mises stress values of approximately 17 MPa were found in the pelvic bone. Comparing this results with a previous study of our group, the new findings show that the introduction of muscular boundary conditions and the flexion-extension movement of the hip reduce the regions of high stress and distributes more uniformly the stress across the pelvic bone. Thus, it is thought that muscle force has a relevant impact in reducing stresses in pelvic bone during walking of the finite element model proposed in this study. Future work will focus on including other deformable structures, such as the femur and the tibia, and subject-specific material properties.

Contribution of passive actions to the lower limb joint moments and powers during gait: A comparison of models

The lower limb passive actions representing the actions of all the passive periarticular structures have been shown to have a significant contribution to the power generation and absorption during gait. However, the respective magnitude of its different components was not established, although models of ligament moment were implemented in some musculoskeletal models. These ligament moments have shown to have an influence on the musculo-tendon forces and contact forces but the models used were never specifically evaluated, that is, compared to the passive and net joint moments. Two models of passive joint moments and three models of ligament moments were selected from the literature. Ten subjects (23–29 years old, 79.8 ± 9.5 kg, 1.85 ± 0.06 m) participated in the study. Each subject performed three gait cycles in a gait laboratory to acquire the kinematics and ground reaction forces and to compute the ligament, passive and net moments of the right lower limb joints. The contributions of the passive joint moments to the net joint moments were in accordance with the literature, although time shifts appeared for peaks in the hip and knee powers. Two of the models of ligament moments seemed, in fact, to represent the passive joint moments as their contributions were very similar while the third model of ligament moments seemed to represent only penalty-based joint limits. As a conclusion, this study showed that the models of ligament moments existing in the literature do not seem reliable. This study also demonstrated that the use of non-subject-specific models of the passive joint moments could be a valid approach for healthy subjects.

Prediction of In Vivo Knee Joint Loads Using a Global Probabilistic Analysis

Musculoskeletal models are powerful tools that allow biomechanical investigations and predictions of muscle forces not accessible with experiments. A core challenge modelers must confront is validation. Measurements of muscle activity and joint loading are used for qualitative and indirect validation of muscle force predictions. Subject-specific models have reached high levels of complexity and can predict contact loads with surprising accuracy. However, every deterministic musculoskeletal model contains an intrinsic uncertainty due to the high number of parameters not identifiable in vivo. The objective of this work is to test the impact of intrinsic uncertainty in a scaled-generic model on estimates of muscle and joint loads. Uncertainties in marker placement, limb coronal alignment, body segment parameters, Hill-type muscle parameters, and muscle geometry were modeled with a global probabilistic approach (multiple uncertainties included in a single analysis). 5-95% confidence bounds and input/output sensitivities of predicted knee compressive loads and varus/valgus contact moments were estimated for a gait activity of three subjects with telemetric knee implants from the "Grand Challenge Competition." Compressive load predicted for the three subjects showed confidence bounds of 333 ± 248 N, 408 ± 333 N, and 379 ± 244 N when all the sources of uncertainty were included. The measured loads lay inside the predicted 5-95% confidence bounds for 77%, 83%, and 76% of the stance phase. Muscle maximum isometric force, muscle geometry, and marker placement uncertainty most impacted the joint load results. This study demonstrated that identification of these parameters is crucial when subject-specific models are developed.

Influence of the Level of Muscular Redundancy on the Validity of a Musculoskeletal Model

While recent literature has clearly demonstrated that an extensive personalization of the musculoskeletal models was necessary to reach high accuracy, several components of the generic models may be further investigated before defining subject-specific parameters. Among others, the choice in muscular geometry and thus the level of muscular redundancy in the model may have a noticeable influence on the predicted musculotendon and joint contact forces. In this context, the aim of this study was to investigate if the level of muscular redundancy can contribute or not to reduce inaccuracies in tibiofemoral contact forces predictions. For that, the dataset disseminated through the Sixth Grand Challenge Competition to Predict In Vivo Knee Loads was applied to a versatile 3D lower limb musculoskeletal model in which two muscular geometries (i.e., two different levels of muscular redundancy) were implemented. This dataset provides tibiofemoral implant measurements for both medial and lateral compartments and thus allows evaluation of the validity of the model predictions. The results suggest that an increase of the level of muscular redundancy corresponds to a better accuracy of total tibiofemoral contact force whatever the gait pattern investigated. However, the medial and lateral contact forces ratio and accuracy were not necessarily improved when increasing the level of muscular redundancy and may thus be attributed to other parameters such as the location of contact points. To conclude, the muscular geometry, among other components of the generic model, has a noticeable impact on joint contact forces predictions and may thus be correctly chosen even before trying to personalize the model.

In Vivo Knee Contact Force Prediction Using Patient-Specific Musculoskeletal Geometry in a Segment-Based Computational Model

Segment-based musculoskeletal models allow the prediction of muscle, ligament, and joint forces without making assumptions regarding joint degrees-of-freedom (DOF). The dataset published for the “Grand Challenge Competition to Predict in vivo Knee Loads” provides directly measured tibiofemoral contact forces for activities of daily living (ADL). For the Sixth Grand Challenge Competition to Predict in vivo Knee Loads, blinded results for “smooth” and “bouncy” gait trials were predicted using a customized patient-specific musculoskeletal model. For an unblinded comparison, the following modifications were made to improve the predictions: further customizations, including modifications to the knee center of rotation; reductions to the maximum allowable muscle forces to represent known loss of strength in knee arthroplasty patients; and a kinematic constraint to the hip joint to address the sensitivity of the segment-based approach to motion tracking artifact. For validation, the improved model was applied to normal gait, squat, and sit-to-stand for three subjects. Comparisons of the predictions with measured contact forces showed that segment-based musculoskeletal models using patient-specific input data can estimate tibiofemoral contact forces with root mean square errors (RMSEs) of 0.48–0.65 times body weight (BW) for normal gait trials. Comparisons between measured and predicted tibiofemoral contact forces yielded an average coefficient of determination of 0.81 and RMSEs of 0.46–1.01 times BW for squatting and 0.70–0.99 times BW for sit-to-stand tasks. This is comparable to the best validations in the literature using alternative models.

Effect of lower-limb joint models on subject-specific musculoskeletal models and simulations of daily motor activities

Understanding the validity of using musculoskeletal models is critical, making important to assess how model parameters affect predictions. In particular, assumptions on joint models can affect predictions from simulations of movement, and the identification of image-based joints is unavoidably affected by uncertainty that can decrease the benefits of increasing model complexity. We evaluated the effect of different lower-limb joint models on muscle and joint contact forces during four motor tasks, and assessed the sensitivity to the uncertainties in the identification of anatomical four-bar-linkage joints. Three MRI-based musculoskeletal models having different knee and ankle joint models were created and used for the purpose. Model predictions were compared against a baseline model including simpler and widely-adopted joints. In addition, a probabilistic analysis was performed by perturbing four-bar-linkage joint parameters according to their uncertainty. The differences between models depended on the motor task analyzed, and there could be marked differences at peak loading (up to 2.40 BW at the knee and 1.54 BW at the ankle), although they were rather small over the motor task cycles (up to 0.59 BW at the knee and 0.31 BW at the ankle). The model including more degrees of freedom showed more discrepancies in predicted muscle activations compared to measured muscle activity. Further, including image-based four-bar-linkages was robust to simulate walking, chair rise and stair ascent, but not stair descent (peak standard deviation of 2.66 BW), suggesting that joint model complexity should be set according to the imaging dataset available and the intended application, performing sensitivity analyses.

The development of a segment-based musculoskeletal model of the lower limb: introducing FreeBody

Traditional approaches to the biomechanical analysis of movement are joint-based; that is the mechanics of the body are described in terms of the forces and moments acting at the joints, and that muscular forces are considered to create moments about the joints. We have recently shown that segment-based approaches, where the mechanics of the body are described by considering the effect of the muscle, ligament and joint contact forces on the segments themselves, can also prove insightful. We have also previously described a simultaneous, optimization-based, musculoskeletal model of the lower limb. However, this prior model incorporates both joint- and segment-based assumptions. The purpose of this study was therefore to develop an entirely segment-based model of the lower limb and to compare its performance to our previous work. The segment-based model was used to estimate the muscle forces found during vertical jumping, which were in turn compared with the muscular activations that have been found in vertical jumping, by using a Geers' metric to quantify the magnitude and phase errors. The segment-based model was shown to have a similar ability to estimate muscle forces as a model based upon our previous work. In the future, we will evaluate the ability of the segment-based model to be used to provide results with clinical relevance, and compare its performance to joint-based approaches. The segment-based model described in this article is publicly available as a GUI-based Matlab® application and in the original source code (at www.msksoftware.org.uk).

VARIATION OF THE ANKLE MOTION WITH THE PIVOT-POINT POSITION AS PREDICTED BY A SPHERICAL MODEL OF THE JOINT

Recent studies proved that the ankle unloaded motion is nearly spherical and can be accurately modeled by one degree-of-freedom spatial spherical mechanisms. Starting from relevant experimental data measured on nine specimens, the global optimal pivot-point (PP) position, i.e., the center of this spherical motion, is obtained in this study and the effect on ankle motion of the changing position of this PP is analyzed.

Potential of the Pseudo-Inverse Method as a Constrained Static Optimization for Musculo-Tendon Forces Prediction

Inverse dynamics combined with a constrained static optimization analysis has often been proposed to solve the muscular redundancy problem. Typically, the optimization problem consists in a cost function to be minimized and some equality and inequality constraints to be fulfilled. Penalty-based and Lagrange multipliers methods are common optimization methods for the equality constraints management. More recently, the pseudo-inverse method has been introduced in the field of biomechanics. The purpose of this paper is to evaluate the ability and the efficiency of this new method to solve the muscular redundancy problem, by comparing respectively the musculo-tendon forces prediction and its cost-effectiveness against common optimization methods. Since algorithm efficiency and equality constraints fulfillment highly belong to the optimization method, a two-phase procedure is proposed in order to identify and compare the complexity of the cost function, the number of iterations needed to find a solution and the computational time of the penalty-based method, the Lagrange multipliers method and pseudo-inverse method. Using a 2D knee musculo-skeletal model in an isometric context, the study of the cost functions isovalue curves shows that the solution space is 2D with the penalty-based method, 3D with the Lagrange multipliers method and 1D with the pseudo-inverse method. The minimal cost function area (defined as the area corresponding to 5% over the minimal cost) obtained for the pseudo-inverse method is very limited and along the solution space line, whereas the minimal cost function area obtained for other methods are larger or more complex. Moreover, when using a 3D lower limb musculo-skeletal model during a gait cycle simulation, the pseudo-inverse method provides the lowest number of iterations while Lagrange multipliers and pseudo-inverse method have almost the same computational time. The pseudo-inverse method, by providing a better suited cost function and an efficient computational framework, seems to be adapted to the muscular redundancy problem resolution in case of linear equality constraints. Moreover, by reducing the solution space, this method could be a unique opportunity to introduce optimization methods for a posteriori articulation of preference in order to provide a palette of solutions rather than a unique solution based on a lot of hypotheses.