Three-dimensional dynamic simulation of total knee replacement motion during a step-up task.

A three-dimensional dynamic model of the tibiofemoral and patellofemoral articulations was developed to predict the motions of knee implants during a step-up activity. Patterns of muscle activity, initial joint angles and velocities, and kinematics of the hip and tinkle were measured experimentally and used as inputs to the simulation. Prosthetic knee kinematics were determined by integration of dynamic equations of motion subject to forces generated by muscles, ligaments, and contact at both the tibiofemoral and patellofemoral articulations. The modeling of contacts between implants did not rely upon explicit constraint equations; thus, changes in the number of contact points were allowed without modification to the model formulation. The simulation reproduced experimentally measured flexion-extension angle of the knee (within one standard deviation), but translations at the tibiofemoral articulations were larger during the simulated step-up task than those reported for patients with total knee replacements.     

Computational wear prediction of a total knee replacement from in vivo kinematics

Wear of ultra-high molecular weight polyethylene bearings in total knee replacements remains a major limitation to the longevity of these clinically successful devices. Few design tools are currently available to predict mild wear in implants based on varying kinematics, loads, and material properties. This paper reports the implementation of a computer modeling approach that uses fluoroscopically measured motions as inputs and predicts patient-specific implant damage using computationally efficient dynamic contact and tribological analyses. Multibody dynamic simulations of two activities (gait and stair) with two loading conditions (70-30 and 50-50 medial-lateral load splits) were generated from fluoroscopic data to predict contact pressure and slip velocity time histories for individual elements on the tibial insert surface. These time histories were used in a computational wear analysis to predict the depth of damage due to wear and creep experienced by each element. Predicted damage areas, volumes, and maximum depths were evaluated against a tibial insert retrieved from the same patient who provided the in vivo motions. Overall, the predicted damage was in close agreement with damage observed on the retrieval. The gait and stair simulations separately predicted the correct location of maximum damage on the lateral side, whereas a combination of gait and stair was required to predict the correct location on the medial side. Predicted maximum damage depths were consistent with the retrieval as well. Total computation time for each damage prediction was less than 30 min. Continuing refinement of this approach will provide a robust tool for accurately predicting clinically relevant wear in total knee replacements.     

Musculo-skeletal loading conditions at the hip during walking and stair climbing.

Musculo-skeletal loading plays an important role in the primary stability of joint replacements and in the biological processes involved in fracture healing. However, current knowledge of musculo-skeletal loading is still limited. In the past, a number of musculo-skeletal models have been developed to estimate loading conditions at the hip. So far, a cycle-to-cycle validation of predicted musculo-skeletal loading by in vivo measurements has not been possible. The aim of this study was to determine the musculo-skeletal loading conditions during walking and climbing stairs for a number of patients and compare these findings to in vivo data.Following total hip arthroplasty, four patients underwent gait analysis during walking and stair climbing. An instrumented femoral prosthesis enabled simultaneous measurement of in vivo hip contact forces. On the basis of CT and X-ray data, individual musculo-skeletal models of the lower extremity were developed for each patient. Muscle and joint contact forces were calculated using an optimization algorithm. The calculated peak hip contact forces both over- and under-estimated the measured forces. They differed by a mean of 12% during walking and 14% during stair climbing.For the first time, a cycle-to-cycle validation of predicted musculo-skeletal loading was possible for walking and climbing stairs in several patients. In all cases, the comparison of in vivo measured and calculated hip contact forces showed good agreement.Thus, the authors consider the presented approach as a useful means to determine valid conditions for the analysis of prosthesis loading, bone modeling or remodeling processes around implants and fracture stability following internal fixation.     

Sensitivity of medial and lateral knee contact force predictions to frontal plane alignment and contact locations

Musculoskeletal models are increasingly used to estimate medial and lateral knee contact forces, which are difficult to measure in vivo. The sensitivity of contact force predictions to modeling parameters is important to the interpretation and implication of results generated by the model. The purpose of this study was to quantify the sensitivity of knee contact force predictions to simultaneous errors in frontal plane knee alignment and contact locations under different dynamic conditions. We scaled a generic musculoskeletal model for N = 23 subjects’ stature and radiographic knee alignment, then perturbed frontal plane alignment and mediolateral contact locations within experimentally-possible ranges of 10° to −10° and 10 to −10 mm, respectively. The sensitivity of first peak, second peak, and mean medial and lateral knee contact forces to knee adduction angle and contact locations was modeled using linear regression. Medial loads increased, and lateral loads decreased, by between 3% and 6% bodyweight for each degree of varus perturbation. Shifting the medial contact point medially increased medial loads and decreased lateral loads by between 1% and 4% bodyweight per millimeter. This study demonstrates that realistic measurement errors of 5 mm (contact distance) or 5° (frontal plane alignment) could result in a combined 50% BW error in subject specific contact force estimates. We also show that model sensitivity varies between subjects as a result of differences in gait dynamics. These results demonstrate that predicted knee joint contact forces should be considered as a range of possible values determined by model uncertainty.     

Change in knee contact force with simulated change in body weight

The relationship between obesity, weight gain and progression of knee osteoarthritis is well supported, suggesting that excessive joint loading may be a mechanism responsible for cartilage deterioration. Examining the influence of weight gain on joint compressive forces is difficult, as both muscles and ground reaction forces can have a significant impact on the forces experienced during gait. While previous studies have examined the relationship between body weight and knee forces, these studies have used models that were not validated using experimental data. Therefore, the objective of this study was to evaluate the relationship between changes in body weight and changes in knee joint contact forces for an individual's gait pattern using musculoskeletal modeling that is validated against known internal compressive forces. Optimal weighting constants were determined for three subjects to generate valid predictions of knee contact forces (KCFs) using in vivo data collection with instrumented total knee arthroplasty. A total of five simulations per walking trial were generated for each subject, from 80% to 120% body weight in 10% increments, resulting in 50 total simulations. The change in peak KCF with respect to body weight was found to be constant and subject-specific, predominantly determined by the peak force during the baseline condition at 100% body weight. This relationship may be further altered by any change in kinematics or body mass distribution that may occur as a result of a change in body weight or exercise program.     

Simultaneous prediction of muscle and contact forces in the knee during gait

Musculoskeletal models are currently the primary means for estimating in vivo muscle and contact forces in the knee during gait. These models typically couple a dynamic skeletal model with individual muscle models but rarely include articular contact models due to their high computational cost. This study evaluates a novel method for predicting muscle and contact forces simultaneously in the knee during gait. The method utilizes a 12 degree-of-freedom knee model (femur, tibia, and patella) combining muscle, articular contact, and dynamic skeletal models. Eight static optimization problems were formulated using two cost functions (one based on muscle activations and one based on contact forces) and four constraints sets (each composed of different combinations of inverse dynamic loads). The estimated muscle and contact forces were evaluated using in vivo tibial contact force data collected from a patient with a force-measuring knee implant. When the eight optimization problems were solved with added constraints to match the in vivo contact force measurements, root-mean-square errors in predicted contact forces were less than 10 N. Furthermore, muscle and patellar contact forces predicted by the two cost functions became more similar as more inverse dynamic loads were used as constraints. When the contact force constraints were removed, estimated medial contact forces were similar and lateral contact forces lower in magnitude compared to measured contact forces, with estimated muscle forces being sensitive and estimated patellar contact forces relatively insensitive to the choice of cost function and constraint set. These results suggest that optimization problem formulation coupled with knee model complexity can significantly affect predicted muscle and contact forces in the knee during gait. Further research using a complete lower limb model is needed to assess the importance of this finding to the muscle and contact force estimation process.     

Hard-on-Hard Lubrication in the Artificial Hip under Dynamic Loading Conditions

The tribological performance of an artificial hip joint has a particularly strong influence on its success. The principle causes for failure are adverse short- and long-term reactions to wear debris and high frictional torque in the case of poor lubrication that may cause loosening of the implant. Therefore, using experimental and theoretical approaches models have been developed to evaluate lubrication under standardized conditions. A steady-state numerical model has been extended with dynamic experimental data for hard-on-hard bearings used in total hip replacements to verify the tribological relevance of the ISO 14242-1 gait cycle in comparison to experimental data from the Orthoload database and instrumented gait analysis for three additional loading conditions: normal walking, climbing stairs and descending stairs. Ceramic-on-ceramic bearing partners show superior lubrication potential compared to hard-on-hard bearings that work with at least one articulating metal component. Lubrication regimes during the investigated activities are shown to strongly depend on the kinematics and loading conditions. The outcome from the ISO gait is not fully confirmed by the normal walking data and more challenging conditions show evidence of inferior lubrication. These findings may help to explain the differences between the in vitro predictions using the ISO gait cycle and the clinical outcome of some hard-on-hard bearings, e.g., using metal-on-metal.     

Calibration of Optical Tweezers for In Vivo Force Measurements: How do Different Approaches Compare?

There is significant interest in quantifying force production inside cells, but since conditions in vivo are less well controlled than those in vitro, in vivo measurements are challenging. In particular, the in vivo environment may vary locally as far as its optical properties, and the organelles manipulated by the optical trap frequently vary in size and shape. Several methods have been proposed to overcome these difficulties. We evaluate the relative merits of these methods and directly compare two of them, a refractive index matching method, and a light-momentum-change method. Since in vivo forces are frequently relatively high (e.g., can exceed 15 pN for lipid droplets), a high-power laser is employed. We discover that this high-powered trap induces local temperature changes, and we develop an approach to compensate for uncertainties in the magnitude of applied force due to such temperature variations.     

Effect of swing phase load on metal-on-metal hip lubrication, friction and wear

There is renewed interest in metal-on-metal (MOM) total hip replacements (THRs), however, variable wear rates have been observed clinically. It is hypothesised that changes in soft tissue tensioning during surgery may alter loading of THRs during the swing phase of gait leading to changes in fluid film lubrication, friction and wear. This study aimed to assess the effect of swing phase load on the lubrication, friction and wear of MOM hip replacements. Theoretical lubrication modelling was carried out using elastohydrodynamic theory. All the governing equations were solved numerically for the lubricant film thickness between the articulating surfaces under the transient dynamic conditions with low and high swing phase loads. Friction testing was completed using a single axis pendulum simulator, simplified loading cycles were applied with low and high swing phase loads. MOM hip replacements were tested in a hip simulator, modified to provide different swing phase loading regimes; a low (100 N) and a high load (as per ISO 14242-1; 280 N). Results demonstrated that the performance of MOM bearings is highly dependent on swing phase load. Hence, changes in the tension of the tissues at surgery and variations in muscle forces may increase swing phase load, reduce lubrication, increase friction and accelerate wear. This may explain some of the variations that have been observed with clinical wear rates.     

Patellofemoral interactions in walking,stair ascent,and stair descent using a virtual patella model

Restoration of normal patella kinematics is an important clinical outcome of total knee arthroplasty. Failure of the patella within total knee systems has been documented and, upon occurrence, often necessitates revision surgery. It is thus important to understand patella mechanics following implantation, subject to load states that are typically realized during walking and other gaits. Here, a computational model of the patella is developed and used to examine the effects of walking, stair ascent, and stair descent on the development of stress and contact pressure in the patella throughout the gait cycle. Motion of the patella was governed by a combination of kinematic and force control, based on knee flexion and patellofemoral joint reaction force data from the literature. Unlike most previous analyses of full gait, quasi-static equilibrium was enforced throughout the cycle. Results indicate that, though peak forces vary greatly between the three gaits, maximum contact pressure and von Mises stress are roughly equivalent. However, contact area is larger in stair ascent and descent than walking, as patellofemoral loading, implant geometry, and polyethylene yield increase conformity between the femoral component and patella. Additionally, maximum contact pressure does not coincide with maximum load except for the case of walking. Though specific to the implant design considered here, this result has important ramifications for patella testing and emphasizes the need to characterize patella mechanics throughout gait.     

Effects of step length and step frequency on lower-limb muscle function in human gait

The aim of this study was to quantify the effects of step length and step frequency on lower-limb muscle function in walking. Three-dimensional gait data were used in conjunction with musculoskeletal modeling techniques to evaluate muscle function over a range of walking speeds using prescribed combinations of step length and step frequency. The body was modeled as a 10-segment, 21-degree-of-freedom skeleton actuated by 54 muscle-tendon units. Lower-limb muscle forces were calculated using inverse dynamics and static optimization. We found that five muscles – GMAX, GMED, VAS, GAS, and SOL – dominated vertical support and forward progression independent of changes made to either step length or step frequency, and that, overall, changes in step length had a greater influence on lower-limb joint motion, net joint moments and muscle function than step frequency. Peak forces developed by the uniarticular hip and knee extensors, as well as the normalized fiber lengths at which these muscles developed their peak forces, correlated more closely with changes in step length than step frequency. Increasing step length resulted in larger contributions from the hip and knee extensors and smaller contributions from gravitational forces (limb posture) to vertical support. These results provide insight into why older people with weak hip and knee extensors walk more slowly by reducing step length rather than step frequency and also help to identify the key muscle groups that ought to be targeted in exercise programs designed to improve gait biomechanics in older adults.     

Computational stability of human knee joint at early stance in Gait: Effects of muscle coactivity and anterior cruciate ligament deficiency

As one of the most complex and vulnerable structures of body, the human knee joint should maintain dynamic equilibrium and stability in occupational and recreational activities. The evaluation of its stability and factors affecting it is vital in performance evaluation/enhancement, injury prevention and treatment managements. Knee stability often manifests itself by pain, hypermobility and giving-way sensations and is usually assessed by the passive joint laxity tests. Mechanical stability of both the human knee joint and the lower extremity at early stance periods of gait (0% and 5%) were quantified here for the first time using a hybrid musculoskeletal model of the lower extremity. The roles of muscle coactivity, simulated by setting minimum muscle activation at 0–10% levels and ACL deficiency, simulated by reducing ACL resistance by up to 85%, on the stability margin as well as joint biomechanics (contact/muscle/ligament forces) were investigated. Dynamic stability was analyzed using both linear buckling and perturbation approaches at the final deformed configurations in gait. The knee joint was much more stable at 0% stance than at 5% due to smaller ground reaction and contact forces. Muscle coactivity, when at lower intensities (<3% of its maximum active force), increased dynamic stability margin. Greater minimum activation levels, however, acted as an ineffective strategy to enhance stability. Coactivation also substantially increased muscle forces, joint loads and ACL force and hence the risk of further injury and degeneration. A deficiency in ACL decreases total ACL force (by 31% at 85% reduced stiffness) and the stability margin of the knee joint at the heel strike. It also markedly diminishes forces in lateral hamstrings (by up to 39%) and contact forces on the lateral plateau (by up to 17%). Current work emphasizes the need for quantification of the lower extremity stability margin in gait.     

Tibio-femoral joint contact forces in sheep

Although the sheep has become a standard model for understanding the mechanical conditions that occur after injury and investigating surgical treatments such as osteochondral defect healing and ligament reconstruction, no study has yet evaluated the contact forces that occur in the sheep tibio-femoral joint in vivo. In this study, bone pins, together with reflective markers, were used to measure the 3D kinematics of three sheep hind limbs, simultaneously with the ground reaction forces during repetitions of gait trials. Joint contact forces were then calculated using inverse dynamics and optimisation techniques. Whilst average peak axial tibio-femoral contact forces of 2.1 body weight (BW) were calculated across the 3 sheep, only small medio-lateral and antero-posterior shear forces, averaging 0.7 BW, were determined. Average knee flexion angles ranging from 49 degrees to 70 degrees were observed. From the forces determined in this study, we have provided a better understanding of the mechanical loading environment that occurs in sheep. This has important implications for the interpretation of knee studies in quadrupeds and their relevance to the clinical situation.     

A Novel Approach for Dynamic Testing of Total Hip Dislocation under Physiological Conditions

Constant high rates of dislocation-related complications of total hip replacements (THRs) show that contributing factors like implant position and design, soft tissue condition and dynamics of physiological motions have not yet been fully understood. As in vivo measurements of excessive motions are not possible due to ethical objections, a comprehensive approach is proposed which is capable of testing THR stability under dynamic, reproducible and physiological conditions. The approach is based on a hardware-in-the-loop (HiL) simulation where a robotic physical setup interacts with a computational musculoskeletal model based on inverse dynamics. A major objective of this work was the validation of the HiL test system against in vivo data derived from patients with instrumented THRs. Moreover, the impact of certain test conditions, such as joint lubrication, implant position, load level in terms of body mass and removal of muscle structures, was evaluated within several HiL simulations. The outcomes for a normal sitting down and standing up maneuver revealed good agreement in trend and magnitude compared with in vivo measured hip joint forces. For a deep maneuver with femoral adduction, lubrication was shown to cause less friction torques than under dry conditions. Similarly, it could be demonstrated that less cup anteversion and inclination lead to earlier impingement in flexion motion including pelvic tilt for selected combinations of cup and stem positions. Reducing body mass did not influence impingement-free range of motion and dislocation behavior; however, higher resisting torques were observed under higher loads. Muscle removal emulating a posterior surgical approach indicated alterations in THR loading and the instability process in contrast to a reference case with intact musculature. Based on the presented data, it can be concluded that the HiL test system is able to reproduce comparable joint dynamics as present in THR patients.     

Forces and moments telemetered from two distal femoral replacements during various activities.

Two distal femoral replacements were instrumented to measure axial force, torque and bending moments in the prosthesis shaft. Data are reported up to 2.5 years for the following activities: uni- and bi-lateral standing, walking, stair climbing and descending, treadmill walking, jogging and jumping. In the first subject the greatest averaged peak shaft forces found were: jogging 3.6Bodyweight (BW), stair descending 3.1BW, walking 2.8BW, treadmill walking 2.75BW, and stair ascending 2.8BW. Bending moments about the antero-posterior axis (varus-valgus) and medio-lateral axis (flexion-extension) peaked in the range 8.5-9.8 and 4.7-7.6BWcm respectively, over the follow-up period. Axial torques peaked in the range 0.2-1.3BWcm, outwardly directed. At most follow-up sessions, forces and moments during jogging were generally greater than those for other gait activities. In the second subject forces and moments were generally only 45-70% of those in the first subject, due to inadequate musculature around the knee. The data can be applied to the design and testing of distal femoral replacements and even to total knee replacements, and contributes to the knowledge of forces acting in the distal femur during activity.     

A musculoskeletal finite element model of rat knee joint for evaluating cartilage biomechanics during gait

Abnormal loading of the knee due to injuries or obesity is thought to contribute to the development of osteoarthritis (OA). Small animal models have been used for studying OA progression mechanisms. However, numerical models to study cartilage responses under dynamic loading in preclinical animal models have not been developed. Here we present a musculoskeletal finite element model of a rat knee joint to evaluate cartilage biomechanical responses during a gait cycle. The rat knee joint geometries were obtained from a 3-D MRI dataset and the boundary conditions regarding loading in the joint were extracted from a musculoskeletal model of the rat hindlimb. The fibril-reinforced poroelastic (FRPE) properties of the rat cartilage were derived from data of mechanical indentation tests. Our numerical results showed the relevance of simulating anatomical and locomotion characteristics in the rat knee joint for estimating tissue responses such as contact pressures, stresses, strains, and fluid pressures. We found that the contact pressure and maximum principal strain were virtually constant in the medial compartment whereas they showed the highest values at the beginning of the gait cycle in the lateral compartment. Furthermore, we found that the maximum principal stress increased during the stance phase of gait, with the greatest values at midstance. We anticipate that our approach serves as a first step towards investigating the effects of gait abnormalities on the adaptation and degeneration of rat knee joint tissues and could be used to evaluate biomechanically-driven mechanisms of the progression of OA as a consequence of joint injury or obesity.     

A practical solution to reduce soft tissue artifact error at the knee using adaptive kinematic constraints

Musculoskeletal modeling and simulations have vast potential in clinical and research fields, but face various challenges in representing the complexities of the human body. Soft tissue artifact from skin-mounted markers may lead to non-physiological representation of joint motions being used as inputs to models in simulations. To address this, we have developed adaptive joint constraints on five of the six degree of freedom of the knee joint based on in vivo tibiofemoral joint motions recorded during walking, hopping and cutting motions from subjects instrumented with intra-cortical pins inserted into their tibia and femur. The constraint boundaries vary as a function of knee flexion angle and were tested on four whole-body models including four to six knee degrees of freedom. A musculoskeletal model developed in OpenSim simulation software was constrained to these in vivo boundaries during level gait and inverse kinematics and dynamics were then resolved. Statistical parametric mapping indicated significant differences (p < 0.05) in kinematics between bone pin constrained and unconstrained model conditions, notably in knee translations, while hip and ankle flexion/extension angles were also affected, indicating the error at the knee propagates to surrounding joints. These changes to hip, knee, and ankle kinematics led to measurable changes in hip and knee transverse plane moments, and knee frontal plane moments and forces. Since knee flexion angle can be validly represented using skin mounted markers, our tool uses this reliable measure to guide the five other degrees of freedom at the knee and provide a more valid representation of the kinematics for these degrees of freedom.     

In vivo contact kinematics and contact forces of the knee after total knee arthroplasty during dynamic weight-bearing activities

Analysis of polyethylene component wear and implant loosening in total knee arthroplasty (TKA) requires precise knowledge of in vivo articular motion and loading conditions. This study presents a simultaneous in vivo measurement of tibiofemoral articular contact forces and contact kinematics in three TKA patients. These measurements were accomplished via a dual fluoroscopic imaging system and instrumented tibial implants, during dynamic single leg lunge and chair rising-sitting. The measured forces and contact locations were also used to determine mediolateral distribution of axial contact forces. Contact kinematics data showed a medial pivot during flexion of the knee, for all patients in the study. Average axial forces were higher for lunge compared to chair rising-sitting (224% vs. 187% body weight). In this study, we measured peak anteroposterior and mediolateral forces averaging 13.3% BW during lunge and 18.5% BW during chair rising-sitting. Mediolateral distributions of axial contact force were both patient and activity specific. All patients showed equitable medial-lateral loading during lunge but greater loads at the lateral compartment during chair rising-sitting. The results of this study may enable more accurate reproduction of in vivo loads and articular motion patterns in wear simulators and finite element models. This in turn may help advance our understanding of factors limiting longevity of TKA implants, such as aseptic loosening and polyethylene component wear, and enable improved TKA designs.     

Pattern of anterior cruciate ligament force in normal walking

The goal of this study was to calculate and explain the pattern of anterior cruciate ligament (ACL) loading during normal level walking. Knee-ligament forces were obtained by a two-step procedure. First, a three-dimensional (3D) model of the whole body was used together with dynamic optimization theory to calculate body-segmental motions, ground reaction forces, and leg-muscle forces for one cycle of gait. Joint angles, ground reaction forces, and muscle forces obtained from the gait simulation were then input into a musculoskeletal model of the lower limb that incorporated a 3D model of the knee. The relative positions of the femur, tibia, and patella and the forces induced in the knee ligaments were found by solving a static equilibrium problem at each instant during the simulated gait cycle. The model simulation predicted that the ACL bears load throughout stance. Peak force in the ACL (303 N) occurred at the beginning of single-leg stance (i.e., contralateral toe off). The pattern of ACL force was explained by the shear forces acting at the knee. The balance of muscle forces, ground reaction forces, and joint contact forces applied to the leg determined the magnitude and direction of the total shear force acting at the knee. The ACL was loaded whenever the total shear force pointed anteriorly. In early stance, the anterior shear force from the patellar tendon dominated the total shear force applied to the leg, and so maximum force was transmitted to the ACL at this time. ACL force was small in late stance because the anterior shear forces supplied by the patellar tendon, gastrocnemius, and tibiofemoral contact were nearly balanced by the posterior component of the ground reaction.     

Direct comparison of measured and calculated total knee replacement force envelopes during walking in the presence of normal and abnormal gait patterns

Knee joint forces measured from instrumented implants provide important information for testing the validity of computational models that predict knee joint forces. The purpose of this study was to validate a parametric numerical model for predicting knee joint contact forces against measurements from four subjects with instrumented TKRs during the stance phase of gait. Model sensitivity to abnormal gait patterns was also investigated. The results demonstrated good agreement for three subjects with relatively normal gait patterns, where the difference between the mean measured and calculated forces ranged from 0.05 to 0.45 body weights, and the envelopes of measured and calculated forces (from three walking trials) overlapped. The fourth subject, who had a "quadriceps avoidance" external moment pattern, initially had little overlap between the measured and calculated force envelopes. When additional constraints were added, tailored to the subject's gait pattern, the model predictions improved to complete force envelope overlap. Coefficient of multiple determination analysis indicated that the shape of the measured and calculated force waveforms were similar for all subjects (adjusted coefficient of multiple correlation values between 0.88 and 0.92). The parametric model was accurate in predicting both the magnitude and waveform of the contact force, and the accuracy of model predictions was affected by deviations from normal gait patterns. Equally important, the envelope of forces generated by the range of solutions substantially overlapped with the corresponding measured envelope from multiple gait trials for a given subject, suggesting that the variable strategic processes of in vivo force generation are covered by the solution range of this parametric model.