A two-muscle, continuum-mechanical forward simulation of the upper limb

Impact of lengthening velocity on the generation of eccentric force by slow-twitch muscle fibers in long stretches

After an initial increase, isovelocity elongation of a muscle fiber can lead to diminishing (referred to as Give in the literature) and subsequently increasing force. How the stretch velocity affects this behavior in slow-twitch fibers remains largely unexplored. Here, we stretched fully activated individual rat soleus muscle fibers from 0.85 to 1.3 optimal fiber length at stretch velocities of 0.01, 0.1, and 1 maximum shortening velocity, vmax, and compared the results with those of rat EDL fast-twitch fibers obtained in similar experimental conditions. In soleus muscle fibers, Give was 7%, 18%, and 44% of maximum isometric force for 0.01, 0.1, and 1 vmax, respectively. As in EDL fibers, the force increased nearly linearly in the second half of the stretch, although the number of crossbridges decreased, and its slope increased with stretch velocity. Our findings are consistent with the concept of a forceful detachment and subsequent crossbridge reattachment in the stretch's first phase and a strong viscoelastic titin contribution to fiber force in the second phase of the stretch. Interestingly, we found interaction effects of stretch velocity and fiber type on force parameters in both stretch phases, hinting at fiber type-specific differences in crossbridge and titin contributions to eccentric force. Whether fiber type-specific combined XB and non-XB models can explain these effects or if they hint at some not fully understood properties of muscle contraction remains to be shown. These results may stimulate new optimization perspectives in sports training and provide a better understanding of structure-function relations of muscle proteins.

Determination of muscle shape deformations of the tibialis anterior during dynamic contractions using 3D ultrasound

Purpose

In this paper, we introduce a novel method for determining 3D deformations of the human tibialis anterior (TA) muscle during dynamic movements using 3D ultrasound.

Materials and Methods

An existing automated 3D ultrasound system is used for data acquisition, which consists of three moveable axes, along which the probe can move. While the subjects perform continuous plantar- and dorsiflexion movements in two different controlled velocities, the ultrasound probe sweeps cyclically from the ankle to the knee along the anterior shin. The ankle joint angle can be determined using reflective motion capture markers. Since we considered the movement direction of the foot, i.e., active or passive TA, four conditions occur: slow active, slow passive, fast active, fast passive. By employing an algorithm which defines ankle joint angle intervals, i.e., intervals of range of motion (ROM), 3D images of the volumes during movement can be reconstructed.

Results

We found constant muscle volumes between different muscle lengths, i.e., ROM intervals. The results show an increase in mean cross-sectional area (CSA) for TA muscle shortening. Furthermore, a shift in maximum CSA towards the proximal side of the muscle could be observed for muscle shortening. We found significantly different maximum CSA values between the fast active and all other conditions, which might be caused by higher muscle activation due to the faster velocity.

Conclusion

In summary, we present a method for determining muscle volume deformation during dynamic contraction using ultrasound, which will enable future empirical studies and 3D computational models of skeletal muscles.

Determining a musculoskeletal system's pre-stretched state using continuum-mechanical forward modelling and joint range optimization

The subject-specific range of motion (RoM) of a musculoskeletal joint system is balanced by pre-tension levels of individual muscles, which affects their contraction capability. Such an inherent pre-tension or pre-stretch of muscles is not measureable with in vivo experiments. Using a 3D continuum mechanical forward simulation approach for motion analysis of the musculoskeletal system of the forearm with 3 flexor and 2 extensor muscles, we developed an optimization process to determine the muscle fibre pre-stretches for an initial arm position, which is given human dataset. We used RoM values of a healthy person to balance the motion in extension and flexion. The performed sensitivity study shows that the fibre pre-stretches of the m. brachialis, m. biceps brachii and m. triceps brachii with 91 % dominate the objective flexion ratio, while m. brachiradialis and m. anconeus amount 7.8 % and 1.2 % . Within the multi-dimensional space of the surrogate model, 3D sub-spaces of primary variables, namely the dominant muscles and the global objective, flexion ratio, exhibit a path of optimal solutions. Within this optimal path, the muscle fibre pre-stretch of two flexors demonstrate a negative correlation, while, in contrast, the primary extensor, m. triceps brachii correlates positively to each of the flexors. Comparing the global optimum with four other designs along the optimal path, we saw large deviations, e.g., up to 15∘ in motion and up to 40% in muscle force. This underlines the importance of accurate determination of fibre pre-stretch in muscles, especially, their role in pathological muscular disorders and surgical applications such as free muscle or tendon transfer.

The role of pennation angle and architectural gearing to rate of force development in dynamic and isometric muscle contractions

BACKGROUND

Associations between muscle architecture and rate of force development (RFD) have been largely studied during fixed-end (isometric) contractions. Fixed-end contractions may, however, limit muscle shape changes and thus alter the relationship between muscle architecture an RFD.

AIM

We compared the correlation between muscle architecture and architectural gearing and knee extensor RFD when assessed during dynamic versus fixed-end contractions.

METHODS

Twenty-two recreationally active male runners performed dynamic knee extensions at constant acceleration (2000°s-2) and isometric contractions at a fixed knee joint angle (fixed-end contractions). Torque, RFD, vastus lateralis muscle thickness, and fascicle dynamics were compared during 0-75 and 75-150 ms after contraction onset.

RESULTS

Resting fascicle angle was moderately and positively correlated with RFD during fixed-end contractions (r = 0.42 and 0.46 from 0-75 and 75-150 ms, respectively; p < 0.05), while more strongly (p < 0.05) correlated with RFD during dynamic contractions (r = 0.69 and 0.73 at 0-75 and 75-150 ms, respectively; p < 0.05). Resting fascicle angle was (very) strongly correlated with architectural gearing (r = 0.51 and 0.73 at 0-75 ms and 0.50 and 0.70 at 75-150 ms; p < 0.05), with gearing in turn also being moderately to strongly correlated with RFD in both contraction conditions (r = 0.38-0.68).

CONCLUSION

Resting fascicle angle was positively correlated with RFD, with a stronger relationship observed in dynamic than isometric contraction conditions. The stronger relationships observed during dynamic muscle actions likely result from different restrictions on the acute changes in muscle shape and architectural gearing imposed by isometric versus dynamic muscle contractions.

Finite Element Model of the Shoulder with Active Rotator Cuff Muscles: Application to Wheelchair Propulsion

The rotator cuff is prone to injury, remarkably so for manual wheelchair users. To understand its pathomechanisms, finite element models incorporating three-dimensional activated muscles are needed to predict soft tissue strains during given tasks. This study aimed to develop such a model to understand pathomechanisms associated with wheelchair propulsion. We developed an active muscle model associating a passive fiber-reinforced isotropic matrix with an activation law linking calcium ion concentration to tissue tension. This model was first evaluated against known physiological muscle behavior; then used to activate the rotator cuff during a wheelchair propulsion cycle. Here, experimental kinematics and electromyography data was used to drive a shoulder finite element model. Finally, we evaluated the importance of muscle activation by comparing the results of activated and non-activated rotator cuff muscles during both propulsion and isometric contractions. Qualitatively, the muscle constitutive law reasonably reproduced the classical Hill model force-length curve and the behavior of a transversally loaded muscle. During wheelchair propulsion, the deformation and fiber stretch of the supraspinatus muscle-tendon unit pointed towards the possibility for this tendon to develop tendinosis due to the multiaxial loading imposed by the kinematics of propulsion. Finally, differences in local stretch and positions of the lines of action between activated and non-activated models were only observed at activation levels higher than 30%. Our novel finite element model with active muscles is a promising tool for understanding the pathomechanisms of the rotator cuff for various dynamic tasks, especially those with high muscle activation levels.

Influence of muscle length on the three-dimensional architecture and aponeurosis dimensions of rabbit calf muscles

The function of a muscle is highly dependent on its architecture, which is characterized by the length, pennation, and curvature of the fascicles, and the geometry of the aponeuroses. During in vivo function, muscles regularly undergo changes in length, thereby altering their architecture. During passive muscle lengthening, fascicle length (FL) generally increases and the angle of fascicle pennation (FP) and the fascicle curvature (FC) decrease, while the aponeuroses increase in length but decrease in width. Muscles are differently structured, making their change during muscle lengthening complex and multifaceted. To obtain comprehensive data on architectural changes in muscles during passive length, the present study determined the three-dimensional fascicle geometry of rabbit M. gastrocnemius medialis (GM), M. gastrocnemius lateralis (GL), and M. plantaris (PLA). For this purpose, the left and right legs of three rabbits were histologically fixed at targeted ankle joint angles of 95° (short muscle length [SML]) and 60° (long muscle length [LML]), respectively, and the fascicles were tracked by manual three-dimensional digitization. In a second set of experiments, the GM aponeurosis dimensions of ten legs from five rabbits were determined at varying muscle lengths via optical marker tracking. The GM consisted of a uni-pennated compartment, whereas the GL and PLA contained multiple compartments of differently pennated fascicles. In the LML compared to the SML, the GM, GL, and PLA had on average a 41%, 29%, and 41% increased fascicle length, and a 30%, 25%, and 33% decrease in fascicle pennation and a 32%, 11%, and 35% decrease in fascicle curvature, respectively. Architectural properties were also differentiated among the different compartments of the PLA and GL, allowing for a more detailed description of their fascicle structure and changes. It was shown that the compartments change differently with muscle length. It was also shown that for each degree of ankle joint angle reduction, the proximal GM aponeurosis length increased by 0.11%, the aponeurosis width decreased by 0.22%, and the area was decreased by 0.20%. The data provided improve our understanding of muscles and can be used to develop and validate muscle models.

3D ultrasound-based determination of skeletal muscle fascicle orientations

Architectural parameters of skeletal muscle such as pennation angle provide valuable information on muscle function, since they can be related to the muscle force generating capacity, fiber packing, and contraction velocity. In this paper, we introduce a 3D ultrasound-based workflow for determining 3D fascicle orientations of skeletal muscles. We used a custom-designed automated motor driven 3D ultrasound scanning system for obtaining 3D ultrasound images. From these, we applied a custom-developed multiscale-vessel enhancement filter-based fascicle detection algorithm and determined muscle volume and pennation angle. We conducted trials on a phantom and on the human tibialis anterior (TA) muscle of 10 healthy subjects in plantarflexion (157 ± 7∘ ), neutral position (109 ± 7∘ , corresponding to neutral standing), and one resting position in between (145 ± 6∘ ). The results of the phantom trials showed a high accuracy with a mean absolute error of 0.92 ± 0.59∘ . TA pennation angles were significantly different between all positions for the deep muscle compartment; for the superficial compartment, angles are significantly increased for neutral position compared to plantarflexion and resting position. Pennation angles were also significantly different between superficial and deep compartment. The results of constant muscle volumes across the 3 ankle joint angles indicate the suitability of the method for capturing 3D muscle geometry. Absolute pennation angles in our study were slightly lower than recent literature. Decreased pennation angles during plantarflexion are consistent with previous studies. The presented method demonstrates the possibility of determining 3D fascicle orientations of the TA muscle in vivo.

How mechanics of individual muscle-tendon units define knee and ankle joint function in health and cerebral palsy-a narrative review

This study reviews the relationship between muscle-tendon biomechanics and joint function, with a particular focus on how cerebral palsy (CP) affects this relationship. In healthy individuals, muscle size is a critical determinant of strength, with muscle volume, cross-sectional area, and moment arm correlating with knee and ankle joint torque for different isometric/isokinetic contractions. However, in CP, impaired muscle growth contributes to joint pathophysiology even though only a limited number of studies have investigated the impact of deficits in muscle size on pathological joint function. As muscles are the primary factors determining joint torque, in this review two main approaches used for muscle force quantification are discussed. The direct quantification of individual muscle forces from their relevant tendons through intraoperative approaches holds a high potential for characterizing healthy and diseased muscles but poses challenges due to the invasive nature of the technique. On the other hand, musculoskeletal models, using an inverse dynamic approach, can predict muscle forces, but rely on several assumptions and have inherent limitations. Neither technique has become established in routine clinical practice. Nevertheless, identifying the relative contribution of each muscle to the overall joint moment would be key for diagnosis and formulating efficient treatment strategies for patients with CP. This review emphasizes the necessity of implementing the intraoperative approach into general surgical practice, particularly for joint correction operations in diverse patient groups. Obtaining in vivo data directly would enhance musculoskeletal models, providing more accurate force estimations. This integrated approach can improve the clinicians' decision-making process and advance treatment strategies by predicting changes at the muscle and joint levels before interventions, thus, holding the potential to significantly enhance clinical outcomes.

Development and verification of a physiologically motivated internal controller for the open-source extended Hill-type muscle model in LS-DYNA

Nowadays, active human body models are becoming essential tools for the development of integrated occupant safety systems. However, their broad application in industry and research is limited due to the complexity of incorporated muscle controllers, the long simulation runtime, and the non-regular use of physiological motor control approaches. The purpose of this study is to address the challenges in all indicated directions by implementing a muscle controller with several physiologically inspired control strategies into an open-source extended Hill-type muscle model formulated as LS-DYNA user-defined umat41 subroutine written in the Fortran programming language. This results in increased usability, runtime performance and physiological accuracy compared to the standard muscle material existing in LS-DYNA. The proposed controller code is verified with extensive experimental data that include findings for arm muscles, the cervical spine region, and the whole body. Selected verification experiments cover three different muscle activation situations: (1) passive state, (2) open-loop and closed-loop muscle activation, and (3) reflexive behaviour. Two whole body finite element models, the 50th percentile female VIVA OpenHBM and the 50th percentile male THUMS v5, are used for simulations, complemented by the simplified arm model extracted from the 50th percentile male THUMS v3. The obtained results are evaluated additionally with the CORrelation and Analysis methodology and the mean squared error method, showing good to excellent biofidelity and sufficient agreement with the experimental data. It was shown additionally how the integrated controller allows simplified mimicking of the movements for similar musculoskeletal models using the parameters transfer method. Furthermore, the Hill-type muscle model presented in this paper shows better kinematic behaviour even in the passive case compared to the existing one in LS-DYNA due to its improved damping and elastic properties. These findings provide a solid evidence base motivating the application of the enhanced muscle material with the internal controller in future studies with Active Human Body Models under different loading conditions.

Linking cortex and contraction-Integrating models along the corticomuscular pathway

Computational models of the neuromusculoskeletal system provide a deterministic approach to investigate input-output relationships in the human motor system. Neuromusculoskeletal models are typically used to estimate muscle activations and forces that are consistent with observed motion under healthy and pathological conditions. However, many movement pathologies originate in the brain, including stroke, cerebral palsy, and Parkinson's disease, while most neuromusculoskeletal models deal exclusively with the peripheral nervous system and do not incorporate models of the motor cortex, cerebellum, or spinal cord. An integrated understanding of motor control is necessary to reveal underlying neural-input and motor-output relationships. To facilitate the development of integrated corticomuscular motor pathway models, we provide an overview of the neuromusculoskeletal modelling landscape with a focus on integrating computational models of the motor cortex, spinal cord circuitry, α-motoneurons  and skeletal muscle in regard to their role in generating voluntary muscle contraction. Further, we highlight the challenges and opportunities associated with an integrated corticomuscular pathway model, such as challenges in defining neuron connectivities, modelling standardisation, and opportunities in applying models to study emergent behaviour. Integrated corticomuscular pathway models have applications in brain-machine-interaction, education, and our understanding of neurological disease.

Eccentric muscle contractions: from single muscle fibre to whole muscle mechanics

Eccentric muscle loading encompasses several unique features compared to other types of contractions. These features include increased force, work, and performance at decreased oxygen consumption, reduced metabolic cost, improved energy efficiency, as well as decreased muscle activity. This review summarises explanatory approaches to long-standing questions in terms of muscular contraction dynamics and molecular and cellular mechanisms underlying eccentric muscle loading. Moreover, this article intends to underscore the functional link between sarcomeric components, emphasising the fundamental role of titin in skeletal muscle. The giant filament titin reveals versatile functions ranging from sarcomere organisation and maintenance, providing passive tension and elasticity, and operates as a mechanosensory and signalling platform. Structurally, titin consists of a viscoelastic spring segment that allows activation-dependent coupling to actin. This titin-actin interaction can explain linear force increases in active lengthening experiments in biological systems. A three-filament model of skeletal muscle force production (mediated by titin) is supposed to overcome significant deviations between experimental observations and predictions by the classic sliding-filament and cross-bridge theories. Taken together, this review intends to contribute to a more detailed understanding of overall muscle behaviour and force generation-from a microscopic sarcomere level to a macroscopic multi-joint muscle level-impacting muscle modelling, the understanding of muscle function, and disease.

Validation of a Finite Element Simulation for Predicting Individual Knee Joint Kinematics

Goal: We introduce an in-vivo validated finite element (FE) simulation approach for predicting individual knee joint kinematics. Our vision is to improve clinicians' understanding of the complex individual anatomy and potential pathologies to improve treatment and restore physiological joint kinematics. Methods: Our 3D FE modeling approach for individual human knee joints is based on segmentation of anatomical structures extracted from routine static magnetic resonance (MR) images. We validate the predictive abilities of our model using static MR images of the knees of eleven healthy volunteers in dedicated knee poses, which are achieved using a customized MR-compatible pneumatic loading device. Results: Our FE simulations reach an average translational accuracy of 2 mm and an average angular accuracy of 1[Formula: see text] compared to the reference knee pose. Conclusions: Reaching high accuracy, our individual FE model can be used in the decision-making process to restore knee joint stability and functionality after various knee injuries.

Shear wave elastography characterizes passive and active mechanical properties of biceps brachii muscle in vivo

Mechanical characterization of individual muscles in their in vivo environment is not well studied. Shear wave elastography (SWE) as a non-invasive technique was shown to be promising in quantifying the local mechanical properties of skeletal muscles. This study aimed to investigate the mechanics of the biceps brachii muscle (BB) derived from SWE in relation to elbow joint position and contraction intensity during isometric contraction. 14 healthy, young subjects participated in the study and five different joint positions (60°-180° elbow angle) were investigated. Shear elastic modulus and surface electromyography (sEMG) of the BB and elbow torque were measured simultaneously, both in passive (i.e., resting) and active states during slow, sub-maximal isometric ramp contractions up to 25%, 50%, and 75% of the maximum voluntary contraction. At passive state, the shear elastic modulus of the BB increased with increasing elbow angle (p < 0.001). Maximum elbow flexion torque was produced at 60° and it decreased with increasing elbow angle (p = 0.001). During sub-maximal contractions, both elbow angle (p < 0.001) and contraction intensity (p < 0.001) had significant effects on the shear elastic modulus but only contraction intensity (p < 0.001) affected sEMG amplitude of the BB. Although torque was decreased at extended elbow positions (150°, 180°), higher active shear elastic modulus of BB muscle was found compared to flexed positions (60°, 90°). Linear regression of the BB sEMG amplitude over elbow torque showed good agreement for all joint positions (R² between 0.69 and 0.89) while the linear agreement between shear elastic modulus of BB and elbow torque differed between flexed (R² = 0.70 at 60° and R² = 0.79 at 90°) and extended positions (with the lowest R² = 0.57 at 150°). We conclude that using SWE, we can detect length-dependent mechanical changes of BB both in passive and active states. More importantly, SWE can be used to characterize active muscle properties in vivo. The present findings have critical importance for developing muscle stiffness as a measure of individual muscle force to validate muscle models and using SWE in clinical diagnostics.

Modelling motor units in 3D: influence on muscle contraction and joint force via a proof of concept simulation

Functional heterogeneity is a skeletal muscle’s ability to generate diverse force vectors through localised motor unit (MU) recruitment. Existing 3D macroscopic continuum-mechanical finite element (FE) muscle models neglect MU anatomy and recruit muscle volume simultaneously, making them unsuitable for studying functional heterogeneity. Here, we develop a method to incorporate MU anatomy and information in 3D models. Virtual fibres in the muscle are grouped into MUs via a novel “virtual innervation” technique, which can control the units’ size, shape, position, and overlap. The discrete MU anatomy is then mapped to the FE mesh via statistical averaging, resulting in a volumetric MU distribution. Mesh dependency is investigated using a 2D idealised model and revealed that the amount of MU overlap is inversely proportional to mesh dependency. Simultaneous recruitment of a MU’s volume implies that action potentials (AP) propagate instantaneously. A 3D idealised model is used to verify this assumption, revealing that neglecting AP propagation results in a slightly less-steady force, advanced in time by approximately 20 ms, at the tendons. Lastly, the method is applied to a 3D, anatomically realistic model of the masticatory system to demonstrate the functional heterogeneity of masseter muscles in producing bite force. We found that the MU anatomy significantly affected bite force direction compared to bite force magnitude. MU position was much more efficacious in bringing about bite force changes than MU overlap. These results highlight the relevance of MU anatomy to muscle function and joint force, particularly for muscles with complex neuromuscular architecture.

A biophysically guided constitutive law of the musculotendon-complex: modelling and numerical implementation in Abaqus

BACKGROUND AND OBJECTIVE

Many biomedical, clinical, and industrial applications may benefit from musculoskeletal simulations. Three-dimensional macroscopic muscle models (3D models) can more accurately represent muscle architecture than their 1D (line-segment) counterparts. Nevertheless, 3D models remain underutilised in academic, clinical, and commercial environments. Among the reasons for this is a lack of modelling and simulation standardisation, verification, and validation. Here, we strive towards a solution by providing an open-access, characterised, constitutive relation (CR) for 3D musculotendon models.

METHODS

The musculotendon complex is modelled following the state-of-the-art active stress approach and is treated as hyperelastic, transversely isotropic, and nearly incompressible. Furthermore, force-length and -velocity relationships are incorporated, and muscle activation is derived from motor-unit information. The CR was implemented within the commercial finite-element software package Abaqus as a user-subroutine. A masticatory system model with left and right masseters was used to demonstrate active and passive movement.

RESULTS

The CR was characterised by various experimental data sets and was able to capture a wide variety of passive and active behaviours. Furthermore, the masticatory simulations revealed that joint movement was sensitive to the muscle's in-fibre passive response.

CONCLUSIONS

This user-material provides a "plug and play" template for 3D neuro-musculoskeletal finite element modelling. We hope that this reduces modelling effort, fosters exchange, and contributes to the standardisation of such models.

How velocity impacts eccentric force generation of fully activated skinned skeletal muscle fibers in long stretches

Eccentric muscle contractions are fundamental to everyday life. They occur markedly in jumping, running, and accidents. Following an initial force rise, stretching of a fully activated muscle can result in a phase of decreasing force ('Give') followed by force redevelopment. However, how the stretch velocity affects 'Give' and force redevelopment remains largely unknown. We investigated the force produced by fully activated single skinned fibers of rat extensor digitorum longus muscles during long stretches. Fibers were pulled from length .85 to 1.3 optimal fiber length at a rate of 1, 10 and 100% of the estimated maximum shortening velocity. 'Give' was absent in slow stretches. Medium and fast stretches yielded a clear 'Give'. After the initial force peak, forces decreased by 11.2% and 27.8% relative to the initial peak force before rising again. During the last half of the stretch (from 1.07 to 1.3 optimal fiber length, which is within the range of the expected descending limb of the force-length relationship), the linear force slope tripled from slow to medium stretch and increased further by 60% from medium to fast stretch. These results are compatible with forcible cross-bridge detachment and re-development of a cross-bridge distribution, and a viscoelastic titin contribution to fiber force. Accounting for these results can improve muscle models and predictions of multi-body simulations.

A Physiology-Guided Classification of Active-Stress and Active-Strain Approaches for Continuum-Mechanical Modeling of Skeletal Muscle Tissue

The well-established sliding filament and cross-bridge theory explain the major biophysical mechanism responsible for a skeletal muscle's active behavior on a cellular level. However, the biomechanical function of skeletal muscles on the tissue scale, which is caused by the complex interplay of muscle fibers and extracellular connective tissue, is much less understood. Mathematical models provide one possibility to investigate physiological hypotheses. Continuum-mechanical models have hereby proven themselves to be very suitable to study the biomechanical behavior of whole muscles or entire limbs. Existing continuum-mechanical skeletal muscle models use either an active-stress or an active-strain approach to phenomenologically describe the mechanical behavior of active contractions. While any macroscopic constitutive model can be judged by it's ability to accurately replicate experimental data, the evaluation of muscle-specific material descriptions is difficult as suitable data is, unfortunately, currently not available. Thus, the discussions become more philosophical rather than following rigid methodological criteria. Within this work, we provide a extensive discussion on the underlying modeling assumptions of both the active-stress and the active-strain approach in the context of existing hypotheses of skeletal muscle physiology. We conclude that the active-stress approach resolves an idealized tissue transmitting active stresses through an independent pathway. In contrast, the active-strain approach reflects an idealized tissue employing an indirect, coupled pathway for active stress transmission. Finally the physiological hypothesis that skeletal muscles exhibit redundant pathways of intramuscular stress transmission represents the basis for considering a mixed-active-stress-active-strain constitutive framework.

Architectural model for muscle growth during maturation

Muscle architecture, which includes parameters like fascicle length, pennation angle, and physiological cross-sectional area, strongly influences skeletal muscles' mechanical properties. During maturation, the muscle architecture has to adapt to a growing organism. This study aimed to develop an architectural model capable of predicting the complete 3D fascicle architecture for primarily unipennate muscles of an arbitrary age, based on fascicle data for an initial age. For model development, we collected novel data on 3D muscle architecture of the rabbit (Oryctolagus cuniculus) M. plantaris of eight animals ranging in age from 29 to 106 days. Experimental results show that plantaris muscle belly length increases by 73%, whereas mean fascicle length and mean pennation angle increases by 39 and 14%, respectively. Those changes were incorporated into the model. In addition to the data collected for M. plantaris the predictions of the model were compared to existing literature data of rabbit M. soleus and M. gastrocnemius medialis. With an error of −1.0 ± 8.6% for relative differences in aponeurosis length, aponeurosis width, muscle height, and muscle mass, the model delivered good results matching interindividual differences. For future studies, the model could be utilized to generate realistic architectural data sets for simulation studies.

How muscle stiffness affects human body model behavior

BACKGROUND

Active human body models (AHBM) consider musculoskeletal movement and joint stiffness via active muscle truss elements in the finite element (FE) codes in dynamic application. In the latest models, such as THUMS™ Version 5, nearly all human muscle groups are modeled in form of one-dimensional truss elements connecting each joint. While a lot of work has been done to improve the active and passive behavior of this 1D muscle system in the past, the volumetric muscle system of THUMS was modeled in a much more simplified way based on Post Mortem Human Subject (PMHS) test data. The stiffness changing effect of isometric contraction was hardly considered for the volumetric muscle system of whole human body models so far. While previous works considered this aspect for single muscles, the effect of a change in stiffness due to isometric contraction of volumetric muscles on the AHBM behavior and computation time is yet unknown.

METHODS

In this study, a simplified frontal impact using the THUMS Version 5 AM50 occupant model was simulated. Key parameters to regulate muscle tissue stiffness of solid elements in THUMS were identified for the material model MAT_SIMPLIFIED_FOAM and different stiffness states were predefined for the buttock and thigh.

RESULTS

During frontal crash, changes in muscle stiffness had an effect on the overall AHBM behavior including expected injury outcome. Changes in muscle stiffness for the thigh and pelvis, as well as for the entire human body model and for strain-rate-dependent stiffness definitions based on literature data had no significant effect on the computation time.

DISCUSSION

Kinematics, peak impact force and stiffness changes were in general compliance with the literature data. However, different experimental setups had to be considered for comparison, as this topic has not been fully investigated experimentally in automotive applications in the past. Therefore, this study has limitations regarding validation of the frontal impact results.

CONCLUSION

Variations of default THUMS material model parameters allow an efficient change in stiffness of volumetric muscles for whole AHBM applications. The computation time is unaffected by altering muscle stiffness using the method suggested in this work. Due to a lack of validation data, the results of this work can only be validated with certain limitations. In future works, the default material models of THUMS could be replaced with recently published models to achieve a possibly more biofidelic muscle behavior, which would even allow a functional dependency of the 1D and 3D muscle systems. However, the effect on calculation time and model stability of these models is yet unknown and should be considered in future studies for efficient AHBM applications.

Power Amplification Increases With Contraction Velocity During Stretch-Shortening Cycles of Skinned Muscle Fibers

Muscle force, work, and power output during concentric contractions (active muscle shortening) are increased immediately following an eccentric contraction (active muscle lengthening). This increase in performance is known as the stretch-shortening cycle (SSC)-effect. Recent findings demonstrate that the SSC-effect is present in the sarcomere itself. More recently, it has been suggested that cross-bridge (XB) kinetics and non-cross-bridge (non-XB) structures (e.g., titin and nebulin) contribute to the SSC-effect. As XBs and non-XB structures are characterized by a velocity dependence, we investigated the impact of stretch-shortening velocity on the SSC-effect. Accordingly, we performed in vitro isovelocity ramp experiments with varying ramp velocities (30, 60, and 85% of maximum contraction velocity for both stretch and shortening) and constant stretch-shortening magnitudes (17% of the optimum sarcomere length) using single skinned fibers of rat soleus muscles. The different contributions of XB and non-XB structures to force production were identified using the XB-inhibitor Blebbistatin. We show that (i) the SSC-effect is velocity-dependent-since the power output increases with increasing SSC-velocity. (ii) The energy recovery (ratio of elastic energy storage and release in the SSC) is higher in the Blebbistatin condition compared with the control condition. The stored and released energy in the Blebbistatin condition can be explained by the viscoelastic properties of the non-XB structure titin. Consequently, our experimental findings suggest that the energy stored in titin during the eccentric phase contributes to the SSC-effect in a velocity-dependent manner.

Variations in Muscle Activity and Exerted Torque During Temporary Blood Flow Restriction in Healthy Individuals

Recent studies suggest that transitory blood flow restriction (BFR) may improve the outcomes of training from anatomical (hypertrophy) and neural control perspectives. Whilst the chronic consequences of BFR on local metabolism and tissue adaptation have been extensively investigated, its acute effects on motor control are not yet fully understood. In this study, we compared the neuromechanical effects of continuous BFR against non-restricted circulation (atmospheric pressure—AP), during isometric elbow flexions. BFR was achieved applying external pressure either between systolic and diastolic (lower pressure—LP) or 1.3 times the systolic pressure (higher pressure—HP). Three levels of torque (15, 30, and 50% of the maximal voluntary contraction—MVC) were combined with the three levels of pressure for a total of 9 (randomized) test cases. Each condition was repeated 3 times. The protocol was administered to 12 healthy young adults. Neuromechanical measurements (torque and high-density electromyography—HDEMG) and reported discomfort were used to investigate the response of the central nervous system to BFR. The investigated variables were: root mean square (RMS), and area under the curve in the frequency domain—for the torque, and average RMS, median frequency and average muscle fibres conduction velocity—for the EMG. The discomfort caused by BFR was exacerbated by the level of torque and accumulated over time. The torque RMS value did not change across conditions and repetitions. Its spectral content, however, revealed a decrease in power at the tremor band (alpha-band, 5–15 Hz) which was enhanced by the level of pressure and the repetition number. The EMG amplitude showed no differences whilst the median frequency and the conduction velocity decreased over time and across trials, but only for the highest levels of torque and pressure. Taken together, our results show strong yet transitory effects of BFR that are compatible with a motor neuron pool inhibition caused by increased activity of type III and IV afferences, and a decreased activity of spindle afferents. We speculate that a compensation of the central drive may be necessary to maintain the mechanical output unchanged, despite disturbances in the afferent volley to the motor neuron pool.

A Systematic Review and Meta-Analysis on the Longitudinal Effects of Unilateral Knee Extension Exercise on Muscle Strength

The aim of the study was to investigate the time-dependent increase in the knee extensors' isometric strength as a response to voluntary, unilateral, isometric knee extension exercise (UIKEE). To do so, a systematic review was carried out to obtain data for a Bayesian longitudinal model-based meta-analysis (BLMBMA). For the systematic review, PubMed, Web of Science, SCOPUS, Chochrane Library were used as databases. The systematic review included only studies that reported on healthy, young individuals performing UIKEE. Studies utilizing a bilateral training protocol were excluded as the focus of this review lied on unilateral training. Out of the 3,870 studies, which were reviewed, 20 studies fulfilled the selected inclusion criteria. These 20 studies were included in the BLMBMA to investigate the time-dependent effects of UIKEE. If compared to the baseline strength of the trained limb, these data reveal that UKIEE can increase the isometric strength by up to 46%. A meta-analysis based on the last time-point of each available study was employed to support further investigations into UIKEE-induced strength increase. A sensitivity analysis showed that intensity of training (%MVC), fraction of male subjects and the average age of the subject had no significant influence on the strength gain. Convergence of BLMBMA revealed that the peak strength increase is reached after ~4 weeks of UIKEE training.

Surface-based modeling of muscles: Functional simulation of the shoulder

Musculoskeletal simulations are an essential tool for studying functional implications of pathologies and of potential surgical outcomes, e.g., for the complex shoulder anatomy. Most shoulder models rely on line-segment approximation of muscles with potential limitations. Comprehensive shoulder models based on continuum-mechanics are scarce due to their complexity in both modeling and computation. In this paper, we present a surface-based modeling approach for muscles, which simplifies the modeling process and is efficient for computation. We propose to use surface geometries for modeling muscles, and devise an automatic approach to generate such models, given the locations of the origin and insertion of tendons. The surfaces are expressed as higher-order tensor B-splines, which ensure smoothness of the geometrical representation. They are simulated as membrane elements within a finite element simulation. This is demonstrated on a comprehensive model of the upper limb, where muscle activations needed to perform desired motions are obtained by using inverse dynamics. In synthetic examples, we demonstrate our proposed surface elements both to be easy to customize (e.g., with spatially varying material properties) and to be substantially (up to 12 times) faster in simulation compared to their volumetric counterpart. With our presented automatic approach of muscle wrapping around bones, the humeral head is exemplified to be wrapped physiologically consistently with surface elements. Our functional simulation is shown to successfully replicate a tracked shoulder motion during activities of daily living. We demonstrate surface-based models to be a numerically stable and computationally efficient compromise between line-segment and volumetric models, enabling anatomical correctness, subject-specific customization, and fast simulations, for a comprehensive simulation of musculoskeletal motion.

Predicting Perturbed Human Arm Movements in a Neuro-Musculoskeletal Model to Investigate the Muscular Force Response

Human movement is generated by a dynamic interplay between the nervous system, the biomechanical structures, and the environment. To investigate this interaction, we propose a neuro-musculoskeletal model of human goal-directed arm movements. Using this model, we simulated static perturbations of the inertia and damping properties of the arm, as well as dynamic torque perturbations for one-degree-of freedom movements around the elbow joint. The controller consists of a feed-forward motor command and feedback based on muscle fiber length and contraction velocity representing short-latency (25 ms) or long-latency (50 ms) stretch reflexes as the first neuronal responses elicited by an external perturbation. To determine the open-loop control signal, we parameterized the control signal resulting in a piecewise constant stimulation over time for each muscle. Interestingly, such an intermittent open-loop signal results in a smooth movement that is close to experimental observations. So, our model can generate the unperturbed point-to-point movement solely by the feed-forward command. The feedback only contributed to the stimulation in perturbed movements. We found that the relative contribution of this feedback is small compared to the feed-forward control and that the characteristics of the musculoskeletal system create an immediate and beneficial reaction to the investigated perturbations. The novelty of these findings is (1) the reproduction of static as well as dynamic perturbation experiments in one neuro-musculoskeletal model with only one set of basic parameters. This allows to investigate the model's neuro-muscular response to the perturbations that-at least to some degree-represent stereotypical interactions with the environment; (2) the demonstration that in feed-forward driven movements the muscle characteristics generate a mechanical response with zero-time delay which helps to compensate for the perturbations; (3) that this model provides enough biomechanical detail to allow for the prediction of internal forces, including joint loads and muscle-bone contact forces which are relevant in ergonomics and for the development of assistive devices but cannot be observed in experiments.

Occlusal load modelling significantly impacts the predicted tooth stress response during biting: a simulation study

Computational models of the masticatory system can provide estimates of occlusal loading during (static) biting or (dynamic) chewing and therefore can be used to evaluate and optimize functional performance of prosthodontic devices and guide dental surgery planning. The modelling assumptions, however, need to be chosen carefully in order to obtain meaningful predictions. The objectives of this study were two-fold: (i) develop a computational model to calculate the stress response of the first molar during biting of a rubber sample and (ii) evaluate the influence of different occlusal load models on the stress response of dental structures. A three-dimensional finite element model was developed comprising the mandible, first molar, associated dental structures, and the articular fossa and discs. Simulations of a maximum force bite on a rubber sample were performed by applying muscle forces as boundary conditions on the mandible and computing the contact between the rubber and molars (GS case). The molar occlusal force was then modelled as a single point force (CF1 case), four point forces (CF2 case), and as a sphere compressing against the occlusal surface (SL case). The peak enamel stress for the GS case was 110 MPa and 677 MPa, 270 MPa and 305 MPa for the CF1, CF2 and SL cases, respectively. Peak dentin stress for the GS case was 44 MPa and 46 MPa, 50 MPa and 63 MPa for the CF1, CF2 and SL cases, respectively. Furthermore, the enamel stress distribution was also strongly correlated to the occlusal load model. The way in which occlusal load is modelled has a substantial influence on the stress response of enamel during biting, but has relatively little impact on the behavior of dentin. The use of point forces or sphere contact to model occlusal loading during mastication overestimates enamel stress magnitude and also influences enamel stress distribution.

Estimating Biomechanical Time-Series with Wearable Sensors: A Systematic Review of Machine Learning Techniques

Wearable sensors have the potential to enable comprehensive patient characterization and optimized clinical intervention. Critical to realizing this vision is accurate estimation of biomechanical time-series in daily-life, including joint, segment, and muscle kinetics and kinematics, from wearable sensor data. The use of physical models for estimation of these quantities often requires many wearable devices making practical implementation more difficult. However, regression techniques may provide a viable alternative by allowing the use of a reduced number of sensors for estimating biomechanical time-series. Herein, we review 46 articles that used regression algorithms to estimate joint, segment, and muscle kinematics and kinetics. We present a high-level comparison of the many different techniques identified and discuss the implications of our findings concerning practical implementation and further improving estimation accuracy. In particular, we found that several studies report the incorporation of domain knowledge often yielded superior performance. Further, most models were trained on small datasets in which case nonparametric regression often performed best. No models were open-sourced, and most were subject-specific and not validated on impaired populations. Future research should focus on developing open-source algorithms using complementary physics-based and machine learning techniques that are validated in clinically impaired populations. This approach may further improve estimation performance and reduce barriers to clinical adoption.

Characterization of Electromechanical Delay Based on a Biophysical Multi-Scale Skeletal Muscle Model

Skeletal muscles can be voluntary controlled by the somatic nervous system yielding an active contractile stress response. Thereby, the active muscle stresses are transmitted to the skeleton by a cascade of connective tissue and thus enable motion. In the context of joint perturbations as well as the assessment of the complexity of neural control, the initial phase of the muscle-tendon system's stress response has a particular importance and is analyzed by means of electromechanical delay (EMD). EMD is defined as the time lag between the stimulation of a muscle and a measurable change in force output. While EMD is believed to depend on multiple structures / phenomena, it is hard to separate their contributions experimentally. We employ a physiologically detailed, three-dimensional, multi-scale model of an idealized muscle-tendon system to analyze the influence of (i) muscle and tendon length, (ii) the material behavior of skeletal muscle and tendon tissue, (iii) the chemo-electro-mechanical behavior of the muscle fibers and (iv) neural control on EMD. Comparisons with experimental data show that simulated EMD values are within the physiological range, i.e., between 6.1 and 68.6 ms, and that the model is able to reproduce the characteristic EMD-stretch curve, yielding the minimum EMD at optimal length. Simulating consecutive recruitment of motor units increases EMD by more than 20 ms, indicating that during voluntary contractions neural control is the dominant factor determining EMD. In contrast, the muscle fiber action potential conduction velocity is found to influence EMD even of a 27 cm long muscle by not more than 3.7 ms. We further demonstrate that in conditions where only little pre-stretch is applied to a muscle-tendon system, the mechanical behavior of both muscle and tendon tissue considerably impacts EMD. Predicting EMD for different muscle and tendon lengths indicates that the anatomy of a specific muscle-tendon system is optimized for its function, i.e., shorter tendon lengths are beneficial to minimize the neural control effort for muscles primary acting as motor in concentric contractions.

Multiscale modeling of the neuromuscular system: Coupling neurophysiology and skeletal muscle mechanics

Mathematical models and computer simulations have the great potential to substantially increase our understanding of the biophysical behavior of the neuromuscular system. This, however, requires detailed multiscale, and multiphysics models. Once validated, such models allow systematic in silico investigations that are not necessarily feasible within experiments and, therefore, have the ability to provide valuable insights into the complex interrelations within the healthy system and for pathological conditions. Most of the existing models focus on individual parts of the neuromuscular system and do not consider the neuromuscular system as an integrated physiological system. Hence, the aim of this advanced review is to facilitate the prospective development of detailed biophysical models of the entire neuromuscular system. For this purpose, this review is subdivided into three parts. The first part introduces the key anatomical and physiological aspects of the healthy neuromuscular system necessary for modeling the neuromuscular system. The second part provides an overview on state-of-the-art modeling approaches representing all major components of the neuromuscular system on different time and length scales. Within the last part, a specific multiscale neuromuscular system model is introduced. The integrated system model combines existing models of the motor neuron pool, of the sensory system and of a multiscale model describing the mechanical behavior of skeletal muscles. Since many sub-models are based on strictly biophysical modeling approaches, it closely represents the underlying physiological system and thus could be employed as starting point for further improvements and future developments. This article is categorized under: Physiology > Mammalian Physiology in Health and Disease Analytical and Computational Methods > Computational Methods Models of Systems Properties and Processes > Organ, Tissue, and Physiological Models.

A comprehensive and volumetric musculoskeletal model for the dynamic simulation of the shoulder function

We present a volumetric and extensive finite element model of the shoulder usable in the context of inverse control, in which the scapula is left unconstrained on the ribcage. Such a model allows for exploring various shoulder movements, which are essential for making patient-specific decisions. The proposed model consists of 23 volumetric muscles parts modelled using the finite element method. The glenohumeral, acromioclavicular and sternoclavicular joints are modelled with soft ball-socket constraints. The musculoskeletal model can be controlled by a tracking-based algorithm, finding the excitations values in the muscles needed to follow some target points. The moment arms obtained during abduction and rotation are compared with the literature, which includes results from cadaveric data and a fine FE model of the rotator cuff and the deltoid. We simulated the paralysis of serratus anterior, a main reason of scapular winging, and compared it with its physiological counterpart. A deficiency in the range of motion as well as a reduction in upward rotation were observed, which both corroborate clinical observations. This is one of the most comprehensive model of the shoulder, which can be used to study complex pathologies of the shoulder and their impact on functional outcome such as range-of-motion.

A Systematic Review of Continuum Modeling of Skeletal Muscles: Current Trends, Limitations, and Recommendations

Finite elasticity theory has been commonly used to model skeletal muscle. A very large range of heterogeneous constitutive laws has been proposed. In this review, the most widely used continuum models of skeletal muscles were synthetized and discussed. Trends and limitations of these laws were highlighted to propose new recommendations for future researches. A systematic review process was performed using two reliable search engines as PubMed and ScienceDirect. 40 representative studies (13 passive muscle materials and 27 active muscle materials) were included into this review. Note that exclusion criteria include tendon models, analytical models, 1D geometrical models, supplement papers, and indexed conference papers. Trends of current skeletal muscle modeling relate to 3D accurate muscle representation, parameter identification in passive muscle modeling, and the integration of coupled biophysical phenomena. Parameter identification for active materials, assumed fiber distribution, data assumption, and model validation are current drawbacks. New recommendations deal with the incorporation of multimodal data derived from medical imaging, the integration of more biophysical phenomena, and model reproducibility. Accounting for data uncertainty in skeletal muscle modeling will be also a challenging issue. This review provides, for the first time, a holistic view of current continuum models of skeletal muscles to identify potential gaps of current models according to the physiology of skeletal muscle. This opens new avenues for improving skeletal muscle modeling in the framework of in silico medicine.

An Efficient Modelling-Simulation-Analysis Workflow to Investigate Stump-Socket Interaction Using Patient-Specific, Three-Dimensional, Continuum-Mechanical, Finite Element Residual Limb Models

The lack of an efficient modelling-simulation-analysis workflow for creating and utilising detailed subject-specific computational models is one of the key reasons why simulation-based approaches for analysing socket-stump interaction have not yet been successfully established. Herein, we propose a novel and efficient modelling-simulation-analysis workflow that uses commercial software for generating a detailed subject-specific, three-dimensional finite element model of an entire residual limb from Diffusion Tensor MRI images in <20 min. Moreover, to complete the modelling-simulation-analysis workflow, the generated subject-specific residual limb model is used within an implicit dynamic FE simulation of bipedal stance to predict the potential sites of deep tissue injury. For this purpose, a nonlinear hyperelastic, transversely isotropic skeletal muscle constitutive law containing a deep tissue injury model was implemented in LS-DYNA. To demonstrate the feasibility of the entire modelling-simulation-analysis workflow and the fact that detailed, anatomically realistic, multi-muscle models are superior to state-of-the-art, fused-muscle models, an implicit dynamic FE analysis of 2-h bipedal stance is carried out. By analysing the potential volume of damaged muscle tissue after donning an optimally-fitted and a misfitted socket, i.e., a socket whose volume was isotropically shrunk by 10%, we were able to highlight the differences between the detailed individual- and fused-muscle models. The results of the bipedal stance simulation showed that peak stresses in the fused-muscle model were four times lower when compared to the multi-muscle model. The peak interface stress in the individual-muscle model, at the end of bipedal stance analysis, was 2.63 times lower than that in the deep tissues of the stump. At the end of the bipedal stance analysis using the misfitted socket, the fused-muscle model predicted that 7.65% of the residual limb volume was injured, while the detailed-model predicted 16.03%. The proposed approach is not only limited to modelling residual limbs but also has applications in predicting the impact of plastic surgery, for detailed forward-dynamics simulations of normal musculoskeletal systems.

Enabling Detailed, Biophysics-Based Skeletal Muscle Models on HPC Systems

Realistic simulations of detailed, biophysics-based, multi-scale models often require very high resolution and, thus, large-scale compute facilities. Existing simulation environments, especially for biomedical applications, are typically designed to allow for high flexibility and generality in model development. Flexibility and model development, however, are often a limiting factor for large-scale simulations. Therefore, new models are typically tested and run on small-scale compute facilities. By using a detailed biophysics-based, chemo-electromechanical skeletal muscle model and the international open-source software library OpenCMISS as an example, we present an approach to upgrade an existing muscle simulation framework from a moderately parallel version toward a massively parallel one that scales both in terms of problem size and in terms of the number of parallel processes. For this purpose, we investigate different modeling, algorithmic and implementational aspects. We present improvements addressing both numerical and parallel scalability. In addition, our approach includes a novel visualization environment which is based on the MegaMol framework and is capable of handling large amounts of simulated data. We present the results of a number of scaling studies at the Tier-1 supercomputer HazelHen at the High Performance Computing Center Stuttgart (HLRS). We improve the overall runtime by a factor of up to 2.6 and achieve good scalability on up to 768 cores.

Gradient-based optimization with B-splines on sparse grids for solving forward-dynamics simulations of three-dimensional, continuum-mechanical musculoskeletal system models

Investigating the interplay between muscular activity and motion is the basis to improve our understanding of healthy or diseased musculoskeletal systems. To be able to analyze the musculoskeletal systems, computational models are used. Albeit some severe modeling assumptions, almost all existing musculoskeletal system simulations appeal to multibody simulation frameworks. Although continuum-mechanical musculoskeletal system models can compensate for some of these limitations, they are essentially not considered because of their computational complexity and cost. The proposed framework is the first activation-driven musculoskeletal system model, in which the exerted skeletal muscle forces are computed using 3-dimensional, continuum-mechanical skeletal muscle models and in which muscle activations are determined based on a constraint optimization problem. Numerical feasibility is achieved by computing sparse grid surrogates with hierarchical B-splines, and adaptive sparse grid refinement further reduces the computational effort. The choice of B-splines allows the use of all existing gradient-based optimization techniques without further numerical approximation. This paper demonstrates that the resulting surrogates have low relative errors (less than 0.76%) and can be used within forward simulations that are subject to constraint optimization. To demonstrate this, we set up several different test scenarios in which an upper limb model consisting of the elbow joint, the biceps and triceps brachii, and an external load is subjected to different optimization criteria. Even though this novel method has only been demonstrated for a 2-muscle system, it can easily be extended to musculoskeletal systems with 3 or more muscles.

Intramuscular Pressure of Tibialis Anterior Reflects Ankle Torque but Does Not Follow Joint Angle-Torque Relationship

Intramuscular pressure (IMP) is the hydrostatic fluid pressure that is directly related to muscle force production. Electromechanical delay (EMD) provides a link between mechanical and electrophysiological quantities and IMP has potential to detect local electromechanical changes. The goal of this study was to assess the relationship of IMP with the mechanical and electrical characteristics of the tibialis anterior muscle (TA) activity at different ankle positions. We hypothesized that (1) the TA IMP and the surface EMG (sEMG) and fine-wire EMG (fwEMG) correlate to ankle joint torque, (2) the isometric force of TA increases at increased muscle lengths, which were imposed by a change in ankle angle and IMP follows the length-tension relationship characteristics, and (3) the electromechanical delay (EMD) is greater than the EMD of IMP during isometric contractions. Fourteen healthy adults [7 female; mean (SD) age = 26.9 (4.2) years old with 25.9 (5.5) kg/m² body mass index] performed (i) three isometric dorsiflexion (DF) maximum voluntary contraction (MVC) and (ii) three isometric DF ramp contractions from 0 to 80% MVC at rate of 15% MVC/second at DF, Neutral, and plantarflexion (PF) positions. Ankle torque, IMP, TA fwEMG, and TA sEMG were measured simultaneously. The IMP, fwEMG, and sEMG were significantly correlated to the ankle torque during ramp contractions at each ankle position tested. This suggests that IMP captures in vivo mechanical properties of active muscles. The ankle torque changed significantly at different ankle positions however, the IMP did not reflect the change. This is explained with the opposing effects of higher compartmental pressure at DF in contrast to the increased force at PF position. Additionally, the onset of IMP activity is found to be significantly earlier than the onset of force which indicates that IMP can be designed to detect muscular changes in the course of neuromuscular diseases impairing electromechanical transmission.

Elucidating Sensorimotor Control Principles with Myoelectric Musculoskeletal Models

There is an old saying that you must walk a mile in someone's shoes to truly understand them. This mini-review will synthesize and discuss recent research that attempts to make humans "walk a mile" in an artificial musculoskeletal system to gain insight into the principles governing human movement control. In this approach, electromyography (EMG) is used to sample human motor commands; these commands serve as inputs to mathematical models of muscular dynamics, which in turn act on a model of skeletal dynamics to produce a simulated motor action in real-time (i.e., the model's state is updated fast enough produce smooth motion without noticeable transitions; Manal et al., 2002). In this mini-review, these are termed myoelectric musculoskeletal models (MMMs). After a brief overview of typical MMM design and operation principles, the review will highlight how MMMs have been used for understanding human sensorimotor control and learning by evoking apparent alterations in a user's biomechanics, neural control, and sensory feedback experiences.