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.     

A Three-Dimensional Musculoskeletal Model of the Human Knee Joint. Part 1: Theoretical Construction

A three-dimensional model of the knee is developed to study the interactions between the muscles, ligaments, and bones during activity. The geometry of the distal femur, proximal tibia, and patella is based on cadaver data reported for an average-size knee. The shapes of the femoral condyles are represented by high-order polynomials: the tibial plateaux and patellar facets are approximated as flat surfaces. The contacting surfaces of the femur and tibia are modeled as deformable, while those of the femur and patella are assumed to be rigid. Interpenetration of the femur and tibia is taken into account by modeling cartilage as a thin, linear, elastic layer, mounted on rigid bone. Twelve elastic elements describe the geometry and mechanical properties of the cruciate ligaments, the collateral ligaments, and the posterior capsule. The model is actuated by thirteen musculotendinous units, each unit modeled as a three-element muscle in series with tendon. The path of each muscle is approximated as a straight line, except where it contacts and wraps around bone and other muscles; changes in muscle paths are taken into account using data obtained from MRI. In the first part of this paper, the model is used to simulate passive knee flexion. Quantitative comparisons of the model results with experimental data reported in the literature indicate that the relative movements of the bones and the geometry of the ligaments and muscles in the model are similar to those evident in the real knee. In Part II, the model is used to describe knee-ligament function during anterior-posterior draw, axial rotation, and isometric knee-extension and knee-flexion exercises.     

Comparison of an EMG-based and a stress-based method to predict shoulder muscle forces

The estimation of muscle forces in musculoskeletal shoulder models is still controversial. Two different methods are widely used to solve the indeterminacy of the system: electromyography (EMG)-based methods and stress-based methods. The goal of this work was to evaluate the influence of these two methods on the prediction of muscle forces, glenohumeral load and joint stability after total shoulder arthroplasty. An EMG-based and a stress-based method were implemented into the same musculoskeletal shoulder model. The model replicated the glenohumeral joint after total shoulder arthroplasty. It contained the scapula, the humerus, the joint prosthesis, the rotator cuff muscles supraspinatus, subscapularis and infraspinatus and the middle, anterior and posterior deltoid muscles. A movement of abduction was simulated in the plane of the scapula. The EMG-based method replicated muscular activity of experimentally measured EMG. The stress-based method minimised a cost function based on muscle stresses. We compared muscle forces, joint reaction force, articular contact pressure and translation of the humeral head. The stress-based method predicted a lower force of the rotator cuff muscles. This was partly counter-balanced by a higher force of the middle part of the deltoid muscle. As a consequence, the stress-based method predicted a lower joint load (16% reduced) and a higher superior–inferior translation of the humeral head (increased by 1.2 mm). The EMG-based method has the advantage of replicating the observed cocontraction of stabilising muscles of the rotator cuff. This method is, however, limited to available EMG measurements. The stress-based method has thus an advantage of flexibility, but may overestimate glenohumeral subluxation.     

A muscle-path-plane method for representing muscle contraction during joint movement

Traditional muscle paths (the straight-line model and the viapoint-line model) emphasise either the mechanical properties that arouse joint movement or the morphological characteristics of the muscles. To consider both the factors, a muscle-path-plane (MPP) method is introduced to model the paths of muscles during joint movement. This method is based on the premise that there is a MPP, constructed by origin, insertion and MPP control point, which represents the major direction of the muscle contraction for an arbitrary joint configuration at any time. Then, we can calculate the functions and the lengths of the muscle paths during instantaneous joint movement in MPP by mathematical approaches. Taking the triceps brachii as an example, the lengths of its paths during elbow flexion are calculated and compared with the relative studies reported in the literature. It is concluded that this method can provide an insight into the simulation of the muscle contraction.     

Subject-specific geometry affects acetabular contact pressure during gait more than subject-specific loading patterns

Abstract

Finite element modeling (FEM) can predict hip cartilage contact mechanics. This study investigated how subject-specific boundary conditions and joint geometry affect acetabular cartilage contact mechanics using a multi-scale workflow. For two healthy subjects, musculoskeletal models calculated subject-specific hip kinematics and loading, which were used as boundary conditions for FEM. Cartilage contact mechanics were predicted using either generic or subject-specific FEM and boundary conditions. A subject-specific mesh resulted in a more lateral contact. Effects of subject-specific boundary conditions varied between both subjects. Results highlight the complex interplay between loading and kinematics and their effect on cartilage contact mechanics.     

A Framework for Structured Modeling of Skeletal Muscle

The aim of this study is to present a detailed continuum mechanics formulation, and the corresponding algorithms, to predict the deformation of skeletal muscle at different structural levels, starting from the muscle fiber level. The model is used to investigate force production and structural changes during isometric and dynamic contractions of the cat medial gastrocnemius. From a comparison with experimental data obtained in our own laboratories, we conclude that the model faithfully predicts all of the observations pertaining to force production, fascicle length and angle of pennation under various test conditions.     

EMG muscle scanning: Comparison to attached surface electrodes

Electromyographic (EMG) muscle scanning measures brief samples of integrated muscle action potentials from individual muscles using a hand-held scanner with post-style electrodes. This scanning technique is widely used by biofeedback practitioners to quickly assess muscle activity in the diagnosis of musculoskeletal disorders. In an effort to compare muscle scanning with the established technique using attached surface electrodes, ten healthy subjects (25–35 years old) were scanned using 2-second sampling at five bilateral muscle sites while simultaneously monitoring the same sites with surface electrodes. This was repeated using 10-second scanning samplings. Pearson's product-moment correlations between scanning for 2 seconds and prolonged surface recording at all sites were 0.54–0.89. Scanning for 10 seconds improved the correlations to 0.68–0.91. EMG scanning for 2 seconds compares favorably with attached surface electrode recording. Comparisons are further improved by 10-second scans.     

Skeletal muscle design to meet functional demands

Skeletal muscles are length- and velocity-sensitive force producers, constructed of a vast array of sarcomeres. Muscles come in a variety of sizes and shapes to accomplish a wide variety of tasks. How does muscle design match task performance? In this review, we outline muscle''s basic properties and strategies that are used to produce movement. Several examples are provided, primarily for human muscles, in which skeletal muscle architecture and moment arms are tailored to a particular performance requirement. In addition, the concept that muscles may have a preferred sarcomere length operating range is also introduced. Taken together, the case is made that muscles can be fine-tuned to perform specific tasks that require actuators with a wide range of properties.     

Functional differentiation of trailing and leading forelimbs during locomotion on the ground and on a horizontal branch in the European red squirrel (Sciurus vulgaris, Rodentia)

Mammalian locomotion is characterized by the frequent use of in-phase gaits in which the footfalls of the left and right fore- or hindlimbs are unevenly spaced in time. Although previous studies have identified a functional differentiation between the first limb (trailing limb) and the second limb (leading limb) to touch the ground during terrestrial locomotion, the influence of a horizontal branch on limb function has never been explored. To determine the functional differences between trailing and leading forelimbs during locomotion on the ground and on a horizontal branch, X-ray motion analysis and force measurements were carried out in two European red squirrels (Sciurus vulgaris, Rodentia). The differences observed between trailing and leading forelimbs were minimal during terrestrial locomotion, where both limbs fulfill two functions and go through a shock-absorbing phase followed by a generating phase. During locomotion on a horizontal branch, European red squirrels reduce speed and all substrate reaction forces transmitted may be due to the reduction of vertical oscillation of the center of mass. Further adjustments during locomotion on a horizontal branch differ significantly between trailing and leading forelimbs and include limb flexion, lead intervals, limb protraction and vertical displacement of the scapular pivot. Consequently, trailing and leading forelimbs perform different functions. Trailing forelimbs function primarily as shock-absorbing elements, whereas leading forelimbs are characterized by a high level of stiffness. This functional differentiation indicates that European red squirrels ‘test’ the substrate for stability with the trailing forelimb, while the leading forelimb responds to or counteracts swinging or snapping branches.     

Performance criteria for generating predictive optimal control simulations of bicycle pedaling

The purpose of this study was to identify one or more performance-based criteria that may be used to generate predictive optimal control simulations of submaximal pedaling. Two-legged pedaling simulations were generated based on minimizing muscle activation, muscle stress, metabolic energy, time derivative of muscle force, and minimizing metabolic energy while pedaling smoothly. The simulations based on minimizing muscle activation and muscle stress most closely matched experimental pedaling data, with the activation criterion better matching experimental muscle activation timing. We conclude that predictive simulations of submaximal pedaling may be generated using a cost function based on minimizing muscle activation.