The role of viscoelasticity of collagen fibers in articular cartilage: theory and numerical formulation

The relative importance of fluid-dependent and fluid-independent transient mechanical behavior in articular cartilage was examined for tensile and unconfined compression testing using a fibril reinforced model. The collagen matrix of articular cartilage was modeled as viscoelastic using a quasi-linear viscoelastic formulation with strain-dependent elastic modulus, while the proteoglycan matrix was considered as linearly elastic. The collagen viscoelastic properties were obtained by fitting experimental data from a tensile test. These properties were used to investigate unconfined compression testing, and the sensitivity of the properties was also explored. It was predicted that the stress relaxation observed in tensile tests was not caused by fluid pressurization at the macroscopic level. A multi-step tensile stress relaxation test could be approximated using a hereditary integral in which the elastic fibrillar modulus was taken to be a linear function of the fibrillar strain. Applying the same formulation to the radial fibers in unconfined compression, stress relaxation could not be simulated if fluid pressurization were absent. Collagen viscoelasticity was found to slightly weaken fluid pressurization in unconfined compression, and this effect was relatively more significant at moderate strain rates. Therefore, collagen viscoelasticity appears to play an import role in articular cartilage in tensile testing, while fluid pressurization dominates the transient mechanical behavior in compression. Collagen viscoelasticity plays a minor role in the mechanical response of cartilage in unconfined compression if significant fluid flow is present.     

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.     

Compressive fatigue and endurance of juvenile bovine articular cartilage explants

Articular cartilage is an enduring tissue. For most individuals, articular cartilage facilitates a lifetime of pain-free ambulation, supporting millions of loading cycles from activities of daily living. Although early studies into osteoarthritis focused on the role of mechanical fatigue in articular cartilage degeneration, much is still unknown regarding its strength and endurance characteristics. The compressive strength of juvenile, bovine articular cartilage explants was determined to be loading rate-dependent, reaching a maximum strength of 29.5 ± 4.8 MPa at a strain rate of 0.10 %/sec. The fatigue and endurance properties of articular cartilage were characterized utilizing a material testing system, as well as a custom, validated instrument termed the two degrees-of-freedom endurance meter (endurometer). These instruments characterized fatigue in articular cartilage explants at loading levels ranging from 10 to 80 % strength (%S), up to 100,000 cycles. Cartilage explants displayed characteristics of fatigue – fatigue life increased as the loading magnitude decreased. All explants failed within 14,000 cycles at loading levels between 50 and 80 %S. At 10 and 20 %S, all explants endured 100,000 loading cycles. There was no significant difference in equilibrium compressive modulus between run-out explants and unloaded controls, although the pooled modulus increased in response to testing. Histological staining and biochemical assays revealed no material changes in collagen, sulfated glycosaminoglycan, or hydration content between unloaded controls and explants cyclically loaded to run-out. These results suggest articular cartilage may have a putative endurance limit of 20 %S (5.86 MPa), with implications for articular cartilage biomechanics and mechanobiology.     

Time-dependent mechanical behavior of sheep digital tendons, including the effects of preconditioning

The time-dependent mechanical properties of sheep digital extensor tendons were studied by sequences of stress-relaxation tests. The results exhibited irreversible preconditioning and reversible viscoelasticity. Preconditioning effects were manifested by stress decay during consecutive stretch cycles to the same strain level, accompanied by elongation of the tendon's reference length. They intensified with increased strain level, and were reduced or became negligible as the strain decreased. The significance of intrinsic response mechanisms was studied via a structural model that includes viscoelasticity, preconditioning, and morphology of the tendon's collagen fibers. Model/data comparisons showed good agreement and good predictive power, suggesting that preconditioning can be integrated into comprehensive material characterization of tendons.     

Experimental verification of the roles of intrinsic matrix viscoelasticity and tension-compression nonlinearity in the biphasic response of cartilage

A biphasic-CLE-QLV model proposed in our recent study [2001, J. Biomech. Eng., 123, pp. 410-417] extended the biphasic theory of Mow et al. [1980, J. Biomech. Eng., 102, pp. 73-84] to include both tension-compression nonlinearity and intrinsic viscoelasticity of the cartilage solid matrix by incorporating it with the conewise linear elasticity (CLE) model [1995, J. Elasticity, 37, pp. 1-38] and the quasi-linear viscoelasticity (QLV) model [Biomechanics: Its foundations and objectives, Prentice Hall, Englewood Cliffs, 1972]. This model demonstrates that a simultaneous prediction of compression and tension experiments of articular cartilage, under stress-relaxation and dynamic loading, can be achieved when properly taking into account both flow-dependent and flow-independent viscoelastic effects, as well as tension-compression nonlinearity. The objective of this study is to directly test this biphasic-CLE-QLV model against experimental data from unconfined compression stress-relaxation tests at slow and fast strain rates as well as dynamic loading. Twelve full-thickness cartilage cylindrical plugs were harvested from six bovine glenohumeral joints and multiple confined and unconfined compression stress-relaxation tests were performed on each specimen. The material properties of specimens were determined by curve-fitting the experimental results from the confined and unconfined compression stress relaxation tests. The findings of this study demonstrate that the biphasic-CLE-QLV model is able to describe the strain-rate-dependent mechanical behaviors of articular cartilage in unconfined compression as attested by good agreements between experimental and theoretical curvefits (r2 = 0.966 +/- 0.032 for testing at slow strain rate; r2 = 0.998 +/- 0.002 for testing at fast strain rate) and predictions of the dynamic response (r2 = 0.91 +/- 0.06). This experimental study also provides supporting evidence for the hypothesis that both tension-compression nonlinearity and intrinsic viscoelasticity of the solid matrix of cartilage are necessary for modeling the transient and equilibrium responses of this tissue in tension and compression. Furthermore, the biphasic-CLE-QLV model can produce better predictions of the dynamic modulus of cartilage in unconfined dynamic compression than the biphasic-CLE and biphasic poroviscoelastic models, indicating that intrinsic viscoelasticity and tension-compression nonlinearity of articular cartilage may play important roles in the load-support mechanism of cartilage under physiologic loading.     

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.     

Biphasic poroviscoelastic simulation of the unconfined compression of articular cartilage: II--Effect of variable strain rates

This study investigated the abilities of the linear biphasic poroviscoelastic (BPVE) model and the linear biphasic poroelastic (BPE) model to simulate the effect of variable ramp strain rates on the unconfined compression stress relaxation response of articular cartilage. Curve fitting of experimental data showed that the BPVE model was able to successfully account for the ramp strain rate-dependent viscoelastic behavior of articular cartilage under unconfined compression, while the BPE model was able to account for the complete viscoelastic response at a slow strain rate, but only the long-term viscoelastic response at faster strain rates. We concluded that the short-term viscoelastic behavior of articular cartilage, when subjected to a fast ramp strain rate, is primarily governed by a fluid flow-independent (intrinsic) viscoelastic mechanism, whereas the long-term viscoelastic behavior is governed by a fluid flow-dependent (biphasic) viscoelastic mechanism. Furthermore, a linear viscoelastic representation of the solid stress was found to be a valid model assumption for the simulation of ramp strain rate-dependent relaxation behaviors of articular cartilage within the range of ramp strain rates investigated.     

Preconditioning is correlated with altered collagen fiber alignment in ligament

Although the mechanical phenomena associated with preconditioning are well-established, the underlying mechanisms responsible for this behavior are still not fully understood. Using quantitative polarized light imaging, this study assessed whether preconditioning alters the collagen fiber alignment of ligament tissue, and determined whether changes in fiber organization are associated with the reduced force and stiffness observed during loading. Collagen fiber alignment maps of facet capsular ligaments (n?=?8) were generated before and after 30 cycles of cyclic tensile loading, and alignment vectors were correlated between the maps to identify altered fiber organization. The change in peak force and tangent stiffness between the 1st and 30th cycle were determined from the force-displacement response, and the principal strain field of the capsular ligament after preconditioning was calculated from the fiber alignment images. The decreases in peak ligament force and tangent stiffness between the 1st and 30th cycles of preconditioning were significantly correlated (R ≥ 0.976, p?相似文献   

Alterations in mechanical behaviour of articular cartilage due to changes in depth varying material properties--a nonhomogeneous poroelastic model study

The depth dependence of the material properties is present in normal adult cartilage and is believed to have significant implications in its normal mechanical function. Cartilage pathology may alter the depth dependence, e.g. a reduced depth dependence of the fibril stiffness has been observed in osteoarthritic cartilage. The objective of the present study is to investigate the alterations in the mechanical response of articular cartilage when the depth dependence of the material properties is varied to simulate healthy and pathological situations. This study is made possible by a recently developed nonhomogeneous poroelastic model. Depth variations of the strains and stresses for individual material phases (collagen, proteoglycan and fluid) are obtained for cartilage disks in unconfined compression using the finite element method. The mean nominal axial strain considered is up to 15%, while the axial strain at the articular surface can reach 33%. This paper demonstrates how the mechanical behaviours of cartilage are affected by individual depth dependent cartilage properties, while such observations are not fully available in experimental investigations. This study suggests the possibility of diagnosing cartilage health by analysing its mechanical behaviours.     

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.     

Determination of the Poisson's ratio of the cell: recovery properties of chondrocytes after release from complete micropipette aspiration

Chondrocytes in articular cartilage are regularly subjected to compression and recovery due to dynamic loading of the joint. Previous studies have investigated the elastic and viscoelastic properties of chondrocytes using micropipette aspiration techniques, but in order to calculate cell properties, these studies have generally assumed that cells are incompressible with a Poisson's ratio of 0.5. The goal of this study was to measure the Poisson's ratio and recovery properties of the chondrocyte by combining theoretical modeling with experimental measures of complete cellular aspiration and release from a micropipette. Chondrocytes isolated from non-osteoarthritic and osteoarthritic cartilage were fully aspirated into a micropipette and allowed to reach mechanical equilibrium. Cells were then extruded from the micropipette and cell volume and morphology were measured throughout the experiment. This experimental procedure was simulated with finite element analysis, modeling the chondrocyte as either a compressible two-mode viscoelastic solid, or as a biphasic viscoelastic material. By fitting the experimental data to the theoretically predicted cell response, the Poisson's ratio and the viscoelastic recovery properties of the cell were determined. The Poisson's ratio of chondrocytes was found to be 0.38 for non-osteoarthritic cartilage and 0.36 for osteoarthritic chondrocytes (no significant difference). Osteoarthritic chondrocytes showed an increased recovery time following full aspiration. In contrast to previous assumptions, these findings suggest that chondrocytes are compressible, consistent with previous studies showing cell volume changes with compression of the extracellular matrix.     

The role of flow-independent viscoelasticity in the biphasic tensile and compressive responses of articular cartilage

A long-standing challenge in the biomechanics of connective tissues (e.g., articular cartilage, ligament, tendon) has been the reported disparities between their tensile and compressive properties. In general, the intrinsic tensile properties of the solid matrices of these tissues are dictated by the collagen content and microstructural architecture, and the intrinsic compressive properties are dictated by their proteoglycan content and molecular organization as well as water content. These distinct materials give rise to a pronounced and experimentally well-documented nonlinear tension-compression stress-strain responses, as well as biphasic or intrinsic extracellular matrix viscoelastic responses. While many constitutive models of articular cartilage have captured one or more of these experimental responses, no single constitutive law has successfully described the uniaxial tensile and compressive responses of cartilage within the same framework. The objective of this study was to combine two previously proposed extensions of the biphasic theory of Mow et al. [1980, ASME J. Biomech. Eng., 102, pp. 73-84] to incorporate tension-compression nonlinearity as well as intrinsic viscoelasticity of the solid matrix of cartilage. The biphasic-conewise linear elastic model proposed by Soltz and Ateshian [2000, ASME J. Biomech. Eng., 122, pp. 576-586] and based on the bimodular stress-strain constitutive law introduced by Curnier et al. [1995, J. Elasticity, 37, pp. 1-38], as well as the biphasic poroviscoelastic model of Mak [1986, ASME J. Biomech. Eng., 108, pp. 123-130], which employs the quasi-linear viscoelastic model of Fung [1981, Biomechanics: Mechanical Properties of Living Tissues, Springer-Verlag, New York], were combined in a single model to analyze the response of cartilage to standard testing configurations. Results were compared to experimental data from the literature and it was found that a simultaneous prediction of compression and tension experiments of articular cartilage, under stress-relaxation and dynamic loading, can be achieved when properly taking into account both flow-dependent and flow-independent viscoelasticity effects, as well as tension-compression nonlinearity.     

Modeling of constrained articular cartilage growth in an intact knee with focal knee resurfacing metal implant

The purpose of the present study was to develop a model to simulate the articular cartilage growth in an intact knee model with a metal implant replacing a degenerated portion of the femoral cartilage. The human knee joint was approximated with a simplified axisymmetric shape of the femoral condyle along with the cartilage, meniscus and bones. Two individually growing constituents (proteoglycans and collagen) bound to solid matrix were considered in the solid phase of the cartilage. The cartilage behavior was modeled with a nonlinear biphasic porohyperelastic material model, and meniscus with a transversely isotropic linear biphasic poroelastic material model. Two criteria (permeation and shear), both driven by mechanical loading, were considered to trigger the growth in the solid constituents. Mechanical loading with sixty heavy cycles was considered to represent daily walking activity. The growth algorithm was implemented for 90 days after implantation. The results from simulations show that both cartilage layers were more stimulated near the implant which lead to more growth of the cartilage near the defect. The method developed in the present work could be a powerful technique if more accurate material data and growth laws were available.     

Tendon exhibits complex poroelastic behavior at the nanoscale as revealed by high-frequency AFM-based rheology

Tendons transmit load from muscle to bone by utilizing their unique static and viscoelastic tensile properties. These properties are highly dependent on the composition and structure of the tissue matrix, including the collagen I hierarchy, proteoglycans, and water. While the role of matrix constituents in the tensile response has been studied, their role in compression, particularly in matrix pressurization via regulation of fluid flow, is not well understood. Injured or diseased tendons and tendon regions that naturally experience compression are known to have alterations in glycosaminoglycan content, which could modulate fluid flow and ultimately mechanical function. While recent theoretical studies have predicted tendon mechanics using poroelastic theory, no experimental data have directly demonstrated such behavior. In this study, we use high-bandwidth AFM-based rheology to determine the dynamic response of tendons to compressive loading at the nanoscale and to determine the presence of poroelastic behavior. Tendons are found to have significant characteristic dynamic relaxation behavior occurring at both low and high frequencies. Classic poroelastic behavior is observed, although we hypothesize that the full dynamic response is caused by a combination of flow-dependent poroelasticity as well as flow-independent viscoelasticity. Tendons also demonstrate regional dependence in their dynamic response, particularly near the junction of tendon and bone, suggesting that the structural and compositional heterogeneity in tendon may be responsible for regional poroelastic behavior. Overall, these experiments provide the foundation for understanding fluid-flow-dependent poroelastic mechanics of tendon, and the methodology is valuable for assessing changes in tendon matrix compressive behavior at the nanoscale.     

The possible role of poroelasticity in the apparent viscoelastic behavior of passive cardiac muscle

This paper investigates the contribution of extracellular fluid flow to the apparent viscoelastic behavior of passive cardiac muscle. The muscle is modeled as an incompressible, isotropic, poroelastic solid saturated by an incompressible viscous fluid. Based on Biot's linear and nonlinear consolidation theories, solutions are presented for general time-dependent uniaxial loading of unconfined cylindrical muscle specimens. The nonlinear analysis includes the effects of large strain, material nonlinearity, and strain-dependent permeability. The computed results show that, for axial stretch ratios greater than 1.1, the changing permeability and the loading rate strongly affect the total stress relaxation and the short-time relaxation rate. Comparisons of theoretical and published experimental results show that extracellular fluid flow can account for several observed biomechanical features of passive myocardium, including the insensitivity of stress-strain curves to loading rate and of stress-relaxation curves to the amount of stretch. Theoretical hysteresis loops, however, are too small. Thus, both poroelastic and tissue viscoelastic effects must be considered in studies of passive cardiac muscle.     

Contribution of tissue composition and structure to mechanical response of articular cartilage under different loading geometries and strain rates

Mechanical function of articular cartilage in joints between articulating bones is dependent on the composition and structure of the tissue. The mechanical properties of articular cartilage are traditionally tested in compression using one of the three loading geometries, i.e., confined compression, unconfined compression or indentation. The aim of this study was to utilize a composition-based finite element model in combination with a fractional factorial design to determine the importance of different cartilage constituents in the mechanical response of the tissue, and to compare the importance of the tissue constituents with different loading geometries and loading rates. The evaluated parameters included water and collagen fraction as well as fixed charge density on cartilage surface and their slope over the tissue thickness. The thicknesses of superficial and middle zones, as based on the collagen orientation, were also included in the evaluated parameters. A three-level resolution V fractional factorial design was used. The model results showed that inhomogeneous composition plays only a minor role in indentation, though that role becomes more significant in confined compression and unconfined compression. In contrast, the collagen architecture and content had a more profound role in indentation than with two other loading geometries. These differences in the mechanical role of composition and structure between the loading geometries were emphasized at higher loading rates. These findings highlight how the results from mechanical tests of articular cartilage under different loading conditions are dependent upon tissue composition and structure.     

Collagen Network of Articular Cartilage Modulates Fluid Flow and Mechanical Stresses in Chondrocyte

The extracellular matrix of articular cartilage modulates the mechanical signals sensed by the chondrocytes. In the present study, a finite element model (FEM) of the chondrocyte and its microenvironment was reconstructed using the information from fourier transform infrared imaging spectroscopy. This environment consisted of pericellular, territorial (mainly proteoglycans), and inter-territorial (mainly collagen) matrices. The chondrocyte, pericellular, and territorial matrix were assumedto be mechanically isotropic and poroelastic, whereas the inter-territorial matrix, due to its high collagen content, was assumed to be transversely isotropic and poroelastic. Under instantaneous strain-controlled compression, the FEM indicated that the fluid pressure within the chondrocyte increased nonlinearly as a function of the in-plane Young’s modulus of the collagen network. Under instantaneous force-controlled compression, the chondrocyte experienced the highest fluid pressure when the in-plane Young’s modulus of the collagen network was ~4 MPa. Based on the present results, the mechanical characteristics of the collagen network of articular cartilage can modify fluid flow and stresses in chondrocytes. Therefore, the integrity of the collagen network may be an important determinant in cell stimulation and in the control of the matrix maintenance.