Investigation of the Peculiarities of Two-joint Muscles Using a 3 DOF Model of the Human Upper Limb in the Sagittal Plane: an Optimization Approach

Abstract

The purpose of this paper is an investigation of the peculiarities of biarticular muscles by means of modelling and analytical solution of the indeterminate problem. The basic model includes 10 muscle elements performing flexio/extensio in the shoulder, elbow and wrist. Four of them are biarticular muscles. Two modifications of the model with only monoarticular muscles are developed. The indeterminate problem is solved analytically using the objective criterion σc_iF_i ² where F( is the module of the i-th muscle force and Cj is a weight coefficient. The predicted muscle forces, joint reactions and moments are compared in-between the basic model and its two modifications for different joint angles, external loading and weight coefficients. The main conclusions are: it is impossible to formulate strict advantages of the biarticular muscles under quasistatical conditions, their peculiarities depend on limb position, external loading and neural control; in general, monoarticular muscles are more powerful than biarticular ones; the biarticular muscles fine tune muscle coordination, their control is more precise and graceful; the presence of biarticular muscles leads to an increase of the joint reactions and moments, thus stabilizing the limb.     

The Roles of Monoarticular and Biarticular Muscles of the Lower Limbs in Terrestrial Locomotion

Specific features of the functioning of mono- and biarticular muscles were studied using a multijoint movement (a high jump) as an example. The powers of the knee and ankle joint extensors are insufficient for a strong and quick movement such as a high jump. Biarticular muscles (m. rectus femoris) transfer forces/powers from one joint to another, thereby compensating for the physiological shortcoming of monoarticular muscles, that is, a decrease in the tractive force with increasing contraction rate. In a high jump, a power of 300 W may be transferred from the hip to the knee joint via the m. rectus femoris; 230 W, from the knee to the hip joint via the hamstring muscle; 210 W, from the knee joint to the ankle via the m. gastrocnemius; and 15 W, from the metatarsophalangeal joint to the ankle via the mm. flexors.     

Determination of muscle loading at the hip joint for use in pre-clinical testing

The stability of joint endoprostheses depends on the loading conditions to which the implant-bone complex is exposed. Due to a lack of appropriate muscle force data, less complex loading conditions tend to be considered in vitro. The goal of this study was to develop a load profile that better simulates the in vivo loading conditions of a "typical" total hip replacement patient and considers the interdependence of muscle and joint forces. The development of the load profile was based on a computer model of the lower extremities that has been validated against in vivo data. This model was simplified by grouping functionally similar hip muscles. Muscle and joint contact forces were computed for an average data set of up to four patients throughout walking and stair climbing. The calculated hip contact forces were compared to the average of the in vivo measured forces. The final derived load profile included the forces of up to four muscles at the instances of maximum in vivo hip joint loading during both walking and stair climbing. The hip contact forces differed by less than 10% from the peak in vivo value for a "typical" patient. The derived load profile presented here is the first that is based on validated musculoskeletal analyses and seems achievable in an in vitro test set-up. It should therefore form the basis for further standardisation of pre-clinical testing by providing a more realistic approximation of physiological loading conditions.     

Contributions of individual muscles to hip joint contact force in normal walking

The human hip joint withstands high contact forces during daily activity and is therefore susceptible to injury and structural deterioration over time. Knowledge of muscle-force contributions to hip joint loading may assist in the development of strategies to prevent and manage conditions such as osteoarthritis, femoro-acetabular impingement and fracture. The main aim of this study was to determine the contributions of individual muscles to hip contact force in normal walking. Muscle contributions to hip contact force were calculated based on a previously published dynamic optimization solution for normal walking, which provided the time histories of joint motion, ground reaction forces, and muscle forces during the stance and swing phases of gait. The force developed by each muscle plus its contribution to the ground reaction force were used to determine the muscle’s contribution to hip contact force. Muscles were the major contributors to hip contact force, with gravitational and centrifugal forces combined contributing less than 5% of the total contact force. Four muscles that span the hip – gluteus medius, gluteus maximus, iliopsoas, and hamstrings – contributed most significantly to the three components of the hip contact force and hip contact impulse (integral of hip contact force over time). Three muscles that do not span the hip – vasti, soleus, and gastrocnemius – also contributed substantially to hip joint loading. These results provide additional insight into lower-limb muscle function during walking and may also be relevant to studies of cartilage degeneration and bone remodelling at the hip.     

A comparison of optimisation methods and knee joint degrees of freedom on muscle force predictions during single-leg hop landings

The aim of this paper was to compare the effect of different optimisation methods and different knee joint degrees of freedom (DOF) on muscle force predictions during a single legged hop. Nineteen subjects performed single-legged hopping manoeuvres and subject-specific musculoskeletal models were developed to predict muscle forces during the movement. Muscle forces were predicted using static optimisation (SO) and computed muscle control (CMC) methods using either 1 or 3 DOF knee joint models. All sagittal and transverse plane joint angles calculated using inverse kinematics or CMC in a 1 DOF or 3 DOF knee were well-matched (RMS error<3°). Biarticular muscles (hamstrings, rectus femoris and gastrocnemius) showed more differences in muscle force profiles when comparing between the different muscle prediction approaches where these muscles showed larger time delays for many of the comparisons. The muscle force magnitudes of vasti, gluteus maximus and gluteus medius were not greatly influenced by the choice of muscle force prediction method with low normalised root mean squared errors (<48%) observed in most comparisons. We conclude that SO and CMC can be used to predict lower-limb muscle co-contraction during hopping movements. However, care must be taken in interpreting the magnitude of force predicted in the biarticular muscles and the soleus, especially when using a 1 DOF knee. Despite this limitation, given that SO is a more robust and computationally efficient method for predicting muscle forces than CMC, we suggest that SO can be used in conjunction with musculoskeletal models that have a 1 or 3 DOF knee joint to study the relative differences and the role of muscles during hopping activities in future studies.     

Influence of muscle activity on the forces in the femur: An in vivo study

Experiments were performed on two patients with custom-made instrumented massive proximal femoral prostheses implanted after tumour resection. In vivo axial forces transmitted along the prostheses were telemetered during level walking, single- and double-leg stance, and isometric exercises of the hip muscles. These activities varied the lever arms available to the external loads: minimum for double-leg stance and maximum for hip isometric exercises. Kinematic, force plate, EMG and telemetered force data were recorded simultaneously. The force magnification ration (FMR; the ratio of the telemetered axial force to the external force) was calculated. The FMRs ranged from 1.3 (during double-leg stance) to 29.8 (during abductors test), indicating that a major part of the axial force in the long bones is a response to muscle activity, the strength of which depends on the lever arms available to the external loads. From these results, it was shown that the bulk of the bending moment along limbs is transmitted by a combination of tensile forces in muscles and compressive forces in bones, so moments transmitted by the bones are smaller than the limb moments. It was concluded that appropriate simulation of muscle forces is important in experimental or theoretical studies of load transmission along bones.     

Activity of mono- and biarticular leg muscles during sprint running

A cinematographic recording of the movements of the lower limbs together with simultaneous emg tracings from nine lower limb muscles were obtained from two male track sprinters during three phases of a 100 m sprint run. The extensor muscles of the hip joint were found to be the primary movers by acceleration of the body's center of gravity (C.G.) during the ground phase of the running cycle. The extensors of the knee joint were also important in this, but to a minor extent, while the plantar flexors of the ankle joint showed the least contribution. The biarticular muscles functioned in a way different from the monoarticular muscles in the sense that they perform eccentric work during the flight and recovery phases and concentric work during the whole ground phase (support), whereas the monoarticular muscles are restricted first to eccentric work and then to concentric work during the ground phase. Furthermore, the biarticular muscles show variation (and rate of variation) in muscle length to a larger extent than the monoarticular muscles. Paradoxical muscle actions appear to take place around the knee joint, where the hamstring muscles, m. gastrocnemius, m. vastus laterialis and m. vastus medialis act as synergists by extending the knee joint during the last part of the ground phase.     

Contributions of muscle forces and toe-off kinematics to peak knee flexion during the swing phase of normal gait: an induced position analysis

A three-dimensional dynamic simulation of walking was used together with induced position analysis to determine how kinematic conditions at toe-off and muscle forces following toe-off affect peak knee flexion during the swing phase of normal gait. The flexion velocity of the swing-limb knee at toe-off contributed 30 degrees to the peak knee flexion angle; this was larger than any contribution from an individual muscle or joint moment. Swing-limb muscles individually made large contributions to knee angle (i.e., as large as 22 degrees), but their actions tended to balance one another, so that the combined contribution from all swing-limb muscles was small (i.e., less than 3 degrees of flexion). The uniarticular muscles of the swing limb made contributions to knee flexion that were an order of magnitude larger than the biarticular muscles of the swing limb. The results of the induced position analysis make clear the importance of knee flexion velocity at toe-off relative to the effects of muscle forces exerted after toe-off in generating peak knee flexion angle. In addition to improving our understanding of normal gait, this study provides a basis for analyzing stiff-knee gait, a movement abnormality in which knee flexion in swing is diminished.     

Quantitative functional anatomy of the lower limb with application to human gait

A scheme was developed to classify muscles according to their primary, secondary and tertiary functions, e.g. a muscle which produces primarily a flexion moment may also produce secondary abduction and tertiary internal rotation moments. The functions of muscles crossing the hip and knee joints were computed based upon the changing relative positions of joint centers and muscle origins and insertions during one gait cycle. The function of several of the major muscles crossing the hip and knee joints is reported for the different limb positions corresponding to normal gait. It was found that the amount of force necessary to produce a given moment about a joint was dependent upon the limb position. In addition, the muscle functions changed significantly with limb position. Electrical stimulation of muscles of a paralyzed subject gave qualitative support to the results.     

The effects of sloped surfaces on locomotion: backward walking as a perturbation

The purpose of this study was to examine lower extremity kinetics and muscle activity during backward slope walking to clarify the relationship between joint moments and powers and muscle activity patterns observed in forward slope walking. Nine healthy volunteers walked backward on an instrumented ramp at three grades (-39% (-21 degrees ), 0% (level), +39% (+21 degrees )). EMG activity was recorded from major lower extremity muscles. Joint kinetics were obtained from kinematic and force platform data. The knee joint moment and power generation increased significantly during upslope walking; hip joint moment and power absorption increased significantly during downslope walking. When compared to data from forward slope walking, these backward walking data suggest that power requirements of a task dictate the muscle activity pattern needed to accomplish that movement. During downslope walking tasks, power absorption increased and changes in muscle activity patterns were directly related to the changes in the joint moment patterns. In contrast, during upslope walking tasks, power generation increased and changes in the muscle activity were related to the changes in the joint moments only at the 'primary' joint; at adjacent joints the changes in muscle activity were unrelated to the joint moment pattern. The 'paradoxical' changes in the muscle activity at the adjacent joints are possibly related to the activation of biarticular muscles required by the increased power generation at the primary joint. In total, these data suggest that changing power requirements at a joint impact the control of muscle activity at that and adjacent joints.     

The force resulting from the action of mono- and biarticular muscles in a limb.

Human and animal limbs can be modelled as a chain of segments connected at joints. For a static limb, the force exerted at the endpoint due to the force of a single muscle has been calculated. It turns out that there are marked differences in the action of mono- vs. biarticular muscles. Monoarticular muscles produce an endpoint force that is directed in the lengthwise direction of the limb, i.e. in the direction of one of the segments. The force from biarticular muscles can have a marked transverse component. The 'principal direction' of this endpoint force is also the movement direction of the endpoint which is the most favourable for the muscle to do work. The reasoning presented can explain e.g. the differences in the activity of mono- and biarticular muscles in cycling.     

Contribution of mono- and biarticular muscles to extending knee joint moments in runners and cyclists.

Motor actions are governed by coordinated activation of mono- and biarticular muscles. This study considered differences in mono- and biarticular knee extensors between runners and cyclists in the context of adaptations to task-specific movement requirements. Two hypotheses were tested: 1) the length-at-use hypothesis, which is that muscle adapts to have it operate around optimal length; and 2) the contraction-mode hypothesis, which is that eccentrically active muscles prefer to operate on the ascending limb of the length-force curve. Ten runners and ten cyclists performed maximal, isometric knee extensions on a dynamometer at five knee and four hip joint angles. This approach allowed the separation of the contribution of mono- and biarticular extensors. Three major differences occurred: 1) compared with runners, monoarticular extensors of cyclists reach optimal length at larger muscle length; 2) in runners, optimal length of the biarticular extensor is shifted to larger lengths; and 3) the moment generated by monoarticular extensor was larger in cyclists. Mono- and biarticular extensors respond to different adaptation triggers in runners and cyclists. Monoarticular muscles seem to adapt to the length-at-use, whereas biarticular muscles were found to be sensitive to the contraction-mode hypothesis.     

A noninvasive technique for in vivo measurement of joint torques of biarticular muscles.

The objective of this work was to develop a noninvasive method to measure the joint torques produced by biarticular muscles at two joints simultaneously. During intramuscular stimulation of the cat medial gastrocnemius (MG) muscle, torques at the ankle and knee joints were calculated from forces measured in two dimensions at the end point of the cat paw under isometric conditions. The method was verified by the known anatomical properties of cat MG muscle and the tibialis anterior (TA) muscle. The MG muscle was shown to produce a significant flexion torque at the knee, besides an extension torque at the ankle. This was in agreement with its anatomical arrangement. The TA muscle produced primarily an ankle flexion torque. The small knee torque, due to measurement errors, yielded an estimate of measurement accuracy of 3.0 +/- 2.1% (n = 52). The coupling ratio of the MG muscle, defined as T(ankle)/T(knee), varied significantly with both knee and ankle angles. The profile of MG mechanical coupling agreed qualitatively with changes in limb configuration. The method can be used to measure recruitment properties of electrically stimulated biarticular muscles, and may potentially be used to study the biomechanics of biarticular coupling.     

Muscle contributions to mediolateral and anteroposterior foot placement during walking

Foot placement is critical to balance control during walking and is primarily controlled by muscle force generation. Although gluteus medius activity has been associated with mediolateral foot placement, how other muscles contribute to foot placement is not clear. Furthermore, although dynamic walking models have suggested that anteroposterior foot placement can be passively controlled, the extent to which muscles actively contribute to anteroposterior foot placement has not been determined. The objective of this study was to identify individual muscle contributions to mediolateral and anteroposterior foot placement during walking in healthy adults. Dynamic simulations of walking were developed for six older adults and a segmental power analysis was performed to determine the individual muscle contributions to the mediolateral and anteroposterior power delivered to the foot segment. The simulations revealed the ipsilateral swing limb gluteus medius, iliopsoas, rectus femoris and hamstrings and the contralateral stance limb gluteus medius and ankle plantarflexors were primary contributors to both mediolateral and anteroposterior foot placement. Muscle contributions to foot placement were found to be highly influenced by their contributions to pelvis power, which was dominated by those muscles crossing the hip joint. Thus, impaired balance control may be improved by focusing rehabilitation interventions on optimizing the coordination of those muscles crossing the hip joint and the ankle plantarflexors.     

Joint moments and contact forces in the foot during walking

The net force and moment of a joint have been widely used to understand joint disease in the foot. Meanwhile, it does not reflect the physiological forces on muscles and contact surfaces. The objective of the study is to estimate active moments by muscles, passive moments by connective tissues and joint contact forces in the foot joints during walking. Joint kinematics and external forces of ten healthy subjects (all males, 24.7 ± 1.2 years) were acquired during walking. The data were entered into the five-segment musculoskeletal foot model to calculate muscle forces and joint contact forces of the foot joints using an inverse dynamics-based optimization. Joint reaction forces and active, passive and net moments of each joint were calculated from muscle and ligament forces. The maximum joint reaction forces were 8.72, 4.31, 2.65, and 3.41 body weight (BW) for the ankle, Chopart’s, Lisfranc and metatarsophalangeal joints, respectively. Active and passive moments along with net moments were also obtained. The maximum net moments were 8.6, 8.4, 5.4 and 0.8%BW∙HT, respectively. While the trend of net moment was very similar between the four joints, the magnitudes and directions of the active and passive moments varied between joints. The active and passive moments during walking could reveal the roles of muscles and ligaments in each of the foot joints, which was not obvious in the net moment. This method may help narrow down the source of joint problems if applied to clinical studies.     

In vivo measurement of dynamic rectus femoris function at postures representative of early swing phase

Forward dynamic models suggest that muscle-induced joint motions depend on dynamic coupling between body segments. As a result, biarticular muscles may exhibit non-intuitive behavior in which the induced joint motion is opposite to that assumed based on anatomy. Empirical validation of such predictions is important for models to be relied upon to characterize muscle function. In this study, we measured, in vivo, the hip and knee accelerations induced by electrical stimulation of the rectus femoris (RF) and the vastus medialis (VM) at postures representatives of the toe-off and early swing phases of the gait cycle. Seven healthy young subjects were positioned side-lying with their lower limb supported on air bearings while a 90 ms pulse train stimulated each muscle separately or simultaneously. Lower limb kinematics were measured and compared to predictions from a similarly configured dynamic model of the lower limb. We found that both RF and VM, when stimulated independently, accelerated the hip and knee into extension at these postures, consistent with model predictions. Predicted ratios of hip acceleration to knee acceleration were generally within 1 s.d. of average values. In addition, measured responses to simultaneous RF and VM stimulation were within 13% of predictions based on the assumption that joint accelerations induced by activating two muscles simultaneously can be found by adding the joint accelerations induced by activating the same muscles independently. These results provide empirical evidence of the importance of considering dynamic effects when interpreting the role of muscles in generating movement.     

Multi-functionality of the cat medical gastrocnemius during locomotion

The functional role of biarticular muscles was investigated based on direct force measurement in the cat medial gastrocnemius (MG) and analysis of hindlimb kinematics and kinetics for the stance phase of level, uphill, and downhill walking. Four primary functional roles of biarticular muscles have been proposed in the past. These functional roles have typically been discussed independently of each other, and biarticular muscles have rarely been assigned more than one functional roles for different phases of the work cycle. The purpose of this study was to elucidate the functional role of the biarticular cat MG during locomotion. It was found that MG forces were primarily associated with the moment requirements at the ankle for most of the stance phase, but also helped to satisfy the moments at the knee in the initial phase of stance. In the second half of stance, MG transferred mechanical energy from the knee to the ankle from the knee to the ankle, while simultaneously producing a substantial amount of mechanical work. Based on these results, we hypothesize that MG's primary function is that of an ankle extensor. However, because of the coupling of the ankle extensor moment with a knee flexor moment in the initial, and a knee extensor moment in the final phase of stance, MG satisfies two joint moments in early stance, and transfers mechanical energy from the knee to the ankle in late stance. We conclude that cat MG has multiple functional roles during the stance phase of locomotion, and speculate that such multi-functionality also exists in other bi- and multi-articular muscles.     

Individual muscle contributions to the axial knee joint contact force during normal walking

Muscles are significant contributors to the high joint forces developed in the knee during human walking. Not only do muscles contribute to the knee joint forces by acting to compress the joint, but they also develop joint forces indirectly through their contributions to the ground reaction forces via dynamic coupling. Thus, muscles can have significant contributions to forces at joints they do not span. However, few studies have investigated how the major lower-limb muscles contribute to the knee joint contact forces during walking. The goal of this study was to use a muscle-actuated forward dynamics simulation of walking to identify how individual muscles contribute to the axial tibio-femoral joint force. The simulation results showed that the vastii muscles are the primary contributors to the axial joint force in early stance while the gastrocnemius is the primary contributor in late stance. The tibio-femoral joint force generated by these muscles was at times greater than the muscle forces themselves. Muscles that do not cross the knee joint (e.g., the gluteus maximus and soleus) also have significant contributions to the tibio-femoral joint force through their contributions to the ground reaction forces. Further, small changes in walking kinematics (e.g., knee flexion angle) can have a significant effect on the magnitude of the knee joint forces. Thus, altering walking mechanics and muscle coordination patterns to utilize muscle groups that perform the same biomechanical function, yet contribute less to the knee joint forces may be an effective way to reduce knee joint loading during walking.     

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.     

Hip joint contact force in the emu (Dromaius novaehollandiae) during normal level walking

The emu is a large, (bipedal) flightless bird that potentially can be used to study various orthopaedic disorders in which load protection of the experimental limb is a limitation of quadrupedal models. An anatomy-based analysis of normal emu walking gait was undertaken to determine hip contact forces for comparison with human data. Kinematic and kinetic data captured for two laboratory-habituated emus were used to drive the model. Muscle attachment data were obtained by dissection, and bony geometries were obtained by CT scan. Inverse dynamics calculations at all major lower-limb joints were used in conjunction with optimization of muscle forces to determine hip contact forces. Like human walking gait, emu ground reaction forces showed a bimodal distribution over the course of the stance phase. Two-bird averaged maximum hip contact force was approximately 5.5 times body weight, directed nominally axially along the femur. This value is only modestly larger than optimization-based hip contact forces reported in literature for humans. The interspecies similarity in hip contact forces makes the emu a biomechanically attractive animal in which to model loading-dependent human orthopaedic hip disorders.