Biphasic poroviscoelastic simulation of the unconfined compression of articular cartilage: I--Simultaneous prediction of reaction force and lateral displacement

This study investigated the ability of the linear biphasic poroelastic (BPE) model and the linear biphasic poroviscoelastic (BPVE) model to simultaneously predict the reaction force and lateral displacement exhibited by articular cartilage during stress relaxation in unconfined compression. Both models consider articular cartilage as a binary mixture of a porous incompressible solid phase and an incompressible inviscid fluid phase. The BPE model assumes the solid phase is elastic, while the BPVE model assumes the solid phase is viscoelastic. In addition, the efficacy of two additional models was also examined, i.e., the transversely isotropic BPE (TIBPE) model, which considers transverse isotropy of the solid matrix within the framework of the linear BPE model assumptions, and a linear viscoelastic solid (LVE) model, which assumes that the viscoelastic behavior of articular cartilage is solely governed by the intrinsic viscoelastic nature of the solid matrix, independent of the interstitial fluid flow. It was found that the BPE model was able to accurately account for the lateral displacement, but unable to fit the short-term reaction force data of all specimens tested. The TIBPE model was able to account for either the lateral displacement or the reaction force, but not both simultaneously. The LVE model was able to account for the complete reaction force, but unable to fit the lateral displacement measured experimentally. The BPVE model was able to completely account for both lateral displacement and reaction force for all specimens tested. These results suggest that both the fluid flow-dependent and fluid flow-independent viscoelastic mechanisms are essential for a complete simulation of the viscoelastic phenomena of 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.     

Ultrasonic measurement of depth-dependent transient behaviors of articular cartilage under compression

We previously reported an ultrasound method for measuring the depth-dependent equilibrium mechanical properties of articular cartilage using quasi-static compression. The objective of this paper was to introduce our recent development for nondestructively measuring the transient depth-dependent strains of full-thickness articular cartilage specimens prepared from bovine patellae. A 50 MHz focused ultrasound transducer was used to collect ultrasound echoes from articular cartilage specimens (n=8) and sponge phantoms with open pores (n=10) during tests of compression and subsequent stress-relaxation. The transient displacements of the tissues at different depths along the compression direction were calculated from the ultrasound echoes using a cross-correlation tracking technique. An LVDT sensor and a load cell were used to measure the overall deformation of the tissue and the applied force, respectively. Results showed that the tissues inside the cartilage layer continued to move during the stress-relaxation phase after the compression was completed. In the equilibrium state, the displacements of the cartilage tissues at the depths of 1/4, 1/2, and 3/4 of the full-thickness reduced by 51%+/-22%, 54%+/-17%, and 50+/-17%, respectively, in comparison with its peak value. However, no similar phenomenon was observed in the sponge phantoms. Our preliminary results demonstrated that this ultrasound method may provide a potential tool for the nondestructive measurement of the transient depth-dependent processes involved in biological and bioengineered soft tissues as well as soft biomaterials under dynamic loading.     

Importance of collagen orientation and depth-dependent fixed charge densities of cartilage on mechanical behavior of chondrocytes

The collagen network and proteoglycan matrix of articular cartilage are thought to play an important role in controlling the stresses and strains in and around chondrocytes, in regulating the biosynthesis of the solid matrix, and consequently in maintaining the health of diarthrodial joints. Understanding the detailed effects of the mechanical environment of chondrocytes on cell behavior is therefore essential for the study of the development, adaptation, and degeneration of articular cartilage. Recent progress in macroscopic models has improved our understanding of depth-dependent properties of cartilage. However, none of the previous works considered the effect of realistic collagen orientation or depth-dependent negative charges in microscopic models of chondrocyte mechanics. The aim of this study was to investigate the effects of the collagen network and fixed charge densities of cartilage on the mechanical environment of the chondrocytes in a depth-dependent manner. We developed an anisotropic, inhomogeneous, microstructural fibril-reinforced finite element model of articular cartilage for application in unconfined compression. The model consisted of the extracellular matrix and chondrocytes located in the superficial, middle, and deep zones. Chondrocytes were surrounded by a pericellular matrix and were assumed spherical prior to tissue swelling and load application. Material properties of the chondrocytes, pericellular matrix, and extracellular matrix were obtained from the literature. The loading protocol included a free swelling step followed by a stress-relaxation step. Results from traditional isotropic and transversely isotropic biphasic models were used for comparison with predictions from the current model. In the superficial zone, cell shapes changed from rounded to elliptic after free swelling. The stresses and strains as well as fluid flow in cells were greatly affected by the modulus of the collagen network. The fixed charge density of the chondrocytes, pericellular matrix, and extracellular matrix primarily affected the aspect ratios (height/width) and the solid matrix stresses of cells. The mechanical responses of the cells were strongly location and time dependent. The current model highlights that the collagen orientation and the depth-dependent negative fixed charge densities of articular cartilage have a great effect in modulating the mechanical environment in the vicinity of chondrocytes, and it provides an important improvement over earlier models in describing the possible pathways from loading of articular cartilage to the mechanical and biological responses of chondrocytes.     

Unconfined compression of articular cartilage: nonlinear behavior and comparison with a fibril-reinforced biphasic model

Mechanical behavior of articular cartilage was characterized in unconfined compression to delineate regimes of linear and nonlinear behavior, to investigate the ability of a fibril-reinforced biphasic model to describe measurements, and to test the prediction of biphasic and poroelastic models that tissue dimensions alter tissue stiffness through a specific scaling law for time and frequency. Disks of full-thickness adult articular cartilage from bovine humeral heads were subjected to successive applications of small-amplitude ramp compressions cumulating to a 10 percent compression offset where a series of sinusoidal and ramp compression and ramp release displacements were superposed. We found all equilibrium behavior (up to 10 percent axial compression offset) to be linear, while most nonequilibrium behavior was nonlinear, with the exception of small-amplitude ramp compressions applied from the same compression offset. Observed nonlinear behavior included compression-offset-dependent stiffening of the transient response to ramp compression, nonlinear maintenance of compressive stress during release from a prescribed offset, and a nonlinear reduction in dynamic stiffness with increasing amplitudes of sinusoidal compression. The fibril-reinforced biphasic model was able to describe stress relaxation response to ramp compression, including the high ratio of peak to equilibrium load. However, compression offset-dependent stiffening appeared to suggest strain-dependent parameters involving strain-dependent fibril network stiffness and strain-dependent hydraulic permeability. Finally, testing of disks of different diameters and rescaling of the frequency according to the rule prescribed by current biphasic and poroelastic models (rescaling with respect to the sample's radius squared) reasonably confirmed the validity of that scaling rule. The overall results of this study support several aspects of current theoretical models of articular cartilage mechanical behavior, motivate further experimental characterization, and suggest the inclusion of specific nonlinear behaviors to models.     

Stress distribution and consolidation in cartilage constituents is influenced by cyclic loading and osteoarthritic degeneration

The understanding of load support mechanisms in cartilage has evolved with computational models that better mimic the tissue ultrastructure. Fibril-reinforced poroelastic models can reproduce cartilage behaviour in a variety of test conditions and can be used to model tissue anisotropy as well as assess stress and pressure partitioning to the tissue constituents. The goal of this study was to examine the stress distribution in the fibrillar and non-fibrillar solid phase and pressure in the fluid phase of cartilage in axisymmetric models of a healthy and osteoarthritic hip joint. Material properties, based on values from the literature, were assigned to the fibrillar and poroelastic components of cartilage and cancellous and subchondral compact bone regions. A cyclic load representing walking was applied for 25 cycles. Contact stresses in the fibrillar and non-fibrillar solid phase supported less than 1% of the contact force and increased only minimally with load cycles. Simulated proteoglycan depletion increased stresses in the radial and tangential collagen fibrils, whereas fibrillation of the tangential fibrils resulted in increased compressive stress in the non-fibrillar component and tensile stress in the radial fibrils. However neither had an effect on fluid pressure. Subchondral sclerosis was found to have the largest effect, resulting in increased fluid pressure, non-fibrillar compressive stress, tangential fibril stress and greater cartilage consolidation. Subchondral bone stiffening may play an important role in the degenerative cascade and may adversely affect tissue repair and regeneration treatments.     

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.     

A poroelastic finite element model of the bone–cartilage unit to determine the effects of changes in permeability with osteoarthritis

The changes experienced in synovial joints with osteoarthritis involve coupled chemical, biological, and mechanical processes. The aim of this study was to investigate the consequences of increasing permeability in articular cartilage (AC), calcified cartilage (CC), subchondral cortical bone (SCB), and subchondral trabecular bone (STB) as observed with osteoarthritis. Two poroelastic finite element models were developed using a depth-dependent anisotropic model of AC with strain-dependent permeability and poroelastic models of calcified tissues (CC, SCB, and STB). The first model simulated a bone-cartilage unit (BCU) in uniaxial unconfined compression, while the second model simulated spherical indentation of the AC surface. Results indicate that the permeability of AC is the primary determinant of the BCU’s poromechanical response while the permeability of calcified tissues exerts no appreciable effect on the force-indentation response of the BCU. In spherical indentation simulations with osteoarthritic permeability properties, fluid velocities were larger in magnitude and distributed over a smaller area compared to normal tissues. In vivo, this phenomenon would likely lead to chondrocyte death, tissue remodeling, alterations in joint lubrication, and the progression of osteoarthritis. For osteoarthritic and normal tissue permeability values, fluid flow was predicted to occur across the osteochondral interface. These results help elucidate the consequences of increases in the permeability of the BCU that occur with osteoarthritis. Furthermore, this study may guide future treatments to counteract osteoarthritis.     

Depth-dependent Compressive Equilibrium Properties of Articular Cartilage Explained by its Composition

For this study, we hypothesized that the depth-dependent compressive equilibrium properties of articular cartilage are the inherent consequence of its depth-dependent composition, and not the result of depth-dependent material properties. To test this hypothesis, our recently developed fibril-reinforced poroviscoelastic swelling model was expanded to include the influence of intra- and extra-fibrillar water content, and the influence of the solid fraction on the compressive properties of the tissue. With this model, the depth-dependent compressive equilibrium properties of articular cartilage were determined, and compared with experimental data from the literature. The typical depth-dependent behavior of articular cartilage was predicted by this model. The effective aggregate modulus was highly strain-dependent. It decreased with increasing strain for low strains, and increases with increasing strain for high strains. This effect was more pronounced with increasing distance from the articular surface. The main insight from this study is that the depth-dependent material behavior of articular cartilage can be obtained from its depth-dependent composition only. This eliminates the need for the assumption that the material properties of the different constituents themselves vary with depth. Such insights are important for understanding cartilage mechanical behavior, cartilage damage mechanisms and tissue engineering studies.     

Contribution of postnatal collagen reorientation to depth-dependent mechanical properties of articular cartilage

The collagen fibril network is an important factor for the depth-dependent mechanical behaviour of adult articular cartilage (AC). Recent studies show that collagen orientation is parallel to the articular surface throughout the tissue depth in perinatal animals, and that the collagen orientations transform to a depth-dependent arcade-like structure in adult animals. Current understanding on the mechanobiology of postnatal AC development is incomplete. In the current paper, we investigate the contribution of collagen fibril orientation changes to the depth-dependent mechanical properties of AC. We use a composition-based finite element model to simulate in a 1-D confined compression geometry the effects of ten different collagen orientation patterns that were measured in developing sheep. In initial postnatal life, AC is mostly subject to growth and we observe only small changes in depth-dependent mechanical behaviour. Functional adaptation of depth-dependent mechanical behaviour of AC takes place in the second half of life before puberty. Changes in fibril orientation alone increase cartilage stiffness during development through the modulation of swelling strains and osmotic pressures. Changes in stiffness are most pronounced for small stresses and for cartilage adjacent to the bone. We hypothesize that postnatal changes in collagen fibril orientation induce mechanical effects that in turn promote these changes. We further hypothesize that a part of the depth-dependent postnatal increase in collagen content in literature is initiated by the depth-dependent postnatal increase in fibril strain due to collagen fibril reorientation.     

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.     

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.     

Modelling cartilage mechanobiology

The growth, maintenance and ossification of cartilage are fundamental to skeletal development and are regulated throughout life by the mechanical cues that are imposed by physical activities. Finite element computer analyses have been used to study the role of local tissue mechanics on endochondral ossification patterns, skeletal morphology and articular cartilage thickness distributions. Using single-phase continuum material representations of cartilage, the results have indicated that local intermittent hydrostatic pressure promotes cartilage maintenance. Cyclic tensile strains (or shear), however, promote cartilage growth and ossification. Because single-phase material models cannot capture fluid exudation in articular cartilage, poroelastic (or biphasic) solid/fluid models are often implemented to study joint mechanics. In the middle and deep layers of articular cartilage where poroelastic analyses predict little fluid exudation, the cartilage phenotype is maintained by cyclic fluid pressure (consistent with the single-phase theory). In superficial articular layers the chondrocytes are exposed to tangential tensile strain in addition to the high fluid pressure. Furthermore, there is fluid exudation and matrix consolidation, leading to cell 'flattening'. As a result, the superficial layer assumes an altered, more fibrous phenotype. These computer model predictions of cartilage mechanobiology are consistent with results of in vitro cell and tissue and molecular biology experiments.     

A cross-validation of the biphasic poroviscoelastic model of articular cartilage in unconfined compression, indentation, and confined compression

The biphasic poroviscoelastic (BPVE) model was curve fit to the simultaneous relaxation of reaction force and lateral displacement exhibited by articular cartilage in unconfined compression (n=18). Model predictions were also made for the relaxation observed in reaction force during indentation with a porous plane-ended metal indenter (n=4), indentation with a nonporous plane ended metal indenter (n=4), and during confined compression (n=4). Each prediction was made using material parameters resulting from curve fits of the unconfined compression response of the same tissue. The BPVE model was able to account for both the reaction force and the lateral displacement during unconfined compression very well. Furthermore, model predictions for both indentation and confined compression also followed the experimental data well. These results provide substantial evidence for the efficacy of the biphasic poroviscoelastic model for articular cartilage, as no successful cross-validation of a model simulation has been demonstrated using other mathematical models.     

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

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

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

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

Fibril reinforced poroelastic model predicts specifically mechanical behavior of normal,proteoglycan depleted and collagen degraded articular cartilage

Degradation of collagen network and proteoglycan (PG) macromolecules are signs of articular cartilage degeneration. These changes impair cartilage mechanical function. Effects of collagen degradation and PG depletion on the time-dependent mechanical behavior of cartilage are different. In this study, numerical analyses, which take the compression-tension nonlinearity of the tissue into account, were carried out using a fibril reinforced poroelastic finite element model. The study aimed at improving our understanding of the stress-relaxation behavior of normal and degenerated cartilage in unconfined compression. PG and collagen degradations were simulated by decreasing the Young's modulus of the drained porous (nonfibrillar) matrix and the fibril network, respectively. Numerical analyses were compared to results from experimental tests with chondroitinase ABC (PG depletion) or collagenase (collagen degradation) digested samples. Fibril reinforced poroelastic model predicted the experimental behavior of cartilage after chondroitinase ABC digestion by a major decrease of the drained porous matrix modulus (-64+/-28%) and a minor decrease of the fibril network modulus (-11+/-9%). After collagenase digestion, in contrast, the numerical analyses predicted the experimental behavior of cartilage by a major decrease of the fibril network modulus (-69+/-5%) and a decrease of the drained porous matrix modulus (-44+/-18%). The reduction of the drained porous matrix modulus after collagenase digestion was consistent with the microscopically observed secondary PG loss from the tissue. The present results indicate that the fibril reinforced poroelastic model is able to predict specifically characteristic alterations in the stress-relaxation behavior of cartilage after enzymatic modifications of the tissue. We conclude that the compression-tension nonlinearity of the tissue is needed to capture realistically the mechanical behavior of normal and degenerated articular cartilage.     

The role of fibril reinforcement in the mechanical behavior of cartilage

Collagen fibril reinforcement was incorporated into a nonlinear poroelastic model for articular cartilage in unconfined compression. It was found that the radial fibrils play a predominant role in the transient mechanical behavior but a less important role in the equilibrium response of cartilage. The radial fibrils are in tension and can be highly stressed during compression, in contrast to low compressive stresses in all directions for the proteoglycan matrix after a small initial compression. The strain dependent fibril stiffening produces strong nonlinear transient response; the fibrils provide extra stiffness to balance a rising fluid pressure and to restrain stress increase in the proteoglycans. The fibril reinforcement, induced by the fluid pressure and flow, also accounts for a complex pattern of strain-magnitude and strain-rate dependence of cartilage stiffness.     

Acoustic properties of articular cartilage under mechanical stress

Mechano-acoustic and elastographic techniques may provide quantitative means for the in vivo diagnostics of articular cartilage. These techniques assume that sound speed does not change during tissue loading. As articular cartilage shows volumetric changes during compression, acoustic properties of cartilage may change affecting the validity of mechano-acoustic measurements. In this study, we examined the ultrasound propagation through human, bovine and porcine articular cartilage during stress-relaxation in unconfined compression. The time of flight (TOF) technique with known cartilage thickness (true sound speed) as well as in situ calibration method [Suh, Youn, Fu, J. Biomech. 34 (2001), 1347-1353] were used for the determination of sound speed. Ultrasound speed and attenuation decreased in articular cartilage during ramp compression, but returned towards the level of original values during relaxation. Variations in ultrasound speed induced an error in strain and compressive moduli provided that constant ultrasound speed and time-of-flight data was used to determine the tissue thickness. Highest errors in strain (-11.8 +/- 12.0%) and dynamic modulus (15.4 +/- 17.9%) were recorded in bovine cartilage. TOF and in situ calibration methods yielded different results for changes in sound speed during compression. We speculate that the variations in acoustic properties in loaded cartilage are related to rearrangement of the interstitial matrix, especially to that of collagen fibers. In human cartilage the changes, are, however relatively small and, according to the numerical simulations, mechano-acoustic techniques that assume constant acoustic properties for the cartilage will not be significantly impaired by this phenomenon.     

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.