Modeling Skeletal Muscle Stress and Intramuscular Pressure: A Whole Muscle Active-Passive Approach

Effects of static exercises on hip muscle fatigue and knee wobble assessed by surface electromyography and inertial measurement unit data

Hip muscle weakness can be a precursor to or a result of lower limb injuries. Assessment of hip muscle strength and muscle motor fatigue in the clinic is important for diagnosing and treating hip-related impairments. Muscle motor fatigue can be assessed with surface electromyography (sEMG), however sEMG requires specialized equipment and training. Inertial measurement units (IMUs) are wearable devices used to measure human motion, yet it remains unclear if they can be used as a low-cost alternative method to measure hip muscle fatigue. The goals of this work were to (1) identify which of five pre-selected exercises most consistently and effectively elicited muscle fatigue in the gluteus maximus, gluteus medius, and rectus femoris muscles and (2) determine the relationship between muscle fatigue using sEMG sensors and knee wobble using an IMU device. This work suggests that a wall sit and single leg knee raise activity fatigue the gluteus medius, gluteus maximus, and rectus femoris muscles most reliably (p < 0.05) and that the gluteus medius and gluteus maximus muscles were fatigued to a greater extent than the rectus femoris (p = 0.031 and p = 0.0023, respectively). Additionally, while acceleration data from a single IMU placed on the knee suggested that more knee wobble may be an indicator of muscle fatigue, this single IMU is not capable of reliably assessing fatigue level. These results suggest the wall sit activity could be used as simple, static exercise to elicit hip muscle fatigue in the clinic, and that assessment of knee wobble in addition to other IMU measures could potentially be used to infer muscle fatigue under controlled conditions. Future work examining the relationship between IMU data, muscle fatigue, and multi-limb dynamics should be explored to develop an accessible, low-cost, fast and standardized method to measure fatiguability of the hip muscles in the clinic.

Male and female syringeal muscles exhibit superfast shortening velocities in zebra finches

Vocalisations play a key role in the communication behaviour of many vertebrates. Vocal production requires extremely precise motor control, which is executed by superfast vocal muscles that can operate at cycle frequencies over 100 Hz and up to 250 Hz. The mechanical performance of these muscles has been quantified with isometric performance and the workloop technique, but owing to methodological limitations we lack a key muscle property characterising muscle performance, the force-velocity relationship. Here, we quantified the force-velocity relationship in zebra finch superfast syringeal muscles using the isovelocity technique and tested whether the maximal shortening velocity is different between males and females. We show that syringeal muscles exhibit high maximal shortening velocities of 25L0 s-1 at 30°C. Using Q10-based extrapolation, we estimate they can reach 37-42L0 s-1 on average at body temperature, exceeding other vocal and non-avian skeletal muscles. The increased speed does not adequately compensate for reduced force, which results in low power output. This further highlights the importance of high-frequency operation in these muscles. Furthermore, we show that isometric properties positively correlate with maximal shortening velocities. Although male and female muscles differ in isometric force development rates, maximal shortening velocity is not sex dependent. We also show that cyclical methods to measure force-length properties used in laryngeal studies give the same result as conventional stepwise methodologies, suggesting either approach is appropriate. We argue that vocal behaviour may be affected by the high thermal dependence of superfast vocal muscle performance.

Three-dimensional mapping of ultrasound-derived skeletal muscle shear wave velocity

Introduction: The mechanical properties of skeletal muscle are indicative of its capacity to perform physical work, state of disease, or risk of injury. Ultrasound shear wave elastography conducts a quantitative analysis of a tissue's shear stiffness, but current implementations only provide two-dimensional measurements with limited spatial extent. We propose and assess a framework to overcome this inherent limitation by acquiring numerous and contiguous measurements while tracking the probe position to create a volumetric scan of the muscle. This volume reconstruction is then mapped into a parameterized representation in reference to geometric and anatomical properties of the muscle. Such an approach allows to quantify regional differences in muscle stiffness to be identified across the entire muscle volume assessed, which could be linked to functional implications. Methods: We performed shear wave elastography measurements on the vastus lateralis (VL) and the biceps femoris long head (BFlh) muscle of 16 healthy volunteers. We assessed test-retest reliability, explored the potential of the proposed framework in aggregating measurements of multiple subjects, and studied the acute effects of muscular contraction on the regional shear wave velocity post-measured at rest. Results: The proposed approach yielded moderate to good reliability (ICC between 0.578 and 0.801). Aggregation of multiple subject measurements revealed considerable but consistent regional variations in shear wave velocity. As a result of muscle contraction, the shear wave velocity was elevated in various regions of the muscle; showing pre-to-post regional differences for the radial assessement of VL and longitudinally for BFlh. Post-contraction shear wave velocity was associated with maximum eccentric hamstring strength produced during six Nordic hamstring exercise repetitions. Discussion and Conclusion: The presented approach provides reliable, spatially resolved representations of skeletal muscle shear wave velocity and is capable of detecting changes in three-dimensional shear wave velocity patterns, such as those induced by muscle contraction. The observed systematic inter-subject variations in shear wave velocity throughout skeletal muscle additionally underline the necessity of accurate spatial referencing of measurements. Short high-effort exercise bouts increase muscle shear wave velocity. Further studies should investigate the potential of shear wave elastography in predicting the muscle's capacity to perform work.

Aponeurosis structure-function properties: Evidence of heterogeneity and implications for muscle function

Aponeurosis is a sheath-like connective tissue that aids in force transmission from muscle to tendon and can be found throughout the musculoskeletal system. The key role of aponeurosis in muscle-tendon unit mechanics is clouded by a lack of understanding of aponeurosis structure-function properties. This work aimed to determine the heterogeneous material properties of porcine triceps brachii aponeurosis tissue with materials testing and evaluate heterogeneous aponeurosis microstructure with scanning electron microscopy. We found that aponeurosis may exhibit more microstructural collagen waviness in the insertion region (near the tendon) compared to the transition region (near the muscle midbelly) (1.20 versus 1.12, p = 0.055), which and a less stiff stress-strain response in the insertion versus transition regions (p < 0.05). We also showed that different assumptions of aponeurosis heterogeneity, specifically variations in elastic modulus with location can alter the stiffness (by more than 10x) and strain (by approximately 10% muscle fiber strain) of a finite element model of muscle and aponeurosis. Collectively, these results suggest that aponeurosis heterogeneity could be due to variations in tissue microstructure and that different approaches to modeling tissue heterogeneity alters the behavior of computational models of muscle-tendon units. STATEMENT OF SIGNIFICANCE: Aponeurosis is a connective tissue found in many muscle tendon units that aids in force transmission, yet little is known about the specific material properties of aponeurosis. This work aimed to determine how the properties of aponeurosis tissue varied with location. We found that aponeurosis exhibits more microstructural waviness near the tendon compared to near the muscle midbelly, which was associated with differences in tissue stiffness. We also showed that different variations in aponeurosis modulus (stiffness) can alter the stiffness and stretch of a computer model of muscle tissue. These results show that assuming uniform aponeurosis structure and modulus, which is common, may lead to inaccurate models of the musculoskeletal system.

Exploring the role of intratendinous pressure in the pathogenesis of tendon pathology: a narrative review and conceptual framework

Despite the high prevalence of tendon pathology in athletes, the underlying pathogenesis is still poorly understood. Various aetiological theories have been presented and rejected in the past, but the tendon cell response model still holds true. This model describes how the tendon cell is the key regulator of the extracellular matrix and how pathology is induced by a failed adaptation to a disturbance of tissue homeostasis. Such failure has been attributed to various kinds of stressors (eg, mechanical, thermal and ischaemic), but crucial elements seem to be missing to fully understand the pathogenesis. Importantly, a disturbance of tissue pressure homeostasis has not yet been considered a possible factor, despite it being associated with numerous pathologies. Therefore, we conducted an extensive narrative literature review on the possible role of intratendinous pressure in the pathogenesis of tendon pathology. This review explores the current understanding of pressure dynamics and the role of tissue pressure in the pathogenesis of other disorders with structural similarities to tendons. By bridging these insights with known structural changes that occur in tendon pathology, a conceptual model was constituted. This model provides an overview of the possible mechanism of how an increase in intratendinous pressure might be involved in the development and progression of tendon pathology and contribute to tendon pain. In addition, some therapies that could reduce intratendinous pressure and accelerate tendon healing are proposed. Further experimental research is encouraged to investigate our hypotheses and to initiate debate on the relevance of intratendinous pressure in tendon pathology.

Modeling of active skeletal muscles: a 3D continuum approach incorporating multiple muscle interactions

Skeletal muscles have a highly organized hierarchical structure, whose main function is to generate forces for movement and stability. To understand the complex heterogeneous behaviors of muscles, computational modeling has advanced as a non-invasive approach to evaluate relevant mechanical quantities. Aiming to improve musculoskeletal predictions, this paper presents a framework for modeling 3D deformable muscles that includes continuum constitutive representation, parametric determination, model validation, fiber distribution estimation, and integration of multiple muscles into a system level for joint motion simulation. The passive and active muscle properties were modeled based on the strain energy approach with Hill-type hyperelastic constitutive laws. A parametric study was conducted to validate the model using experimental datasets of passive and active rabbit leg muscles. The active muscle model with calibrated material parameters was then implemented to simulate knee bending during a squat with multiple quadriceps muscles. A computational fluid dynamics (CFD) fiber simulation approach was utilized to estimate the fiber arrangements for each muscle, and a cohesive contact approach was applied to simulate the interactions among muscles. The single muscle simulation results showed that both passive and active muscle elongation responses matched the range of the testing data. The dynamic simulation of knee flexion and extension showed the predictive capability of the model for estimating the active quadriceps responses, which indicates that the presented modeling pipeline is effective and stable for simulating multiple muscle configurations. This work provided an effective framework of a 3D continuum muscle model for complex muscle behavior simulation, which will facilitate additional computational and experimental studies of skeletal muscle mechanics. This study will offer valuable insight into the future development of multiscale neuromuscular models and applications of these models to a wide variety of relevant areas such as biomechanics and clinical research.

Formulation and exploration of novel, intramuscular pressure based, muscle activation strategies in a spine model

Optimization models are often devised to assess spinal stability via estimating individual muscle forces. However, neglecting muscles' fluidic behavior remains an approximation due to the role of muscle pressure in force transmission. The purpose of this study was to leverage a validated Finite Element (FE) model of the spine, inclusive of Intra-Muscular Pressure (IMP), to explore muscle activation strategies towards maintaining equilibrium spinal stability. Three conventional strategies governing minimizing muscle effort, minimizing IVD compressive forces, and maintaining stability at all costs were first investigated to explore model's validity. Thereafter, two novel IMP-based strategies were devised and explored, specifically minimizing and maximizing IMP. The model was previously shown valid in light of in vivo and in silico observations with an average discrepancy of 6%. This being the case, the conventional strategies dictated efficacy in muscular activations whilst maintaining an equilibrium stable position, as quantified in the present paper, with a difference of 9.8% from documented data. In addition, the explored novel IMP-based strategies suggested the presence of a threshold individual muscles IMP, approximately 272 mmHg for the longissimus muscle for example, beyond which muscles potentially start to share radial loads with surrounding tissues, whilst limiting the contraction of the underlying muscles. In conclusion, this study theoretically supports the possibility of activation strategies based on muscular pressure, which the developed, verified, and validated FE spine model was leveraged to investigate. The explored novel IMP-based strategies may have significance in informing clinical applications such as motion analysis and functional electrical stimulation of muscles.

The effects of gravity and compression on interstitial fluid transport in the lower limb

Edema in the limbs can arise from pathologies such as elevated capillary pressures due to failure of venous valves, elevated capillary permeability from local inflammation, and insufficient fluid clearance by the lymphatic system. The most common treatments include elevation of the limb, compression wraps and manual lymphatic drainage therapy. To better understand these clinical situations, we have developed a comprehensive model of the solid and fluid mechanics of a lower limb that includes the effects of gravity. The local fluid balance in the interstitial space includes a source from the capillaries, a sink due to lymphatic clearance, and movement through the interstitial space due to both gravity and gradients in interstitial fluid pressure (IFP). From dimensional analysis and numerical solutions of the governing equations we have identified several parameter groups that determine the essential length and time scales involved. We find that gravity can have dramatic effects on the fluid balance in the limb with the possibility that a positive feedback loop can develop that facilitates chronic edema. This process involves localized tissue swelling which increases the hydraulic conductivity, thus allowing the movement of interstitial fluid vertically throughout the limb due to gravity and causing further swelling. The presence of a compression wrap can interrupt this feedback loop. We find that only by modeling the complex interplay between the solid and fluid mechanics can we adequately investigate edema development and treatment in a gravity dependent limb.

Systematic review of skeletal muscle passive mechanics experimental methodology

Understanding passive skeletal muscle mechanics is critical in defining structure-function relationships in skeletal muscle and ultimately understanding pathologically impaired muscle. In this systematic review, we performed an exhaustive literature search using PRISMA guidelines to quantify passive muscle mechanical properties, summarized the methods used to create these data, and make recommendations to standardize future studies. We screened over 7500 papers and found 80 papers that met the inclusion criteria. These papers reported passive muscle mechanics from single muscle fiber to whole muscle across 16 species and 54 distinct muscles. We found a wide range of methodological differences in sample selection, preparation, testing, and analysis. The systematic review revealed that passive muscle mechanics is species and scale dependent-specifically within mammals, the passive mechanics increases non-linearly with scale. However, a detailed understanding of passive mechanics is still unclear because the varied methodologies impede comparisons across studies, scales, species, and muscles. Therefore, we recommend the following: smaller scales may be maintained within storage solution prior to testing, when samples are tested statically use 2-3 min of relaxation time, stress normalization at the whole muscle level be to physiologic cross-sectional area, strain normalization be to sarcomere length when possible, and an exponential equation be used to fit the data. Additional studies using these recommendations will allow exploration of the multiscale relationship of passive force within and across species to provide the fundamental knowledge needed to improve our understanding of passive muscle mechanics.

Optimization of traditional Chinese medicine rolling manipulation and pressure attenuation

BACKGROUND

Traditional Chinese Medicine (TCM) massage utilizes mechanical force stimulation, and the amount of mechanical force influences therapeutic outcome. This amount is determined by pressure, frequency, and duration; however, there are no standard definitions for these measures.

METHODS

An orthogonal design was used to evaluate massage efficacy using muscle tension as an index. Pressure (2, 4, 6 kg), duration (5, 10, 15 min), frequency (60, 120, 180 repetitions/min), pain (mild, medium, severe), weight (<60, 60-75, >75 kg), and sex (male, female) were evaluated. Additionally, a porcine model of muscle tension was used to construct pressure-time curves for muscle tissues under static and dynamic pressure.

RESULTS

We identified an interaction among the six massage measures (P<0.05). Of these measures, only two were individually significant: manipulation frequency and patient pain level (P<0.05). Specifically, 120 repetitions/min improved muscle tension significantly more than 60 or 180 repetitions/min (P<0.05), and patients with severe pain had significantly improved muscle tension compared to those with medium or mild pain (P<0.05). In the porcine muscle model, both static and dynamic pressure were attenuated by approximately 12.5% per cm. This attenuation dropped to 10% per cm when the pressure sensor was placed below tissues with different thicknesses instead of being inserted into tissues at different levels.

CONCLUSION

Manipulation frequency and patient pain level were primarily responsible for the therapeutic effects of TCM massage. Mechanistically, pressure was attenuated by nearly 75% at a depth of 2 cm from the muscle surface during TCM massage.

Skeletal muscle: Modeling the mechanical behavior by taking the hierarchical microstructure into account

Skeletal muscles ensure the mobility of mammals and are complex natural fiber-matrix-composites with a hierarchical microstructure. In this work, we analyze the muscle's mechanical behavior on the level of fascicles and muscle fibers. We introduce continuum mechanics hyperelastic material models for the connective tissue endomysium and the embedded muscle fibers. The coupled electrical, chemical and mechanical processes taking place in activated contracting muscle fibers are captured including the temporal change of the activation level and the spatial propagation of the activation potential in fibers. In our model, we investigate the material behavior of fascicle, fiber and endomysium in the fiber direction and examine interactions between muscle fiber and endomysium by considering the temporal and spatial change of muscle fiber activation. In addition, a loading case of normal and shear forces is applied to analyze the fiber lifting force and the lifting height of unipennate muscles with different pennation angles. Moreover, the development of local stresses and strains in fibers and endomysium for different strains are studied. The simulation results allow to identify regions in high risk of damage. Optimal arrangements of unipennate muscle microstructure are found for either very small or very large pennation angles.

Longer electromechanical delay in paretic triceps surae muscles during voluntary isometric plantarflexion torque generation in chronic hemispheric stroke survivors

Electromechanical delay (EMD) is the time delay between the onset of muscle activity and the onset of force/joint torque. This delay appears to be linked to muscular contraction efficiency. However, to our knowledge, limited evidence is available regarding the magnitude of the EMD in stroke-impaired muscles. Accordingly, this study aims to quantify the EMD in both paretic and non-paretic triceps surae muscles of chronic hemispheric stroke survivors, and to investigate whether the EMD is related to voluntary force-generating capacity in this muscle group. Nine male chronic stroke survivors were asked to perform isometric plantarflexion contractions at different force levels and at different ankle joint angles ranging from maximum plantarflexion to maximum dorsiflexion. The surface electromyograms were recorded from triceps surae muscles. The longest EMD among triceps surae muscles was chosen as the EMD for each side. Our results revealed that the EMD in paretic muscles was significantly longer than in non-paretic muscles. Moreover, both paretic and non-paretic muscles showed a negative correlation between the EMD and maximum torque-generating capacity. In addition, there was a strong positive relationship between the EMD and shear wave speed in paretic muscles as well as a negative relationship between the EMD and passive ankle joint range of motion. These findings imply that the EMD may be a useful biomarker, in part, associated with contractile and material properties in stroke-impaired muscles.

Anisotropic and viscoelastic tensile mechanical properties of aponeurosis: Experimentation, modeling, and tissue microstructure

Aponeuroses are stiff sheath-like components of the muscle-tendon unit that play a vital role in force transmission and thus locomotion. There is clear importance of the aponeurosis in musculoskeletal function, but there have been relatively few studies of aponeurosis material properties to date. The goals of this work were to: 1) perform tensile stress-relaxation tests, 2) perform planar biaxial tests, 3) employ computational modeling to the data from 1 to 2, and 4) perform scanning electron microscopy to determine collagen fibril organization for aponeurosis tissue. Viscoelastic modeling and statistical analysis of stress-relaxation data showed that while relaxation rate differed statistically between strain levels (p = 0.044), functionally the relaxation behavior was nearly the same. Biaxial testing and associated modeling highlighted the nonlinear (toe region of ~2-3% strain) and anisotropic (longitudinal direction linear modulus ~50 MPa, transverse ~2.5 MPa) tensile mechanical behavior of aponeurosis tissue. Comparisons of various constitutive formulations showed that a transversely isotropic Ogden approach balanced strong fitting (goodness of fit 0.984) with a limited number of parameters (five), while damage modeling parameters were also provided. Scanning electron microscopy showed a composite structure of highly aligned, partially wavy collagen fibrils with more random collagen cables for aponeurosis microstructure. Future work to expand microstructural analysis and use these data to inform computational modeling would benefit this work and the field.

Investigating Passive Muscle Mechanics With Biaxial Stretch

Introduction

The passive stiffness of skeletal muscle can drastically affect muscle function in vivo, such as the case for fibrotic tissue or patients with cerebral palsy. The two constituents of skeletal muscle that dominate passive stiffness are the intracellular protein titin and the collagenous extracellular matrix (ECM). However, efforts to correlate stiffness and measurements of specific muscle constituents have been mixed, and thus the complete mechanisms for changes to muscle stiffness remain unknown. We hypothesize that biaxial stretch can provide an improved approach to evaluating passive muscle stiffness.

Methods

We performed planar biaxial materials testing of passively stretched skeletal muscle and identified three previously published datasets of uniaxial materials testing. We developed and employed a constitutive model of passive skeletal muscle that includes aligned muscle fibers and dispersed ECM collagen fibers with a bimodal von Mises distribution. Parametric modeling studies and fits to experimental data (both biaxial and previously published) were completed.

Results

Biaxial data exhibited differences in time dependent behavior based on orientation (p < 0.0001), suggesting different mechanisms supporting load in the direction of muscle fibers (longitudinal) and in the perpendicular (transverse) directions. Model parametric studies and fits to experimental data exhibited the robustness of the model (<20% error) and how differences in tissue stiffness may not be observed in uniaxial longitudinal stretch, but are apparent in biaxial stretch.

Conclusion

This work presents novel materials testing data of passively stretched skeletal muscle and use of constitutive modeling and finite element analysis to explore the interaction between stiffness, constituent variability, and applied deformation in passive skeletal muscle. The results highlight the importance of biaxial stretch in evaluating muscle stiffness and in further considering the role of ECM collagen in modulating passive muscle stiffness.

An experimental and computational investigation of the effects of volumetric boundary conditions on the compressive mechanics of passive skeletal muscle

Computational modeling, such as finite element analysis, is employed in a range of biomechanics specialties, including impact biomechanics and surgical planning. These models rely on accurate material properties for skeletal muscle, which comprises roughly 40% of the human body. Due to surrounding tissues, compressed skeletal muscle in vivo likely experiences a semi-confined state. Nearly all previous studies investigating passively compressed muscle at the tissue level have focused on muscle in unconfined compression. The goals of this study were to (1) examine the stiffness and time-dependent material properties of skeletal muscle subjected to both confined and unconfined compression (2) develop a model that captures passive muscle mechanics under both conditions and (3) determine the extent to which different assumptions of volumetric behavior affect model results. Muscle in confined compression exhibited stiffer behavior, agreeing with previous assumptions of near-incompressibility. Stress relaxation was found to be faster under unconfined compression, suggesting there may be different mechanisms that support load these two conditions. Finite element calibration was achieved through nonlinear optimization (normalized root mean square error <6%) and model validation was strong (normalized root mean square error <17%). Comparisons to commonly employed assumptions of bulk behavior showed that a simple one parameter approach does not accurately simulate confined compression. We thus recommend the use of a properly calibrated, nonlinear bulk constitutive model for modeling of skeletal muscle in vivo. Future work to determine mechanisms of passive muscle stiffness would enhance the efforts presented here.