Skip Navigation

This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (11)
Right arrowRequest Permissions
Right arrow Disclaimer
Google Scholar
Right arrow Articles by Barker, M. K.
Right arrow Articles by Seedhom, B. B.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Barker, M. K.
Right arrow Articles by Seedhom, B. B.
Related Collections
Right arrow Osteoarthritis and Cartilage
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

Rheumatology 2001; 40: 274-284
© 2001 British Society for Rheumatology

The relationship of the compressive modulus of articular cartilage with its deformation response to cyclic loading: does cartilage optimize its modulus so as to minimize the strains arising in it due to the prevalent loading regime?

M. K. Barker and B. B. Seedhom

Rheumatology and Rehabilitation Research Unit, University of Leeds, 36 Clarendon Road, Leeds LS2 9NZ, UK


    Abstract
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 Conclusions
 References
 
Aim. To investigate the relationship of the instantaneous compressive modulus with its deformation response to cyclic loading typical of that encountered at the knee joint during level walking.

Method. The study was performed on 24 osteochondral plugs taken from three unembalmed cadaveric knees. As the compressive modulus of cartilage has been shown to vary topographically across the knee in an established manner, the specimens were taken from specific sites on the femur and tibia of each knee. All the cartilage specimens were immersed in Hanks’ salt solution at 37°C and were subjected to the same cyclic loading regimen that was representative of a typical walking cycle in a specialized indentation apparatus, for over 1 h.

Results and conclusion. The viscous and elastic components of matrix strain, the creep rate and the cartilage compressive modulus were measured. The latter was found to be significantly related to the strain response of cartilage to cyclic loading. Elastic strain varied exponentially with the compressive modulus; specimens with a modulus less than 4 MPa experienced elastic strains in the range 0.18–0.36, whereas stiffer specimens experienced strains between 0.05 and 0.13. Viscous strain varied linearly with cartilage stiffness and was as low as 0.02 at the lower values of the compressive modulus but increased to 0.22 for a compressive modulus of 18 MN/m2. The rate of creep under cyclic load was inversely linearly related to cartilage stiffness. The strain response of soft specimens approached steady state by 200 cycles but that of stiff specimens did not approach it until 1300 cycles. It was hypothesized that the viscous strain response of cartilage can be explained in terms of differences in permeability between specimens of different compressive modulus, stiffer cartilage having a lower permeability than soft cartilage.

KEY WORDS: Articular cartilage, Compressive modulus, Cyclic loading, Strain, Fluid flow.


    Introduction
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 Conclusions
 References
 
The deformation of cartilage has been studied for the most part under static loading, and the various studies have shown that cartilage exhibits viscoelastic behaviour [15]. The viscoelasticity of cartilage tissue has been explained in terms of interstitial fluid flow within the cartilage matrix and of the intrinsic viscoelasticity of the matrix itself [69].

During locomotion and exercise, cartilage is subjected to cyclic loading and the loads are applied at fast rates, the rise times being 10–150 ms [10], and act for short durations (10–400 ms) [11, 12, 34]. The deformation response of cartilage under cyclic loading conditions has received limited attention. Basic experimental measurements have been attempted [1419], while various models have been used to predict cyclic strain behaviour [2022]. Maroudas [23] made predictions as to the likely deformation response of cartilage under cyclic loading conditions. She postulated that recovery between load cycles is likely to be incomplete because the recovery rate would be governed by the swelling pressure, which would be much lower than the physiological stress caused by the cyclic load that is acting.

The deformation of cartilage and rate of flow of interstitial fluid through its matrix must each be a function of the parameters of the duty cycle as well as of cartilage properties. Parameters of the duty cycle include: (i) the loading rate; (ii) the frequency of load application, (iii) the amplitude of the applied stress; (iv) the ratio of the loading period to the recovery period within a loading cycle; and (iv) the number of loading cycles. The property of cartilage most relevant to its deformation would be its compressive modulus.

This study focuses on the effect of the compressive modulus on cartilage deformation under cyclic loading conditions—particularly those arising during level walking, which is the predominant human activity. One important reason for focusing this investigation on the effect of the compressive modulus on cartilage deformation rather than on the other parameters mentioned above is the greater variation in the magnitude of the compressive modulus of cartilage compared with the other parameters. The modulus ranges from 1 to 20 MN/m2, whereas the frequency of loading falls within a narrow band, between 1 Hz (during slow walking) and 2.25 Hz (during sprinting). Furthermore, the average stress in the knees during walking is 1–1.5 MPa.

It is interesting that the relationship between cartilage deformation and its modulus has not been subject to any detailed experimental study. This might well be due to an expectation that such a relationship would obviously be of an inverse, perhaps monotonic nature and therefore in itself would not be of interest. However, since cartilage viscoelastic behaviour, as mentioned earlier, is contributed to by the intrinsic viscoelastic nature of the cartilage matrix and by interstitial fluid flow in undefined proportions, cartilage behaviour under cyclic loading is complex and may not be as predictable as was at first thought.


    Materials and methods
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 Conclusions
 References
 
Materials
We used eight cartilage specimens from each of three unembalmed human cadaveric knee joints. When these had been procured they were immediately stored at -20°C, and before use they were thawed overnight at room temperature. Each knee was dissected to remove all soft tissue and to expose the cartilage surfaces of the femoral condyles and tibial plateaux. Meachim's Indian ink test [24] was performed to identify any areas of fibrillated cartilage, and only healthy areas of cartilage were selected for testing. Osteochondral plugs 12 mm in diameter and 12 mm long were removed for testing by the use of a specially designed reamer. Eight such specimens were harvested from each knee joint from specific sites on the femur and tibia, as illustrated in Fig. 1Go. Cartilage specimens from these sites have been shown to have very different stiffnesses [25, 26], which allowed the effect of cartilage stiffness on its deformation response to be studied.



View larger version (12K):
[in this window]
[in a new window]
 
FIG. 1. Locations and labels of the eight specimens taken from each knee joint. The prefixes L and M signify lateral and medial respectively. PS, patellar surface of the femur; FC, femoral condyle; TC, areas of cartilage on the tibial plateau covered by the menisci; TU, areas of cartilage on the tibial plateau which come into direct contact with the femur. Previous studies have shown that FC sites are significantly stiffer than PS sites, and TC sites are significantly stiffer than TU sites.

 

Methods
Specimen alignment.
The osteochondral specimen, which was 12 mm in height, was secured in a specimen- holder with dental cement (Fig. 2aGo), leaving the cartilage layer exposed and protruding above the upper surface of the specimen-holder (the cartilage was thus unconstrained at its edges). The cartilage specimen was then positioned beneath the tip of the indenter in the test apparatus. Using both the naked eye and an arthroscopic camera (and monitor) to view the indenter tip and cartilage surface, fine adjustment of the alignment was made using three screws on the base of the specimen-holder (Fig. 2bGo). Throughout the alignment procedure, the cartilage specimen was kept moist by irrigating it with Hanks’ balanced salt solution (HBSS). The accuracy of the alignment technique was better than 1°. Once perpendicular alignment of the cartilage surface to the indenter had been achieved, a heating stage was fixed to the specimen-holder. This heating stage also formed a chamber around the cartilage specimen, so that it could be immersed in HBSS maintained at 37°C by circulating heated water in the heating element.



View larger version (55K):
[in this window]
[in a new window]
 
FIG. 2. Section diagrams of the specimen-holder and the alignment procedure. (1) Adjustment screw; (2) base plate; (3) specimen holder; (4) dental cement; (5) attachment to cyclic indentation rig; (6) plane-ended impervious indenter; (7) osteochondral specimen; (8) arthroscope; (9) specimen plate; (10) hole for pushing specimen to remove it; (11) clamp and nut; (12) heating element; (13) heating element fluid; (14) HBSS; (15) O-ring seal; (16) inlet and outlet ports for heating fluid; (17) specimen bone substrate. (a) The cartilage is secured in dental cement so that it protrudes entirely above the upper surface of the specimen-holder, thus allowing visualization of the indenter tip and cartilage surface for the purpose of alignment. (b) After alignment of the specimen, the heating stage is fitted; this provides a chamber for saline. After testing, hole 10 is used as access to drive the specimen out of the holder, while plate 9 prevents the tool from damaging the specimen.

 

Measurement of cartilage modulus.
The stiffness of each specimen was measured in terms of its 50 ms compressive modulus, which is representative of the physiological loading rate encountered during walking. It is based on the deformation of cartilage 50 ms after initial contact of the indenter tip with the cartilage surface. The calculation of the modulus required three measurements to be obtained: the cartilage deformation after 50 ms of loading, the applied load after 50 ms of loading and the cartilage thickness. The indentation and load measurements were performed twice on each specimen, once before and once after performing the cyclic loading experiment. This was done in order to check if the repetitive loading had damaged the specimens. As the rate of loading in these modulus measurements was high (rise time 50 ms or less) and fluid flow (if any) would not be significant within such intervals (the measurement was undertaken in a single indentation), it was justified to treat cartilage as an elastic material and to use an impervious indenter. This latter was hemispherical and had a radius of 1.5 mm. It is appropriate to comment on the effect of the curvature of the cartilage of the various osteochondral plugs tested on the calculated values of the cartilage modulus. The equation from which such values are calculated assumes contact between a hemispherical rigid indenter and a flat cartilage surface. The curvature of the specimens will vary from 20 mm on the femoral condyles to 70 mm on the tibial plateaux (the medial being concave and the lateral convex). At the most, the equivalent radius based on the measurements would be lower by 7.5% for a surface of 20 mm radius, such as the femoral condyle, and by 2% for a surface of a 70 mm radius, such as that of a tibial plateau. Such variations in contact geometry will result in a small error in the calculated value of the compressive modulus. The test parameters were identical to those used in two previous studies from our laboratory [26, 27], which induced physiological stresses.

The cartilage thickness was measured using a technique that employed a needle, which caused disruption of the surface of the cartilage specimen. For this reason, the thickness measurement was left until all other testing had been completed so as not to violate the cartilage surface. This technique is described by Swann and Seedhom [25].

Deformation response under cyclic loading.
This was investigated immediately after we had measured the indentation for the determination of the modulus. Therefore, before the cyclic load was applied, the specimen was left for 30 min to ensure full recovery from the indentation caused during the modulus measurement. The hemispherical indenter used for this measurement was then replaced with a plane-ended impervious indenter of the same diameter (3 mm), in order to ensure that all cartilage specimens were subjected to the same contact stress. The duty cycle parameters were as follows: amplitude of applied stress 1.4 MPa (typical of the stress occurring during walking); frequency of load application 1 Hz; load rise time 20 ms; and ratio of loading duration to recovery duration 1:2 (loading period 330 ms, recovery period 670 ms) (Fig. 3Go). Load and displacement traces were recorded simultaneously during each of the first 64 consecutive cycles and thereafter during one complete cycle every 64 cycles, for a little over 1 h (these numbers were dictated by the electronics of the data acquisition system). In total, data were collected from 4160 loading cycles.



View larger version (16K):
[in this window]
[in a new window]
 
FIG. 3. Profile of a single load cycle applied to the cartilage to simulate physiological aspects of the walking cycle. The total cycle duration was 1 s. The duration of the loading phase was 330 ms and that of the recovery phase (during which the indenter was lifted clear of the cartilage surface) was 670 ms. This gave a loading-to-recovery ratio of 1:2, which is typical of the loading experienced by many sites on the articular surface of the femur. The load rise time was 20 ms (typical of the fast loading rates experienced physiologically), and the peak stress was 1.4 MPa (a stress level typically encountered at the knee joint during walking). A total of 4160 of these cycles were applied consecutively to each cartilage specimen in each experiment.

 

Data analysis
Strain components.
The various load and displacement traces were then analysed using a Fortran programme producing two sets of data; one set consisted of the position of the cartilage surface at the instant of contact between indenter and cartilage at the beginning of each loading cycle. The other data set consisted of the position of the cartilage surface prior to the instant of cessation of contact between cartilage and the indenter, for each cycle, as described by Barker and Seedhom [13]. A typical set of these data points is shown in Fig. 4aGo and is described here in order to define the variables analysed. The curves exhibited two characteristic phases of the cartilage deformation under cyclic loading. The first phase occurs within the first 1000 cycles and the data show that the cartilage surface does not fully recover its deformation between individual loading cycles. This implies that more water is being exuded during the loading phase than is being imbibed during the load-free recovery phase. The curve is steepest during the first 200 loading cycles but the gradient reduces with the number of cycles, indicating a reduction in the amount of fluid being lost. The second phase occurs from 1000 cycles onwards. The curve heads asymptotically towards a steady state, during which cartilage is being depressed and then allowed to recover fully during each cycle. Analysis of the load and displacement during this steady state shows that there are two processes occurring. The first process, which accounts for most of the deformation, is that of instantaneous deformation of the solid matrix during the loading phase and an equal, instantaneous recovery during the recovery phase. These deformations are equal and opposite after the specimen has reached a steady state. The second process, which accounts for a small fraction of the total deformation, is that of fluid exudation during loading and fluid uptake during recovery. These small volumes of fluid flow are also equal. Hence, the overall deformation response at this steady state is termed ‘elastic’, because the cartilage surface position at the beginning of each cycle (i.e. when the indenter comes into contact with the cartilage surface) remains unchanged.



View larger version (25K):
[in this window]
[in a new window]
 
FIG. 4. (a) Definition of viscous, elastic and total strain data for specimen A_LPS. This specimen was 2.8 mm thick and had a 50 ms compressive modulus of 9.4 MN/m2. (b) Exponential curve fitted to the total strain data in Fig. 4aGo, together with the differential of the strain curve with respect to the number of loading cycles. The diagram illustrates how the gradient value of -0.003 was used to define the number of cycles to approach the steady-state condition.

 
The cyclic strain thus comprised two components: a viscous component and an elastic one (Fig. 4aGo). The higher curve represents the strain that is purely due to the cumulative loss of fluid from the loaded region, which has been termed ‘viscous strain’ ({varepsilon}viscous). The difference between the upper and lower curves is the elastic strain of the solid cartilage matrix ({varepsilon}elastic). The total matrix strain ({varepsilon}total) experienced during a particular loading cycle, which is the sum of the viscous and elastic strain components, is represented by the lower of the two curves.

Cartilage creep.
The creep rate of each specimen throughout the test was also investigated. A two-term exponential curve (equation 1) was fitted to the total strain curve ({varepsilon}total) of each specimen (Fig. 4bGo). The two data curves were very similar in shape. However, the lower one was used for the curve fit because the points generally displayed a smoother curve. The slight variation in the points of the upper curve was due to the inaccuracies that arise in measuring the surface position of the cartilage at such high approach velocities. These inaccuracies were not present on the lower curve as the indenter is virtually static at this point. The form of the curve fit was:


(1)
where {varepsilon}total is total tissue strain, x is the number of cycles, and C0, C1, C2, t1 and t2 are constants. For each specimen, equation 1 was differentiated to investigate the gradient of the strain curve and hence the rate of cumulative strain:


(2)

To obtain a measure of the number of load cycles applied to the specimen before reaching steady state, equation 1 was differentiated to yield equation 2. Solving equation 2 for C3=0, the value of x, provides a measure of the number of cycles applied until the gradient of the fitted creep curve (1) became zero (or steady state). Some samples were still undergoing creep after the 1 h of loading applied during this test series, so C3=0 could not be used. Instead, a value of C3=-0.003 was found to provide a solution across all curve fits and represented a gradient of creep curve very nearly approaching the steady state. The solution x was in fact the cycle number (Nss) at which this gradient occurred and was therefore calculated for each specimen. Figure 4bGo shows the exponential curve, which was fitted to the {varepsilon}total data in Fig. 4aGo. The differential of this curve is also displayed. In this case it took 600 cycles before the steady state was reached.


    Results
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 Conclusions
 References
 
Compressive modulus values before and after cyclic loading tests
The compressive modulus for each of these specimens had a value that was consistent with the general pattern of topographical variation in knee cartilage modulus reported previously [25, 26, 28]. Thus, specimens from the femoral condyles were stiffer than specimens from the patellar surface of the femur. Likewise, specimens from the areas of the tibial plateau covered by the meniscus were stiffer than those from the areas that come into direct contact with the femur. In all cases the specimens from the area of the tibial plateau, which come into direct contact with the femur, were the softest. The modulus had a range of values between 1 and 19.5 MN/m2.

Figure 5Go and Table 1Go show the 50 ms compressive modulus values for each specimen before and after completion of the cyclic loading regimen. In one case the specimen had a significantly higher compressive modulus after testing, indicating that some tissue damage may have occurred to it during the cyclic loading, and this specimen was therefore excluded from further analysis. The sample was probably damaged as a result of the abnormal stress gradients imposed at the edge of the flat-ended indenter. The ‘before’ and ‘after’ 50 ms compressive modulus values determined for all the other specimens fell within the bounds of accuracy of the apparatus, and these specimens were therefore assumed to have been undamaged by the cyclic loading. Therefore 23 specimens were included in this study.



View larger version (22K):
[in this window]
[in a new window]
 
FIG. 5. The 50 ms compressive modulus value after completion of the cyclic loading testing regimen plotted against its value before cyclic loading. All points fell within the measurement accuracy of the apparatus (dotted line), except for sample A_MFC (asterisk), which was significantly stiffer after cyclic loading. This specimen was assumed to be damaged after testing and was excluded from further analysis.

 

View this table:
[in this window]
[in a new window]
 
TABLE 1. Data obtained after cyclic load testing for each of the 24 specimens

 

Cartilage response to cyclic loading
Stiff and soft cartilage.
A comparison between typical deformation responses of soft and stiff cartilage to cyclic loading is shown in Fig. 6Go. A specimen with a high compressive modulus of 16.8 MN/m2 (Fig. 6aGo) is compared with a specimen of low compressive modulus of 3.1 MN/m2 (Fig. 6bGo). The viscous and total strain curves are plotted for each of these two specimens, together with their respective curve fits. The elastic strain is the difference between these two curves. The diagram highlights the significant differences in the cyclic deformation responses of the two specimens. The viscous strain of the softer cartilage is less than that of the stiffer cartilage, whereas its elastic strain during each cycle is much greater. Furthermore, the softer cartilage reaches its steady-state deformation sooner than the stiff cartilage. The results comparing the strain responses of all specimens demonstrate this behaviour consistently.



View larger version (12K):
[in this window]
[in a new window]
 
FIG. 6. Comparison of typical cyclic strain responses of (a) a specimen with a high 50 ms compressive modulus of 16.8 MN/m2, and (b) a softer specimen with a 50 ms compressive modulus of 3.1 MN/m2. Compared with the stiffer cartilage, the softer specimen had higher elastic strain and less viscous strain, and reached steady state sooner. Equations of the best fit curves to these data are as follows: (a) {varepsilon}visc=0.17-0.051e(-x/50)-0.11e(-x/847) and {varepsilon}total=0.24-0.049e(-x/80)-0.1e(-x/949); (b) {varepsilon}visc=0.096-0.051e(-x/110)-0.042e(-x/1405) and {varepsilon}total=0.32-0.07e(-x/85)-0.039e(-x/984).

 

Strain components and creep.
Best-fit curves were applied to the elastic strain, viscous strain and creep rate data. The elastic strain data were best fitted by an exponential curve, while the viscous strain and creep rate data were best fitted by linear regression. The equations of all the curve fits, together with their associated Pearson coefficients and significance values, are shown in Table 2Go. The Pearson coefficient and significance value associated with the exponential curve were determined by plotting the natural logarithm of the strain values against the compressive modulus to give a linear relationship. As shown in Table 2Go, the P values were all less than 0.001, indicating strong, significant correlations. The stiffness of the cartilage, expressed in terms of its 50 ms (instantaneous) compressive modulus, was thus significantly correlated to both the elastic and viscous strains, and also to the creep rate. As the compressive modulus increased, the elastic strain decreased exponentially, heading towards an asymptotic value of elastic strain between 0.05 and 0.1. Softer specimens with a compressive modulus less than 4 MN/m2 experienced elastic strains in the range 0.18–0.36, while stiffer specimens experienced strains between 0.05 and 0.13. The viscous strains increased linearly with the 50 ms compressive modulus. Soft specimens exhibited low viscous strains (0.02), while those of stiff specimens were up to 0.22. The creep rate was inversely and linearly related to the 50 ms compressive modulus. Soft specimens approached steady state in as few as 200 cycles, while stiff specimens took up to 1300 cycles.


View this table:
[in this window]
[in a new window]
 
TABLE 2. Regression analysis data for elastic strain, viscous strain and Nss

 
Data on the compressive modulus, elastic strain, viscous strain and creep rate are summarized in Table 1Go and plotted in Fig. 7Go for each of the 23 specimens.



View larger version (17K):
[in this window]
[in a new window]
 
FIG. 7. Results and best-fit curves for all 23 specimens. (a) Variation of elastic strain with cartilage compressive modulus can be expressed by an exponential relationship. (b) The viscous strain of cartilage has a linear correlation with its compressive modulus. (c) The number of load cycles at which steady-state behaviour is approached also varies linearly with the cartilage compressive modulus. Nss corresponds to the number of cycles taken to reach a particular gradient d{varepsilon}/dx=-0.003, {varepsilon} being the total strain and x the number of load cycles.

 

Total strain.
The total strain (sum of the elastic and viscous contributions) incurred by each specimen at the final load cycle is plotted in Fig. 8Go. The best-fit curve to the data was a second-order polynomial with a minimum at a compressive modulus value of approximately 10 MN/m2.



View larger version (17K):
[in this window]
[in a new window]
 
FIG. 8. Total strain data plotted against specimen 50 ms compressive modulus, showing a second-order polynomial best fit and a region of optimum stiffness.

 


    Discussion
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 Conclusions
 References
 
Few experimental studies have measured the deformation and fluid flow of articular cartilage during prolonged cyclic loading [1419]. Hence, the response of cartilage to cyclic loading has been subject to conjecture and prediction based on experimental data from the application of a single load in uniaxial compression/ recovery experiments, and also on various mathematical models.

In his experiments on intact canine joints, Linn [14] showed that static loading caused an initial deformation followed by a creep response that took over 24 h to reach equilibrium. In contrast, he showed that when the load oscillated the deformation became constant, arriving at a load-specific value in 5–6 min. Simon [15] quoted few results but concluded that for a loading interval of 1 s the cumulative deformation was less than that resulting after continuous loading over the same period. Johnson et al. [16] and Higginson and Snaith [17] performed very similar experiments, and both concluded that the tissue response was elastic, fluid flow playing no part in determining the material properties or the behaviour of the tissue under short-term loading. Lee et al. [18] performed sinusoidal loading from 0.001 to 20 Hz and concluded that fluid flow was significant across the whole range of test frequencies, and showed that, at 1 Hz in confined compression, cartilage did not behave as a linear elastic solid. Torzilli [19] was in agreement with Lee et al. [18], and proposed that the deformation behaviour of articular cartilage due to an oscillatory load would be governed by the oscillatory movement of the interstitial fluid during each load cycle. Few results were quoted, but the observations made were similar to those of Linn [14] and Simon [29], namely that the deformation of the tissue reached a cyclic steady state faster than under a static load. In all cases, the accuracy of the measurement techniques was limited and few results were quoted. Varied loading cycles were employed, none of which closely resembled those experienced physiologically.

Maroudas [23] performed no experimental work in this area, but she postulated that recovery between load cycles was likely to be incomplete because the recovery rate would be governed by the swelling pressure. This pressure is always much lower than the physiological stress arising due to the loads acting.

This study measured for the first time the strain response, to prolonged cyclic loading, of knee cartilage specimens possessing a wide range of compressive modulus values. The data obtained show that the technique used is sensitive enough to measure the small residual deformations occurring over individual cycles. Two aspects of the results call for comment and explanation. The first of these is the marked difference in behaviour of stiff and soft cartilage. The data characterize the history of strain components (viscous and elastic) obtained through the entire test period and which revealed a steady state achieved after periods of loading that varied with the stiffness of the cartilage specimen. The second is the variation in the total strain of cartilage over the range of its compressive modulus. The total strain has a minimum value that corresponds to the mid-range value of the compressive modulus. This is different from what would be predicted intuitively; as the stiffness of a material increases its deformation under load is expected to decrease monotonically.

Strain response of soft and stiff cartilage
Considering the rates of fluid movement during the loading and recovery phases of the loading cycle can lead to an explanation of the observations made in the present experiments. During the loading phase, interstitial fluid moves away from the loaded region towards the unloaded region under the action of the pressure gradient caused by the difference between the applied stress and the cartilage swelling pressure. During the recovery phase, the direction of interstitial fluid flow is reversed, but the flow is driven by a much lower pressure gradient as the swelling pressure during recovery is much lower than the applied stress during the loading phase. Consequently, the recovery of cartilage between successive cycles is not complete, and a residual strain is observed. In the present experiment it was observed that, for a stiffer cartilage specimen, larger residual strains occurred (between individual cycles) than those observed for a softer one. Maroudas found that the stiffness and swelling pressure of cartilage are both directly related to its proteoglycan content [30]. Hence, after cartilage is subjected to compression, the swelling pressure of stiff cartilage during the recovery phase of the loading cycle must be higher than that of softer cartilage. Were both tissues (i.e. the stiff and soft cartilage) to have similar permeability, it would be expected that stiff cartilage should recover faster between successive load cycles than does soft cartilage. However, it was shown in the present study that the converse is true. Furthermore, permeability is directly related to cartilage stiffness, stiff cartilage having lower permeability [3133]. Therefore, it is more likely that matrix permeability, which also controls interstitial fluid flow, will be dominant in controlling tissue recovery, and may thereby account for the much slower recovery of the stiffer cartilage.

The above explanation is based on the assumption that the solid matrix of cartilage is elastic. However, in reality it may be viscoelastic in part, and the effect of this on the total strain of cartilage is difficult to measure [6]. If this contribution were large it might instead be the dominant factor in the observed behaviour of cartilage. To clarify this ambiguity, an experiment was undertaken in which cartilage was subjected to the same cyclic loading regime, but a small tare load that was 1% of the maximum load was maintained throughout the test. Thus, the surface of the cartilage was not exposed to the fluid during the period of recovery but was in contact throughout the test with the flat-ended indenter. This indenter was of a weight that supplied the tare load, and transmitted the cyclic load to the cartilage surface via its free end. This continuous contact between the indenter and the cartilage surface drastically slowed the recovery of the deformation between cycles and so increased the time to attain the steady state by a factor of at least three [34]. We postulate that the contact between the indenter and cartilage, even under such a small tare load, has almost blocked the path of the fluid being imbibed by the cartilage matrix during the recovery phase of the loading cycle. Furthermore, the fluid path is further restricted by the decreased permeability of the surface layer caused by its significant compression, even under the action of such a small tare load [19]. Were the inherent viscoelasticity of the matrix a major contributor to the observed viscous behaviour, the recovery time would not have been greatly influenced by the presence of the tare load. It appears, therefore, that the interstitial fluid flow, which is governed by the permeability of cartilage, is the dominant factor in its viscous response to load.

The steady-state results similarly reflect the permeability effect. Steady state is attained when sufficient interstitial fluid has been lost from the loaded regions of cartilage, such that the stiffness of the resulting compacted cartilage is sufficient to carry the applied dynamic stresses. Tests showed that all cartilage specimens lost interstitial fluid to attain this steady state, but the softer cartilage reached it sooner than the stiffer. As would be expected of a more permeable tissue, the softer cartilage lost its excess interstitial fluid more readily than did the stiffer, less permeable cartilage.

Optimum stiffness: does cartilage adapt to the loading regime?
The total strain data gave rise to the observation of an optimum range of cartilage stiffness (8–12 MN/m2) within which cartilage incurred minimum strain—about 23%. Specimens with stiffness outside this range experienced higher total strains. Stiffer cartilage outside the range underwent higher total strains that were contributed to by large viscous strains, whereas softer cartilage outside the range underwent higher total strains due to the contribution of the large elastic strains. On the basis of this observation, it may be hypothesized that cartilage adapts its matrix constituents to be least susceptible to damage, as it minimizes the total matrix strain by optimally balancing the viscous and elastic strain contributions. Intuitively, softer cartilage would be prone to early failure because of large elastic strains, while, perhaps more surprisingly, stiffer cartilage may also be prone to failure after many load cycles due to excessive strains caused by viscous losses.

The results may thus have implications for chondrocyte biosynthesis. In softer cartilage the chondrocyte is subject to large elastic matrix strains coupled with little fluid flow, whereas chondrocytes in stiffer cartilage are subject to smaller elastic strains of the surrounding matrix coupled with greater local fluid flow. These differing combinations of matrix strain and fluid flow could be important mechano-transduction factors affecting the the local structure of tissue synthesized by the chondrocyte.


    Conclusions
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 Conclusions
 References
 
The deformation of cartilage under compressive cyclic loading conditions has a complex relationship with its compressive modulus. The relationship is not an inverse but a bimodal one, and the total strain of cartilage (sum of the elastic and viscous components) had a minimum at the mid-range of the modulus determined at eight predetermined sites on the surfaces of the knee joint. The elastic strain component had an inverse exponential relationship with the modulus, whereas the viscous component increased linearly with increase in the modulus.

The viscoelastic behaviour of cartilage observed under the cyclic loading regime in this study was attributed to the interstitial fluid flow within the cartilage matrix and was explained primarily in terms of the permeability of cartilage rather than the intrinsic viscoelasticity of the matrix itself. This explanation is supported by data from a further experiment in which the total strain of cartilage was measured under the same cyclic loading conditions but in the presence of a small tare load on the cartilage. This delayed the recovery of cartilage between cycles by a factor of at least three, which was attributed to the presence of the tare load, which was transmitted via the flat-ended indenter and which presented the increased resistance in the path of the fluid being imbibed by the cartilage matrix during the recovery period between consecutive load cycles.


    Acknowledgments
 
The authors would like to thank to their technicians Brian Whitham and Michael Pullan for their assistance during the course of this work. This work was supported by a scholarship awarded by the University of Leeds.


    Notes
 
Correspondence to: B. B. Seedhom. Back


    References
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 Conclusions
 References
 

  1. Hayes WC, Mockros, LF. Viscoelastic properties of human articular cartilage. J Appl Physiol1971;31:562–8.[Free Full Text]
  2. Coletti JM, Akeson WH, Woo SL. A comparison of the physical behaviour of normal articular cartilage and the arthroplasty surface. J Bone Joint Surg1972; 54:147–60.[Abstract/Free Full Text]
  3. Elmore SM, Sokoloff L, Norris G, Carmeci P. Nature of ‘imperfect’ elasticity of articular cartilage. J Appl Physiol1963;18:393–6.[Abstract/Free Full Text]
  4. Kempson GE, Freeman MA, Swanson SA. The determination of a creep modulus for articular cartilage from indentation tests of the human femoral head. J Biomech1971a;4:239–50.[Medline]
  5. Kempson GE, Spivey CJ, Swanson SA, Freeman, MA. Patterns of cartilage stiffness on normal and degenerate human femoral heads. J Biomech1971b;4:597–609.[Web of Science][Medline]
  6. Hayes WC, Bodine AJ. Flow-independent viscoelastic properties of articular cartilage matrix. J Biomech1978;11:407–19.[Web of Science][Medline]
  7. Mow VC, Kuei SC, Lai WM, Armstrong CG. Biphasic creep and stress relaxation of articular cartilage in compression. Theory and experiments. J Biomech Eng1980;102:73–84.[Web of Science][Medline]
  8. Mak AF, Lai WM, Mow VC. Biphasic indentation of cartilage: i. Theoretical analysis. J Biomech1987; 20:703–14.[Web of Science][Medline]
  9. Setton LA, Zhu W, Mow VC. The biphasic porovisco-elastic behaviour of articular cartilage: role of the surface zone in governing the compressive behaviour. J Biomech1993;26:581–92.[Web of Science][Medline]
  10. Radin EL, Whittle MW, Yang KH et al. The heelstrike transient, its relationship with the angular velocity of the shank and the effects of quadriceps paralysis. Adv Bioeng1986;BED:121–3.
  11. Mann RA, Hagy J. Biomechanics of walking, running and sprinting. Am J Sports Med1980;8:345–50.[Abstract/Free Full Text]
  12. Nilsson J, Thorstensson A. Adaptability in frequency and amplitude of leg movements during human locomotion at different speeds. Acta Physiol Scand1987; 129:107–14.[Web of Science][Medline]
  13. Barker MK, Seedhom BB. Articular cartilage deformation under physiological cyclic loading—apparatus and measurement technique. J Biomech1997;30:377–81.[Web of Science][Medline]
  14. Linn FC. Lubrication of animal joints. J Bone Joint Surg1967;49A:1079–98.
  15. Simon WH. Scale effects in animal joints. II. Thickness and elasticity in the deformability of articular cartilage. Arthritis Rheum1971;14:493–502.[Medline]
  16. Johnson GR, Dowson D, Wright V. The elastic behaviour of articular cartilage under a sinusoidally varying compressive load. Int J Mech Sci1977;19:301–8.
  17. Higginson GR, Snaith JE. The mechanical stiffness of articular cartilage in confined oscillating compression. Eng Med1979;8:11–4.
  18. Lee RC, Frank EH, Grodzinsky AJ, Roylance DK. Oscillatory compressional behavior of articular cartilage and its associated electromechanical properties. J Biomech Eng1981;103:280–92.[Web of Science][Medline]
  19. Torzilli PA. Mechanical response of articular cartilage to an oscillating load. Mech Res Commun1984; 11:75–82.
  20. Simon BR, Coats RS, Woo SLY. Relaxation and creep quasilinear visco-elastic models for normal articular cartilage. J Biomech Eng1984;106:159–65.[Medline]
  21. Guilak F, Spilker RL, Mow VC. A finite element model of cartilage extracellular matrix response to static and cyclic compressive loading. Adv Biomech1990, BED-IM:225–8.
  22. Suh J, Li Z, Wood S. Dynamic behaviour of a biphasic cartilage model under cyclic compressive loading. J Biomech1995;28:357–64.[Medline]
  23. Maroudas A. Physico-chemical properties of articular cartilage. In: Freeman MAR, ed. Adult articular cartilage. London: Pitman Medical, 1973:163.
  24. Meachim G. Light microscopy of Indian ink preparations of fibrillated cartilage. Ann Rheum Dis1972;31:457–64.[Free Full Text]
  25. Swann AC, Seedhom BB. The stiffness of normal articular cartilage and the predominant acting stress levels: implications for the aetiology of osteoarthrosis. Br J Rheumatol1993;32:16–25.[Abstract/Free Full Text]
  26. Shepherd DET, Seedhom BB. A technique for measuring the compressive modulus of articular cartilage under physiological loading rates, with preliminary results. Eng Med1997;211:155–65.
  27. Shepherd DET, Seedhom BB. The ‘instantaneous’ compressive modulus of human articular cartilage in joints of the lower limb. Rheumatology1999;38:124–32.[Free Full Text]
  28. Yao JQ, Seedhom BB. Mechanical conditioning of articular cartilage to prevalent stresses. Br J Rheumatol1993;32:956–65.[Abstract/Free Full Text]
  29. Simon WH. Scale effects in animal joints. I. Articular cartilage thickness and compressive stress. Arthritis Rheum1970;13:244–55.[Medline]
  30. Maroudas A. A balance between swelling pressure and collagen tension in normal and degenerate cartilage. Nature1976;260:808–9.[Medline]
  31. Jurvelin J, Kiviranta I, Arokoski J, Tammi M, Helminen HJ. Indentation study of the biochemical properties of articular cartilage in the canine knee. Eng Med1987, 16:15–22.[Medline]
  32. Armstrong CG, Mow VC. Variations in the intrinsic mechanical properties of human articular cartilage with age, degeneration and water content. J Bone Joint Surg1982;64A:88–94.
  33. Athanasiou KA, Rosenwasser MP, Buckwalter JA, Malinin TI, Mow VC. Interspecies comparisons on in-situ intrinsic mechanical properties of distal femoral cartilage. J Orthop Res1991;12:340–9.
  34. Barker MK. Aspects of articular cartilage response to physiological cyclic loading with special reference to cumulative deformation [PhD thesis]. Research School of Medicine, University of Leeds, 1997.
Submitted 19 June 2000; Accepted 18 September 2000


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?



This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (11)
Right arrowRequest Permissions
Right arrow Disclaimer
Google Scholar
Right arrow Articles by Barker, M. K.
Right arrow Articles by Seedhom, B. B.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Barker, M. K.
Right arrow Articles by Seedhom, B. B.
Related Collections
Right arrow Osteoarthritis and Cartilage
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?