- 1Tianjin Key Laboratory for Advanced Mechatronic System Design and Intelligent Control, School of Mechanical Engineering, Tianjin University of Technology, Tianjin, China
- 2National Demonstration Center for Experimental Mechanical and Electrical Engineering Education, Tianjin University of Technology, Tianjin, China
- 3Tianjin Stomatological Hospital, The Affiliated Stomatological Hospital of Nankai University, Tianjin Key Laboratory of Oral and Maxillofacial Function Reconstruction, Tianjin, China
- 4School of Chemical Engineering and Technology, Tianjin University, Tianjin, China
Creep deformation in cartilage can be observed under physiological loads in daily activities such as standing, single-leg lunge, the stance phase of gait. If not fully recovered in time, it may induce irreversible damage in cartilage and further lead to early osteoarthritis. In this study, 36 cruciform-shape samples in total from 18 bulls were employed to conduct the uniaxial and biaxial creep-recovery tests by using a biaxial cyclic testing system. Effects of stress level (σ = .5, 1.0, 1.5 MPa) and biaxial stress ratio (B = 0, .3, .5, 1.0) on creep-recovery behaviors of cartilage were characterized. And then, a viscoelastic constitutive model was employed to predict its creep-recovery behaviors. The results showed that the creep strain and its three components, namely instantaneous elastic strain, delayed elastic strain and viscous flow strain, increase with the increasing stress level or with the decreasing biaxial stress ratio. Compared with uniaxial creep-recovery, biaxial creep-recovery exhibits a smaller creep strain, a faster recovery rate of creep strain and a smaller residual strain. Besides, the built viscoelastic model can be used to describe the uniaxial creep-recovery behaviors of cartilage as a good correlation between the fitted results and test results is achieved. The findings are expected to provide new insights into understanding normal joint function and cartilage pathology.
Highlights
• Creep-recovery behaviors are explored at different stress levels and biaxial stress ratios.
• Creep strain and its components increase with stress level.
• A biaxial stress ratio of 1 imposes the smallest creep strain.
• Increasing biaxial stress ratio accelerates the recovery of creep strain and decreases the residual strain.
• The higher the biaxial stress ratio, the smaller the viscous flow strain and the delayed elastic strain.
1 Introduction
Articular cartilage has the functions to absorb shocks, transmit loads and sustain complex mechanical loading histories, thus maintaining dynamic mechanical equilibrium of knee joint in daily activities. As a kind of fluid saturated porous material, articular cartilage exhibits viscoelastic responses under mechanical loads. And creep is observed in cartilage when the mechanical load is held constant during standing, single-leg lunge (Hosseini et al., 2010; Choi et al., 2016; Qiu et al., 2016), the stance phase of gait (Liu et al., 2010; Li et al., 2014) and so on. Creep deformation under physiological loads may serve as a functional indicator for cartilage healthy. The majority of creep deformation in cartilage can recover timely. However, long term overloading may cause irreversible deformation, further induce damage in cartilage and lead to early osteoarthritis (Roos et al., 1998; Felson, 2013; Martel-Pelletier et al., 2016; Georgiev and Angelov, 2019; Peters et al., 2022). Studying the creep and creep-recovery behaviors of articular cartilage are helpful for understanding the normal joint function and cartilage pathology.
In vivo creep tests were usually conducted by applying 50%–150% body weight on the joint for 300–900 s (Herberhold et al., 1998; Herberhold et al., 1999; Hosseini et al., 2010; Choi et al., 2016; Uzuner et al., 2020), and the creep deformations of cartilage were measured in vivo by using imaging techniques like computed tomography (CT), magnetic resonance (MR), dual fluoroscopic imaging (DF). Herberhold et al. found that the creep deformation had high inter-individual variability (Herberhold et al., 1998; 1999). Hosseini et al. found the cartilage deformation of human tibiofemoral increases sharply in the first 20 s and almost held constant beyond 50 s when a constant full body weight was quickly applied and maintained for 300 s (Hosseini et al., 2010). Choi et al. studied the time-dependent creep behavior of tibial cartilage by applying a load of 75% body weight on limb and they reported that the creep deformation increased rapidly over the first minute and kept stable after 5 min (Choi et al., 2016). Uzuner et al. reported a substantial increase in cartilage deformation when joint loading increased from nil to 75% body weight and a continued small increase over time when the joint loading was held constant (Uzuner et al., 2020). In vivo experimental data help to understand the deformation characters of cartilage under physiological loads. However, technical limitations impede the accurate measurement of creep recovery in the physiological range and it non-etheless has significant difficulty to quantitatively analyze the creep-recovery deformation law of cartilage.
To cover the shortage of in vivo tests, in vitro creep tests are widely used to investigate the creep and creep-recovery behaviors of cartilage. Athanasiou et al. and Stolberg et al. employed the indentation test system to characterize the creep-recovery behavior of cartilage and they found that the initial creep rate is larger than the initial recovery rate and more than 95% of the creep deformation could recovered within 90 min (Athanasiou et al., 1991; Stolberg-Stolberg et al., 2018). Boettcher et al. found that the cartilage deformation increased with the increasing stress. (Boettcher et al., 2016). Cutclife and DeFrate validated the use of the recovery response for mechanical characterization of cartilage in a controlled, ex vivo environment (Cutcliffe and DeFrate, 2020). Reuter and Hurschler used a biphasic 3D-FE-based method to determine the biomechanical properties of equine articular cartilage (Reuter and Hurschler., 2021). In our previous study, the depth-dependent creep strain of cartilage under uniaxial compressive load was investigated by an optimized digital image correlation (DIC) technique (Si et al., 2022). And the creep deformation could accelerate the fatigue damage of cartilage in knee joint (Gao et al., 2019). It should be mentioned that these studies on the creep-recovery behaviors of articular cartilage were probed under uniaxial loading.
The cartilage with freely squeezing and stretching in vivo is subjected to biaxial or even multiaxial loads, and however literature provides limited data on these. Thus, it is of great significance to study the biaxial or multiaxial creep-recovery behaviors of cartilage. The goal of this study is to investigate the creep-recovery behaviors of articular cartilage under biaxial tensile load. The hypothesize is that changing the stress state from uniaxial loading to biaxial loading has a significant effect on creep deformation of cartilage. To verify the hypothesize, both the uniaxial and biaxial creep-recovery experiments were conducted by using a biaxial cyclic testing system. Effects of stress level and biaxial stress ratio on cartilage creep-recovery were characterized. Creep strain components were analyzed and compared at different loading conditions. Additionally, the creep-recovery behaviors of cartilage were predicted by mathematical models. This study is original because it provides the first insight into the effect of biaxial constraint on tensile creep behavior of cartilage. The results could help to explore intrinsic mechanical properties of articular cartilage and to prevent cartilage damage.
2 Materials and methods
2.1 Sample preparation
In total, 36 cruciform-shape samples from 18 bulls (around 8 months) were prepared for mechanical tests. The cruciform-shape samples were selected as test samples since they were widely used in biaxial tests to measure in plane mechanical properties of anisotropic material (Leotoing et al., 2013; Lamkanfi et al., 2015; Chen et al., 2022) and the sample preparation process in this study is similar to literature (Kamalanathan and Broom, 1993). As shown in Figure 1A , firstly, the cartilage with subchondral bone was harvested from the medial and lateral femoral trochleae of each femur. Secondly, the cruciform shape was marked on the cartilage with subchondral bone. Thirdly, the cartilage with subchondral bone was cut into the cruciform shape by a band saw. Finally, the subchondral bone was removed from cartilage by cutting with the band saw and abrading with an electric eraser (Si et al., 2022). Thus, the cruciform-shape cartilage sample with the thickness of 2.00 mm ± .2 mm was obtained. Figure 1B shows the specific structural dimensions of the cruciform specimen. It should be mentioned that since the curvature structure of cartilage induced by the orientation of collagen fibrils may affect the results, all the samples were cut from the same position of the femur. The pre-prepared samples were wrapped in gauze moistened with phosphate buffered saline and placed in a refrigerator at −20°C to preserve the collagen fibers intact of cartilage samples. Prior to tests, the samples were defrosted by being immersed in phosphate Buffered Saline (PBS, pH 7.4) for 12 h at room temperature.
2.2 Experimental description
Three kinds of mechanical tests, namely the uniaxial tensile tests, the uniaxial creep recovery tests and the biaxial creep recovery tests were conducted on the biaxial dynamic mechanical testing system (IPBF-300, CARE Measurement & Control Co., Ltd., China). For uniaxial tensile and uniaxial creep-recovery tests, only one side of the specimen was loaded, leaving the other side of the specimen freely shrinking; for biaxial creep-recovery tests, both loading axes were loaded and recovered at the same time. The uniaxial tensile tests were performed at different strain rates of .01%/s, .1%/s and 1%/s considering the limit of physiological strain rate of articular cartilage (Huang et al., 2003; Li et al., 2003; Li and Herzog, 2004). As reported that the deformation of cartilage under physiological loads is no larger than 30% with the average deformation value lower than 20%. Thus, the rate-dependent tensile tests were stopped when the strain reached 20% (Liu et al., 2010; Wang et al., 2015).
The stress on the cartilage was .535 MPa in double leg stance and .953 MPa in single leg stance (Qiu et al., 2016). In order to approximate the in vivo stress state, both uniaxial and biaxial creep-recovery tests were carried out at different stress levels σ of .5, 1.0 and 1.5 MPa. In addition, for the biaxial creep-recovery tests, the different biaxial stress ratios (denoted as B, indicates the ratio of stress in one loading direction to that in the other direction) of .3, .5 and 1.0 were applied. It can be seen that the biaxial stress ratio B) of uniaxial creep tests is 0 and B = 1 means that the cartilage is in an equi-biaxial stress state. For all the creep-recovery test conditions mentioned above, the creep time was 2,400 s and followed by the strain recovery (at 0 MPa) time of 3,600 s to make sure that most creep deformation of cartilage has recovered (Eckstein et al., 1999; Erisken et al., 2010).
All the tests were conducted at room temperature with a preload of 0.1 N being applied before tests to keep the samples align with the axis of loading. Besides, to prevent cartilage samples from drying (Erisken et al., 2008), an air humidifier was placed near the sample and a sprayer was employed to spray the PBS solution on the surface of sample every 5 min during tests (Si et al., 2022).
2.3 Image and data analysis
During mechanical tests, the load and displacement of the central area (10 mm × 10 mm) were recorded by the mechanical sensor and the digital image correlation (DIC) system equipped by the biaxial dynamic mechanical testing machine. The tensile stress of cartilage is the load divided by the cross-sectional area. The measurement of strain is similar with that reported in literatures (Kamalanathan and Broom, 1993; Lin et al., 2017; Lin et al., 2018; Zhang et al., 2019). As shown in Figure 2C, two pairs of thin lines were drawn in the middle local mark section of cruciform sample. A optical camera was applied to measure the variation of distance between fine lines. The distances between two pairs of marked lines were recorded as length
Only the strain in x-axis is recorded for uniaxial mechanical tests. For the creep-recovery tests, as shown in Figure 1D, the creep strain is divided into three components, namely instantaneous elastic strain
2.4 Statistical analysis
Considering the random error of experiments, three samples were repeated by the same person under the same loading condition, and test data were used in data analysis. The figures were plotted in Origin 2018 and the statistical analyses were performed using SPSS 22.0. R-squared value was used to indicate the agreement between the test result and its prediction. A one-way analysis of variance (ANOVAs) was employed to detect significant differences between different testing conditions. The statistical results were significant if p < .05.
3 Theoretical model
The viscoelastic constitutive model Eq. 2 introduced by Lou and Schapery is used to describe the unconfined mechanical behavior of cartilage (Lou and Schapery, 1971). The constitutive model under uniaxial loading can be expressed as
Where
Here
When the step-stress input
And the recovery strain can be expressed as
Where in Eq. (5)
The creep strain at t1 (just before unloading) can be described as follows in Eq. (7). At the moment after unloading the initial strain in recovery phase can be described as follows in Eq. (8).
Therefore, the strain jump value at
The transient component of linear viscoelastic creep compliance is expressed as
Where C is creep coefficient, n is creep index.
By substituting Eq. (10) into Eq. (4), we obtain the non-linear viscoelastic creep model.
By substituting Eq. (10) into Eq. (5), we obtain the expression of recovery strain:
Where in Eq. (12)
and
It is noted that the parameters in creep model Eq. 11 and recovery model (Eqs. 12–14 were obtained by fitting the test data with the non-linear least square method in the software origin 2018. Thus, the uniaxial creep-recovery strains of cartilage can be fitted by the built model.
4 Results
4.1 Tensile properties of articular cartilage at different strain rates
Figure 3A shows the stress-strain curves of articular cartilage at different strain rates. The stress-strain curves are not coincident, indicating that the tensile behavior of articular cartilage is strain rate-dependent. According to literature (Korhonen et al., 2002), and considering the error in data at the very beginning of the experiments, the Young’s modulus of cartilage was determined by calculating the slope of linear range of stress-strain curve when the strain value is 1%–5%. Effect of strain rate on Young’s modulus is shown in Figure 3B. Young’s modulus increases obviously with the increase of strain rate.
FIGURE 3. (A) Tensile stress-strain curves of articular cartilage; (B) Young’s modulus of articular cartilage.
4.2 Uniaxial creep-recovery behaviors of articular cartilage at different stress levels
The uniaxial creep-recovery curves of cartilage at different stress levels are shown in Figure 4. In the phase of creep, the creep strain at each stress level increases rapidly at first, and followed by the decrease of growth rate with creep time. The creep strain increases with the increase of stress level (p < .05). In the phase of recovery (σ = 0 MPa), the strain decreases rapidly at the beginning, and then the declining rate decreases with recovery time. When the recovery time is 2,400 s, the strain decreases to 45.1%, 50.3% and 53.5% of initial recovery strain at stress level of .5, 1.0, 1.5 MPa. When the recovery time is 3,600 s, the recovery strain decreases to 39.5%, 46.4% and 50.5% of initial recovery strain. Enough unloading time is beneficial for creep recovery. And longer recovery time is required for creep recovery at higher stress level.
Simultaneously, the creep strains and recovery strains at different stress levels have also been fitted by the theoretical creep-recovery model (Eq. (11) and Eq, 12. Figure 4 demonstrates the comparisons of test results and fitting results. There are good agreements between them (R > .8).
Based on the creep strain classification in Figure 2D, each component at different stress levels is compared in Figure 5. It is noted that the increase of stress level results in an increase in all strain components. As shown in Figure 5D, by comparing the percentage of each strain component, the contribution of each strain component at different stress levels was evaluated. Decrease in proportion of instantaneous elastic strain and increase in proportions of delayed elastic strain and viscous flow strain are observed with the increase of stress level.
FIGURE 5. (A–C) Instantaneous elastic strain, delayed elastic strain and viscous flow strain components of uniaxial creep strain, respectively; (D) percentage of each strain component at different stress.
4.3 Biaxial creep-recovery behaviors of articular cartilage at different stress levels
Biaxial creep-recovery tests with different stress levels were carried out at
FIGURE 7. (A–C) Instantaneous elastic strain, delayed elastic strain and viscous flow strain components of biaxial creep strain, respectively; (D) percentage of each biaxial strain component at different stress levels.
To further quantify the effect of biaxial stress level on creep strain, the strain ratio R was proposed and it is defined as uniaxial creep strain component divided by corresponding biaxial creep strain component. Figure 8 compares the strain ratios of instantaneous elastic strain, delayed elastic strain and viscous flow strain at different stress levels. It is found that all strain ratios at different stress levels are larger than 1, which means that each uniaxial creep strain component is larger than the corresponding biaxial creep strain component. At the same stress level, the strain ratio of viscous flow strain is larger than that of the other two components. The strain ratio increases with the increase of stress level by focusing on certain creep strain component. Compared with .5 MPa and 1.0 MPa, 1.5 MPa exhibits the most remarkable increase in all three creep strain components. Besides, the strain ratios of viscous flow strain and delayed elastic strain experience greater increase than that of instantaneous, regardless the stress level. Consequently, biaxial loading reduces the creep strain by retarding viscous flow strain and delayed elastic strain, especially at high stress of 1.5 MPa.
FIGURE 8. The strain ratio as a function of stress level for instantaneous elastic strain, delayed elastic strain and viscous flow strain.
4.4 Effect of biaxial stress ratio on biaxial creep-recovery behavior
Figure 9 shows the creep-recovery responses of cartilage at different biaxial stress ratios. The biaxial creep strain decreases with increase of biaxial stress ratio (p < .05), indicating that the stronger the biaxial constraint, the smaller the creep strain. Besides, the larger biaxial creep strain shows the smaller residual strain. As shown in Figures 10A–C, increasing biaxial stress ratio B from .3 to 1.0 brings a large reduction of each creep component, leading to the decreasing of creep and recovered strain. In particular, the delayed elastic strain decreases from 4.74 to 1.81 with biaxial stress ratio from B = .5 to B = 1.0, corresponding to a 72% reduction.
FIGURE 10. (A–C) Instantaneous elastic strain, delayed elastic strain and viscous flow strain components of biaxial creep strain, respectively.
In addition, the strain ratios of instantaneous elastic strain, delayed elastic strain and viscous flow strain at different stress ratios were compared in Figure 11. At the same stress level, the strain ratio of instantaneous elastic strain is smaller than that of the other two components. The strain ratio of instantaneous elastic strain rises slightly with the rise of biaxial stress ratio. By comparison, the strain ratios of delayed elastic strain and viscous flow strain show the more remarkable increase with the rise of biaxial stress ratio, especially when the biaxial stress ratio increase from B = .5 to B = 1.0. The results indicate that biaxial loading limits creep deformation of cartilage by inhibiting its delayed elastic strain and viscous flow strain, and the equi-biaxial stress state imposes the largest constriction.
FIGURE 11. The strain ratio as a function of biaxial stress ratio for instantaneous elastic strain, delayed elastic strain and viscous flow strain.
5 Discussions
In daily activities, articular cartilage is subjected to long time biaxial loads or multiaxial loads. The creep deformation is induced in cartilage under mechanical loads due to the alternation of unequal fluid exudation and imbibition within its matrix. Once unloaded, most of the creep deformation will gradually recover with time. Uniaxial creep-recovery of cartilage have been well described in previous studies (Gao et al., 2019; Cutcliffe and DeFrate, 2020). However, seldom has any study focused on biaxial or multiaxial creep-recovery of cartilage. Thus, we probed the creep-recovery behaviors of articular cartilage by comparing uniaxial and biaxial tensile loading results.
The creep strain evolution of cartilage consists of two parts, namely initial creep stage and the steady creep stage. In initial creep stage, the uniaxial and biaxial creep strain showed a similar evolution that the creep strain increased rapidly at the beginning, and gradually tended to be stable with time (Figure 4; Figure 6; Figure 9). The rapid increase of strain in initial creep stage for cartilage is due to the rapid outflow of interstitial fluid. When the unbound fluid in cartilage is squeezed out, the burden of supporting tensile load is gradually shifted to the collagen network. The cartilage becomes more resistant to deformation and the next creep stage starts. In steady creep stage, the increase rate of creep strain stays constant. The creep strain evolution observed in this study is similar with that reported by Hosseini et al. and Choi et al. as they found that the creep deformation in vivo kept stable after a rapid growth (Hosseini et al., 2010; Choi et al., 2016). However, it should be mentioned that there is no accurate cut-off point between the first creep stage and the steady creep stage since the creep deformation of cartilage varies with different physiological loads (Schoenbauer et al., 2015; Wang et al., 2015; Chan et al., 2016; Choi et al., 2016). Compared with uniaxial creep, biaxial creep exhibits smaller strain at the same stress level, especially at high stress level of 1.5 MPa. It is found that biaxial loading reduces the creep strain through retarding all creep strain components (Figure 5; Figure 7), especially through retarding delayed elastic strain and viscous flow strain (Figure 8). The smaller strain under biaxial creep is ascribed to the higher constraint imposed by the equi-biaxial stress state, as higher constraints can increase the stiffness of the cartilage and restrict the rapid outflow of interstitial fluid.
Generally, the deformation of articular cartilage can be fully recovered within a sufficient period of unloading time due to the viscoelastic properties of articular cartilage (Torzilli, 1984; Athanasiou et al., 1994). Eckstein et al. suggested at least 90 min of rest might be required to allow articular cartilage to fully recover according to the imaging based studies (Eckstein et al., 1999). Erisken et al. reported that 1,000 s was sufficient for the compressive stress to relax in bovine cartilage (Erisken et al., 2010). In this study, it is found that both uniaxial and biaxial strain have not been fully recovered with the residual strains left when the recovery time is 3,600 s. Compared with uniaxial loading, biaxial loading shows a smaller residual strain (Figure 4; Figure 6; Figure 9), due to the significant decrease in the viscous flow strain at equi-biaxial stress state (Figure 5C; Figure 7C). The residual strain, as viscous flow strain, is slow to recover. Thus, more time may be needed for the residual creep deformation to recover.
It is very likely for cartilage to experience non-uniform stress distribution in daily activities due to its curvature structure, thickness variation, diffusion kinetics and other factors (Kamalanathan and Broom, 1993). Kamalanathan and Broom reported that the tissue was significantly stiffer, i.e. more strain limiting, along the split-line direction than across it. They found that the tissue was significantly stiffer, i.e. more strain limiting, along the split-line direction than across it and the anisotropic strain-limiting response of the fibrillar array will result in fibrils acting more collectively along the less extensible direction than across it. Obviously, the stress state of cartilage under physiological load is not equi-biaxial. To better understand how biaxial constraint influences biaxial creep behavior, we studied the effect of biaxial stress ratio on creep-recovery behavior of cartilage (Figure 8). Both creep strain and recovered strain decrease with the increase of biaxial stress ratio due to the reason that the larger the biaxial stress ratio, the stronger the biaxial constraint. The equi-biaxial stress state (B = 1) imposes the largest biaxial constraint. Thus, when the biaxial stress ratio B is 1.0, the lowest creep and recovery strain is induced in cartilage. Increasing the biaxial stress ratio from .3 to 1.0 causes a decrease in all creep strain components, especially in delayed elastic strain (Figure 9; Figure 10). Furthermore, less proportion of delayed elastic strain results in less time to reach steady recovery state (Figure 8) since delayed elastic strain is slow to recover. In addition, the higher constraint was imposed by the larger biaxial stress ratio, which in turn increases the stiffness of cartilage (Lin et al., 2017; Zhang et al., 2019) resulting in the smaller creep strain and residual strain.
Viscoelastic (Leipzig and Athanasiou, 2005; Si et al., 2022), biphasic (Huang et al., 2003) and poroviscoelastic models (Wang et al., 2001) are employed to describe mechanical behaviors of cartilage in literature. Leipzig and Athanasiou used the viscoelastic model to describe the mechanical behaviors of cartilage cells (Leipzig and Athanasiou, 2005). Gao et al. employed the non-linear viscoelastic constitutive model to predict the depth-dependent creep behaviors of cartilage under compression (Gao et al., 2014). The viscoelastic constitutive model employed in this study is proposed by Lou and sharper (Lou and Schapery, 1971). It has the outstanding advantage of retaining the single time integral form even in the non-linear range. And the non-linear effects are expressed by means of stress-dependent material functions determined experimentally. It has been applied to predict the mechanical behaviors of cracked cartilage and exhibits good results (Si et al., 2022). In this study, the proposed viscoelastic model with the consideration of irrecoverable deformation was developed to describe the uniaxial creep-recovery behaviors of cartilage and a good consistence between test data and fitted data is obtained.
6 Conclusion
The study probed the uniaxial and biaxial creep-recovery behaviors of articular cartilage. Changing the stress state from uniaxial loading to biaxial loading, creep strain decreases significantly under the same stress level. Increasing biaxial stress ratio (from B = .3–1) causes notable decrease in creep strain. During the recovery process of creep, it takes less time for the biaxial strain than the uniaxial strain to reach the steady recovery state. The cartilage with higher stress level or smaller biaxial stress ratio shows larger residual strain, indicating that the cartilage damage may be accelerated at higher stress level or at smaller biaxial stress ratio. These results provide a new sight for cartilage creep and are significant to improve the durability of cartilage in daily activities.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors.
Author contributions
LG: Writing—original draft, Methodology. GL: Writing—original draft, Investigation. YT: Writing—review and editing. RL: Software. CZ: Methodology, Supervision. HG: Supervision. BZ: Supervision.
Acknowledgments
The authors gratefully acknowledge the financial support from National Natural Science Foundation of China (Nos. 32271371, 12072235), Tianjin Natural Science Foundation (Nos. 21JCQNJC01310 and 21JCYBJC00940) and Technical Innovation Guidance Project of Tianjin Science and Technology Bureau (No. 22YDTPJC00440).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Athanasiou, K. A., Agarwal, A., and Dzida, F. J. (1994). Comparative study of the intrinsic mechanical properties of the human acetabular and femoral head cartilage. J. Orthop. Res. 12, 340–349. doi:10.1002/jor.1100120306
Athanasiou, K. A., Rosenwasser, M. P., Buckwalter, J. A., Malinin, T. I., and Mow, V. C. (1991). Interspecies comparisons of in situ intrinsic mechanical properties of distal femoral cartilage. J. Orthop. Res. 9, 330–340. doi:10.1002/jor.1100090304
Boettcher, K., Kienle, S., Nachtsheim, J., Burgkart, R., Hugel, T., and Lieleg, O. (2016). The structure and mechanical properties of articular cartilage are highly resilient towards transient dehydration. Acta Biomater. 29, 180–187. doi:10.1016/j.actbio.2015.09.034
Chan, D. D., Cai, L. Y., Butz, K. D., Trippel, S. B., Nauman, E. A., and Neu, C. P. (2016). In vivo articular cartilage deformation: Noninvasive quantification of intratissue strain during joint contact in the human knee. Sci. Rep. 6, 19220. doi:10.1038/srep19220
Chen, J. X., Zhang, J. H., and Zhao, H. W. (2022). Quantifying alignment deviations for the in-plane biaxial test system via a shape-optimised cruciform specimen. Materials 15, 4949. doi:10.3390/ma15144949
Choi, J. H., Mcwalter, E. J., Datta, S., Mueller, K., Maier, A., Fahrig, R., et al. (2016). Tibial cartilage creep during weight bearing:in vivo 3D CT assessment. Osteoarthr. Cartil. 24, S104. doi:10.1016/j.joca.2016.01.210
Cutcliffe, H. C., and Defrate, L. E. (2020). Comparison of cartilage mechanical properties measured during creep and recovery. Sci. Rep. 10, 1547. doi:10.1038/s41598-020-58220-2
Eckstein, F., Tieschky, M., Faber, S., Englmeier, K. H., and Reiser, M. (1999). Functional analysis of articular cartilage deformation, recovery, and fluid flow following dynamic exercise in vivo. Anat. Embryology 200, 419–424. doi:10.1007/s004290050291
Erisken, C., Kalyon, D. M., and Wang, H. J. (2008). Functionally graded electrospun polycaprolactone and beta-tricalcium phosphate nanocomposites for tissue engineering applications. Biomaterials 29, 4065–4073. doi:10.1016/j.biomaterials.2008.06.022
Erisken, C., Kalyon, D. M., and Wang, H. J. (2010). Viscoelastic and biomechanical properties of osteochondral tissue constructs generated from graded polycaprolactone and beta-tricalcium phosphate composites. J. Biomechanical Engineering-Transactions Asme 132, 091013. doi:10.1115/1.4001884
Felson, D. T. (2013). Osteoarthritis as a disease of mechanics. Osteoarthr. Cartil. 21, 10–15. doi:10.1016/j.joca.2012.09.012
Gao, L. L., Liu, D. D., Gao, H., Lv, L. W., and Zhang, C. Q. (2019). Effects of creep and creep-recovery on ratcheting strain of articular cartilage under cyclic compression. Mater. Sci. Eng. C-Materials Biol. Appl. 94, 988–997. doi:10.1016/j.msec.2018.10.047
Gao, L. L., Zhang, C. Q., Gao, H., Liu, Z. D., and Xiao, P. P. (2014). Depth and rate dependent mechanical behaviors for articular cartilage: Experiments and theoretical predictions. Mater. Sci. Eng. C-Materials Biol. Appl. 38, 244–251. doi:10.1016/j.msec.2014.02.009
Georgiev, T., and Angelov, A. K. (2019). Modifiable risk factors in knee osteoarthritis: Treatment implications. Rheumatol. Int. 39, 1145–1157. doi:10.1007/s00296-019-04290-z
Herberhold, C., Faber, S., Stammberger, T., Steinlechner, M., Putz, R., Englmeier, K. H., et al. (1999). In situ measurement of articular cartilage deformation in intact femoropatellar joints under static loading. J. Biomechanics 32, 1287–1295. doi:10.1016/s0021-9290(99)00130-x
Herberhold, C., Stammberger, T., Faber, S., Putz, R., Englmeier, K. H., Reiser, M., et al. (1998). An MR-based technique for quantifying the deformation of articular cartilage during mechanical loading in an intact cadaver joint. Magnetic Reson. Med. 39, 843–850. doi:10.1002/mrm.1910390522
Hosseini, A., Van De Velde, S. K., Kozanek, M., Gill, T. J., Grodzinsky, A. J., Rubash, H. E., et al. (2010). In-vivo time-dependent articular cartilage contact behavior of the tibiofemoral joint. Osteoarthr. Cartil. 18, 909–916. doi:10.1016/j.joca.2010.04.011
Huang, C. Y., Soltz, M. A., Kopacz, M., Mow, V. C., and Ateshian, G. A. (2003). Experimental verification of the roles of intrinsic matrix viscoelasticity and tension-compression nonlinearity in the biphasic response of cartilage. J. Biomechanical Engineering-Transactions Asme 125, 84–93. doi:10.1115/1.1531656
Kamalanathan, S., and Broom, N. D. (1993). The biomechanical ambiguity of the articular surface. J. Anat. 183, 567–578.
Korhonen, R. K., Laasanen, M. S., Toyras, J., Rieppo, J., Hirvonen, J., Helminen, H. J., et al. (2002). Comparison of the equilibrium response of articular cartilage in unconfined compression, confined compression and indentation. J. Biomechanics 35, 903–909. doi:10.1016/s0021-9290(02)00052-0
Lamkanfi, E., Van Paepegem, W., and Degrieck, J. (2015). Shape optimization of a cruciform geometry for biaxial testing of polymers. Polym. Test. 41, 7–16. doi:10.1016/j.polymertesting.2014.09.016
Leipzig, N. D., and Athanasiou, K. A. (2005). Unconfined creep compression of chondrocytes. J. Biomechanics 38, 77–85. doi:10.1016/j.jbiomech.2004.03.013
Leotoing, L., Guines, D., Zidane, I., and Ragneau, E. (2013). Cruciform shape benefits for experimental and numerical evaluation of sheet metal formability. J. Mater. Process. Technol. 213, 856–863. doi:10.1016/j.jmatprotec.2012.12.013
Li, J. S., Tsai, T. Y., Wang, S. B., Li, P. Y., Kwon, Y. M., Freiberg, A., et al. (2014). Prediction of in vivo knee joint kinematics using a combined dual fluoroscopy imaging and statistical shape modeling technique. J. Biomechanical Engineering-Transactions Asme 136, 124503. doi:10.1115/1.4028819
Li, L. P., Buschmann, M. D., and Shirazi-Adl, A. (2003). Strain-rate dependent stiffness of articular cartilage in unconfined compression. J. Biomechanical Engineering-Transactions Asme 125, 161–168. doi:10.1115/1.1560142
Li, L. P., and Herzog, W. (2004). Strain-rate dependence of cartilage stiffness in unconfined compression: The role of fibril reinforcement versus tissue volume change in fluid pressurization. J. Biomechanics 37, 375–382. doi:10.1016/s0021-9290(03)00263-x
Lin, Q., Shi, S. W., Wang, L., Chen, S., Chen, X., and Chen, G. (2017). In-plane biaxial cyclic mechanical behavior of proton exchange membranes. J. Power Sources 360, 495–503. doi:10.1016/j.jpowsour.2017.06.040
Lin, Q., Shi, S. W., Wang, L., Chen, X., and Chen, G. (2018). Biaxial fatigue crack propagation behavior of perfluorosulfonic-acid membranes. J. Power Sources 384, 58–65. doi:10.1016/j.jpowsour.2018.02.002
Liu, F., Kozanek, M., Hosseini, A., Van De Velde, S. K., Gill, T. J., Rubash, H. E., et al. (2010). In vivo tibiofemoral cartilage deformation during the stance phase of gait. J. Biomechanics 43, 658–665. doi:10.1016/j.jbiomech.2009.10.028
Lou, Y. C., and Schapery, R. A. (1971). Viscoelastic characterization of a nonlinear fiber-reinforced plastic. J. Compos. Mater. 5, 208–234. doi:10.1177/002199837100500206
Martel-Pelletier, J., Barr, A. J., Cicuttini, F. M., Conaghan, P. G., Cooper, C., Goldring, M. B., et al. (2016). Osteoarthritis. Nat. Rev. Dis. Prim. 2, 16072. doi:10.1038/nrdp.2016.72
Peters, A. E., Geraghty, B., Bates, K. T., Akhtar, R., Readioff, R., and Comerford, E. (2022). Ligament mechanics of ageing and osteoarthritic human knees. Front. Bioeng. Biotechnol. 10, 954837. doi:10.3389/fbioe.2022.954837
Qiu, L. L., Ma, X. M., Gao, L. L., Men, Y. T., and Zhang, C. Q. (2016). Analysis of the mechanical state of the human joint with defect cartilage in standing. J. Mech. Med. Biol. 16, 1640021. doi:10.1142/s0219519416400212
Reuter, T., and Hurschler, C. (2021). Comparison of biphasic material properties of equine articular cartilage estimated from stress relaxation and creep indentation tests. Curr. Dir. Biomed. Eng. 7, 363–366. doi:10.1515/cdbme-2021-2092
Roos, E. M., Roos, H. P., Ekdahl, C., and Lohmander, L. S. (1998). Knee injury and osteoarthritis outcome score (KOOS) - validation of a Swedish version. Scand. J. Med. Sci. Sports 8, 439–448. doi:10.1111/j.1600-0838.1998.tb00465.x
Schoenbauer, E., Szomolanyi, P., Shiomi, T., Juras, V., Zbyn, S., Zak, L., et al. (2015). Cartilage evaluation with biochemical MR imaging using in vivo Knee compression at 3 T-comparison of patients after cartilage repair with healthy volunteers. J. Biomechanics 48, 3349–3355. doi:10.1016/j.jbiomech.2015.06.016
Si, Y. P., Tan, Y. S., Gao, L. L., Li, R. X., Zhang, C. Q., Gao, H., et al. (2022). Mechanical properties of cracked articular cartilage under uniaxial creep and cyclic tensile loading. J. Biomechanics 134, 110988. doi:10.1016/j.jbiomech.2022.110988
Stolberg-Stolberg, J., Foehr, P., Pflieger, I., Kuntz, L., Von Deimling, C., Obermeier, A., et al. (2018). Analysis of cartilage creep recovery using a highly dynamic closed-loop test system. J. Bionic Eng. 15, 1057–1066. doi:10.1007/s42235-018-0093-x
Torzilli, P. A. (1984). Mechanical response of articular-cartilage to an oscillating load. Mech. Res. Commun. 11, 75–82. doi:10.1016/0093-6413(84)90100-9
Uzuner, S., Li, L. P., Kucuk, S., and Memisoglu, K. (2020). Changes in knee joint mechanics after medial meniscectomy determined with a poromechanical model. J. Biomechanical Engineering-Transactions Asme 142, 101006. doi:10.1115/1.4047343
Wang, C. C. B., Hung, C. T., and Mow, V. C. (2001). An analysis of the effects of depth-dependent aggregate modulus on articular cartilage stress-relaxation behavior in compression. J. Biomechanics 34, 75–84. doi:10.1016/s0021-9290(00)00137-8
Wang, H. S., Koff, M. F., Potter, H. G., Warren, R. F., Rodeo, S. A., and Maher, S. A. (2015). An MRI-compatible loading device to assess knee joint cartilage deformation: Effect of preloading and inter-test repeatability. J. Biomechanics 48, 2934–2940. doi:10.1016/j.jbiomech.2015.08.006
Keywords: articular cartilage, biaxial loading, stress ratio, creep-recovery, strain ratio
Citation: Gao L, Liu G, Tan Y, Li R, Zhang C, Gao H and Zhao B (2023) Creep-recovery behaviors of articular cartilage under uniaxial and biaxial tensile loadings. Front. Bioeng. Biotechnol. 10:1085062. doi: 10.3389/fbioe.2022.1085062
Received: 31 October 2022; Accepted: 28 December 2022;
Published: 10 January 2023.
Edited by:
Massimiliano Zingales, University of Palermo, ItalyReviewed by:
Cevat Erisken, Nazarbayev University, KazakhstanMatteo Berni, Rizzoli Orthopedic Institute (IRCCS), Italy
Copyright © 2023 Gao, Liu, Tan, Li, Zhang, Gao and Zhao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yansong Tan, dG9tb3Jyb3cyMDEyQDE2My5jb20=; Chunqiu Zhang, emhhbmdfY2h1bnFpdUAxMjYuY29t