The low permeability of healthy meniscus and labrum limit articular cartilage consolidation and maintain fluid load support in the knee and hip

The knee meniscus and hip labrum appear to be important for joint health, but the mechanisms by which these structures perform their functions are not fully understood. The fluid phase of articular cartilage provides compressive stiffness and aids in maintaining a low friction articulation. Healthy fibrocartilage, the tissue of meniscus and labrum, has a lower fluid permeability than articular cartilage. In this study we hypothesized that an important function of the knee meniscus and the hip labrum is to augment fluid retention in the articular cartilage of a mechanically loaded joint. Axisymmetric hyperporoelastic finite element models were analyzed for an idealized knee and an idealized hip. The results indicate that the meniscus maintained fluid pressure and inhibited fluid exudation in knee articular cartilage. Similar, but smaller, effects were seen with the labrum in the hip. Increasing the fibrocartilage permeability relative to that of articular cartilage gave a consolidation rate and loss of fluid load support comparable to that predicted by meniscectomy or labrectomy. The reduced articular cartilage fluid pressure that was calculated for the joint periphery is consistent with patterns of endochondral ossification and osteophyte formation in knee and hip osteoarthritis. High articular central strains and loss of fluid load support after meniscectomy could lead to fibrillation. An intact low-permeability fibrocartilage is important for limiting fluid exudation from articular cartilage in the hip and knee. This may be an important aspect of the role of fibrocartilage in protecting these joints from osteoarthritis.     

Partial meniscectomy changes fluid pressurization in articular cartilage in human knees

Partial meniscectomy is believed to change the biomechanics of the knee joint through alterations in the contact of articular cartilages and menisci. Although fluid pressure plays an important role in the load support mechanism of the knee, the fluid pressurization in the cartilages and menisci has been ignored in the finite element studies of the mechanics of meniscectomy. In the present study, a 3D fibril-reinforced poromechanical model of the knee joint was used to explore the fluid flow dependent changes in articular cartilage following partial medial and lateral meniscectomies. Six partial longitudinal meniscectomies were considered under relaxation, simple creep, and combined creep loading conditions. In comparison to the intact knee, partial meniscectomy not only caused a substantial increase in the maximum fluid pressure but also shifted the location of this pressure in the femoral cartilage. Furthermore, these changes were positively correlated to the size of meniscal resection. While in the intact joint, the location of the maximum fluid pressure was dependent on the loading conditions, in the meniscectomized joint the location was predominantly determined by the site of meniscal resection. The partial meniscectomy also reduced the rate of the pressure dissipation, resulting in even larger difference between creep and relaxation times as compared to the case of the intact knee. The knee joint became stiffer after meniscectomy because of higher fluid pressure at knee compression followed by slower pressure dissipation. The present study indicated the role of fluid pressurization in the altered mechanics of meniscectomized knees.     

3D finite element model of meniscectomy: changes in joint contact behavior

The goal of this study is to quantify changes in knee joint contact behavior following varying degrees of the medial partial meniscectomy. A previously validated 3D finite element model was used to simulate 11 different meniscectomies. The accompanying changes in the contact pressure on the superior surface of the menisci and tibial plateau were quantified as was the axial strain in the menisci and articular cartilage. The percentage of medial meniscus removed was linearly correlated with maximum contact pressure, mean contact pressure, and contact area. The lateral hemi-joint was minimally affected by the simulated medial meniscectomies. The location of maximum strain and location of maximum contact pressure did not change with varying degrees of partial medial meniscectomy. When 60% of the medial meniscus was removed, contact pressures increased 65% on the remaining medial meniscus and 55% on the medial tibial plateau. These data will be helpful for assessing potential complications with the surgical treatment of meniscal tears. Additionally, these data provide insight into the role of mechanical loading in the etiology of post-meniscectomy osteoarthritis.     

Effect of superficial collagen patterns and fibrillation of femoral articular cartilage on knee joint mechanics-a 3D finite element analysis

Collagen fibrils of articular cartilage have specific depth-dependent orientations and the fibrils bend in the cartilage surface to exhibit split-lines. Fibrillation of superficial collagen takes place in osteoarthritis. We aimed to investigate the effect of superficial collagen fibril patterns and collagen fibrillation of cartilage on stresses and strains within a knee joint. A 3D finite element model of a knee joint with cartilage and menisci was constructed based on magnetic resonance imaging. The fibril-reinforced poroviscoelastic material properties with depth-dependent collagen orientations and split-line patterns were included in the model. The effects of joint loading on stresses and strains in cartilage with various split-line patterns and medial collagen fibrillation were simulated under axial impact loading of 1000 N. In the model, the collagen fibrils resisted strains along the split-line directions. This increased also stresses along the split-lines. On the contrary, contact and pore pressures were not affected by split-line patterns. Simulated medial osteoarthritis increased tissue strains in both medial and lateral femoral condyles, and contact and pore pressures in the lateral femoral condyle. This study highlights the importance of the collagen fibril organization, especially that indicated by split-line patterns, for the weight-bearing properties of articular cartilage. Osteoarthritic changes of cartilage in the medial femoral condyle created a possible failure point in the lateral femoral condyle. This study provides further evidence on the importance of the collagen fibril organization for the optimal function of articular cartilage.     

Pathways of load-induced cartilage damage causing cartilage degeneration in the knee after meniscectomy

Results of both clinical and animal studies show that meniscectomy often leads to osteoarthritic degenerative changes in articular cartilage. It is generally assumed that this process of cartilage degeneration is due to changes in mechanical loading after meniscectomy. It is, however, not known why and where this cartilage degeneration starts. Load induced cartilage damage is characterized as either type (1)--damage without disruption of the underlying bone or calcified cartilage layer--or type (2), subchondral fracture with or without damage to the overlying cartilage. We asked the question whether cartilage degeneration after meniscectomy is likely to be initiated by type (1) and/or type (2) cartilage damage. To investigate that we applied an axisymmetric biphasic finite element analysis model of the knee joint. In this model the articular cartilage layers of the tibial and the femoral condyles, the meniscus and the bone underlying the articular cartilage of the tibia plateau were included. The model was validated with data from clinical studies, in which the effects of meniscectomy on contact areas and pressures were measured. It was found that both the maximal values and the distributions of the shear stress in the articular cartilage changed after meniscectomy, and that these changes could lead to both type (1) and type (2) cartilage damage. Hence it likely that the cartilage degeneration seen after meniscectomy is initiated by both type (1) and type (2) cartilage damage.     

In situ measurement of articular cartilage deformation in intact femoropatellar joints under static loading

The deformational behavior of articular cartilage has been investigated in confined and unconfined compression experiments and indentation tests, but to date there exist no reliable data on the in situ deformation of the cartilage during static loading. The objective of the current study was to perform a systematic study into cartilage compression of intact human femoro-patellar joints under short- and long-term static loading with MR imaging. A non-metallic pneumatic pressure device was used to apply loads of 150% body weight to six joints within the extremity coil of an MRI scanner. The cartilage was delineated during the compression experiment with previously validated 2D and 3D fat-suppressed gradient echo sequences. We observed a mean (maximal) in situ deformation of 44% (57%) in patellar cartilage after 32 h of loading (mean contact pressure 3.6 MPa), the femoral cartilage showing a smaller amount of deformation than the patella. However, only around 7% of the final deformation (3% absolute deformation) occurred during the first minute of loading. A 43% fluid loss from the interstitial patellar matrix was recorded, the initial fluid flux being 0.217 +/- 0.083 microm/s, and a high inter-individual variability of the deformational behavior (coefficients of variation 11-38%). In conjunction with finite-element analyses, these data may be used to compute the load partitioning between the solid matrix and fluid phase, and to elucidate the etiologic factors relevant in mechanically induced osteoarthritis. They can also provide direct estimates of the mechanical strain to be encountered by cartilage transplants.     

The effect of meniscal tears and resultant partial meniscectomies on the knee contact stresses: a finite element analysis

Knee osteoarthritis (OA) is believed to result from high levels of contact stresses on the articular cartilage and meniscus after meniscal damage. This study investigated the effect of meniscal tears and partial meniscectomies on the peak compressive and shear stresses in the human knee joint. An elaborate three-dimensional finite element model of knee joint including bones, articular cartilages, menisci and main ligaments was developed from computed tomography and magnetic resonance imaging images. This model was used to model four types of meniscal tears and their resultant partial meniscectomies and analysed under an axial 1150 N load at 0° flexion. Three different conditions were compared: a healthy knee joint, a knee joint with medial meniscal tears and a knee joint following partial meniscectomies. The numerical results showed that each meniscal tear and its resultant partial meniscectomy led to an increase in the peak compressive and shear stresses on the articular cartilages and meniscus in the medial knee compartment, especially for partial meniscectomy. Among the four types of meniscal tears, the oblique tear resulted in the highest values of the peak compressive and shear stresses. For the four partial meniscectomies, longitudinal meniscectomy led to the largest increase in these two stresses. The lateral compartment was minimally affected by all the simulations. The results of this study demonstrate meniscal tear and its resultant partial meniscectomy has a positive impact on the maintenance of high levels of contact stresses, which may improve the progression of knee OA, especially for partial meniscectomy. Surgeons should adopt a prudent strategy to preserve the greatest amount of meniscus possible.     

Proteoglycan 4 downregulation in a sheep meniscectomy model of early osteoarthritis

Osteoarthritis is a disease of multifactorial aetiology characterised by progressive breakdown of articular cartilage. In the early stages of the disease, changes become apparent in the superficial zone of articular cartilage, including fibrillation and fissuring. Normally, a monolayer of lubricating molecules is adsorbed on the surface of cartilage and contributes to the minimal friction and wear properties of synovial joints. Proteoglycan 4 is the lubricating glycoprotein believed to be primarily responsible for this boundary lubrication. Here we have used an established ovine meniscectomy model of osteoarthritis, in which typical degenerative changes are observed in the operated knee joints at three months after surgery, to evaluate alterations in proteoglycan 4 expression and localisation in the early phases of the disease. In normal control joints, proteoglycan 4 was immunolocalised in the superficial zone of cartilage, particularly in those regions of the knee joint covered by a meniscus. After the onset of early osteoarthritis, we demonstrated a loss of cellular proteoglycan 4 immunostaining in degenerative articular cartilage, accompanied by a significant (p < 0.01) decrease in corresponding mRNA levels. Early loss of proteoglycan 4 from the cartilage surface in association with a decrease in its expression by superficial-zone chondrocytes might have a role in the pathogenesis of osteoarthritis.     

An augmented lagrangian method for sliding contact of soft tissue

Despite the importance of sliding contact in diarthrodial joints, only a limited number of studies have addressed this type of problem, with the result that the mechanical behavior of articular cartilage in daily life remains poorly understood. In this paper, a finite element formulation is developed for the sliding contact of biphasic soft tissues. The augmented Lagrangian method is used to enforce the continuity of contact traction and fluid pressure across the contact interface. The resulting method is implemented in the commercial software COMSOL Multiphysics. The accuracy of the new implementation is verified using an example problem of sliding contact between a rigid, impermeable indenter and a cartilage layer for which analytical solutions have been obtained. The new implementation's capability to handle a complex loading regime is verified by modeling plowing tests of the temporomandibular joint (TMJ) disc.     

Effects of inserting a pressensor film into articular joints on the actual contact mechanics.

Fuji film has been widely used in studies aimed at obtaining the contact mechanics of articular joints. Once sealed for practical use in biological joints, Fuji Pressensor film has a total effective thickness of 0.30 mm, which is comparable to the cartilage thickness in the joints of many small animals. The average effective elastic modulus of Fuji film is approximately 100 MPa in compression, which is larger by a factor of 100-300 compared to that of normal articular cartilage. Therefore, inserting a Pressensor film into an articular joint will change the contact mechanics of the joint. The measurement precision of the Pressensor film has been determined systematically; however, the changes in contact mechanics associated with inserting the film into joints have not been investigated. This study was aimed at quantifying the changes in the contact mechanics associated with inserting sealed Fuji Pressensor film into joints. Spherical and cylindrical articular joint contact mechanics with and without Pressensor film and for varying degrees of surface congruency were analyzed and compared by using finite element models. The Pressensor film was taken as linearly elastic and the cartilage was assumed to be biphasic, composed of a linear elastic solid phase and an inviscid fluid phase. The present analyses showed that measurements of the joint contact pressures with Fuji Pressensor film will change the maximum true contact pressures by 10-26 percent depending on the loading, geometry of the joints, and the mechanical properties of cartilage. Considering this effect plus the measurement precision of the film (approximately 10 percent), the measured joint contact pressures in a joint may contain errors as large as 14-28 percent.     

Application of the u-p finite element method to the study of articular cartilage.

The finite element method using the principle of virtual work was applied to the biphasic theory to establish a numerical routine for analyses of articular cartilage behavior. The matrix equations that resulted contained displacements of the solid matrix (mu) and true fluid pressure (p) as the unknown variables at the element nodes. Both small and large strain conditions were considered. The algorithms and computer code for the analysis of two-dimensional plane strain, plane stress, and axially symmetric cases were developed. The u-p finite element numerical procedure demonstrated excellent agreement with available closed-form and numerical solutions for the configurations of confined compression and unconfined compression under small strains, and for confined compression under large strains. The model was also used to examine the behavior of a repaired articular surface. The differences in material properties between the repair tissue and normal cartilage resulted in significant deformation gradients across the repair interface as well as increased fluid efflux from the tissue.     

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.     

Biphasic finite element modeling of hydrated soft tissue contact using an augmented Lagrangian method

A study of biphasic soft tissues contact is fundamental to understanding the biomechanical behavior of human diarthrodial joints. To date, biphasic-biphasic contact has been developed for idealized geometries and not been accessible for more general geometries. In this paper a finite element formulation is developed for contact of biphasic tissues. The augmented Lagrangian method is used to enforce the continuity of contact traction and fluid pressure across the contact interface, and the resulting method is implemented in the commercial software COMSOL Multiphysics. The accuracy of the implementation is verified using 2D axisymmetric problems, including indentation with a flat-ended indenter, indentation with spherical-ended indenter, and contact of glenohumeral cartilage layers. The biphasic finite element contact formulation and its implementation are shown to be robust and able to handle physiologically relevant problems.     

An in vitro model for the pathological degradation of articular cartilage in osteoarthritis

The objective of this study was to develop an in vitro cartilage degradation model that emulates the damage seen in early-stage osteoarthritis. To this end, cartilage explants were collagenase-treated to induce enzymatic degradation of collagen fibers and proteoglycans at the articular surface. To assess changes in mechanical properties, intact and degraded cartilage explants were subjected to a series of confined compression creep tests. Changes in extracellular matrix structure and composition were determined using biochemical and histological approaches. Our results show that collagenase-induced degradation increased the amount of deformation experienced by the cartilage explants under compression. An increase in apparent permeability as well as a decrease in instantaneous and aggregate moduli was measured following collagenase treatment. Histological analysis of degraded explants revealed the presence of surface fibrillation, proteoglycan depletion in the superficial and intermediate zones and loss of the lamina splendens. Collagen cleavage was confirmed by the Col II–3/4C_short antibody. Degraded specimens experienced a significant decrease in proteoglycan content but maintained total collagen content. Repetitive testing of degraded samples resulted in the gradual collapse of the articular surface and the compaction of the superficial zone. Taken together, our data demonstrates that enzymatic degradation with collagenase can be used to emulate changes seen in early-stage osteoarthritis. Further, our in vitro model provides information on cartilage mechanics and insights on how matrix changes can affect cartilage's functional properties. More importantly, our model can be applied to develop and test treatment options for tissue repair.     

High-resolution measurements of the multilayer ultra-structure of articular cartilage and their translational potential

Current musculoskeletal imaging techniques usually target the macro-morphology of articular cartilage or use histological analysis. These techniques are able to reveal advanced osteoarthritic changes in articular cartilage but fail to give detailed information to distinguish early osteoarthritis from healthy cartilage, and this necessitates high-resolution imaging techniques measuring cells and the extracellular matrix within the multilayer structure of articular cartilage. This review provides a comprehensive exploration of the cellular components and extracellular matrix of articular cartilage as well as high-resolution imaging techniques, including magnetic resonance image, electron microscopy, confocal laser scanning microscopy, second harmonic generation microscopy, and laser scanning confocal arthroscopy, in the measurement of multilayer ultra-structures of articular cartilage. This review also provides an overview for micro-structural analysis of the main components of normal or osteoarthritic cartilage and discusses the potential and challenges associated with developing non-invasive high-resolution imaging techniques for both research and clinical diagnosis of early to late osteoarthritis.     

A finite element exploration of cartilage stress near an articular incongruity during unstable motion

Both instability and residual articular incongruity are implicated in the development of post-traumatic osteoarthritis (OA) following intra-articular fracture, but currently no information exists regarding cartilage stresses for unstable residual incongruities. In this study, a transversely isotropic poroelastic cartilage finite element model was implemented and validated within physiologically relevant loading ranges. This material model was then used to simulate the loading of cartilage during stable and unstable motion accompanying a step-off incongruity residual from intra-articular fracture, using load data from previous cadaver tests of ankle instability. Peak solid-phase stresses and fluid pressure were found to increase markedly in the presence of instability. Solid-phase transients of normal stress increased from 2.00 to 13.8 MPa/s for stable compared to unstable motion, and tangential stress transients increased from 17.1 to 118.1 MPa/s. Corresponding fluid pressure transients increased from 15.1 to 117.9 MPa/s for unstable motion. In the most rapidly loaded sections of cartilage, the fluid was found to carry nearly all of the normal load, with the pressurization of the fluid resulting in high solid matrix tangential stresses.     

A finite element formulation and program to study transient swelling and load-carriage in healthy and degenerate articular cartilage

The theory of poroelasticity is extended to include physico-chemical swelling and used to predict the transient responses of normal and degenerate articular cartilage to both chemical and mechanical loading; with emphasis on isolating the influence of the major parameters which govern its deformation. Using a new hybrid element, our mathematical relationships were implemented in a purpose-built poroelastic finite element analysis algorithm (u-pi-c program) which was used to resolve the nature of the coupling between the mechanical and chemical responses of cartilage when subjected to ionic transport across its membranous skeleton. Our results demonstrate that one of the roles of the strain-dependent matrix permeability is to limit the rate of transmission of stresses from the fluid to the collagen-proteoglycan solid skeleton in the incipient stages of loading, and that the major contribution of the swelling pressure is that of preventing any excessive deformation of the matrix.     

The mechanics of focal chondral defects in the hip

There is a mean incidence of osteoarthritis (OA) of the hip in 8% of the overall population. In the presence of focal chondral defects, defined as localized damage to the articular cartilage, there is an increased risk of symptomatic progression toward OA. This relationship between chondral defects and subsequent development of OA has led to substantial efforts to develop effective procedures for surgical cartilage repair. This study examined the effects of chondral defects and labral delamination on cartilage mechanics in the dysplastic hip during the gait cycle using subject-specific finite element analysis. Models were generated from volumetric CT data and analyzed with simulated chondral defects at the chondrolabral junction on the posterior acetabulum during five distinct points in the gait cycle. Focal chondral defects increased maximum shear stress on the osteochondral surface of the acetabular cartilage, when compared to the intact case. This effect was amplified with labral delamination. Additionally, chondral defects increased the first principal Lagrange strain on the articular surface of the acetabular cartilage and labrum. Labral delamination relieved some of this tensile strain. As defect size was increased, contact stress increased in the medial zone of the acetabulum, while it decreased anteriorly. The results suggest that in the presence of chondral defects and labral delamination the cartilage experiences elevated tensile strains and shear and contact stress, which could lead to further damage of the cartilage, and subsequent arthritic progression. The framework presented here will serve as the procedure for future finite element studies on cartilage mechanics in hips with varying disease states with simulated chondral defects and labral tears.     

The influence and biomechanical role of cartilage split line pattern on tibiofemoral cartilage stress distribution during the stance phase of gait

This study presents an evaluation of the role that cartilage fibre ‘split line’ orientation plays in informing femoral cartilage stress patterns. A two-stage model is presented consisting of a whole knee joint coupled to a tissue-level cartilage model for computational efficiency. The whole joint model may be easily customised to any MRI or CT geometry using free-form deformation. Three ‘split line’ patterns (medial–lateral, anterior–posterior and random) were implemented in a finite element model with constitutive properties referring to this ‘split line’ orientation as a finite element fibre field. The medial–lateral orientation was similar to anatomy and was derived from imaging studies. Model predictions showed that ‘split lines’ are formed along the line of maximum principal strains and may have a biomechanical role of protecting the cartilage by limiting the cartilage deformation to the area of higher cartilage thickness.     

Fluid load support during localized indentation of cartilage with a spherical probe

Interstitial fluid pressurization, a consequence of a biphasic tissue structure, is essential to the load bearing and lubrication properties of articular cartilage. Focal tissue degradation may interfere with this protective mechanism, eventually leading to gross degeneration and osteoarthritis. Our long-term goal is to determine whether local contacts can be used as a means to probe local tissue integrity and functionality. In the present work, Hertzian rate-controlled microindentation was used as a model of the more complicated sliding system to directly determine the effects of contact radius and deformation rate on interstitial load support. During localized contact between a steel spherical probe and bovine articular cartilage, the equilibrium and non-equilibrium responses were well-fit by the Hertz model (R(2)>0.998) with a mean equilibrium contact modulus of 0.93 MPa. The effective contact modulus and fluid load fraction were independent of indentation depth, contact radius, and normal force; both increased monotonically with indentation rate. At 21 μm/s indentation rate, the cartilage was effectively stiffened by 6-fold with the fluid pressure supporting 85% of the contact force. The results motivated a simple analytical model that directly links the tribomechanical response (including fluid load support) and the Peclet number to measurable material properties and controllable experimental variables. This paper demonstrates that tribological contacts can be used to probe local functional properties. Such measurements can add important insights into the roles of focal tissue damage and impaired local functionality in the pathogenesis of osteoarthritis.