Role of the acetabular labrum on articular cartilage consolidation patterns

Biomechanics and Modeling in Mechanobiology - Damage to the acetabular labrum has been associated with cartilage degeneration. Because conventional pressure measurement devices were unable to...     

Role of the acetabular labrum in load support across the hip joint

The relatively high incidence of labral tears among patients presenting with hip pain suggests that the acetabular labrum is often subjected to injurious loading in vivo. However, it is unclear whether the labrum participates in load transfer across the joint during activities of daily living. This study examined the role of the acetabular labrum in load transfer for hips with normal acetabular geometry and acetabular dysplasia using subject-specific finite element analysis. Models were generated from volumetric CT data and analyzed with and without the labrum during activities of daily living. The labrum in the dysplastic model supported 4-11% of the total load transferred across the joint, while the labrum in the normal model supported only 1-2% of the total load. Despite the increased load transferred to the acetabular cartilage in simulations without the labrum, there were minimal differences in cartilage contact stresses. This was because the load supported by the cartilage correlated with the cartilage contact area. A higher percentage of load was transferred to the labrum in the dysplastic model because the femoral head achieved equilibrium near the lateral edge of the acetabulum. The results of this study suggest that the labrum plays a larger role in load transfer and joint stability in hips with acetabular dysplasia than in hips with normal acetabular geometry.     

Articular cartilage friction increases in hip joints after the removal of acetabular labrum

The acetabular labrum is believed to have a sealing function. However, a torn labrum may not effectively prevent joint fluid from escaping a compressed joint, resulting in impaired lubrication. We aimed to understand the role of the acetabular labrum in maintaining a low friction environment in the hip joint. We did this by measuring the resistance to rotation (RTR) of the hip, which reflects the friction of the articular cartilage surface, following focal and complete labrectomy. Five cadaveric hips without evidence of osteoarthritis and impingement were tested. We measured resistance to rotation of the hip joint during 0.5, 1, 2, and 3 times body weight (BW) cyclic loading in the intact hip, and after focal and complete labrectomy. Resistance to rotation, which reflects articular cartilage friction in an intact hip was significantly increased following focal labrectomy at 1-3 BW loading, and following complete labrectomy at all load levels. The acetabular labrum appears to maintain a low friction environment, possibly by sealing the joint from fluid exudation. Even focal labrectomy may result in increased joint friction, a condition that may be detrimental to articular cartilage and lead to osteoarthritis.     

Effects of idealized joint geometry on finite element predictions of cartilage contact stresses in the hip

Computational models may have the ability to quantify the relationship between hip morphology, cartilage mechanics and osteoarthritis. Most models have assumed the hip joint to be a perfect ball and socket joint and have neglected deformation at the bone-cartilage interface. The objective of this study was to analyze finite element (FE) models of hip cartilage mechanics with varying degrees of simplified geometry and a model with a rigid bone material assumption to elucidate the effects on predictions of cartilage stress. A previously validated subject-specific FE model of a cadaveric hip joint was used as the basis for the models. Geometry for the bone-cartilage interface was either: (1) subject-specific (i.e. irregular), (2) spherical, or (3) a rotational conchoid. Cartilage was assigned either a varying (irregular) or constant thickness (smoothed). Loading conditions simulated walking, stair-climbing and descending stairs. FE predictions of contact stress for the simplified models were compared with predictions from the subject-specific model. Both spheres and conchoids provided a good approximation of native hip joint geometry (average fitting error ～0.5 mm). However, models with spherical/conchoid bone geometry and smoothed articulating cartilage surfaces grossly underestimated peak and average contact pressures (50% and 25% lower, respectively) and overestimated contact area when compared to the subject-specific FE model. Models incorporating subject-specific bone geometry with smoothed articulating cartilage also underestimated pressures and predicted evenly distributed patterns of contact. The model with rigid bones predicted much higher pressures than the subject-specific model with deformable bones. The results demonstrate that simplifications to the geometry of the bone-cartilage interface, cartilage surface and bone material properties can have a dramatic effect on the predicted magnitude and distribution of cartilage contact pressures in the hip joint.     

Biomechanical influence of disk properties on the load transfer of healthy and degenerated disks using a poroelastic finite element model

Spine degeneration is a pathology that will affect 80% of the population. Since the intervertebral disks play an important role in transmitting loads through the spine, the aim of this study was to evaluate the biomechanical impact of disk properties on the load carried by healthy (Thompson grade I) and degenerated (Thompson grades III and IV) disks. A three-dimensional parametric poroelastic finite element model of the L4/L5 motion segment was developed. Grade I, grade II, and grade IV disks were modeled by altering the biomechanical properties of both the annulus and nucleus. Models were validated using published creep experiments, in which a constant compressive axial stress of 0.35 MPa was applied for 4 h. Pore pressure (PP) and effective stress (S(E)) were analyzed as a function of time following loading application (1 min, 5 min, 45 min, 125 min, and 245 min) and discal region along the midsagittal profile for each disk grade. A design of experiments was further implemented to analyze the influence of six disk parameters (disk height (H), fiber proportion (%F), drained Young's modulus (E(a),E(n)), and initial permeability (k(a),k(n)) of both the annulus and nucleus) on load-sharing for disk grades I and IV. Simulations of grade I, grade III, and grade IV disks agreed well with the available published experimental data. Disk height (H) had a significant influence (p<0.05) on the PP and S(E) during the entire loading history for both healthy and degenerated disk models. Young's modulus of the annulus (E(a)) significantly affected not only S(E) in the annular region for both disk grades in the initial creep response but also S(E) in the nucleus zone for degenerated disks with further creep response. The nucleus and annulus permeabilities had a significant influence on the PP distribution for both disk grades, but this effect occurred at earlier stages of loading for degenerated than for healthy disk models. This is the first study that investigates the biomechanical influence of both geometrical and material disk properties on the load transfer of healthy and degenerated disks. Disk height is a significant parameter for both healthy and degenerated disks during the entire loading. Changes in the annulus stiffness, as well as in the annulus and nucleus permeability, control load-sharing in different ways for healthy and degenerated disks.     

The contribution of the glenoid labrum to glenohumeral stability under physiological joint loading using finite element analysis

The labrum contributes to passive glenohumeral joint stability. Cadaveric studies have demonstrated that this has position and load dependency, which has not been quantified under physiological loads. This study aims to validate subject-specific finite element (FE) models against in vitro measurements of joint stability and to utilise the FE models to predict joint stability under physiological loads. The predicted stability values were within ± one standard deviation of experimental data and the FE models showed a reduction in stability of 10–15% with high, physiological, loads. The developed regression equations provide the first representation of passive glenohumeral stability and will aid surgical decision-making.     

A nonlinear finite element model of cartilage growth

The long range objective of this work is to develop a cartilage growth finite element model (CGFEM), based on the theories of growing mixtures that has the capability to depict the evolution of the anisotropic and inhomogeneous mechanical properties, residual stresses, and nonhomogeneities that are attained by native adult cartilage. The CGFEM developed here simulates isotropic in vitro growth of cartilage with and without mechanical stimulation. To accomplish this analysis a commercial finite element code (ABAQUS) is combined with an external program (MATLAB) to solve an incremental equilibrium boundary value problem representing one increment of growth. This procedure is repeated for as many increments as needed to simulate the desired growth protocol. A case study is presented utilizing a growth law dependent on the magnitude of the diffusive fluid velocity to simulate an in vitro dynamic confined compression loading protocol run for 2 weeks. The results include changes in tissue size and shape, nonhomogeneities that develop in the tissue, as well as the variation that occurs in the tissue constitutive behavior from growth.     

The influence of the acetabular labrum seal,intact articular superficial zone and synovial fluid thixotropy on squeeze-film lubrication of a spherical synovial joint

A model of synovial fluid (SF) filtration by articular cartilage (AC) in a step-loaded spherical synovial joint at rest is presented. The effects of joint pathology (such as a depleted acetabular labrum, a depleted cartilage superficial zone consistent with early osteoarthritis and an inflammatory SF) on the squeezed synovial film are also investigated. Biphasic mixture models for AC (ideal fluid and elastic porous transversely isotropic two-layer matrix) and for SF (ideal and thixotropic fluids) are applied and the following results are obtained. If the acetabular labrum is able to seal the pressurised SF between the articular surfaces (as in the normal hip joint), the fluid in the synovial film and in the cartilage within the labral ring is homogeneously pressurised. The articular surfaces remain separated by a fluid film for minutes. If the labrum is destroyed or absent and the SF can escape across the contact edge, the fluid pressure is non-homogeneous and with a small jump at the articular surface at the very moment of load application. The ensuing synovial film filtration by porous cartilage is lower for the normal cartilage (with the intact superficial zone) than if this zone is already depleted or rubbed off as in the early stage of primary osteoarthritis. Compared with the inflammatory (Newtonian) SF, the normal (thixotropic) fluid applies favourably in the squeezed film near the contact centre only, yielding a thicker SF film there, but not affecting the minimum thickness in the fluid film profile at a fixed time. For all that, in the unsealed case for both the normal and pathological joint, the macromolecular concentration of the hyaluronic acid-protein complex in the synovial film quickly increases due to the filtration in the greater part of the contact. A stable synovial gel film, thick on the order of 10(-7)m, protecting the articular surfaces from the intimate contact, is formed within a couple of seconds. Boundary lubrication by the synovial gel is established if sliding motion follows until a fresh SF is entrained into the contact. This theoretical prediction is open for experimental verifications.     

A musculoskeletal finite element model of rat knee joint for evaluating cartilage biomechanics during gait

Abnormal loading of the knee due to injuries or obesity is thought to contribute to the development of osteoarthritis (OA). Small animal models have been used for studying OA progression mechanisms. However, numerical models to study cartilage responses under dynamic loading in preclinical animal models have not been developed. Here we present a musculoskeletal finite element model of a rat knee joint to evaluate cartilage biomechanical responses during a gait cycle. The rat knee joint geometries were obtained from a 3-D MRI dataset and the boundary conditions regarding loading in the joint were extracted from a musculoskeletal model of the rat hindlimb. The fibril-reinforced poroelastic (FRPE) properties of the rat cartilage were derived from data of mechanical indentation tests. Our numerical results showed the relevance of simulating anatomical and locomotion characteristics in the rat knee joint for estimating tissue responses such as contact pressures, stresses, strains, and fluid pressures. We found that the contact pressure and maximum principal strain were virtually constant in the medial compartment whereas they showed the highest values at the beginning of the gait cycle in the lateral compartment. Furthermore, we found that the maximum principal stress increased during the stance phase of gait, with the greatest values at midstance. We anticipate that our approach serves as a first step towards investigating the effects of gait abnormalities on the adaptation and degeneration of rat knee joint tissues and could be used to evaluate biomechanically-driven mechanisms of the progression of OA as a consequence of joint injury or obesity.     

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.     

The mechanical environment of the chondrocyte: a biphasic finite element model of cell-matrix interactions in articular cartilage

Mechanical compression of the cartilage extracellular matrix has a significant effect on the metabolic activity of the chondrocytes. However, the relationship between the stress–strain and fluid-flow fields at the macroscopic “tissue” level and those at the microscopic “cellular” level are not fully understood. Based on the existing experimental data on the deformation behavior and biomechanical properties of articular cartilage and chondrocytes, a multi-scale biphasic finite element model was developed of the chondrocyte as a spheroidal inclusion embedded within the extracellular matrix of a cartilage explant. The mechanical environment at the cellular level was found to be time-varying and inhomogeneous, and the large difference (3 orders of magnitude) in the elastic properties of the chondrocyte and those of the extracellular matrix results in stress concentrations at the cell–matrix border and a nearly two-fold increase in strain and dilatation (volume change) at the cellular level, as compared to the macroscopic level. The presence of a narrow “pericellular matrix” with different properties than that of the chondrocyte or extracellular matrix significantly altered the principal stress and strain magnitudes within the chondrocyte, suggesting a functional biomechanical role for the pericellular matrix. These findings suggest that even under simple compressive loading conditions, chondrocytes are subjected to a complex local mechanical environment consisting of tension, compression, shear, and fluid pressure. Knowledge of the local stress and strain fields in the extracellular matrix is an important step in the interpretation of studies of mechanical signal transduction in cartilage explant culture models.     

Determining the location of hip joint centre: application of a conchoid's shape to the acetabular cartilage surface of magnetic resonance images

Preoperative planning, or intraoperative navigation of hip surgery, including joint-preserving procedures such as osteotomy or joint-replacing procedures such as total arthroplasty, needs to be performed with a high degree of accuracy to ensure a successful outcome. The ability to precisely localise the hip joint rotation centre may prove to be very useful in this context. The human hip joint has been shown to be a conchoid shape, and therefore the accurate location of the hip joint centre (HJC) cannot be computed simply as the centre of a sphere. This study describes a method for determining the HJC by applying a conchoid shape to the acetabular cartilage surface of magnetic resonance images, in order to increase the accuracy of the HJC location which had previously been calculated by a functional method using reconstructed three-dimensional surface bony models. By approximating a conchoid shape to the acetabulum, it was possible to compensate for HJC calculation errors.     

The effect of degenerative morphological changes of the intervertebral disc on the lumbar spine biomechanics: a poroelastic finite element investigation

Intervertebral disc degeneration involves changes in the spinal anatomical structures. The mechanical relevance of the following changes was investigated: disc height, endplate sclerosis, disc water content, permeability and depressurisation. A poroelastic nonlinear finite element model of the L4-L5 human spine segments was employed. Loads represented a daily cycle (500 N compression combined with flexion-extension motion for 16 h followed by 200 N compression for 8 h). In non-degenerative conditions, the model predicted a diurnal axial displacement of 1.32 mm and a peak intradiscal pressure of 0.47 MPa. Axial displacement, facet force and range of motion in flexion-extension are decreased by decreasing disc height. By decreasing the initial water content, axial displacement, facet force and fluid loss were all reduced. Endplate sclerosis did not have a significant influence on the calculated results. Depressurisation determined an increase of the disc effective stress, possibly inducing failure. Degenerative instability was not calculated in any simulations.     

The effect of degenerative morphological changes of the intervertebral disc on the lumbar spine biomechanics: a poroelastic finite element investigation

Intervertebral disc degeneration involves changes in the spinal anatomical structures. The mechanical relevance of the following changes was investigated: disc height, endplate sclerosis, disc water content, permeability and depressurisation. A poroelastic nonlinear finite element model of the L4–L5 human spine segments was employed. Loads represented a daily cycle (500 N compression combined with flexion–extension motion for 16 h followed by 200 N compression for 8 h). In non-degenerative conditions, the model predicted a diurnal axial displacement of 1.32 mm and a peak intradiscal pressure of 0.47 MPa. Axial displacement, facet force and range of motion in flexion–extension are decreased by decreasing disc height. By decreasing the initial water content, axial displacement, facet force and fluid loss were all reduced. Endplate sclerosis did not have a significant influence on the calculated results. Depressurisation determined an increase of the disc effective stress, possibly inducing failure. Degenerative instability was not calculated in any simulations.     

Configuration of anchorage holes affects fixation of the acetabular component in cemented total hip replacement--a finite element study

Our survey of current practice among UK orthopaedic surgeons shows wide variations in fixation techniques. The aim of this study, is to investigate the effect of drilling different configurations of anchorage holes in the acetabulum on implant stability. To avoid variables that could incur during in vitro testing, we used commercially available COSMOS finite element analysis package to investigate the stress distributions, deformations, and strains on the cement mantle when drilling three large anchorage holes and six smaller ones, with straight and rounded cement pegs. The results, which are in line with our in vitro studies on simulated reconstructed acetabulae, indicate better stability of the acetabular component when three larger holes than six smaller holes are drilled and when the necks of the anchorage holes are rounded. The longevity of total hip replacements could be improved by drilling three large anchorage holes, rather than many smaller ones, as initially proposed by Charnley.     

Discrete element analysis is a valid method for computing joint contact stress in the hip before and after acetabular fracture

Evaluation of abnormalities in joint contact stress that develop after inaccurate reduction of an acetabular fracture may provide a potential means for predicting the risk of developing post-traumatic osteoarthritis. Discrete element analysis (DEA) is a computational technique for calculating intra-articular contact stress distributions in a fraction of the time required to obtain the same information using the more commonly employed finite element analysis technique. The goal of this work was to validate the accuracy of DEA-computed contact stress against physical measurements of contact stress made in cadaveric hips using Tekscan sensors. Four static loading tests in a variety of poses from heel-strike to toe-off were performed in two different cadaveric hip specimens with the acetabulum intact and again with an intentionally malreduced posterior wall acetabular fracture. DEA-computed contact stress was compared on a point-by-point basis to stress measured from the physical experiments. There was good agreement between computed and measured contact stress over the entire contact area (correlation coefficients ranged from 0.88 to 0.99). DEA-computed peak contact stress was within an average of 0.5 MPa (range 0.2–0.8 MPa) of the Tekscan peak stress for intact hips, and within an average of 0.6 MPa (range 0–1.6 MPa) for fractured cases. DEA-computed contact areas were within an average of 33% of the Tekscan-measured areas (range: 1.4–60%). These results indicate that the DEA methodology is a valid method for accurately estimating contact stress in both intact and fractured hips.     

The influence of ligament modelling strategies on the predictive capability of finite element models of the human knee joint

In finite element (FE) models knee ligaments can represented either by a group of one-dimensional springs, or by three-dimensional continuum elements based on segmentations. Continuum models closer approximate the anatomy, and facilitate ligament wrapping, while spring models are computationally less expensive. The mechanical properties of ligaments can be based on literature, or adjusted specifically for the subject. In the current study we investigated the effect of ligament modelling strategy on the predictive capability of FE models of the human knee joint. The effect of literature-based versus specimen-specific optimized material parameters was evaluated. Experiments were performed on three human cadaver knees, which were modelled in FE models with ligaments represented either using springs, or using continuum representations. In spring representation collateral ligaments were each modelled with three and cruciate ligaments with two single-element bundles. Stiffness parameters and pre-strains were optimized based on laxity tests for both approaches. Validation experiments were conducted to evaluate the outcomes of the FE models.Models (both spring and continuum) with subject-specific properties improved the predicted kinematics and contact outcome parameters. Models incorporating literature-based parameters, and particularly the spring models (with the representations implemented in this study), led to relatively high errors in kinematics and contact pressures. Using a continuum modelling approach resulted in more accurate contact outcome variables than the spring representation with two (cruciate ligaments) and three (collateral ligaments) single-element-bundle representations. However, when the prediction of joint kinematics is of main interest, spring ligament models provide a faster option with acceptable outcome.