Determination of muscle orientations and moment arms

The moment arms and lines of action of subscapularis after the Latarjet procedure

The Latarjet procedure is an established surgical treatment for recurrent glenohumeral joint instability with glenoid bone loss. Intraoperatively, the conjoint tendon and its attachement on the coracoid bone graft is routed through a split in subscapularis where the graft is fixed to and augments the anteroinferior glenoid. The objective of this in vitro study was to quantify the influence of glenohumeral joint position and conjoint tendon force on the lines of action and moment arms of subscapularis muscle sub-regions after Latarjet surgery. Eight fresh-frozen, entire upper extremities were mounted onto a testing apparatus, and a cable-pulley system was used to apply physiological muscle loading to the major shoulder muscles. The lines of action and moment arms of four subregions of subscapularis (superior, mid-superior, mid-inferior, and inferior) were quantified radiographically with the conjoint tendon unloaded and loaded while the shoulder was in (i) 0° abduction (ii) 90° abduction (iii) 90° abduction and full external rotation (ABER), and (iv) the apprehension position, defined as ABER with 30° horizontal extension. Conjoint tendon loading after Latarjet surgery significantly increased the inferior inclination of the lines of action of the mid-inferior and inferior subregions of subscapularis in the scapular plane in ABER and apprehension positions (p < 0.001), as well as decreased the horizontal flexion moment arm of the inferior subscapularis (p = 0.040). Increased subscapularis inferior inclination may ultimately increase inferior joint shear potential, while smaller horizontal flexion leverage may reduce joint flexion capacity. The findings have implications for Latarjet surgical planning and postoperative rehabilitation prescription.

Effect of the temporal coordination and volume of cyclic mechanical loading on human Achilles tendon adaptation in men

Human tendons adapt to mechanical loading, yet there is little information on the effect of the temporal coordination of loading and recovery or the dose-response relationship. For this reason, we assigned adult men to either a control or intervention group. In the intervention group, the two legs were randomly assigned to one of five high-intensity Achilles tendon (AT) loading protocols (i.e., 90% maximum voluntary contraction and approximately 4.5 to 6.5% tendon strain) that were systematically modified in terms of loading frequency (i.e., sessions per week) and overall loading volume (i.e., total time under loading). Before, at mid-term (8 weeks) and after completion of the 16 weeks intervention, AT mechanical properties were determined using a combination of inverse dynamics and ultrasonography. The cross-sectional area (CSA) and length of the free AT were measured using magnetic resonance imaging pre- and post-intervention. The data analysis with a linear mixed model showed significant increases in muscle strength, rest length-normalized AT stiffness, and CSA of the free AT in the intervention group (p < 0.05), yet with no marked differences between protocols. No systematic effects were found considering the temporal coordination of loading and overall loading volume. In all protocols, the major changes in normalized AT stiffness occurred within the first 8 weeks and were mostly due to material rather than morphological changes. Our findings suggest that-in the range of 2.5-5 sessions per week and 180-300 s total high strain loading-the temporal coordination of loading and recovery and overall loading volume is rather secondary for tendon adaptation.

Comparison of Synergy Extrapolation and Static Optimization for Estimating Multiple Unmeasured Muscle Activations during Walking

Background

Calibrated electromyography (EMG)-driven musculoskeletal models can provide great insight into internal quantities (e.g., muscle forces) that are difficult or impossible to measure experimentally. However, the need for EMG data from all involved muscles presents a significant barrier to the widespread application of EMG-driven modeling methods. Synergy extrapolation (SynX) is a computational method that can estimate a single missing EMG signal with reasonable accuracy during the EMG-driven model calibration process, yet its performance in estimating a larger number of missing EMG signals remains unclear.

Methods

This study assessed the accuracy with which SynX can use eight measured EMG signals to estimate muscle activations and forces associated with eight missing EMG signals in the same leg during walking while simultaneously performing EMG-driven model calibration. Experimental gait data collected from two individuals post-stroke, including 16 channels of EMG data per leg, were used to calibrate an EMG-driven musculoskeletal model, providing "gold standard" muscle activations and forces for evaluation purposes. SynX was then used to predict the muscle activations and forces associated with the eight missing EMG signals while simultaneously calibrating EMG-driven model parameter values. Due to its widespread use, static optimization (SO) was also utilized to estimate the same muscle activations and forces. Estimation accuracy for SynX and SO was evaluated using root mean square errors (RMSE) to quantify amplitude errors and correlation coefficient r values to quantify shape similarity, each calculated with respect to "gold standard" muscle activations and forces.

Results

On average, SynX produced significantly more accurate amplitude and shape estimates for unmeasured muscle activations (RMSE 0.08 vs. 0.15 , r value 0.55 vs. 0.12) and forces (RMSE 101.3 N vs. 174.4   N , r value 0.53 vs. 0.07) compared to SO. SynX yielded calibrated Hill-type muscle-tendon model parameter values for all muscles and activation dynamics model parameter values for measured muscles that were similar to "gold standard" calibrated model parameter values.

Conclusions

These findings suggest that SynX could make it possible to calibrate EMG-driven musculoskeletal models for all important lower-extremity muscles with as few as eight carefully chosen EMG signals and eventually contribute to the design of personalized rehabilitation and surgical interventions for mobility impairments.

A Dynamic Elbow Testing Apparatus for Simulating Elbow Joint Motion in Varying Shoulder Positions

Purpose

To develop and evaluate the capabilities of a dynamic elbow testing apparatus that simulates unconstrained elbow motion throughout the range of humerothoracic (HTA) abduction.

Methods

Elbow flexion was generated by six computer-controlled electromechanical actuators that simulated muscle action, while six degree-of-freedom joint motion was measured using an optical tracking device. Repeatability of joint kinematics was assessed at four HTA angles (0°, 45°, 90°, 135°) and with two muscle force combinations (A1-biceps brachialis, brachioradialis and A2-biceps, brachioradialis). Repeatability was determined by comparing kinematics at every 10° of flexion over five flexion-extension cycles (0° to 100°).

Results

Multiple muscle force combinations can be used at each HTA angle to generate elbow flexion. Trials showed that the testing apparatus produced highly repeatable joint motion at each HTA angle and with varying muscle force combinations. The intraclass correlation coefficient was greater than 0.95 for all conditions.

Conclusions

Repeatable smooth cadaveric elbow motion was created that mimicked the in vivo situation.

Clinical relevance

These results suggest that the dynamic elbow testing apparatus can be used to characterize elbow biomechanics in cadaver upper extremities.

EMG-Driven Musculoskeletal Model Calibration With Wrapping Surface Personalization

Muscle forces and joint moments estimated by electromyography (EMG)-driven musculoskeletal models are sensitive to the wrapping surface geometry defining muscle-tendon lengths and moment arms. Despite this sensitivity, wrapping surface properties are typically not personalized to subject movement data. This study developed a novel method for personalizing OpenSim cylindrical wrapping surfaces during EMG-driven model calibration. To avoid the high computational cost of repeated OpenSim muscle analyses, the method uses two-level polynomial surrogate models. Outer-level models specify time-varying muscle-tendon lengths and moment arms as functions of joint angles, while inner-level models specify time-invariant outer-level polynomial coefficients as functions of wrapping surface parameters. To evaluate the method, we used walking data collected from two individuals post-stroke and performed four variations of EMG-driven lower extremity model calibration: 1) no calibration of scaled generic wrapping surfaces (NGA), 2) calibration of outer-level polynomial coefficients for all muscles (SGA), 3) calibration of outer-level polynomial coefficients only for muscles with wrapping surfaces (LSGA), and 4) calibration of cylindrical wrapping surface parameters for muscles with wrapping surfaces (PGA). On average compared to NGA, SGA reduced lower extremity joint moment matching errors by 31%, LSGA by 24%, and PGA by 12%, with the largest reductions occurring at the hip. Furthermore, PGA reduced peak hip joint contact force by 47% bodyweight, which was the most consistent with published in vivo measurements. The proposed method for EMG-driven model calibration with wrapping surface personalization produces physically realistic OpenSim models that reduce joint moment matching errors while improving prediction of hip joint contact force.

Advances in EMG measurement techniques, analysis procedures, and the impact of muscle mechanics on future requirements for the methodology

Muscular coordination enables locomotion and interaction with the environment. For more than 50 years electromyography (EMG) has provided insights into the central nervous system control of individual muscles or muscle groups, enabling both fine and gross motor functions. This information is available either at individual motor units (Mus) level or on a more global level from the coordination of different muscles or muscle groups. In particular, non-invasive EMG methods such as surface EMG (sEMG) or, more recently, spatial mapping methods (High-Density EMG - HDsEMG) have found their place in research into biomechanics, sport and exercise, ergonomics, rehabilitation, diagnostics, and increasingly for the control of technical devices. With further technical advances and a growing understanding of the relationship between EMG and movement task execution, it is expected that with time, especially non-invasive EMG methods will become increasingly important in movement sciences. However, while the total number of publications per year on non-invasive EMG methods is growing exponentially, the number of publications on this topic in journals with a scope in movement sciences has stagnated in the last decade. This review paper contextualizes non-invasive EMG development over the last 50 years, highlighting methodological progress. Changes in research topics related to non-invasive EMG were identified. Today non-invasive EMG procedures are increasingly used to control technical devices, where muscle mechanics have a minor influence. In movement science, however, the effect of muscle mechanics on the EMG signal cannot be neglected. This explains why non-invasive EMG's relevance in movement sciences has not developed as expected.

Three-dimensional skeletal muscle architecture in the lower legs of living human infants

Little is known about the skeletal muscle architecture of living humans at birth. In this study, we used magnetic resonance imaging (MRI) to measure the volumes of ten muscle groups in the lower legs of eight human infants aged less than three months. We then combined MRI and diffusion tensor imaging (DTI) to provide detailed, high-resolution reconstructions and measurements of moment arms, fascicle lengths, physiological cross-sectional areas (PCSAs), pennation angles and diffusion parameters of the medial (MG) and lateral gastrocnemius (LG) muscles. On average, the total lower leg muscle volume was 29.2 cm³. The largest muscle was the soleus muscle with a mean volume of 6.5 cm³. Compared to the LG muscles, the MG muscles had, on average, greater volumes (by ∼35%) and greater PCSAs (by ∼63%) but similar ankle-to-knee moment arm ratios (∼0.1 difference), fascicle lengths (∼5.7 mm difference) and pennation angles (∼2.7° difference). The MG data were compared with data previously collected from adults. The MG muscles of adults had, on average, a 63-fold greater volume, a 36-fold greater PCSA, and 1.7-fold greater fascicle length. This study demonstrates the feasibility of using MRI and DTI to reconstruct the three-dimensional architecture of skeletal muscles in living human infants. It is shown that, between infancy and adulthood, MG muscle fascicles grow primarily in cross-section rather than in length.

Modern three-dimensional digital methods for studying locomotor biomechanics in tetrapods

Here, we review the modern interface of three-dimensional (3D) empirical (e.g. motion capture) and theoretical (e.g. modelling and simulation) approaches to the study of terrestrial locomotion using appendages in tetrapod vertebrates. These tools span a spectrum from more empirical approaches such as XROMM, to potentially more intermediate approaches such as finite element analysis, to more theoretical approaches such as dynamic musculoskeletal simulations or conceptual models. These methods have much in common beyond the importance of 3D digital technologies, and are powerfully synergistic when integrated, opening a wide range of hypotheses that can be tested. We discuss the pitfalls and challenges of these 3D methods, leading to consideration of the problems and potential in their current and future usage. The tools (hardware and software) and approaches (e.g. methods for using hardware and software) in the 3D analysis of tetrapod locomotion have matured to the point where now we can use this integration to answer questions we could never have tackled 20 years ago, and apply insights gleaned from them to other fields.

Towards modern understanding of the Achilles tendon properties in human movement research

The Achilles tendon (AT) is the strongest tendon in humans, yet it often suffers from injury. The mechanical properties of the AT afford efficient movement, power amplification and power attenuation during locomotor tasks. The properties and the unique structure of the AT as a common tendon for three muscles have been studied frequently in humans using in vivo methods since 1990's. As a part of the celebration of 50 years history of the International Society of Biomechanics, this paper reviews the history of the AT research focusing on its mechanical properties in humans. The questions addressed are: What are the most important mechanical properties of the Achilles tendon, how are they studied, what is their significance to human movement, and how do they adapt? We foresee that the ongoing developments in experimental methods and modeling can provide ways to advance knowledge of the complex three-dimensional structure and properties of the Achilles tendon in vivo, and to enable monitoring of the loading and recovery for optimizing individual adaptations.

A Shoulder Musculoskeletal Model with Three-Dimensional Complex Muscle Geometries

Muscle structure is an essential component in typical computational models of the musculoskeletal system. Almost all musculoskeletal models represent muscle geometry using a set of line segments. The straight-line approach limits models' ability to accurately predict the paths of muscles with complex geometry. This approach needs knowledge of how the muscle changes shape and interacts with fundamental structures like muscles, bones, and joints that move. Moreover, the moment arms are supposed to be equivalent to all the fibers in the muscle. This study aims to create a shoulder musculoskeletal model that includes complex muscle geometries. We reconstructed the shape of fibers in the entire volume of six muscles adjacent to the shoulder using an automated technique. This method generates many fibers from the surface geometry of the skeletal muscle and its attachment areas. Highly discretized muscle representations for all muscles were created and used to simulate different shoulder movements. The moment arms of each muscle were calculated and validated against cadaveric measurements and models of the same muscles from the literature. We found that simulations using the developed musculoskeletal models generated more realistic geometries, which expands the physical representation of muscles compared to line segments. The shoulder musculoskeletal model with complex muscle geometry is created to increase the anatomical reality of models and the lines action of muscle fibers, and to be used for finite element investigations.

The load borne by the Achilles tendon during exercise: A systematic review of normative values

The Achilles tendon (AT) can be exposed to considerable stress during athletic activities and is often subject to pathologies such as tendinopathies. When designing a prevention or rehabilitation protocol, mechanical loading is a key factor to consider. This implies being able to accurately determine the load applied to the AT when performing exercises that stress this tendon. A systematic review was performed to synthesize the load borne by the AT during exercises/activities. Three databases (Pubmed, Embase and Cochrane) were searched for articles up to May 2021, and only the studies assessing the AT load in newtons relative to body-weight (BW) on humans during activities or exercises were included. Most of the 11 included studies assessed AT load when running or walking (N = 10), and only three tested exercises were usually performed during rehabilitation. The load on the tendon ranged from 2.7 to 3.95 BW when walking, from 4.15 to 7.71 BW when running, and from 0.41 to 7.3 BW according to the strengthening exercise performed. From the collected data, a progression of exercises progressively loading the Achilles tendon, as well as the possible connections with walking and running activities, could be defined. However, the trends highlighted in the relationship between tendon loading and walking or running speeds present some inconsistencies. Further research is still needed to clarify them, but also to complete the data set in healthy and injured people.

In Achilles Tendinopathy the Symptomatic Tendon Differs from the Asymptomatic Tendon While Exercise Therapy Has Little Effect on Asymmetries-An Ancillary Analysis of Data from a Controlled Clinical Trial

BACKGROUND

As inter-limb asymmetries can be associated with higher injury risk, we aimed to investigate their role in Achilles tendinopathy patients.

METHODS

In Achilles tendinopathy patients (n = 41), we assessed inter-limb asymmetries of mechanical, material, and morphological musculoskeletal properties and function and how those were affected by 12 weeks of exercise intervention (high-load protocol, n = 13; Alfredson protocol, n = 11). Moreover, we assessed whether asymmetry reductions correlated with improved Patient-Reported Outcomes (VISA-A score).

RESULTS

At baseline, tendinopathic tendons demonstrated lower tendon force (p = 0.017), lower tendon stress (p < 0.0001), larger tendon cross-sectional area (CSA) (p < 0.001), and increased intratendinous (p = 0.042) and tendon overall (p = 0.021) vascularization. For the high-load group, PRE-to-POST asymmetry comparisons revealed an asymmetry increase for the counter-movement jump (CMJ) (p = 0.034) and PRE-to-POST VISA-A score improvements correlated with CSA asymmetry reductions (p = 0.024). Within the Alfredson group, PRE-to-POST VISA-A score improvements correlated with CMJ asymmetry reductions (p = 0.044) and tendon stiffness asymmetry increases (p = 0.037). POST-to-POST in-between group comparisons revealed lower asymmetry in the high-load group for tendon elongation (p = 0.021) and tendon strain (p = 0.026).

CONCLUSIONS

The tendinopathic limb differs from the asymptomatic limb while therapeutic exercise interventions have little effect on asymmetries. Asymmetry reductions are not necessarily associated with tendon health improvements.

Evidence-Based High-Loading Tendon Exercise for 12 Weeks Leads to Increased Tendon Stiffness and Cross-Sectional Area in Achilles Tendinopathy: A Controlled Clinical Trial

BACKGROUND

Assuming that the mechanisms inducing adaptation in healthy tendons yield similar responses in tendinopathic tendons, we hypothesized that a high-loading exercise protocol that increases tendon stiffness and cross-sectional area in male healthy Achilles tendons may also induce comparable beneficial adaptations in male tendinopathic Achilles tendons in addition to improving pain and function.

OBJECTIVES

We investigated the effectiveness of high-loading exercise in Achilles tendinopathy in terms of inducing mechanical (tendon stiffness, maximum strain), material (Young's modulus), morphological (tendon cross-sectional area (CSA)), maximum voluntary isometric plantar flexor strength (MVC) as well as clinical adaptations (Victorian Institute of Sports Assessment-Achilles (VISA-A) score and pain (numerical rating scale (NRS))) as the primary outcomes. As secondary outcomes, drop (DJ) and counter-movement jump (CMJ) height and intratendinous vascularity were assessed.

METHODS

We conducted a controlled clinical trial with a 3-month intervention phase. Eligibility criteria were assessed by researchers and medical doctors. Inclusion criteria were male sex, aged between 20 and 55 years, chronic Achilles tendinopathy confirmed by a medical doctor via ultrasound-assisted assessment, and a severity level of less than 80 points on the VISA-A score. Thirty-nine patients were assigned by sequential allocation to one of three parallel arms: a high-loading intervention (training at ~ 90% of the MVC) (n = 15), eccentric exercise (according to the Alfredson protocol) as the standard therapy (n = 15) and passive therapy (n = 14). Parameters were assessed pre- and-post-intervention. Data analysis was blinded.

RESULTS

Primary outcomes: Plantar flexor MVC, tendon stiffness, mean CSA and maximum tendon strain improved only in the high-loading intervention group by 7.2 ± 9.9% (p = 0.045), 20.1 ± 20.5% (p = 0.049), 8.98 ± 5.8% (p < 0.001) and -12.4 ± 10.3% (p = 0.001), respectively. Stiffness decreased in the passive therapy group (-7.7 ± 21.2%; p = 0.042). There was no change in Young's modulus in either group (p > 0.05). The VISA-A score increased in all groups on average by 19.8 ± 15.3 points (p < 0.001), while pain (NRS) dropped by -0.55 ± 0.9 points (p < 0.001).

SECONDARY OUTCOMES

CMJ height decreased for all groups (-0.63 ± 4.07 cm; p = 0.005). There was no change in DJ height and vascularity (p > 0.05) in either group.

CONCLUSION

Despite an overall clinical improvement, it was exclusively the high-loading intervention that induced significant mechanical and morphological adaptations of the plantar flexor muscle-tendon unit. This might contribute to protecting the tendon from strain-induced injury. Thus, we recommend the high-loading intervention as an effective (alternative) therapeutic protocol in Achilles tendinopathy rehabilitation management in males.

CLINICAL TRIALS REGISTRATION NUMBER

NCT02732782.

EMG-driven musculoskeletal model calibration with estimation of unmeasured muscle excitations via synergy extrapolation

Subject-specific electromyography (EMG)-driven musculoskeletal models that predict muscle forces have the potential to enhance our knowledge of internal biomechanics and neural control of normal and pathological movements. However, technical gaps in experimental EMG measurement, such as inaccessibility of deep muscles using surface electrodes or an insufficient number of EMG channels, can cause difficulties in collecting EMG data from muscles that contribute substantially to joint moments, thereby hindering the ability of EMG-driven models to predict muscle forces and joint moments reliably. This study presents a novel computational approach to address the problem of a small number of missing EMG signals during EMG-driven model calibration. The approach (henceforth called "synergy extrapolation" or SynX) linearly combines time-varying synergy excitations extracted from measured muscle excitations to estimate 1) unmeasured muscle excitations and 2) residual muscle excitations added to measured muscle excitations. Time-invariant synergy vector weights defining the contribution of each measured synergy excitation to all unmeasured and residual muscle excitations were calibrated simultaneously with EMG-driven model parameters through a multi-objective optimization. The cost function was formulated as a trade-off between minimizing joint moment tracking errors and minimizing unmeasured and residual muscle activation magnitudes. We developed and evaluated the approach by treating a measured fine wire EMG signal (iliopsoas) as though it were "unmeasured" for walking datasets collected from two individuals post-stroke-one high functioning and one low functioning. How well unmeasured muscle excitations and activations could be predicted with SynX was assessed quantitatively for different combinations of SynX methodological choices, including the number of synergies and categories of variability in unmeasured and residual synergy vector weights across trials. The two best methodological combinations were identified, one for analyzing experimental walking trials used for calibration and another for analyzing experimental walking trials not used for calibration or for predicting new walking motions computationally. Both methodological combinations consistently provided reliable and efficient estimates of unmeasured muscle excitations and activations, muscle forces, and joint moments across both subjects. This approach broadens the possibilities for EMG-driven calibration of muscle-tendon properties in personalized neuromusculoskeletal models and may eventually contribute to the design of personalized treatments for mobility impairments.

Effective stretching position for the posterior deltoid muscle evaluated by shear wave elastography

BACKGROUND

Deteriorated extensibility of the posterior deltoid muscle is one of the factors of posterior shoulder tightness, and improvement in its extensibility is needed. However, no study has investigated which shoulder positions effectively stretch the posterior deltoid muscle in vivo. The aim of this study was to verify the effective stretching position of the posterior deltoid muscle in vivo by shear wave elastography.

METHODS

Fifteen healthy men participated in this study. The shear modulus of the posterior deltoid was measured at resting and 13 stretching positions: 60°, 90°, and 120° shoulder flexion; maximum shoulder flexion, horizontal adductions at 60°, 90°, and 120° shoulder flexion; internal rotations at 60°, 90°, and 120° shoulder flexion; and combinations of horizontal adduction with internal rotation at 60°, 90°, and 120° shoulder flexion. The shear moduli of each stretching position were compared to those of the rest. Then, among the stretching positions for which the shear modulus was significantly different from the rest, the shear moduli were compared using a three-way analysis of variance with repeated measures of the 3 factors-flexion, horizontal adduction, and internal rotation.

RESULTS

The shear moduli in all stretching positions were significantly higher than those of the rest, except for maximum shoulder flexion. The three-way analysis of variance with repeated measures revealed significant main effects in flexion and horizontal adduction. Comparing the flexion angles, the shear modulus was significantly higher at 90° than that at 60° and 120°. The shear modulus with horizontal adduction was significantly higher than that without horizontal adduction. Moreover, a significant two-way interaction was found only at flexion and horizontal adduction. The shear modulus with horizontal adduction was significantly higher at all angles than that without horizontal adduction at each flexion angle. Comparing the flexion angles with horizontal adduction, the shear modulus was significantly higher at 90° than that at 60° and 120°. No significant three-way interactions were found.

CONCLUSION

Shoulder flexion and horizontal adduction affected the extensibility of the posterior deltoid muscle, whereas the effect of shoulder internal rotation was limited. More precisely, maximal horizontal adduction at 90° shoulder flexion was the most effective stretching position for the posterior deltoid muscle.

Achilles tendon mechanical properties during walking and running are underestimated when its curvature is not accounted for

Achilles tendon (AT) mechanical properties can be estimated using an inverse dynamic approach, taking into account the tendon internal moment arm (IMA) and its kinematic behavior. Although AT presents a curvilinear line of action, a straight-line function to estimate IMA and AT length is often utilized in the literature. In this study, we combined kinetic, kinematic and ultrasound data to understand the impact of two different approaches (straight-line vs. curvilinear) in determining AT mechanical properties in vivo (during walking and running at the self-selected speed). AT force and power were calculated based on data of AT IMA and AT length derived by both respective methods. All investigated parameters were significantly affected by the method utilized (paired t-test; p < 0.05): when using the curvilinear method IMA was about 5% lower and AT length about 1.2% higher, whereas peak and mean values of AT force and power were 5% higher when compared to the straight-line method (both in walking and running). Statistic-parametric mapping (SMP) analysis revealed significant differences in IMA during the early and the late stance phase of walking and during the late stance phase of running (p < 0.01); SPM revealed significant differences also in AT length during the entire stance phase in both locomotion modes (p < 0.01). These results confirm and extend previous findings to human locomotion: neglecting the AT curvature might be a source of error, resulting in underestimates not only of internal moment arm and tendon length, but also of tendon force and power.

Percutaneous Sonographically Guided Release of Carpal Tunnel and Trigger Finger: Biomechanics, Clinical Results, Technical Developments

The interventional use of sonography is growing fast, and percutaneous sonographically guided release is more and more used as minimally invasive treatment of carpal tunnel syndrome as well as trigger finger digit. The benefits of these procedures seem promising in clinical studies, but biomechanical studies comparing these procedures with open classical surgery are scarce. Minimally invasive releases of carpal tunnel and trigger finger could limit the phenomenon of tendon bowstringing observed after open surgery. A new model is presented to compare the biomechanical effects of open and sono-guided carpal tunnel and trigger finger releases.

An Automatic and Simplified Approach to Muscle Path Modeling

This paper aims at proposing an automatic method to design and adjust simplified muscle paths of a musculoskeletal model. These muscle paths are composed of straight lines described by a limited set of fixed active via points and an optimization routine is developed to place these via points on the model in order to fit moment arms and musculotendon lengths input data. The method has been applied to a forearm musculoskeletal model extracted from the literature, using theoretical input data as an example. Results showed that for 75% of the muscle set, the relative root-mean-square error between literature theoretical data and the results from optimized muscle path was under 29.23% for moment arms and of 1.09% for musculotendon lengths. These results confirm the ability of the method to automatically generate computationally efficient muscle paths for musculoskeletal simulations. Using only via points lowers computational expense compared to paths exhibiting wrapping objects. A proper balance between computational time and anatomical realism should be found to help those models being interpreted by practitioners.

Control of Mammalian Locomotion by Somatosensory Feedback

When animals walk overground, mechanical stimuli activate various receptors located in muscles, joints, and skin. Afferents from these mechanoreceptors project to neuronal networks controlling locomotion in the spinal cord and brain. The dynamic interactions between the control systems at different levels of the neuraxis ensure that locomotion adjusts to its environment and meets task demands. In this article, we describe and discuss the essential contribution of somatosensory feedback to locomotion. We start with a discussion of how biomechanical properties of the body affect somatosensory feedback. We follow with the different types of mechanoreceptors and somatosensory afferents and their activity during locomotion. We then describe central projections to locomotor networks and the modulation of somatosensory feedback during locomotion and its mechanisms. We then discuss experimental approaches and animal models used to investigate the control of locomotion by somatosensory feedback before providing an overview of the different functional roles of somatosensory feedback for locomotion. Lastly, we briefly describe the role of somatosensory feedback in the recovery of locomotion after neurological injury. We highlight the fact that somatosensory feedback is an essential component of a highly integrated system for locomotor control. © 2021 American Physiological Society. Compr Physiol 11:1-71, 2021.

Validation of an Echidna Forelimb Musculoskeletal Model Using XROMM and diceCT

In evolutionary biomechanics, musculoskeletal computer models of extant and extinct taxa are often used to estimate joint range of motion (ROM) and muscle moment arms (MMAs), two parameters which form the basis of functional inferences. However, relatively few experimental studies have been performed to validate model outputs. Previously, we built a model of the short-beaked echidna (Tachyglossus aculeatus) forelimb using a traditional modelling workflow, and in this study we evaluate its behaviour and outputs using experimental data. The echidna is an unusual animal representing an edge-case for model validation: it uses a unique form of sprawling locomotion, and possesses a suite of derived anatomical features, in addition to other features reminiscent of extinct early relatives of mammals. Here we use diffusible iodine-based contrast-enhanced computed tomography (diceCT) alongside digital and traditional dissection to evaluate muscle attachments, modelled muscle paths, and the effects of model alterations on the MMA outputs. We use X-ray Reconstruction of Moving Morphology (XROMM) to compare ex vivo joint ROM to model estimates based on osteological limits predicted via single-axis rotation, and to calculate experimental MMAs from implanted muscles using a novel geometric method. We also add additional levels of model detail, in the form of muscle architecture, to evaluate how muscle torque might alter the inferences made from MMAs alone, as is typical in evolutionary studies. Our study identifies several key findings that can be applied to future models. 1) A light-touch approach to model building can generate reasonably accurate muscle paths, and small alterations in attachment site seem to have minimal effects on model output. 2) Simultaneous movement through multiple degrees of freedom, including rotations and translation at joints, are necessary to ensure full joint ROM is captured; however, single-axis ROM can provide a reasonable approximation of mobility depending on the modelling objectives. 3) Our geometric method of calculating MMAs is consistent with model-predicted MMAs calculated via partial velocity, and is a potentially useful tool for others to create and validate musculoskeletal models. 4) Inclusion of muscle architecture data can change some functional inferences, but in many cases reinforced conclusions based on MMA alone.

Muscle-specific economy of force generation and efficiency of work production during human running

Human running features a spring-like interaction of body and ground, enabled by elastic tendons that store mechanical energy and facilitate muscle operating conditions to minimize the metabolic cost. By experimentally assessing the operating conditions of two important muscles for running, the soleus and vastus lateralis, we investigated physiological mechanisms of muscle work production and muscle force generation. We found that the soleus continuously shortened throughout the stance phase, operating as work generator under conditions that are considered optimal for work production: high force-length potential and high enthalpy efficiency. The vastus lateralis promoted tendon energy storage and contracted nearly isometrically close to optimal length, resulting in a high force-length-velocity potential beneficial for economical force generation. The favorable operating conditions of both muscles were a result of an effective length and velocity-decoupling of fascicles and muscle-tendon unit, mostly due to tendon compliance and, in the soleus, marginally by fascicle rotation.

Calculation of flexor pollicis longus moment arm for wrist motion in a cadaver model validates the tenodesis effect for therapy

INTRODUCTION

Synergies of fingers and wrist motion have been incorporated into therapies for finger flexor tendon injuries to improve repair outcomes. Similar synergistic therapy strategies have not been well documented for the thumb.

PURPOSE OF THE STUDY

The purpose of this study was to investigate the extent to which wrist motion enables a synergistic effect at the thumb in a cadaveric model by measuring flexor pollicis longus excursion and calculating the moment arm of this tendon at the wrist joint.

STUDY DESIGN

This is a basic science research.

METHODS

Eight fresh-frozen cadaveric arms were obtained from our anatomical bequest program. The proximal arm was fixed in neutral pronation/supination position, and motion of the wrist was guided through either flexion/extension or radial/ulnar deviation. Fingers were fixed in extension, thumb interphalangeal and metacarpophalangeal joints were fixed in neutral extension, and the carpometacarpal joint was fixed at 30° palmar abduction. The flexor pollicis longus tendon was exposed proximal to the wrist crease and connected to a rotary potentiometer to measure tendon excursion. Optical markers were attached to the hand to capture kinematics. Wrists were moved from a neutral position over the range of flexion and extension and then from the neutral position through the range of radial to ulnar deviation. Moment arms were calculated.

RESULTS

Moment arm calculation indicated that the flexor pollicis longus acts as a wrist flexor over the entire motion range and as a weak radial deviator at ulnarly-deviated positions.

CONCLUSIONS

This study provides a mechanistic rationale for passive interphalangeal joint motion in varying wrist positions when treating thumb flexor tendon injuries, with benefits seen primarily for wrist extension.

Recovery from Achilles Tendon Repair: A Combination of Postsurgery Outcomes and Insufficient Remodeling of Muscle and Tendon

INTRODUCTION

Achilles tendon rupture (ATR) patients have persistent functional deficits in the triceps surae muscle-tendon unit (MTU). The complex remodeling of the MTU accompanying these deficits remains poorly understood. The purpose of the present study was to associate in vivo and in silico data to investigate the relations between changes in MTU properties and strength deficits in ATR patients.

METHODS

Eleven male subjects who had undergone surgical repair of complete unilateral ATR were examined 4.6 ± 2.0 (mean ± SD) yr after rupture. Gastrocnemius medialis (GM) tendon stiffness, morphology, and muscle architecture were determined using ultrasonography. The force-length relation of the plantar flexor muscles was assessed at five ankle joint angles. In addition, simulations (OpenSim) of the GM MTU force-length properties were performed with various iterations of MTU properties found between the unaffected and the affected side.

RESULTS

The affected side of the patients displayed a longer, larger, and stiffer GM tendon (13% ± 10%, 105% ± 28%, and 54% ± 24%, respectively) compared with the unaffected side. The GM muscle fascicles of the affected side were shorter (32% ± 12%) and with greater pennation angles (31% ± 26%). A mean deficit in plantarflexion moment of 31% ± 10% was measured. Simulations indicate that pairing an intact muscle with a longer tendon shifts the optimal angular range of peak force outside physiological angular ranges, whereas the shorter muscle fascicles and tendon stiffening seen in the affected side decrease this shift, albeit incompletely.

CONCLUSIONS

These results suggest that the substantial changes in MTU properties found in ATR patients may partly result from compensatory remodeling, although this process appears insufficient to fully restore muscle function.

A Home-Based Bilateral Rehabilitation System With sEMG-based Real-Time Variable Stiffness

Bilateral rehabilitation allows patients with hemiparesis to exploit the cooperative capabilities of both arms to promote the recovery process. Although various approaches have been proposed to facilitate synchronized robot-assisted bilateral movements, few studies have focused on addressing the varying joint stiffness resulting from dynamic motions. This paper presents a novel bilateral rehabilitation system that implements a surface electromyography (sEMG)-based stiffness control to achieve real-time stiffness adjustment based on the user's dynamic motion. An sEMG-driven musculoskeletal model that incorporates the muscle activation and muscular contraction dynamics is developed to provide reference signals for the robot's real-time stiffness control. Preliminary experiments were conducted to evaluate the system performance in tracking accuracy and comfortability, which showed the proposed rehabilitation system with sEMG-based real-time stiffness variation achieved fast adaption to the patient's dynamic movement as well as improving the comfort in robot-assisted bilateral training.

Regional differential stretching of the pectoralis major muscle: An ultrasound elastography study

Pectoralis major (PMa) muscle injuries are becoming more prevalent, and their incidence differs among the PMa regions, i.e., the clavicular, sternal, and abdominal regions. Therefore, identifying the position for effectively lengthening each PMa region is critical in preventing PMa injuries. The purpose of this study was to determine the effective stretching position for each PMa region through shear wave elastography, which can indirectly assess individual muscle lengthening. Fifteen men participated in this study. Twelve stretching positions were compounded with shoulder abductions (45°, 90°, and 135°), pelvic rotation (with or without), shoulder external rotation (with or without), and shoulder horizontal abductions. The shear modulus of each PMa region was measured through shear wave elastography in the stretching positions mentioned above. At the clavicular region, the shear modulus was higher for three stretching positions: shoulder horizontal abduction at 45° abduction during pelvic rotation and shoulder external rotation, shoulder horizontal abduction at 90° abduction, and shoulder horizontal abduction at 90° abduction while considering shoulder external rotation. For the sternal region, the shear modulus was higher in two stretching positions: shoulder horizontal abduction at 90° abduction while adding external rotation, and combination of pelvic rotation and external rotation. For the abdominal region, the shear modulus was higher in the shoulder horizontal abduction at 135° abduction with pelvic and external rotation. These results indicated that the effective stretching position was different for each PMa region.

The biomechanical role of the lacertus fibrosus of the biceps brachii Muscle

PURPOSE

The lacertus fibrosus (LF) is involved in various surgeries. However, the biomechanical contribution of the LF remains unclear. The aim of this study was to determine the role of the lacertus fibrosus on the elbow and forearm kinematics and on the biceps brachii muscle lever arms.

METHODS

This biomechanical study was performed on seven fresh-frozen upper limbs of cadavers. Elbow flexion, forearm supination, and biceps brachii muscle lever arms were analyzed in the intact conditions (I) and after superficial (R) and deep part (R2) of the lacertus fibrosus release, respectively.

RESULTS

Elbow flexion shows a significant difference (p < 0.0001) between I, R, R2. Abduction/adduction shows a significant difference between I-R (p < 0.0001) and I-R2 (p < 0.0001). Supination does not show a significant difference in mean maximum amplitude, but between 40 and 70%, there are significant differences. There is a significant mean decrease of lever arm in flexion (28%) and supination (50%) after superficial and deep part of the lacertus fibrosus release.

CONCLUSION

The results of this study show that the lacertus fibrosus increases the lever arm during flexion and supination. It limits the flexion and abduction of the elbow and supination of the forearm. Lacertus fibrosus maintains the rhythmicity between the elbow flexion and supination of the forearm.

LEVEL OF EVIDENCE

Basic science study, biomechanics.

Gracilis and semitendinosus moment arm decreased by fascial tissue release after hamstring harvesting surgery: a key parameter to understand the peak torque obtained to a shallow angle of the knee

PURPOSE

Semitendinosus and gracilis muscles whose tendons are used in surgical reconstruction of the anterior cruciate ligament maintain their contractile ability, but the peak torque angle of hamstring muscles shifted to a shallow angle postoperatively. The goal was to quantify the influence of the myofascial structures on instantaneous moment arms of knee muscles to attempt explaining the above-mentioned post-surgical observations.

METHODS

Hamstring harvesting procedures were performed by a senior orthopaedic surgeon on seven lower limbs from fresh-frozen specimens. Femoro-tibial kinematics and tendons excursion were simultaneously recorded at each step of the surgery.

RESULTS

No significant difference was demonstrated for instantaneous moment arm of gracilis during anterior cruciate ligament surgery (84% of the maximum intact values; P ≥ 0.05). The first significant semitendinosus moment arm decrease was observed after tendon harvesting (61% of the maximum intact values; p ≤ 0.005). After hamstring harvesting, the maximum and minimum moment arm (both gracilis and semi tendinosus) shifted to a shallow angle and 90°, respectively.

CONCLUSION

Moment arm modifications by paratenons and the loose connective tissue release are essential to understand the peak torque obtained to a shallow angle.

LEVEL OF EVIDENCE

Basic science study, biomechanics.

Quantifying mechanical loading and elastic strain energy of the human Achilles tendon during walking and running

The purpose of the current study was to assess in vivo Achilles tendon (AT) mechanical loading and strain energy during locomotion. We measured AT length considering its curve-path shape. Eleven participants walked at 1.4 m/s and ran at 2.5 m/s and 3.5 m/s on a treadmill. The AT length was defined as the distance between its origin at the gastrocnemius medialis myotendinous junction (MTJ) and the calcaneal insertion. The MTJ was tracked using ultrasonography and projected to the reconstructed skin surface to account for its misalignment. Skin-to-bone displacements were assessed during a passive rotation (5°/s) of the ankle joint. Force and strain energy of the AT during locomotion were calculated by fitting a quadratic function to the experimentally measured tendon force-length curve obtained from maximum voluntary isometric contractions. The maximum AT strain and force were affected by speed (p < 0.05, ranging from 4.0 to 4.9% strain and 1.989 to 2.556 kN), yet insufficient in magnitude to be considered as an effective stimulus for tendon adaptation. Besides the important tendon energy recoil during the propulsion phase (7.8 to 11.3 J), we found a recoil of elastic strain energy at the beginning of the stance phase of running (70-77 ms after touch down) between 1.7 ± 0.6 and 1.9 ± 1.1 J, which might be functionally relevant for running efficiency.

Exploring the functional morphology of the Gorilla shoulder through musculoskeletal modelling

Musculoskeletal computer models allow us to quantitatively relate morphological features to biomechanical performance. In non‐human apes, certain morphological features have long been linked to greater arm abduction potential and increased arm‐raising performance, compared to humans. Here, we present the first musculoskeletal model of a western lowland gorilla shoulder to test some of these long‐standing proposals. Estimates of moment arms and moments of the glenohumeral abductors (deltoid, supraspinatus and infraspinatus muscles) over arm abduction were conducted for the gorilla model and a previously published human shoulder model. Contrary to previous assumptions, we found that overall glenohumeral abduction potential is similar between Gorilla and Homo. However, gorillas differ by maintaining high abduction moment capacity with the arm raised above horizontal. This difference is linked to a disparity in soft tissue properties, indicating that scapular morphological features like a cranially oriented scapular spine and glenoid do not enhance the abductor function of the gorilla glenohumeral muscles. A functional enhancement due to differences in skeletal morphology was only demonstrated in the gorilla supraspinatus muscle. Contrary to earlier ideas linking a more obliquely oriented scapular spine to greater supraspinatus leverage, our results suggest that increased lateral projection of the greater tubercle of the humerus accounts for the greater biomechanical performance in Gorilla. This study enhances our understanding of the evolution of gorilla locomotion, as well as providing greater insight into the general interaction between anatomy, function and locomotor biomechanics.

The moment arms and leverage of the human finger muscles

The moment arm of a muscle's force represents the muscle's leverage or mechanical advantage in producing a joint moment. It is indicative of the muscle's potential to contribute to actuation of a joint in a particular joint motion direction and defines the role of the muscle, for example, as a joint flexor or abductor. The aims of this study were, firstly, to measure the moment arms of the flexor and extensor muscles of the metacarpophalangeal (MCP), proximal interphalangeal (PIP) and distal interphalangeal (DIP) joints, and the moment arms of the major abductor and adductor muscles of the metacarpophalangeal (MCP) joint of each finger in the hand; secondly, to assess the effect of change in joint angle on these moment arms; and thirdly, to determine if there are differences in a given flexor or extensor's muscle moment arms between the joints it spans on a given finger, and across its tendon slips to multiple fingers. The tendon-excursion method was used to measure instantaneous muscle moment arms in nine fresh-frozen entire forearm cadaver specimens. Joint flexion angle was found to have significant effects on the moment arms of the extensor muscles at the MCP and PIP joints (p < 0.05). In contrast, the digital flexor muscles maintained relatively constant moment arms through the range of joint flexion. The moment arms of the digital flexors and extensors spanning multiple joints in a finger were largest at the MCP joints and smallest at the DIP joints. The findings demonstrate greater torque generating capacity for tasks such as grasping at the proximal interphalangeal joints, and smaller torque capacity for finer movement control at the distal interphalangeal joints. The dataset generated in this study may be useful in the development and validation of computational models used in surgical planning, and rehabilitation.

Approximating complex musculoskeletal biomechanics using multidimensional autogenerating polynomials

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

The community in the fields of biomechanics, neural engineering, and neuroscience has the need to understand and simulate realistic muscle actions in real-time. In biomechanics, the models of muscle structure have been of paramount importance for understanding the mechanical demands of movements. In neural engineering, the use of biomimetic control schemes require realistic and real-time computations with low latencies to achieve an intuitive interface with high-dimensional active prostheses or orthoses. In neuroscience, the new realization of the close relationship between neural computations and body mechanics has been promoted under the concept of neuromechanics. This concept has been instrumental in the understanding of neural computations for movement planning and execution. To enable the theoretical framework of neuromechanical computations embedded within musculoskeletal organization we propose a novel method for calculating muscle biomechanics in real-time with objective approximations that embed structural and functional attributes of simulated muscles. This description offers a scalable solution that accurately computes muscle kinematic states with real-time latencies surpassing the previous results by an order of magnitude.

Surface Electromyography Meets Biomechanics: Correct Interpretation of sEMG-Signals in Neuro-Rehabilitation Needs Biomechanical Input

Coordinated activation of muscles is the basis for human locomotion. Impaired muscular activation is related to poor movement performance and disability. To restore movement performance, information about the subject's individual muscular activation is of high relevance. Surface electromyography (sEMG) allows the pain-free assessment of muscular activation and many ready-to-use technologies are available. They enable the usage of sEMG measurements in several applications. However, due to the fact that in most rehabilitation applications dynamic conditions are analyzed, the correct interpretation of sEMG signals remains difficult which hinders the spread of sEMG in clinical applications. From biomechanics it is well-known that the sEMG signal depends on muscle fiber length, contraction velocity, contraction type and on the muscle's biomechanical moment. In non-isometric conditions these biomechanical factors have to be considered when analyzing sEMG signals. Additionally, the central nervous system control strategies used to activate synergistic and antagonistic muscles have to be taken into consideration. These central nervous system activation strategies are rarely known in physiology and are hard to manage in pathology. In this perspective report we discuss how the consideration of biomechanical factors leads to more reliable information extraction from sEMG signals and how the limitations of sEMG can be overcome in dynamic conditions. This is a prerequisite if the use of sEMG in rehabilitation applications is to extend. Examples will be given showing how the integration of biomechanical knowledge into the interpretation of sEMG helps to identify the central nervous system activation strategies involved and leads to relevant clinical information.

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.

Mechanisms of reduced plantarflexor function in Cerebral palsy: smaller triceps surae moment arm and reduced muscle force

Both muscle forces, and moment arm (MA) could contribute to reduced muscle moment in people with Cerebral Palsy (CP). Current reports in CP are conflicting. The tendon travel method of estimating MA requires constant force, but passive force is high and variable in CP, and range of motion is limited. Therefore, the purpose of this study was to examine triceps surae muscle MA in 12 subjects with mild to moderate CP (15-32 years) and 10 typically developing peers (TD, 17-26 years) by tendon travel and by visually measuring the apparent MA. MA was calculated at 90° and at a reference angle (∼106°) with zero net passive moment. The tendon travel (28.8 ± 5.6 mm) and visual methods (29.1 ± 5.5 mm) yielded similar MA in CP (p = 0.94) at the reference angle. TD had significantly larger triceps surae muscle MA than CP subjects (p = 0.002), 35.4 ± 4.1 mm at the reference angle for tendon travel and 35.4 ± 3.6 mm by the visual method. Test/retest revealed less bias (0.8 mm) using the visual method. Calculated active peak isometric force was significantly less in CP (1983.8 ± 887.0 N) than TD (4104.9 ± 1154.9 N, p < 0.001). There are challenges in estimating MA in CP, but the visual method is more reliable. Although a shorter moment arm would reduce the joint moment, joint angular velocity for a given velocity of muscle shortening would be enhanced. Strength training may mitigate the effects of the smaller moment arm and reduced joint moment generated in those with CP.

Influence of joint angle on muscle fascicle dynamics and rate of torque development during isometric explosive contractions

Ankle angle influences the operating muscle fascicle lengths of gastrocnemius medialis and the rate of torque development during explosive isometric plantar flexions. The rate of torque development peaks in neutral angles where muscle fascicles shorten over the plateau of the force-length relationship. When fascicles operate over the plateau of the force-length relationship (neutral ankle positions), the force-velocity properties represent a limiting factor for the rapid force-generating capacity from 100 ms after the onset of explosive contractions.

Torus obstacle method as a wrapping approach of the deltoid muscle group for humeral abduction in musculoskeletal simulation

Musculoskeletal models of the shoulder complex are valuable research aids to investigate tears of the supraspinatus and the resulting mechanical impact during abduction of the humerus. One of the major contributors to this motion is the deltoid muscle group and for this, an accurate modeling of the lines of action is indispensable. The aim of this work was to utilize a torus obstacle wrapping approach for the deltoids of an existing shoulder model and assess the feasibility of the approach during humeral abduction. The shoulder model from the AnyBody™ modeling system was used as a platform. The size of the tori is based on a magnetic resonance imaging (MRI) approach and several kinematic couplings are implemented to determine the trajectories of the tori during abduction. To assess the model behavior, the moment arms of the virtual muscle elements and the resultant glenohumeral joint reaction force (GHJF) were compared with reference data from the literature during abduction of the humerus in the range 20°-120°. The root mean square error for the anterior, lateral and posterior part between the simulated muscle elements and reference data from the literature was 3.9, 1.7 and 5.8 mm, respectively. The largest deviation occurred on the outer elements of the muscle groups, with 12.6, 10.4 and 20.5 mm, respectively. During abduction, there is no overlapping of the muscle elements and these are in continuous contact with the torus obstacles, thus enabling a continuous force transmission. This results in a rising trend of the resultant GHJF. The torus obstacle approach as a wrapping method for the deltoid muscles provides a guided muscle pathing by simultaneously approximating the curvature of the deltoid muscle. The results from the comparison of the simulated moment arms and the resultant GHJF are in accordance with those in the literature in the range 20°-120° of abduction. Although this study shows the strength of the torus obstacle as a wrapping approach, the method of fitting the tori according to MRI data was not suitable. A cadaver study is recommended to better validate and mathematically describe the torus approach.

Voluntary Control of an Ankle Joint Exoskeleton by Able-Bodied Individuals and Stroke Survivors Using EMG-Based Admittance Control Scheme

Control schemes based on electromyography (EMG) have demonstrated their superiority in human-robot cooperation due to the fact that motion intention can be well estimated by EMG signals. However, there are several limitations due to the noisy nature of EMG signals and the inaccuracy of EMG-force/torque estimation, which might deteriorate the stability of human-robot cooperation movement. To improve the movement stability, an EMG-based admittance control scheme (EACS) was proposed, comprised of an EMG-driven musculoskeletal model (EDMM), an admittance filter and an inner position controller. To investigate the performance of EACS, a series of sinusoidal tracking tasks were conducted with 12 healthy participants and 4 stroke survivors in an ankle exoskeleton in comparison with the EMG-based open-loop control scheme (EOCS). The experimental results indicated that both EACS and EOCS could improve stroke survivors' ankle range of motion (ROM). The experimental results of both healthy participants and stroke survivors showed that the assistance torque, tracking error and jerk values of EACS were lower than those of EOCS. The interaction torque of EACS decreased towards the increasing assistance ratio while that of EOCS increased. Moreover, the EMG levels of tibialis anterior (TA) decreased towards the increasing assistance ratio but were higher than those of EOCS. EACS was effective in improving movements stability, and had the potential to be applied in robot-assisted rehabilitation training to address the foot-drop problem.

Automated Generation of Three-Dimensional Complex Muscle Geometries for Use in Personalised Musculoskeletal Models

The geometrical representation of muscles in computational models of the musculoskeletal system typically consists of a series of line segments. These muscle anatomies are based on measurements from a limited number of cadaveric studies that recently have been used as atlases for creating subject-specific models from medical images, so potentially restricting the options for personalisation and assessment of muscle geometrical models. To overcome this methodological limitation, we propose a novel, completely automated technique that, from a surface geometry of a skeletal muscle and its attachment areas, can generate an arbitrary number of lines of action (fibres) composed by a user-defined number of straight-line segments. These fibres can be included in standard musculoskeletal models and used in biomechanical simulations. This methodology was applied to the surfaces of four muscles surrounding the hip joint (iliacus, psoas, gluteus maximus and gluteus medius), segmented on magnetic resonance imaging scans from a cadaveric dataset, for which highly discretised muscle representations were created and used to simulate functional tasks. The fibres’ moment arms were validated against measurements and models of the same muscles from the literature with promising outcomes. The proposed approach is expected to improve the anatomical representation of skeletal muscles in personalised biomechanical models and finite element applications.

Direct Validation of Model-Predicted Muscle Forces in the Cat Hindlimb During Locomotion

Various methods are available for simulating the movement patterns of musculoskeletal systems and determining individual muscle forces, but the results obtained from these methods have not been rigorously validated against experiment. The aim of this study was to compare model predictions of muscle force derived for a cat hindlimb during locomotion against direct measurements of muscle force obtained in vivo. The cat hindlimb was represented as a 5-segment, 13-degrees-of-freedom (DOF), articulated linkage actuated by 25 Hill-type muscle-tendon units (MTUs). Individual muscle forces were determined by combining gait data with two widely used computational methods—static optimization and computed muscle control (CMC)—available in opensim, an open-source musculoskeletal modeling and simulation environment. The forces developed by the soleus, medial gastrocnemius (MG), and tibialis anterior muscles during free locomotion were measured using buckle transducers attached to the tendons. Muscle electromyographic activity and MTU length changes were also measured and compared against the corresponding data predicted by the model. Model-predicted muscle forces, activation levels, and MTU length changes were consistent with the corresponding quantities obtained from experiment. The calculated values of muscle force obtained from static optimization agreed more closely with experiment than those derived from CMC.

Area 2 of primary somatosensory cortex encodes kinematics of the whole arm

Proprioception, the sense of body position, movement, and associated forces, remains poorly understood, despite its critical role in movement. Most studies of area 2, a proprioceptive area of somatosensory cortex, have simply compared neurons' activities to the movement of the hand through space. Using motion tracking, we sought to elaborate this relationship by characterizing how area 2 activity relates to whole arm movements. We found that a whole-arm model, unlike classic models, successfully predicted how features of neural activity changed as monkeys reached to targets in two workspaces. However, when we then evaluated this whole-arm model across active and passive movements, we found that many neurons did not consistently represent the whole arm over both conditions. These results suggest that 1) neural activity in area 2 includes representation of the whole arm during reaching and 2) many of these neurons represented limb state differently during active and passive movements.

A Multi-Mode Rehabilitation Robot With Magnetorheological Actuators Based on Human Motion Intention Estimation

Lower extremity paralysis has become common in recent years, and robots have been developed to help patients recover from it. This paper presents such a robotic system that allows for two working modes, the robot-active mode and human-active mode. The robot is designed to be equipped with magnetorheological (MR) actuators that have the advantages of high torque, fast response, flexible controllability, low power consumption and safety guarantee. The design and characteristics of the MR actuator are introduced. In the robot-active mode, the MR actuator works as a clutch to transfer the torque to the robotic joint safely. In the human-active mode, the MR actuator functions as a brake to provide resistance to help strengthen muscles. The working mode is determined by the human motion intention, which is detected via the skin surface electromyography (EMG) signals. The human-robot interaction torques are estimated using the EMG-driven impedance model. The biomechanical analysis based on AnyBody Modeling System (AMS) is used to help optimization. Then, an adaptive control method is proposed to realize the assist-as-needed (AAN) training strategy, where the robot can switch between these two modes. Experiments are conducted to validate the proposed design.

Posterior and inferior glenosphere position in reverse total shoulder arthroplasty supports deltoid efficiency for shoulder flexion and elevation

BACKGROUND

For humeral flexion and elevation, most relevant for daily activities with reverse total shoulder arthroplasty, the anterior and lateral deltoid muscles are most important. However, how this direction of movement is best supported with the glenosphere position is not fully understood. We hypothesized that both inferior positioning and posterior positioning of the glenosphere may best support this direction of movement.

METHODS

A validated, anatomic biomechanical shoulder model was modified to host a reverse shoulder prosthesis. The glenoid baseplate was altered to allow inferior, lateral, and posterior center-of-rotation (COR) offsets. An optical tracking system was used to track the excursion of ropes simulating portions of various shoulder muscles during humeral abduction, elevation, and flexion.

RESULTS

The inferior COR offset resulted in a significant increase in the deltoid moment arm in all 3 planes of motion. The lateral COR offset showed a significantly lower posterior deltoid moment arm during humeral abduction and a significantly lower lateral deltoid moment arm during humeral elevation. The posterior offset showed significantly larger anterior and lateral deltoid moment arms during humeral flexion.

DISCUSSION AND CONCLUSION

Owing to the oblique direction of the deltoid muscle across the shoulder joint, an inferior offset of the COR in reverse total shoulder arthroplasty increases the deltoid moment arm during abduction, elevation, and flexion, whereas it mainly supports humeral flexion at a posterior offset. For humeral elevation and flexion, favorable positioning of the glenosphere may, therefore, be defined by a more inferior and posterior placement compared with the non-offset position.

Capturing 3D large-strain Euler-bending filament dynamics in fibrous media simulations; sample case of compression collapse in dendritic actin network

Cytoskeletal networks to transmission towers are comprised of slender elements. Slender filaments bend and buckle more easily than stretch. Therefore a deforming network is expected to exhaust all possible bending-based modes before engaging filament stretch. While the large-strain bending critically determines fibrous-media response, simulations use small-strain and jointed approximations. At low resolution, these approximations inflate bending resistance and delay buckling onset. The proposed string-of-continuous-beams (SOCB) approach captures 3D nonlinear Euler bending of filaments with high fidelity at low cost. Bending geometry (i.e. angles and its differentials) is solved as primary variables, to fit a 5th order polynomial of the contour angle. Displacement, solved simultaneously as length conservation, is predicted with C3 and C6 smoothness between and within segments, using only 2 nodes. In the chosen analysis frame, in-plane and out-plane moments can be decoupled for arbitrarily-curved segments. Complex crosslink force-transfers can be specified. Simulations show that when a daughter branch is appended, the buckling resistance of a filament changes from linear to nonlinear before reversible collapse. An actin outcrop with 8 generations of mother-daughter branching produced the linear, nonlinear, and collapse regimes observed in compression experiments. 'Collapse' was a redistribution of outcrop forces following the buckling of few strands.

Integrated System for Monitoring Muscular States during Elbow Flexor Resistance Training in Bedridden Patients

To improve or maintain the physical function of bedridden patients, appropriate and effective exercises are required during the patient's bed rest. Resistance training (RT) is an effective exercise for improving the physical function of bedridden patients, and the improvement of the physical function is caused by mechanical stimuli associated with RT. Currently, the measured mechanical stimuli are external variables which represent the synthetic effect of multiple muscles and body movements. Important features of stimuli experienced by muscles are of crucial importance in explaining muscular strength and power adaptation. This study describes an integrated system for assessing muscular states during elbow flexor resistance training in bedridden patients, and some experiments were carried out to test and evaluate this system. The integrated system incorporates an elbow joint angle estimation model (EJAEM), a musculoskeletal model (MSM), and a muscle-tendon model. The EJAEM enables real-time interaction between patient and MSM. The MSM is a three-dimensional model of the upper extremity, including major muscles that make up the elbow flexor and extensor, and was built based on public data. One set of concentric and eccentric contraction was performed by a healthy subject, and the results of the calculations were analyzed to show important features of mechanical stimuli experienced by muscles during the training. The integrated system provides a considerable method to monitor the body-level and muscle-level mechanical stimuli during elbow flexor resistance training in bedridden patients.

Analyzing Moment Arm Profiles in a Full-Muscle Rat Hindlimb Model

Understanding the kinematics of a hindlimb model is a fundamental aspect of modeling coordinated locomotion. This work describes the development process of a rat hindlimb model that contains a complete muscular system and incorporates physiological walking data to examine realistic muscle movements during a step cycle. Moment arm profiles for selected muscles are analyzed and presented as the first steps to calculating torque generation at hindlimb joints. A technique for calculating muscle moment arms from muscle attachment points in a three-dimensional (3D) space has been established. This model accounts for the configuration of adjacent joints, a critical aspect of biarticular moment arm analysis that must be considered when calculating joint torque. Moment arm profiles from isolated muscle motions are compared to two existing models. The dependence of biarticular muscle's moment arms on the configuration of the adjacent joint is a critical aspect of moment arm analysis that must be considered when calculating joint torque. The variability in moment arm profiles suggests changes in muscle function during a step.