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

Clinical treatments of skeletal muscle weakness are hindered by a lack of an approach to evaluate individual muscle force. Intramuscular pressure (IMP) has shown a correlation to muscle force in vivo, but patient to patient and muscle to muscle variability results in difficulty of utilizing IMP to estimate muscle force. The goal of this work was to develop a finite element model of whole skeletal muscle that can predict IMP under passive and active conditions to further investigate the mechanisms of IMP variability. A previously validated hypervisco-poroelastic constitutive approach was modified to incorporate muscle activation through an inhomogeneous geometry. Model parameters were optimized to fit model stress to experimental data, and the resulting model fluid pressurization data were utilized for validation. Model fitting was excellent (root-mean-square error or RMSE <1.5 kPa for passive and active conditions), and IMP predictive capability was strong for both passive (RMSE 3.5 mmHg) and active (RMSE 10 mmHg at in vivo lengths) conditions. Additionally, model fluid pressure was affected by length under isometric conditions, as increases in stretch yielded decreases in fluid pressurization following a contraction, resulting from counteracting Poisson effects. Model pressure also varied spatially, with the highest gradients located near aponeuroses. These findings may explain variability of in vivo IMP measurements in the clinic, and thus help reduce this variability in future studies. Further development of this model to include isotonic contractions and muscle weakness would greatly benefit this work.