Musculotendon lengths and moment arms for a three-dimensional upper-extremity model

Generating muscle-driven forward dynamics simulations of human movement using detailed musculoskeletal models can be computationally expensive. This is due in part to the time required to calculate musculotendon geometry (e.g., musculotendon lengths and moment arms), which is necessary to determine and apply individual musculotendon forces during the simulation. Modeling upper-extremity musculotendon geometry can be especially challenging due to the large number of multi-articular muscles and complex muscle paths. To accurately represent this geometry, wrapping surface algorithms and/or other computationally expensive techniques (e.g., phantom segments) are used. This paper provides a set of computationally efficient polynomial regression equations that estimate musculotendon length and moment arms for thirty-two (32) upper-extremity musculotendon actuators representing the major muscles crossing the shoulder, elbow and wrist joints. Equations were developed using a least squares fitting technique based on geometry values obtained from a validated public-domain upper-extremity musculoskeletal model that used wrapping surface elements (Holzbaur et al., 2005). In general, the regression equations fit well the original model values, with an average root mean square difference for all musculotendon actuators over the represented joint space of 0.39 mm (1.1% of peak value). In addition, the equations reduced the computational time required to simulate a representative upper-extremity movement (i.e., wheelchair propulsion) by more than two orders of magnitude (315 versus 2.3 s). Thus, these equations can assist in generating computationally efficient forward dynamics simulations of a wide range of upper-extremity movements.     

Estimation of musculotendon kinematics in large musculoskeletal models using multidimensional B-splines

We present a robust and computationally inexpensive method to estimate the lengths and three-dimensional moment arms for a large number of musculotendon actuators of the human lower limb. Using a musculoskeletal model of the lower extremity, a set of values was established for the length of each musculotendon actuator for different lower limb generalized coordinates (joint angles). A multidimensional spline function was then used to fit these data. Muscle moment arms were obtained by differentiating the musculotendon length spline function with respect to the generalized coordinate of interest. This new method was then compared to a previously used polynomial regression method. Compared to the polynomial regression method, the multidimensional spline method produced lower errors for estimating musculotendon lengths and moment arms throughout the whole generalized coordinate workspace. The fitting accuracy was also less affected by the number of dependent degrees of freedom and by the amount of experimental data available. The spline method only required information on musculotendon lengths to estimate both musculotendon lengths and moment arms, thus relaxing data input requirements, whereas the polynomial regression requires different equations to be used for both musculotendon lengths and moment arms. Finally, we used the spline method in conjunction with an electromyography driven musculoskeletal model to estimate muscle forces under different contractile conditions, which showed that the method is suitable for the integration into large scale neuromusculoskeletal models.     

Moment arms and musculotendon lengths estimation for a three-dimensional lower-limb model

This paper presents a set of polynomial expressions that can be used as regression equations to estimate length and three-dimensional moment arms of 43 lower-limb musculotendon actuators. These equations allow one to find, at a low computational cost, the musculotendon geometric parameters required for numerical simulation of large musculoskeletal models. Nominal values for these biomechanical parameters were established using a public-domain musculoskeletal model of the lower limb (IEEE Trans. Biomed. Eng. 37 (1990) 757). To fit these nominal values, regression equations with different levels of complexity were generated, based on the number of generalized coordinates of the joints spanned by each musculotendon actuator. Least squares fitting was used to identify regression equation coefficients. The goodness of the fit and confidence intervals were assessed, and the best fitting equations selected.     

Biarticular hip extensor and knee flexor muscle moment arms of the feline hindlimb

Moment arms are important for understanding muscular behavior and for calculating internal muscle forces in musculoskeletal simulations. Biarticular muscles cross two joints and have moment arms that depend on the angle of both joints the muscles cross. The tendon excursion method was used to measure the joint angle-dependence of hamstring (biceps femoris, semimembranosus and semitendinosus) moment arm magnitudes of the feline hindlimb at the knee and hip joints. Knee angle influenced hamstring moment arm magnitudes at the hip joint; compared to a flexed knee joint, the moment arm for semimembranosus posterior at the hip was at most 7.4 mm (25%) larger when the knee was extended. On average, hamstring moment arms at the hip increased by 4.9 mm when the knee was more extended. In contrast, moment arm magnitudes at the knee varied by less than 2.8 mm (mean=1.6 mm) for all hamstring muscles at the two hip joint angles tested. Thus, hamstring moment arms at the hip were dependent on knee position, while hamstring moment arms at the knee were not as strongly associated with relative hip position. Additionally, the feline hamstring muscle group had a larger mechanical advantage at the hip than at the knee joint.     

Evaluation of different analytical methods for subject-specific scaling of musculotendon parameters

Musculoskeletal models are often used to estimate internal muscle forces and the effects of those forces on the development of human movement. The Hill-type muscle model is an important component of many of these models, yet it requires specific knowledge of several muscle and tendon properties. These include the optimal muscle fibre length, the length at which the muscle can generate maximum force, and the tendon slack length, the length at which the tendon starts to generate a resistive force to stretch. Both of these parameters greatly influence the force-generating behaviour of a musculotendon unit and vary with the size of the person. However, these are difficult to measure directly and are often estimated using the results of cadaver studies, which do not account for differences in subject size. This paper examined several different techniques that can be used to scale the optimal muscle fibre length and tendon slack length of a musculotendon unit according to subject size. The techniques were divided into three categories corresponding to linear scaling, scaling by maintaining a constant tendon slack length throughout the range of joint motion, and scaling by maintaining muscle operating range throughout the range of joint motion. We suggest that a good rationale for scaling muscle properties should be to maintain the same force-generating characteristics of a musculotendon unit for all subjects, which is best achieved by scaling that preserves the muscle operating range when the muscle is maximally activated.     

Architectural properties of distal forelimb muscles in horses,Equus caballus

Articular injuries in athletic horses are associated with large forces from ground impact and from muscular contraction. To accurately and noninvasively predict muscle and joint contact forces, a detailed model of musculoskeletal geometry and muscle architecture is required. Moreover, muscle architectural data can increase our understanding of the relationship between muscle structure and function in the equine distal forelimb. Muscle architectural data were collected from seven limbs obtained from five thoroughbred and thoroughbred-cross horses. Muscle belly rest length, tendon rest length, muscle volume, muscle fiber length, and pennation angle were measured for nine distal forelimb muscles. Physiological cross-sectional area (PCSA) was determined from muscle volume and muscle fiber length. The superficial and deep digital flexor muscles displayed markedly different muscle volumes (227 and 656 cm3, respectively), but their PCSAs were very similar due to a significant difference in muscle fiber length (i.e., the superficial digital flexor muscle had very short fibers, while those of the deep digital flexor muscle were relatively long). The ulnaris lateralis and flexor carpi ulnaris muscles had short fibers (17.4 and 18.3 mm, respectively). These actuators were strong (peak isometric force, Fmax=5,814 and 4,017 N, respectively) and stiff (tendon rest length to muscle fiber length, LT:LMF=5.3 and 2.1, respectively), and are probably well adapted to stabilizing the carpus during the stance phase of gait. In contrast, the flexor carpi radialis muscle displayed long fibers (89.7 mm), low peak isometric force (Fmax=555 N), and high stiffness (LT:LMF=1.6). Due to its long fibers and low Fmax, flexor carpi radialis appears to be better adapted to flexion and extension of the limb during the swing phase of gait than to stabilization of the carpus during stance. Including muscle architectural parameters in a musculoskeletal model of the equine distal forelimb may lead to more realistic estimates not only of the magnitudes of muscle forces, but also of the distribution of forces among the muscles crossing any given joint.     

Fibre operating lengths of human lower limb muscles during walking

Muscles actuate movement by generating forces. The forces generated by muscles are highly dependent on their fibre lengths, yet it is difficult to measure the lengths over which muscle fibres operate during movement. We combined experimental measurements of joint angles and muscle activation patterns during walking with a musculoskeletal model that captures the relationships between muscle fibre lengths, joint angles and muscle activations for muscles of the lower limb. We used this musculoskeletal model to produce a simulation of muscle-tendon dynamics during walking and calculated fibre operating lengths (i.e. the length of muscle fibres relative to their optimal fibre length) for 17 lower limb muscles. Our results indicate that when musculotendon compliance is low, the muscle fibre operating length is determined predominantly by the joint angles and muscle moment arms. If musculotendon compliance is high, muscle fibre operating length is more dependent on activation level and force-length-velocity effects. We found that muscles operate on multiple limbs of the force-length curve (i.e. ascending, plateau and descending limbs) during the gait cycle, but are active within a smaller portion of their total operating range.     

Longer moment arm results in smaller joint moment development, power and work outputs in fast motions

Effects of moment arm length on kinetic outputs of a musculoskeletal system (muscle force development, joint moment development, joint power output and joint work output) were evaluated using computer simulation. A skeletal system of the human ankle joint was constructed: a lower leg segment and a foot segment were connected with a hinge joint. A Hill-type model of the musculus soleus (m. soleus), consisting of a contractile element and a series elastic element, was attached to the skeletal system. The model of the m. soleus was maximally activated, while the ankle joint was plantarflexed/dorsiflexed at a variation of constant angular velocities, simulating isokinetic exercises on a muscle testing machine. Profiles of the kinetic outputs (muscle force development, joint moment development, joint power output and joint work output) were obtained. Thereafter, the location of the insertion of the m. soleus was shifted toward the dorsal/ventral direction by 1cm, which had an effect of lengthening/shortening the moment arm length, respectively. The kinetic outputs of the musculoskeletal system during the simulated isokinetic exercises were evaluated with these longer/shorter moment arm lengths. It was found that longer moment arm resulted in smaller joint moment development, smaller joint power output and smaller joint work output in the larger plantarflexion angular velocity region (>120 degrees/s). This is because larger muscle shortening velocity was required with longer moment arm to achieve a certain joint angular velocity. Larger muscle shortening velocity resulted in smaller muscle force development because of the force-velocity relation of the muscle. It was suggested that this phenomenon should be taken into consideration when investigating the joint moment-joint angle and/or joint moment-joint angular velocity characteristics of experimental data.     

A Neuro-Muscular Elasto-Dynamic Model of the Human Arm Part 2: Musculotendon Dynamics and Related Stress Effects

In this paper we develop an elasto-dynamic model of the human arm that includes effects of neuro-muscular control uponelastic deformation in the limb.The elasto-dynamic model of the arm is based on hybrid parameter multiple body systemvariational projection principles presented in the companion paper.Though the technique is suitable for detailed bone and jointmodeling,we present simulations for simplified geometry of the bones,discretized as Rayleigh beams with elongation,whileallowing for large deflections.Motion of the upper extremity is simulated by incorporating muscle forces derived from aHill-type model of musculotendon dynamics.The effects of muscle force are modeled in two ways.In one approach,aneffective joint torque is calculated by multiplying the muscle force by a joint moment ann.A second approach models themuscle as acting along a straight line between the origin and insertion sites of the tendon.Simple arm motion is simulated byutilizing neural feedback and feedforward control.Simulations illustrate the combined effects of neural control strategies,models of muscle force inclusion,and elastic assumptions on joint trajectories and stress and strain development in the bone andtendon.     

Approximating complex musculoskeletal biomechanics using multidimensional autogenerating polynomials

Computational models of the musculoskeletal system are scientific tools used to study human movement, quantify the effects of injury and disease, plan surgical interventions, or control realistic high-dimensional articulated prosthetic limbs. If the models are sufficiently accurate, they may embed complex relationships within the sensorimotor system. These potential benefits are limited by the challenge of implementing fast and accurate musculoskeletal computations. A typical hand muscle spans over 3 degrees of freedom (DOF), wrapping over complex geometrical constraints that change its moment arms and lead to complex posture-dependent variation in torque generation. Here, we report a method to accurately and efficiently calculate musculotendon length and moment arms across all physiological postures of the forearm muscles that actuate the hand and wrist. Then, we use this model to test the hypothesis that the functional similarities of muscle actions are embedded in muscle structure. The posture dependent muscle geometry, moment arms and lengths of modeled muscles were captured using autogenerating polynomials that expanded their optimal selection of terms using information measurements. The iterative process approximated 33 musculotendon actuators, each spanning up to 6 DOFs in an 18 DOF model of the human arm and hand, defined over the full physiological range of motion. Using these polynomials, the entire forearm anatomy could be computed in <10 μs, which is far better than what is required for real-time performance, and with low errors in moment arms (below 5%) and lengths (below 0.4%). Moreover, we demonstrate that the number of elements in these autogenerating polynomials does not increase exponentially with increasing muscle complexity; complexity increases linearly instead. Dimensionality reduction using the polynomial terms alone resulted in clusters comprised of muscles with similar functions, indicating the high accuracy of approximating models. We propose that this novel method of describing musculoskeletal biomechanics might further improve the applications of detailed and scalable models to describe human movement.     

Inevitable joint angular rotation affects muscle architecture during isometric contraction

The purpose of this study was to quantify the influence of inevitable ankle joint motion during an isometric contraction on the measured change of the gastrocnemius medialis muscle (GM) architecture in vivo during the loading and the unloading phase. Sitting on a dynamometer subjects performed isometric maximal voluntary contractions as well as contractions induced by electrostimulation. Synchronous joint angular motion, plantarflexion moment, foot’s centre of pressure and real-time ultrasonography of muscle architecture changes of the GM were obtained. During the contraction the ankle joint position altered and significantly affected the change in muscle architecture. At maximal tendon force (1094 ± 323 N), the measured fascicle length overestimated the change in fascicle length due to the tendon force by 1.53 cm, while the measured pennation angle overestimated the change in pennation angle due to the tendon force by 5.5°. At the same tendon force the measured fascicle length and pennation angle were significantly different between loading and unloading conditions. After correcting the values for the change in ankle joint angle no differences between the loading and the unloading phase at the same tendon force were found. Concerning the estimation of GM fascicle length–force and pennation angle–force curves during the loading and unloading phase of an isometric contraction, these findings indicate that not accounting for ankle joint motion will produce unreliable results.     

Control of a one-link arm by burst signal generators

 The focus of this paper is the study of stability and point-to-point movement of a one-link arm. The sagittal arm has two musculotendon actuators, two neural oscillators that generate burst signals as motoneuron inputs, and spindles and Golgi tendon organs for extrinsic reflex feedbacks. It is shown that coactivation leads to intrinsic position and velocity feedback, and that the tendons introduce intrinsic force and rate of force feedback. In addition, the integrating effects of the tendons are studied when the actuator is constructed from a large number of identical fibers that are excited by alpha signals whose arrival times at the fiber are randomly distributed. Each of the musculotendon actuators receives two input signals – a burst signal analogous to alpha inputs and a conventional analogue signal that represents fusimotor input to the spindles. The process of combining burst signals and conventional analogue signals is studied. Simulation results show that the movement of the system with burst signals as inputs has overshoot and speed similar to the system with analogue signals. Received: 30 May 1994/Accepted in revised form: 13 January 1995     

Measurement of human Gracilis muscle isometric forces as a function of knee angle,intraoperatively

The goals of the present study were (1) to measure the previously unstudied isometric forces of activated human Gracilis (G) muscle as a function of knee joint angle and (2) to test whether length history effects are important also for human muscle. Experiments were conducted intraoperatively during anterior cruciate ligament (ACL) reconstruction surgery (n=8). Mean peak G muscle force, mean peak G tendon stress and mean optimal knee angle equals 178.5±270.3 N, 24.4±20.6 MPa and 67.5±41.7°, respectively. The substantial inter-subject variability found (e.g., peak G force ranges between 17.2 and 490.5 N) indicate that the contribution of the G muscle to knee flexion moment may vary considerably among subjects. Moreover, typical subject anthropometrics did not appear to provide a sound estimate of the peak G force: only a limited insignificant correlation was found between peak G force and subject mass as well as mid-thigh perimeter and no correlation was found between peak G force and thigh length. The functional joint range of motion for human G muscle was determined to be at least as wide as full knee extension to 120° of knee flexion. However; the portion of the knee angle–muscle force relationship operationalized is not unique but individual specific: our data suggest for most subjects that G muscle operates in both ascending and descending limbs of its length–force characteristics whereas, for the remainder of the subjects, its function is limited to the descending limb, exclusively. Previous activity of G muscle at high muscle length attained during collection of a complete set of knee angle–force data showed for the first time that such length history effects are important also for human muscles: a significant correlation was found between optimal knee angle and absolute value of % force change. Except for two of the subjects, G muscle force measured at low length was lower than that measured during collection of knee joint–force data (maximally by 42.3%).     

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.     

Influence of muscle anatomical cross-sectional area on the moment arm length of the triceps brachii muscle at the elbow joint

The purpose of this study was to test the hypothesis that the musculotendon moment arm length is affected by the muscle anatomical cross-sectional area. The moment arm length of the triceps brachii (TB) muscle at 30°, 50°, 70°, 90°, 110° elbow flexion positions was measured in sagittal magnetic resonance images (MRI) of 18 subjects as the perpendicular distance between the center of the pulley of the humerus to the line through the center of the TB tendon. The moment arm increased as the elbow flexion angle decreased, from 1.74±0.13 cm at 110° to 2.39±0.14 cm at 30°. The maximal anatomical cross-sectional area of the TB muscle was significantly correlated with the moment arms at all joint positions (r=0.545–0.803, p<0.05). Furthermore, the circumference of the upper arm was also significantly correlated with the moment arms at all joint positions, except for 70° (r=0.504–0.702, p<0.05). These results indicate that the moment arm length of the TB muscle is affected by the muscle anatomical cross-sectional area.     

Impact of ankle foot orthosis stiffness on Achilles tendon and gastrocnemius function during unimpaired gait

Ankle foot orthoses (AFOs) are designed to improve gait for individuals with neuromuscular conditions and have also been used to reduce energy costs of walking for unimpaired individuals. AFOs influence joint motion and metabolic cost, but how they impact muscle function remains unclear. This study investigated the impact of different stiffness AFOs on medial gastrocnemius muscle (MG) and Achilles tendon (AT) function during two walking speeds. We performed gait analyses for eight unimpaired individuals. Each individual walked at slow and very slow speeds with a 3D printed AFO with no resistance (free hinge condition) and four levels of ankle dorsiflexion stiffness: 0.25 Nm/°, 1 Nm/°, 2 Nm/°, and 3.7 Nm/°. Motion capture, ultrasound, and musculoskeletal modeling were used to quantify MG and AT lengths with each AFO condition. Increasing AFO stiffness increased peak AFO dorsiflexion moment with decreased peak knee extension and peak ankle dorsiflexion angles. Overall musculotendon length and peak AT length decreased, while peak MG length increased with increasing AFO stiffness. Peak MG activity, length, and velocity significantly decreased with slower walking speed. This study provides experimental evidence of the impact of AFO stiffness and walking speed on joint kinematics and musculotendon function. These methods can provide insight to improve AFO designs and optimize musculotendon function for rehabilitation, performance, or other goals.     

Chimpanzee thumb muscle cross sections, moment arms and potential torques, and comparisons with humans.

This study investigates the morphological basis of differences between humans and chimpanzees in the kinematical and dynamical parameters of the musculature of the thumb. It is partly intended to test an hypothesis that human thumb muscles can exert significantly greater torques, due to larger muscle cross-sectional areas or to longer tendon moment arms or to both. We focus on the estimation of the potentials of thumb muscles to exert torques about joint axes in a sample of eight chimpanzee cadaver hands. The potential torque of a muscle is estimated by taking the product of a muscle's physiological cross-sectional area (an estimator of force) with its dynamical moment arm (derived from the slope of tendon excursion versus joint angular displacement, obtained during passive movements of cadaver thumb joints). Comparison of our results with similar data obtained for humans at the same Mayo Clinic laboratory shows significant differences between humans and chimpanzees in potential torque of most thumb muscles, those of humans generally exhibiting larger values. The primary reason for the larger torques in humans is that their average moment arms are significantly longer, permitting greater torque for a given muscle size. An additional finding is that chimpanzees and humans differ in the direction of secondary thumb metacarpal movements elicited by contraction of some muscles, as shown by differences in moment arm signs for a given movement in the same muscle. The differences appear to be related to differences in the musculo-skeletal structures of the trapeziometacarpal joint.     

In vivo changes in the human patellar tendon moment arm length with different modes and intensities of muscle contraction

The purpose of this study was to examine the effect of different muscle contraction modes and intensities on patellar tendon moment arm length (d(PT)). Five men performed isokinetic concentric, eccentric and passive knee extensions at an angular velocity of 60 deg/s and six men performed gradually increasing to maximum effort isometric muscle contractions at 90( composite function) and 20( composite function) of knee flexion. During the tests, lateral X-ray fluoroscopy imaging was used to scan the knee joint. The d(PT) differences between the passive state and the isokinetic concentric and extension were quantified at 15( composite function) intervals of knee joint flexion angle. Furthermore, the changes of the d(PT) as a function of the isometric muscle contraction intensities were determined during the isometric knee extension at 90( composite function) and 20( composite function) of knee joint flexion. Muscle contraction-induced changes in knee joint flexion angle during the isometric muscle contraction were also taken into account for the d(PT) measurements. During the two isometric knee extensions, d(PT) increased from rest to maximum voluntary muscle contraction (MVC) by 14-15%. However, when changes in knee joint flexion angle induced by the muscle contraction were taken into account, d(PT) during MVC increased by 6-26% compared with rest. Moreover, d(PT) increased during concentric and eccentric knee extension by 3-15%, depending on knee flexion angle, compared with passive knee extension. These findings have important implications for estimating musculoskeletal loads using modelling under static and dynamic conditions.