Error propagation from kinematic data to modeled muscle-tendon lengths during walking

Kinematic data from 3D gait analysis together with musculoskeletal modeling techniques allow the derivation of muscle-tendon lengths during walking. However, kinematic data are subject to soft tissue artifacts (STA), referring to skin marker displacements during movement. STA are known to significantly affect the computation of joint kinematics, and would therefore also have an effect on muscle-tendon lengths which are derived from the segmental positions. The present study aimed to introduce an analytical approach to calculate the error propagation from STA to modeled muscle-tendon lengths. Skin marker coordinates were assigned uncorrelated, isotropic error functions with given standard deviations accounting for STA. Two different musculoskeletal models were specified; one with the joints moving freely in all directions, and one with the joints constrained to rotation but no translation. Using reference kinematic data from two healthy boys (mean age 9 y 5 m), the propagation of STA to muscle-tendon lengths was quantified for semimembranosus, gastrocnemius and soleus. The resulting average SD ranged from 6% to 50% of the normalized muscle-tendon lengths during gait depending on the muscle, the STA magnitudes and the musculoskeletal model. These results highlight the potential impact STA has on the biomechanical analysis of modeled muscle-tendon lengths during walking, and suggest the need for caution in the clinical interpretation of muscle-tendon lengths derived from joint kinematics.     

Finger Muscle Attachments for an OpenSim Upper-Extremity Model

We determined muscle attachment points for the index, middle, ring and little fingers in an OpenSim upper-extremity model. Attachment points were selected to match both experimentally measured locations and mechanical function (moment arms). Although experimental measurements of finger muscle attachments have been made, models differ from specimens in many respects such as bone segment ratio, joint kinematics and coordinate system. Likewise, moment arms are not available for all intrinsic finger muscles. Therefore, it was necessary to scale and translate muscle attachments from one experimental or model environment to another while preserving mechanical function. We used a two-step process. First, we estimated muscle function by calculating moment arms for all intrinsic and extrinsic muscles using the partial velocity method. Second, optimization using Simulated Annealing and Hooke-Jeeves algorithms found muscle-tendon paths that minimized root mean square (RMS) differences between experimental and modeled moment arms. The partial velocity method resulted in variance accounted for (VAF) between measured and calculated moment arms of 75.5% on average (range from 48.5% to 99.5%) for intrinsic and extrinsic index finger muscles where measured data were available. RMS error between experimental and optimized values was within one standard deviation (S.D) of measured moment arm (mean RMS error = 1.5 mm < measured S.D = 2.5 mm). Validation of both steps of the technique allowed for estimation of muscle attachment points for muscles whose moment arms have not been measured. Differences between modeled and experimentally measured muscle attachments, averaged over all finger joints, were less than 4.9 mm (within 7.1% of the average length of the muscle-tendon paths). The resulting non-proprietary musculoskeletal model of the human fingers could be useful for many applications, including better understanding of complex multi-touch and gestural movements.     

Talocrural and talocalcaneal joint kinematics and kinetics during the stance phase of walking

The purpose of this investigation was to study the kinematics and kinetics of the joints between the leg and calcaneus during the stance phase of walking. The talocrural and talocalcaneal joints were each assumed to act as monocentric single degree of freedom hinge joints. Motion at one joint was defined by the relative rotation of a point on the opposing joint. The results, based upon the gait of three subjects, showed that the hinge joint assumption may be reasonable. A discrepancy in the kinematics was shown between the talocrural joint rotation and its commonly assumed sagittal plane representation, especially during initial flatfoot. This discrepancy is due to the fact that the sagittal plane rotation is created by the combined rotations of the talocrural and talocalcaneal joints. The talocalcaneal joint showed a peak 25-30 Nm supinatory moment at 80% of stance. The talocrural joint moment was qualitatively similar to the commonly measured sagittal plane moment, but the present results show that the sagittal plane moment overpredicted the true moment by 6-22% due to the two-dimensional assumption.     

Estimation of cat medial gastrocnemius fascicle lengths during dynamic contractions

In typical muscle models, it is often assumed that the contractile element (fascicle) length depends exclusively on the instantaneous muscle-tendon length and the instantaneous muscle force. In order to test whether the instantaneous fascicle length during dynamic contractions can be predicted from muscle-tendon length and force, fascicle lengths, muscle-tendon lengths, and muscle forces were directly measured in cat medial gastrocnemii during isometric and dynamic contractions. Two theoretical muscle models were developed: model A was based on force-time data obtained during the activation phase and model D on force-time data obtained during the deactivation phase of isometric contractions. To test the models, instantaneous fascicle lengths were predicted from muscle-tendon lengths and forces during dynamic contractions that simulated cat locomotion for speeds ranging from 0.4 to 1.6m/s. The theoretically predicted fascicle lengths were compared with the experimentally measured fascicle lengths. It was found that fascicle lengths were not uniquely associated with muscle-tendon lengths and forces; that is, for a given muscle-tendon length and force, fascicle lengths varied depending on the contractile history. Consequently, models A and D differed in fascicle length predictions; model D (maximum average error=8.5%) was considerably better than model A (maximum average error=22.3%). We conclude from this study that it is not possible to predict the exact fascicle lengths from muscle-tendon lengths and forces alone, however, adequate predictions seem possible based on such a model. The relationship between fascicle length and muscle force and muscle-tendon length is complex and highly non-linear, thus, it appears unlikely that accurate fascicle length predictions can be made without some reference contractions in which fascicle length, muscle-tendon length, and force are measured simultaneously.     

Mechanical behavior of the quadriceps femoris muscle tendon unit during low-load contractions.

We examined the relationships between morphology and muscle-tendon dynamics of the quadriceps femoris muscle of 11 men using velocity-encoded phase-contrast magnetic resonance imaging (MRI). Thigh muscle electromyography and joint range of motion were first measured outside the MRI scanner during knee extension-flexion tasks that were performed at a rate of 40 times/min with elastic bands providing peak resistance of 5.2 kp (SD 0.4) to the extension. The same movement was repeated inside the MRI scanner bore where tissue velocities and muscle morphology were recorded. The average displacement in the proximal and distal halves of the rectus femoris and vastus intermedius aponeuroses was different (P = 0.049), reflecting shortening (1.6%), but the tensile strain along the length of the aponeuroses was uniform. The aponeurosis behavior varied among individuals, and these individual patterns were best explained by the differences in relative cross-sectional area of rectus femoris to vastus muscles (r = 0.71, P = 0.014). During dynamic contraction, considerable deformation of muscles in the axial plane caused an anatomic measure such as muscle thickness to change differently (decrease or increase) in different sites of measurement. For example, when analyzed from the axial images, the vastus lateralis thickness did not change (P = 0.946) in the frontal plane through femur but increased in a 45 degrees oblique plane between the frontal and sagittal planes (P = 0.004). The present observations of the heterogeneity and individual behavior emphasize the fact that single-point measurements do not always reflect the overall behavior of muscle-tendon unit.     

Shoulder joint kinematics during elevation measured by ultrasound-based measuring system.

In order to analyze shoulder joint movements, the authors use a ZEBRIS CMS-HS ultrasound-based movement analysis system. In essence, the measurement involves the determination of the spatial position of the 16 anatomical points, which are specified on the basis of the coordinates of ultrasound-based triplets positioned on the upper limb, the scapula, and the thorax; their spatial position is measured in the course of motion. Kinematic characteristics of 74 shoulder joints of 50 healthy persons were identified during elevation in the plane of the scapula. Kinematic characteristics of motion were identified by scapulothoracic, glenohumeral, and humeral elevation angles; range of angles; scapulothoracis and glenohumeral rhythm; scapulothoracic, glenohumeral, and scapuloglenoid ratios; and the relative displacement between the rotation centers of the humerus and the scapula. Motion of the humerus and the scapula relative to each other was characterized by their rotation as well as the relative displacement between the rotation centers of scapula and humerus. The biomechanical model of the shoulder joint during elevation can be described by analyzing the results of the measurements performed.     

Moving muscle points provide accurate curved muscle paths in a model of the cervical spine

Muscle paths in musculoskeletal models have been modeled using several different methods; however, deformation of soft tissue with changes in posture is rarely accounted for, and often only the neutral posture is used to define a muscle path. The objective of this study was to model curved muscle paths in the cervical spine that take into consideration soft tissue deformation with changes in neck posture. Two subject-specific models were created from magnetic resonance images (MRI) in 5 different sagittal plane neck postures. Curved paths of flexor and extensor muscles were modeled using piecewise linear lines-of-action in two ways; (1) using fixed via points determined from muscle paths in the neutral posture and (2) using moving muscle points that moved relative to the bones determined from muscle paths in all 5 postures. Accuracy of each curved modeled muscle path was evaluated by an error metric, the distance from the anatomic (centroid) muscle path determined from the MRI. Error metric was compared among three modeled muscle path types (straight, fixed via and moving muscle point) using a repeated measures one-way ANOVA (α=0.05). Moving muscle point paths had 21% lower error metric than fixed via point paths over all 15 pairs of neck muscles examined over 5 postures (3.86 mm vs. 4.88 mm). This study highlights the importance of defining muscle paths in multiple postures in order to properly define the changing curvature of a muscle path due to soft tissue deformation with posture.     

A light weight compliant hand mechanism with high degrees of freedom

This paper presents the design and prototyping of an inherently compliant lightweight hand mechanism. The hand mechanism itself has 15 degrees of freedom and five fingers. Although the degrees of freedom in each finger are coupled, reducing the number of independent degrees of freedom to 5, the 15 degrees of freedom of the hand could potentially be individually actuated. Each joint consists of a novel flexing mechanism that is based on the loading of a compression spring in the axial and transverse direction via a cable and conduit system. Currently, a bench top version of the prototype is being developed; the three joints of each finger are coupled together to simplify the control system. The current control scheme under investigation simulates a control scheme where myoelectric signals in the wrist flexor and extensor muscles are converted in to x and y coordinates on a control scheme chart. Static load-deformation analysis of finger segments is studied based on a 3-dimensional model without taking the stiffener into account, and the experiment validates the simulation.     

Three-dimensional model of the feline hindlimb

This article describes a three-dimensional musculoskeletal model of the feline hindlimb based on digitized musculoskeletal anatomy. The model consists of seven degrees of freedom: three at the hip and two each at the knee and ankle. Lines of action and via points for 32 major muscles of the limb are described. Interspecimen variability of muscle paths was surprisingly low; most via points displayed a scatter of only a few millimeters. Joint axes identified by mechanical techniques as noncoincident and nonorthogonal were further honed to yield moment arms consistent with previous reports. Interspecimen variability in joint axes was greater than that of muscle paths and highlights the importance of joint axes in kinematic models. The contribution of specific muscles to the direction of endpoint force generation is discussed.     

Biomechanical behavior of muscle-tendon complex during dynamic human movements

This paper reviews the research findings regarding the force and length changes of the muscle-tendon complex during dynamic human movements, especially those using ultrasonography and computer simulation. The use of ultrasonography demonstrated that the tendinous structures of the muscle-tendon complex are compliant enough to influence the biomechanical behavior (length change, shortening velocity, and so on) of fascicles substantially. It was discussed that the fascicles are a force generator rather than a work generator; the tendinous structures function not only as an energy re-distributor but also as a power amplifier, and the interaction between fascicles and tendinous structures is essential for generating higher joint power outputs during the late pushoff phase in human vertical jumping. This phenomenon could be explained based on the force-length/velocity relationships of each element (contractile and series elastic elements) in the muscle-tendon complex during movements. Through computer simulation using a Hill-type muscle-tendon complex model, the benefit of making a countermovement was examined in relation to the compliance of the muscle-tendon complex and the length ratio between the contractile and series elastic elements. Also, the integral roles of the series elastic element were simulated in a cyclic human heel-raise exercise. It was suggested that the storage and reutilization of elastic energy by the tendinous structures play an important role in enhancing work output and movement efficiency in many sorts of human movements.     

Human leg model predicts ankle muscle-tendon morphology, state, roles and energetics in walking

A common feature in biological neuromuscular systems is the redundancy in joint actuation. Understanding how these redundancies are resolved in typical joint movements has been a long-standing problem in biomechanics, neuroscience and prosthetics. Many empirical studies have uncovered neural, mechanical and energetic aspects of how humans resolve these degrees of freedom to actuate leg joints for common tasks like walking. However, a unifying theoretical framework that explains the many independent empirical observations and predicts individual muscle and tendon contributions to joint actuation is yet to be established. Here we develop a computational framework to address how the ankle joint actuation problem is resolved by the neuromuscular system in walking. Our framework is founded upon the proposal that a consideration of both neural control and leg muscle-tendon morphology is critical to obtain predictive, mechanistic insight into individual muscle and tendon contributions to joint actuation. We examine kinetic, kinematic and electromyographic data from healthy walking subjects to find that human leg muscle-tendon morphology and neural activations enable a metabolically optimal realization of biological ankle mechanics in walking. This optimal realization (a) corresponds to independent empirical observations of operation and performance of the soleus and gastrocnemius muscles, (b) gives rise to an efficient load-sharing amongst ankle muscle-tendon units and (c) causes soleus and gastrocnemius muscle fibers to take on distinct mechanical roles of force generation and power production at the end of stance phase in walking. The framework outlined here suggests that the dynamical interplay between leg structure and neural control may be key to the high walking economy of humans, and has implications as a means to obtain insight into empirically inaccessible features of individual muscle and tendons in biomechanical tasks.     

Robust estimation of dominant axis of rotation

A simple method is developed for robustly estimating a fixed dominant axis of rotation (AoR) of anatomical joints from surface marker data. Previous approaches which assume a model of circular marker trajectories use plane fitting to estimate the direction of the AoR. However, when there is limited joint range of motion and rotation due to a second degree of freedom, minimizing only the planar error can give poor estimates of the AoR direction. Optimizing a cost function which includes the error component within a plane, instead of only the component orthogonal to a plane, leads to improved estimates of the AoR direction for joints which exhibit additional rotational motion from a second degree of freedom. Results from synthetic data validation show the ranges of motion where the new method has lower estimation error compared to plane-fitting techniques. Estimates of the flexion-extension AoR from empirical motion capture data of the knee and index finger joints were also more anatomically plausible.     

Electronic methods for the analysis of the human post-canine dentition

The human post-canine dentitions are described and analyzed, in terms of the relative plane surface area of the occlusal surface, by electronic methods. Dental casts are duplicated and the marked duplicates are then photographed in a 1:1 relationship. The photographic negatives are digitized, i.e., the points marked on the casts are converted to X and Y coordinates which are recorded on punched cards. Computer facilities are used to convert the X and Y coordinates into mensurational data and to calculate the relative plane surface areas of the crown components and the total relative plane surface of the crown. The accuracy of this method was verified by the statistical analyses.     

Tracking the motion of hidden segments using kinematic constraints and Kalman filtering

Motion capture for biomechanical applications involves in almost all cases sensors or markers that are applied to the skin of the body segments of interest. This paper deals with the problem of estimating the movement of connected skeletal segments from 3D position data of markers attached to the skin. The use of kinematic constraints has been shown previously to reduce the error in estimated segment movement that are due to skin and muscles moving with respect to the underlying segment. A kinematic constraint reduces the number of degrees of freedom between two articulating segments. Moreover, kinematic constraints can help reveal the movement of some segments when the 3D marker data otherwise are insufficient. Important cases include the human ankle complex and the phalangeal segments of the horse, where the movement of small segments is almost completely hidden from external observation by joint capsules and ligaments. This paper discusses the use of an extended Kalman filter for tracking a system of connected segments. The system is modeled using rigid segments connected by simplified joint models. The position and orientation of the mechanism are specified by a set of generalized coordinates corresponding to the mechanism's degrees of motion. The generalized coordinates together with their first time derivatives can be used as the state vector of a state space model governing the kinematics of the mechanism. The data collected are marker trajectories from skin-mounted markers, and the state vector is related to the position of the markers through a nonlinear function. The Jacobian of this function is derived. The practical use of the method is demonstrated on a model of the distal part of the limb of the horse. Monte Carlo simulations of marker data for a two-segment system connected by a joint with three degrees of freedom indicate that the proposed method gives significant improvement over a method, which does not make use of the joint constraint, but the method requires that the model is a good approximation of the true mechanism. Applying the method to data on the movement of the four distal-most segments of the horse's limb shows good between trial consistency and small differences between measured marker positions and marker positions predicted by the model.     

An inertial and magnetic sensor based technique for joint angle measurement

This paper describes the design and evaluation of a miniature kinematic sensor based three dimensional (3D) joint angle measurement technique. The technique uses a combination of rate gyroscope, accelerometer and magnetometer sensor signals. The technique enables 3D inter-segment joint angle measurement and could be of benefit in a variety of applications which require monitoring of joint angles. The technique is not dependent on a fixed reference coordinate system and thus may be suitable for use in a dynamic system such as a moving vehicle. The technique was evaluated by applying it to joint angle measurement of the ankle joint. Experimental results show that accurate measurement of ankle joint angles is achieved by the technique during a variety of lower leg exercises including walking.     

Isokinetic plantar flexion: experimental results and model calculations

In isokinetic experiments on human subjects, conducted to determine moments that can be exerted about a joint at different angular velocities, joint rotation starts as soon as the moment increases above the resting level. This contraction history differs from the one in experiments on isolated muscle, where the force is allowed to increase to an isometric level before shortening is initiated. The purpose of the present study was to determine the influence of contraction history on plantar flexing moments found during maximal voluntary plantar flexion on an isokinetic dynamometer. In ten subjects, plantar flexing moments were measured as a function of ankle angle at different angular velocities. They were also calculated using a model of the muscle-tendon complex of the human triceps surae. The model incorporates elastic tendinous tissue in series with muscle fibers. The input of the model consists of time histories of active state (the force generating capacity of contractile elements) and shortening velocity of the muscle-tendon complex. Different time courses of active state were offered at fixed length of the muscle-tendon complex. The time course yielding a close match between the calculated rise of plantar flexing moment and the rise measured during fixed angle contractions was used to calculate moment-angle curves for isokinetic plantar flexion. The active state value reached when a peak occurred in calculated moment-angle curves was found to be lower if the angular velocity was made higher. Comparing measured and calculated results, it was concluded that moment-angular velocity diagrams determined in studies of isokinetic plantar flexion in human subjects reflect not only the influence of shortening velocity of contractile elements on the force which can be produced by plantar flexors.     

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.     

Fusion of radiostereometric analysis data into computed tomography space: application to the elbow joint

Improvement of joint prostheses is dependent upon information concerning the biomechanical properties of the joint. Radiostereometric analysis (RSA) and electromagnetic techniques have been applied in previous cadaver and in vivo studies on the elbow joint to provide valuable information concerning joint motion axes. However, such information is limited to mathematically calculated positions of the axes according to an orthogonal coordinate system and is difficult to relate to individual skeletal anatomy. The aim of this study was to evaluate the in vivo application of a new fusion method to provide three-dimensional (3D) visualization of flexion axes according to bony landmarks. In vivo RSA data of the elbow joint's flexion axes was combined with data obtained by 3D computed tomography (CT). Results were obtained from five healthy subjects after one was excluded due to an instable RSA marker. The median error between imported and transformed RSA marker coordinates and those obtained in the CT volume was 0.22 mm. Median maximal rotation error after transformation of the rigid RSA body to the CT volume was 0.003 degrees . Points of interception with a plane calculated in the RSA orthogonal coordinate system were imported into the CT volume, facilitating the 3D visualization of the flexion axes. This study demonstrates a successful fusion of RSA and CT data, without significant loss of RSA accuracy. The method could be used for relating individual motion axes to a 3D representation of relevant joint anatomy, thus providing important information for clinical applications such as the development of joint prostheses.     

Influence of the 3D inverse dynamic method on the joint forces and moments during gait

The joint forces and moments are commonly used in gait analysis. They can be computed by four different 3D inverse dynamic methods proposed in the literature, either based on vectors and Euler angles, wrenches and quaternions, homogeneous matrices, or generalized coordinates and forces. In order to analyze the influence of the inverse dynamic method, the joint forces and moments were computed during gait on nine healthy subjects. A ratio was computed between the relative dispersions (due to the method) and the absolute amplitudes of the gait curves. The influence of the inverse dynamic method was negligible at the ankle (2%) but major at the knee and the hip joints (40%). This influence seems to be due to the dynamic computation rather than the kinematic computation. Compared to the influence of the joint center location, the body segment inertial parameter estimation, and more, the influence of the inverse dynamic method is at least of equivalent importance. This point should be confirmed with other subjects, possibly pathologic, and other movements.