Investigation into the biphasic properties of a hydrogel for use in a cushion form replacement joint.

A hydrogel with potential applications in the role of a cushion form replacement joint bearing surface material has been investigated. The material properties are required for further development and design studies and have not previously been quantified. Creep indentation experiments were therefore performed on samples of the hydrogel. The biphasic model developed by Mow and co-workers (Mak et al., 1987; Mow et al., 1989a) was used to curve-fit the experimental data to theoretical solutions in order to extract the three intrinsic biphasic material properties of the hydrogel (aggregate modulus, HA, Poisson's ratio, Vs, and permeability, k). Ranges of material properties were determined: aggregate modulus was calculated to be between 18.4 and 27.5 MPa, Poisson's ratio 0.0-0.307, and permeability 0.012-7.27 x 10(-17) m4/Ns. The hydrogel thus had a higher aggregate modulus than values published for natural normal articular cartilage, the Poisson's ratios were similar to articular cartilage, and finally the hydrogel was found to be less permeable than articular cartilage. The determination of these values will facilitate further numerical analysis of the stress distribution in a cushion form replacement joint.     

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.     

Congruency effects on load bearing in diarthrodial joints

Modelling load bearing in diarthrodial joints is challenging, due to the complexity of the materials, the boundary and interface conditions and the geometry. The articulating surfaces are covered with cartilage layers that are filled with a fluid that plays a major role in load bearing [Mow, V.C., Holmes, M.H., Lai, W.M. (1984) "Survey article: fluid transport and mechanical properties of articular cartilage: a review", Journal of Biomechanics 17(5), 377-394]. Researchers have tended to approximate joint geometry using axisymmetry [Donzelli, P.S., Spilker, R.L., Ateshian, G.A., Mow, V.C. (1999) "Contact analysis of biphasic transversely isotropic cartilage layers and correlations with tissue failure", Journal of Biomechanics 32, 1037-1047], often with a rounded upper articulating surface, creating a form of Hertz problem [Donzelli, P.S., Spilker, R.L., Ateshian, G.A., Mow, V.C. (1999) "Contact analysis of biphasic transversely isotropic cartilage layers and correlations with tissue failure", Journal of Biomechanics 32, 1037-1047]. However, diarthrodial joints (shoulder, hip and knee) are equipped with peripheral structures (glenoid labrum, acetabular labrum and meniscus, respectively) that tend to deepen the joint contact and thus cause initial contact to be established at the periphery of the joint rather than "centrally". The surface geometries are purposefully incongruent, and the incongruency has a significant effect on the stresses, pressures and pressure gradients inside the tissue. The models show the importance of the peripheral structures and the incongruency from a load-bearing perspective. Joint shapes must provide a compromise between demands for load-bearing, lubrication and the supply of nutrients to the chondrocytes of the cartilage and cells of the peripheral structures. Retention and repair of the functionality of these peripheral structures should be a prime consideration in any surgical treatment of an injured joint.     

Measuring fixed charge density of goat articular cartilage using indentation methods and biochemical analysis

An important indicator of osteoarthritis (OA) progression is the loss of proteoglycan (PG) aggregates from the cartilage tissue. Using the indentation creep test, two analytical methods, as previously developed by Lu et al. [Lu, X. L., Miller, C., Chen, F. H., Guo, X. E., Mow, V. C., 2007. The generalized triphasic correspondence principle for simultaneous determination of the mechanical properties and proteoglycan content of articular cartilage by indentation. Journal of Biomechanics 40, 2434-2441 (EPub).], for predicting the fixed charge density (FCD) of goat knee articular cartilage in the normal (control) and degenerated states were compared: (1) a "dual-stage" method to calculate FCD from the mechanical properties of the tissue when tested in isotonic and hypertonic solutions; and (2) a "single-stage" method to predict FCD (as in (1)) assuming an intrinsic Poisson's ratio of 0.05 in the hypertonic state. A biochemical analysis using 1,9-dimethylmethylene blue (DMMB) assay was conducted to directly measure PG content, and hence FCD. The association between the FCD and the aggregate modulus of the tissue was also explored. The mean (+/-S.D.) FCD values measured using the dual-stage method were the closest (control: 0.129+/-0.039, degenerated: 0.046+/-029) to the DMMB results (control: 0.125+/-0.034, degenerated: 0.057+/-0.024) as compared to those of the single-stage method (control: 0.147+/-0.035, degenerated: 0.063+/-0.026). The single-stage method was more reliable (r(2)=0.81) when compared to the dual-stage method (r(2)=0.79). A prediction of FCD from the aggregate modulus generated the least reliable FCD prediction (r(2)=0.68). Because both the dual- and single-stage methods provided reliable FCD estimates for normal and degenerated tissue, the less time-consuming single-stage method was concluded to be the ideal technique for predicting FCD and hence PG content of the tissue.     

Congruency Effects on Load Bearing in Diarthrodial Joints

Modelling load bearing in diarthrodial joints is challenging, due to the complexity of the materials, the boundary and interface conditions and the geometry. The articulating surfaces are covered with cartilage layers that are filled with a fluid that plays a major role in load bearing [Mow, V.C., Holmes, M.H., Lai, W.M. (1984) “Survey article: fluid transport and mechanical properties of articular cartilage: a review”, Journal of Biomechanics 17(5), 377–394]. Researchers have tended to approximate joint geometry using axisymmetry [Donzelli, P.S., Spilker, R.L., Ateshian, G.A., Mow, V.C. (1999) “Contact analysis of biphasic transversely isotropic cartilage layers and correlations with tissue failure”, Journal of Biomechanics 32, 1037–1047], often with a rounded upper articulating surface, creating a form of Hertz problem [Donzelli, P.S., Spilker, R.L., Ateshian, G.A., Mow, V.C. (1999) “Contact analysis of biphasic transversely isotropic cartilage layers and correlations with tissue failure”, Journal of Biomechanics 32, 1037–1047]. However, diarthrodial joints (shoulder, hip and knee) are equipped with peripheral structures (glenoid labrum, acetabular labrum and meniscus, respectively) that tend to deepen the joint contact and thus cause initial contact to be established at the periphery of the joint rather than “centrally”. The surface geometries are purposefully incongruent, and the incongruency has a significant effect on the stresses, pressures and pressure gradients inside the tissue. The models show the importance of the peripheral structures and the incongruency from a load-bearing perspective. Joint shapes must provide a compromise between demands for load-bearing, lubrication and the supply of nutrients to the chondrocytes of the cartilage and cells of the peripheral structures. Retention and repair of the functionality of these peripheral structures should be a prime consideration in any surgical treatment of an injured joint.     

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.     

Time and depth dependent poisson’s ratio of cartilage explained by an inhomogeneous orthotropic fiber embedded biphasic model

A time- and depth-dependent Poisson’s ratio has been observed during unconfined compression experiments on articular cartilage, but existing cartilage models have not fully addressed these phenomena. The goal of this study was to develop a model which is able to predict and explain these phenomena, while also being able to fit other experimental scenarios on full depth cartilage specimens such as confined and unconfined compressions. A biphasic (poroelastic), fiber-embedded cartilage model was developed. The heterogeneous material properties of the cartilage (aggregate modulus, void ratio tensile modulus) were extracted from reported experiments on individual layers of bovine articular cartilage. The nonlinear permeability material constants were found by fitting the overall response to published experimental data from confined compression. The matrix of the cartilage was modelled as an inhomogeneous isotropic biphasic material with nonlinear strain dependent permeability. Orthotropic layers were added as embedded elements to represent collagen fibers. Material parameters for these layers were derived from tensile tests of different layers of cartilage. With these predefined tensile parameters, the model showed a good fit with multi-step confined and unconfined compression experiments (R²=0.984 and 0.977, respectively) and could also predict the depth-dependent Poisson’s ratio (R²=0.981). The highlight of the model is the ability to explain the time-depth dependent Poisson's ratio and, by association, the strong effect of material inhomogeneity on local stress and strain patterns within the cartilage layer. This material model’s response may provide valuable new insight into potential initiation of cartilage fibrillation or delamination in whole-joint simulations.     

A biphasic multiscale study of the mechanical microenvironment of chondrocytes within articular cartilage under unconfined compression

Computational analyses have been used to study the biomechanical microenvironment of the chondrocyte that cannot be assessed by in vitro experimental studies; yet all computational studies thus far have focused on the effect of zonal location (superficial, middle, and deep) on the mechanical microenvironment of chondrocytes. The aim of this paper was to study the effect of both zonal and radial locations on the biomechanical microenvironment of chondrocytes in inhomogeneous cartilage under unconfined stress relaxation. A biphasic multiscale approach was employed and nine chondrocytes in different locations were studied. Hyperelastic biphasic theory and depth-dependent aggregate modulus and permeability of articular cartilage were included in the models. It was found that both zonal and radial locations affected the biomechanical stresses and strains of the chondrocytes. Chondrocytes in the mid-radial location had increased volume during the early stage of the loading process. Maximum principal shear stress at the interface between the chondrocyte and the extracellular matrix (ECM) increased with depth, yet that at the ECM–pericellular matrix (PCM) interface had an inverse trend. Fluid pressure decreased with depth, while the fluid pressure difference between the top and bottom boundaries of the microscale model increased with depth. Regardless of location, fluid was exchanged between the chondrocyte, PCM, and ECM. These findings suggested that even under simple compressive loading conditions, the biomechanical microenvironment of the chondrocytes, PCM and ECM was spatially dependent. The current study provides new insight on chondrocyte biomechanics.     

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.     

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 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.     

High-resolution spatial mapping of shear properties in cartilage

Structural properties of articular cartilage such as proteoglycan content, collagen content and collagen alignment are known to vary over length scales as small as a few microns (Bullough and Goodfellow, 1968; Bi et al., 2006). Characterizing the resulting variation in mechanical properties is critical for understanding how the inhomogeneous architecture of this tissue gives rise to its function. Previous studies have measured the depth-dependent shear modulus of articular cartilage using methods such as particle image velocimetry (PIV) that rely on cells and cell nuclei as fiducial markers to track tissue deformation (Buckley et al., 2008; Wong et al., 2008a). However, such techniques are limited by the density of trackable markers, which may be too low to take full advantage of optical microscopy. This limitation leads to noise in the acquired data, which is often exacerbated when the data is manipulated. In this study, we report on two techniques for increasing the accuracy of tissue deformation measurements. In the first technique, deformations were tracked in a grid that was photobleached on each tissue sample (Bruehlmann et al., 2004). In the second, a numerical technique was implemented that allowed for accurate differentiation of optical displacement measurements by minimizing the propagated experimental error while ensuring that truncation error associated with local averaging of the data remained small. To test their efficacy, we employed these techniques to compare the depth-dependent shear moduli of neonatal bovine and adult human articular cartilage. Using a photobleached grid and numerical optimization to gather and analyze data led to results consistent with those reported previously (Buckley et al., 2008; Wong et al., 2008a), but with increased spatial resolution and characteristic coefficients of variation that were reduced up to a factor of 3. This increased resolution allowed us to determine that the shear modulus of neonatal bovine and adult human tissue both exhibit a global minimum at a depth z of around 100 μm and plateau at large depths. The consistency of the depth dependence of |G*|(Z) for adult human and neonatal bovine tissue suggests a functional advantage resulting from this behavior.     

Tensile and compressive properties of healthy and osteoarthritic human articular cartilage

Osteoarthritis (OA) is a disease affecting articular cartilage and the underlying bone, resulting from many biological and mechanical interacting factors which change the extracellular matrix (ECM) and cells and lead to increasing levels of cartilage degeneration, like softening, fibrillation, ulceration and cartilage loss. The early diagnosis of the disease is fundamental to prevent pain, further tissue degeneration and reduce hospital costs. Although morphological modifications can be detected by modern non-invasive diagnostic techniques, they may not be evident in the early stages of OA. The mechanical properties of articular cartilage are related to its composition and structure and are sensitive to even small changes in the ECM that could occur in early OA. The aim of the present study was to compare the mechanical properties of healthy and OA cartilage using a combined experimental-numerical approach. Experimental assessments consisted of step wise confined and unconfined compression and tension stress relaxation tests on disks (for compression) or strips (for tension) of cartilage obtained from human femoral heads discarded from the operating room after total hip replacement. The numerical model was based on the biphasic theory and included the tension-compression non-linearity. Considering OA samples vs normal samples, the static compressive modulus was 55-68% lower, the permeability was 60-80% higher, the dynamic compressive modulus was 59-64% lower, the static tension modulus was 72-83% lower. The model successfully simulated the experimental tests performed on healthy and OA cartilage and was used in combination with the experimental tests to evaluate the role of different ECM components in the mechanical response of normal and OA cartilage.     

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.     

Topographical variation of the elastic properties of articular cartilage in the canine knee

Equilibrium response of articular cartilage to indentation loading is controlled by the thickness (h) and elastic properties (shear modulus, mu, and Poisson's ratio, nu) of the tissue. In this study, we characterized topographical variation of Poisson's ratio of the articular cartilage in the canine knee joint (N=6). Poisson's ratio was measured using a microscopic technique. In this technique, the shape change of the cartilage disk was visualized while the cartilage was immersed in physiological solution and compressed in unconfined geometry. After a constant 5% axial strain, the lateral strain was measured during stress relaxation. At equilibrium, the lateral-to-axial strain ratio indicates the Poisson's ratio of the tissue. Indentation (equilibrium) data from our prior study (Arokoski et al., 1994. International Journal of Sports Medicine 15, 254-260) was re-analyzed using the Poisson's ratio results at the test site to derive values for shear and aggregate moduli. The lowest Poisson's ratio (0.070+/-0.016) located at the patellar surface of femur (FPI) and the highest (0.236+/-0.026) at the medial tibial plateau (TMI). The stiffest cartilage was found at the patellar groove of femur (micro=0.964+/-0.189MPa, H(a)=2.084+/-0. 409MPa) and the softest at the tibial plateaus (micro=0.385+/-0. 062MPa, H(a)=1.113+/-0.141MPa). Comparison of the mechanical results and the biochemical composition of the tissue (Jurvelin et al., 1988. Engineering in Medicine 17, 157-162) at the matched sites of the canine knee joint indicated a negative correlation between the Poisson's ratio and collagen-to-PG content ratio. This is in harmony with our previous findings which suggested that, in unconfined compression, the degree of lateral expansion in different tissue zones is related to collagen-to-PG ratio of the zone.     

Effect of sulfated glycosaminoglycan digestion on the transverse permeability of medial collateral ligament

Dermatan and chondroitin sulfate glycosaminoglycans (GAGs) comprise over 90% of the GAG content in ligament. Studies of their mechanical contribution to soft tissues have reported conflicting results. Measuring the transient compressive response and biphasic material parameters of the tissue may elucidate the contributions of GAGs to the viscoelastic response to deformation. The hypotheses of the current study were that digestion of sulfated GAGs would decrease compressive stress and aggregate modulus while increasing the permeability of porcine medial collateral ligament (MCL). Confined compression stress relaxation experiments were carried out on porcine MCL and tissue treated with chondroitinase ABC (ChABC). Results were fit to a biphasic constitutive model to derive permeability and aggregate modulus. Bovine articular cartilage was used as a benchmark tissue to verify that the apparatus provided reliable results. GAG digestion removed up to 88% of sulfated GAGs from the ligament. Removal of sulfated GAGs increased the permeability of porcine MCL nearly 6-fold versus control tissues. Peak stress decreased significantly. Bovine articular cartilage exhibited the typical reduction of GAG content and resultant decreases in stress and modulus and increases in permeability with ChABC digestion. Given the relatively small amount of GAG in ligament (<1% of tissue dry weight) and the significant change in peak stress and permeability upon removal of GAGs, sulfated GAGs may play a significant role in maintaining the apposition of collagen fibrils in the transverse direction, thus supporting dynamic compressive loads experienced by the ligament during complex joint motion.     

Lubrication of the human ankle joint in walking with the synovial fluid filtrated by the cartilage with the surface zone worn out: steady pure sliding motion.

A mixture model of synovial fluid filtration by cartilage in the human ankle joint during walking is presented for steady sliding motion of the articular surfaces. In the paper the cartilage surface zone is assumed worn out. The same model has been recently applied to the squeeze-film problem for the human hip joint loaded by the body weight during standing (Hlavácek, Journal of Biomechanics 26, 1145-1150, 1151-1160, 1993; Hlavácek and Novák, Journal of Biomechanics 28, 1193-1198, 1199-1205, 1995). The linear biphasic model for cartilage (elastic porous matrix + ideal fluid) due to Prof. V. C. Mow and his co-workers and the biphasic model for synovial fluid (viscous fluid + ideal fluid), as used in the above-mentioned squeeze-film problem, are applied. For the physiologic parameters of the ankle joint during walking, a continuous synovial fluid film about 1 microm thick is maintained under steady entraining motion according to the classical model without the fluid transport across the articular surface. This is not the case in the filtration model with the cartilage surface zones worn out. On the contrary, this filtration model indicates that synovial fluid is intensively filtrated by such cartilage, so that no continuous fluid film is maintained and a synovial gel layer, about 10(-8) m thick, develops over the majority of the contact. Thus, if the cartilage surface zones are worn out, boundary lubrication should prevail in the ankle joint under steady sliding motion for the mean values of loading and the sliding velocity encountered in walking cycle.     

Biphasic indentation of articular cartilage--II. A numerical algorithm and an experimental study

Part I (Mak et al., 1987, J. Biomechanics 20, 703-714) presented the theoretical solutions for the biphasic indentation of articular cartilage under creep and stress-relaxation conditions. In this study, using the creep solution, we developed an efficient numerical algorithm to compute all three material coefficients of cartilage in situ on the joint surface from the indentation creep experiment. With this method we determined the average values of the aggregate modulus. Poisson's ratio and permeability for young bovine femoral condylar cartilage in situ to be HA = 0.90 MPa, vs = 0.39 and k = 0.44 x 10(-15) m4/Ns respectively, and those for patellar groove cartilage to be HA = 0.47 MPa, vs = 0.24, k = 1.42 x 10(-15) m4/Ns. One surprising finding from this study is that the in situ Poisson's ratio of cartilage (0.13-0.45) may be much less than those determined from measurements performed on excised osteochondral plugs (0.40-0.49) reported in the literature. We also found the permeability of patellar groove cartilage to be several times higher than femoral condyle cartilage. These findings may have important implications on understanding the functional behavior of cartilage in situ and on methods used to determine the elastic moduli of cartilage using the indentation experiments.     

Osteoarthritic changes in the biphasic mechanical properties of the chondrocyte pericellular matrix in articular cartilage

The pericellular matrix (PCM) is a narrow region of cartilaginous tissue that surrounds chondrocytes in articular cartilage. Previous modeling studies indicate that the mechanical properties of the PCM relative to those of the extracellular matrix (ECM) can significantly affect the stress-strain, fluid flow, and physicochemical environments of the chondrocyte, suggesting that the PCM plays a biomechanical role in articular cartilage. The goals of this study were to measure the mechanical properties of the PCM using micropipette aspiration coupled with a linear biphasic finite element model, and to determine the alterations in the mechanical properties of the PCM with osteoarthritis (OA). Using a recently developed isolation technique, chondrons (the chondrocyte and its PCM) were mechanically extracted from non-degenerate and osteoarthritic human cartilage. The transient mechanical behavior of the PCM was well-described by a biphasic model, suggesting that the viscoelastic response of the PCM is attributable to flow-dependent effects, similar to that of the ECM. With OA, the mean Young's modulus of the PCM was significantly decreased (38.7+/-16.2 kPa vs. 23.5+/-12.9 kPa, p < 0.001), and the permeability was significantly elevated (4.19+/-3.78 x10(-17) m(4)/Ns vs. 10.2+/-9.38 x 10(-17) m(4)/Ns, p < 0.01). The Poisson's ratio was similar for both non-degenerate and OA PCM (0.044+/-0.063 vs. 0.030+/-0.068, p > 0.6). These findings suggest that the PCM may undergo degenerative processes with OA, similar to those occurring in the ECM. In combination with previous theoretical models of cell-matrix interactions in cartilage, our findings suggest that changes in the properties of the PCM with OA may have an important influence on the biomechanical environment of the chondrocyte.