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.     

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.     

Purification, cloning, and functional expression of phenylalanine aminomutase: the first committed step in Taxol side-chain biosynthesis

The biomechanical functions of articular cartilage are governed largely by the composition and density of its specialized extracellular matrix. Relationships between matrix density and functional indices such as mechanical properties or interstitial solute diffusivities have been previously explored. However, direct correlations between mechanical properties and solute transport parameters have received less attention, despite potential application of this information for cartilage functional assessment both in vivo and in vitro. The objective of this study was therefore to examine relationships among solute diffusivities, mechanical properties, and matrix density of compressed articular cartilage. Matrix density varied due to natural variation among explants and due to applied static compression. Matrix density of statically compressed cartilage explants was characterized by glycoaminoglycan (GAG) weight fraction and fluid volume fraction, while diffusion coefficients of a wide range of solutes were measured to characterize the transport environment. Explant mechanical properties were characterized by a non-linear Young's modulus (axial stress-strain ratio) and a non-linear Poisson's ratio (radial-to-axial strain ratio). Solute diffusivities were consistently correlated with Young's modulus, as well as with explant GAG weight and fluid volume fractions. Therefore, in vitro mechanical tests may provide a means of assessing transport environments in cartilage-like materials, while in vivo measurements of solute transport (for example with magnetic resonance imaging) may be a useful complement in identifying localized differences in matrix density and mechanical properties.     

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.     

Confined compression and torsion experiments on a pHEMA gel in various bath concentrations

The constitutive behaviour of cartilaginous tissue is the result of complex interaction between electrical, chemical and mechanical forces. Electrostatic interactions between fixed charges and mobile ions are usually accounted for by means of Donnan osmotic pressure. Recent experimental data show, however, that the shear modulus of articular cartilage depends on ionic concentration even if the strain is kept constant. Poisson–Boltzmann simulations suggest that this dependence is intrinsic to the double-layer around the proteoglycan chains. In order to verify this premise, this study measures whether—at a given strain—this ionic concentration–dependent shear modulus is present in a polymerized hydroxy-ethyl-methacrylate gel or not. A combined 1D confined compression and torque experiment is performed on a thin cylindrical hydrogel sample, which is brought in equilibrium with, respectively, 1, 0.1 and 0.03 M NaCl. The sample was placed in a chamber that consists of a stainless steel ring placed on a sintered glass filter, and on top a sintered glass piston. Stepwise ionic loading was cascaded by stepwise 1D compression, measuring the total stress after equilibration of the sample. In addition, a torque experiment was interweaved by applying a harmonic angular displacement and measuring the torque, revealing the relation between aggregate shear modulus and salt concentration at a given strain.     

Direct measurement of the Poisson's ratio of human patella cartilage in tension

Articular cartilage has been shown to exhibit large transverse contractions when loaded in tension, suggesting the existence of large values for the Poisson's ratio. Previous studies have suggested that this effect is dependent on amplitude of applied strain, so that a single Poisson's ratio may not be sufficient to describe cartilage behavior. In this study, the Poisson's ratio (v), toe region modulus (Eo), and linear region modulus (E) of human patellar articular cartilage were calculated in simple tension tests from optical analysis of the two-dimensional strain fields at equilibrium. The Poisson's ratio was found to be independent of strain due to the absence of viscoelastic effects during testing. The Poisson's ratio was found to be significantly higher in the surface zone (1.87 +/- 1.11, p<0.01) than in the middle zone (0.62 +/- 0.23), with no significant correlation of v with age of the cartilage. In general, values for Poisson's ratio were greater than 0.5, suggesting cartilage behavior in tension deviates from isotropy. Reported values for the Poisson's ratio of cartilage in compression have been much lower than values measured here in tension, reflecting a mechanical contribution of the collagen fibers to anisotropy in tension but not compression. The toe-region modulus (Eo) was significantly higher in the surface zone (4.51 +/- 2.78 MPa, n=8) compared to the middle zone (2.51 +/- 1.93 MPa, n=10). In addition, the linear-region modulus (E) in the surface zone, but not middle zone (3.42 +/- 2.17 MPa, n=10), was found to correlate with age (R=0.97, p<0.02) with values of surface zone E equal to 23.92 +/- 12.29 MPa (n=5) for subjects under 70 yr of age, and 4.27 +/- 2.89 MPa (n=3) for subjects over 70 yr. Moduli values and trends with depth were consistent with previous studies of human and animal cartilage. From direct measures of two independent material properties, v and E, we calculated a shear modulus, G, which had not been previously reported for cartilage from tensile testing. Calculated values for surface zone G were 3.64 +/- 1.80 MPa for subjects under 70 yr old and 0.96 +/- 0.69 MPa for subjects over 70 yr old, and were significantly higher in the surface zone than in the middle zone (1.10 +/- 0.78 MPa). This study provides an intrinsic measure for the Poisson's ratio of articular cartilage and its dependence on depth which will be important in understanding the nonlinear tension-compression and anisotropic behaviors of articular cartilage.     

Mechano-acoustic determination of Young's modulus of articular cartilage

The compressive stiffness of an elastic material is traditionally characterized by its Young's modulus. Young's modulus of articular cartilage can be directly measured using unconfined compression geometry by assuming the cartilage to be homogeneous and isotropic. In isotropic materials, Young's modulus can also be determined acoustically by the measurement of sound speed and density of the material. In the present study, acoustic and mechanical techniques, feasible for in vivo measurements, were investigated to quantify the static and dynamic compressive stiffness of bovine articular cartilage in situ. Ultrasound reflection from the cartilage surface, as well as the dynamic modulus were determined with the recently developed ultrasound indentation instrument and compared with the reference mechanical and ultrasound speed measurements in unconfined compression (n=72). In addition, the applicability of manual creep measurements with the ultrasound indentation instrument was evaluated both experimentally and numerically. Our experimental results indicated that the sound speed could predict 47% and 53% of the variation in the Young's modulus and dynamic modulus of cartilage, respectively. The dynamic modulus, as determined manually with the ultrasound indentation instrument, showed significant linear correlations with the reference Young's modulus (r(2)=0.445, p<0.01, n=70) and dynamic modulus (r(2)=0.779, p<0.01, n=70) of the cartilage. Numerical analyses indicated that the creep measurements, conducted manually with the ultrasound indentation instrument, were sensitive to changes in Young's modulus and permeability of the tissue, and were significantly influenced by the tissue thickness. We conclude that acoustic parameters, i.e. ultrasound speed and reflection, are indicative to the intrinsic mechanical properties of the articular cartilage. Ultrasound indentation instrument, when further developed, provides an applicable tool for the in vivo detection of cartilage mechano-acoustic properties. These techniques could promote the diagnostics of osteoarthrosis.     

An analysis of the unconfined compression of articular cartilage

Analytical solutions have been obtained for the internal deformation and fluid-flow fields and the externally observable creep, stress relaxation, and constant strain-rate behaviors which occur during the unconfined compression of a cylindrical specimen of a fluid-filled, porous, elastic solid, such as articular cartilage, between smooth, impermeable plates. Instantaneously, the "biphasic" continuum deforms without change in volume and behaves like an incompressible elastic solid of the same shear modulus. Radial fluid flow then allows the internal fluid pressure to equilibrate with the external environment. The equilibrium response is controlled by the Young's modulus and Poisson's ratio of the solid matrix.     

Prediction of biomechanical properties of articular cartilage with quantitative magnetic resonance imaging

Quantitative magnetic resonance imaging (MRI) is the most potential non-invasive means for revealing the structure, composition and pathology of articular cartilage. Here we hypothesize that cartilage mechanical properties as determined by the macromolecular framework and their interactions can be accessed by quantitative MRI. To test this, adjacent cartilage disk pairs (n=32) were prepared from bovine proximal humerus and patellofemoral surfaces. For one sample, the tissue Young's modulus, aggregate modulus, dynamic modulus and Poisson's ratio were determined in unconfined compression. The adjacent disk was studied at 9.4T to determine the tissue T(2) relaxation time, sensitive to the integrity of the collagen network, and T(1) relaxation time in the presence of Gd-DTPA, a technique developed for the estimation of cartilage proteoglycan (PG) content. Quantitative MRI parameters were able to explain up to 87% of the variations in certain biomechanical parameters. Correlations were further improved when data from the proximal humerus was assessed separately. MRI parameters revealed a topographical variation similar to that of mechanical parameters. Linear regression analysis revealed that Young's modulus of cartilage may be characterized more completely by combining both collagen- and PG-sensitive MRI parameters. The present results suggest that quantitative MRI can provide important information on the mechanical properties of articular cartilage. The results are encouraging with respect to functional imaging of cartilage, although in vivo applicability may be limited by the inferior resolution of clinical MRI instruments.     

Determination of Poisson's ratio of articular cartilage by indentation using different-sized indenters

Articular cartilage is often characterized as an isotropic elastic material with no interstitial fluid flow during instantaneous and equilibrium conditions, and indentation testing commonly used to deduce material properties of Young's modulus and Poisson's ratio. Since only one elastic parameter can be deduced from a single indentation test, some other test method is often used to allow separate measurement of both parameters. In this study, a new method is introduced by which the two material parameters can be obtained using indentation tests alone, without requiring a secondary different type of test. This feature makes the method more suitable for testing small samples in situ. The method takes advantages of the finite layer effect. By indenting the sample twice with different-sized indenters, a nonlinear equation with the Poisson's ratio as the only unknown can be formed and Poisson's ratio obtained by solving the nonlinear equation. The method was validated by comparing the predicted Poisson's ratio for urethane rubber with the manufacturer's supplied value, and comparing the predicted Young's modulus for urethane rubber and an elastic foam material with modulii measured by unconfined compression. Anisotropic and nonhomogeneous finite-element (FE) models of the indentation were developed to aid in data interpretation. Applying the method to bovine patellar cartilage, the tissue Young's modulus was found to be 1.79 +/- 0.59 MPa in instantaneous response and 0.45 +/- 0.26 MPa in equilibrium, and the Poisson's ratio 0.503 +/- 0.028 and 0.463 +/- 0.073 in instantaneous and equilibrium, respectively. The equilibrium Poisson's ratio obtained in our work was substantially higher than those derived from biphasic indentation theory and those optically measured in an unconfined compression test. The finite element model results and examination of viscoelastic-biphasic models suggest this could be due to viscoelastic, inhomogeneity, and anisotropy effects.     

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.     

Direct measurement of osmotic pressure of glycosaminoglycan solutions by membrane osmometry at room temperature

Articular cartilage is a hydrated soft tissue composed of negatively charged proteoglycans fixed within a collagen matrix. This charge gradient causes the tissue to imbibe water and swell, creating a net osmotic pressure that enhances the tissue's ability to bear load. In this study we designed and utilized an apparatus for directly measuring the osmotic pressure of chondroitin sulfate, the primary glycosaminoglycan found in articular cartilage, in solution with varying bathing ionic strength (0.015 M, 0.15 M, 0.5 M, 1 M, and 2 M NaCl) at room temperature. The osmotic pressure (pi) was found to increase nonlinearly with increasing chondroitin sulfate concentration and decreasing NaCl ionic bath environment. Above 1 M NaCl, pi changes negligibly with further increases in salt concentration, suggesting that Donnan osmotic pressure is negligible above this threshold, and the resulting pressure is attributed to configurational entropy. Results of the current study were also used to estimate the contribution of osmotic pressure to the stiffness of cartilage based on theoretical and experimental considerations. Our findings indicate that the osmotic pressure resulting from configurational entropy is much smaller in cartilage (based on an earlier study on bovine articular cartilage) than in free solution. The rate of change of osmotic pressure with compressive strain is found to contribute approximately one-third of the compressive modulus (H(A)(eff)) of cartilage (Pi approximately H(A)(eff)/3), with the balance contributed by the intrinsic structural modulus of the solid matrix (i.e., H(A) approximately 2H(A)(eff)/3). A strong dependence of this intrinsic modulus on salt concentration was found; therefore, it appears that proteoglycans contribute structurally to the magnitude of H(A), in a manner independent of osmotic pressure.     

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.     

Volumetric changes of articular cartilage during stress relaxation in unconfined compression

The time-dependent lateral expansion and load relaxation of cartilage cylinders subjected to unconfined compression were simultaneously recorded. These measurements were used to (1) test the assumption of incompressibility for articular cartilage, (2) measure the Poisson's ratio of articular cartilage in compression and (3) investigate the relationship between stress relaxation and volumetric change. Mechanical tests were performed on fetal, calf, and adult humeral head articular cartilage. The instantaneous Poisson's ratio of adult cartilage was 0.49+/-0.08 (mean+S.D.), thus confirming the assumption of incompressibility for this tissue. The instantaneous Poisson's ratio was significantly lower for calf (0. 38+/-0.04) and fetal cartilage (0.36+/-0.04). The equilibrium Poisson's ratio, i.e. true Poisson's ratio of the solid matrix, was significantly higher for the adult tissue (0.26+/-0.11) compared to both the fetal (0.09+/-0.02) and calf (0.11+/-0.03) cartilage. A linear relationship between time-matched load and lateral expansion after the first minute of stress relaxation was observed.     

The functional environment of chondrocytes within cartilage subjected to compressive loading: a theoretical and experimental approach

A non-invasive methodology (based on video microscopy, optimized digital image correlation and thin plate spline smoothing technique) has been developed to determine the intrinsic tissue stiffness (H(a)) and the intrinsic fixed charge density (c(0)(F)) distribution for hydrated soft tissues such as articular cartilage. Using this technique, the depth-dependent inhomogeneous parameters H(a)(z) and c(0)(F)(z) were determined for young bovine cartilage and incorporated into a triphasic mixture model. This model was then used to predict the mechanical and electrochemical events (stress, strain, fluid/osmotic pressure, and electrical potentials) inside the tissue specimen under a confined compression stress relaxation test. The integration of experimental measurements with theoretical analyses can help to understand the unique material behaviors of articular cartilage. Coupled with biological assays of cell-scale biosynthesis, there is also a great potential in the future to study chondrocyte mechanotransduction in situ with a new level of specificity.     

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.     

Relative contribution of articular cartilage’s constitutive components to load support depending on strain rate

Cartilage is considered a biphasic material in which the solid is composed of proteoglycans and collagen. In biphasic tissue, the hydraulic pressure is believed to bear most of the load under higher strain rates and its dissipation due to fluid flow determines creep and relaxation behavior. In equilibrium, hydraulic pressure is zero and load bearing is transferred to the solid matrix. The viscoelasticity of the collagen network also contributes to its time-dependent behavior, and the osmotic pressure to load bearing in equilibrium. The aim of the present study was to determine the relative contributions of hydraulic pressure, viscoelastic collagen stress, solid matrix stiffness and osmotic pressure to load carriage in cartilage under transient and equilibrium conditions. Unconfined compression experiments were simulated using a fibril-reinforced poroviscoelastic model of articular cartilage, including water, fibrillar viscoelastic collagen and non-fibrillar charged glycosaminoglycans. The relative contributions of hydraulic and osmotic pressures and stresses in the fibrillar and non-fibrillar network were evaluated in the superficial, middle and deep zone of cartilage under five different strain rates and after relaxation. Initially upon loading, the hydraulic pressure carried most of the load in all three zones. The osmotic swelling pressure carried most of the equilibrium load. In the surface zone, where the fibers were loaded in tension, the collagen network carried 20 % of the load for all strain rates. The importance of these fibers was illustrated by artificially modifying the fiber architecture, which reduced the overall stiffness of cartilage in all conditions. In conclusion, although hydraulic pressure dominates the transient behavior during cartilage loading, due to its viscoelastic nature the superficial zone collagen fibers carry a substantial part of the load under transient conditions. This becomes increasingly important with higher strain rates. The interesting and striking new insight from this study suggests that under equilibrium conditions, the swelling pressure generated by the combination of proteoglycans and collagen reinforcement accounts cartilage stiffness for more than 90 % of the loads carried by articular cartilage. This finding is different from the common thought that load is transferred from fluid to solid and is carried by the aggregate modulus of the solid. Rather, it is transformed from hydraulic to osmotic swelling pressure. These results show the importance of considering both (viscoelastic) collagen fibers as well as swelling pressure in studies of the (transient) mechanical behavior of cartilage.     

Fixed Electrical Charges and Mobile Ions Affect the Measurable Mechano-Electrochemical Properties of Charged-Hydrated Biological Tissues: The Articular Cartilage Paradigm

The triphasic constitutive law [Lai, Hou and Mow (1991)] has been shown in some special 1D cases to successfully model the deformational and transport behaviors of charged-hydrated, porous-permeable, soft biological tissues, as typified by articular cartilage. Due to nonlinearities and other mathematical complexities of these equations, few problems for the deformation of such materials have ever been solved analytically. Using a perturbation procedure, we have linearized the triphasic equations with respect to a small imposed axial compressive strain, and obtained an equilibrium solution, as well as a short-time boundary layer solution for the mechano- electrochemical (MEC) fields for such a material under a 2D unconfined compression test. The present results show that the key physical parameter determining the deformational behaviors is the ratio of the perturbation of osmotic pressure to elastic stress, which leads to changes of the measurable elastic coefficients. From the short-time boundary layer solution, both the lateral expansion and the applied load are found to decrease with the square root of time. The predicted deformations, flow fields and stresses are consistent with the analysis of the short time and equilibrium biphasic (i.e., the solid matrix has no attached electric charges) [Armstrong, Lai and Mow (1984)]. These results provide a better understanding of the manner in which fixed electric charges and mobile ions within the tissue contribute to the observed material responses.     

An automated approach for direct measurement of two-dimensional strain distributions within articular cartilage under unconfined compression

An automated approachfor measuring in situ two-dimensional strain fields was developed and validated for its application to cartilage mechanics. This approach combines video microscopy, optimized digital image correlation (DIC), thin-plate spline smoothing (TPSS) and generalized cross-validation (GCV) techniques to achieve the desired efficiency and accuracy. Results demonstrate that sub-pixel accuracies can be achieved for measuring tissue displacements with this methodology with a measurement uncertainty ranging from 0.25 to 0.30 pixels. The deformational gradients (from which the strains are determined) can be evaluated directly using the optimized DIC, with a measurement uncertainty of 0.017 to approximately 0.032. In actual measurements of strain in cartilage, TPSS and differentiation can be used to achieve a more accurate measurement of the gradients from the displacement data. Using this automated approach, the two-dimensional strain fields inside immature bovine carpometacarpal joint cartilage specimens under unconfined compression were characterized (n=21). The depth-dependent apparent elastic modulus and Poisson's ratio were also determined and found to be smallest at the articular surface and increasing with depth. The apparent Poisson's ratio is found to decrease with increasing compressive strain, with values as low as 0.01 observed near the articular surface at 25% compression. The variation of the apparent Poisson's ratio with depth is found to be consistent with a theoretical model of cartilage which accounts for the disparity in its tensile and compressive moduli.     

Characterization of articular cartilage by combining microscopic analysis with a fibril-reinforced finite-element model

Load-bearing characteristics of articular cartilage are impaired during tissue degeneration. Quantitative microscopy enables in vitro investigation of cartilage structure but determination of tissue functional properties necessitates experimental mechanical testing. The fibril-reinforced poroviscoelastic (FRPVE) model has been used successfully for estimation of cartilage mechanical properties. The model includes realistic collagen network architecture, as shown by microscopic imaging techniques. The aim of the present study was to investigate the relationships between the cartilage proteoglycan (PG) and collagen content as assessed by quantitative microscopic findings, and model-based mechanical parameters of the tissue. Site-specific variation of the collagen network moduli, PG matrix modulus and permeability was analyzed. Cylindrical cartilage samples (n=22) were harvested from various sites of the bovine knee and shoulder joints. Collagen orientation, as quantitated by polarized light microscopy, was incorporated into the finite-element model. Stepwise stress-relaxation experiments in unconfined compression were conducted for the samples, and sample-specific models were fitted to the experimental data in order to determine values of the model parameters. For comparison, Fourier transform infrared imaging and digital densitometry were used for the determination of collagen and PG content in the same samples, respectively. The initial and strain-dependent fibril network moduli as well as the initial permeability correlated significantly with the tissue collagen content. The equilibrium Young's modulus of the nonfibrillar matrix and the strain dependency of permeability were significantly associated with the tissue PG content. The present study demonstrates that modern quantitative microscopic methods in combination with the FRPVE model are feasible methods to characterize the structure-function relationships of articular cartilage.