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.     

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.     

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.     

Material properties of articular cartilage in the rabbit tibial plateau

The material properties of articular cartilage in the rabbit tibial plateau were determined using biphasic indentation creep tests. Cartilage specimens from matched-pair hind limbs of rabbits approximately 4 months of age and greater than 12 months of age were tested on two locations within each compartment using a custom built materials testing apparatus. A three-way ANOVA was used to determine the effect of leg, compartment, and test location on the material properties (aggregate modulus, permeability, and Poisson's ratio) and thickness of the cartilage for each set of specimens. While no differences were observed in cartilage properties between the left and right legs, differences between compartments were found in each set of specimens. For cartilage from the adolescent group, values for aggregate modulus were 40% less in the medial compartment compared to the lateral compartment, while values for permeability and thickness were greater in the medial compartment compared to the lateral compartment (57% and 30%, respectively). Values for Poisson's ratio were 19% less in the medial compartment compared to the lateral compartment. There was also a strong trend for thickness to differ between test locations. Similar findings were observed for cartilage from the mature group with values for permeability and thickness being greater in the medial compartment compared to the lateral compartment (66% and 34%, respectively). Values for Poisson's ratio were 22% less in the medial compartment compared to the lateral compartment.     

Compressive properties of mouse articular cartilage determined in a novel micro-indentation test method and biphasic finite element model

The mechanical properties of articular cartilage serve as important measures of tissue function or degeneration, and are known to change significantly with osteoarthritis. Interest in small animal and mouse models of osteoarthritis has increased as studies reveal the importance of genetic background in determining predisposition to osteoarthritis. While indentation testing provides a method of determining cartilage mechanical properties in situ, it has been of limited value in studying mouse joints due to the relatively small size of the joint and thickness of the cartilage layer. In this study, we developed a micro-indentation testing system to determine the compressive and biphasic mechanical properties of cartilage in the small joints of the mouse. A nonlinear optimization program employing a genetic algorithm for parameter estimation, combined with a biphasic finite element model of the micro-indentation test, was developed to obtain the biphasic, compressive material properties of articular cartilage. The creep response and material properties of lateral tibial plateau cartilage were obtained for wild-type mouse knee joints, by the micro-indentation testing and optimization algorithm. The newly developed genetic algorithm was found to be efficient and accurate when used with the finite element simulations for nonlinear optimization to the experimental creep data. The biphasic mechanical properties of mouse cartilage in compression (average values: Young's modulus, 2.0 MPa; Poisson's ratio, 0.20; and hydraulic permeability, 1.1 x 10(-16) m4/N-s) were found to be of similar orders of magnitude as previous findings for other animal cartilages, including human, bovine, rat, and rabbit and demonstrate the utility of the new test methods. This study provides the first available data for biphasic compressive properties in mouse cartilage and suggests a promising method for detecting altered cartilage mechanics in small animal models of osteoarthritis.     

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.     

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.     

Comparison of the equilibrium response of articular cartilage in unconfined compression,confined compression and indentation

At mechanical equilibrium, articular cartilage is usually characterized as an isotropic elastic material with no interstitial fluid flow. In this study, the equilibrium properties (Young's modulus, aggregate modulus and Poisson's ratio) of bovine humeral, patellar and femoral cartilage specimens (n=26) were investigated using unconfined compression, confined compression, and indentation tests. Optical measurements of the Poisson's ratio of cartilage were also carried out. Mean values of the Young's modulus (assessed from the unconfined compression test) were 0.80+/-0.33, 0.57+/-0.17 and 0.31+/-0.18MPa and of the Poisson's ratio (assessed from the optical test) 0.15+/-0.06, 0.16+/-0.05 and 0.21+/-0.05 for humeral, patellar, and femoral cartilages, respectively. The indentation tests showed 30-79% (p<0.01) higher Young's modulus values than the unconfined compression tests. In indentation, values of the Young's modulus were independent of the indenter diameter only in the humeral cartilage. The mean values of the Poisson's ratio, obtained indirectly using the mathematical relation between the Young's modulus and the aggregate modulus in isotropic material, were 0.16+/-0.06, 0.21+/-0.05, and 0.26+/-0.08 for humeral, patellar, and femoral cartilages, respectively. We conclude that the values of the elastic parameters of the cartilage are dependent on the measurement technique in use. Based on the similar values of Poisson's ratios, as determined directly or indirectly, the equilibrium response of articular cartilage under unconfined and confined compression is satisfactorily described by the isotropic elastic model. However, values of the isotropic Young's modulus obtained from the in situ indentation tests are higher than those obtained from the in vitro unconfined or confined compression tests and may depend on the indenter size in use.     

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.     

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.     

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.     

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.     

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.     

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.     

Unconfined creep compression of chondrocytes

The study of single cell mechanics offers a valuable tool for understanding cellular milieus. Specific knowledge of chondrocyte biomechanics could lead to elucidation of disease etiologies and the biomechanical factors most critical to stimulating regenerative processes in articular cartilage. Recent studies in our laboratory have suggested that it may be acceptable to approximate the shape of a single chondrocyte as a disc. This geometry is easily utilized for generating models of unconfined compression. In this study, three continuum mechanics models of increasing complexity were formulated and used to fit unconfined compression creep data. Creep curves were obtained from middle/deep zone chondrocytes (n = 15) and separately fit using the three continuum models. The linear elastic solid model yielded a Young's modulus of 2.55+/-0.85 kPa. The viscoelastic model (adapted from the Kelvin model) generated an instantaneous modulus of 2.47+/-0.85 kPa, a relaxed modulus of 1.48+/-0.35 kPa, and an apparent viscosity of 1.92+/-1.80 kPa-s. Finally, a linear biphasic model produced an aggregate modulus of 2.58+/-0.87 kPa, a permeability of 2.57 x 10(-12)+/-3.09 m(4)/N-s, and a Poisson's ratio of 0.069+/-0.021. The results of this study demonstrate that similar values for the cell modulus can be obtained from three models of increasing complexity. The elastic model provides an easy method for determining the cell modulus, however, the viscoelastic and biphasic models generate additional material properties that are important for characterizing the transient response of compressed chondrocytes.     

Material and functional properties of articular cartilage and patellofemoral contact mechanics in an experimental model of osteoarthritis

The purposes of this study were to determine the in situ functional and material properties of articular cartilage in an experimental model of joint injury, and to quantify the corresponding in situ joint contact mechanics. Experiments were performed in the anterior cruciate ligament (ACL) transected knee of the cat and the corresponding, intact contralateral knee, 16 weeks following intervention. Cartilage thickness, stiffness, effective Young’s modulus, and permeability were measured and derived from six locations of the knee. The total contact area and peak pressures in the patellofemoral joint were obtained in situ using Fuji Pressensor film, and comparisons between experimental and contralateral joint were made for corresponding loading conditions. Total joint contact area and peak pressure were increased and decreased significantly (=0.01), respectively, in the experimental compared to the contralateral joint. Articular cartilage thickness and stiffness were increased and decreased significantly (=0.01), respectively, in the experimental compared to the contralateral joint in the four femoral and patellar test locations. Articular cartilage material properties (effective Young’s modulus and permeability) were the same in the ACL-transected and intact joints. These results demonstrate for the first time the effect of changes in articular cartilage properties on the load transmission across a joint. They further demonstrate a substantial change in the joint contact mechanics within 16 weeks of ACL transection. The results were corroborated by theoretical analysis of the contact mechanics in the intact and ACL-transected knee using biphasic contact analysis and direct input of cartilage properties and joint surface geometry from the experimental animals. We conclude that the joint contact mechanics in the ACL-transected cat change within 16 weeks of experimental intervention.     

Experimental determination of the linear biphasic constitutive coefficients of human fetal proximal femoral chondroepiphysis

The mechanical properties of the cartilaginous regions of the proximal femoral epiphysis are an important factor in load transmission through the hip joint of young children. Cylindrical test specimens excised from the chondroepiphysis of human stillborn femoral heads were subjected to uniaxial loading in peripherally-unconfined compression, using a ramp/plateau input strain history. The corresponding load vs time curves were analyzed in terms of a recent analytical solution for a linear biphasic material (the well-known KLM model), allowing calculation of that model's three fundamental constitutive coefficients (permeability, equilibrium modulus and solid-phase Poisson ratio) for this material. The numerical algorithm developed to evaluate the biphasic solution yielded very precise replication of previously published KLM parametric plots. When fitted to experimental load histories, however, the model provided only a rather loose approximation of specimen behavior, due apparently to a substantial underestimation of the transient response component associated with interstitial fluid transport. Averaged over the series, the best-fit values for permeability (2.51 X 10(-15) m4 Ns-1) and equilibrium modulus (0.699 MPa) were in the range of values accepted for human adult articular cartilage. A consequence of the coarseness of the analytical curve fits was that a solid-phase Poisson ratio of 0.0 was inferred for all specimens. The permeability vs equilibrium modulus exhibited a nearly linear (r = 0.74) inverse relationship similar to that reported for adult articular cartilage.     

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.     

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.