Erratum: Uncertainty in muscle-tendon parameters can greatly influence the accuracy of knee contact force estimates of musculoskeletal models
- 1Laboratory for Movement Biomechanics, ETH Zürich, Zürich, Switzerland
- 2Aesculap AG, Tuttlingen, Germany
- 3Department of Orthopaedic and Trauma Surgery, Ludwig Maximilians University Munich, Musculoskeletal University Center Munich (MUM), Campus Grosshadern, Munich, Germany
- 4Julius Wolff Institute, Berlin Institute of Health at Charité—Universitätsmedizin Berlin, Berlin, Germany
Understanding the sources of error is critical before models of the musculoskeletal system can be usefully translated. Using in vivo measured tibiofemoral forces, the impact of uncertainty in muscle–tendon parameters on the accuracy of knee contact force estimates of a generic musculoskeletal model was investigated following a probabilistic approach. Population variability was introduced to the routine musculoskeletal modeling framework by perturbing input parameters of the lower limb muscles around their baseline values. Using ground reaction force and skin marker trajectory data collected from six subjects performing body-weight squat, the knee contact force was calculated for the perturbed models. The combined impact of input uncertainties resulted in a considerable variation in the knee contact force estimates (up to 2.1 BW change in the predicted force), especially at larger knee flexion angles, hence explaining up to 70% of the simulation error. Although individual muscle groups exhibited different contributions to the overall error, variation in the maximum isometric force and pathway of the muscles showed the highest impacts on the model outcomes. Importantly, this study highlights parameters that should be personalized in order to achieve the best possible predictions when using generic musculoskeletal models for activities involving deep knee flexion.
Introduction
Musculoskeletal models have been widely used to estimate in vivo loading conditions within the knee (Richards and Higginson, 2010; Worsley et al., 2011; Trepczynski et al., 2012; Gerus et al., 2013; Stylianou et al., 2013; Chen et al., 2016; Trepczynski et al., 2018; van Rossom et al., 2018; Imani Nejad et al., 2020). Outputs of musculoskeletal simulations can be used to predict postoperative functional outcomes of different surgeries (Barry et al., 2010; Chen et al., 2016), optimize rehabilitation protocols (Barry et al., 2010; Jamwal et al., 2020; Li et al., 2020), and enhance athletic performance (Heron, 2015; Langholz et al., 2016; Ataei et al., 2020; Seow et al., 2020). However, when musculoskeletal predictions of knee loads are compared against in vivo measurements, substantial errors are common (Lundberg et al., 2012; Valente et al., 2014; Charles et al., 2020; Koller et al., 2021), especially when generic models are used (e.g., errors of up to 150% for body-weight squat (Schellenberg et al., 2018; Imani Nejad et al., 2020)).
Generic musculoskeletal models are generally developed based on averaged anatomical data obtained from cadavers or living subjects (Delp et al., 1990; Arnold et al., 2010; Rajagopal et al., 2016). To represent a specific subject, certain model parameters are usually scaled based on the individual’s weight and skin marker locations in a static pose (Kainz et al., 2017; Imani Nejad et al., 2020; Ma et al., 2021). However, the scaled parameters possess inherent uncertainty due to marker placement inaccuracy, on top of considerable inter-subject variability due to skeletal morphology, for example, femoral anteversion angle (Kainz et al., 2021; Modenese et al., 2021) and muscle parameters, for example, tendon slack length (Winby et al., 2008; Heinen et al., 2016). As a result, the importance of using subject-specific muscle parameters to reduce uncertainty in modeling outcomes has been consistently emphasized in the literature (Barry et al., 2010; Gerus et al., 2013; Myers et al., 2015; Chen et al., 2016; Koller et al., 2021), but such parameters are generally difficult to measure. For example, to determine the exact muscle pathway, multiple MR images should be captured along the length of the muscle and preferably during different body poses (Schmid et al., 2009; Fernandez et al., 2016). Muscle force generation parameters also play an important role: the maximum isometric force (MIF) of muscles can be estimated during maximum voluntary contraction trials using either force or torque sensors (Be Groote et al., 2010; Ivanovic and Dopsaj, 2013; Kainz et al., 2018; Higa et al., 2019) or electromyography (EMG) data (Pizzolato et al., 2015; Hoang et al., 2019); however, these assessments are generally not sufficiently repeatable or sensitive to determine the MIF of a single muscle (Bolliger et al., 2008). It is also known that EMG signals measured during some dynamic activities can exceed the values measured during maximum voluntary contraction trials (Burden, 2010), thus the measured MIF may not accurately represent the real muscle strength. As a result, measurement and implementation of all subject-specific muscle parameters for all muscles modeled is extremely difficult, expensive, and consequently impracticable. Therefore, it is crucial to limit subject-specific measurements to the muscles and parameters that have a considerable impact on the modeling outputs.
A probabilistic analysis can provide important insights into the relative impact of different sources of uncertainty on the muscle and joint reaction force estimates of generic musculoskeletal models (Barry et al., 2010; Valente et al., 2014; Myers et al., 2015; Lamberto et al., 2017; Zuk et al., 2018). Using experimental gait data collected from a single subject as inputs to a generic musculoskeletal model, Myers et al. (2015) found that measurement inaccuracies including marker placement and skin movement artifact have a relatively small contribution toward the overall uncertainty in muscle force estimates for level walking (e.g., 20 N variation in estimated gastrocnemius force due to marker placement errors, compared to 80 N variation due to errors in muscle input parameters). Nevertheless, they found 1.7 times greater impact on muscle force outputs due to uncertainty in muscle parameters, particularly maximum isometric force and tendon slack length. Similar findings were reported by Navacchia et al. (2016), who used a global probabilistic modeling approach to assess uncertainty in knee contact forces (KCFs). Although these investigations provide an underlying understanding of how uncertainties propagate through different stages of musculoskeletal simulations, they generally focused on level walking with experimental data collected from a single subject and without reliable in vivo data for validation. As a result, to date, there is no comprehensive probabilistic analysis that provides a reference for guiding the selection of influential model parameters that need to be personalized, especially for activities involving deep knee flexion.
In this study, a probabilistic musculoskeletal modeling approach was exploited to quantify the impact of uncertainty in lower limb muscle parameters on KCF estimates of a generic musculoskeletal model used to simulate squat activity. We used skin marker trajectories, ground reaction forces (GRF), and in vivo measured tibiofemoral joint contact forces reported within the CAMS-Knee data sets (Taylor et al., 2017) to perform the simulation and validate the results. The findings of this study aim to provide an improved understanding of the influence of model uncertainty on musculoskeletal simulation outcomes, as well as highlight parameters that should be personalized in order to achieve the best possible predictions of knee contact force.
Methods
Skin marker trajectories, ground reaction forces, and in vivo KCF data, measured in six subjects (74 ± 6 years, 89 ± 13 kg, and 172 ± 4 cm) with an instrumented knee implant (INNEX FIXUC, Zimmer, Switzerland) were obtained from the CAMS-Knee data sets (Taylor et al., 2017).
A generic musculoskeletal model (Rajagopal et al., 2016) developed within the OpenSim modeling environment (Delp et al., 1990) with 37 degrees of freedom (DoFs) and 80 muscle–tendon units actuating the lower limb was selected for this study. To obtain the so-called baseline models representing each subject’s anthropometry, the position of selected bone landmarks and anatomic joint centers (3D position of these landmarks are extracted from the segmented CT images and available within the CAMS-Knee data sets) were used to scale each bone in the lower limb. More specifically, the 3D distances between six bone landmarks on the pelvis (left and right anterior and left and right posterior iliac spine as well as the two hip joint centers) were calculated based on the data obtained from CT images and compared against the distance between corresponding landmarks on the generic bone to scale the pelvis in three dimensions. The distance between proximal and distal joint centers was used to scale the lengths of the tibiae and femora. An average scale factor was calculated based on the distance between the medial and lateral epicondyles as well as the distance between the medial and lateral malleoli to scale the width and thickness of each femur, tibia, and patella. The feet and upper limb segments were scaled based on skin marker positions in a static pose. Optimal fiber length and tendon slack length were scaled based on the standard methods built into the OpenSim Scale Tool, while maximum isometric forces of the muscles were taken directly from the generic model.
Muscles in the generic model were represented by an elastic tendon Hill-type muscle model (Millard2012EquilibriumMuscle (Millard et al., 2013)) with MIF, tendon slack length (TSL), and pennation angle (PEN) obtained from cadavers (Ward et al., 2009; Handsfield et al., 2014). The pathway of each muscle–tendon unit was specified using its origin and insertion points (OIPs), as well as via points (VIA) and/or wrapping objects, which were additionally used to represent the curvilinear path of the muscles. Using the OpenSim/MATLAB application programming interface, values for MIF, PEN, TSL, OIP, and VIA of seven major groups of the lower limb muscles (Table 1) were perturbed around their baseline values to perform a series of Monte Carlo (MC) simulations (Figure 1). Here, in addition to a general MC analysis where all muscles of the lower limb were perturbed simultaneously, individual MC analyses were performed to assess the relative contribution of the knee (flexors and extensors), ankle, and hip muscles toward the overall uncertainty in KCF estimates. In addition, the biarticular muscles were investigated separately to understand how much uncertainty propagates from neighboring joints to the knee.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/808027/fbioe-10-808027-HTML-r2/image_m/fbioe-10-808027-g001.jpg)
FIGURE 1. Probabilistic modeling flowchart. GRF, ground reaction force; IK, inverse kinematics; SO, static optimization; JRF, joint reaction force; KCF, knee contact force.
The muscle–tendon model parameters were probabilistically represented as Gaussian distributions based on the inter-individual variability reported in the literature (Friederich and Brand, 1990; Fukunaga et al., 1996; Lloyd and Besier, 2003; Ward et al., 2009), Supplementary Table S1). In case different coefficients of variation were reported for the same parameter, the weighted average coefficient was used, where the weight was determined based on the number of studied subjects. If the parameter variability for a specific muscle was not available in the literature, the mean of the reported coefficients of variation for other muscles was used to sample random input parameters. Regarding lack of information about the range of variability in muscle pathway parameters, similar to a previous probabilistic study (Navacchia et al., 2016), 5 mm was used as the standard deviation for OIP and VIA distributions.
To assess the overall impact of uncertainty in muscle model parameters on predicted muscle and knee contact forces, we performed one general MC simulation with 2000 iterations for each subject and muscle group, where all the five input parameters of the included muscles were perturbed. In addition, to assess the isolated effect of each source of uncertainty, five individual MC simulations each with 500 iterations were performed where only a single muscle parameter was perturbed at a time. The number of iterations was chosen based on convergence at approximately 1,600 and 350 simulations during early model testing. Here, similar to previous probabilistic musculoskeletal modeling studies (Valente et al., 2014; Myers et al., 2015), convergence of the Monte Carlo simulations was ensured when adding 10% more iterations resulted in less than 1% change in the mean confidence bound of the outcome KCF. In order to ensure that convergence could be reasonably reached for all models within the analysis pipeline, higher values of 2000 and 500 iterations were chosen. During post-processing, all simulations were checked and confirmed that convergence had indeed been reached within these provided analysis conditions. As a result, the designed MC had a total number of 189,000 iterations.
For each subject, skin marker trajectory and GRF data of a single representative squat trial were input to the standard OpenSim tools. Inverse kinematics (IK), static optimization (SO), and joint reaction force (JRF) analysis tools were used to calculate muscle activation and total KCFs in the baseline as well as in the perturbed models by minimizing the sum of the squared muscle activations. To assess the influence of using more simplified muscle models, baseline model simulations were repeated without muscle–tendon force–length characteristics in the static optimization tool.
The measured EMG signals were processed within MATLAB (R2017b, MathWorks, United States) where the raw signal was bandpass filtered (fourth order Butterworth, lowpass 10 Hz, and highpass 300 Hz), offset corrected, rectified, and finally lowpass filtered using a moving average filter. The maximum activation value recorded across all trials of all activities performed for each subject was used to normalize the processed EMG signals. Regarding measurement errors/artefacts and uncertainty introduced by normalization technique, EMG signals were only used for qualitative validation of the muscle activation patterns estimated by the baseline models. In addition, the 5th–95th percentile range of the muscle activation outputs was used to understand how model parameter uncertainty influenced muscle redundancy solutions. The calculated 5th–95th percentile range of the KCF error distributions obtained from general MC was compared between different muscle groups to identify the muscle group with the highest contribution toward the overall KCF uncertainty. In addition, for each muscle group, the ratio between 5th–95th percentile range of individual and general MC outputs was calculated to assess the relative contribution of each muscle parameter to the overall uncertainty observed in the KCFs. Finally, perturbed models with the least root mean squared error (RMSE) between predicted and measured KCFs were selected to understand whether subject-specific measurement of all muscle parameters can guarantee fully accurate KCF estimates when squat activity is simulated using generic models.
Results
Activation patterns of the knee extensor muscles predicted by the baseline models displayed similar trends to the measured EMG signals, especially during the descending phase of the squat (Figure 2). The flexion-dependent increase in EMG patterns of the rectus femoris and vastii was also present in the simulation outcomes. However, while the simulation results showed minimal activation of almost all muscles at the end of the squat cycle, the measured EMG signals indicated some residual co-contraction of the extensor and flexor muscles. The estimated muscle activation for the biceps femoris long head displayed a distinct peak at deep knee flexion, where the EMG signals showed minimal activation of this muscle. The models predicted very minimal activation of the gastrocnemii throughout the entire squat cycle, which was also observed in the EMG signals. However in some subjects, EMG sensors measured a relatively high activation in these muscles at the beginning or end of the activity cycle.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/808027/fbioe-10-808027-HTML-r2/image_m/fbioe-10-808027-g002.jpg)
FIGURE 2. In vivo measured EMG (dotted lines) vs. predicted muscle activation levels for selected lower limb muscles during the studied squat trials. Solid lines represent average (baseline) activation patterns obtained from baseline models, whereas shaded areas display 5th–95th percentile range of the muscle activations obtained from general Monte Carlo simulations. The downward and upward arrows represent the descending and ascending phases of the squat activity.
Simultaneous perturbation of the MIF, PEN, TSL, OIP, and VIA parameters of the lower limb muscles resulted in variations of up to 80% in the activation levels of the rectus femoris, vastus medialis, and biceps femoris long head, as well as up to 50% variation in activation of the vastus lateralis and semimembranosus (Figure 2). These variations were generally subject-specific and flexion-dependent. However, the impact of uncertainty in input parameters on the activation of tibialis anterior and gastrocnemii muscles was generally very small (Figure 2).
As required to balance the external knee moment, predicted muscle forces from the baseline models generally showed the greatest forces in the knee extensor muscles (Figure 3). For example, vastus lateralis force reached 1.8–2.5 BW at the instance of deepest squat. However, antagonist forces were also predicted in the biarticular knee flexors to generate the necessary hip and ankle moments (Figure 4). Here, the greatest force in the knee flexors was about 0.6 BW generated by the biceps femoris long head. Interestingly, forces in the vastus muscles were clearly flexion-dependent, while other knee muscles showed subject-specific force patterns. During the flexion phase of squat, force in the gastrocnemius muscles exhibited a rapid decline, directly after knee flexion was initiated, followed by a gradual decrease with increasing knee flexion angle (Figure 4). During the extension phase, this pattern was reversed.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/808027/fbioe-10-808027-HTML-r2/image_m/fbioe-10-808027-g003.jpg)
FIGURE 3. Estimated forces in the knee extensor muscles obtained from the baseline musculoskeletal models (solid lines) and 5th–95th percentile range of the general Monte Carlo simulation outputs (shaded area) for the studied squat trials.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/808027/fbioe-10-808027-HTML-r2/image_m/fbioe-10-808027-g004.jpg)
FIGURE 4. Estimated forces in the knee flexor muscles obtained from the baseline musculoskeletal models (solid lines) and 5th–95th percentile range of the general Monte Carlo simulation outputs (shaded area) for the studied squat trials.
Simultaneous perturbation of the input parameters of all lower limb muscles resulted in considerable variations in the estimated muscle forces (Figures 3, 4 and Supplementary Figures S1, S2). These variations were generally subject-specific (e.g., 3.4 BW variation in vastus lateralis force for K2L compared with 2.2 BW for K8L, both values occurring at deepest squat) and more pronounced at larger knee flexion angles.
Baseline simulations showed relatively large inter-subject variability in KCF estimates, with peak forces ranging from 3 to 5 times body weight. However, the KCFs showed a consistent flexion-dependent pattern across the studied subjects. Simulation error exhibited a similar flexion-dependent trend with the largest errors of up to 100% occurring at deepest squat when compared to the in vivo measured forces (Figure 5). Using simplified muscle models (no muscle force-length characteristics) resulted in a general increase in KCF estimates (up to 50% BW change in the KCF predictions of the baseline OpenSim models) and up to 50% change in the corresponding simulation errors (Supplementary Figure S3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/808027/fbioe-10-808027-HTML-r2/image_m/fbioe-10-808027-g005.jpg)
FIGURE 5. (A) In vivo measured (dotted lines) vs. predicted knee contact forces (solid lines) for the studied squat trials. (B) Simulation errors of the baseline musculoskeletal models. (C) Mean (solid lines) and 5th–95th percentile range (shaded area) of the knee contact force obtained from general Monte Carlo simulation outputs. (D) Mean (solid lines) and 5th–95th percentile range (shaded area) of the modeling error calculated from general Monte Carlo simulations.
The MC simulations revealed that uncertainty in the lower limb muscle parameters can result in considerable variation of the KCF, especially at larger knee flexion angles (e.g., 2.1 BW at 88° knee flexion for K5R). Similarly, the KCF simulation error was greatly influenced by the input uncertainty and showed large variation at deep knee flexion that could explain up to 70% of the simulation error (Figure 5). Here, the knee flexor muscles showed the highest contribution (50–100%) to the overall KCF uncertainty at small flexion angles, whereas knee extensors were the main contributors (70–80%) at larger flexion angles (Figure 6). Variation of the hip muscle parameters had a very small impact until around 50° knee flexion but became more influential thereafter. Ankle muscles had a considerable contribution toward the overall uncertainty when the knee was close to full extension; however, their influence declined with increasing knee flexion angle.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/808027/fbioe-10-808027-HTML-r2/image_m/fbioe-10-808027-g006.jpg)
FIGURE 6. Contribution of different muscle groups toward the overall variability in KCF estimates. Solid lines represent average contributions obtained from the general Monte Carlo analysis with a total number of 72,000 joint reaction force simulations (2000 iterations per subject and muscle group). KH (Knee-Hip) and KA (Knee-Ankle).
Among the studied parameters, MIF and PEN were the most and least influential factors, explaining around 70% and 10% of the overall uncertainty (Figure 7; Table 2). Variability of the TSL had a considerable impact on the KCF estimates until around 50° knee flexion, but its influence declined thereafter. Interestingly, the contribution of the VIA parameter toward the overall KCF uncertainty was almost linearly dependent on the knee flexion angle, with around 60% contribution at 90° knee flexion. Here, perturbation of the VIA of the knee extensor muscles introduced up to 15 mm variation in the muscle moment arms (Supplementary Figure S4) that consequently resulted in around 1.2BW variation in the KCF estimates (also explaining up to 40% of the simulation error, Supplementary Figure S5). Similar to MIF and PEN, variability induced by OIP variation remained consistent over the entire range of knee flexion, explaining about 35% of the overall KCF variability.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/808027/fbioe-10-808027-HTML-r2/image_m/fbioe-10-808027-g007.jpg)
FIGURE 7. Contribution of different muscle parameters toward the overall variability in KCF estimates. Solid lines represent average contribution factors obtained from the general and individual MCs with a total number of 24,000 JRF simulations (4,000 iterations per subject). MIF, maximum isometric force; TSL, tendon slack length; PEN, pennation angle; OIP, origin and insertion point; VIA, VIA points.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/808027/fbioe-10-808027-HTML-r2/image_m/fbioe-10-808027-t002.jpg)
TABLE 2. Contribution factor (in percent) of the input parameters of the studied muscle groups toward the overall variability in the knee contact force estimates at deepest knee flexion angle during squat.
Perturbed models with the lowest RMSE of the simulation results showed the greatest improvements in the KCF predictions (27% decrease in average RMSE, Supplementary Table S2 and Figure S6). However, these improvements were mostly highlighted at larger flexion angles where baseline model errors of up to 100% were decreased to less than 10% for models with optimized muscle parameters reported in Supplementary Table S3–S8. It should be noted, however, that even in these models with minimal KCF RMSE over the complete squat cycle, KCF predictions were not substantially better than those obtained from baseline models at the beginning and end of the cycle. Regarding the discrepancy between the estimated muscle activation and measured EMG signals (Figure 2), errors at the start and end of the squat cycle mainly originated from the inability of the static optimization method to account for the observed muscle co-contractions.
Discussion
This study assessed the influence of muscle model parameter uncertainty on simulation predictions of a generic musculoskeletal model used to estimate KCFs in six instrumented TKA subjects performing a body-weight squat. The results indicate uncertainties in model parameters propagated through the simulation workflow, resulting in a significant impact on the muscle and knee joint contact force predictions. Simultaneous variation of the muscle force-length and geometric pathway parameters of the lower limb muscles resulted in more than 2 BW variations in the KCF and 70% variation in the modeling error. We also found that the contribution of different muscle groups and input parameters toward the overall uncertainty in the KCF outcome may vary over the range of knee flexion angle. These results indicate that, to reduce uncertainty in the predicted KCFs, subject-specific MIF and VIA of the lower limb muscles in generic models must be determined through measurement or calibration, especially when activities with larger flexion angles are considered.
Similar to the in vivo measured tibiofemoral forces, KCFs predicted by the baseline models indicated flexion-dependent patterns, with peaks occurring at the largest flexion angles during squat. While at small flexion angles, KCFs were both over- and under-predicted by the baseline models, at larger knee flexion angles, model estimates were substantially greater than the measured outputs. Other studies that used the CAMS-knee data sets have reported similar KCF errors obtained not only from the same generic model used in the current study but also when using a different generic model (Gait 2,392) (Schellenberg et al., 2018; Imani Nejad et al., 2020). In an another modeling study (Ding et al., 2016), when KCF estimates from squat simulation were compared against the in vivo data reported in the sixth Grand Challenge Competition (Fregly et al., 2012), errors of up to 100% BW were identified. van Rossom et al. (2018) used a six DoF knee model integrated in a full body musculoskeletal model to estimate KCFs in 15 young healthy subjects performing squats. They found an average maximum KCF of around 4 times BW, which is comparable with the results reported by Shelburne and Pandy (2002), who used a simpler musculoskeletal model. Despite using a different generic model, the results of our baseline model simulations concur well with those reported from these aforementioned studies. It is important to note, however, that the majority of predictions presented to date are generally higher than the tibiofemoral forces measured in living subjects performing squat (2.2-3BW, (Mundermann et al., 2008; Fregly et al., 2012; Trepczynski et al., 2012; Bergmann et al., 2014; Mizu-uchi et al., 2015; Catelli et al., 2019)).
We found varying contributions between the studied muscle groups toward the overall uncertainty in KCF estimates. Here, the knee flexor muscles played a greater role at lower flexion angles, while knee extensors became considerably more influential at larger flexion angles (Figure 6). We also found a slightly higher contribution of the knee extensor muscles during the extension phase of the squat compared with the flexion phase. These findings can be partially explained by the different activation patterns of the leg muscles during squat. Escamilla et al. (2001) measured EMG patterns of different muscle groups in 10 subjects and found that the peak quadriceps activity occurs at 80–90° knee flexion and that the activity is 25–50% greater in the knee extension phase. Conversely, knee flexor muscles are more active at smaller flexion angles during squat. In particular, at the beginning of the flexion phase, the biceps femoris, gastrocnemius, and semitendinosus all act together to initiate flexion of the knee (Robertson et al., 2008). Although knee flexor muscles remain less active after their initial burst (Wilk et al., 1996; Isear et al., 1997), biceps femoris long head, semimembranosus, and semitendinosus continue lengthening and thus apply a considerable passive force until around 60° flexion (Sinclair et al., 2017). From 60° to 90°, the hamstring muscles remain nearly isometric and minimally activated (Wilk et al., 1996; Isear et al., 1997; Sinclair et al., 2017). This can explain our results indicating a decline in contribution of knee flexors toward uncertainty in the KCF after 60° knee flexion. The increasing contribution of the hip muscles after around 50° knee flexion is likely due to the increases in hip moment magnitudes (in all 3 DoF) and increased inaccuracies in muscle pathway representations and resulting moment arms of these muscles.
MIF was found to be the most influential parameter on the resultant KCF throughout the entire range of knee flexion (Figure 7), which concurs with findings of previous probabilistic studies performed on level walking (Navacchia et al., 2016; Zuk et al., 2018). However, it should be mentioned that, due to the nature of SO, changing the MIF of all muscles surrounding the joint by the same amount (e.g., 30%) is unlikely to greatly alter the KCF estimates. Regarding the substantial inter-subject variability of MIF of the lower limb muscles (Friederich and Brand, 1990; Ward et al., 2009; Yucesoy et al., 2010; Bakenecker et al., 2020; Montefiori et al., 2020), using the baseline values from generic models may result in sizable KCF errors. Medical imaging techniques (e.g., MR and ultrasound imaging) and isokinetic dynamometers can be used to estimate subject-specific MIF of major muscle groups (Maganaris et al., 2001; Correa and Pandy, 2011; Ivanovic and Dopsaj, 2013; Kainz et al., 2018; Higa et al., 2019). Although, these methods have limited accuracies and are not applicable to all individual muscles, subject-specific determination of MIF inputs for major muscle groups seems to be able to considerably reduce uncertainty in KCF estimates of generic musculoskeletal models.
After MIF, VIA showed the largest impact on KCF predictions (Figure 7). Variation of the muscle VIA influences the muscle moment arm and thereby the recruitment strategy selected by the static optimization approach to balance the external joint moment. In the current study, variation in VIA was found to have a flexion-dependent contribution toward the overall KCF variability, with a maximum contribution of around 70% at deepest squat. Using a global probabilistic approach, Navacchia et al. (2016) found that the uncertainty in VIA can explain around 50% of the overall variability in the estimated KCFs for level walking. However, they did not provide details on how VIA contribution changes over the range of knee flexion. The large contribution of VIA uncertainty in our study is due to the high range of loaded knee flexion during squat, whereas the flexion-dependent pattern is likely to be associated with the decreased accuracy in muscle pathway representation in generic models. Musculoskeletal models have been mainly developed to assess walking tasks, and hence their muscle moment arms are validated for a small range of knee and hip flexion (Delp et al., 1990; Arnold et al., 2010; Rajagopal et al., 2016). Therefore, as the knee flexion angle increases, variation of the VIA points of the highly loaded muscles (mainly knee extensors) has a more highlighted effect on the KCF output. This is not limited to the knee muscles. For example, variation of the VIA of the hip muscles (mainly hip adductors) can change the entire force distribution between the hip muscles, which consequently affects the KCF due to the changes induced in biarticular muscle forces. Regarding the increasing contribution of VIA and extensor muscles with increasing knee flexion angle, the observed flexion-dependent errors in our baseline simulations most likely originate from inaccurate representation of the knee extensor muscle pathways in the scaled generic models. Our results, therefore, highlight the importance of subject-specific modeling of the muscle pathway geometry when activities with large knee flexions are investigated using generic musculoskeletal models. To reduce uncertainty in muscle moment arms, more accurate VIA pathways can be obtained from MRI images (Blemker et al., 2007; Scheys et al., 2008); however, implementation of these parameters in the model necessitates the use of simplified wrapping objects. As a result, even personalized models have an inherent uncertainty in their muscle geometry representation. Moreover, other sources of inter-subject variability such as bone geometry (e.g., femoral anteversion angle or tibial torsion) may also influence the muscle attachment points and pathway geometry, thereby impacting the predicted KCFs (Kainz et al., 2020; Modenese et al., 2021).
Our data suggest that parameters with low impacts can be directly taken from the generic model (or scaled based on subject’s anthropometry). The pennation angle of the muscle fibers (PEN) is a clear example of such a parameter, considering its relatively small impact on the KCF estimates (Figure 7, also reported by Navacchia et al. (2016)). For the studied squat cycles, TSL had a considerable impact on KCF variability, especially between 20° and 50° knee flexion. The large sensitivity of predicted KCFs to uncertainty in TSL has been also reported in previous studies on level walking (Myers et al., 2015; Carbone et al., 2016; Navacchia et al., 2016). However, plausibly due to the very large ±20% range of variation of the TSL (compared to 2–9% used in the current study), Carbone and co-workers found TSL to be more influential than MIF in prediction of the muscle forces. Regarding technical challenges in subject-specific measurement of TSL in living subjects, the baseline values from the generic models are either directly used or adjusted such that the differences between experimental and model-based joint moments is minimized (Van Campen et al., 2014; Heinen et al., 2019). It is also possible to ignore the force-length characteristics of the muscle-tendon units in the static optimization procedure. This approach excludes the role of the tendon slack length and optimal fiber length, while the maximum isometric force and muscle pathway geometry (which determines the muscle moment arm) remain influential. In the current study ignoring muscle force–length relationships did not play a dominant role in the knee contact force estimates; however, it did result in a general increase in KCF estimates and changes of up to 50% in simulation error (Supplementary Figure S3). Similar results were reported when muscle models were simplified during KCF estimation in level walking simulations (Modenese et al., 2018). Therefore, inclusion of force-length characteristics of the muscles in static optimization seems to be necessary in order to reduce uncertainty in modeling outcomes. Hence, a probabilistic modeling approach, as used in the current study, can provide an estimate of the range of possible errors that might originate from TSL uncertainty. Nevertheless, a complete understanding of the real influence of PEN and TSL on the model outcomes is limited by the 1D representation of the muscles in the selected musculoskeletal model. Therefore, volumetric modeling of muscles (Murai et al., 2016; Modenese and Kohout, 2020) may better clarify interaction between uncertainty in muscle parameters and variation in the joint contact force estimates, but such modeling approaches are extremely expensive and technically challenging.
Our results indicate that uncertainty in muscle parameters cannot completely explain KCF errors when generic models are used to simulate activities involving high knee joint flexion (Supplementary Table S2 and Figure S6). As a result, considerable errors, especially at small flexion angles, can still be expected even with subject-specific muscle parameters. A comparison between the predicted muscle activation and measured EMG signals (Figure 2) confirms the inability of utilized static optimization technique to account for co-contraction of the flexor and extensor muscles at the beginning and end of the squat cycle. Here, it is entirely plausible that subjects change neuromuscular control from stability and balance during standing to a more economical/preservation mechanism to avoid muscle overloading and/or damage during higher loading scenarios. Such control programs might nicely be accounted for through adapting the cost-function within the SO to account for, for example, muscle co-contractions. Thresholds or optimal approaches to implement and validate such novel methods would, however, need to be fully verified in future investigations.
Several limitations need to be considered when interpreting the results of this study. First, we chose the generic model developed by Rajagopal et al. (2016) because it has been widely used to assess KCFs during different activities. It is likely that uncertainty in input parameters may propagate differently into the outcomes of other generic models with different numbers of muscles and joint DoFs (Cleather and Bull, 2011; Kainz et al., 2016). Moreover, while CAMS-Knee data sets provide access to the fluoroscopically measured tibiofemoral implant kinematics, lack of subject-specific patellofemoral kinematics and bone geometries precluded prescription of the tibiofemoral kinematics within the IK process. As a result, similar to many other studies in the field, the model joints were driven by skin marker trajectories, and their kinematics may, therefore, differ from the subject-specific joint movements. In addition, while a large number of factors may affect the accuracy of musculoskeletal modeling knee contact force predictions, we only perturbed the muscle parameters in the selected generic musculoskeletal model. Although inclusion of other sources of uncertainty (e.g., skin marker location and body segment parameters) may also change the predicted results, previous probabilistic simulations reported no major influence of such factors (Valente et al., 2014; Myers et al., 2015; Navacchia et al., 2016). It is also important to mention that our probabilistic approach did not account for inter-subject variability of bone geometry. Such variations can impact the 3D pathways of the muscles and thereby the joint contact force estimation (Kainz et al., 2020; Modenese et al., 2021). More importantly, we only used SO to solve the muscle redundancy problem through minimizing the sum of the squared muscle activation levels. It is known that different optimization approaches may lead to different modeling outcomes (Li et al., 1998; Mokhtarzadeh et al., 2014). In addition, the unknown level of antagonistic muscle co-contraction (Trepczynski et al., 2018) and inability of the routine muscle optimization approaches to accurately predict such contributions (Mokhtarzadeh et al., 2014) may also hinder a clear understanding of uncertainty propagation in musculoskeletal modeling. Moreover, the assumed distributions for the uncertain muscle parameters also influence the variability in the predicted KCF, and since the range of variation in input parameters for some muscles was missing in the literature, the average of the reported values for other muscles was used to generate random inputs for the MC simulations. Finally, we assessed KCFs during squat cycles with a limited range of flexion performed by only six TKA subjects. Uncertainty propagation through the modeling framework may be different at larger knee flexion angles during squats performed by young healthy population.
Conclusion
Personalization of musculoskeletal models remains technically and financially challenging. This study provides a more complete understanding of the impact of uncertainty in the input parameters on the estimated KCFs and can help to identify the most influential parameters that need to be accurately measured and implemented. The results suggest that attention should be paid to the subject-specific measurement of MIF and VIA in order to ensure errors in musculoskeletal modeling remain minimal. In addition, adaptive muscle optimization techniques that enable accurate prediction of the muscle co-contractions at the beginning and end of the squat cycle might well present improvements in KCF prediction deficits revealed in this study.
Data Availability Statement
Requests to access the data sets used in this study should be directed to: https://cams-knee.orthoload.com/data/data-download.
Ethics Statement
The studies involving human participants were reviewed and approved by Charité Ethikkommission (EA4/069/06) and ETH Ethikkommission (EK 2013-N-90). The patients/participants provided their written informed consent to participate in this study.
Author Contributions
SH, CS, AM, and WT contributed to the conception and design of the study. SH and AV developed the scripts, performed simulations, and conducted the statistical analysis. PS, PD, AT, and JD collected and processed the experimental data and contributed to the design of the results section. All authors contributed to preparation of the manuscript and approved the submission.
Funding
This work was financially supported by B. Braun Aesculap, Tuttlingen. Partial funding was provided by the German Research Foundation (DFG) (TR 1657/1-1, EH 422-2-1/MO 3865-1-1, DA 1786/5-1, BE 804/18-1, and SFB 1444), Federal Ministry of Education and Research (OVERLOAD-PrevOP (01EC1408A) and workHEALTH (01 EC 1905D)), European Regional Development Fund (ERDF/EFRE 16409608, OrthoLoadLab), and the OrthoLoadClub. Open access funding provided by ETH Zurich.
Conflict of Interest
AM was employed by Aesculap AG.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We gratefully acknowledge the financial support of B. Braun Aesculap, Tuttlingen. We would also like to thank RMS Foundation for the collection of the original CAMS-Knee data sets.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe.2022.808027/full#supplementary-material
References
Arnold, E. M., Ward, S. R., Lieber, R. L., and Delp, S. L. (2010). A Model of the Lower Limb for Analysis of Human Movement. Ann. Biomed. Eng. 38, 269–279. doi:10.1007/s10439-009-9852-5
Ataei, G., Abedi, R., Mohammadi, Y., and Fatouraee, N. (2020). Analysing the Effect of Wearable Lift-Assist Vest in Squat Lifting Task Using Back Muscle EMG Data and Musculoskeletal Model. Phys. Eng. Sci. Med. 43, 651–658. doi:10.1007/s13246-020-00872-5
Bakenecker, P., Raiteri, B. J., and Hahn, D. (2020). Force Enhancement in the Human Vastus Lateralis Is Muscle-length-dependent Following Stretch but Not during Stretch. Eur. J. Appl. Physiol. 120, 2597–2610. doi:10.1007/s00421-020-04488-1
Barry, M. J., Kwon, T. H., and Dhaher, Y. Y. (2010). Probabilistic Musculoskeletal Modeling of the Knee: A Preliminary Examination of an ACL-Reconstruction. Annu. Int. Conf. IEEE Eng. Med. Biol. Soc. 2010, 5440–5443. doi:10.1109/IEMBS.2010.5626511
Bergmann, G., Bender, A., Graichen, F., Dymke, J., Rohlmann, A., Trepczynski, A., et al. (2014). Standardized Loads Acting in Knee Implants. Plos One 9 (1), e86035. doi:10.1371/journal.pone.0086035
Blemker, S. S., Asakawa, D. S., Gold, G. E., and Delp, S. L. (2007). Image-based Musculoskeletal Modeling: Applications, Advances, and Future Opportunities. J. Magn. Reson. Imaging 25, 441–451. doi:10.1002/jmri.20805
Bolliger, M., Banz, R., Dietz, V., and Lünenburger, L. (2008). Standardized Voluntary Force Measurement in a Lower Extremity Rehabilitation Robot. J. NeuroEngineering Rehabil. 5, 23. doi:10.1186/1743-0003-5-23
Burden, A. (2010). How Should We Normalize Electromyograms Obtained from Healthy Participants? what We Have Learned from over 25years of Research. J. Electromyogr. Kinesiol. 20, 1023–1035. doi:10.1016/j.jelekin.2010.07.004
Carbone, V., Van der Krogt, M. M., Koopman, H. F. J. M., and Verdonschot, N. (2016). Sensitivity of Subject-specific Models to Hill Muscle-Tendon Model Parameters in Simulations of Gait. J. Biomechanics 49, 1953–1960. doi:10.1016/j.jbiomech.2016.04.008
Catelli, D. S., Wesseling, M., Jonkers, I., and Lamontagne, M. (2019). A Musculoskeletal Model Customized for Squatting Task. Comput. Methods Biomechanics Biomed. Eng. 22, 21–24. doi:10.1080/10255842.2018.1523396
Charles, J. P., Grant, B., D’Août, K., and Bates, K. T. (2020). Subject-specific Muscle Properties from Diffusion Tensor Imaging Significantly Improve the Accuracy of Musculoskeletal Models. J. Anat. 237, 941–959. doi:10.1111/joa.13261
Chen, Z., Zhang, Z., Wang, L., Li, D., Zhang, Y., and Jin, Z. (2016). Evaluation of a Subject-specific Musculoskeletal Modelling Framework for Load Prediction in Total Knee Arthroplasty. Med. Eng. Phys. 38, 708–716. doi:10.1016/j.medengphy.2016.04.010
Cleather, D. J., and Bull, A. M. J. (2011). Knee and Hip Joint Forces - Sensitivity to the Degrees of Freedom Classification at the Knee. Proc. Inst. Mech. Eng. H. 225, 621–626. doi:10.1177/0954411911399975
Correa, T. A., and Pandy, M. G. (2011). A Mass-Length Scaling Law for Modeling Muscle Strength in the Lower Limb. J. Biomechanics 44, 2782–2789. doi:10.1016/j.jbiomech.2011.08.024
De Groote, F., Van Campen, A., Jonkers, I., and De Schutter, J. (2010). Sensitivity of Dynamic Simulations of Gait and Dynamometer Experiments to Hill Muscle Model Parameters of Knee Flexors and Extensors. J. Biomech. 43, 1876–1883. doi:10.1016/j.jbiomech.2010.03.022
Delp, S. L., Loan, J. P., Hoy, M. G., Zajac, F. E., Topp, E. L., and Rosen, J. M. (1990). An Interactive Graphics-Based Model of the Lower Extremity to Study Orthopaedic Surgical Procedures. IEEE Trans. Biomed. Eng. 37, 757–767. doi:10.1109/10.102791
Ding, Z., Nolte, D., Kit Tsang, C., Cleather, D. J., Kedgley, A. E., and Bull, A. M. (2016). In Vivo Knee Contact Force Prediction Using Patient-specific Musculoskeletal Geometry in a Segment-Based Computational Model. J. Biomech. Eng. 138, 021018. doi:10.1115/1.4032412
Escamilla, R. F., Fleisig, G. S., Zheng, N., Lander, J. E., Barrentine, S. W., Andrews, J. R., et al. (2001). Effects of Technique Variations on Knee Biomechanics during the Squat and Leg Press. Med. Sci. Sports Exerc. 33, 1552–1566. doi:10.1097/00005768-200109000-00020
Fernandez, J., Zhang, J., Heidlauf, T., Sartori, M., Besier, T., Röhrle, O., et al. (2016). Multiscale Musculoskeletal Modelling, Data-Model Fusion and Electromyography-Informed Modelling. Interface Focus. 6, 20150084. doi:10.1098/rsfs.2015.0084
Fregly, B. J., Besier, T. F., Lloyd, D. G., Delp, S. L., Banks, S. A., Pandy, M. G., et al. (2012). Grand Challenge Competition to Predict In Vivo Knee Loads. J. Orthop. Res. 30, 503–513. doi:10.1002/jor.22023
Friederich, J. A., and Brand, R. A. (1990). Muscle Fiber Architecture in the Human Lower Limb. J. Biomechanics 23, 91–95. doi:10.1016/0021-9290(90)90373-b
Fukunaga, T., Roy, R. R., Shellock, F. G., Hodgson, J. A., and Edgerton, V. R. (1996). Specific Tension of Human Plantar Flexors and Dorsiflexors. J. Appl. Physiology 80, 158–165. doi:10.1152/jappl.1996.80.1.158
Gerus, P., Sartori, M., Besier, T. F., Fregly, B. J., Delp, S. L., Banks, S. A., et al. (2013). Subject-specific Knee Joint Geometry Improves Predictions of Medial Tibiofemoral Contact Forces. J. Biomechanics 46, 2778–2786. doi:10.1016/j.jbiomech.2013.09.005
Handsfield, G. G., Meyer, C. H., Hart, J. M., Abel, M. F., and Blemker, S. S. (2014). Relationships of 35 Lower Limb Muscles to Height and Body Mass Quantified Using MRI. J. Biomechanics 47, 631–638. doi:10.1016/j.jbiomech.2013.12.002
Heinen, F., S Rensen, S. R., King, M., Lewis, M., Lund, M. E., Rasmussen, J., et al. (2019). Muscle-tendon Unit Parameter Estimation of a Hill-type Musculoskeletal Model Based on Experimentally Obtained Subject-specific Torque Profiles. J. Biomech. Eng. 141 (6), 061005. doi:10.1115/1.4043356
Heinen, F., Lund, M. E., Rasmussen, J., and De Zee, M. (2016). Muscle-tendon Unit Scaling Methods of Hill-type Musculoskeletal Models: An Overview. Proc. Inst. Mech. Eng. H. 230, 976–984. doi:10.1177/0954411916659894
Heron, N. (2015). A 'novel' Model for Integrating Sport and Exercise Medicine (SEM) and Musculoskeletal (MSK) Management into Primary Care in the UK. BMJ Open Sport Exerc Med. 1, e000027. doi:10.1136/bmjsem-2015-000027
Higa, M., Araki, N., Nakatani, S., and Toji, H. (2019). Measurements of Isometric Strength and Electromyography of Elbow Flexors under Active and Passive Conditions. Abe 8, 124–129. doi:10.14326/abe.8.124
Hoang, H. X., Diamond, L. E., Lloyd, D. G., and Pizzolato, C. (2019). A Calibrated EMG-Informed Neuromusculoskeletal Model Can Appropriately Account for Muscle Co-contraction in the Estimation of Hip Joint Contact Forces in People with Hip Osteoarthritis. J. Biomechanics 83, 134–142. doi:10.1016/j.jbiomech.2018.11.042
Imani nejad, Z., Khalili, K., Hosseini Nasab, S. H., Schütz, P., Damm, P., Trepczynski, A., et al. (2020). The Capacity of Generic Musculoskeletal Simulations to Predict Knee Joint Loading Using the CAMS-Knee Datasets. Ann. Biomed. Eng. 48, 1430–1440. doi:10.1007/s10439-020-02465-5
Isear, J. A., Erickson, J. C., and Worrell, T. W. (1997). EMG Analysis of Lower Extremity Muscle Recruitment Patterns during an Unloaded Squat. Med. & Sci. Sports & Exerc. 29, 532–539. doi:10.1097/00005768-199704000-00016
Ivanovic, J., and Dopsaj, M. (2013). Reliability of Force-Time Curve Characteristics during Maximal Isometric Leg Press in Differently Trained High-Level Athletes. Measurement 46, 2146–2154. doi:10.1016/j.measurement.2013.03.008
Jamwal, P. K., Hussain, S., Tsoi, Y. H., and Xie, S. Q. (2020). Musculoskeletal Model for Path Generation and Modification of an Ankle Rehabilitation Robot. IEEE Trans. Human-Mach. Syst. 50, 373–383. doi:10.1109/thms.2020.2989688
Kainz, H., Goudriaan, M., Falisse, A., Huenaerts, C., Desloovere, K., De Groote, F., et al. (2018). The Influence of Maximum Isometric Muscle Force Scaling on Estimated Muscle Forces from Musculoskeletal Models of Children with Cerebral Palsy. Gait Posture 65, 213–220. doi:10.1016/j.gaitpost.2018.07.172
Kainz, H., Hoang, H. X., Stockton, C., Boyd, R. R., Lloyd, D. G., and Carty, C. P. (2017). Accuracy and Reliability of Marker-Based Approaches to Scale the Pelvis, Thigh, and Shank Segments in Musculoskeletal Models. J. Appl. Biomechanics 33, 354–360. doi:10.1123/jab.2016-0282
Kainz, H., Killen, B. A., Wesseling, M., Perez-Boerema, F., Pitto, L., Garcia Aznar, J. M., et al. (2020). A Multi-Scale Modelling Framework Combining Musculoskeletal Rigid-Body Simulations with Adaptive Finite Element Analyses, to Evaluate the Impact of Femoral Geometry on Hip Joint Contact Forces and Femoral Bone Growth. PLoS One 15, e0235966. doi:10.1371/journal.pone.0235966
Kainz, H., Modenese, L., Lloyd, D. G., Maine, S., Walsh, H. P. J., and Carty, C. P. (2016). Joint Kinematic Calculation Based on Clinical Direct Kinematic versus Inverse Kinematic Gait Models. J. Biomechanics 49, 1658–1669. doi:10.1016/j.jbiomech.2016.03.052
Kainz, H., Wesseling, M., and Jonkers, I. (2021). Generic Scaled versus Subject-specific Models for the Calculation of Musculoskeletal Loading in Cerebral Palsy Gait: Effect of Personalized Musculoskeletal Geometry Outweighs the Effect of Personalized Neural Control. Clin. Biomech. 87, 105402. doi:10.1016/j.clinbiomech.2021.105402
Koller, W., Baca, A., and Kainz, H. (2021). Impact of Scaling Errors of the Thigh and Shank Segments on Musculoskeletal Simulation Results. Gait Posture 87, 65–74. doi:10.1016/j.gaitpost.2021.02.016
Lamberto, G., Martelli, S., Cappozzo, A., and Mazzà, C. (2017). To what Extent Is Joint and Muscle Mechanics Predicted by Musculoskeletal Models Sensitive to Soft Tissue Artefacts? J. Biomechanics 62, 68–76. doi:10.1016/j.jbiomech.2016.07.042
Langholz, J. B., Westman, G., and Karlsteen, M. (2016). Musculoskeletal Modelling in Sports - Evaluation of Different Software Tools with Focus on Swimming. Procedia Eng. 147 (147), 281–287. doi:10.1016/j.proeng.2016.06.278
Li, G., Kawamura, K., Barrance, P., Chao, E. Y. S., and Kaufman, K. (1998). Prediction of Muscle Recruitment and its Effect on Joint Reaction Forces during Knee Exercises. Ann. Biomed. Eng. 26, 725–733. doi:10.1114/1.104
Li, N., Yang, T., Yang, Y., Yu, P., Xue, X., Zhao, X., et al. (2020). Bioinspired Musculoskeletal Model-Based Soft Wrist Exoskeleton for Stroke Rehabilitation. J. Bionic Eng. 17, 1163–1174. doi:10.1007/s42235-020-0101-9
Lloyd, D. G., and Besier, T. F. (2003). An EMG-Driven Musculoskeletal Model to Estimate Muscle Forces and Knee Joint Moments In Vivo. J. Biomechanics 36, 765–776. doi:10.1016/s0021-9290(03)00010-1
Lundberg, H. J., Foucher, K. C., Andriacchi, T. P., and Wimmer, M. A. (2012). Direct Comparison of Measured and Calculated Total Knee Replacement Force Envelopes during Walking in the Presence of Normal and Abnormal Gait Patterns. J. Biomechanics 45, 990–996. doi:10.1016/j.jbiomech.2012.01.015
Ma, Y., Jiang, S., Mithraratne, K., Wilson, N., Yu, Y., and Zhang, Y. (2021). The Effect of Musculoskeletal Model Scaling Methods on Ankle Joint Kinematics and Muscle Force Prediction during Gait for Children with Cerebral Palsy and Equinus Gait. Comput. Biol. Med. 134, 104436. doi:10.1016/j.compbiomed.2021.104436
Maganaris, C. N., Baltzopoulos, V., Ball, D., and Sargeant, A. J. (2001). In Vivo specific Tension of Human Skeletal Muscle. J. Appl. Physiology 90, 865–872. doi:10.1152/jappl.2001.90.3.865
Millard, M., Uchida, T., Seth, A., and Delp, S. L. (2013). Flexing Computational Muscle: Modeling and Simulation of Musculotendon Dynamics. J. Biomech. Eng. 135, 021005. doi:10.1115/1.4023390
Mizu-uchi, H., Colwell, C. W., Flores-Hernandez, C., Fregly, B. J., Matsuda, S., and D’Lima, D. D. (2015). Patient-Specific Computer Model of Dynamic Squatting after Total Knee Arthroplasty. J. Arthroplasty 30, 870–874. doi:10.1016/j.arth.2014.12.021
Modenese, L., Barzan, M., and Carty, C. P. (2021). Dependency of Lower Limb Joint Reaction Forces on Femoral Version. Gait Posture 88, 318–321. doi:10.1016/j.gaitpost.2021.06.014
Modenese, L., and Kohout, J. (2020). Automated Generation of Three-Dimensional Complex Muscle Geometries for Use in Personalised Musculoskeletal Models. Ann. Biomed. Eng. 48, 1793–1804. doi:10.1007/s10439-020-02490-4
Modenese, L., Montefiori, E., Wang, A., Wesarg, S., Viceconti, M., and Mazzà, C. (2018). Investigation of the Dependence of Joint Contact Forces on Musculotendon Parameters Using a Codified Workflow for Image-Based Modelling. J. Biomechanics 73, 108–118. doi:10.1016/j.jbiomech.2018.03.039
Mokhtarzadeh, H., Perraton, L., Fok, L., Muñoz, M. A., Clark, R., Pivonka, P., et al. (2014). A Comparison of Optimisation Methods and Knee Joint Degrees of Freedom on Muscle Force Predictions during Single-Leg Hop Landings. J. Biomechanics 47, 2863–2868. doi:10.1016/j.jbiomech.2014.07.027
Montefiori, E., Kalkman, B. M., Henson, W. H., Paggiosi, M. A., Mccloskey, E. V., and Mazzà, C. (2020). MRI-based Anatomical Characterisation of Lower-Limb Muscles in Older Women. PLoS One 15, e0242973. doi:10.1371/journal.pone.0242973
Mündermann, A., Dyrby, C. O., D'Lima, D. D., Colwell, C. W., and Andriacchi, T. P. (2008). In Vivo knee Loading Characteristics during Activities of Daily Living as Measured by an Instrumented Total Knee Replacement. J. Orthop. Res. 26, 1167–1172. doi:10.1002/jor.20655
Murai, A., Endo, Y., and Tada, M. (2016). Anatomographic Volumetric Skin-Musculoskeletal Model and its Kinematic Deformation with Surface-Based SSD. IEEE Robot. Autom. Lett. 1 (2), 1103–1109. doi:10.1109/lra.2016.2524069
Myers, C. A., Laz, P. J., Shelburne, K. B., and Davidson, B. S. (2015). A Probabilistic Approach to Quantify the Impact of Uncertainty Propagation in Musculoskeletal Simulations. Ann. Biomed. Eng. 43, 1098–1111. doi:10.1007/s10439-014-1181-7
Navacchia, A., Myers, C. A., Rullkoetter, P. J., and Shelburne, K. B. (2016). Prediction of In Vivo Knee Joint Loads Using a Global Probabilistic Analysis. J. Biomech. Eng. 138, 4032379. doi:10.1115/1.4032379
Pizzolato, C., Lloyd, D. G., Sartori, M., Ceseracciu, E., Besier, T. F., Fregly, B. J., et al. (2015). CEINMS: A Toolbox to Investigate the Influence of Different Neural Control Solutions on the Prediction of Muscle Excitation and Joint Moments during Dynamic Motor Tasks. J. Biomechanics 48, 3929–3936. doi:10.1016/j.jbiomech.2015.09.021
Rajagopal, A., Dembia, C. L., Demers, M. S., Delp, D. D., Hicks, J. L., and Delp, S. L. (2016). Full-Body Musculoskeletal Model for Muscle-Driven Simulation of Human Gait. IEEE Trans. Biomed. Eng. 63, 2068–2079. doi:10.1109/tbme.2016.2586891
Richards, C., and Higginson, J. S. (2010). Knee Contact Force in Subjects with Symmetrical OA Grades: Differences between OA Severities. J. Biomechanics 43, 2595–2600. doi:10.1016/j.jbiomech.2010.05.006
Robertson, D. G. E., Wilson, J.-M. J., and Pierre, T. A. S. (2008). Lower Extremity Muscle Functions during Full Squats. J. Appl. Biomechanics 24, 333–339. doi:10.1123/jab.24.4.333
Schellenberg, F., Taylor, W. R., Trepczynski, A., List, R., Kutzner, I., Schütz, P., et al. (2018). Evaluation of the Accuracy of Musculoskeletal Simulation during Squats by Means of Instrumented Knee Prostheses. Med. Eng. Phys. 61, 95–99. doi:10.1016/j.medengphy.2018.09.004
Scheys, L., Spaepen, A., Suetens, P., and Jonkers, I. (2008). Calculated Moment-Arm and Muscle-Tendon Lengths during Gait Differ Substantially Using MR Based versus Rescaled Generic Lower-Limb Musculoskeletal Models. Gait Posture 28, 640–648. doi:10.1016/j.gaitpost.2008.04.010
Schmid, J., Sandholm, A., Chung, F., Thalmann, D., Delingette, H., and Magnenat-Thalmann, N. (2009). “Musculoskeletal Simulation Model Generation from MRI Data Sets and Motion Capture Data,” in Recent Advances in the 3d Physiological Human. Editors N. Magnenat-Thalmann, J. Zhang, and D. Feng (London: Springer). doi:10.1007/978-1-84882-565-9_1
Seow, D., Graham, I., and Massey, A. (2020). Prediction Models for Musculoskeletal Injuries in Professional Sporting Activities: A Systematic Review. Transl. Sports Med. 3, 505–517. doi:10.1002/tsm2.181
Shelburne, K. B., and Pandy, M. G. (2002). A Dynamic Model of the Knee and Lower Limb for Simulating Rising Movements. Comput. Methods Biomechanics Biomed. Eng. 5, 149–159. doi:10.1080/10255840290010265
Sinclair, J., Brooks, D., and Atkins, S. (2017). An Examination of the Hamstring and the Quadriceps Muscle Kinematics during the Front and Back Squat in Males. Bjhpa 9, 37–45. doi:10.29359/bjhpa.09.1.04
Stylianou, A. P., Guess, T. M., and Kia, M. (2013). Multibody Muscle Driven Model of an Instrumented Prosthetic Knee during Squat and Toe Rise Motions. J. Biomech. Eng. 135, 041008. doi:10.1115/1.4023982
Taylor, W. R., Schütz, P., Bergmann, G., List, R., Postolka, B., Hitz, M., et al. (2017). A Comprehensive Assessment of the Musculoskeletal System: The CAMS-Knee Data Set. J. Biomechanics 65, 32–39. doi:10.1016/j.jbiomech.2017.09.022
Trepczynski, A., Kutzner, I., Kornaropoulos, E., Taylor, W. R., Duda, G. N., Bergmann, G., et al. (2012). Patellofemoral Joint Contact Forces during Activities with High Knee Flexion. J. Orthop. Res. 30, 408–415. doi:10.1002/jor.21540
Trepczynski, A., Kutzner, I., Schwachmeyer, V., Heller, M. O., Pfitzner, T., and Duda, G. N. (2018). Impact of Antagonistic Muscle Co-contraction on In Vivo Knee Contact Forces. J. NeuroEngineering Rehabil. 15, 101. doi:10.1186/s12984-018-0434-3
Valente, G., Pitto, L., Testi, D., Seth, A., Delp, S. L., Stagni, R., et al. (2014). Are Subject-specific Musculoskeletal Models Robust to the Uncertainties in Parameter Identification? PLoS One 9, e112625. doi:10.1371/journal.pone.0112625
Van campen, A., Pipeleers, G., De Groote, F., Jonkers, I., and De Schutter, J. (2014). A New Method for Estimating Subject-specific Muscle-Tendon Parameters of the Knee Joint Actuators: a Simulation Study. Int. J. Numer. Meth. Biomed. Engng. 30, 969–987. doi:10.1002/cnm.2639
Van rossom, S., Smith, C. R., Thelen, D. G., Vanwanseele, B., Van Assche, D., and Jonkers, I. (2018). Knee Joint Loading in Healthy Adults during Functional Exercises: Implications for Rehabilitation Guidelines. J. Orthop. Sports Phys. Ther. 48, 162–173. doi:10.2519/jospt.2018.7459
Ward, S. R., Eng, C. M., Smallwood, L. H., and Lieber, R. L. (2009). Are Current Measurements of Lower Extremity Muscle Architecture Accurate? Clin. Orthop. Relat. Res. 467, 1074–1082. doi:10.1007/s11999-008-0594-8
Wilk, K. E., Escamilla, R. F., Fleisig, G. S., Barrentine, S. W., Andrews, J. R., and Boyd, M. L. (1996). A Comparison of Tibiofemoral Joint Forces and Electromyographic Activit during Open and Closed Kinetic Chain Exercises. Am. J. Sports Med. 24, 518–527. doi:10.1177/036354659602400418
Winby, C. R., Lloyd, D. G., and Kirk, T. B. (2008). Evaluation of Different Analytical Methods for Subject-specific Scaling of Musculotendon Parameters. J. Biomechanics 41, 1682–1688. doi:10.1016/j.jbiomech.2008.03.008
Worsley, P., Stokes, M., and Taylor, M. (2011). Predicted Knee Kinematics and Kinetics during Functional Activities Using Motion Capture and Musculoskeletal Modelling in Healthy Older People. Gait Posture 33, 268–273. doi:10.1016/j.gaitpost.2010.11.018
Yucesoy, C. A., Ateş, F., Akgün, U., and Karahan, M. (2010). Measurement of Human Gracilis Muscle Isometric Forces as a Function of Knee Angle, Intraoperatively. J. Biomechanics 43, 2665–2671. doi:10.1016/j.jbiomech.2010.06.002
Keywords: probabilistic analysis, musculoskeletal modeling, uncertainty, muscle parameters, knee contact force
Citation: Hosseini Nasab SH, Smith CR, Maas A, Vollenweider A, Dymke J, Schütz P, Damm P, Trepczynski A and Taylor WR (2022) Uncertainty in Muscle–Tendon Parameters can Greatly Influence the Accuracy of Knee Contact Force Estimates of Musculoskeletal Models. Front. Bioeng. Biotechnol. 10:808027. doi: 10.3389/fbioe.2022.808027
Received: 02 November 2021; Accepted: 03 May 2022;
Published: 03 June 2022.
Edited by:
Marwan El-Rich, Khalifa University, United Arab EmiratesReviewed by:
Hans Kainz, University of Vienna, AustriaNathalie Alexander, Children’s Hospital of Eastern Switzerland, Switzerland
Copyright © 2022 Hosseini Nasab, Smith, Maas, Vollenweider, Dymke, Schütz, Damm, Trepczynski and Taylor. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Seyyed Hamed Hosseini Nasab, seyyed.hosseini@hest.ethz.ch; William R. Taylor, bt@ethz.ch