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.     

Effect of fiber orientation and strain rate on the nonlinear uniaxial tensile material properties of tendon

Tendons are exposed to complex loading scenarios that can only be quantified by mathematical models, requiring a full knowledge of tendon mechanical properties. This study measured the anisotropic, nonlinear, elastic material properties of tendon. Previous studies have primarily used constant strain-rate tensile tests to determine elastic modulus in the fiber direction. Data for Poisson's ratio aligned with the fiber direction and all material properties transverse to the fiber direction are sparse. Additionally, it is not known whether quasi-static constant strain-rate tests represent equilibrium elastic tissue behavior. Incremental stress-relaxation and constant strain-rate tensile tests were performed on sheep flexor tendon samples aligned with the tendon fiber direction or transverse to the fiber direction to determine the anisotropic properties of toe-region modulus (E0), linear-region modulus (E), and Poisson's ratio (v). Among the modulus values calculated, only fiber-aligned linear-region modulus (E1) was found to be strain-rate dependent. The E1 calculated from the constant strain-rate tests were significantly greater than the value calculated from incremental stress-relaxation testing. Fiber-aligned toe-region modulus (E(1)0 = 10.5 +/- 4.7 MPa) and linear-region modulus (E1 = 34.0 +/- 15.5 MPa) were consistently 2 orders of magnitude greater than transverse moduli (E(2)0 = 0.055 +/- 0.044 MPa, E2 = 0.157 +/- 0.154 MPa). Poisson's ratio values were not found to be rate-dependent in either the fiber-aligned (v12 = 2.98 +/- 2.59, n = 24) or transverse (v21 = 0.488 +/- 0.653, n = 22) directions, and average Poisson's ratio values in the fiber-aligned direction were six times greater than in the transverse direction. The lack of strain-rate dependence of transverse properties demonstrates that slow constant strain-rate tests represent elastic properties in the transverse direction. However, the strain-rate dependence demonstrated by the fiber-aligned linear-region modulus suggests that incremental stress-relaxation tests are necessary to determine the equilibrium elastic properties of tendon, and may be more appropriate for determining the properties to be used in elastic mathematical models.     

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.     

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.     

Anisotropic strain-dependent material properties of bovine articular cartilage in the transitional range from tension to compression

Articular cartilage exhibits complex mechanical properties such as anisotropy, inhomogeneity and tension-compression nonlinearity. This study proposes and demonstrates that the application of compressive loading in the presence of osmotic swelling can be used to acquire a spectrum of incremental cartilage moduli (EYi) and Poisson's ratios (upsilon ij) from tension to compression. Furthermore, the anisotropy of the tissue can be characterized in both tension and compression by conducting these experiments along three mutually perpendicular loading directions: parallel to split-line (1-direction), perpendicular to split-line (2-direction) and along the depth direction (3-direction, perpendicular to articular surface), accounting for tissue inhomogeneity between the surface and deep layers in the latter direction. Tensile moduli were found to be strain-dependent while compressive moduli were nearly constant. The peak tensile (+) Young's moduli in 0.15M NaCl were E+Y1=3.1+/-2.3, E+Y2=1.3+/-0.3, E+Y3(Surface)=0.65+/-0.29 and E+Y3(Deep)=2.1+/-1.2 MPa. The corresponding compressive (-) Young's moduli were E-Y1=0.23+/-0.07, E-Y2=0.22+/-0.07, E-Y3(Surface)=0.18+/-0.07 and E-Y3(Deep)=0.35+/-0.11 MPa. Peak tensile Poisson's ratios were upsilon+12=0.22+/-0.06, upsilon+21=0.13+/-0.07, upsilon+31(Surface)=0.10+/-0.03 and upsilon+31(Deep)=0.20+/-0.05 while compressive Poisson's ratios were upsilon-12=0.027+/-0.012, upsilon-21=0.017+/-0.07, upsilon-31(Surface)=0.034+/-0.009 and upsilon-31(Deep)=0.065+/-0.024. Similar measurements were also performed at 0.015 M and 2 M NaCl, showing strong variations with ionic strength. Results indicate that (a) a smooth transition occurs in the stress-strain and modulus-strain responses between the tensile and compressive regimes, and (b) cartilage exhibits orthotropic symmetry within the framework of tension-compression nonlinearity. The strain-softening behavior of cartilage (the initial decrease in EYi with increasing compressive strain) can be interpreted in the context of osmotic swelling and tension-compression nonlinearity.     

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.     

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.     

Ultrasound indentation of bovine knee articular cartilage in situ

We have earlier developed a handheld ultrasound indentation instrument for the diagnosis of articular cartilage degeneration. In ultrasound indentation, cartilage is compressed with the ultrasound transducer. Tissue thickness and deformation are calculated from the A-mode ultrasound signal and the stress applied is registered with the strain gauges. In this study, the applicability of the ultrasound indentation instrument to quantify site-dependent variation in the mechano-acoustic properties of bovine knee cartilage was investigated. Osteochondral blocks (n=6 per site) were prepared from the femoral medial condyle (FMC), the lateral facet of the patello-femoral groove (LPG) and the medial tibial plateau (MTP). Cartilage stiffness (dynamic modulus, E(dyn)), as obtained with the ultrasound indentation instrument in situ, correlated highly linearly (r=0.913, p<0.01) with the values obtained using the reference material-testing device in vitro. Reproducibility (standardized coefficient of variation) of the ultrasound indentation measurements was 5.2%, 1.7% and 3.1% for E(dyn), ultrasound reflection coefficient of articular surface (R) and thickness, respectively. E(dyn) and R were site dependent (p<0.05, Kruskall-Wallis H test). E(dyn) was significantly higher (p<0.05, Kruskall-Wallis Post Hoc test) in LPG (mean+/-SD: 10.1+/-3.1MPa) than in MTP (2.9+/-1.4MPa). In FMC, E(dyn) was 4.6+/-1.3MPa. R was significantly (p<0.05) lower at MTP (2.0+/-0.7%) than at other sites (FMC: 4.2+/-0.9%; LPG: 4.4+/-0.8%). Cartilage glycosaminoglycan concentration, as quantified with the digital densitometry, correlated positively with E(dyn) (r=0.678, p<0.01) and especially with the equilibrium Young's modulus (reference device, r=0.874, p<0.01) but it was not associated with R (r=0.294, p=0.24). We conclude that manual measurements are reproducible and the instrument may be used for detection of cartilage quality in situ. Especially, combined measurement of thickness, E(dyn) and R provides valuable diagnostic information on cartilage status.     

Indentation stiffness of young canine knee articular cartilage—Influence of strenuous joint loading

The indentation stiffness of knee articular cartilage subjected to strenuous physical training (SPT: treadmill running 20 km day⁻¹ for 15 weeks, n = 6) of young Beagles was tested and compared to that obtained from age-matched (55 weeks, n = 9) controls. The mathematical solution for the shear modulus, as determined from indentation of an elastic layer bonded to a rigid half space, was extended to small Poisson's ratios and applied to the analysis of cartilage response after a step stress (0.39 MPa) application. In these measurements with an impervious, plane-ended indenter, the equilibrium deformation was systematically greater than values predicted from the instant response by the linear biphasic theory. Therefore, the accurate determination of Poisson's ratio from the creep curves was not possible. The mean shear modulus (calculated by using the deformation at 900 s after load application and assuming a constant Poisson's ratio of 0.40 for the matrix) of canine knee articular cartilage was 0.37 MPa. While the cartilage thickness was not affected by SPT, the cartilage of the lateral tibial plateau was stiffer (13.3%, p<0.05) than that in controls. However, in the femoral condyles, the stiffness was at the control level or even below. Our results on cartilage structure and properties suggest that SPT, in contrast to our previous findings with moderate training, does not necessarily improve the biological properties of articular cartilage in young animals.     

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.     

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.     

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.     

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.     

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

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

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.     

Influence of stress magnitude on water loss and chondrocyte viability in impacted articular cartilage

The objective of this study was to assess mechano-biological response of articular cartilage when subjected to a single impact stress. Mature bovine cartilage explants were impacted with peak stresses ranging from 10 to 60 MPa at a stress rate of 350 MPa/s. Water loss, matrix axial deformation, dynamic impact modulus (DIM), and cell viability were measured immediately after impaction. The water loss through the articular surface (AS) was small and ranged from 1% to 6% with increasing peak stress. The corresponding axial strains ranged from 2.5% to 25%, respectively, while the DIM was 455.9 +/- 111.9 MPa. Chondrocyte death started at the articular surface and increased in depth to a maximum of 6% (70 microns) of the cartilage thickness at the highest stress. We found that the volumetric (axial) strain was more than twice the amount of water loss at the highest peak stress. Furthermore, specimens impacted such that the interstitial water was forced through the deep zone (DZ) had less water loss, a higher DIM, and no cell death. These findings appear to be due to matrix compaction in the superficial region causing higher compressive strains to occur at the surface rather than in the deeper zones.     

Anisotropic and inhomogeneous tensile behavior of the human anulus fibrosus: experimental measurement and material model predictions

The anulus fibrosus (AF) of the intervertebral disc exhibits spatial variations in structure and composition that give rise to both anisotropy and inhomogeneity in its material behaviors in tension. In this study, the tensile moduli and Poisson's ratios were measured in samples of human AF along circumferential, axial, and radial directions at inner and outer sites. There was evidence of significant inhomogeneity in the linear-region circumferential tensile modulus (17.4+/-14.3 MPa versus 5.6+/-4.7 MPa, outer versus inner sites) and the Poisson's ratio v21 (0.67+/-0.22 versus 1.6+/-0.7, outer versus inner), but not in the axial modulus (0.8+/-0.9 MPa) or the Poisson's ratios V12 (1.8+/-1.4) or v13 (0.6+/-0.7). These properties were implemented in a linear an isotropic material model of the AF to determine a complete set of model properties and to predict material behaviors for the AF under idealized kinematic states. These predictions demonstrate that interactions between fiber populations in the multilamellae AF significantly contribute to the material behavior, suggesting that a model for th     

Alterations in the mechanical properties of the human chondrocyte pericellular matrix with osteoarthritis

In articular cartilage, chondrocytes are surrounded by a pericellular matrix (PCM), which together with the chondrocyte have been termed the "chondron." While the precise function of the PCM is not know there has been considerable speculation that it plays a role in regulating the biomechanical environment of the chondrocyte. In this study, we measured the Young's modulus of the PCM from normal and osteoarthritic cartilage using the micropipette aspiration technique, coupled with a newly developed axisymmetric elastic layered half-space model of the experimental configuration. Viable, intact chondrons were extracted from human articular cartilage using a new microaspiration-based isolation technique. In normal cartilage, the Young's modulus of the PCM was similar in chondrons isolated from the surface zone (68.9 +/- 18.9 kPa) as compared to the middle and deep layers (62.0 +/- 30.5 kPa). However, the mean Young's modulus of the PCM (pooled for the two zones) was significantly decreased in osteoarthritic cartilage (66.5 +/- 23.3 kPa versus 41.3 +/- 21.1 kPa, p < 0.001). In combination with previous theoretical models of cell-matrix interactions in cartilage, these findings suggest that the PCM has an important influence on the stress-strain environment of the chondrocyte that potentially varies with depth from the cartilage surface. Furthermore, the significant loss of PCM stiffness that was observed in osteoarthritic cartilage may affect the magnitude and distribution of biomechanical signals perceived by the chondrocytes.     

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.