Vol.:(0123456789)1 3 Biomechanics and Modeling in Mechanobiology (2023) 22:1499–1514 https://doi.org/10.1007/s10237-022-01664-4 ORIGINAL PAPER On a three–dimensional model for the description of the passive characteristics of skeletal muscle tissue Fabian Walter1 · Robert Seydewitz1 · Philipp Mitterbach2 · Tobias Siebert3 · Markus Böl1 Received: 7 September 2022 / Accepted: 29 November 2022 / Published online: 22 December 2022 © The Author(s) 2022 Abstract In this work, a three–dimensional model was developed to describe the passive mechanical behaviour of anisotropic skeletal muscle tissue. To validate the model, orientation–dependent axial ( 0◦ , 45◦ , 90◦ ) and semi–confined compression experiments (mode I, II, III) were performed on soleus muscle tissue from rabbits. In the latter experiments, specimen deformation is prescribed in the loading direction and prevented in an additional spatial direction, fibre compression at 0◦ (mode I), fibre elongation at 90◦ (mode II) and a neutral state of the fibres at 90◦ where their length is kept constant (mode III). Overall, the model can adequately describe the mechanical behaviour with a relatively small number of model parameters. The stiff- est tissue response during orientation–dependent axial compression ( − 7.7 ± 1.3 kPa) occurs when the fibres are oriented perpendicular to the loading direction ( 90◦ ) and are thus stretched during loading. Semi–confined compression experiments yielded the stiffest tissue ( − 36.7 ± 11.2 kPa) in mode II when the muscle fibres are stretched. The extensive data set col- lected in this study allows to study the different error measures depending on the deformation state or the combination of deformation states. Keywords Skeletal muscle tissue model · Axial compression experiment · Semi–confined compression experiment · Parameter identification · Passive skeletal muscle behaviour 1 Introduction The mechanical behaviour of passive skeletal muscle tis- sue exhibits some unique characteristics that distinguish it from ordinary composite materials. Although the mechani- cal behaviour is well documented in the literature through numerous experimental investigations, the load–bearing mechanism of the complicated three–dimensional and hier- archical structure remains largely unknown. The interaction of muscle fibres and the collagenous, extracellular matrix (ECM) is responsible for the generation and transmission of active and passive forces and to resist various external loads. While the main function of muscle fibres is to gener- ate active forces, the ECM serves as a connector to transfer load across different length scales (Böl et al. 2014, 2022). Experiments were performed on proteins (Linke et al. 1998; Baker et al. 2002; Pertici et al. 2018; Li et al. 2020), sar- comeres/myofibrils (Shalabi et al. 2017; de Souza Leite et al. 2017; Haeger and Rassier 2020; Lee et al. 2020; Hae- ger et al. 2020; Swist et al. 2020; Ward et al. 2020; Scel- lini et al. 2021; Marston 2022), fibres (Bartoo et al. 1997; Toursel et al. 2002; Fridén and Lieber 2003; Noonan et al. 2020; Böl et al. 2019, 2020, 2022), fibre bundles (Mutungi and Ranatunga 1996; Lieber et al. 2003; Meyer and Lieber 2011; Brown et al. 2012; Mathewson et al. 2014; Wood et al. 2014; Tamura et al. 2016; Wu et al. 2016; Noonan et al. 2020; Ward et al. 2020), and tissue level (Calvo et al. 2010; Morrow et al. 2010; Böl et al. 2012; Gras et al. 2012; Takaza et al. 2013b, a; Böl et al. 2014, 2016; Hashemi et al. 2020; Wheatley 2020; Leichsenring et al. 2021; Kohn et al. 2021; Kuravi et al. 2021; Böl et al. 2022) to characterise the mechanical properties at different scales. As previous studies at fibre (Böl et al. 2019, 2020) and tissue scale (Böl et al. 2014, 2016) have shown, individual fibres as well as * Markus Böl m.boel@tu-braunschweig.de 1 Institute of Mechanics and Adaptronics, Technische Universität Braunschweig, D–38106 Braunschweig, Germany 2 Mechanical Engineering, Eindhoven University of Technology, NLD–5612 Eindhoven, The Netherlands 3 Institute of Sport and Motion Science, University of Stuttgart, D–70569 Stuttgart, Germany http://crossmark.crossref.org/dialog/?doi=10.1007/s10237-022-01664-4&domain=pdf 1500 F. Walter et al. 1 3 fibres in combination with the ECM are capable of carrying compressive loads, suggesting a more comprehensive theory for the mechanical behaviour of muscle fibres and their con- tribution to the overall load response. Constitutive models originally developed for tissue materials such as the arterial wall (Holzapfel and Ogden 2010; Gasser 2018) or the heart muscle (Humphrey et al. 1990; Ambrosi and Pezzuto 2012), skeletal muscle (Martins et al. 1998; Johannson et al. 2000; Blemker et al. 2005; van Loocke et al. 2006; Röhrle and Pullan 2007; Odegard et al. 2008; Tang et al. 2009; Calvo et al. 2010; Ehret et al. 2011; Moerman et al. 2016), but also general formulations for fibre–reinforced soft tissue materi- als (Balzani et al. 2006; Ehret and Itskov 2007; Federico and Gasser 2010) use phenomenological approaches based on a hyperelastic strain energy function to predict the character- istic anisotropic load behaviour of soft tissue materials. Due to their relatively simple structure, they are easy to imple- ment and adapt to different tissue types. A comprehensive overview of constitutive equations for various types of soft tissue materials are presented in the review article of Chag- non et al. (2015). To account for the anisotropic behaviour, additional so–called pseudo–invariants are often used in the constitutive relations, which take into account the deforma- tion of the oriented fibres. Neglecting the ability of the fibre material to carry compressive loads, but also to account for different loading responses in tension and compression, the anisotropic contribution of the strain energy function during fibre compression is suppressed by using an additional on/ off or heaviside function (Holzapfel et al. 2000; Federico and Gasser 2010). Model calibration often uses unloaded uniaxial compression tests in which the fibres are either compressed or stretched as a result of lateral expansion. There are only a few studies with so–called semi–confined compression experiments (Böl et al. 2014, 2016; Leichsen- ring et al. 2021). Here, the experimental setup allows, among other situations, the fibre length to be kept constant in order to investigate the anisotropic material behaviour under pres- sure more closely. In this sense, Böl et al. (2014) defined the fibre states under compression, tension, and at constant length as mode I, II, and III, respectively. Despite the dif- ferences, the experimental results of Böl et al. (2014) and Leichsenring et al. (2021) show different material responses for mode I and III, which illustrates the contribution of the muscle fibres to the load transfer of the entire tissue. It is easy to see that all models with the anisotropic con- tribution disabled during compression are not able to distin- guish between mode I and III , as the load transfer behaviour of the muscle fibres is not taken into account in these for- mulations. Therefore, it is not possible to find a set of mate- rial parameters that can predict the material response for all three fibre states. But also continuous approaches (Gasser et al. 2006) have difficulties in predicting the material behav- iour. Consequently, the aim of the present work is to develop a simple phenomenological modelling approach at tissue level to predict the material response under different loading conditions and in different fibre directions within the limits of compression. Avoiding an on/off relationship allows a continuous transition between stress and compression for different fibre orientations. To verify the performance of the model, we apply the model to experimental studies under axial compression and semi–confined compression. 2 Materials and methods 2.1 Ethical approval The study was exempted from ethical committee review according to national regulations (German Animal Wel- fare Act), as healthy, female domestic rabbits (Oryctolagus cuniculus) were obtained from a rabbit farm immediately after animal sacrifice. 2.2 Tissue sample dissection and processing Following Leichsenring et al. (2021), experiments were per- formed on muscle cubes taken from the soleus muscles of 18 female, white New Zealand rabbits (weight: 2700 ± 598 g, age: 75.2 ± 16.1 days). After transporting the animals in a cool box (at a temperature of 4◦ C) to the laboratory, the soleus muscles were dissected out. The rabbit soleus has a simple unipennate muscle architecture with a mean optimal muscle fibre length of about 22 mm (Siebert et al. 2015), which allows cutting of tissue samples with fibre orienta- tions parallel to the cube edge and at 45◦ . To stabilise the samples during cutting, the tissue was embedded in algi- nate, allowing the cutting of cubic samples with a mean edge length of 4.36 ± 0.56 mm. After preparation, the tis- sue samples were wrapped in cloths soaked with Dulbecco’s phosphate–buffered saline (DPBS) and stored in a climate chamber at 4◦ C and 50% humidity. To warm the samples to a physiological temperature of 38.5◦ C for the mechani- cal tests (Cooper et al. 1965; Siebert et al. 2017), they were removed from the climate chamber before the mechanical tests and placed in a bath of DPBS solution for 30 seconds. 2.3 Experimental investigation on skeletal muscle tissue As shown in previous work (e.g. Böl et al. 2013, 2014, 2016), the deformation state during muscle activation is highly inhomogeneous, indicating a complex load distri- bution within the muscle tissue where tension and com- pression are present. To gain a better understanding of the load–bearing mechanisms, in this study, we look at the material response of muscle tissue samples under axial and 1501On a three–dimensional model for the description of the passive characteristics of skeletal… 1 3 semi–confined compression for different muscle fibre ori- entations, see Fig. 1. In the axial compression experiment (a), specimens are tested with their fibres oriented 0◦ , 45◦ , or 90◦ with respect to the loading direction. In semi–confined compression (b), where specimen deformation is prescribed in the load direc- tion and prevented in an additional spatial direction, fibre compression at 0◦ (mode I), fibre elongation at 90◦ (mode II) and a neutral state of the fibres at 90◦(mode III) are exam- ined. Table 1 shows an overview of the experiments per- formed in this study. 2.3.1 Axial compression experiments Following our previous investigation by Leichsenring et al. (2021), forty–two ( n = 42 ) cubic samples were subjected to orientation–dependent, axial compression (index: ac) testing up to 40% nominal compressive strain � at a strain rate of �̇� = 0.5 %s−1 . Specifically, three fibre orientations ( 0◦ , 45◦ , 90◦ ), to be defined as the angle between the fibre and the applied loading direction, were considered, see also Fig. 1a. While the displacement u of the upper platen was predefined, the resulting force F was measured and con- verted to mean engineering stress via P = F∕Aac by divi- sion through the axial cross–sectional area Aac , measured from a digital image which was recorded before testing. The compressive strain � = u∕h0 was calculated from u and the undeformed sample height h0 . To diminish friction effects, the polytetrafluoroethylene (PTFE) platens were coated with silicone oil for lubrication (Böl et al. 2012). Additionally, a camera was positioned in front of the samples recording the sample’s transversal deformation. 2.3.2 Semi–confined compression experiments Fourty–four ( n = 44 ) semi–confined (index: sc), orienta- tion–dependent compression tests were performed using a custom–made testing device (Leichsenring et al. 2021). Three orientations, in which the fibres are aligned in the x–, y–, and z–directions, were considered, see Fig. 1b. These specimen orientations cause the fibre to be compressed (mode I), stretched (mode II), and held constant (mode III) during deformation. For details on the experimental setup, we refer to previous investigations (Böl et al. 2014, 2015a, 2016; Leichsenring et al. 2021). To sample the specimens, they were carefully placed in the device, taking into account the fibre orientation. To reduce friction effects between the tissue and the device, the contact surfaces were sprayed with low viscosity sili- cone oil prior to the test, which had proven to be very effective in previous studies (Böl et al. 2014, 2015a, 2016; Leichsenring et al. 2021). Similar to the axial compression experiments, the tissue was compressed at a constant strain rate of �̇� = 0.5 %s−1 to a maximum compression level of 40%. Again, the displacement u of the upper plunger was predefined, while the resulting force F was measured and (a) (b) u h0 u h0 plungertissue sample force sensor force sensorplunger PTFE platens 3 cm 3 cm 10 mm 10 mm tissue sample y x z Fig. 1 Deformation state, realised within this study: a Axial compression and b semi–confined compression experiments. All experiments were conducted with an axial testing machine (Zwick Z0.5, Zwick GmbH & Co. Ulm, Germany), equipped with a 5 N or 10 N load cell Table 1 Overview of the number of experiments realised in this study Axial compression Semi–confined compression Orientation 0◦ 45◦ 90◦ I II III No. of samples 12 16 14 14 15 15 Total 86 1502 F. Walter et al. 1 3 converted to mean engineering stress P = F∕Asc , where Asc defines the axial cross–sectional area of the undeformed specimen. The axial strain � = u∕h0 was calculated from u and the undeformed sample height h0. 3 Continuum mechanical framework of passive muscle tissue modelling 3.1 Kinematics Following general principles of continuum mechanics, we introduce B0 and B as the reference and current configura- tions of a body, in which the position of the particle is given by the vectors X and x = �(X) , respectively. Herein, the non–linear mapping of a material particle from the reference to the current configuration is represented by �(X, t) and the corresponding deformation gradient F = ∇X�(X, t) . The Jacobian J = detF > 0 represents the volume change of the material particle and ∇X(∙) indicates the spatial derivative with respect to the reference coordinates. Accounting for the anisotropic behaviour of the tissue material, we introduce the fibre direction vector M and the corresponding structural tensor in the reference configuration. The corresponding muscle fibre stretch depends on the right Cauchy-Green tensor C = FTF and is in relation to the fourth invariant I4. Due to the high amount of water in living tissues, it is justified to assume incompressible material behaviour (Baskin and Paolini 1967). In computational mechanics it is often convenient to allow for slight volumetric changes and accordingly, to consider nearly incompressible material. The basic step for this procedure is a multiplicative split of the total deformation gradient into volumetric and isochoric parts according to Flory (1961) as where F̄ is the isochoric part of F . Accordingly, we obtain the isochoric counterpart of the right Cauchy–Green tensor as as well as the corresponding scalar valued first invariant (1)Z = M ⊗M (2)�2 = I4 = C ∶ Z (3)F = J1∕3F̄, (4)C̄ = J−2∕3C = F̄ T F̄ (5)Ī1 = trC̄, representing an isotropic measure for the distortion of the matrix–fibre material. 3.2 Constitutive equation for fibrous material under compression Based on the results of recent experimental investigations, see Sects. 4.1 and 4.2, we propose a phenomenological modelling approach to describe the stress/strain behaviour during axial and semi–confined compression for different fibre orienta- tions. For this purpose, we assume homogeneous deforma- tions with perfectly aligned fibres in the direction of the main axes. In contrast to axial compression tests, the movement in semi–confined experiments is restricted in two directions, resulting in fibre being compressed (mode I, 𝜆 < 1 ), stretched (mode II, 𝜆 > 1 ), and held at a constant length (mode III � = 1 ). It is assumed that in mode I the compression counter- acts the structural integrity of the fibre material, resulting in a weaker material response than in mode III, where the fibre length remains constant. To allow a smooth transition between the different modes and to predict the material behaviour for realistic deformations containing a combination of the three modes, a continuous strain energy function is proposed in the framework of hyperelasticity as Herein, the material parameters c1 and c3 , c2 , and K vol con- trol the general load response, the anisotropic behaviour, and volume preservation, respectively. The bulk modulus K vol was set to 5000 times the shear modulus of the tissue, resulting in nearly incompressible behaviour. Further, Equa- tion (6) depends on the fourth invariant ( I4 ) and the isochoric part ( ̄I1 ) of the first invariant, see also Remark 1. The corresponding stress in terms of the 1st Piola-Kirchhoff stress follows from the derivation of the strain energy function with respect to the right Cauchy-Green tensor as The variable load behaviour for mode  I, II, and III is mainly controlled by g(I4) , since any deformation F ≠ I with detF > 0 leads to f (Ī1) > 0 , where the particular (6) 𝜓(I4, Ī1) = c1 e c2(I4−1) ��� g(I4) (ec3(Ī1−3) − 1) ������������� f (Ī1) +K vol (J2 − 1 − 2 ln J). (7) P = 2F 𝜕𝜓 𝜕C = 2c1 ( exp[c2(I4 − 1) + c3(Ī1 − 3)] [ c2FM + c3J −2∕3 ( F − 1 3 trCF −T )] − c2 exp[c2(I4 − 1)]FM ) + 2K vol (J 2 − 1)F−T . 1503On a three–dimensional model for the description of the passive characteristics of skeletal… 1 3 formulation of f (Ī1) satisfies the stress–free reference con- dition of � . A common approach to modelling the aniso- tropic material behaviour is an additive split of the strain energy function into an isotropic and an anisotropic con- tribution that includes the fibre strain I4 (e.g. Balzani et al. 2006; Moerman et al. 2016). However, deactivating the anisotropic component for I4 < 1 excludes the fibres during compression, limiting the ability of the model formulation to predict the material response for all three modes and thus general deformations. For the current approach, g(I4) < 1 in mode I reduces the overall load response in the strain energy function � . In contrast, in mode II, g(I4) > 1 leads to an increase in � , while in mode III, the length of the muscle fibre remains constant and g(I4) = 1 , leading to 𝜓 = 𝜓(Ī1, J) . Note, Equation (6) is not convex due to fibre compression, which can lead to undesirable effects in a finite element sim- ulation depending on the relationship between the material parameters c2 and c3 and the deformation present. However, the main objective of the present work was to find a simple structure for a strain energy function capable of capturing the different states of fibre compression. Remark 1 (Volumetric/isochoric split of F ) In Equation (6) and consequently in (7) we use the combination of Ī1 and I4 to implement one–dimensional stresses in the fibre structure with respect to the reference configuration, which is con- sistent with the simplified representation that fibres act as one–dimensional truss elements. As already mentioned by Sansour (2008), the use of Ī4 does not imply a one–dimen- sional stress condition. Therefore, the use of Ī4 is limited to purely incompressible materials where the incompressibility condition is enforced by Lagrange multipliers (Nolan et al. 2014). In the present work, an alternative approach has been used that allows for slight volume changes using a penalty function and a static parameter. 3.3 Optimisation scheme for material parameter identification To achieve successful identification of the material param- eters, an inverse numerical optimisation procedure was introduced. In this process, a forward finite element analy- sis is coupled with an optimisation algorithm that deter- mines the optimal set of parameters by adjusting the mate- rial parameters step by step to minimise the error between the measured and simulated responses, further denoted as inverse finite element method (iFEM). The basis for such parameter identification is, on the numerical side, the model approach for describing the mechanical behaviour of the anisotropic muscle tissue shown in Equation (6) and, on the experimental side, two deformation states (axial compression, semi–confined compression) described in Sect. 2.3 for three fibre orientations (axial compression: 0◦ , 45◦ , 90◦ and semi–confined compression: mode  I, mode II, mode III) each. Thus, the agreement between simulated and experimental data is evaluated using the objective function which represents the combined relative square error of the included compression experiments. Herein, iPsim j and iPexp j are the orientation dependent simulated and measured stress values for every deformation increment, respectively. The use of silicone oil as a lubricant reduces frictional effects to a minimum and the resulting friction has been shown to be negligible (Böl et al. 2012). Therefore, friction effects are neglected in the FEM analyses within the scope of this study. Considering parallel ( 0◦ ) and perpendicular ( 90◦ ) oriented fibres, the axial and semi–confined compression tests within the friction limit lead to a homogeneous stress state, which allows the use of a single element for the FEM analysis. To achieve a compression of 40%, the FEM procedure performs 50 continuous load steps. For the subsequent adjustment of the strain measuring points of the calculation data, linear interpolation within measuring points is applied to both data sets. In addition, m takes into account the number of experi- ments (inform of stress–stretch relations) used simultane- ously in the optimisation and n accounts for the numer of data sets used in the respective experiment. The following combinations were used in this study, see also Sect. 4.3: First, each of the 86 experiments was considered individu- ally in the optimisation, resulting in m = 1 . In addition, the 3 directions of the axial compression and the semi–confined compression experiments (each in the form of mean value curves, see Figs. 2 and 3) were taken into account in the optimisation, resulting in a value of m = 3 in each case. Finally, the mean stress–stretch relations of all directions and deformation states were considered simultaneously, so that m = 6 . To do this, mean value curves are created by first performing an interpolation for each data set and then cal- culating the mean stress at each sample point. To find an optimal parameter set in terms of the parameter vector p = (c1, c2, c3) T the objective function was minimised, using the Nelder–Mead simplex algorithm (fminsearch, matlabⓇ R2018a) (Nelder and Mead 1965; Lagarias et al. 1998). Based on this algorithm, parameter sets were considered optimal when the objective function changed by less than the tolerance between two consecutive fitting steps. (8)O(p) = m ∑ i=1 n ∑ j=1 1 mn √ √ √ √ (iPsim j −iP exp j )2 iP exp j , 1504 F. Walter et al. 1 3 4 Results 4.1 Axial compression experiments Overall, the results of the axial compression tests show a non–linear, exponential material behaviour, independent of loading direction, see Fig. 2. The maximum standard deviations consistently occur at the maximum compressive stress value. Expressed as a percentage of the corresponding mean value, i.e. s.d./ mean ×100, the maximum standard deviations are less than 44%. Furthermore, it can be seen that the stiffest tissue response occurs for samples in which the fibres are ori- ented perpendicular to the loading direction, i.e., 90◦ . The softest response, on the other hand, is seen for tissue that is loaded parallel to the fibre orientation, i.e., 0◦ . There is a tendency for significant differences to exist between the 0◦ and 90◦ as well as between the 45◦ and 90◦ samples. It is particularly noticeable that the mean stress at 40% compression is –7.7±1.3 kPa for the sample set at 90◦ , and therefore about 3 times higher than for the sample set at 0◦ , to be –2.5±1.1 kPa, indicating strong anisotropic characteristics. 4.2 Semi–confined compression experiments Figure 3 illustrates the results of the orientation–dependent semi–confined compression experiments in form of mean stress–stretch curves and standard deviations. Independent of the fibre orientation all curves are charac- terised by a clearly non–linear material response. For all ori- entations, the highest standard deviations are found at maxi- mum compression. In general, these standard deviations are smaller than 43% of the corresponding mean value. Addi- tionally, from Fig. 3 it is clear to see, that mode II, where the muscle fibres are stretched, shows with –36.7±11.2 kPa the stiffest response at 40% of compression. Here, mode Fig. 2 Stress–strain curves of the axial compression tests in dependence on the fibre orientation: a 0◦ , b 45◦ , and c 90◦ . Black curves indicate mean values, the shaded areas depict the standard deviation, and grey curves identify the single measurements Fig. 3 Stress–strain curves of the semi–confined compression tests in dependence on the different modes: a Mode I, b mode II , and c mode III. Black curves indicate mean values, the shaded areas depict the standard deviation, and grey curves identify the single measurements 1505On a three–dimensional model for the description of the passive characteristics of skeletal… 1 3 III features the second stiffest response (–16.5±7.1 kPa), here the fibre length remains constant during deformation, followed by mode I (–10.4±3.5 kPa), where the fibres are compressed. 4.3 Optimised material parameters Using the parameter identification scheme described in Sect. 3.3, a wide variety of parameter sets were identified. In a first step, all experiments were optimised individually, resulting in 86 parameter sets, see Table 2. Within a deformation state, these parameter sets show similar parameter values overall. Only in mode II parameter values deviate decisively. This is also evident in the error measures. While these are very low overall, they increase significantly in mode II. Since the identified individual parameter sets as described above only inadequately describe the mechanical behaviour of all experiments, parameters from further experimental set (combinations) were identified. Thus, material parameters were identified only from the (i) axial compression experi- ments ( 0◦ , 45◦ , 90◦ , i.e. related to Equation (8) m = 3 ), (ii) from the semi–confined compression experiments (mode I, mode II, mode III, i.e. m = 3 ), related to Equation (8) m = 3 ), and (iii) from all axial compression experiments ( 0◦ , 45◦ , 90◦ ) and semi–confined compression experiments (mode I, mode II, mode III). Thus m = 6 in Equation (8). The latter case (iii) represents the most realistic situation, since on average the parameters for all experiments should be adequately represented here. The results of the three cases are shown in Table 3. Experimental set (i) yields the lowest error with O = 0.122 . The highest error, with O = 0.512 , is calcu- lated for (ii) and especially the parameters c1 and c3 deviate strongly from those in parameter set (i). Experimental set (iii), with an error measure of O = 0.473 , lies between cases (i) and (ii) and the parameters are similar in magnitude to those of case (i). Figure 4 shows the model response for experimental set (iii), providing the most meaningful parameters, compared to the axial compression experiments. Basically, the model can reproduce all three directions of the deformation state well. While the model for the 0◦ direction shows a slightly too stiff behaviour beyond a stretching value of approximately 0.75, the model for the 45◦ specimens generally shows a slightly stiffer material response. The 90◦ model response, on the other hand, shows very good agreement with the experimen- tal response over the entire stretching range. Figure 5 shows the model prediction for the semi–con- fined compression experiments using the same material parameters as for the axial compression simulations shown in Figure 4. While the model can reproduce the tissue response for mode I and III very well, the model for mode II deviates sig- nificantly by overestimating the progressive increase in the stress response. The model shows a too soft response up to an elongation value of approximately 0.72. Compared to the experimental data, a stronger curvature of the stress–strain relationship can also be observed. 5 Discussion 5.1 Mechanical muscle tissue characteristics To mechanically characterise the passive behaviour of the muscle tissue, in a first step orientation–dependent axial compression tests were performed, cf. Fig. 2. The stiffest tissue response (–7.7 ± 1.3 kPa) occurs when the fibres are oriented perpendicular to the loading direction ( 90◦ ) and are thus stretched during loading. Comparing the different fibre orientations, the mean stress for the 90◦ orientation is 188% and 303% larger compared to the samples with a 45◦ fibre orientation (–4.1 ± 2.0 kPa) and with a 0◦ fibre orien- tation (–2.5 ± 1.1 kPa), respectively. This behaviour agrees well with some studies carried out on muscle tissue from pigs (van Loocke et al. 2006), rabbits (Böl et al. 2012), or chickens (Mohammadkhah et al. 2016), in which the passive properties of muscle tissue were determined as a function of fibre orientation. The anisotropic material behaviour of the muscle tissue was also reported by Pietsch et al. (2014) and Simms et al. (2016), investigating compressive loading of porcine muscle tissue in two different directions ( 0◦ and 90◦ ). Both studies also determined lower stresses when the fibres are oriented in line to the loading direction. In a second step, semi–confined compression experiments were performed to describe the mechanical properties of the muscle tissue. The stiffest tissue response occurs with –36.7 ± 11.2 kPa when the muscle fibres are stretched, i.e. mode II, see Fig. 3. This observation is consistent with tis- sue measurements performed on biceps femoris muscles of female domestic pigs (Böl et al. 2014, 2016) and soleus muscle of rabbits (Leichsenring et al. 2021). For mode I, the tissue response in the present study is softest at –10.4 ± 3.5 kPa, which is consistent with the study by Leichsenring et al. (2021), which also sampled soleus muscle tissue from rabbits. In contrast, Böl et al. (2014) and Böl et al. (2016) found the lowest stresses when the fibre lengths remained constant during the compression test, i.e., at mode III. The reason for the different mechanical results between the aforementioned studies is thought to be the ratio of fibres and ECM in the specific muscle tissue. For example, it is known that skeletal muscles within an animal or between different species differ both in their ECM (Gillies and Lieber 2011) and in the intracellular passive muscle structures (e.g. 1506 F. Walter et al. 1 3 Ta bl e 2 M at er ia l p ar am et er s a nd e rr or m ea su re s d et er m in ed fr om si ng le , i nd iv id ua l e xp er im en ts A xi al c om pr es si on Se m i– co nfi ne d co m pr es si on Se t n o. 0 ◦ 4 5 ◦ 9 0 ◦ M od e I M od e II M od e II I O c 1 [k Pa ] c 2 [– ] c 3 [– ] c 1 [k Pa ] c 2 [– ] c 3 [– ] c 1 [k Pa ] c 2 [– ] c 3 [– ] c 1 [k Pa ] c 2 [– ] c 3 [– ] c 1 [k Pa ] c 2 [– ] c 3 [– ] c 1 [k Pa ] c 2 [– ] c 3 [– ] [– ] I 0. 47 1. 34 1. 41 – – – – – – – – – – – – – – – 0. 04 0 II 1. 54 1. 72 0. 91 – – – – – – – – – – – – – – – 0. 00 8 II I 0. 85 1. 29 1. 18 – – – – – – – – – – – – – – – 0. 04 0 IV 0. 99 0. 75 0. 82 – – – – – – – – – – – – – – – 0. 02 7 V 0. 62 1. 16 0. 84 – – – – – – – – – – – – – – – 0. 00 4 V I 2. 31 2. 69 0. 84 – – – – – – – – – – – – – – – 0. 01 2 V II 3. 40 0. 01 0. 13 – – – – – – – – – – – – – – – 0. 05 8 V II I 6. 13 1. 26 0. 20 – – – – – – – – – – – – – – – 0. 03 8 IX 1. 45 0. 39 0. 69 – – – – – – – – – – – – – – – 0. 01 4 X 0. 15 1. 35 1. 73 – – – – – – – – – – – – – – – 0. 01 1 X I 0. 83 1. 47 0. 71 – – – – – – – – – – – – – – – 0. 05 1 X II 0. 82 1. 12 0. 79 – – – – – – – – – – – – – – – 0. 00 7 X II I – – – 0. 16 0. 00 1. 63 – – – – – – – – – – – – 0. 04 9 X IV – – – 0. 19 0. 00 1. 32 – – – – – – – – – – – – 0. 05 3 X V – – – 0. 14 0. 00 2. 33 – – – – – – – – – – – – 0. 08 7 X V I – – – 0. 12 1. 09 1. 66 – – – – – – – – – – – – 0. 01 7 X V II – – – 0. 88 2. 70 0. 95 – – – – – – – – – – – – 0. 05 2 X V II I – – – 0. 12 0. 00 1. 46 – – – – – – – – – – – – 0. 02 2 X IX – – – 0. 36 1. 52 1. 17 – – – – – – – – – – – – 0. 02 1 X X – – – 0. 20 0. 00 1. 84 – – – – – – – – – – – – 0. 03 4 X X I – – – 0. 14 0. 00 2. 51 – – – – – – – – – – – – 0. 06 3 X X II – – – 0. 64 0. 23 1. 26 – – – – – – – – – – – – 0. 04 0 X X II I – – – 0. 10 0. 27 2. 27 – – – – – – – – – – – – 0. 00 9 X X IV – – – 0. 16 0. 00 1. 95 – – – – – – – – – – – – 0. 04 0 X X V – – – 0. 11 0. 00 2. 33 – – – – – – – – – – – – 0. 05 2 X X V I – – – 0. 36 1. 16 0. 95 – – – – – – – – – – – – 0. 00 8 X X V II – – – 0. 11 0. 00 1. 84 – – – – – – – – – – – – 0. 01 6 X X V II I – – – 0. 32 0. 00 1. 43 – – – – – – – – – – – – 0. 21 9 X X IX – – – – – – 0. 12 1. 54 1. 81 – – – – – – – – – 0. 08 6 X X X – – – – – – 0. 16 1. 56 1. 59 – – – – – – – – – 0. 10 7 X X X I – – – – – – 0. 33 1. 57 0. 96 – – – – – – – – – 0. 13 0 X X X II – – – – – – 0. 20 1. 55 1. 54 – – – – – – – – – 0. 08 5 X X X II I – – – – – – 0. 69 1. 58 0. 67 – – – – – – – – – 0. 05 5 X X X IV – – – – – – 0. 07 0. 01 3. 13 – – – – – – – – – 0. 11 3 X X X V – – – – – – 0. 40 1. 58 1. 07 – – – – – – – – – 0. 15 2 1507On a three–dimensional model for the description of the passive characteristics of skeletal… 1 3 Ta bl e 2 (c on tin ue d) A xi al c om pr es si on Se m i– co nfi ne d co m pr es si on Se t n o. 0 ◦ 4 5 ◦ 9 0 ◦ M od e I M od e II M od e II I O c 1 [k Pa ] c 2 [– ] c 3 [– ] c 1 [k Pa ] c 2 [– ] c 3 [– ] c 1 [k Pa ] c 2 [– ] c 3 [– ] c 1 [k Pa ] c 2 [– ] c 3 [– ] c 1 [k Pa ] c 2 [– ] c 3 [– ] c 1 [k Pa ] c 2 [– ] c 3 [– ] [– ] X X X V I – – – – – – 0. 19 1. 56 1. 49 – – – – – – – – – 0. 12 5 X X X V II – – – – – – 0. 11 1. 62 1. 97 – – – – – – – – – 0. 09 2 X X X V II I – – – – – – 0. 10 0. 22 2. 62 – – – – – – – – – 0. 06 8 X X X IX – – – – – – 0. 11 1. 55 1. 92 – – – – – – – – – 0. 04 0 X L – – – – – – 0. 37 1. 54 0. 90 – – – – – – – – – 0. 04 1 X LI – – – – – – 0. 22 1. 56 1. 56 – – – – – – – – – 0. 10 3 X LI I – – – – – – 0. 52 1. 57 0. 98 – – – – – – – – – 0. 12 1 X LI II – – – – – – – – – 0. 80 0. 14 0. 90 – – – – – – 0. 05 8 X LI V – – – – – – – – – 1. 56 0. 92 0. 84 – – – – – – 0. 07 8 X LV – – – – – – – – – 0. 90 0. 27 0. 72 – – – – – – 0. 07 5 X LV I – – – – – – – – – 1. 31 0. 00 0. 44 – – – – – – 0. 06 9 X LV II – – – – – – – – – 4. 60 1. 09 0. 40 – – – – – – 0. 03 7 X LV II I – – – – – – – – – 9. 33 0. 14 0. 07 – – – – – – 0. 09 2 X LI X – – – – – – – – – 2. 74 0. 75 0. 45 – – – – – – 0. 03 1 L – – – – – – – – – 0. 58 0. 00 1. 16 – – – – – – 0. 09 9 LI – – – – – – – – – 1. 50 0. 00 0. 37 – – – – – – 0. 07 0 LI I – – – – – – – – – 2. 09 1. 41 0. 61 – – – – – – 0. 05 4 LI II – – – – – – – – – 1. 24 1. 02 0. 87 – – – – – – 0. 04 6 LI V – – – – – – – – – 1. 48 1. 17 1. 02 – – – – – – 0. 06 9 LV – – – – – – – – – 1. 70 0. 86 0. 62 – – – – – – 0. 04 5 LV I – – – – – – – – – 1. 82 0. 00 0. 40 – – – – – – 0. 10 0 LV II – – – – – – – – – – – – 33 .6 8 0. 01 0. 06 – – – 0. 11 9 LV II I – – – – – – – – – – – – 14 .9 0 0. 00 0. 27 – – – 0. 31 2 LI X – – – – – – – – – – – – 22 1. 16 0. 00 0. 02 – – – 0. 77 2 LX – – – – – – – – – – – – 21 4. 93 0. 18 0. 02 – – – 1. 05 8 LX I – – – – – – – – – – – – 2. 07 0. 06 0. 89 – – – 0. 56 7 LX II – – – – – – – – – – – – 20 2. 29 0. 00 0. 02 – – – 0. 47 8 LX II I – – – – – – – – – – – – 73 0. 50 0. 00 0. 01 – – – 0. 47 0 LX IV – – – – – – – – – – – – 22 3. 62 0. 42 0. 01 – – – 1. 22 8 LX V – – – – – – – – – – – – 4. 12 0. 38 0. 31 – – – 0. 82 9 LX V I – – – – – – – – – – – – 15 96 .5 7 0. 00 0. 00 – – – 0. 81 4 LX V II – – – – – – – – – – – – 4. 44 0. 32 0. 29 – – – 0. 47 0 LX V II I – – – – – – – – – – – – 22 3. 56 0. 17 0. 02 – – – 0. 86 9 1508 F. Walter et al. 1 3 Ta bl e 2 (c on tin ue d) A xi al c om pr es si on Se m i– co nfi ne d co m pr es si on Se t n o. 0 ◦ 4 5 ◦ 9 0 ◦ M od e I M od e II M od e II I O c 1 [k Pa ] c 2 [– ] c 3 [– ] c 1 [k Pa ] c 2 [– ] c 3 [– ] c 1 [k Pa ] c 2 [– ] c 3 [– ] c 1 [k Pa ] c 2 [– ] c 3 [– ] c 1 [k Pa ] c 2 [– ] c 3 [– ] c 1 [k Pa ] c 2 [– ] c 3 [– ] [– ] LX IX – – – – – – – – – – – – 6. 71 0. 19 0. 38 – – – 0. 40 3 LX X – – – – – – – – – – – – 11 .2 0 0. 01 0. 26 – – – 0. 22 1 LX X I – – – – – – – – – – – – 2. 48 0. 24 0. 57 – – – 0. 17 6 LX X II – – – – – – – – – – – – – – – 4. 07 2. 95 0. 16 0. 11 5 LX X II I – – – – – – – – – – – – – – – 0. 61 0. 99 1. 53 0. 25 3 LX X IV – – – – – – – – – – – – – – – 0. 42 0. 67 1. 28 0. 09 0 LX X V – – – – – – – – – – – – – – – 0. 39 1. 42 1. 44 0. 14 4 LX X V I – – – – – – – – – – – – – – – 1. 22 1. 36 0. 81 0. 22 4 LX X V II – – – – – – – – – – – – – – – 0. 58 0. 35 0. 82 0. 07 8 LX X V II – – – – – – – – – – – – – – – 0. 74 1. 34 1. 12 0. 13 4 LX X IX – – – – – – – – – – – – – – – 0. 63 0. 47 0. 86 0. 05 0 LX X X – – – – – – – – – – – – – – – 0. 13 2. 12 1. 58 0. 06 7 LX X X I – – – – – – – – – – – – – – – 0. 75 0. 30 1. 00 0. 10 2 LX X X II – – – – – – – – – – – – – – – 1. 13 0. 73 1. 07 0. 25 1 LX X X II I – – – – – – – – – – – – – – – 0. 23 0. 10 1. 76 0. 17 2 LX X X IV – – – – – – – – – – – – – – – 0. 17 1. 58 1. 41 0. 08 5 LX X X V – – – – – – – – – – – – – – – 0. 24 1. 11 1. 49 0. 24 7 LX X X V I – – – – – – – – – – – – – – – 0. 57 1. 02 1. 18 0. 12 5 M ea n 1. 63 1. 21 0. 85 0. 26 0. 44 1. 68 0. 26 1. 36 1. 59 2. 26 0. 56 0. 63 23 2. 82 0. 13 0. 21 0. 79 1. 10 1. 17 0. 17 S. D 1. 60 0. 64 0. 43 0. 21 0. 76 0. 48 0. 18 0. 51 0. 66 2. 18 0. 50 0. 29 40 8. 84 0. 15 0. 25 0. 93 0. 72 0. 39 0. 24 1509On a three–dimensional model for the description of the passive characteristics of skeletal… 1 3 titin) (Prado et al. 2005). In addition, the ECM structure and amount changes with age (Ramaswamy et al. 2011; Leich- senring et al. 2021) and disease (Lieber and Ward 2013). These factors can lead to different directional behaviour of the muscle tissue under compression. The main aim of Böl et al. (2014) was to analyse the load transfer mechanisms between fibres and ECM. Since the muscles of different animals cited here may differ in their structure and thus in the relationship between fibres and ECM, the mechanical responses also vary. Furthermore, Böl et al. (2014) and Böl et al. (2016) describe significant differences between all three modes, which can also be confirmed in the present study and thus reveals information about the load transfer behaviour of the muscle fibres. 5.2 Comparison with existing muscle models As described in Sect. 4.3, the model proposed in this study is able to satisfactorily reproduce the deformation states shown here. However, to place this result, and thus the present mus- cle model, in the landscape of existing muscle models, the proposed material model formulation is compared with other skeletal muscle models. For the comparison given here, the typical models published in the field of passive muscle mod- elling were used, assuming ideal incompressibility for all Table 3 Material parameters and error measures determined from various experiments/experiment combinations Parameter set p Experi- mental set Considered deformation state(s) c 1 [kPa] c 2 [–] c 3 [–] O(p) (i) Axial compression ( 0◦ , 45◦ , 90◦) 0.24 0.70 1.56 0.122 (ii) Semi–confined compres- sion (modes I, II, III) 1.39 0.52 0.66 0.512 (iii) (i) & (ii) 0.39 0.53 1.27 0.473 Fig. 4 Comparison of the data from the axial compression experi- ments with the model response (with parameters from the experimen- tal set (iii), see Table 3) for a 0◦ , b 45◦ , and c 90◦ . Green curves illus- trate the model response, black curves indicate the mean value of the experimental data, and the shaded areas depict the standard deviation Fig. 5 Comparison of the data from the semi–confined compres- sion experiments with the model response (with parameters from the experimental set (iii), see Table  3) for a mode I, b mode II, and c mode III. Green curves illustrate the model response, black curves indicate the mean value of the experimental data, and the shaded areas depict the standard deviation 1510 F. Walter et al. 1 3 models for simplicity. To evaluate the overall performance of the model, Table 4 lists the error values with respect to the deformation states, the model’s ability to produce higher stresses in mode III than in mode I, and the number of material parameters used in the corresponding model (last column). As shown in the penultimate column, only the present approach, as well as the models of Blemker et al. (2005), Odegard et al. (2008), and Ehret et al. (2011) are basically able to generate higher stresses in mode III than I. Thereby, the stress responses of the model of Ehret et al. (2011) agree well with the experimental data. Figure 6 shows smaller deviations for axial compression and larger deviations for mode II in case of the semi–con- fined experiments, see Fig. 7b. In particular, the stresses at 40% compression are significantly higher than the max- imum stresses produced by the present model in mode II. Overall, both models can reproduce the stress–stretch behaviour of the semi–confined deformation state. How- ever, the present approach shows a slightly better overall performance. Although the model of Odegard et al. (2008) produces higher stresses in mode III than in mode I, the accuracy of the model prediction is limited. E.g., the stress response of the model fits very well for axial compres- sion in the case of 45◦ oriented fibres, but clearly deviates in case of the 0◦ fibre orientation, which becomes even clearer for 90◦ fibre orientation, see Fig. 6. Also for mode II the model shows clear deviations from the experimen- tal data, see Fig. 7b. Within the models able to generate higher stresses in mode III than in mode I, the approach by Blemker et  al. (2005) generates the highest error (mean value O(p) = 0.68 ). The model is able to reproduce mode II very well compared to the other models. For all other deformation states, rather larger deviations can be seen. In particular, it should be mentioned that the model does not provide any stress in the case of the axial com- pression test ( 0◦ fibre orientation), which is related to the model formulation. For more detailed information, the interested reader is referred to Blemker et al. (2005). Table 4 Error O(p) of the different material models when compared with the experimental data, model ability to produce higher stresses in mode III than in mode I (✓: distinction between mode III and mode I, ✗: no distinction between mode III and mode I), and the number of material parameters used in the corresponding model Error O(p) Reference 0◦ 45◦ 90◦ Mode I Mode II Mode III Mean PI < PIII No. model parameters Blemker et al. (2005) 0.89 0.24 0.61 1.21 0.45 0.67 0.68 ✓ 6 Calvo et al. (2010) 0.15 0.31 0.51 0.50 0.59 0.69 0.46 ✗ 5 Ehret et al. (2011) 0.33 0.37 0.16 0.36 1.74 0.39 0.56 ✓ 5 Johannson et al. (2000) 0.60 0.24 0.37 0.04 1.46 0.54 0.54 ✗ 4 Martins et al. (1998) 0.74 0.31 0.50 0.16 1.63 0.70 0.67 ✗ 4 Moerman et al. (2016) 0.23 0.28 0.46 0.60 0.59 0.77 0.49 ✗ 3 Odegard et al. (2008) 0.34 0.03 0.88 0.26 1.71 0.20 0.57 ✓ 4 Röhrle and Pullan (2007) 0.15 0.30 0.51 0.51 0.44 0.70 0.43 ✗ 6 Mean 0.43 0.26 0.50 0.46 1.08 0.58 Present model 0.17 0.29 0.20 0.38 1.58 0.23 0.48 ✓ 3 Fig. 6 Comparison of the data from the axial compression experiments with existing models for a 0◦ , b 45◦ , and c 90◦ . Coloured curves illustrate the model responses, black curves indicate the mean value of the experimental data, and the shaded areas depict the standard deviation 1511On a three–dimensional model for the description of the passive characteristics of skeletal… 1 3 While the stress responses, whose deviations from the experimental data in the form of error measures, are a char- acteristic, the number of parameters used by a model is also an important indicator of its quality. In the last column of the Table 4 the parameter numbers are given. While the model presented here uses 3 parameters, the approaches of Blemker et al. (2005), Odegard et al. (2008) and Ehret et al. (2011) use 6, 4, and 5 parameters respectively, which underlines the quality of the model presented here. The models of Röhrle and Pullan (2007), Moerman et al. (2016), and Calvo et al. (2010) produce the lowest mean error values. However, due to the Heaviside formulation to switch off the fibre contribu- tion for I4 < 1 , the models do not distinguish between mode III and mode I. This further results in the model responses of these three models being identical for modes I and III, as they only provide the tissue matrix response. Overall, the presented model shows a good performance between the prediction of the stress and stretch behaviour of the fibres for different fibre orientations under compressive loading. Also in comparison to the modelling concepts by Blemker et al. (2005), Odegard et al. (2008), and Ehret et al. (2011), the model shows good performance resulting in the lowest mean error of O(p) = 0.48. 6 Conclusions The experimental basis for this study is an earlier experi- mental campaign (Leichsenring et al. 2021), providing com- prehensive data on the muscle tissue of rabbit soleus mus- cles. Two orientation–dependent deformation states, namely axial and semi–confined compression, were realised. The orientation–dependent compression experiments show that the generated stresses increase from 0◦ over 45◦ up to 90◦ fibre orientation. In the case of the semi–confined compres- sion experiments, mode I shows the lowest and mode II the stiffest tissue response, with mode III in between. In this work we developed a phenomenologically motivated mate- rial formulation which can predict the stresses of the previ- ously described deformation states with good agreement. A special feature here is the fact that mode III is stiffer than mode I in the semi–confined experiments. The proposed model contributes to a better understanding of the load trans- fer mechanisms between fibres and ECM (Böl et al. 2014, 2015b). Furthermore, it could lead to improved predictions of models of human muscles under compression, which are used e.g. in rehabilitation engineering, impact biomechan- ics, or for the simulation of surgical interventions (Guc- cione et al. 2001; van Rooij et al. 2003; Linder–Ganz et al. 2009; Siebert et al. 2018). While the model developed here can describe this behaviour, not all models in the literature are able to do so. Overall, the model presented here shows a good performance and describes the deformation states shown here with only three model parameters in a good way. Author contributions FW Formal analysis, Writing—Original Draft, Writing— Review & Editing; RS Conceptualization, Formal analysis, Writing —Original Draft, Writing—Review & Editing; PM Formal analysis, Review & Editing; TS Funding acquisition, Writing— Review & Editing; MB Funding acquisition, Conceptualization, Visualiza- tion, Formal analysis, Writing— Original Draft, Writing—Review & Editing Funding Open Access funding enabled and organized by Projekt DEAL. This work was supported by the Deutsche Forschungsgemein- schaft (DFG) under Grants 127753742 and 316485047. Data availability Data and materials are available within this manuscript. Declarations Conflict of interests Not applicable Ethical approval See Sect. 2.1. Fig. 7 Comparison of the data from the semi–confined compres- sion experiments with existing models for a mode I, b mode II, and c mode III. Coloured curves illustrate the model responses, black curves indicate the mean value of the experimental data, and the shaded areas depict the standard deviation 1512 F. Walter et al. 1 3 Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/. References Ambrosi D, Pezzuto S (2012) Active stress versus active strain in mechanobiology: constitutive issues. J Elast 107(2):199–212 Baker JE, Brosseau C, Joel PB, Warshaw DM (2002) The biochemical kinetics underlying actin movement generated by one and many skeletal muscle myosin molecules. Biophys J 82(4):2134–2147 Balzani D, Neff P, Schröder J, Holzapfel GA (2006) A polyconvex framework for soft biological tissues. adjustment to experimental data. Int J Solids Struct 43(20):6052–6070 Bartoo ML, Linke WA, Pollack GH (1997) Basis of passive tension and stiffness in isolated rabbit myofibrils. Am J Physiol–Cell Physiol 273(1):C266–C276 Baskin RJ, Paolini PJ (1967) Volume change and pressure develop- ment in muscle during contraction. Am J Phys Legacy Content 213(4):1025–1030 Blemker SS, Pinsky PM, Delp SL (2005) A 3D model of muscle reveals the causes of nonuniform strains in the biceps brachii. J Biomech 38(4):657–665 Böl M, Kruse R, Ehret AE, Leichsenring K, Siebert T (2012) Compres- sive properties of passive skeletal muscle – The impact of precise sample geometry on parameter identification in inverse finite ele- ment analysis. J Biomech 45(15):2673–2679 Böl M, Leichsenring K, Weichert C, Sturmat M, Schenk P, Blickhan R, Siebert T (2013) Three–dimensional surface geometries of the rabbit soleus muscle during contraction: input for biomechani- cal modelling and its validation. Biomech Model Mechanobiol 12(6):1205–1220 Böl M, Ehret AE, Leichsenring K, Weichert C, Kruse R (2014) On the anisotropy of skeletal muscle tissue under compression. Acta Biomater 10(7):3225–3234 Böl M, Ehret AE, Leichsenring K, Ernst M (2015) Tissue–scale anisot- ropy and compressibility of tendon in semi–confined compression tests. J Biomech 48(6):1092–1098 Böl M, Leichsenring K, Ernst M, Wick C, Blickhan R, Siebert T (2015) Novel microstructural findings in M. plantaris and their impact during active and passive loading at the macro level. J Mech Behav Biomed Mater 51:25–39 Böl M, Leichsenring K, Ernst M, Ehret AE (2016) Long–term mechan- ical behaviour of skeletal muscle tissue in semi–confined com- pression experiments. J Mech Behav Biomed Mater 63:115–124 Böl M, Iyer R, Dittmann J, Garcés–Schröder M, Dietzel A (2019) Investigating the passive mechanical behaviour of skeletal muscle fibres: micromechanical experiments and Bayesian hierarchical modelling. Acta Biomater 92:277–289 Böl M, Iyer R, Garcés–Schröder M, Kohn S, Dietzel A (2020) Mech- ano–geometrical skeletal muscle fibre characterisation under cyclic and relaxation loading. J Mech Behav Biomed Mater 110:104001 Böl M, Kohn S, Leichsenring K, Morales–Orcajo E, Ehret AE (2022) On multiscale tension–compression asymmetry in skeletal muscle. Acta Biomater 96:1–2 Brown SHM, Carr JA, Ward SR, Lieber RL (2012) Passive mechanical properties of rat abdominal wall muscles suggest an important role of the extracellular connective tissue matrix. J Orthop Res 30(8):1321–1326 Calvo B, Ramírez A, Alonso A, Grasa J, Soteras F, Osta R, Muñoz MJ (2010) Passive nonlinear elastic behaviour of skeletal mus- cle: experimental results and model formulation. J Biomech 43(2):318–325 Chagnon G, Rebouah M, Favier D (2015) Hyperelastic energy densi- ties for soft biological tissues: a review. J Elast 120(2):129–160 Cooper KE, Cranston WI, Honour AJ (1965) Effects of intraventricular and intrahypothalamic injection of noradrenaline and 5–HT on body temperature in conscious rabbits. J Physiol 181(4):852–864 de Souza Leite F, Minozzo FC, Altman D, Rassier DE (2017) Microflu- idic perfusion shows intersarcomere dynamics within single skel- etal muscle myofibrils. Proc Natl Acad Sci 114(33):8794–8799 Ehret AE, Itskov M (2007) A polyconvex hyperelastic model for fiber– reinforced materials in application to soft tissues. J Mater Sci 42(21):8853–8863 Ehret AE, Böl M, Itskov M (2011) A continuum constitutive model for the active behaviour of skeletal muscle. J Mech Phys Solids 59(3):625–636 Federico S, Gasser TC (2010) Nonlinear elasticity of biological tissues with statistical fibre orientation. J R Soc Interface 7(47):955–966 Flory PJ (1961) Thermodynamic relations for high elastic materials. Trans Faraday Soc 57:829–838 Fridén J, Lieber RL (2003) Spastic muscle cells are shorter and stiffer than normal cells. Muscle Nerve Off J Am Assoc Electrodiagn Med 27(2):157–164 Gasser TC (2018) Vascular tissue biomechanics : constitutive modeling of the arterial wall. Elsevier, Heidelberg, pp 265–274 Gasser TC, Ogden RW, Holzapfel GA (2006) Hyperelastic modelling of arterial layers with distributed collagen fibre orientations. J R Soc Interface 3(6):15–35 Gillies AR, Lieber RL (2011) Structure and function of the skeletal muscle extracellular matrix. Muscle Nerve 44(3):318–331 Gras LL, Mitton D, Viot P, Laporte S (2012) Hyper–elastic properties of the human sternocleidomastoideus muscle in tension. J Mech Behav Biomed Mater 15:131–140 Guccione JM, Moonly SM, Wallace AW, Ratcliffe MB (2001) Residual stress produced by ventricular volume reduction surgery has lit- tle effect on ventricular function and mechanics: a finite element model study. J Thorac Cardiovasc Surg 122(3):592–599 Haeger RM, Rassier DE (2020) Force enhancement after stretch of iso- lated myofibrils is increased by sarcomere length non–uniformi- ties. Sci Rep 10(1):21590 Haeger R, de Souza Leite F, Rassier DE (2020) Sarcomere length non–uniformities dictate force production along the descend- ing limb of the force–length relation. Proc Royal Soc B Biol Sci 287(1937):20202133 Hashemi SS, Asgari M, Rasoulian A (2020) An experimental study of nonlinear rate–dependent behaviour of skeletal muscle to obtain passive mechanical properties. Proc Inst Mech Eng 234(6):590–602 Holzapfel GA, Ogden RW (2010) Constitutive modelling of arteries. Proc Royal Soc A Math Phys Eng Sci 466(2118):1551–1597 Holzapfel GA, Gasser TC, Ogden RW (2000) A new constitutive framework for arterial wall mechanics and a comparative study of material models. J Elast Phys Sci Solids 61(1):1–48 Humphrey JD, Strumpf RK, Yin FCP (1990) Determination of a consti- tutive relation for passive myocardium: i. a new functional form. J Biomech Eng 112(3):333–339 http://creativecommons.org/licenses/by/4.0/ http://creativecommons.org/licenses/by/4.0/ 1513On a three–dimensional model for the description of the passive characteristics of skeletal… 1 3 Johannson T, Meier P, Blickhan R (2000) A finite–element model for the mechanical analysis of skeletal muscles. J Theor Biol 206(1):131–149 Kohn S, Leichsenring K, Kuravi R, Ehret AE, Böl M (2021) Direct measurement of the direction–dependent mechanical behaviour of skeletal muscle extracellular matrix. Acta Biomater 122:249–262 Kuravi R, Leichsenring K, Trostorf R, Morales–Orcajo E, Böl M, Ehret AE (2021) Predicting muscle tissue response from calibrated com- ponent models and histology–based finite element models. J Mech Behav Biomed Mater 117:104375 Lagarias JC, Reeds JA, Wright MH, Wright PE (1998) Convergence properties of the Nelder–Mead simplex method in low dimen- sions. SIAM J Optim 9(1):112–147 Lee EJ, Jang HC, Koo KH, Kim HK, Lim JY (2020) Mechanical properties of single muscle fibers: understanding poor muscle quality in older adults with diabetes. Ann Geriatric Med Res 24(4):267–273 Leichsenring K, Viswanathan A, Kutschke S, Siebert T, Böl M (2021) Age–dependent mechanical and microstructural prop- erties of the rabbit soleus muscle. Acta Biomater 134:453–465 Li Y, Hessel AL, Unger A, Ing D, Recker J, Koser F, Freundt JK, Linke WA (2020) Graded titin cleavage progressively reduces tension and uncovers the source of A–band stability in contract- ing muscle. Elife 9:e64107 Lieber RL, Ward SR (2013) Cellular mechanisms of tissue fibrosis. 4. structural and functional consequences of skeletal muscle fibrosis. Am J Physiol Cell Physiol 305(3):C241–C252 Lieber RL, Runesson E, Einarsson F, Fridén J (2003) Inferior mechanical properties of spastic muscle bundles due to hyper- trophic but compromised extracellular matrix material. Muscle Nerve 28(4):464–471 Linder–Ganz E, Yarnitzky G, Yizhar Z, Siev–Ner I, Gefen A (2009) Real–time finite element monitoring of sub–dermal tissue stresses in individuals with spinal cord injury: toward preven- tion of pressure ulcers. Ann Biomed Eng 37(2):387–400 Linke WA, Ivemeyer M, Mundel P, Stockmeier MR, Kolmerer B (1998) Nature of PEVK–titin elasticity in skeletal muscle. Proc Natl Acad Sci 95(14):8052–8057 Marston S (2022) Force measurements from myofibril to filament. Front Physiol. https:// doi. org/ 10. 3389/ fphys. 2021. 817036 Martins JA, Pires EB, Salvado R, Dinis PB (1998) A numerical model of passive and active behavior of skeletal muscles. Com- put Methods Appl Mech Eng 151(3):419–433 Mathewson MA, Chambers HG, Girard PJ, Tenenhaus M, Schwartz AK, Lieber RL (2014) Stiff muscle fibers in calf muscles of patients with cerebral palsy lead to high passive muscle stiff- ness. J Orthop Res 32(12):1667–1674 Meyer GA, Lieber RL (2011) Elucidation of extracellular matrix mechanics from muscle fibers and fiber bundles. J Biomech 44(4):771–773 Moerman KM, Simms CK, Nagel T (2016) Control of tension– compression asymmetry in Ogden hyperelasticity with appli- cation to soft tissue modelling. J Mech Behav Biomed Mater 56:218–228 Mohammadkhah M, Murphy P, Simms CK (2016) The in vitro pas- sive elastic response of chicken pectoralis muscle to applied ten- sile and compressive deformation. J Mech Behav Biomed Mater 62:468–480 Morrow DA, Haut Donahue TL, Odegard GM, Kaufman KR (2010) Transversely isotropic tensile material properties of skeletal mus- cle tissue. J Mech Behav Biomed Mater 3(1):124–129 Mutungi G, Ranatunga KW (1996) The viscous, viscoelastic and elastic characteristics of resting fast and slow mammalian (rat) muscle fibres. J Physiol 496(3):827–836 Nelder JA, Mead R (1965) A simplex method for function minimiza- tion. Comput J 7(4):308–313 Nolan DR, Gower AL, Destrade M, Ogden RW, McGarry JP (2014) A robust anisotropic hyperelastic formulation for the modelling of soft tissue. J Mech Behav Biomed Mater 39:48–60 Noonan AM, Mazara N, Zwambag DP, Weersink E, Power GA, Brown SHM (2020) Age–related changes in human single muscle fibre passive elastic properties are sarcomere length dependent. Exp Gerontol 137:110968 Odegard GM, Haut Donahue TL, Morrow DA, Kaufman KR (2008) Constitutive modeling of skeletal muscle tissue with an explicit strain–energy function. J Biomech Eng 130(6):061017 Pertici I, Bongini L, Melli L, Bianchi G, Salvi L, Falorsi G, Squarci C, Bozó T, Cojoc D, Kellermayer MSZ, Lombardi V, Bianco P (2018) A myosin II nanomachine mimicking the striated muscle. Nat Commun 9(1):3532 Pietsch R, Wheatley BB, Haut Donahue TL, Gilbrech R, Prabhu R, Liao J, Williams LN (2014) Anisotropic compressive properties of passive porcine muscle tissue. J Biomech Eng 136(11):111003 Prado LG, Makarenko I, Andresen C, Krüger M, Opitz CA, Linke WA (2005) Isoform diversity of giant proteins in relation to pas- sive and active contractile properties of rabbit skeletal muscles. J Gen Physiol 126(5):461–480 Ramaswamy KS, Palmer ML, van der Meulen JH, Renoux A, Kostro- minova TY, Michele DE, Faulkner JA (2011) Lateral transmis- sion of force is impaired in skeletal muscles of dystrophic mice and very old rats. J Physiol 589(5):1195–1208 Röhrle O, Pullan AJ (2007) Three–dimensional finite element modelling of muscle forces during mastication. J Biomech 40(15):3363–3372 Sansour C (2008) On the physical assumptions underlying the volu- metric–isochoric split and the case of anisotropy. Eur J Mech A Solids 27(1):28–39 Scellini B, Piroddi N, Dente M, Vitale G, Pioner JM, Coppini R, Ferrantini C, Poggesi C, Tesi C (2021) Mavacamten has a differ- ential impact on force generation in myofibrils from rabbit psoas and human cardiac muscle. J Gen Physiol 153(7):e202012789 Shalabi N, Cornachione A, de Souza Leite F, Vengallatore S, Rassier DE (2017) Residual force enhancement is regulated by titin in skeletal and cardiac myofibrils. J Physiol 595(6):2085–2098 Siebert T, Leichsenring K, Rode C, Wick C, Stutzig N, Schubert H, Blickhan R, Böl M (2015) Three–dimensional muscle architec- ture and comprehensive dynamic properties of rabbit gastrocne- mius, plantaris and soleus: input for simulation studies. PLoS ONE 10(6):e0130985 Siebert T, Tomalka A, Stutzig N, Leichsenring K, Böl M (2017) Changes in three–dimensional muscle structure of rabbit gas- trocnemius, flexor digitorum longus, and tibialis anterior during growth. J Mech Behav Biomed Mater 74:507–519 Siebert T, Stutzig N, Rode C (2018) A hill–type muscle model expansion accounting for effects of varying transverse muscle load. J Biomech 66:57–62 Simms C, Kilroy H, Blackburn G, Takaza M (2016) The influence of physical dimension on apparent stress–strain behaviour of in vitro passive skeletal muscle samples. J Strain Anal Eng Des 52(1):3–11 Swist S, Unger A, Li Y, Vöge A, von Frieling–Salewsky M, Skärlén A, Cacciani N, Braun T, Larsson L, Linke WA (2020) Main- tenance of sarcomeric integrity in adult muscle cells crucially depends on Z–disc anchored titin. Nat Commun 11(1):4479 Takaza M, Moerman KM, Gindre J, Lyons G, Simms CK (2013) The anisotropic mechanical behaviour of passive skeletal muscle tissue subjected to large tensile strain. J Mech Behav Biomed Mater 17:209–220 Takaza M, Moerman KM, Simms CK (2013) Passive skeletal muscle response to impact loading: experimental testing and inverse modelling. J Mech Behav Biomed Mater 27:214–225 https://doi.org/10.3389/fphys.2021.817036 1514 F. Walter et al. 1 3 Tamura A, Hayashi S, Matsumoto T (2016) Effect of loading rate on viscoelastic properties and local mechanical heterogeneity of freshly isolated muscle fiber bundles subjected to uniaxial stretching. J Mech Med Biol 16(06):1650086 Tang CY, Zhang G, Tsui CP (2009) A 3D skeletal muscle model cou- pled with active contraction of muscle fibres and hyperelastic behaviour. J Biomech 42(7):865–872 Toursel T, Stevens L, Granzier H, Mounier Y (2002) Passive tension of rat skeletal soleus muscle fibers: effects of unloading condi- tions. J Appl Physiol 92(4):1465–1472 van Loocke M, Lyons CG, Simms CK (2006) A validated model of passive muscle in compression. J Biomech 39(16):2999–3009 van Rooij L, Bours R, van Hoof J, Mihm JJ, Ridella SA, Bass CR, Crandall JR (2003) The development, validation and application of a finite element upper extremity model subjected to air bag loading. Stapp Car Crash J 47:55–78 Ward SR, Winters TM, O’Connor SM, Lieber RL (2020) Non–linear scaling of passive mechanical properties in fibers, bundles, fasci- cles and whole rabbit muscles. Front Physiol 11:211 Wheatley BB (2020) Investigating passive muscle mechanics with biaxial stretch. Front Physiol. https:// doi. org/ 10. 3389/ fphys. 2020. 01021 Wood LK, Kayupov E, Gumucio JP, Mendias CL, Claflin DR, Brooks SV (2014) Intrinsic stiffness of extracellular matrix increases with age in skeletal muscles of mice. J Appl Physiol 117(4):363–369 Wu YN, Ren Y, Tsai LC, Gao F, Zhang LQ (2016) In vivo simultane- ous evaluations of sarcomere imaging and muscle fiber tension. J Biomech 49(5):797–801 Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. https://doi.org/10.3389/fphys.2020.01021 https://doi.org/10.3389/fphys.2020.01021 On a three–dimensional model for the description of the passive characteristics of skeletal muscle tissue Abstract 1 Introduction 2 Materials and methods 2.1 Ethical approval 2.2 Tissue sample dissection and processing 2.3 Experimental investigation on skeletal muscle tissue 2.3.1 Axial compression experiments 2.3.2 Semi–confined compression experiments 3 Continuum mechanical framework of passive muscle tissue modelling 3.1 Kinematics 3.2 Constitutive equation for fibrous material under compression 3.3 Optimisation scheme for material parameter identification 4 Results 4.1 Axial compression experiments 4.2 Semi–confined compression experiments 4.3 Optimised material parameters 5 Discussion 5.1 Mechanical muscle tissue characteristics 5.2 Comparison with existing muscle models 6 Conclusions References