Velocity-dependent cost function for the prediction of force sharing among synergistic muscles in a one degree of freedom model

Prediction of accurate and meaningful force sharing among synergistic muscles is a major problem in biomechanics research. Given a resultant joint moment, a unique set of muscle forces can be obtained from this mathematically redundant system using nonlinear optimization. The classical cost functions for optimization involve a normalization of the muscle forces to the absolute force capacity of the target muscles, usually by the cross-sectional area or the maximal isometric force. In a one degree of freedom model this leads to a functional relationship between moment arms and the predicted muscle forces, such that for constant moment arms, or constant ratios of moment arms, agonistic muscle forces increase or decrease in unison. Experimental studies have shown however that the relationship between muscle forces is highly task-dependent often causing forces to increase in one muscle while decreasing in a functional agonist, likely because of the contractile conditions and contractile properties of the involved muscles. We therefore, suggest a modified cost function that accounts for the instantaneous contraction velocity of the muscles and its effect on the instantaneous maximal force. With this novel objective function, a task-dependent prediction of muscle force distribution is obtained that allows, even in a one degree of freedom system, the prediction of force sharing loops, and simultaneously increasing and decreasing forces for agonist pairs of muscles.     

Muscular synergism--I. On criteria for load sharing between synergistic muscles

The use of optimization techniques to predict individual muscle forces in redundant biomechanical systems implies the formulation of a criterion for load sharing between the muscles. In part I of this paper, the characteristics and performance of several linear and non-linear criteria reported in the literature have been compared for static-isometric knee flexion. The results show that linear criteria inherently predict discrete muscle action (orderly recruitment of muscles) whereas non-linear criteria can predict synergistic action. All criteria predict that relatively more force is allocated to muscles with large moment arms. When muscle stresses (or ratios of muscle force to maximum muscle force) are used as the decision variables in the objective function, then relatively more force is allocated to muscles with large maximum possible force as well. Future formulations of the optimization should consider the differences in fiber type composition among the muscles. Such an approach is presented in part II of the paper.     

Electromyography and mechanics of mastication in the albino rat.

The masticatory apparatus in the albino rat was studied by means of electromyography and subsequent estimation of muscular forces. The activity patterns of the trigeminal and suprahyoid musculature and the mandibular movements were recorded simultaneously during feeding. The relative forces of the individual muscles in the different stages of chewing cycles and biting were estimated on the basis of their physiological cross sections and their activity levels, as measured from integrated electromyograms. Workinglines and moment arms of these muscles were determined for different jaw positions. In the anteriorly directed masticatory grinding stroke the resultants of the muscle forces at each side are identical; they direct anteriorly, dorsally and slightly lingually and pass along the lateral side of the second molar. Almost the entire muscular resultant force is transmitted to the molars while the temporo-mandibular joint remains unloaded. A small transverse force, produced by the tense symphyseal cruciate ligaments balances the couple of muscle resultant and molar reaction force in the transverse plane. After each grinding stroke the mandible is repositioned for the next stroke by the overlapping actions of three muscle groups: the pterygoids and suprahyoids produce depression and forward shift, the suprahyoids and temporal backward shift and elevation of the mandible while the subsequent co-operation of the temporal and masseter causes final closure of the mouth and starting of the forward grinding movement. All muscles act in a bilaterally symmetrical fashion. The pterygoids contract more strongly, the masseter more weakly during biting than during chewing. The wide gape shifts the resultant of the muscle forces more vertically and moreposteriorly. The joint then becomes strongly loaded because the reaction forces are applied far anteriorly on the incisors. The charateristic angle between the almost horizontal biting force and the surface of the food pellet indicates that the lower incisors produce a chisel-like action. Tooth structure reflects chewing and biting forces. The transverse molar lamellae lie about parallel to the chewing forces whereas perpendicular loading of the occlusal surfaces is achieved by their inclination in the transverse plane. The incisors are loaded approximately parallel to their longitudinal axis, placement that avoids bending forces during biting. It is suggested that a predominantly protrusive musculature favors the effective force transmission to the lower incisors, required for gnawing. By grinding food across transversely oriented molar ridges the protrusive components of the muscles would be utilized best. From the relative weights of the masticatory muscles in their topographical relations with joints, molars and incisors it may be concluded that the masticatory apparatus is a construction adapted to optimal transmission of force from muscles to teeth.     

Relationships between the size and spatial morphology of human masseter and medial pterygoid muscles, the craniofacial skeleton, and jaw biomechanics

The relationship between human craniofacial morphology and the biomechanical efficiency of bite force generation in widely varying muscular and skeletal types is unknown. To address this problem, we selected 22 subjects with different facial morphologies and used magnetic resonance imaging, cephalometric radiography, and data from dental casts to reconstruct their craniofacial tissues in three dimensions. Conventional cephalometric analyses were carried out, and the cross-sectional sizes of the masseter and medial pterygoid muscles were measured from reconstituted sections. The potential abilities of the muscles to generate bite forces at the molar teeth and mandibular condyles were calculated according to static equilibrium theory using muscle, first molar, and condylar moment arms. On average, the masseter muscle was about 66% larger in cross section than the medial pterygoid and was inclined more anteriorly relative to the functional occlusal plane. There was a significant positive correlation (P less than 0.01) between the cross-sectional areas of the masseter and medial pterygoid muscles (r = 0.75) and between the bizygomatic arch width and masseter cross-sectional area (r = 0.56) and medial pterygoid cross-sectional area (r = 0.69). The masseter muscle was always a more efficient producer of vertically oriented bite force than the medial pterygoid. Putative bite force from the medial pterygoid muscle alone correlated positively with mandibular length and inversely with upper face height. When muscle and tooth moment arms were considered together, a system efficient at producing force on the first molar was statistically associated with a face having a large intergonial width, small intercondylar width, narrow dental arch, forward maxilla, and forward mandible. There was no significant correlation between muscle cross-sectional areas and their respective putative bite forces. This suggests that there is no simple relationship between the tension-generating capacity of the muscles and their mechanical efficiency as described by their spatial arrangement. The study shows that in a modern human population so many combinations of biomechanically relevant variables are possible that subjects cannot easily be placed into ideal or nonideal categories for producing molar force. Our findings also confirm the impression that similar bite-force efficiencies can be found in subjects with disparate facial features.     

Individual muscle force parameters and fiber operating ranges for elbow flexion-extension and forearm pronation-supination

We have quantified individual muscle force and moment contributions to net joint moments and estimated the operating ranges of the individual muscle fibers over the full range of motion for elbow flexion/extension and forearm pronation/supination. A three dimensional computer graphics model was developed in order to estimate individual muscle contributions in each degree of freedom over the full range of motion generated by 17 muscles crossing the elbow and forearm. Optimal fiber length, tendon slack length, and muscle specific tension values were adjusted within the literature range from cadaver studies such that the net isometric joint moments of the model approximated experimental joint moments within one standard deviation. Analysis of the model revealed that the muscles operate on varying portions of the ascending limb, plateau region, and descending limb of the force-length curve. This model can be used to further understand isometric force and moment contributions of individual muscles to net joint moments of the arm and forearm and can serve as a comprehensive reference for the forces and moments generated by 17 major muscles crossing the elbow and wrist.     

Effect of bone quality on the forces generated by compression screws.

Internal fixation of the fractured scaphoid bone is used to promote union between bone fragments and to decrease wrist immobilization. Headless screws are commonly used because they minimize interference with articular surfaces and reduce tissue irritation and immobilization. In the present experiment, compressive force was measured as a function of bone quality for two headless screw types, the Herbert and the Acutrak. Forty-seven cylindrical samples of cancellous bone were prepared from fresh, previously frozen human cadaveric distal femora. The compressive forces generated as the screws were advanced into the specimens were measured and correlated to the specimens' bone mineral density (BMD) and density. Over the range tested, the average compressive force for the Acutrack screw was approximately 42% higher than that of the Herbert. Statistical significance, however, could not established because of the low statistical power of the test due to the inherent spread in the data. For the Acutrak screw, force was best fit to BMD and to density by second-order polynomials. Regression analysis indicated that similar correlations did not exist between force and BMD or between force and density for the Herbert screw. The correlation shown by the Acutrak screw indicates that it may be a more predictable as well as more effective system and therefore there may be some advantage in selecting this system. Furthermore, results suggest that the Acutrak screw generates greater forces with increasing bone density and could be more effective for a younger population.     

Static optimization underestimates antagonist muscle activity at the glenohumeral joint: A musculoskeletal modeling study

Static optimization is commonly employed in musculoskeletal modeling to estimate muscle and joint loading; however, the ability of this approach to predict antagonist muscle activity at the shoulder is poorly understood. Antagonist muscles, which contribute negatively to a net joint moment, are known to be important for maintaining glenohumeral joint stability. This study aimed to compare muscle and joint force predictions from a subject-specific neuromusculoskeletal model of the shoulder driven entirely by measured muscle electromyography (EMG) data with those from a musculoskeletal model employing static optimization. Four healthy adults performed six sub-maximal upper-limb contractions including shoulder abduction, adduction, flexion, extension, internal rotation and external rotation. EMG data were simultaneously measured from 16 shoulder muscles using surface and intramuscular electrodes, and joint motion evaluated using video motion analysis. Muscle and joint forces were calculated using both a calibrated EMG-driven neuromusculoskeletal modeling framework, and musculoskeletal model simulations that employed static optimization. The EMG-driven model predicted antagonistic muscle function for pectoralis major, latissimus dorsi and teres major during abduction and flexion; supraspinatus during adduction; middle deltoid during extension; and subscapularis, pectoralis major and latissimus dorsi during external rotation. In contrast, static optimization neural solutions showed little or no recruitment of these muscles, and preferentially activated agonistic prime movers with large moment arms. As a consequence, glenohumeral joint force calculations varied substantially between models. The findings suggest that static optimization may under-estimate the activity of muscle antagonists, and therefore, their contribution to glenohumeral joint stability.     

Biomechanics of the masticatory apparatus of Pteropus giganteus (Megachiroptera)

Jaw mechanics in Pteropus were studied by means of a three-dimensional model. The model included several parameters of muscle architecture, combined with quantified movement and electromyographical data. Estimates of the nature of the applied forces that act upon the mandible during a chewing cycle, and subsequent estimates of reaction forces at the bite point and joints during the powerstroke, were thus obtained for different food consistencies. The resultant muscle force (relative to the palate) shifts from upward and slightly backward at large gapes to upward and markedly backward at the end of closing. The resultant simultaneously moves anteriorly. During the powerstroke it retains a constant position and orientation along the thickened anterior edge of the coronoid process. The early stages of opening are guided by the slope of the teeth and mandibular fossa; during the remaining part of opening the working line of the resultant crosses the skull behind the joint and thus acquires an opening moment. The bite force has downward and forward components, and a slight transverse component. For a given applied muscular force its magnitude is larger in more posteriorly positioned bite points. Both joints are loaded, the contralateral one more than the ipsilateral. Food consistency affects magnitude and orientation of the applied force, and hence, magnitude and orientation of the bite force and magnitude of the joint reaction forces. The magnitude of masseter activity relative to temporalis activity appears to be the key factor for the orientation of the bite force, and hence for the mechanical optimal position of the food. The adaptive value of the general topography of the masticatory muscles in Pteropus is discussed.     

Balancing function of the masticatory muscles during incisal biting in two murid rodents,Apodemus speciosus and Clethrionomys rufocanus

The functional significance of masticatory muscle direction was estimated using a mechanical model in two murid rodents: the Japanese field mouse (Apodemus speciosus) and the gray red-backed vole (Clethrionomys rufocanus). Theoretical analyses of the data suggest that a balancing mechanism among the muscle forces occurs during incisal power stroke. The activation of the large deep masseter in both murids results in marked tensile separation of two hemimandibles at the flexible mandibular symphysis. Activation of the internal pterygoid decreases this large tensile force at the symphysis more efficiently than other muscles. The lines of action of the deep masseter and internal pterygoid are aligned to produce such a balancing function in both species studied here. The resultant force generated by the deep masseter on both sides is opposite in direction to the reaction force at the lower incisor tip. Therefore, the large deep masseter forms an effective mandibular support mechanism when the reaction forces during biting push the mandible downward. Because of the area of insertion and the line of action, the posterior temporalis appears to have an important role in stabilizing the position of the mandibular condyle in the glenoid fossa during incisal biting. J. Morphol. 236:49–56, 1998. © 1998 Wiley-Liss, Inc.     

Sensitivity of model predictions of muscle function to changes in moment arms and muscle-tendon properties: a Monte-Carlo analysis

Hill-type muscle models are commonly used in musculoskeletal models to estimate muscle forces during human movement. However, the sensitivity of model predictions of muscle function to changes in muscle moment arms and muscle-tendon properties is not well understood. In the present study, a three-dimensional muscle-actuated model of the body was used to evaluate the sensitivity of the function of the major lower limb muscles in accelerating the whole-body center of mass during gait. Monte-Carlo analyses were used to quantify the effects of entire distributions of perturbations in the moment arms and architectural properties of muscles. In most cases, varying the moment arm and architectural properties of a muscle affected the torque generated by that muscle about the joint(s) it spanned as well as the torques generated by adjacent muscles. Muscle function was most sensitive to changes in tendon slack length and least sensitive to changes in muscle moment arm. However, the sensitivity of muscle function to changes in moment arms and architectural properties was highly muscle-specific; muscle function was most sensitive in the cases of gastrocnemius and rectus femoris and insensitive in the cases of hamstrings and the medial sub-region of gluteus maximus. The sensitivity of a muscle's function was influenced by the magnitude of the muscle's force as well as the operating region of the muscle on its force-length curve. These findings have implications for the development of subject-specific models of the human musculoskeletal system.     

Predicted pattern of human muscle activity during clenching derived from a computer assisted model: symmetric vertical bite forces

A computer assisted three-dimensional model of the jaw, based on linear programming, is presented. The upper and lower attachments of the muscles of mastication have been measured on a single human skull and divided into thirteen independent units on each side--a total of 26 muscle elements. The direction (in three dimensions) and maximum forces that could be developed by each muscle element, the bite reaction and two joint reactions are included in the model. It is shown for symmetrical biting that a model which minimizes the sum of the muscle forces used to produce a given bite force activates muscles in a way which corresponds well with previous observations on human subjects. A model which minimizes the joint reactions behaves differently and is rejected. An analysis of the way the chosen model operates suggests that there are two types of jaw muscles, power muscles and control muscles. Power muscles (superficial masseter, medial pterygoid and some of temporalis) produce the bite force but tend to displace the condyle up or down the articular eminence. This displacement is prevented by control muscles (oblique temporalis and lateral pterygoid) which have very poor moment arms for generating usual bite forces, but are efficient for preventing condylar slide. The model incorporates the concept that muscles consist of elements which can contract independently. It predicts that those muscle elements with longer moment arms relative to the joint are the first to be activated and, as the bite force increases, a ripple of activity spreads into elements with shorter moment arms. In general, the model can be used to study the three-dimensional activity in any system of joints and muscles.     

Prediction of pedal forces in bicycling using optimization methods

The bicycle-rider system is modeled as a planar five-bar linkage with pedal forces and pedal dynamics as input. The pedal force profile input is varied, maintaining constant average bicycle power, in order to obtain the optimal pedal force profile that minimizes two cost functions. One cost function is based on joint moments and the other is based on muscle stresses. Predicted (optimal) pedal profiles as well as joint moment time histories are compared to representative real data to examine cost function appropriateness. Both cost functions offer reasonable predictions of pedal forces. The muscle stress cost function, however, better predicts joint moments. Predicted muscle activity also correlates well with myoelectric data. The factors that lead to effective (i.e. low cost) pedalling are examined. Pedalling effectiveness is found to be a complex function of pedal force vector orientation and muscle mechanics.     

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.     

Sensitivity of predicted muscle forces to parameters of the optimization-based human leg model revealed by analytical and numerical analyses

There are different opinions in the literature on whether the cost functions: the sum of muscle stresses squared and the sum of muscle stresses cubed, can reasonably predict muscle forces in humans. One potential reason for the discrepancy in the results could be that different authors use different sets of model parameters which could substantially affect forces predicted by optimization-based models. In this study, the sensitivity of the optimal solution obtained by minimizing the above cost functions for a planar three degrees-of-freedom (DOF) model of the leg with nine muscles was investigated analytically for the quadratic function and numerically for the cubic function. Analytical results revealed that, generally, the non-zero optimal force of each muscle depends in a very complex non-linear way on moments at all three joints and moment arms and physiological cross-sectional areas (PCSAs) of all muscles. Deviations of the model parameters (moment arms and PCSAs) from their nominal values within a physiologically feasible range affected not only the magnitude of the forces predicted by both criteria, but also the number of non-zero forces in the optimal solution and the combination of muscles with non-zero predicted forces. Muscle force magnitudes calculated by both criteria were similar. They could change several times as model parameters changed, whereas patterns of muscle forces were typically not as sensitive. It is concluded that different opinions in the literature about the behavior of optimization-based models can be potentially explained by differences in employed model parameters.     

About weight factors in the non-linear objective functions used for solving indeterminate problems in biomechanics.

A lot of non-linear objective criteria are applied for solving the indeterminate problems formulated for different biomechanical models--most of them can be covered by the expression [formula in text]. It might be noted, however, that most of the suggested criteria are not applicable if considerable antagonistic co-contractions exist. This could be an effect of treating the agonistic muscles and their respective antagonists in one and the same manner in the objective function. Using a completely inverse approach (the muscle forces are supposed to be known quantities) and a simple 1DOF model (actuated by three agonistic muscles and one corresponding antagonist) it has been shown which values of the weight factors c(i) may predict different levels of muscle forces from the two antagonistic groups. Three hypothetical border variants for magnitudes of the muscle forces are considered (flexor muscles are only active, extensor muscles are only active, considerable co-contraction of flexors and extensors exists). The main conclusions are: the signs of c(i) at agonistic muscles have to be opposite to the c(i) signs at their antagonists; the signs of the weight factors depend on the direction of the net external joint moment; the closer c(i) to zero, the bigger force will be predicted in the ith muscle.     

Evaluation of the influence of muscle deactivation on other muscles and joints during gait motion

When any muscle in the human musculoskeletal system is damaged, other muscles and ligaments tend to compensate for the role of the damaged muscle by exerting extra effort. It is beneficial to clarify how the roles of the damaged muscles are compensated by other parts of the musculoskeletal system from the following points of view: From a clinical point of view, it will be possible to know how the abnormal muscle and joint forces caused by the acute compensations lead to further physical damage to the musculoskeletal system. From the viewpoint of rehabilitation, it will be possible to know how the role of the damaged muscle can be compensated by extra training of the other muscles. A method to evaluate the influence of muscle deactivation on other muscles and joints is proposed in this report. Methodology based on inverse dynamics and static optimization, which is applicable to arbitrary motion was used in this study. The evaluation method was applied to gait motion to obtain matrices representing (1) the dependence of muscle force compensation and (2) the change to bone-on-bone contact forces. These matrices make it possible to evaluate the effects of deactivation of one of the muscles of the musculoskeletal system on the forces exerted by other muscles as well as the change to the bone-on-bone forces when the musculoskeletal system is performing the same motion. Through observation of this matrix, it was found that deactivation of a muscle often results in increment/decrement of force developed by muscles with completely different primary functions and bone-on-bone contact force in different parts of the body. For example, deactivation of the iliopsoas leads to a large reduction in force by the soleus. The results suggest that acute deactivation of a muscle can result in damage to another part of the body. The results also suggest that the whole musculoskeletal system must go through extra retraining in the case of damage to certain muscles.     

Critical analysis of musculoskeletal modelling complexity in multibody biomechanical models of the upper limb

The inverse dynamics technique applied to musculoskeletal models, and supported by optimisation techniques, is used extensively to estimate muscle and joint reaction forces. However, the solutions of the redundant muscle force sharing problem are sensitive to the detail and modelling assumptions of the models used. This study presents four alternative biomechanical models of the upper limb with different levels of discretisation of muscles by bundles and muscle paths, and their consequences on the estimation of the muscle and joint reaction forces. The muscle force sharing problem is solved for the motions of abduction and anterior flexion, acquired using video imaging, through the minimisation of an objective function describing muscle metabolic energy consumption. While looking for the optimal solution, not only the equations of motion are satisfied but also the stability of the glenohumeral and scapulothoracic joints is preserved. The results show that a lower level of muscle discretisation provides worse estimations regarding the muscle forces. Moreover, the poor discretisation of muscles relevant to the joint in analysis limits the applicability of the biomechanical model. In this study, the biomechanical model of the upper limb describing the infraspinatus by a single bundle could not solve the complete motion of anterior flexion. Despite the small differences in the magnitude of the forces predicted by the biomechanical models with more complex muscular systems, in general, there are no significant variations in the muscular activity of equivalent muscles.     

Simple and complex models for studying muscle function in walking

While simple models can be helpful in identifying basic features of muscle function, more complex models are needed to discern the functional roles of specific muscles in movement. In this paper, two very different models of walking, one simple and one complex, are used to study how muscle forces, gravitational forces and centrifugal forces (i.e. forces arising from motion of the joints) combine to produce the pattern of force exerted on the ground. Both the simple model and the complex one predict that muscles contribute significantly to the ground force pattern generated in walking; indeed, both models show that muscle action is responsible for the appearance of the two peaks in the vertical force. The simple model, an inverted double pendulum, suggests further that the first and second peaks are due to net extensor muscle moments exerted about the knee and ankle, respectively. Analyses based on a much more complex, muscle-actuated simulation of walking are in general agreement with these results; however, the more detailed model also reveals that both the hip extensor and hip abductor muscles contribute significantly to vertical motion of the centre of mass, and therefore to the appearance of the first peak in the vertical ground force, in early single-leg stance. This discrepancy in the model predictions is most probably explained by the difference in model complexity. First, movements of the upper body in the sagittal plane are not represented properly in the double-pendulum model, which may explain the anomalous result obtained for the contribution of a hip-extensor torque to the vertical ground force. Second, the double-pendulum model incorporates only three of the six major elements of walking, whereas the complex model is fully 3D and incorporates all six gait determinants. In particular, pelvic list occurs primarily in the frontal plane, so there is the potential for this mechanism to contribute significantly to the vertical ground force, especially during early single-leg stance when the hip abductors are activated with considerable force.     

A comparison of the triceps surae and residual muscle moments at the ankle during cycling

The rigid linked system model and principles of inverse dynamics have been widely used to calculate residual muscle moments during various activities. EMG driven models and optimization algorithms have also been presented in the literature in efforts to estimate skeletal muscle forces and evaluate their possible contribution to the residual muscle moment. Additionally, skeletal muscle-tendon forces have been measured, directly, in both animals and humans. The purpose of this investigation was to calculate the moment produced by the triceps surae muscles and compare it to the residual muscle moment at the ankle during cycling at three power outputs (90, 180 and 270 W). Inferences were made regarding the potential contribution made by each triceps surae component to the tendon force using EMG and muscle-tendon length changes. A buckle-type transducer was surgically implanted on the right Achilles tendon of one male subject. Achilles tendon forces measured in vivo were multiplied by their corresponding moment arms to yield the triceps surae moment during the three working conditions. Moment arm lengths were obtained in a separate experiment using magnetic resonance imaging (MRI). Pedal reaction forces, body segment accelerations (determined from high speed film), and appropriate mass parameters served as input to the inverse solution. The triceps surae moment was temporally in phase with and consistently represented approximately 65% of the residual muscle moment at the ankle. These data demonstrate the feasibility of using implanted transducers in human subjects and provide a greater understanding of musculoskeletal mechanics during normal human movements.