The role of flow-independent viscoelasticity in the biphasic tensile and compressive responses of articular cartilage

A long-standing challenge in the biomechanics of connective tissues (e.g., articular cartilage, ligament, tendon) has been the reported disparities between their tensile and compressive properties. In general, the intrinsic tensile properties of the solid matrices of these tissues are dictated by the collagen content and microstructural architecture, and the intrinsic compressive properties are dictated by their proteoglycan content and molecular organization as well as water content. These distinct materials give rise to a pronounced and experimentally well-documented nonlinear tension-compression stress-strain responses, as well as biphasic or intrinsic extracellular matrix viscoelastic responses. While many constitutive models of articular cartilage have captured one or more of these experimental responses, no single constitutive law has successfully described the uniaxial tensile and compressive responses of cartilage within the same framework. The objective of this study was to combine two previously proposed extensions of the biphasic theory of Mow et al. [1980, ASME J. Biomech. Eng., 102, pp. 73-84] to incorporate tension-compression nonlinearity as well as intrinsic viscoelasticity of the solid matrix of cartilage. The biphasic-conewise linear elastic model proposed by Soltz and Ateshian [2000, ASME J. Biomech. Eng., 122, pp. 576-586] and based on the bimodular stress-strain constitutive law introduced by Curnier et al. [1995, J. Elasticity, 37, pp. 1-38], as well as the biphasic poroviscoelastic model of Mak [1986, ASME J. Biomech. Eng., 108, pp. 123-130], which employs the quasi-linear viscoelastic model of Fung [1981, Biomechanics: Mechanical Properties of Living Tissues, Springer-Verlag, New York], were combined in a single model to analyze the response of cartilage to standard testing configurations. Results were compared to experimental data from the literature and it was found that a simultaneous prediction of compression and tension experiments of articular cartilage, under stress-relaxation and dynamic loading, can be achieved when properly taking into account both flow-dependent and flow-independent viscoelasticity effects, as well as tension-compression nonlinearity.     

Effects of tension-compression nonlinearity on solute transport in charged hydrated fibrous tissues under dynamic unconfined compression

Cartilage is a charged hydrated fibrous tissue exhibiting a high degree of tension-compression nonlinearity (i.e., tissue anisotropy). The effect of tension-compression nonlinearity on solute transport has not been investigated in cartilaginous tissue under dynamic loading conditions. In this study, a new model was developed based on the mechano-electrochemical mixture model [Yao and Gu, 2007, J. Biomech. Model Mechanobiol., 6, pp. 63-72, Lai et al., 1991, J. Biomech. Eng., 113, pp. 245-258], and conewise linear elasticity model [Soltz and Ateshian, 2000, J. Biomech. Eng., 122, pp. 576-586; Curnier et al., 1995, J. Elasticity, 37, pp. 1-38]. The solute desorption in cartilage under unconfined dynamic compression was investigated numerically using this new model. Analyses and results demonstrated that a high degree of tissue tension-compression nonlinearity could enhance the transport of large solutes considerably in the cartilage sample under dynamic unconfined compression, whereas it had little effect on the transport of small solutes (at 5% dynamic strain level). The loading-induced convection is an important mechanism for enhancing the transport of large solutes in the cartilage sample with tension-compression nonlinearity. The dynamic compression also promoted diffusion of large solutes in both tissues with and without tension-compression nonlinearity. These findings provide a new insight into the mechanisms of solute transport in hydrated, fibrous soft tissues.     

Dynamic response of immature bovine articular cartilage in tension and compression, and nonlinear viscoelastic modeling of the tensile response

Very limited information is currently available on the constitutive modeling of the tensile response of articular cartilage and its dynamic modulus at various loading frequencies. The objectives of this study were to (1) formulate and experimentally validate a constitutive model for the intrinsic viscoelasticity of cartilage in tension, (2) confirm the hypothesis that energy dissipation in tension is less than in compression at various loading frequencies, and (3) test the hypothesis that the dynamic modulus of cartilage in unconfined compression is dependent upon the dynamic tensile modulus. Experiment 1: Immature bovine articular cartilage samples were tested in tensile stress relaxation and cyclical loading. A proposed reduced relaxation function was fitted to the stress-relaxation response and the resulting material coefficients were used to predict the response to cyclical loading. Adjoining tissue samples were tested in unconfined compression stress relaxation and cyclical loading. Experiment 2: Tensile stress relaxation experiments were performed at varying strains to explore the strain-dependence of the viscoelastic response. The proposed relaxation function successfully fit the experimental tensile stress-relaxation response, with R2 = 0.970+/-0.019 at 1% strain and R2 = 0.992+/-0.007 at 2% strain. The predicted cyclical response agreed well with experimental measurements, particularly for the dynamic modulus at various frequencies. The relaxation function, measured from 2% to 10% strain, was found to be strain dependent, indicating that cartilage is nonlinearly viscoelastic in tension. Under dynamic loading, the tensile modulus at 10 Hz was approximately 2.3 times the value of the equilibrium modulus. In contrast, the dynamic stiffening ratio in unconfined compression was approximately 24. The energy dissipation in tension was found to be significantly smaller than in compression (dynamic phase angle of 16.7+/-7.4 deg versus 53.5+/-12.8 deg at 10(-3) Hz). A very strong linear correlation was observed between the dynamic tensile and dynamic compressive moduli at various frequencies (R2 = 0.908+/-0.100). The tensile response of cartilage is nonlinearly viscoelastic, with the relaxation response varying with strain. A proposed constitutive relation for the tensile response was successfully validated. The frequency response of the tensile modulus of cartilage was reported for the first time. Results emphasize that fluid-flow dependent viscoelasticity dominates the compressive response of cartilage, whereas intrinsic solid matrix viscoelasticity dominates the tensile response. Yet the dynamic compressive modulus of cartilage is critically dependent upon elevated values of the dynamic tensile modulus.     

Acoustic properties of articular cartilage under mechanical stress

Mechano-acoustic and elastographic techniques may provide quantitative means for the in vivo diagnostics of articular cartilage. These techniques assume that sound speed does not change during tissue loading. As articular cartilage shows volumetric changes during compression, acoustic properties of cartilage may change affecting the validity of mechano-acoustic measurements. In this study, we examined the ultrasound propagation through human, bovine and porcine articular cartilage during stress-relaxation in unconfined compression. The time of flight (TOF) technique with known cartilage thickness (true sound speed) as well as in situ calibration method [Suh, Youn, Fu, J. Biomech. 34 (2001), 1347-1353] were used for the determination of sound speed. Ultrasound speed and attenuation decreased in articular cartilage during ramp compression, but returned towards the level of original values during relaxation. Variations in ultrasound speed induced an error in strain and compressive moduli provided that constant ultrasound speed and time-of-flight data was used to determine the tissue thickness. Highest errors in strain (-11.8 +/- 12.0%) and dynamic modulus (15.4 +/- 17.9%) were recorded in bovine cartilage. TOF and in situ calibration methods yielded different results for changes in sound speed during compression. We speculate that the variations in acoustic properties in loaded cartilage are related to rearrangement of the interstitial matrix, especially to that of collagen fibers. In human cartilage the changes, are, however relatively small and, according to the numerical simulations, mechano-acoustic techniques that assume constant acoustic properties for the cartilage will not be significantly impaired by this phenomenon.     

A Conewise Linear Elasticity mixture model for the analysis of tension-compression nonlinearity in articular cartilage

A biphasic mixture model is developed that can account for the observed tension-compression nonlinearity of cartilage by employing the continuum-based Conewise Linear Elasticity (CLE) model of Curnier et al. (J. Elasticity, 37, 1-38, 1995) to describe the solid phase of the mixture. In this first investigation, the orthotropic octantwise linear elasticity model was reduced to the more specialized case of cubic symmetry, to reduce the number of elastic constants from twelve to four. Confined and unconfined compression stress-relaxation, and torsional shear testing were performed on each of nine bovine humeral head articular cartilage cylindrical plugs from 6 month old calves. Using the CLE model with cubic symmetry, the aggregate modulus in compression and axial permeability were obtained from confined compression (H-A = 0.64 +/- 0.22 MPa, k2 = 3.62 +/- 0.97 x 10(-16) m4/N.s, r2 = 0.95 +/- 0.03), the tensile modulus, compressive Poisson ratio, and radial permeability were obtained from unconfined compression (E+Y = 12.75 +/- 1.56 MPa, v- = 0.03 +/- 0.01, kr = 6.06 +/- 2.10 x 10(-16) m4/N.s, r2 = 0.99 +/- 0.00), and the shear modulus was obtained from torsional shear (mu = 0.17 +/- 0.06 MPa). The model was also employed to predict the interstitial fluid pressure successfully at the center of the cartilage plug in unconfined compression (r2 = 0.98 +/- 0.01). The results of this study demonstrate that the integration of the CLE model with the biphasic mixture theory can provide a model of cartilage that can successfully curve-fit three distinct testing configurations while producing material parameters consistent with previous reports in the literature.     

The correspondence between equilibrium biphasic and triphasic material properties in mixture models of articular cartilage

Mixture models have been successfully used to describe the response of articular cartilage to various loading conditions. Mow et al. (J. Biomech. Eng. 102 (1980) 73) formulated a biphasic mixture model of articular cartilage where the collagen-proteoglycan matrix is modeled as an intrinsically incompressible porous-permeable solid matrix, and the interstitial fluid is modeled as an incompressible fluid. Lai et al. (J. Biomech. Eng. 113 (1991) 245) proposed a triphasic model of articular cartilage as an extension of their biphasic theory, where negatively charged proteoglycans are modeled to be fixed to the solid matrix, and monovalent ions in the interstitial fluid are modeled as additional fluid phases. Since both models co-exist in the cartilage literature, it is useful to show how the measured properties of articular cartilage (the confined and unconfined compressive and tensile moduli, the compressive and tensile Poisson's ratios, and the shear modulus) relate to both theories. In this study, closed-form expressions are presented that relate biphasic and triphasic material properties in tension, compression and shear. These expressions are then compared to experimental findings in the literature to provide greater insight into the measured properties of articular cartilage as a function of bathing solutions salt concentrations and proteoglycan fixed-charge density.     

New insight into deformation-dependent hydraulic permeability of gels and cartilage,and dynamic behavior of agarose gels in confined compression

Equilibrium, creep, and dynamic behaviors of agarose gels (2.0-14.8%) in confined compression were investigated in this study. The hydraulic permeabilities of gels were determined by curve-fitting creep data to the biphasic model (J. Biomech. Eng. 102 (1980) 73) and found to be similar in value to those published in the literature (AIChE J. 42 (1996) 1220). A new relationship between intrinsic permeability and volume fraction of water was found for agarose gel, capable of predicting deformation-dependent permeabilities of bovine articular cartilage and 2% agarose gel published in literature. This relationship is accurate for gels and cartilage over a wide range of permeabilities (four orders of magnitude variation). The dynamic stiffness of the gels increases with gel concentration and loading frequency (0.01-1.0Hz). The increase in dynamic stiffness with loading frequency is less pronounced for gels with higher concentrations. The results of this study provide a new insight into deformation-dependent permeability behavior of agarose gel and cartilage, and are important for understanding biological responses of cells to interstitial fluid flow in gel or in cartilage under dynamic mechanical loading.     

Fibril reinforced poroelastic model predicts specifically mechanical behavior of normal,proteoglycan depleted and collagen degraded articular cartilage

Degradation of collagen network and proteoglycan (PG) macromolecules are signs of articular cartilage degeneration. These changes impair cartilage mechanical function. Effects of collagen degradation and PG depletion on the time-dependent mechanical behavior of cartilage are different. In this study, numerical analyses, which take the compression-tension nonlinearity of the tissue into account, were carried out using a fibril reinforced poroelastic finite element model. The study aimed at improving our understanding of the stress-relaxation behavior of normal and degenerated cartilage in unconfined compression. PG and collagen degradations were simulated by decreasing the Young's modulus of the drained porous (nonfibrillar) matrix and the fibril network, respectively. Numerical analyses were compared to results from experimental tests with chondroitinase ABC (PG depletion) or collagenase (collagen degradation) digested samples. Fibril reinforced poroelastic model predicted the experimental behavior of cartilage after chondroitinase ABC digestion by a major decrease of the drained porous matrix modulus (-64+/-28%) and a minor decrease of the fibril network modulus (-11+/-9%). After collagenase digestion, in contrast, the numerical analyses predicted the experimental behavior of cartilage by a major decrease of the fibril network modulus (-69+/-5%) and a decrease of the drained porous matrix modulus (-44+/-18%). The reduction of the drained porous matrix modulus after collagenase digestion was consistent with the microscopically observed secondary PG loss from the tissue. The present results indicate that the fibril reinforced poroelastic model is able to predict specifically characteristic alterations in the stress-relaxation behavior of cartilage after enzymatic modifications of the tissue. We conclude that the compression-tension nonlinearity of the tissue is needed to capture realistically the mechanical behavior of normal and degenerated articular cartilage.     

On the electric potentials inside a charged soft hydrated biological tissue: streaming potential versus diffusion potential

The main objective of this study is to determine the nature of electric fields inside articular cartilage while accounting for the effects of both streaming potential and diffusion potential. Specifically, we solve two tissue mechano-electrochemical problems using the triphasic theories developed by Lai et al. (1991, ASME J. Biomech Eng., 113, pp. 245-258) and Gu et al. (1998, ASME J. Biomech. Eng., 120, pp. 169-180) (1) the steady one-dimensional permeation problem; and (2) the transient one-dimensional ramped-displacement, confined-compression, stress-relaxation problem (both in an open circuit condition) so as to be able to calculate the compressive strain, the electric potential, and the fixed charged density (FCD) inside cartilage. Our calculations show that in these two technically important problems, the diffusion potential effects compete against the flow-induced kinetic effects (streaming potential) for dominance of the electric potential inside the tissue. For softer tissues of similar FCD (i.e., lower aggregate modulus), the diffusion potential effects are enhanced when the tissue is being compressed (i.e., increasing its FCD in a nonuniform manner) either by direct compression or by drag-induced compaction; indeed, the diffusion potential effect may dominate over the streaming potential effect. The polarity of the electric potential field is in the same direction of interstitial fluid flow when streaming potential dominates, and in the opposite direction of fluid flow when diffusion potential dominates. For physiologically realistic articular cartilage material parameters, the polarity of electric potential across the tissue on the outside (surface to surface) may be opposite to the polarity across the tissue on the inside (surface to surface). Since the electromechanical signals that chondrocytes perceive in situ are the stresses, strains, pressures and the electric field generated inside the extracellular matrix when the tissue is deformed, the results from this study offer new challenges for the understanding of possible mechanisms that control chondrocyte biosyntheses.     

Cartilage interstitial fluid load support in unconfined compression

Under physiological conditions of loading, articular cartilage is subjected to both compressive strains, normal to the articular surface, and tensile strains, tangential to the articular surface. Previous studies have shown that articular cartilage exhibits a much higher modulus in tension than in compression, and theoretical analyses have suggested that this tension–compression nonlinearity enhances the magnitude of interstitial fluid pressurization during loading in unconfined compression, above a theoretical threshold of 33% of the average applied stress. The first hypothesis of this experimental study is that the peak fluid load support in unconfined compression is significantly greater than the 33% theoretical limit predicted for porous permeable tissues modeled with equal moduli in tension and compression. The second hypothesis is that the peak fluid load support is higher at the articular surface side of the tissue samples than near the deep zone, because the disparity between the tensile and compressive moduli is greater at the surface zone. Ten human cartilage samples from six patellofemoral joints, and 10 bovine cartilage specimens from three calf patellofemoral joints were tested in unconfined compression. The peak fluid load support was measured at 79±11% and 69±15% at the articular surface and deep zone of human cartilage, respectively, and at 94±4% and 71±8% at the articular surface and deep zone of bovine calf cartilage, respectively. Statistical analyses confirmed both hypotheses of this study. These experimental results suggest that the tension–compression nonlinearity of cartilage is an essential functional property of the tissue which makes interstitial fluid pressurization the dominant mechanism of load support in articular cartilage.     

Stress-relaxation of human patellar articular cartilage in unconfined compression: prediction of mechanical response by tissue composition and structure

Mechanical properties of articular cartilage are controlled by tissue composition and structure. Cartilage function is sensitively altered during tissue degeneration, in osteoarthritis (OA). However, mechanical properties of the tissue cannot be determined non-invasively. In the present study, we evaluate the feasibility to predict, without mechanical testing, the stress-relaxation response of human articular cartilage under unconfined compression. This is carried out by combining microscopic and biochemical analyses with composition-based mathematical modeling. Cartilage samples from five cadaver patellae were mechanically tested under unconfined compression. Depth-dependent collagen content and fibril orientation, as well as proteoglycan and water content were derived by combining Fourier transform infrared imaging, biochemical analyses and polarized light microscopy. Finite element models were constructed for each sample in unconfined compression geometry. First, composition-based fibril-reinforced poroviscoelastic swelling models, including composition and structure obtained from microscopical and biochemical analyses were fitted to experimental stress-relaxation responses of three samples. Subsequently, optimized values of model constants, as well as compositional and structural parameters were implemented in the models of two additional samples to validate the optimization. Theoretical stress-relaxation curves agreed with the experimental tests (R=0.95-0.99). Using the optimized values of mechanical parameters, as well as composition and structure of additional samples, we were able to predict their mechanical behavior in unconfined compression, without mechanical testing (R=0.98). Our results suggest that specific information on tissue composition and structure might enable assessment of cartilage mechanics without mechanical testing.     

A nonlinear biphasic viscohyperelastic model for articular cartilage

Experiments on articular cartilage have shown nonlinear stress-strain curves under finite deformations as well as intrinsic viscous effects of the solid phase. The aim of this study was to propose a nonlinear biphasic viscohyperelastic model that combines the intrinsic viscous effects of the proteoglycan matrix with a nonlinear hyperelastic constitutive equation. The proposed equation satisfies objectivity and reduces for uniaxial loading to a solid type viscous model in which the actions of the springs are represented by the hyperelastic function proposed by Holmes and Mow [1990. J. Biomechanics 23, 1145-1156.]. Results of the model, that were efficiently implemented in an updated Lagrangian algorithm, were compared with experimental infinitesimal data reported by DiSilverstro and Suh [2001. J. Biomechanics 34, 519-525.] and showed acceptable fitting for the axial force (R(2)=0.991) and lateral displacement (R(2)=0.914) curves in unconfined compression as well as a good fitting of the axial indentation force curve (R(2)=0.982). In addition, the model showed an excellent fitting of finite-deformation confined compression stress relaxation data reported by Ateshian et al. [1997. J. Biomechanics 30, 1157-1164.] and Huang et al. [2005. J. Biomechanics 38, 799-809.] (R(2)=0.993 and R(2)=0.995, respectively). The constitutive equation may be used to represent the mechanical behavior of the proteoglycan matrix in a fiber reinforced model of articular cartilage.     

The role of viscoelasticity of collagen fibers in articular cartilage: theory and numerical formulation

The relative importance of fluid-dependent and fluid-independent transient mechanical behavior in articular cartilage was examined for tensile and unconfined compression testing using a fibril reinforced model. The collagen matrix of articular cartilage was modeled as viscoelastic using a quasi-linear viscoelastic formulation with strain-dependent elastic modulus, while the proteoglycan matrix was considered as linearly elastic. The collagen viscoelastic properties were obtained by fitting experimental data from a tensile test. These properties were used to investigate unconfined compression testing, and the sensitivity of the properties was also explored. It was predicted that the stress relaxation observed in tensile tests was not caused by fluid pressurization at the macroscopic level. A multi-step tensile stress relaxation test could be approximated using a hereditary integral in which the elastic fibrillar modulus was taken to be a linear function of the fibrillar strain. Applying the same formulation to the radial fibers in unconfined compression, stress relaxation could not be simulated if fluid pressurization were absent. Collagen viscoelasticity was found to slightly weaken fluid pressurization in unconfined compression, and this effect was relatively more significant at moderate strain rates. Therefore, collagen viscoelasticity appears to play an import role in articular cartilage in tensile testing, while fluid pressurization dominates the transient mechanical behavior in compression. Collagen viscoelasticity plays a minor role in the mechanical response of cartilage in unconfined compression if significant fluid flow is present.     

Optical determination of anisotropic material properties of bovine articular cartilage in compression

The precise nature of the material symmetry of articular cartilage in compression remains to be elucidated. The primary objective of this study was to determine the equilibrium compressive Young's moduli and Poisson's ratios of bovine cartilage along multiple directions (parallel and perpendicular to the split line direction, and normal to the articular surface) by loading small cubic specimens (0.9 x 0.9 x 0.8 mm, n =15) in unconfined compression, with the expectation that the material symmetry of cartilage could be determined more accurately with the help of a more complete set of material properties. The second objective was to investigate how the tension-compression nonlinearity of cartilage might alter the interpretation of material symmetry. Optimized digital image correlation was used to accurately determine the resultant strain fields within the specimens under loading. Experimental results demonstrated that neither the Young's moduli nor the Poisson's ratios exhibit the same values when measured along the three loading directions. The main findings of this study are that the framework of linear orthotropic elasticity (as well as higher symmetries of linear elasticity) is not suitable to describe the equilibrium response of articular cartilage nor characterize its material symmetry; a framework which accounts for the distinctly different responses of cartilage in tension and compression is more suitable for describing the equilibrium response of cartilage; within this framework, cartilage exhibits no lower than orthotropic symmetry.     

Ultrasonic measurement of depth-dependent transient behaviors of articular cartilage under compression

We previously reported an ultrasound method for measuring the depth-dependent equilibrium mechanical properties of articular cartilage using quasi-static compression. The objective of this paper was to introduce our recent development for nondestructively measuring the transient depth-dependent strains of full-thickness articular cartilage specimens prepared from bovine patellae. A 50 MHz focused ultrasound transducer was used to collect ultrasound echoes from articular cartilage specimens (n=8) and sponge phantoms with open pores (n=10) during tests of compression and subsequent stress-relaxation. The transient displacements of the tissues at different depths along the compression direction were calculated from the ultrasound echoes using a cross-correlation tracking technique. An LVDT sensor and a load cell were used to measure the overall deformation of the tissue and the applied force, respectively. Results showed that the tissues inside the cartilage layer continued to move during the stress-relaxation phase after the compression was completed. In the equilibrium state, the displacements of the cartilage tissues at the depths of 1/4, 1/2, and 3/4 of the full-thickness reduced by 51%+/-22%, 54%+/-17%, and 50+/-17%, respectively, in comparison with its peak value. However, no similar phenomenon was observed in the sponge phantoms. Our preliminary results demonstrated that this ultrasound method may provide a potential tool for the nondestructive measurement of the transient depth-dependent processes involved in biological and bioengineered soft tissues as well as soft biomaterials under dynamic loading.     

The apparent viscoelastic behavior of articular cartilage--the contributions from the intrinsic matrix viscoelasticity and interstitial fluid flows

Articular cartilage was modeled rheologically as a biphasic poroviscoelastic material. A specific integral-type linear viscoelastic model was used to describe the constitutive relation of the collagen-proteoglycan matrix in shear. For bulk deformation, the matrix was assumed either to be linearly elastic, or viscoelastic with an identical reduced relaxation spectrum as in shear. The interstitial fluid was considered to be incompressible and inviscid. The creep and the rate-controlled stress-relaxation experiments on articular cartilage under confined compression were analyzed using this model. Using the material data available in the literature, it was concluded that both the interstitial fluid flow and the intrinsic matrix viscoelasticity contribute significantly to the apparent viscoelastic behavior of this tissue under confined compression.     

Biomechanical properties of knee articular cartilage

Structure and properties of knee articular cartilage are adapted to stresses exposed on it during physiological activities. In this study, we describe site- and depth-dependence of the biomechanical properties of bovine knee articular cartilage. We also investigate the effects of tissue structure and composition on the biomechanical parameters as well as characterize experimentally and numerically the compression-tension nonlinearity of the cartilage matrix. In vitro mechano-optical measurements of articular cartilage in unconfined compression geometry are conducted to obtain material parameters, such as thickness, Young's and aggregate modulus or Poisson's ratio of the tissue. The experimental results revealed significant site- and depth-dependent variations in recorded parameters. After enzymatic modification of matrix collagen or proteoglycans our results show that collagen primarily controls the dynamic tissue response while proteoglycans affect more the static properties. Experimental measurements in compression and tension suggest a nonlinear compression-tension behavior of articular cartilage in the direction perpendicular to articular surface. Fibril reinforced poroelastic finite element model was used to capture the experimentally found compression-tension nonlinearity of articular cartilage.     

A comparison of cartilage stress-relaxation models in unconfined compression: QLV and stretched exponential in combination with fluid flow

Cartilage exhibits nonlinear viscoelastic behaviour. Various models have been proposed to explain cartilage stress relaxation, but it is unclear whether explicit modelling of fluid flow in unconfined compression is needed. This study compared Fung's quasi-linear viscoelastic (QLV) model with a stretched-exponential model of cartilage stress relaxation and examined each of these models both alone and in combination with a fluid-flow model in unconfined compression. Cartilage explants were harvested from bovine calf patellofemoral joints and equilibrated in tissue culture for 5 days before stress-relaxation testing in unconfined compression at 5% nominal strain. The stretched exponential models fit as well as the QLV models. Furthermore, the average stretched exponential relaxation time determined by this model lies within the range of experimentally measured relaxation times for extracted proteoglycan aggregates, consistent with the hypothesis that the stretched exponential model represents polymeric mechanisms of cartilage viscoelasticity.