Cartilage interstitial fluid load support in unconfined compression following enzymatic digestion

Interstitial fluid pressurization plays an important role in cartilage biomechanics and is believed to be a primary mechanism of load support in synovial joints. The objective of this study was to investigate the effects of enzymatic degradation on the interstitial fluid load support mechanism of articular cartilage in unconfined compression. Thirty-seven immature bovine cartilage plugs were tested in unconfined compression before and after enzymatic digestion. The peak fluid load support decreased significantly (p < 0.0001) from 84 +/- 10% to 53 +/- 19% and from 80 +/- 10% to 46 +/- 21% after 18-hours digestion with 1.0 u/mg-wet-weight and 0.7 u/mg-wet-weight of collagenase, respectively. Treatment with 0.1 u/ml of chondroitinase ABC for 24 hours also significantly reduced the peak fluid load support from 83 +/- 12% to 48 +/- 16% (p < 0.0001). The drop in interstitial fluid load support following enzymatic treatment is believed to result from a decrease in the ratio of tensile to compressive moduli of the solid matrix.     

Inhomogeneous cartilage properties enhance superficial interstitial fluid support and frictional properties, but do not provide a homogeneous state of stress

It has been well established that articular cartilage is compositionally and mechanically inhomogenous through its depth. To what extent this structural inhomogeneity is a prerequisite for appropriate cartilage function and integrity is not well understood. The first hypothesis to be tested in this study was that the depth-dependent inhomogeneity of the cartilage acts to maximize the interstitial fluid load support at the articular surface, to provide efficient frictional and wear properties. The second hypothesis was that the inhomogeneity produces a more homogeneous state of elastic stress in the matrix than would be achieved with uniform properties. We have, for the first time, simultaneously determined depth-dependent tensile and compressive properties of human patellofemoral cartilage from unconfined compression stress relaxation tests. The results show that the tensile modulus increases significantly from 4.1 +/- 1.9 MPa in the deep zone to 8.3 +/- 3.7 MPa at the superficial zone, while the compressive modulus decreases from 0.73 +/- 0.26 MPa to 0.28 +/- 0.16 MPa. The experimental measurements were then implemented with the finite-element method to compute the response of an inhomogeneous and homogeneous cartilage layer to loading. The finite-element models demonstrate that structural inhomogeneity acts to increase the interstitial fluid load support at the articular surface. However, the state of stress, strain, or strain energy density in the solid matrix remained inhomogeneous through the depth of the articular layer, whether or not inhomogeneous material properties were employed. We suggest that increased fluid load support at the articular surface enhances the frictional and wear properties of articular cartilage, but that the tissue is not functionally adapted to produce homogeneous stress, strain, or strain energy density distributions. Interstitial fluid pressurization, but not a homogeneous elastic stress distribution, appears thus to be a prerequisite for the functional and morphological integrity of the 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.     

Effects of enzymatic degradation on the frictional response of articular cartilage in stress relaxation

It was recently shown experimentally that the friction coefficient of articular cartilage correlates with the interstitial fluid pressurization, supporting the hypothesis that interstitial water pressurization plays a fundamental role in the frictional response by supporting most of the load during the early time response. A recent study showed that enzymatic treatment with chondroitinase ABC causes a decrease in the maximum fluid load support of bovine articular cartilage in unconfined compression. The hypothesis of this study is that treatment with chondroitinase ABC will increase the friction coefficient of articular cartilage in stress relaxation. Articular cartilage samples (n = 34) harvested from the femoral condyles of five bovine knee joints (1-3 months old) were tested in unconfined compression with simultaneous continuous sliding (+/-1.5 mm at 1 mm/s) under stress relaxation. Results showed a significantly higher minimum friction coefficient in specimens treated with 0.1 micro/ml of chondroitinase ABC for 24 h (micro(min) = 0.082+/-0.024) compared to control specimens (micro(min) = 0.047+/-0.014). Treated samples also exhibited higher equilibrium friction coefficient (micro(eq) = 0.232+/-0.049) than control samples (micro(eq) = 0.184+/-0.036), which suggest that the frictional response is greatly influenced by the degree of tissue degradation. The fluid load support was predicted from theory, and the maximum value (as a percentage of the total applied load) was lower in treated specimens (77+/-12%) than in control specimens (85+/-6%). Based on earlier findings, the increase in the ratio micro(min)/micro(eq) may be attributed to the decrease in fluid load support.     

Strain-rate dependence of cartilage stiffness in unconfined compression: the role of fibril reinforcement versus tissue volume change in fluid pressurization

The strain and strain-rate-dependent response of articular cartilage in unconfined compression was studied theoretically. The transient stress and stiffness of cartilage were determined for strain rates ranging from zero to infinity. It is shown, for a given compressive strain, that the axial stress initially increases quickly as a function of strain rate, and then increases progressively more slowly towards the stress corresponding to the instantaneous response. The volume change of the tissue does not give its transient stiffness uniquely, because of the strong strain-rate dependence. The variation of tissue stiffness is primarily determined by the transient stiffness of the radial fibrils. Load sharing between the solid matrix and fluid pressurization also depends on the strain rate. At 15% axial compression, the matrix bears more than 80% of the applied load at a strain rate of 0.005%/s, while the fluid pressurization contributes more than 80% of the load at a strain rate of 0.15%/s. These results show the interplay between fibril reinforcement and fluid pressurization in articular cartilage: the fluid drives fibril stiffening which in turn produces high pore pressure at high strain rates.As a secondary objective of the present work, a fibrillar continuum element was formulated to replace the fibrillar spring element used previously in fibril-reinforced modeling, in order to eliminate the deformation incompatibility between the spring system and the nonfibrillar matrix. The results obtained using the two fibrillar elements were compared with the closed-form solutions for the static and instantaneous responses for the case of large deformation. It was found for unconfined compression that using the spring elements did not generally result in greater numerical errors than using the fibrillar continuum elements.     

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.     

Streaming potentials maps are spatially resolved indicators of amplitude, frequency and ionic strength dependant responses of articular cartilage to load.

Streaming potential distributions were measured on the surface of articular cartilage in uniaxial unconfined compression using a linear array of microelectrodes. Potential profiles were obtained for sinusoidal and ramp/stress-relaxation displacements and exhibited dependencies on radial position, sinusoidal amplitude and frequency, time during stress relaxation, and on ionic strength. The measurements agreed with trends predicted by biphasic and related models. In particular, the absolute potential amplitude was maximal at the disk center, as was the predicted fluid pressure and the potential gradient (the electric field) was seen to be maximal at the disk periphery, as was the predicted fluid velocity. We also observed a similarity between non-linear behavior of streaming potential amplitude and load amplitude with respect to sinusoidal displacement amplitude. Taken together, these results support many of the phenomena concerning relative fluid-solid movement and fluid pressurization predicted by biphasic and related models, and they indicate the general utility of spatially resolved measurements of streaming potentials for the investigation of electromechanical phenomena in tissues. For example, these streaming potential maps could be used to non-destructively diagnose cartilage extracellular matrix composition and function, as well as to quantify spatially and temporally varying physical signals in cartilage that can induce cellular and extracellular biological responses to load.     

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.     

Dynamic unconfined compression of articular cartilage under a cyclic compressive load

To study the effect of dynamic mechanical force on cartilage metabolism, many investigators have applied a cyclic compressive load to cartilage disc explants in vitro. The most frequently used in vitro testing protocol has been the cyclic unconfined compression of articular cartilage in a bath of culture medium. Cyclic compression has been achieved by applying either a prescribed cyclic displacement or a prescribed cyclic force on a loading platen placed on the top surface of a cylindrical cartilage disc. It was found that the separation of the loading platen from the tissue surface was likely when a prescribed cyclic displacement was applied at a high frequency.The purpose of the present study was to simulate mathematically the dynamic behavior of a cylindrical cartilage disc subjected to cyclic unconfined compression under a dynamic force boundary condition protocol, and to provide a parametric analysis of mechanical deformations within the extracellular matrix. The frequency-dependent dynamic characteristics of dilatation, hydrostatic pressure and interstitial fluid velocity were analyzed over a wide range of loading frequencies without the separation of the loading platen. The result predicted that a cyclic compressive force created an oscillating positive-negative hydrostatic pressure together with a forced circulation of interstitial fluid within the tissue matrix. It was also found that the load partitioning mechanism between the solid and fluid phases was a function of loading frequency. At a relatively high loading frequency, a localized dynamic zone was developed near the peripheral free surface of the cartilage disc, where a large dynamic pressure gradient exists, causing vigorous interstitial fluid flow.     

The Mechanical Behaviour of Chondrocytes Predicted with a Micro-structural Model of Articular Cartilage

The integrity of articular cartilage depends on the proper functioning and mechanical stimulation of chondrocytes, the cells that synthesize extracellular matrix and maintain tissue health. The biosynthetic activity of chondrocytes is influenced by genetic factors, environmental influences, extracellular matrix composition, and mechanical factors. The mechanical environment of chondrocytes is believed to be an important determinant for joint health, and chondrocyte deformation in response to mechanical loading is speculated to be an important regulator of metabolic activity. In previous studies of chondrocyte deformation, articular cartilage was described as a biphasic material consisting of a homogeneous, isotropic, linearly elastic solid phase, and an inviscid fluid phase. However, articular cartilage is known to be anisotropic and inhomogeneous across its depth. Therefore, isotropic and homogeneous models cannot make appropriate predictions for tissue and cell stresses and strains. Here, we modelled articular cartilage as a transversely isotropic, inhomogeneous (TI) material in which the anisotropy and inhomogeneity arose naturally from the microstructure of the depth-dependent collagen fibril orientation and volumetric fraction, as well as the chondrocyte shape and volumetric fraction. The purpose of this study was to analyse the deformation behaviour of chondrocytes using the TI model of articular cartilage. In order to evaluate our model against experimental results, we simulated indentation and unconfined compression tests for nominal compressions of 15%. Chondrocyte deformations were analysed as a function of location within the tissue. The TI model predicted a non-uniform behaviour across tissue depth: in indentation testing, cell height decreased by 43% in the superficial zone and between 11 and 29% in the deep zone. In unconfined compression testing, cell height decreased by 32% in the superficial zone, 25% in the middle, and 18% in the deep zones. This predicted non-uniformity is in agreement with experimental studies. The novelty of this study is the use of a cartilage material model accounting for the intrinsic inhomogeneity and anisotropy of cartilage caused by its microstructure.     

Experimental verification of the roles of intrinsic matrix viscoelasticity and tension-compression nonlinearity in the biphasic response of cartilage

A biphasic-CLE-QLV model proposed in our recent study [2001, J. Biomech. Eng., 123, pp. 410-417] extended the biphasic theory of Mow et al. [1980, J. Biomech. Eng., 102, pp. 73-84] to include both tension-compression nonlinearity and intrinsic viscoelasticity of the cartilage solid matrix by incorporating it with the conewise linear elasticity (CLE) model [1995, J. Elasticity, 37, pp. 1-38] and the quasi-linear viscoelasticity (QLV) model [Biomechanics: Its foundations and objectives, Prentice Hall, Englewood Cliffs, 1972]. This model demonstrates 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 viscoelastic effects, as well as tension-compression nonlinearity. The objective of this study is to directly test this biphasic-CLE-QLV model against experimental data from unconfined compression stress-relaxation tests at slow and fast strain rates as well as dynamic loading. Twelve full-thickness cartilage cylindrical plugs were harvested from six bovine glenohumeral joints and multiple confined and unconfined compression stress-relaxation tests were performed on each specimen. The material properties of specimens were determined by curve-fitting the experimental results from the confined and unconfined compression stress relaxation tests. The findings of this study demonstrate that the biphasic-CLE-QLV model is able to describe the strain-rate-dependent mechanical behaviors of articular cartilage in unconfined compression as attested by good agreements between experimental and theoretical curvefits (r2 = 0.966 +/- 0.032 for testing at slow strain rate; r2 = 0.998 +/- 0.002 for testing at fast strain rate) and predictions of the dynamic response (r2 = 0.91 +/- 0.06). This experimental study also provides supporting evidence for the hypothesis that both tension-compression nonlinearity and intrinsic viscoelasticity of the solid matrix of cartilage are necessary for modeling the transient and equilibrium responses of this tissue in tension and compression. Furthermore, the biphasic-CLE-QLV model can produce better predictions of the dynamic modulus of cartilage in unconfined dynamic compression than the biphasic-CLE and biphasic poroviscoelastic models, indicating that intrinsic viscoelasticity and tension-compression nonlinearity of articular cartilage may play important roles in the load-support mechanism of cartilage under physiologic loading.     

The influence of the fixed negative charges on mechanical and electrical behaviors of articular cartilage under unconfined compression

Unconfined compression test has been frequently used to study the mechanical behaviors of articular cartilage, both theoretically and experimentally. It has also been used in explant and gel-cell-complex studies in tissue engineering. In biphasic and poroelastic theories, the effect of charges fixed on the proteoglycan macromolecules in articular cartilage is embodied in the apparent compressive Young's modulus and the apparent Poisson's ratio of the tissue, and the fluid pressure is considered to be the portion above the osmotic pressure. In order to understand how proteoglycan fixed charges might affect the mechanical behaviors of articular cartilage, and in order to predict the osmotic pressure and electric fields inside the tissue in this experimental configuration, it is necessary to use a model that explicitly takes into account the charged nature of the tissue and the flow of ions within its porous interstices. In this paper, we used a finite element model based on the triphasic theory to study how fixed charges in the porous-permeable soft tissue can modulate its mechanical and electrochemical responses under a step displacement in unconfined compression. The results from finite element calculations showed that: 1) A charged tissue always supports a larger load than an uncharged tissue of the same intrinsic elastic moduli. 2) The apparent Young's modulus (the ratio of the equilibrium axial stress to the axial strain) is always greater than the intrinsic Young's modulus of an uncharged tissue. 3) The apparent Poisson's ratio (the negative ratio of the lateral strain to the axial strain) is always larger than the intrinsic Poisson's ratio of an uncharged tissue. 4) Load support derives from three sources: intrinsic matrix stiffness, hydraulic pressure and osmotic pressure. Under the unconfined compression, the Donnan osmotic pressure can constitute between 13%-22% of the total load support at equilibrium. 5) During the stress-relaxation process following the initial instant of loading, the diffusion potential (due to the gradient of the fixed charge density and the associated gradient of ion concentrations) and the streaming potential (due to fluid convection) compete against each other. Within the physiological range of material parameters, the polarity of the electric potential depends on both the mechanical properties and the fixed charge density (FCD) of the tissue. For softer tissues, the diffusion effects dominate the electromechanical response, while for stiffer tissues, the streaming potential dominates this response. 6) Fixed charges do not affect the instantaneous strain field relative to the initial equilibrium state. However, there is a sudden increase in the fluid pressure above the initial equilibrium osmotic pressure. These new findings are relevant and necessary for the understanding of cartilage mechanics, cartilage biosynthesis, electromechanical signal transduction by chondrocytes, and tissue engineering.     

Articular cartilage mechanical and biochemical property relations before and after in vitro growth

The aim of this study was to design in vitro growth protocols that can comprehensively quantify articular cartilage structure–function relations via measurement of mechanical and biochemical properties. Newborn bovine patellofemoral groove articular cartilage explants were tested sequentially in confined compression (CC), unconfined compression (UCC), and torsional shear before (D0, i.e. day zero) and after (D14, i.e. day 14) unstimulated in vitro growth. The contents of collagen (COL), collagen-specific pyridinoline (PYR) crosslinks, glycosaminoglycan, and DNA significantly decreased during in vitro growth; consequently, a wide range of biochemical properties existed for investigating structure–function relations when pooling the D0 and D14 groups. All D0 mechanical properties were independent of compression strain while only Poisson's ratios were dependent on direction (i.e. anisotropic). Select D0 and D14 group mechanical properties were correlated with biochemical measures; including (but not limited to) results that CC/UCC moduli and UCC Poisson's ratios were correlated with COL and PYR. COL network weakening during in vitro growth due to reduced COL and PYR was accompanied by reduced CC/UCC moduli and increased UCC Poisson's ratios.     

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.     

Cyclic compression of articular cartilage explants is associated with progressive consolidation and altered expression pattern of extracellular matrix proteins.

In this study, we investigated the biosynthetic response of full thickness, adult bovine articular cartilage explants to 45 h of static and cyclic unconfined compression. The cyclic compression of articular cartilage resulted in a progressive consolidation of the cartilage matrix. The oscillatory loading increased protein synthesis ([35S]methionine incorporation) by as much as 50% above free swelling control values, but had an inhibitory influence on proteoglycan synthesis ([35SO4] incorporation). As expected, static compression was associated with a dose-dependent decrease in biosynthetic activity. ECM oligomeric proteins which were most affected by mechanical loading were fibronectin and cartilage oligomeric matrix protein (COMP). Static compression at all amplitudes caused a significant increase in fibronectin synthesis over free swelling control levels. Cyclic compression of articular cartilage at 0.1 Hz and higher was consistently associated with a dramatic increase in the synthesis of COMP as well as fibronectin. The biosynthetic activity of chondrocytes appears to be sensitive to both the frequency and amplitude of the applied load. The results of this study support the hypothesis that cartilage tissue can remodel its extracellular matrix in response to alterations in functional demand.     

Hip chondrolabral mechanics during activities of daily living: Role of the labrum and interstitial fluid pressurization

Osteoarthritis of the hip can result from mechanical factors, which can be studied using finite element (FE) analysis. FE studies of the hip often assume there is no significant loss of fluid pressurization in the articular cartilage during simulated activities and approximate the material as incompressible and elastic. This study examined the conditions under which interstitial fluid load support remains sustained during physiological motions, as well as the role of the labrum in maintaining fluid load support and the effect of its presence on the solid phase of the surrounding cartilage. We found that dynamic motions of gait and squatting maintained consistent fluid load support between cycles, while static single-leg stance experienced slight fluid depressurization with significant reduction of solid phase stress and strain. Presence of the labrum did not significantly influence fluid load support within the articular cartilage, but prevented deformation at the cartilage edge, leading to lower stress and strain conditions in the cartilage. A morphologically accurate representation of collagen fibril orientation through the thickness of the articular cartilage was not necessary to predict fluid load support. However, comparison with simplified fibril reinforcement underscored the physiological importance. The results of this study demonstrate that an elastic incompressible material approximation is reasonable for modeling a limited number of cyclic motions of gait and squatting without significant loss of accuracy, but is not appropriate for static motions or numerous repeated motions. Additionally, effects seen from removal of the labrum motivate evaluation of labral reattachment strategies in the context of labral repair.     

Biomechanics of single zonal chondrocytes

Articular cartilage has a distinct zonal architecture, and previous work has shown that chondrocytes from different zones exhibit variations in gene expression and biosynthesis. In this study, the material properties of single chondrocytes from the superficial and middle/deep zones of bovine distal metatarsal articular cartilage were determined using unconfined compression and digital videocapture. To determine the viscoelastic properties of zonal chondrocytes, unconfined creep compression experiments were performed and the resulting creep curves of individual cells were fit using a standard linear viscoelastic solid model. In the model, a fixed value of the Poisson's ratio was used, determined optically from direct compression of middle/deep chondrocytes. The two approaches used in this study yielded the following average material properties of single chondrocytes: Poisson's ratio of 0.26+/-0.08, instantaneous modulus of 1.06+/-0.82 kPa, relaxed modulus of 0.78+/-0.58 kPa, and apparent viscosity of 4.08+/-7.20 kPa s. Superficial zone chondrocytes were found to be significantly stiffer than middle/deep zone chondrocytes. Attachment time did not affect the stiffness of the cells. The zonal variation in viscoelastic properties may result from the distinct mechanical environments experienced by the cells in vivo. Identifying intrinsic differences in the biomechanics of superficial and middle/deep zone chondrocytes is an important component in understanding how biomechanics influence articular cartilage health and disease.     

Biomechanical properties of human articular cartilage under compressive loads

The function of articular cartilage is to support and distribute loads and to provide lubrication in the diarthrodial joints. Cartilage function is described by proper mechanical and rheological properties, strain and depth-dependent, which are not completely assessed. Unconfined and confined compression are commonly used to evaluate the Young's modulus (E) and the aggregate modulus (H(A)), respectively. The Poisson's ratio (nu) can be calculated indirectly from the equilibrium compression data, or using the biphasic indentation technique; it has recently been optically evaluated by using video microscopy during unconfined compression. The transient response of articular cartilage during confined compression depends on its permeability k; a constant value of k can be easily identified by a simple analytical model of confined compression tests, whereas more complex models or direct measurements (permeation tests) are needed to study the permeability dependence on deformation. A poroelastic finite element model of articular cartilage was developed for this purpose. The elastic parameters (E,nu) of the model were evaluated performing unconfined compression creep tests on human articular cartilage disks, whereas k was identified from the confined test response. Our combined experimental and computational method can be used to identify the parameters that define the permeability dependence on deformation, as a function of depth from articular surface.     

Anisotropy, inhomogeneity, and tension-compression nonlinearity of human glenohumeral cartilage in finite deformation

The tensile and compressive properties of human glenohumeral cartilage were determined by testing 120 rectangular strips in uniaxial tension and 70 cylindrical plugs in confined compression, obtained from five human glenohumeral joints. Specimens were harvested from five regions across the articular surface of the humeral head and two regions on the glenoid. Tensile strips were obtained along two orientations, parallel and perpendicular to the split-line directions. Two serial slices through the thickness, corresponding to the superficial and middle zones of the cartilage layers, were prepared from each tensile strip and each compressive plug. The equilibrium tensile modulus and compressive aggregate modulus of cartilage were determined from the uniaxial tensile and confined compression tests, respectively. Significant differences in the tensile moduli were found with depth and orientation relative to the local split-line direction. Articular cartilage of the humeral head was significantly stiffer in tension than that of the glenoid. There were significant differences in the aggregate compressive moduli of articular cartilage between superficial and middle zones in the humeral head. Furthermore, tensile and compressive stress-strain responses exhibited nonlinearity under finite strain, while the tensile modulus differed by up to two orders of magnitude from the compressive aggregate modulus at 0% strain, demonstrating a high degree of tension-compression nonlinearity. The complexity of the mechanical properties of human glenohumeral cartilage was exposed in this study, showing anisotropy, inhomogeneity, and tension-compression nonlinearity within the same joint. The observed differences in the tensile properties of human glenohumeral cartilage suggest that the glenoid may be more susceptible to cartilage degeneration than the humeral head.     

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.