A novel approach for tetrahedral-element-based finite element simulations of anisotropic hyperelastic intervertebral disc behavior.

Intervertebral discs are microstructurally complex spinal tissues that add greatly to the flexibility and mechanical strength of the human spine. Attempting to provide an adjustable basis for capturing a wide range of mechanical characteristics and to better address known challenges of numerical modeling of the disc, we present a robust finite-element-based model formulation for spinal segments in a hyperelastic framework using tetrahedral elements. We evaluate the model stability and accuracy using numerical simulations, with particular attention to the degenerated intervertebral discs and their likely skewed and narrowed geometry. To this end, 1) annulus fibrosus is modeled as a fiber-reinforced Mooney-Rivlin type solid for numerical analysis. 2) An adaptive state-variable dependent explicit time step is proposed and utilized here as a computationally efficient alternative to theoretical estimates. 3) Tetrahedral-element-based FE models for spinal segments under various loading conditions are evaluated for their use in robust numerical simulations. For flexion, extension, lateral bending, and axial rotation load cases, numerical simulations reveal that a suitable framework based on tetrahedral elements can provide greater stability and flexibility concerning geometrical meshing over commonly employed hexahedral-element-based ones for representation and study of spinal segments in various stages of degeneration.

Predicting muscle tissue response from calibrated component models and histology-based finite element models.

Skeletal muscle is an anisotropic soft biological tissue composed of muscle fibres embedded in a structurally complex, hierarchically organised extracellular matrix. In a recent work (Kuravi et al., 2021) we have developed 3D finite element models from series of histological sections. Moreover, based on decellularisation of fresh tissue samples, a novel set of experimental data on the direction dependent mechanical properties of collagenous ECM was established (Kohn et al., 2021). Together with existing information on the material properties of single muscle fibres, the combination of these techniques allows computing predictions of the composite tissue response. To this end, an inverse finite element procedure is proposed in the present work to calibrate a constitutive model of the extracellular matrix, and supplementary biaxial tensile tests on fresh and decellularised tissues are performed for model validation. The results of this rigorously predictive and thus unforgiving strategy suggest that the prediction of the tissue response from the individual characteristics of muscle cells and decellularised tissue is only possible within clear limits. While orders of magnitude are well matched, and the qualitative behaviour in a wide range of load cases is largely captured, the existing deviations point at potentially missing components of the model and highlight the incomplete experimental information in bottom-up multiscale approaches to model skeletal muscle tissue.

Direct measurement of the direction-dependent mechanical behaviour of skeletal muscle extracellular matrix.

This paper reports the first comprehensive data set on the anisotropic mechanical properties of isolated endo- and perimysial extracellular matrix of skeletal muscle, and presents the corresponding protocols for preparing and testing the samples. In particular, decellularisation of porcine skeletal muscle is achieved with caustic soda solution, and mechanical parameters are defined based on compressive and tensile testing in order to identify the optimal treatment time such that muscle fibres are dissolved whereas the extracellular matrix remains largely intact and mechanically functional. At around 18 h, a time window was found and confirmed by histology, in which axial tensile experiments were performed to characterise the direction-dependent mechanical response of the extracellular matrix samples, and the effect of lateral pre-compression was studied. The typical, large variability in the experimental stress response could be largely reduced by varying a single scalar factor, which was attributed to the variation of the fraction of extracellular matrix within the tissue. While experimental results on the mechanical properties of intact muscle tissue and single muscle fibres are increasingly available in literature, there is a lack of information on the properties of the collagenous components of skeletal muscle. The present work aims at closing this gap and thus contributes to an improved understanding of the mechanics of skeletal muscle tissue and provides a missing piece of information for the development of corresponding constitutive and computational models.