12 Biomechanics of Cartilage



10.1055/b-0035-122012

12 Biomechanics of Cartilage

Lutz Dürselen and Andreas Martin Seitz

Diarthrodial joints are covered with hyaline articular cartilage (AC) enabling the transmission of high loads at low friction. Hence, AC plays a very important mechanical role in our musculoskeletal system. Unfortunately, AC can degenerate idiopathically or due to trauma or joint misalignment leading to osteoarthritis with a potential need for joint replacement. For decades, methods have been investigated to heal or regenerate AC with different strategies. To evaluate such methods, it is important to understand the biomechanical way of functioning of AC. Its complex functionality is accomplished by the biphasic composition of this tissue. It contains chondrocytes, which secrete an extracellular matrix with the main constituents collagen and proteoglycans. Besides this solid matrix, AC possesses a high water content as its fluid component. While collagen is mainly responsible for the mechanical strength of the solid matrix, the proteoglycans are able to store water by osmotic swelling providing the resistance to compression and creating the time-dependent behavior, which results in creep and relaxation phenomena. Any change in composition of the extracellular matrix resulting from trauma or degenerative processes leads to changes in biomechanical properties. Therefore, the biomechanical characterization of AC is an important and useful tool to assess its functional quality and can be used in any kind of research on AC such as basic questions on the functionality of AC, degenerative processes, repair, and tissue engineering aiming to replace destroyed cartilage. A variety of different test methods for the biomechanical characterization of AC are available depending on which properties are of interest regarding a special research question. This chapter gives an overview of practical methods to assess the compressive and time-dependent properties including friction of AC.



12.1 Unconfined Compression and Tensile Testing


The simplest way to determine the material properties of AC is unconfined compression or tensile testing of geometrically standardized specimens. From such experiments, parameters can be derived such as stiffness (K in N/mm), deduced from the force-elongation diagram, and Young′s modulus (E in N/mm2 or MPa), determined from the stress-strain diagram. Although both tensile and unconfined compression testing lead to the same parameters, the results are different, depending on the sample and test conditions. For example, the cartilage tensile modulus has been reported to be one to two orders of magnitude greater than the corresponding compressive modulus. 1 Testing AC under compressive loads commonly involves cylindrical specimens, which can be easily harvested by punches or hollow drills, or cubic specimens, whereas for tensile tests mostly rectangular strips or dumbbell-shaped samples are used, which are more complicated to create. A beneficial attribute associated with the compressive testing method is its comparatively simple test setup (Fig. 12.1a) where a flat-ended rod is used to compress the specimen to a previously defined load, stress, deformation, or strain. Due to the varying Poisson′s ratio of different soft tissues, the cross section of the cartilage specimen changes differently during compressive (ν = 0.2) 2 and tensile (ν > 0.5)1 loads. Compared to unconfined compression, tensile testing is more elaborate, because it implies the problem of clamping soft tissue, which is a general pitfall in soft tissue biomechanics. Therefore, the test setup for tensile testing (Fig. 12.1b) commonly consists of two clamps to ensure adequate fixation of the samples during tensile loading.

Fig. 12.1 Unconfined compression of a cubic cartilage specimen (a) leading to a change in its geometric shape and cross-sectional area, and tensile test of a dumbbell-shaped cartilage specimen (b) leading to thinning of its cross-sectional area under tensile loads.

During testing, it should be considered that both test methods do not intend to represent physiological or in vivo loading conditions. This involves the integrity of the cartilage sample with respect to the loading characteristics, as well as the ambient conditions (e.g., temperature, nutrition, and moistness). Due to its anisotropic and biphasic material behavior, AC shows a typical strain-rate–dependent response during unconfined compression and tensile loading where an increased strain rate is associated with strengthening the transient stiffness and Young′s modulus 3, 4 (Fig. 12.2a). Another phenomenon can be observed when loading the cartilage specimen to a certain point with instant unloading, leading to a so-called hysteresis loop, which is common for biphasic soft tissues like cartilage 5 (Fig. 12.2b).

Fig. 12.2 The Young′s modulus of articular cartilage samples is dramatically influenced by the loading rate during compression or tensile testing (a): With increasing loading rate (ν) the initial resistance, and therefore the Young′s modulus, significantly increases. The same phenomenon can be observed when loading and immediately unloading the specimen, resulting in a so-called hysteresis loop (b). Common for biphasic materials like cartilage, the hysteresis loop indicates dissipated energy (∆W), which is defined as the area embedded between the loading and unloading curve. This area is also loading-rate (ν) sensitive.


Jargon Simplified:




  • Load (F in N): Usually results in a deformation of the tested material.



  • Stress (σ in N/mm2 or MPa): Measure of the internal force acting on a material. Defined as the measure of the force that is applied via a certain area.



  • Deformation (∆l in mm): Change in the dimensions of a material to external forces or stress.



  • Strain (ε in %): Deformation of a material related to its initial length or height.



  • Poisson′s Ratio (ν): Tapering or bulging ratio of a material in the direction perpendicular to the applied tensile or compressive stress, relative to its initial state.



12.2 Creep/Relaxation and Dynamic Compression



Jargon Simplified




  • Aggregate Modulus (HA in N/mm2 or MPa): Modulus of a biphasic material in its equilibrium state reflecting the properties of the solid matrix. Can be calculated by dividing the equilibrium stress (σt→∝) by the corresponding strain (ε):



  • Hydraulic Permeability (k in m4/Ns): Parameter to describe the fluid flow through a porous material.


When repeating the previously introduced hysteresis loop several times, their characteristics will change in a way that the differences of the loading and unloading curve start to converge (Fig. 12.3a). Additionally, when using a time-resolving illustration of the stress and strain curves of these loading and unloading ramps, a characteristic phase shift of the material response could be observed (Fig. 12.3b). This phase shift also happens during sinusoidal or other alternating loading schemes.

Fig. 12.3 (a) Alternating hysteresis loop: The initial hysteresis loop indicates a bigger difference between the loading and unloading curve leading to a more distinct amount of dissipated energy, compared to the fifth hysteresis loop. (b) Typical phase shift (α) during cyclic loading of cartilage: For a biphasic material, like cartilage, the strain follows the stress. In contrast, a purely elastic material will not have any phase shift, whereas a purely viscoelastic material will have a phase shift of 90 degrees, which means that the highest deformation would be observed at the lowest applied force level and vice versa.

In addition, the characteristic time-dependent properties of AC could be assessed using two more loading tests:




  • Creep: Loading of a cartilage sample with a constant load. The samples thickness diminishes by fluid efllux until an asymptotic thickness is reached.



  • Stress-relaxation: Compression of a cartilage sample until a desired strain is reached. The applied strain is kept constant over time, whereas the applied load decreases until an asymptotic equilibrium force is noted.


Both load methods are able to reveal the viscoelastic behavior of cartilage over time. However, to gain the aggregate modulus (HA in MPa or kPa) and hydraulic permeability (k in m4/Ns), each test routine could be conducted in two different ways:




  • Indentation test (Fig. 12.4a). Here, a flat indenter is pushed into the cartilaginous surface. The respective decrease of thickness (creep) or load (relaxation) is then measured over time. There is no need for preparing complex, geometrically standardized cartilage samples, because the whole specimen could be mounted into the testing rig, while ensuring adequate positioning. However, the cartilage thickness must be determined in another test (e.g., needle test), because without knowing the sample thickness the application of a certain strain level and relative loading rate is impossible.



  • Confined compression testing (Fig. 12.4b). A coplanar, geometrically defined cartilage disk is compressed to a certain load or strain level and held constant until an equilibrium state is reached. In contrast, the test setup is more complex as it must be ensured to constrain the samples in their radial direction while allowing an adequate fluid flow in the axial direction.

Fig. 12.4 (a) Setup for indentation testing: An indenter is adapted to a load cell applying a defined load perpendicular to the cartilage surface, which is inserted in an adjustable measuring chamber, filled with saline solution. The indentation depth is measured using an accurate displacement transducer. (b) Schematic drawing of a test setup for a compression test under confined conditions: The geometrically defined sample is radially constrained allowing fluid flow only in the axial direction. The load is applied via a punch and a porous cylinder, which is connected to a load cell. The strain is accurately measured using a high-precision laser displacement transducer.

When determining the viscoelastic material parameters, various constitutive mathematical models could be used to solve their according equations. 6, 7, 8 However, it should be kept in mind that each test method and each mathematical model can lead to different results for the time-dependent, viscoelastic parameters, because mathematical models mostly are simplified models and do not necessarily match reality 100%. Therefore, care must be taken when comparing results with already existing results in literature.



12.2.1 Example for the Experimental Characterization of the Viscoelastic Properties of Articular Cartilage: Stress-Relaxation under Confined Compression


After harvesting a coplanar cartilaginous disk (thickness: 1 mm; Ø: 5 mm), the sample could be tested using a special measuring chamber providing confined conditions (Fig. 12.4b). This means that the cartilage sample cannot extend in the radial direction and no fluid can flow through the sides of the cartilage disk. The measurement chamber is filled with saline solution to prevent dehydration. The load is applied via a porous Al2O3 cylinder allowing adequate fluid flow through the surface of the sample.


During testing, it is important to measure the deformation, especially the strain, directly at the surface of the cartilage sample as, for example, measurement errors related to deflection of the crosshead of a materials testing machine or of the load cell should be avoided. A dedicated external displacement sensor with high accuracy can be used for this critical issue.


The viscoelastic properties strongly depend on the extent of the applied compression. Therefore, it is recommended to perform a test routine with incremental steps, for example increasing the strain stepwise from 5 to 20% of its initial thickness (Fig. 12.5a). The decrease of load is recorded for about 10 minutes or at least until the equilibrium state of each strain step is reached. As previously mentioned, the strain rate during the compression phase of the test also considerably influences the peak load occurring directly after reaching the desired strain step. Therefore, the strain rate must be kept constant over the whole test series of the experiment. Thus it is recommended to use a strain-dependent loading rate (e.g., 10% *h 0/min [h 0 = initial height of the specimen]) to load each specimen equally.

Fig. 12.5 (a) Force-time-graph (top) and related compression time graph (bottom) of stress-relaxation test with four defined compression states. (b) Screenshot (JMP V.3.2; SAS Institute GmbH, Germany) of the results of the modulus (H) and hydraulic permeability (k), based on first approximation from fitting the previously measured data to a diffusion equation by using least square fit algorithm.

For data analysis, the modulus (H) indicating a combined modulus of both the fluid and the solid phase of the viscoelastic material and the hydraulic permeability (k) can be extracted in a first approximation from the measured relaxation curves by fitting the data to a diffusion equation by using a least square fitting algorithm using the following equation6,7:


where σ(t) is the present stress in MPa, σ(t=∝) the equilibrium stress in MPa, H the above-mentioned combined modulus in MPa, ∆l/d the applied strain level, d the sample thickness in mm, and k the hydraulic permeability in m4/Ns. The determination of the values for H and k can be performed by any statistical software (Fig. 12.5).


The exemplary results of the nonlinear curve fitting, calculated with JMP (V.3.2; SAS Institute GmbH, Germany) indicate an H of 0.79 MPa and k = 3.05 * 10 to 15 m4/Ns. Furthermore, the graph (Fig. 12.5) shows that the adaptation of the model equation to the measured data is good, but not optimal. This means that the real process does not exactly follow the chosen diffusion model, which does not reflect the dynamic behavior of the cartilage sample directly after application of the load. Therefore, an adequate number of measured data points must be excluded from the fitting process until the correlation coefficient is sufficiently high (e.g., 0.95).



Key Concepts: Mathematical Modeling of Articular Cartilage


Due to the anisotropic behavior of cartilage, there are plenty of mathematical descriptions of the response of cartilage to a various number of external loads and environmental conditions. After measuring the in vitro response of articular cartilage, an adequate mathematical model must be chosen to determine the required material properties. Usually, the model parameters are then varied until the model converges with the experimental data. For this purpose, statistical software can be used providing a variety of fitting algorithms (e.g., least-square, Levenberg-Marquardt, Gauss-Newton). After running the fitting algorithm, the mathematical model results in a best possible representation of the experimental data. This could be checked using the Pearson correlation coefficient and setting its accuracy limit, for example, to r = 0.95.

Only gold members can continue reading. Log In or Register to continue

Stay updated, free articles. Join our Telegram channel

Jun 10, 2020 | Posted by in ORTHOPEDIC | Comments Off on 12 Biomechanics of Cartilage

Full access? Get Clinical Tree

Get Clinical Tree app for offline access