Identification of passive elastic joint moment-angle relationships in the lower extremity

A cadaveric study of wrist-joint moments in chimpanzees and orangutans with implications for the evolution of knuckle-walking

Understanding the mechanism underlying the evolution of knuckle-walking in African great apes but not in humans may provide important implications about the origin and evolution of human bipedal locomotion. In this study, aiming to reveal possible structural adaptations of the chimpanzee's forearm and hand musculature related to knuckle-walking, we measure the passive elastic moment of the chimpanzee's and orangutan's wrist as it was rotated into extension, immobilizing the metacarpophalangeal joint at three different positions: extended (as in knuckle-walking), flexed (as in fist-walking), and an intermediate position. Our findings demonstrate that when the metacarpophalangeal joints are extended, the rigidity of the wrist joint in the extended direction increases. This increased rigidity is attributed to the passive elongation and force generation of digital flexor muscles, which are relatively short in chimpanzees. Consequently, this enhanced wrist-joint rigidity contributes to the stability and energetically efficient transmission of propulsive force to the ground during the stance phase. Overall, our study supports the hypothesis that knuckle-walking is an adaptation to terrestrial locomotion for an ancestor characterized by the restricted capacity for wrist extension owing to the relatively shorter tendons of digital flexor muscles.

Insights into muscle metabolic energetics: Modelling muscle-tendon mechanics and metabolic rates during walking across speeds

The metabolic energy rate of individual muscles is impossible to measure without invasive procedures. Prior studies have produced models to predict metabolic rates based on experimental observations of isolated muscle contraction from various species. Such models can provide reliable predictions of metabolic rates in humans if muscle properties and control are accurately modeled. This study aimed to examine how muscle-tendon model individualization and metabolic energy models influenced estimation of muscle-tendon states and time-series metabolic rates, to evaluate the agreement with empirical data, and to provide predictions of the metabolic rate of muscle groups and gait phases across walking speeds. Three-dimensional musculoskeletal simulations with prescribed kinematics and dynamics were performed. An optimal control formulation was used to compute muscle-tendon states with four levels of individualization, ranging from a scaled generic model and muscle controls based on minimal activations, inclusion of calibrated muscle passive forces, personalization of Achilles and quadriceps tendon stiffnesses, to finally informing muscle controls with electromyography. We computed metabolic rates based on existing models. Simulations with calibrated passive forces and personalized tendon stiffness most accurately estimate muscle excitations and fiber lengths. Interestingly, the inclusion of electromyography did not improve our estimates. The whole-body average metabolic cost was better estimated with a subset of metabolic energy models. We estimated metabolic rate peaks near early stance, pre-swing, and initial swing at all walking speeds. Plantarflexors accounted for the highest cost among muscle groups at the preferred speed and were similar to the cost of hip adductors and abductors combined. Also, the swing phase accounted for slightly more than one-quarter of the total cost in a gait cycle, and its relative cost decreased with walking speed. Our prediction might inform the design of assistive devices and rehabilitation treatment. The code and experimental data are available online.

Springs vs. motors: Ideal assistance in the lower limbs during walking at different speeds

Recent years have witnessed breakthroughs in assistive exoskeletons; both passive and active devices have reduced metabolic costs near preferred walking speed by assisting muscle actions. Metabolic reductions at multiple speeds should thus also be attainable. Musculoskeletal simulation can potentially predict the interaction between assistive moments, muscle-tendon mechanics, and walking energetics. In this study, we simulated devices' optimal assistive moments based on minimal muscle activations during walking with prescribed kinematics and dynamics. We used a generic musculoskeletal model with tuned muscle-tendon parameters and computed metabolic rates from muscle actions. We then simulated walking across multiple speeds and with two ideal actuation modes-motor-based and spring-based-to assist ankle plantarflexion, knee extension, hip flexion, and hip abduction and compared computed metabolic rates. We found that both actuation modes considerably reduced physiological joint moments but did not always reduce metabolic rates. Compared to unassisted conditions, motor-based ankle plantarflexion and hip flexion assistance reduced metabolic rates, and this effect was more pronounced as walking speed increased. Spring-based hip flexion and abduction assistance increased metabolic rates at some walking speeds despite a moderate decrease in some muscle activations. Both modes of knee extension assistance reduced metabolic rates to a small extent, even though the actuation contributed with practically the entire net knee extension moment during stance. Motor-based hip abduction assistance reduced metabolic rates more than spring-based assistance, though this reduction was relatively small. Our study also suggests that an assistive strategy based on minimal muscle activations might result in a suboptimal reduction of metabolic rates. Future work should experimentally validate the effects of assistive moments and refine modeling assumptions accordingly. Our computational workflow is freely available online.

Relevance of instrumented gait analysis in the prediction of the rebound phenomenon after guided growth intervention

Predictors of rebound after correction of coronal plane deformities using temporary hemiepiphysiodesis (TH) are not well defined. The following research questions were tested: (1) Is the dynamic knee joint load useful to improve rebound prediction accuracy? (2) Does a large initial deformity play a critical role in rebound development? (3) Are BMI and a young age risk factors for rebound? Fifty children and adolescents with idiopathic knee valgus malalignment were included. A deviation of the mechanical femorotibial angle (MFA) of ≥ 3° into valgus between explantation and the one-year follow-up period was chosen to classify a rebound. A rebound was detected in 22 of the 50 patients (44%). Two predictors of rebound were identified: 1. reduced peak lateral knee joint contact force in the first half of the stance phase at the time of explantation (72.7% prediction); 2. minor initial deformity according to the MFA (70.5% prediction). The best prediction (75%) was obtained by including both parameters in the binary logistic regression method. A TH should not be advised in patients with a minor initial deformity of the leg axis. Dynamic knee joint loading using gait analysis and musculoskeletal modeling can be used to determine the optimum time to remove the plates.

Experiment-guided tuning of muscle-tendon parameters to estimate muscle fiber lengths and passive forces

The workflow to simulate motion with recorded data usually starts with selecting a generic musculoskeletal model and scaling it to represent subject-specific characteristics. Simulating muscle dynamics with muscle-tendon parameters computed from existing scaling methods in literature, however, yields some inconsistencies compared to measurable outcomes. For instance, simulating fiber lengths and muscle excitations during walking with linearly scaled parameters does not resemble established patterns in the literature. This study presents a tool that leverages reported in vivo experimental observations to tune muscle-tendon parameters and evaluates their influence in estimating muscle excitations and metabolic costs during walking. From a scaled generic musculoskeletal model, we tuned optimal fiber length, tendon slack length, and tendon stiffness to match reported fiber lengths from ultrasound imaging and muscle passive force-length relationships to match reported in vivo joint moment-angle relationships. With tuned parameters, muscle contracted more isometrically, and soleus's operating range was better estimated than with linearly scaled parameters. Also, with tuned parameters, on/off timing of nearly all muscles' excitations in the model agreed with reported electromyographic signals, and metabolic rate trajectories varied significantly throughout the gait cycle compared to linearly scaled parameters. Our tool, freely available online, can customize muscle-tendon parameters easily and be adapted to incorporate more experimental data.

Hamstrings passive and active shear modulus: Implications of conventional static stretching and warmup

PURPOSE

This study compares the acute effects of a static stretching and a warmup protocol on the active and passive shear modulus of the hamstring muscles.

METHODS

Muscle shear modulus was assessed at rest and during isometric contractions at 20 % of maximal voluntary isometric contraction (MVIC).

RESULTS

After stretching, the passive shear modulus pattern was not altered, while at 20 % MVIC the biceps femoris short head (BFsh) and semimembranosus showed a shear modulus increase and decrease, respectively, which resulted on BFsh-SM pair differences (pre: 3.8 ± 16.8 vs. post: 39.3 ± 25.1 kPa; p < 0.001; d = 1.66) which was accompanied by a decrease of 18.3 % on MVIC. Following the warmup protocol, passive shear modulus remained unchanged, while active shear modulus was decreased for the semitendinosus (pre: 65.3 ± 13.5 vs. post: 60.3 ± 12.3 kPa; p = 0.035; d = 0.4). However, this difference was within the standard error of measurement (10.54 kPa), and did not impact the force production, since it increased only 1.4 % after the warmup.

CONCLUSIONS

The results of this study suggest that the passive and active shear modulus responses of the individual hamstring muscles to static stretching are muscle-specific and that passive and active hamstring shear modulus are not changed by a standard warmup intervention.

Comparing spinal loads in individuals with unilateral transtibial amputation with and without chronic low back pain: An EMG-informed approach

Chronic low back pain (cLBP) is highly prevalent after lower limb amputation (LLA), likely due in part to biomechanical factors. Here, three-dimensional full-body kinematics and kinetics during level-ground walking, at a self-selected and three controlled speeds (1.0, 1.3, and 1.6 m/s), were collected from twenty-one persons with unilateral transtibial LLA, with (n = 9) and without cLBP (n = 12). Peak compressive, mediolateral, and anteroposterior L5-S1 spinal loads were estimated from a full-body, transtibial amputation-specific OpenSim model and compared between groups. Predicted lumbar joint torques from muscle activations were compared to inverse dynamics and predicted and measured electromyographic muscle activations were compared for model evaluation and verification. There were no group differences in compressive or anterior shear forces (p > 0.466). During intact stance, peak ipsilateral loads increased with speed to a greater extent in the cLBP group vs. no cLBP group (p=0.023), while during prosthetic stance, peak contralateral loads were larger in the no cLBP group (p=0.047) and increased to a greater extent with walking speed compared to the cLBP group (p=0.008). During intact stance, intact side external obliques had higher activations in the no cLBP group (p=0.039), and internal obliques had higher activations in the cLBP group at faster walking speeds compared to the no cLBP group. Predicted muscle activations demonstrated similar activation patterns to electromyographic-measured activations (r = 0.56-0.96), and error between inverse dynamics and simulated spinal moments was low (0.08 Nm RMS error). Persons with transtibial LLA and cLBP may adopt movement strategies during walking to reduce mediolateral shear forces at the L5-S1 joint, particularly as walking speed increases. However, future work is needed to understand the time course from pain onset to chronification and the cumulative influence of increased spinal loads over time.

American Society of Biomechanics Journal of Biomechanics Award 2022: Computer models do not accurately predict human muscle passive muscle force and fiber length: Evaluating subject-specific modeling impact on musculoskeletal model predictions

Musculoskeletal models are valuable for studying and understanding the human body in a variety of clinical applications that include surgical planning, injury prevention, and prosthetic design. Subject-specific models have proven to be more accurate and useful compared to generic models. Nevertheless, it is important to validate all models when possible. To this end, gracilis muscle-tendon parameters were directly measured intraoperatively and used to test model predictions. The aim of this study was to evaluate the benefits and limitations of systematically incorporating subject-specific variables into muscle models used to predict passive force and fiber length. The results showed that incorporating subject-specific values generally reduced errors, although significant errors still existed. Optimization of the modeling parameter "tendon slack length" was explored in two cases: minimizing fiber length error and minimizing passive force error. The results showed that using all subject-specific values yielded the most favorable outcome in both models and optimization cases. However, the trade-off between fiber length error and passive force error will depend on the specific circumstances and research objectives due to significant individual errors. Notably, individual fiber length and passive force errors were as high as 20% and 37% respectively. Finally, the modeling parameter "tendon slack length" did not correlate with any real-world anatomical length.

Simulating Muscle-Level Energetic Cost Savings When Humans Run with a Passive Assistive Device

Connecting the legs with a spring attached to the shoelaces, called an exotendon, can reduce the energetic cost of running, but how the exotendon reduces the energetic burden of individual muscles remains unknown. We generated muscle-driven simulations of seven individuals running with and without the exotendon to discern whether savings occurred during the stance phase or the swing phase, and to identify which muscles contributed to energy savings. We computed differences in muscle-level energy consumption, muscle activations, and changes in muscle-fiber velocity and force between running with and without the exotendon. The seven of nine participants who reduced energy cost when running with the exotendon reduced their measured energy expenditure rate by 0.9 W/kg (8.3%). Simulations predicted a 1.4 W/kg (12.0%) reduction in the average rate of energy expenditure and correctly identified that the exotendon reduced rates of energy expenditure for all seven individuals. Simulations showed most of the savings occurred during stance (1.5 W/kg), though the rate of energy expenditure was also reduced during swing (0.3 W/kg). The energetic savings were distributed across the quadriceps, hip flexor, hip abductor, hamstring, hip adductor, and hip extensor muscle groups, whereas no changes were observed in the plantarflexor or dorsiflexor muscles. Energetic savings were facilitated by reductions in the rate of mechanical work performed by muscles and their estimated rate of heat production. By modeling muscle-level energetics, this simulation framework accurately captured measured changes in whole-body energetics when using an assistive device. This is a useful first step towards using simulation to accelerate device design by predicting how humans will interact with assistive devices that have yet to be built.

An EMG-to-Force Processing Approach to Estimating Knee Muscle Forces during Adult, Self-Selected Speed Gait

BACKGROUND

The purpose of this study was to determine the force production during self-selected speed normal gait by muscle-tendon units that cross the knee. The force of a single knee muscle is not directly measurable without invasive methods, yet invasive techniques are not appropriate for clinical use. Thus, an EMG-to-force processing (EFP) model was developed which scaled muscle-tendon unit (MTU) force output to gait EMG.

METHODS

An EMG-to-force processing (EFP) model was developed which scaled muscle-tendon unit (MTU) force output to gait EMG. Active muscle force power was defined as the product of MTU forces (derived from EFP) and that muscle's contraction velocity. Net knee EFP moment was determined by summing individual active knee muscle moments. Net knee moments were also calculated for these study participants via inverse dynamics (kinetics plus kinematics, KIN). The inverse dynamics technique used are well accepted and the KIN net moment was used to validate or reject this model. Closeness of fit of the moment power curves for the two methods (during active muscle forces) was used to validate the model.

RESULTS

The correlation between the EFP and KIN methods was sufficiently close, suggesting validation of the model's ability to provide reasonable estimates of knee muscle forces.

CONCLUSIONS

The EMG-to-force processing approach provides reasonable estimates of active individual knee muscle forces in self-selected speed walking in neurologically intact adults.

Effects of repeated sprinting on hamstring shear modulus pattern and knee flexor neuromuscular parameters

The purpose of the present study was to examine the acute effects of a maximum repeated sprint protocol on (1) hamstring shear modulus and (2) knee flexor neuromuscular parameters such as peak torque (PT) and rate of torque development (RTD). Muscle shear modulus was assessed in 18 healthy males using shear wave elastography at rest and during 30° isometric knee flexion at 20% of maximal voluntary isometric contraction, before and after a 10 × 30 m repeated sprint protocol. There was a 9% decrease in average speed between the fastest and slowest sprint (p < 0.001; d = 2.27). A pre-post decrease was observed in PT (p = 0.004; η2p = 0.399) and in the 0-50 ms (p = 0.042; η2p = 0.222), and 50-100 ms (p = 0.028; η2p = 0.254) RTD periods. For the active shear modulus, the only significant change after the sprint task was in the biceps femoris long head (BFlh) with an increase of 10% (Pre: 26.29 ± 8.89 kPa; Post: 28.93 ± 8.31 kPa; p = 0.015; d = 0.31). The present study provides evidence that repeated sprinting leads to significant decreases in average speed, PT, early RTD (0-50 ms; 50-100 ms), and to an increase in BFlh active shear modulus without changing the shear modulus of the other hamstrings muscles.

How Connecting the Legs with a Spring Improves Human Running Economy

Connecting the legs with a spring attached to the shoelaces reduces the energy cost of running, but how the spring reduces the energy burden of individual muscles remains unknown. We generated muscle-driven simulations of seven individuals running with and without the spring to discern whether savings occurred during the stance phase or the swing phase, and to identify which muscles contributed to energy savings. We computed differences in muscle-level energy consumption, muscle activations, and changes in muscle-fiber velocity and force between running with and without the spring. Across participants, running with the spring reduced the measured rate of energy expenditure by 0.9 W/kg (8.3%). Simulations predicted a 1.4 W/kg (12.0%) reduction in the average rate of energy expenditure and correctly identified that the spring reduced rates of energy expenditure for all participants. Simulations showed most of the savings occurred during stance (1.5 W/kg), though the rate of energy expenditure was also reduced during swing (0.3 W/kg). The energetic savings were distributed across the quadriceps, hip flexor, hip abductor, hamstring, hip adductor, and hip extensor muscle groups, whereas no changes in the rate of energy expenditure were observed in the plantarflexor or dorsiflexor muscles. Energetic savings were facilitated by reductions in the rate of mechanical work performed by muscles and their estimated rate of heat production. The simulations provide insight into muscle-level changes that occur when utilizing an assistive device and the mechanisms by which a spring connecting the legs improves running economy.

Peak knee joint moments accurately predict medial and lateral knee contact forces in patients with valgus malalignment

Compressive knee joint contact force during walking is thought to be related to initiation and progression of knee osteoarthritis. However, joint loading is often evaluated with surrogate measures, like the external knee adduction moment, due to the complexity of computing joint contact forces. Statistical models have shown promising correlations between medial knee joint contact forces and knee adduction moments in particularly in individuals with knee osteoarthritis or after total knee replacements (R² = 0.44-0.60). The purpose of this study was to evaluate how accurately model-based predictions of peak medial and lateral knee joint contact forces during walking could be estimated by linear mixed-effects models including joint moments for children and adolescents with and without valgus malalignment. Peak knee joint moments were strongly correlated (R² > 0.85, p < 0.001) with both peak medial and lateral knee joint contact forces. The knee flexion and adduction moments were significant covariates in the models, strengthening the understanding of the statistical relationship between both moments and medial and lateral knee joint contact forces. In the future, these models could be used to evaluate peak knee joint contact forces from musculoskeletal simulations using peak joint moments from motion capture software, obviating the need for time-consuming musculoskeletal simulations.

Do Repeated Sprints Affect the Biceps Femoris Long Head Architecture in Football Players with and without an Injury History?-A Retrospective Study

The aim of this study was to compare the biceps femoris long head (BFlh) architecture between football players with (twelve) and without (twenty) history of BFlh injury before and after a repeated sprint task. Fascicle length (FL), pennation angle (PA) and muscle thickness (MT) were assessed at rest and in the active condition before and after the repeated sprint protocol. Athletes with previous BFlh injury showed shorter FL at rest (p = 0.014; η2p = 0.196) and active state (p < 0.001; η2p = 0.413), and greater PA at rest (p = 0.002; η2p = 0.307) and active state (p < 0.001; η2p = 0.368) before and after the task. Intra-individual comparisons showed that injured limbs have shorter FL at rest (p = 0.012; η2p = 0.519) and in the active state (p = 0.039; η2p = 0.332), and greater PA in passive (p < 0.001; η2p = 0.732) and active conditions (p = 0.018; η2p = 0.412), when compared with contralateral limbs. Injured players, at rest and in the active condition, display shorter BFlh FL and greater PA than contralateral and healthy controls after repeated sprints. Moreover, the BFlh of injured players presented a different architectural response to the protocol compared with the healthy controls.

Muscle coordination retraining inspired by musculoskeletal simulations reduces knee contact force

Humans typically coordinate their muscles to meet movement objectives like minimizing energy expenditure. In the presence of pathology, new objectives gain importance, like reducing loading in an osteoarthritic joint, but people often do not change their muscle coordination patterns to meet these new objectives. Here we use musculoskeletal simulations to identify simple changes in coordination that can be taught using electromyographic biofeedback, achieving the therapeutic goal of reducing joint loading. Our simulations predicted that changing the relative activation of two redundant ankle plantarflexor muscles—the gastrocnemius and soleus—could reduce knee contact force during walking, but it was unclear whether humans could re-coordinate redundant muscles during a complex task like walking. Our experiments showed that after a single session of walking with biofeedback of summary measures of plantarflexor muscle activation, healthy individuals reduced the ratio of gastrocnemius-to-soleus muscle activation by 25 ± 15% (p = 0.004, paired t test, n = 10). Participants who walked with this “gastrocnemius avoidance” gait pattern reduced late-stance knee contact force by 12 ± 12% (p = 0.029, paired t test, n = 8). Simulation-informed coordination retraining could be a promising treatment for knee osteoarthritis and a powerful tool for optimizing coordination for a variety of rehabilitation and performance applications.

Handle With Care: The Anterior Hip Capsule Plays a Key Role in Daily Hip Performance

Background:

Passive energy storage and return has long been recognized as one of the central mechanisms for minimizing the energy cost needed for terrestrial locomotion. Although the iliofemoral ligament (IFL) is the strongest ligament in the body, its potential role in energy-efficient walking remains unexplored.

Purpose:

To identify the contribution of the IFL to the amount of work performed by the hip muscles for normal, straight-level walking.

Study Design:

Controlled laboratory study.

Methods:

Straight-level walking of 50 healthy and injury-free adults was simulated using the AnyBody Modeling System. For each participant, the bone morphology and soft tissue properties were nonuniformly scaled. The superior and inferior parts of the IFL were represented by 2 springs each, and a linear force-strain relation was defined. A parameter study was conducted to account for the uncertainty surrounding the mechanical properties of the IFL. The work required from the gluteus, quadriceps, iliopsoas, and sartorius with and without inclusion of the IFL was calculated. Analysis of variance with subsequent post hoc paired t test was used to test the significance of IFL presence on the required mechanical work.

Results:

During walking, the strain in the IFL reached a median of 18.7% (95% CI, 8.0%-26.5%), with the largest values obtained at toe-off. With the IFL undamaged and fully operational, the effort required by the hip flexor muscles was reduced by a median of 54% (99% CI, 45%-62%) for the iliopsoas and by a median of 41% (99% CI, 27%-54%) for the sartorius muscles. The inclusion of the IFL did not significantly alter the work required by the gluteus and the quadriceps.

Conclusion:

The findings emphasized the key role the IFL plays in hip flexion by working synergistically with the hip musculature.

Clinical Relevance:

The importance of the contribution of the IFL to the hip flexors warrants careful handling and repair of these ligaments in cases of surgery and structural damage.

Coupled exoskeleton assistance simplifies control and maintains metabolic benefits: A simulation study

Assistive exoskeletons can reduce the metabolic cost of walking, and recent advances in exoskeleton device design and control have resulted in large metabolic savings. Most exoskeleton devices provide assistance at either the ankle or hip. Exoskeletons that assist multiple joints have the potential to provide greater metabolic savings, but can require many actuators and complicated controllers, making it difficult to design effective assistance. Coupled assistance, when two or more joints are assisted using one actuator or control signal, could reduce control dimensionality while retaining metabolic benefits. However, it is unknown which combinations of assisted joints are most promising and if there are negative consequences associated with coupled assistance. Since designing assistance with human experiments is expensive and time-consuming, we used musculoskeletal simulation to evaluate metabolic savings from multi-joint assistance and identify promising joint combinations. We generated 2D muscle-driven simulations of walking while simultaneously optimizing control strategies for simulated lower-limb exoskeleton assistive devices to minimize metabolic cost. Each device provided assistance either at a single joint or at multiple joints using massless, ideal actuators. To assess if control could be simplified for multi-joint exoskeletons, we simulated different control strategies in which the torque provided at each joint was either controlled independently or coupled between joints. We compared the predicted optimal torque profiles and changes in muscle and total metabolic power consumption across the single joint and multi-joint assistance strategies. We found multi-joint devices-whether independent or coupled-provided 50% greater metabolic savings than single joint devices. The coupled multi-joint devices were able to achieve most of the metabolic savings produced by independently-controlled multi-joint devices. Our results indicate that device designers could simplify multi-joint exoskeleton designs by reducing the number of torque control parameters through coupling, while still maintaining large reductions in metabolic cost.

Effect of knee joint angle on individual hamstrings morphology quantified using free-hand 3D ultrasonography

Exercise responses and injury rates differ between individual hamstrings and this may be linked with their morphology. The aim of this study was to compare muscle length and tendon dimensions between the individual hamstrings at two knee joint angles using free hand three-dimensional ultrasound (3D US). Muscle-tendon length and distal tendon cross-sectional area (CSA), volume, length and echogenicity of biceps femoris long (BFlh) and short (BFsh) head, semimembranosus (SM) and semitendinosus (ST) of 16 individuals were measured using free-hand 3D US at 0° (full extension) and 45° of knee flexion. ST showed the greatest length than all muscles and BFsh the lowest (p < 0.05). No difference was observed between SM and BFlh length (p > 0.05). Of the four muscles, ST tendon was longer, with less volume and CSA but greater echogenicity than the other tendons. In contrast, SM and BFlh showed shorter tendons and lower echogenicity but a greater volume and CSA than ST (p < 0.05). Muscle and tendon lengthened from 45° to 0° knee flexion angle (p < 0.05) but this change was not statistically different between individual hamstrings (p > 0.05). Freehand 3D US indicated that hamstring muscle length and distal tendon dimensions differ between individual hamstrings. All muscles and tendons lengthened as the knee was extended but this change was similar for all individual hamstrings.

Biomechanical responses to landing strategies of female artistic gymnasts

HIGHLIGHTS

An increase in the range of motion at the knee and hip may support the recommendation of the men's landing style.Gymnasts appear to utilise individual landing strategies to complete the landing objective, supporting the use of a single-subject design.

Predictive simulations of running gait reveal a critical dynamic role for the tail in bipedal dinosaur locomotion

Locomotion has influenced the ecology, evolution, and extinction of species throughout history, yet studying locomotion in the fossil record is challenging. Computational biomechanics can provide novel insight by mechanistically relating observed anatomy to whole-animal function and behavior. Here, we leverage optimal control methods to generate the first fully predictive, three-dimensional, muscle-driven simulations of locomotion in an extinct terrestrial vertebrate, the bipedal non-avian theropod dinosaur Coelophysis. Unexpectedly, our simulations involved pronounced lateroflexion movements of the tail. Rather than just being a static counterbalance, simulations indicate that the tail played a crucial dynamic role, with lateroflexion acting as a passive, physics-based mechanism for regulating angular momentum and improving locomotor economy, analogous to the swinging arms of humans. We infer this mechanism to have existed in many other bipedal non-avian dinosaurs as well, and our methodology provides new avenues for exploring the functional diversity of dinosaur tails in the future.

Reducing the energy cost of running using a lightweight, low-profile elastic exosuit

Background

Human beings can enhance their distance running performance with the help of assistive devices. Although several such devices are available, they are heavy and bulky, which limits their use in everyday activities. In this study, we developed a lightweight running assistive device with a low-profile design. The device applies a flexion moment to the hip according to the hip extension within a specific range of motion to assist running.

Methods

A passive exosuit was fabricated using textile materials and elastic bands. The deformation of the suit was measured and compensated for in the design. The fabricated suit was tested on eight participants (age: 24.4 ± 3.8 y; height: 1.72 ± 0.05 m; weight: 74.5 ± 6.1 kg) who were instructed to run on a treadmill at a speed of 2.5 m/s. Through indirect calorimetry, the metabolic rate was measured for the no-suit condition and three band conditions. Variations in the spatiotemporal parameters were measured using a motion capture system and force-sensing resistors (FSRs).

Results

When using the fabricated device, seven out of the eight participants exhibited a reduced metabolic rate in at least one of the three band conditions. An average reduction of − 4.7 ± 1.4% (mean ± standard error of the mean (s.e.m.), two-sided paired t-test, p = 0.017) was achieved when using the best-fitting bands compared to the average of the two no-suit conditions. No statistically significant changes were observed in the spatiotemporal parameters, except for the stance duration in the medium assistance force condition.

Conclusions

The proposed passive exosuit, which has a low weight of 609 g and small extrusion of 2.5 cm from the body in standing posture, can reduce the metabolic rate during running. The proposed device can potentially be used every day owing to its low-profile design and low weight, thereby overcoming the limitations of existing portable devices targeting the hip joints.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12984-021-00928-x.

Influence of the femoral offset on the muscles passive resistance in total hip arthroplasty

BACKGROUND

Soft tissue tension is treated as a crucial factor influencing the post-THA dislocation. The femoral offset is regarded as one of the major parameters responsible for the stabilization of the prosthesis. It is unclear which soft tissue is mostly affected by the offset changes.

METHODS

A finite element model of the hip was created. The model comprised muscles, bones, a stem, the acetabular component and a liner. The muscles were modelled as a Hill-type musculo-tendon nonlinear springs. Nonlinear analyses of the hip flexion and internal rotation were performed for the two values of the femoral stem offset.

RESULTS

We observed that the quadratus femoris and gluteus medius produce the largest resisting moment opposing the external load excreted by the surgeon during the intraoperative hip dislocation test.

CONCLUSIONS

An increased femoral offset increases the stretching of the quadratus femoris muscle significantly and provides the growth of its initial passive force. This muscle serves as a stiff band, providing stabilisation of the hip prosthesis, measured during the simulated intraoperative test.

Predictive Simulations of Musculoskeletal Function and Jumping Performance in a Generalized Bird

Jumping is a common, but demanding, behavior that many animals employ during everyday activity. In contrast to jump-specialists such as anurans and some primates, jumping biomechanics and the factors that influence performance remains little studied for generalized species that lack marked adaptations for jumping. Computational biomechanical modeling approaches offer a way of addressing this in a rigorous, mechanistic fashion. Here, optimal control theory and musculoskeletal modeling are integrated to generate predictive simulations of maximal height jumping in a small ground-dwelling bird, a tinamou. A three-dimensional musculoskeletal model with 36 actuators per leg is used, and direct collocation is employed to formulate a rapidly solvable optimal control problem involving both liftoff and landing phases. The resulting simulation raises the whole-body center of mass to over double its standing height, and key aspects of the simulated behavior qualitatively replicate empirical observations for other jumping birds. However, quantitative performance is lower, with reduced ground forces, jump heights, and muscle–tendon power. A pronounced countermovement maneuver is used during launch. The use of a countermovement is demonstrated to be critical to the achievement of greater jump heights, and this phenomenon may only need to exploit physical principles alone to be successful; amplification of muscle performance may not necessarily be a proximate reason for the use of this maneuver. Increasing muscle strength or contractile velocity above nominal values greatly improves jump performance, and interestingly has the greatest effect on more distal limb extensor muscles (i.e., those of the ankle), suggesting that the distal limb may be a critical link for jumping behavior. These results warrant a re-evaluation of previous inferences of jumping ability in some extinct species with foreshortened distal limb segments, such as dromaeosaurid dinosaurs.

Simulations prédictives de la fonction musculo-squelettique et des performances de saut chez un oiseau généralisé Sauter est un comportement commun, mais exigeant, que de nombreux animaux utilisent au cours de leurs activités quotidiennes. Contrairement aux spécialistes du saut tels que les anoures et certains primates, la biomécanique du saut et les facteurs qui influencent la performance restent peu étudiés pour les espèces généralisées qui n’ont pas d’adaptations marquées pour le saut. Les approches de modélisation biomécanique computationnelle offrent un moyen d’aborder cette question de manière rigoureuse et mécaniste. Ici, la théorie du contrôle optimal et la modélisation musculo-squelettique sont intégrées pour générer des simulations prédictives du saut en hauteur maximal chez un petit oiseau terrestre, le tinamou. Un modèle musculo-squelettique tridimensionnel avec 36 actionneurs par patte est utilisé, et une méthode numérique nommée “direct collocation” est employée pour formuler un problème de contrôle optimal rapidement résoluble impliquant les phases de décollage et d’atterrissage. La simulation qui en résulte élève le centre de masse du corps entier à plus du double de sa hauteur debout, et les aspects clés du comportement simulé reproduisent qualitativement les observations empiriques d’autres oiseaux sauteurs. Cependant, les performances quantitatives sont moindres, avec une réduction des forces au sol, des hauteurs de saut et de la puissance musculo-tendineuse. Une manœuvre de contre-mouvement prononcée est utilisée pendant le lancement. Il a été démontré que l’utilisation d’un contre-mouvement est essentielle à l’obtention de hauteurs de saut plus importantes, et il se peut que ce phénomène doive exploiter uniquement des principes physiques pour réussir; l’amplification de la performance musculaire n’est pas nécessairement une raison immédiate de l’utilisation de cette manœuvre. L’augmentation de la force musculaire ou de la vitesse de contraction au-dessus des valeurs nominales améliore grandement la performance de saut et, fait intéressant, a le plus grand effet sur les muscles extenseurs des membres plus distaux (c'est-à-dire ceux de la cheville), ce qui suggère que le membre distal peut être un lien critique pour le comportement de saut. Ces résultats justifient une réévaluation des déductions précédentes de la capacité de sauter chez certaines espèces éteintes avec des segments de membres distaux raccourcis, comme les dinosaures droméosauridés.

Voorspellende simulaties van musculoskeletale functie en springprestaties bij een gegeneraliseerde vogel Springen is een veel voorkomend, maar veeleisend, gedrag dat veel dieren toepassen tijdens hun dagelijkse bezigheden. In tegenstelling tot de springspecialisten zoals de anura en sommige primaten, is de biomechanica van het springen en de factoren die de prestaties beïnvloeden nog weinig bestudeerd voor algemene soorten die geen uitgesproken adaptaties voor het springen hebben. Computationele biomechanische modelbenaderingen bieden een manier om dit op een rigoureuze, mechanistische manier aan te pakken. Hier worden optimale controle theorie en musculoskeletale modellering geïntegreerd om voorspellende simulaties te genereren van maximale hoogtesprong bij een kleine grondbewonende vogel, een tinamou. Een driedimensionaal musculoskeletaal model met 36 actuatoren per poot wordt gebruikt, en directe collocatie wordt toegepast om een snel oplosbaar optimaal controleprobleem te formuleren dat zowel de opstijg-als de landingsfase omvat. De resulterende simulatie verhoogt het lichaamszwaartepunt tot meer dan het dubbele van de stahoogte, en belangrijke aspecten van het gesimuleerde gedrag komen kwalitatief overeen met empirische waarnemingen voor andere springende vogels. De kwantitatieve prestaties zijn echter minder, met verminderde grondkrachten, spronghoogtes en spierpeeskracht. Tijdens de lancering wordt een uitgesproken tegenbewegingsmanoeuvre gebruikt. Aangetoond is dat het gebruik van een tegenbeweging van cruciaal belang is voor het bereiken van grotere spronghoogten, en dit fenomeen hoeft alleen op fysische principes te berusten om succesvol te zijn; versterking van de spierprestaties hoeft niet noodzakelijk een proximate reden te zijn voor het gebruik van deze manoeuvre. Het verhogen van de spierkracht of van de contractiesnelheid boven de nominale waarden verbetert de sprongprestatie aanzienlijk, en heeft interessant genoeg het grootste effect op de meer distale extensoren van de ledematen (d.w.z. die van de enkel), wat suggereert dat de distale ledematen een kritieke schakel kunnen zijn voor het springgedrag. Deze resultaten rechtvaardigen een herevaluatie van eerdere conclusies over springvermogen bij sommige uitgestorven soorten met voorgekorte distale ledematen, zoals dromaeosauride dinosauriërs.

Prädiktive Simulationen der muskuloskelettalen Funktion und Sprungleistung bei einem generalisierten Vogel Springen ist ein übliches jedoch anstrengendes Verhalten, das viele Tiere bei ihren täglichen Aktivitäten einsetzen. Im Gegensatz zu Springspezialisten, wie Fröschen und einigen Primaten, sind bei allgemeinen Arten, welche keine ausgeprägten Anpassung für Sprungverhalten aufweisen, die Biomechanik beim Springen und die Faktoren, welche die Leistungsfähigkeit beeinflussen, noch wenig untersucht. Computergestützte biomechanische Modellierungsverfahren bieten hier eine Möglichkeit, dies in einer gründlichen, mechanistischen Weise anzugehen. In dieser Arbeit werden die optimale Steuerungstheorie und Muskel-Skelett-Modellierung zusammen eingesetzt, um die maximale Sprunghöhe eines kleinen bodenlebenden Vogels, eines Perlsteisshuhns, zu simulieren und zu prognostizieren. Es wird ein dreidimensionales Muskel-Skelett-Modell mit 36 Aktuatoren pro Bein verwendet, und durch direkte Kollokation wird ein schnell lösbares optimales Steuerungsproblem formuliert, das sowohl die Abstoss- als auch die Landephase umfasst. Die daraus folgende Simulation bringt den Ganzkörperschwerpunkt auf mehr als das Doppelte seiner Standhöhe und entscheidende Aspekte des simulierten Verhaltens entsprechen qualitativ empirischen Beobachtungen für andere springende Vögel. Allerdings ist die quantitative Leistungsfähigkeit geringer, mit reduzierten Bodenkräften, Sprunghöhen und Muskel-Sehnen-Kräften. Beim Abstossen wird ein ausgeprägtes Gegenbewegungsmanöver durchgeführt. Die Durchführung einer Gegenbewegung ist nachweislich entscheidend für das Erreichen grösserer Sprunghöhen, wobei dieses Phänomen möglicherweise nur physikalische Prinzipien auszuschöpfen braucht, um erfolgreich zu sein. Die Verstärkung der Muskelleistung ist daher möglicherweise nicht zwingend ein unmittelbarer Grund für die Verwendung dieses Manövers. Eine Erhöhung der Muskelkraft oder der Kontraktionsgeschwindigkeit über die Nominalwerte hinaus führt zu einer erheblichen Zunahme der Sprungleistung und hat interessanterweise den grössten Effekt bei den weiter distal gelegenen Streckmuskeln der Beine (d.h. bei denjenigen des Sprunggelenks), was darauf hindeutet, dass die distale Gliedmasse ein entscheidendes Element für das Sprungverhalten sein könnte. Diese Ergebnisse geben Anlass zur Überprüfung früherer Schlussfolgerungen hinsichtlich der Sprungfähigkeit einiger ausgestorbener Arten mit verkürzten distalen Gliedmassen, wie beispielsweise bei dromaeosauriden Dinosauriern.

Tendon release reduced joint stiffness with unaltered leg stiffness during gait in spastic diplegic cerebral palsy

Biomechanical deviations at individual joints are often identified by gait analysis of patients with cerebral palsy (CP). Analysis of the control of joint and leg stiffness of the locomotor system during gait in children with spastic diplegic CP has been used to reveal their control strategy, but the differences between before and after surgery remain unknown. The current study aimed to bridge the gap by comparing the leg stiffness—both skeletal and muscular components—and associated joint stiffness during gait in 12 healthy controls and 12 children with spastic diplegic CP before and after tendon release surgery (TRS). Each subject walked at a self-selected pace on a 10-meter walkway while their kinematic and forceplate data were measured to calculate the stiffness-related variables during loading response, mid-stance, terminal stance, and pre-swing. The CP group altered the stiffness of the lower limb joints and decreased the demand on the muscular components while maintaining an unaltered leg stiffness during stance phase after the TRS. The TRS surgery improved the joint and leg stiffness control during gait, although residual deficits and associated deviations still remained. It is suggested that the stiffness-related variables be included in future clinical gait analysis for a more complete assessment of gait in children with CP.

Verification of Finger Joint Stiffness Estimation Method With Soft Robotic Actuator

Stroke has been the leading cause of disability due to the induced spasticity in the upper extremity. The constant flexion of spastic fingers following stroke has not been well described. Accurate measurements for joint stiffness help clinicians have a better access to the level of impairment after stroke. Previously, we conducted a method for quantifying the passive finger joint stiffness based on the pressure-angle relationship between the spastic fingers and the soft-elastic composite actuator (SECA). However, it lacks a ground-truth to demonstrate the compatibility between the SECA-facilitated stiffness estimation and standard joint stiffness quantification procedure. In this study, we compare the passive metacarpophalangeal (MCP) joint stiffness measured using the SECA with the results from our designed standalone mechatronics device, which measures the passive metacarpophalangeal joint torque and angle during passive finger rotation. Results obtained from the fitting model that concludes the stiffness characteristic are further compared with the results obtained from SECA-Finger model, as well as the clinical score of Modified Ashworth Scale (MAS) for grading spasticity. These findings suggest the possibility of passive MCP joint stiffness quantification using the soft robotic actuator during the performance of different tasks in hand rehabilitation.

Musculoskeletal Model Personalization Affects Metabolic Cost Estimates for Walking

Assessment of metabolic cost as a metric for human performance has expanded across various fields within the scientific, clinical, and engineering communities. As an alternative to measuring metabolic cost experimentally, musculoskeletal models incorporating metabolic cost models have been developed. However, to utilize these models for practical applications, the accuracy of their metabolic cost predictions requires improvement. Previous studies have reported the benefits of using personalized musculoskeletal models for various applications, yet no study has evaluated how model personalization affects metabolic cost estimation. This study investigated the effect of musculoskeletal model personalization on estimates of metabolic cost of transport (CoT) during post-stroke walking using three commonly used metabolic cost models. We analyzed walking data previously collected from two male stroke survivors with right-sided hemiparesis. The three metabolic cost models were implemented within three musculoskeletal modeling approaches involving different levels of personalization. The first approach used a scaled generic OpenSim model and found muscle activations via static optimization (SOGen). The second approach used a personalized electromyographic (EMG)-driven musculoskeletal model with personalized functional axes but found muscle activations via static optimization (SOCal). The third approach used the same personalized EMG-driven model but calculated muscle activations directly from EMG data (EMGCal). For each approach, the muscle activation estimates were used to calculate each subject's CoT at different gait speeds using three metabolic cost models (Umberger et al., 2003; Bhargava et al., 2004; Umberger, 2010). The calculated CoT values were compared with published CoT data as a function of walking speed, step length asymmetry, stance time asymmetry, double support time asymmetry, and severity of motor impairment (i.e., Fugl-Meyer score). Overall, only SOCal and EMGCal with the Bhargava metabolic cost model were able to reproduce accurately published experimental trends between CoT and various clinical measures of walking asymmetry post-stroke. Tuning of the parameters in the different metabolic cost models could potentially resolve the observed CoT magnitude differences between model predictions and experimental measurements. Realistic CoT predictions may allow researchers to predict human performance, surgical outcomes, and rehabilitation outcomes reliably using computational simulations.

Differences between joint-space and musculoskeletal estimations of metabolic rate time profiles

Motion capture laboratories can measure multiple variables at high frame rates, but we can only measure the average metabolic rate of a stride using respiratory measurements. Biomechanical simulations with equations for calculating metabolic rate can estimate the time profile of metabolic rate within the stride cycle. A variety of methods and metabolic equations have been proposed, including metabolic time profile estimations based on joint parameters. It is unclear whether differences in estimations are due to differences in experimental data or due to methodological differences. This study aimed to compare two methods for estimating the time profile of metabolic rate, within a single dataset. Knowledge about the consistency of different methods could be useful for applications such as detecting which part of the gait cycle causes increased metabolic cost in patients. Here we compare estimations of metabolic rate time profiles using a musculoskeletal and a joint-space method. The musculoskeletal method was driven by kinematics and electromyography data and used muscle metabolic rate equations, whereas the joint-space method used metabolic rate equations based on joint parameters. Both estimations of changes in stride average metabolic rate correlated significantly with large changes in indirect calorimetry from walking on different grades showing that both methods accurately track changes. However, estimations of changes in stride average metabolic rate did not correlate significantly with more subtle changes in indirect calorimetry due to walking with different shoe inclinations, and both the musculoskeletal and joint-space time profile estimations did not correlate significantly with each other except in the most downward shoe inclination. Estimations of the relative cost of stance and swing matched well with previous simulations with similar methods and estimations from experimental perturbations. Rich experimental datasets could further advance time profile estimations. This knowledge could be useful to develop therapies and assistive devices that target the least metabolically economic part of the gait cycle.

Leg and lower limb dynamic joint stiffness during different walking speeds in healthy adults

BACKGROUND

The differences and relationship between joint stiffness and leg stiffness can be used to characterize the lower limb behavior during different walking speeds.

RESEARCH QUESTION

This study aimed to investigate the differences in whole leg and lower limb joint stiffness at different walking speeds and the interactions between leg and lower limb joint stiffness.

METHODS

Twenty-seven healthy adults, seventeen males (age: 19.6 ± 2.2 years, height: 176.0 ± 6.0 cm, mass: 69.7 ± 8.9 kg), and ten females (age: 19.1 ± 1.9 years, height: 164.0 ± 3.0 cm, mass: 59.6 ± 3.8 kg), were recruited. Dynamic leg and joint stiffness were calculated during eccentric loading from data recorded using 3D infrared motion analysis and force plates at slow, normal, and fast walking speeds. Differences in dynamic stiffness, joint angles and moments were explored between the walking speeds using Repeated Measures ANOVA with Sidak post-hoc tests. Correlations between leg, joint stiffness, and walking speed were also explored.

RESULTS

The results indicated that the leg dynamic stiffness is decreased by walking speed, however, hip and ankle joint stiffness were increased (p < 0.001) and knee stiffness was unaffected. Leg stiffness showed no correlation with hip, knee, or ankle stiffness. A positive significant correlation was seen between hip and ankle stiffness (p < 0.01) and between knee and ankle stiffness (p < 0.001), however, no correlation was seen between hip and knee stiffness.

SIGNIFICANCE

These results suggest leg stiffness is not associated with lower limb joint stiffness during eccentric loading. This provides new information on the responses of ankle, knee and hip joint stiffness to walking speed.

Leg Joint Stiffness Affects Dynamics of Backward Falling From Standing Height: A Simulation Work

Falling backward can lead to injuries including hip fracture, back injury, and traumatic brain impact among older adults. A loss of consciousness is associated with falling backward and accounts for about 13% of all falls among older adults. Little is known about the dynamics of backward falls, such as the falling duration, the impact severity, and how the fall dynamics are affected by the biomechanical properties of the lower limb joints, particularly the rotational stiffness. The purpose of this study was to investigate the influence of the stiffness of individual leg joints on the dynamics of backward falls after losing consciousness in terms of the falling duration and impact velocities. Based on a 15-segment human model, we simulated the process of falling backward by sweeping the parameter space of ankle, knee, and hip's stiffnesses varying from 0 to 8.73 N·m·deg-1 (or 500 N·m·rad-1). The results revealed that the falling duration and impact speeds of the head and hip ranged from 0.27 to 0.63 s, 2.65 to 7.88 m·s-1, and 0.35 to 3.36 m·s-1, respectively, when the stiffness of the leg joints changed within their limits. Overall, the influence of the joint stiffness on the falling dynamics (falling duration and impact speed) is comparable between hip and knee joints, whereas ankle stiffness showed little influence on the backward falling dynamics. Our findings could provide references for designing protective devices to prevent impact-induced injuries after a backward fall.

Effect of viscoelastic properties on passive torque variations at different velocities of the knee joint extension and flexion movements

This study aimed to investigate the rate of passive torque variations of human knee joint in the different velocities of knee flexion and extension movements. Ten healthy men were invited to participate in the tests. All passive torque tests were performed for the knee joint extension and flexion on the sagittal plane in three different angular velocities of 15, 45, and 120°/s; in 5 consecutive cycles; and within 0° to 100° range of motion. The electrical activity of knee joint extensor and flexor muscles was recorded until there was no muscle activity signal. A Three-element Solid Model (SLS) was used to obtain the viscose and elastic coefficients. As the velocity increases, the stretch rate in velocity-independent tissues increases, and the stretch rate in velocity-dependent tissues decreases. By increasing the velocity, the resistance of velocity-dependent parts increases, and the velocity-independent parts are not affected by velocity. Since the first torque that resists the joint movement is passive torque, the elastic and viscous torques should be simultaneously used. It is better to perform the movement at a low velocity so that less energy is lost. The viscoelastic resistance of tissues diminishes. Graphical abstract.

A novel method for prediction of postoperative global sagittal alignment based on full-body musculoskeletal modeling and posture optimization

Associations between spinal sagittal balance and pain and disability are well documented. Reciprocal changes after spinal surgery might be critical for the outcomes, but assessing their extent remains a challenge. This paper proposes a method for predicting full-body sagittal alignment including reciprocal changes in response to spinal fusion, based on musculoskeletal modeling and inverse-inverse dynamics approach. An established body model (AnyBody) was used, with fused segments modeled as rigid. Posture was optimized based on muscle expenditure minimization, following the concept of the cone of economy. The data of adult spinal fusion patients were obtained retrospectively from an ongoing clinical study. Patient spino-pelvic alignment, body weight and height, age- and pathology-related muscle deterioration, and underwent treatment details were represented in the model. Predicted postural changes were compared to follow-up radiographs to evaluate method validity. Twenty-one cases were analyzed in this preliminary study (age range = 48-74; number of fused segments 1-14). The model predictions correlated well with the radiographic measures at follow-up: TPA, r = 0.83; ΔPILL, r = 0.90; LL, r = 0.90; TK, r = 0.77. The model demonstrated high accuracy in predicting sagittal imbalance (positive predictive value = 1.00, negative predictive value = 0.75). The presented method for patient- and treatment-specific postoperative posture prediction can be used to guide preoperative planning of spinal fusion, but more extensive validation is needed.

Computational Design of FastFES Treatment to Improve Propulsive Force Symmetry During Post-stroke Gait: A Feasibility Study

Stroke is a leading cause of long-term disability worldwide and often impairs walking ability. To improve recovery of walking function post-stroke, researchers have investigated the use of treatments such as fast functional electrical stimulation (FastFES). During FastFES treatments, individuals post-stroke walk on a treadmill at their fastest comfortable speed while electrical stimulation is delivered to two muscles of the paretic ankle, ideally to improve paretic leg propulsion and toe clearance. However, muscle selection and stimulation timing are currently standardized based on clinical intuition and a one-size-fits-all approach, which may explain in part why some patients respond to FastFES training while others do not. This study explores how personalized neuromusculoskeletal models could potentially be used to enable individual-specific selection of target muscles and stimulation timing to address unique functional limitations of individual patients post-stroke. Treadmill gait data, including EMG, surface marker positions, and ground reactions, were collected from an individual post-stroke who was a non-responder to FastFES treatment. The patient's gait data were used to personalize key aspects of a full-body neuromusculoskeletal walking model, including lower-body joint functional axes, lower-body muscle force generating properties, deformable foot-ground contact properties, and paretic and non-paretic leg neural control properties. The personalized model was utilized within a direct collocation optimal control framework to reproduce the patient's unstimulated treadmill gait data (verification problem) and to generate three stimulated walking predictions that sought to minimize inter-limb propulsive force asymmetry (prediction problems). The three predictions used: (1) Standard muscle selection (gastrocnemius and tibialis anterior) with standard stimulation timing, (2) Standard muscle selection with optimized stimulation timing, and (3) Optimized muscle selection (soleus and semimembranosus) with optimized stimulation timing. Relative to unstimulated walking, the optimal control problems predicted a 41% reduction in propulsive force asymmetry for scenario (1), a 45% reduction for scenario (2), and a 64% reduction for scenario (3), suggesting that non-standard muscle selection may be superior for this patient. Despite these predicted improvements, kinematic symmetry was not noticeably improved for any of the walking predictions. These results suggest that personalized neuromusculoskeletal models may be able to predict personalized FastFES training prescriptions that could improve propulsive force symmetry, though inclusion of kinematic requirements would be necessary to improve kinematic symmetry as well.

An application of scissored-pair control moment gyroscopes in a design of wearable balance assistance device for the elderly

Impaired balance control ability and degraded functional mobility increases the risk of falling in elderly people. The elderly show more postural sway when standing compared with young people. A sway fall occurs when the center of gravity moves outside the limit of stability. In order to reduce the fall risk from the excessive sway, this study presents the design of wearable balance assistance device for the elderly. Scissored-pair control moment gyroscopes were selected as a torque actuator. A two-axis inclination sensor was used to detect the inclined angle of the wearer's body. The direction of sway was calculated from the detected inclined angle. The designed device weighs 8.2 kg with a height of 32 cm × width of 40 cm × depth of 22 cm. A multi-segment model of a standing human was used to investigate the device's performance for balance recovery. According to the simulations, balance recovery in any direction was successfully accomplished with the appropriate initial angle. The relationship between the effective initial angle and detected inclined angle was subsequently established. The stability provided by activation of the device was able to limit the unstable user's sway boundary. The designed device shows promise for use as a balance assistance device for the elderly.

Intra- and Inter-Muscular Variations in Hamstring Architecture and Mechanics and Their Implications for Injury: A Narrative Review

Understanding the architecture, anatomy, and biomechanics of the hamstrings may assist in explaining the mechanisms that affect and improve their function. The aim of this review is to specifically examine intra- and inter-muscular variations in architecture and mechanical properties of the hamstrings. Of the hamstrings, the long head of the biceps femoris shows the shortest and more pennated fibers. The semimembranosus has a similar muscle architecture with a long head of the biceps femoris but it has a different proximal attachment as well as a different moment arm compared with the long head of the biceps femoris. For the same joint motion, the semitendinosus displays less relative strain than the other hamstrings probably owing to a greater length, longer fascicles and, possibly, a longer tendon. Intra-muscular variations in architecture are documented but their implications are currently unclear. Proximally, the long head of the biceps femoris has shorter and more pennated fibers coupled with a narrower aponeurosis than distally, while the semitendinosus is the only muscle that entails a tendinous inscription. In conclusion, some of the identified intra- and inter-muscular variations in architecture may help explain why some muscles sustain injuries more than others. In the same line, exercises designed for the hamstrings may not provide the same stimulus for all components of this muscle group. Future research could examine whether intervention strategies that target specific muscles or specific areas of the hamstrings may offer additional benefits for injury prevention or rehabilitation of their function.

On the biological mechanics and energetics of the hip joint muscle-tendon system assisted by passive hip exoskeleton

Passive exoskeletons have potential advantages in reducing metabolic energy cost. We consider a passive elastic exoskeleton (peEXO) providing hip flexion moment to assist hip flexors during walking, our goal is to use a biomechanical model to explore the biological mechanics and energetics of the hip joint muscle-tendon-exotendon system for obtaining the optimum stiffness of this peEXO at the muscle-level. Based on our developed hip musculoskeletal model capable of replicating human-like behaviors, the hip peEXO is firstly abstracted as a spring (i.e. exotendon), we then simulate the peEXO assisted human walking over a series of stiffnesses, the biological muscle-tendon dynamics is optimally solved by minimizing the total metabolic cost of muscles. The simulation results are consistent with the experimental data of walking with an exoskeleton prototype. We find peEXO of minor stiffness helps reducing the muscle force, activation, and metabolic energy cost of hip flexors, especially the iliopsoas; while stiffer peEXO causes extra metabolic energy cost of antagonist muscles especially the gluteus maximus. With an optimum rotational stiffness of 350 Nm rad^-1, the peEXO can reduce the metabolic energy cost of walking by ~7.1% and the hip joint muscles simultaneously have a 3% muscle efficiency promotion. The changes in muscle-tendon dynamics indicate it is more economical to assist the hip joint compared with the ankle joint, and periods of high muscle activation are ideal assistance phases for hip joint muscles. The modeling framework provides deep insights into the potential muscle-level mechanisms which are difficult to study via experiments alone, which is helpful to recover the mechanism of how energy cost is reduced under the external passive assistance and guide the design of passive hip EXOs.

Passive elastic contribution of hip extensors to joint moments during walking in people with low back pain

BACKGROUND

It has been found that alterations in passive muscle properties may be associated with low back pain, and these may be responsible for the altered gait parameters often observed in subjects with back pain. The purpose of the present study was to assess total hip and passive hip extensor moments in people with or without low back pain during the hip flexion component of walking.

METHODS

52 subjects volunteered for this study (low back pain group, n = 25 (male n = 13, female n = 12), control group, n = 27 (male n = 15, female n = 12)). Passive hip moments were calculated using an adapted force transducer during supine testing. A biomechanical model and predictive equation were used to calculate passive hip moments during walking. Total hip moments were calculated with the use of a 9 camera, 3-D motion-capture system.

FINDINGS

Independent samples t-tests demonstrated no significant differences between groups for gait parameters or hip or knee angles. Results of the ANOVAs demonstrated significant differences in passive hip flexor moments during the second half of hip flexion (P < 0.05).There were also significant differences in hip power and work done during peaks of power absorption and the second peak of power generation (P < 0.05).

INTERPRETATION

The present data demonstrates that subjects with low back pain have altered passive hip extensor and total power and work done during walking compared with healthy controls. Biomechanical models should include individual measurements of passive joint moments.

Developing reliable measures of the passive torque-angle relationship for shoulder internal and external rotation: Implications for overhead athletics

OBJECTIVES

1) Thoroughly assess shoulder flexibility by establishing the passive torque-angle relationship for internal and external rotation with the arm in an overhead athletics position (abducted 90°) and 2) test the reliability of four passive torque-angle measures.

DESIGN

Reliability study.

SETTING

Data were collected in a university biomechanics laboratory.

PARTICIPANTS

Bilateral shoulder flexibility of 15 male college students (20.7 ± 1.1 y) was evaluated twice in two sessions over 7-10 days.

MAIN OUTCOME MEASURES

For both ER and IR, reliability was assessed bilaterally (intra-session, inter-session, and inter-tester) for the traditional range of motion measure and three novel kinetic measures: torque at end ROM, resistance onset angle, rotational stiffness. This resulted in 48 total assessments.

RESULTS

Thirty-four assessments had good to excellent reliability (ICC ≥ 0.8), 10 had fair reliability (0.7 ≤ ICC < 0.8), and 4 had poor reliability (ICC< 0.7). Three of the four flexibility measures had a good overall ICC score: ROM (0.83), torque at end ROM (0.84), and resistance onset angle (0.81). The fourth, stiffness, had a fair overall reliability score (0.74).

CONCLUSIONS

The passive torque-angle measures should be assimilated into clinical and research settings to determine the relevance to injury, rehabilitation, and performance.

Contribution of passive actions to the lower limb joint moments and powers during gait: A comparison of models

The lower limb passive actions representing the actions of all the passive periarticular structures have been shown to have a significant contribution to the power generation and absorption during gait. However, the respective magnitude of its different components was not established, although models of ligament moment were implemented in some musculoskeletal models. These ligament moments have shown to have an influence on the musculo-tendon forces and contact forces but the models used were never specifically evaluated, that is, compared to the passive and net joint moments. Two models of passive joint moments and three models of ligament moments were selected from the literature. Ten subjects (23–29 years old, 79.8 ± 9.5 kg, 1.85 ± 0.06 m) participated in the study. Each subject performed three gait cycles in a gait laboratory to acquire the kinematics and ground reaction forces and to compute the ligament, passive and net moments of the right lower limb joints. The contributions of the passive joint moments to the net joint moments were in accordance with the literature, although time shifts appeared for peaks in the hip and knee powers. Two of the models of ligament moments seemed, in fact, to represent the passive joint moments as their contributions were very similar while the third model of ligament moments seemed to represent only penalty-based joint limits. As a conclusion, this study showed that the models of ligament moments existing in the literature do not seem reliable. This study also demonstrated that the use of non-subject-specific models of the passive joint moments could be a valid approach for healthy subjects.

Inferring muscle functional roles of the ostrich pelvic limb during walking and running using computer optimization

Owing to their cursorial background, ostriches (Struthio camelus) walk and run with high metabolic economy, can reach very fast running speeds and quickly execute cutting manoeuvres. These capabilities are believed to be a result of their ability to coordinate muscles to take advantage of specialized passive limb structures. This study aimed to infer the functional roles of ostrich pelvic limb muscles during gait. Existing gait data were combined with a newly developed musculoskeletal model to generate simulations of ostrich walking and running that predict muscle excitations, force and mechanical work. Consistent with previous avian electromyography studies, predicted excitation patterns showed that individual muscles tended to be excited primarily during only stance or swing. Work and force estimates show that ostrich gaits are partially hip-driven with the bi-articular hip–knee muscles driving stance mechanics. Conversely, the knee extensors acted as brakes, absorbing energy. The digital extensors generated large amounts of both negative and positive mechanical work, with increased magnitudes during running, providing further evidence that ostriches make extensive use of tendinous elastic energy storage to improve economy. The simulations also highlight the need to carefully consider non-muscular soft tissues that may play a role in ostrich gait.

Reliability of Sonographic Assessment of Biceps Femoris Distal Tendon Strain during Passive Stretching

The purpose of this study was to determine the intra-rater, inter-examiner and inter-observer reliability of biceps femoris long head (BFlh) tendon strain using ultrasound imaging. Nineteen patients (age: 20.4 ± 0.35 y) were tested twice with a 1-wk interval. Each session included passive stretching from three different hip positions. Tests were performed independently by two examiners while BFlh tendon displacement (mm) and strain (%) were manually extracted from ultrasound video footages by two observers. Intra-rater comparisons revealed an intra-class correlation coefficient (ICC_2,1) range of 0.87 to 0.98 and a variability less than 4.74%. Inter-examiner comparisons revealed an ICC_2,1 range of 0.83 to 0.99 and less than 4.69% variability. Inter-observer ICCs ranged from 0.93 to 0.97 with variability less than 4.89%. Using a well-defined scanning protocol, two experienced examiners attained high levels of intra-rater agreement, with similarly excellent results for inter-rater and inter-observer reliability for BFlh tendon displacement and strain.

Lower extremity EMG-driven modeling of walking with automated adjustment of musculoskeletal geometry

Neuromusculoskeletal disorders affecting walking ability are often difficult to manage, in part due to limited understanding of how a patient’s lower extremity muscle excitations contribute to the patient’s lower extremity joint moments. To assist in the study of these disorders, researchers have developed electromyography (EMG) driven neuromusculoskeletal models utilizing scaled generic musculoskeletal geometry. While these models can predict individual muscle contributions to lower extremity joint moments during walking, the accuracy of the predictions can be hindered by errors in the scaled geometry. This study presents a novel EMG-driven modeling method that automatically adjusts surrogate representations of the patient’s musculoskeletal geometry to improve prediction of lower extremity joint moments during walking. In addition to commonly adjusted neuromusculoskeletal model parameters, the proposed method adjusts model parameters defining muscle-tendon lengths, velocities, and moment arms. We evaluated our EMG-driven modeling method using data collected from a high-functioning hemiparetic subject walking on an instrumented treadmill at speeds ranging from 0.4 to 0.8 m/s. EMG-driven model parameter values were calibrated to match inverse dynamic moments for five degrees of freedom in each leg while keeping musculoskeletal geometry close to that of an initial scaled musculoskeletal model. We found that our EMG-driven modeling method incorporating automated adjustment of musculoskeletal geometry predicted net joint moments during walking more accurately than did the same method without geometric adjustments. Geometric adjustments improved moment prediction errors by 25% on average and up to 52%, with the largest improvements occurring at the hip. Predicted adjustments to musculoskeletal geometry were comparable to errors reported in the literature between scaled generic geometric models and measurements made from imaging data. Our results demonstrate that with appropriate experimental data, joint moment predictions for walking generated by an EMG-driven model can be improved significantly when automated adjustment of musculoskeletal geometry is included in the model calibration process.

Musculoskeletal modelling of human ankle complex: Estimation of ankle joint moments

BACKGROUND

A musculoskeletal model for the ankle complex is vital in order to enhance the understanding of neuro-mechanical control of ankle motions, diagnose ankle disorders and assess subsequent treatments. Motions at the human ankle and foot, however, are complex due to simultaneous movements at the two joints namely, the ankle joint and the subtalar joint. The musculoskeletal elements at the ankle complex, such as ligaments, muscles and tendons, have intricate arrangements and exhibit transient and nonlinear behaviour.

METHODS

This paper develops a musculoskeletal model of the ankle complex considering the biaxial ankle structure. The model provides estimates of overall mechanical characteristics (motion and moments) of ankle complex through consideration of forces applied along ligaments and muscle-tendon units. The dynamics of the ankle complex and its surrounding ligaments and muscle-tendon units is modelled and formulated into a state space model to facilitate simulations. A graphical user interface is also developed during this research in order to include the visual anatomical information by converting it to quantitative information on coordinates.

FINDINGS

Validation of the ankle model was carried out by comparing its outputs with those published in literature as well as with experimental data obtained from an existing parallel ankle rehabilitation robot.

INTERPRETATION

Qualitative agreement was observed between the model and measured data for both, the passive and active ankle motions during trials in terms of displacements and moments.

Dynamic Flexion Stiffness of Foot Joints During Walking

BACKGROUND

Dynamic stiffness can be used for studying foot pathologic abnormalities and for developing prostheses and orthoses. Although previous works have studied the role of ankle joint stiffness during gait, other foot joints have not yet been analyzed. We sought to characterize the dynamic stiffness of the ankle, midtarsal, and metatarsophalangeal joints during normal walking.

METHODS

Kinematics and contact data from four healthy individuals during walking were registered with a three-dimensional motion analysis system and a pressure platform. Stance phases with flexion moment-angle linear relationships were identified, and dynamic stiffnesses were calculated from the slope of their linear regressions. Intraparticipant repeatability was analyzed using analyses of variance, and interparticipant variability was checked through the SD of averaged participant stiffnesses.

RESULTS

Flexion moment-angle linear relationships were identified (R(2) > 0.98) during the early and late midstance phases and the propulsion phase at the ankle (2.76, 5.23, and 3.42 N·m/kg/rad, respectively) and midtarsal (15.88, 3.90, and 4.64 N·m/kg/rad, respectively) joints. At the metatarsophalangeal joint, a linear relationship (R(2) > 0.96) occurred only during the propulsion phase (0.11 N·m/kg/rad). High dynamic stiffness variability was observed during the late and early midstance phases at the ankle and midtarsal joints, respectively.

CONCLUSIONS

These results may serve as a basis for future studies aimed at investigating the role of dynamic stiffness identified herein in different foot disorders. The importance of properly controlling the samples in such studies is highlighted. Study of the dynamic stiffnesses identified might be used in the design of prostheses, orthoses, and other assistive devices.

Instantaneous Metabolic Cost of Walking: Joint-Space Dynamic Model with Subject-Specific Heat Rate

A subject-specific model of instantaneous cost of transport (ICOT) is introduced from the joint-space formulation of metabolic energy expenditure using the laws of thermodynamics and the principles of multibody system dynamics. Work and heat are formulated in generalized coordinates as functions of joint kinematic and dynamic variables. Generalized heat rates mapped from muscle energetics are estimated from experimental walking metabolic data for the whole body, including upper-body and bilateral data synchronization. Identified subject-specific energetic parameters-mass, height, (estimated) maximum oxygen uptake, and (estimated) maximum joint torques-are incorporated into the heat rate, as opposed to the traditional in vitro and subject-invariant muscle parameters. The total model metabolic energy expenditure values are within 5.7 ± 4.6% error of the measured values with strong (R2 > 0.90) inter- and intra-subject correlations. The model reliably predicts the characteristic convexity and magnitudes (0.326-0.348) of the experimental total COT (0.311-0.358) across different subjects and speeds. The ICOT as a function of time provides insights into gait energetic causes and effects (e.g., normalized comparison and sensitivity with respect to walking speed) and phase-specific COT, which are unavailable from conventional metabolic measurements or muscle models. Using the joint-space variables from commonly measured or simulated data, the models enable real-time and phase-specific evaluations of transient or non-periodic general tasks that use a range of (aerobic) energy pathway similar to that of steady-state walking.

Muscle Synergies Facilitate Computational Prediction of Subject-Specific Walking Motions

Researchers have explored a variety of neurorehabilitation approaches to restore normal walking function following a stroke. However, there is currently no objective means for prescribing and implementing treatments that are likely to maximize recovery of walking function for any particular patient. As a first step toward optimizing neurorehabilitation effectiveness, this study develops and evaluates a patient-specific synergy-controlled neuromusculoskeletal simulation framework that can predict walking motions for an individual post-stroke. The main question we addressed was whether driving a subject-specific neuromusculoskeletal model with muscle synergy controls (5 per leg) facilitates generation of accurate walking predictions compared to a model driven by muscle activation controls (35 per leg) or joint torque controls (5 per leg). To explore this question, we developed a subject-specific neuromusculoskeletal model of a single high-functioning hemiparetic subject using instrumented treadmill walking data collected at the subject's self-selected speed of 0.5 m/s. The model included subject-specific representations of lower-body kinematic structure, foot-ground contact behavior, electromyography-driven muscle force generation, and neural control limitations and remaining capabilities. Using direct collocation optimal control and the subject-specific model, we evaluated the ability of the three control approaches to predict the subject's walking kinematics and kinetics at two speeds (0.5 and 0.8 m/s) for which experimental data were available from the subject. We also evaluated whether synergy controls could predict a physically realistic gait period at one speed (1.1 m/s) for which no experimental data were available. All three control approaches predicted the subject's walking kinematics and kinetics (including ground reaction forces) well for the model calibration speed of 0.5 m/s. However, only activation and synergy controls could predict the subject's walking kinematics and kinetics well for the faster non-calibration speed of 0.8 m/s, with synergy controls predicting the new gait period the most accurately. When used to predict how the subject would walk at 1.1 m/s, synergy controls predicted a gait period close to that estimated from the linear relationship between gait speed and stride length. These findings suggest that our neuromusculoskeletal simulation framework may be able to bridge the gap between patient-specific muscle synergy information and resulting functional capabilities and limitations.

Muscle size explains low passive skeletal muscle force in heart failure patients

Background

Alterations in skeletal muscle function and architecture have been linked to the compromised exercise capacity characterizing chronic heart failure (CHF). However, how passive skeletal muscle force is affected in CHF is not clear. Understanding passive force characteristics in CHF can help further elucidate the extent to which altered contractile properties and/or architecture might affect muscle and locomotor function. Therefore, the aim of this study was to investigate passive force in a single muscle for which non-invasive measures of muscle size and estimates of fiber force are possible, the soleus (SOL), both in CHF patients and age- and physical activity-matched control participants.

Methods

Passive SOL muscle force and size were obtained by means of a novel approach combining experimental data (dynamometry, electromyography, ultrasound imaging) with a musculoskeletal model.

Results

We found reduced passive SOL forces (∼30%) (at the same relative levels of muscle stretch) in CHF vs. healthy individuals. This difference was eliminated when force was normalized by physiological cross sectional area, indicating that reduced force output may be most strongly associated with muscle size. Nevertheless, passive force was significantly higher in CHF at a given absolute muscle length (non length-normalized) and likely explained by the shorter muscle slack lengths and optimal muscle lengths measured in CHF compared to the control participants. This later factor may lead to altered performance of the SOL in functional tasks such gait.

Discussion

These findings suggest introducing exercise rehabilitation targeting muscle hypertrophy and, specifically for the calf muscles, exercise that promotes muscle lengthening.