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

In vivo mechanical response of thigh soft tissues under compression: A two-layer model allows an improved representation of the local tissue kinematics

Biomechanical parameters have the potential to be used as physical markers for prevention and diagnosis. Finite Element Analysis (FEA) is a widely used tool to evaluate these parameters in vivo. However, the development of clinically relevant FEA requires personalisation of the geometry, boundary conditions, and constitutive parameters. This contribution focuses on the characterisation of mechanical properties in vivo which remains a significant challenge for the community. The aim of this retrospective study is to evaluate the sensitivity of the computed elastic parameters (shear modulus of fat and muscle tissues) derived by inverse analysis as a function of the geometrical modelling assumption (homogenised monolayer vs bilayer) and the formulation of the cost function. The methodology presented here proposes to extract the experimental force-displacement response for each tissue layer (muscle and fat) and construct the associated Finite Element Model for each volunteer, based on data previously collected in our group (N = 7 volunteers) as reported in (Fougeron et al., 2020). The sensitivity analysis indicates that the choice of the cost function has minimal impact on the topology of the response surface in the parametric space. Each surface displays a valley of parameters that minimises the cost function. The constitutive properties of the thigh (reported as median ± interquartile range) were determined to be (μ=198±322Pa,α=37) for the monolayer and (μmuscle=1675±1127Pa,αmuscle=22±14,μfat=537±1131Pa,αfat=32±7) for the bilayer. A comparison of the homogenised monolayer and bilayer models showed that adding a layer reduces the error on the local force displacement curves, increasing the accuracy of the local kinematics of soft tissues during indentation. This allows for an increased understanding of load transmission in soft tissue. The comparison of the two models in terms of strains indicates that the modelling choice significantly influences the localization of maximal compressive strains. These results support the idea that the biomechanical community should conduct further work to develop reliable methodologies for estimating in vivo strain in soft tissue.

Single and bi-compartment poro-elastic model of perfused biological soft tissues: FEniCSx implementation and tutorial

Soft biological tissues demonstrate strong time-dependent and strain-rate mechanical behavior, arising from their intrinsic visco-elasticity and fluid-solid interactions. The time-dependent mechanical properties of soft tissues influence their physiological functions and are related to several pathological processes. Poro-elastic modeling represents a promising approach because it allows the integration of multiscale/multiphysics data to probe biologically relevant phenomena at a smaller scale and embeds the relevant mechanisms at the larger scale. The implementation of multiphase flow poro-elastic models however is a complex undertaking, requiring extensive knowledge. The open-source software FEniCSx Project provides a novel tool for the automated solution of partial differential equations by the finite element method. This paper aims to provide the required tools to model the mixed formulation of poro-elasticity, from the theory to the implementation, within FEniCSx. Several benchmark cases are studied. A column under confined compression conditions is compared to the Terzaghi analytical solution, using the L2-norm. An implementation of poro-hyper-elasticity is proposed. A bi-compartment column is compared to previously published results (Cast3m implementation). For all cases, accurate results are obtained in terms of a normalized Root Mean Square Error (RMSE). Furthermore, the FEniCSx computation is found three times faster than the legacy FEniCS one. The benefits of parallel computation are also highlighted.

Short-term vancomycin and buffer soaking does not change rabbit achilles tendon tensile material properties

BACKGROUND

Allograft tendons are commonly used during orthopedic surgery to reconstruct tissue that is severely damaged. Soaking the tendon in an antibiotic solution, specifically vancomycin, has been shown to lower the risk of post-operative infections. While some material properties of tendon and ligament after antibiotic soaking have previously been characterized, extensive sub-failure allograft tendon material properties after soaking in antibiotic solutions have not.

METHODS

Forty tendons were dissected from rabbits and soaked in either a phosphate buffered saline (PBS) only solution or vancomycin and PBS solution for five or 30 min. Immediately after soaking, quasi-static tensile experiments were performed in a materials testing system.

FINDINGS

Tissue nominal stress, Lagrange strain, toe-region properties and elastic modulus were characterized. For all forty tendons, the average elastic modulus was found to be 455 ± 37 MPa, the average transition strain (from toe-region to linear elastic region) was 0.0487 ± 0.0035, and the average transition stress was 9.71 ± 0.79 MPa. No statistically significant differences in any of these material properties were found across soaking medium or soaking time.

INTERPRETATION

From these results, we conclude that soaking an allograft tendon in antibiotic solution for up to 30 min prior to implantation does not change the tensile material properties of tendons, supporting current clinical practice.

Role of multi-layer tissue composition of musculoskeletal extremities for prediction of in vivo surface indentation response and layer deformations

Emergent mechanics of musculoskeletal extremities (surface indentation stiffness and tissue deformation characteristics) depend on the underlying composition and mechanics of each soft tissue layer (i.e. skin, fat, and muscle). Limited experimental studies have been performed to explore the layer specific relationships that contribute to the surface indentation response. The goal of this study was to examine through statistical modeling how the soft tissue architecture contributed to the aggregate mechanical surface response across 8 different sites of the upper and lower extremities. A publicly available dataset was used to examine the relationship of soft tissue thickness (fat and muscle) to bulk tissue surface compliance. Models required only initial tissue layer thicknesses, making them usable in the future with only a static ultrasound image. Two physics inspired models (series of linear springs), which allowed reduced statistical representations (combined locations and location specific), were explored to determine the best predictability of surface compliance and later individual layer deformations. When considering the predictability of the experimental surface compliance, the physics inspired combined locations model showed an improvement over the location specific model (percent difference of 25.4 +/- 27.9% and 29.7 +/- 31.8% for the combined locations and location specific models, respectively). While the statistical models presented in this study show that tissue compliance relies on the individual layer thicknesses, it is clear that there are other variables that need to be accounted for to improve the model. In addition, the individual layer deformations of fat and muscle tissues can be predicted reasonably well with the physics inspired models, however additional parameters may improve the robustness of the model outcomes, specifically in regard to capturing subject specificity.

Société de Biomécanique Young Investigator Award 2021: Numerical investigation of the time-dependent stress-strain mechanical behaviour of skeletal muscle tissue in the context of pressure ulcer prevention

BACKGROUND

Pressure-induced tissue strain is one major pathway for Pressure Ulcer development and, especially, Deep Tissue Injury. Biomechanical investigation of the time-dependent stress-strain mechanical behaviour of skeletal muscle tissue is therefore essential. In the literature, a viscoelastic formulation is generally assumed for the experimental characterization of skeletal muscles, with the limitation that the underlying physical mechanisms that give rise to the time dependent stress-strain behaviour are not known. The objective of this study is to explore the capability of poroelasticity to reproduce the apparent viscoelastic behaviour of passive muscle tissue under confined compression.

METHODS

Experimental stress-relaxation response of 31 cylindrical porcine samples tested under fast and slow confined compression by Vaidya and collaborators were used. An axisymmetric Finite Element model was developed in ABAQUS and, for each sample a one-to-one inverse analysis was performed to calibrate the specimen-specific constitutive parameters, namely, the drained Young's modulus, the void ratio, hydraulic permeability, the Poisson's ratio, the solid grain's and fluid's bulk moduli.

FINDINGS

The peak stress and consolidation were recovered for most of the samples (N=25) by the poroelastic model (normalised root-mean-square error ≤0.03 for fast and slow confined compression conditions).

INTERPRETATION

The strength of the proposed model is its fewer number of variables (N=6 for the proposed poroelastic model versus N=18 for the viscohyperelastic model proposed by Vaidya and collaborators). The incorporation of poroelasticity to clinical models of Pessure Ulcer formation could lead to more precise and mechanistic explorations of soft tissue injury risk factors.

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.