Combined finite-element and rigid-body analysis of human jaw joint dynamics

The jaw joint plays a crucial role in human mastication. It acts as a guidance for jaw movements and as a fulcrum for force generation. The joint is subjected to loading which causes tensions and deformations in its cartilaginous structures. These are assumed to be a major determinant for development, maintenance and also degeneration of the joint. To analyze the distribution of tensions and deformations in the cartilaginous structures of the jaw joint during jaw movement, a dynamical model of the human masticatory system has been constructed. Its movements are controlled by muscle activation. The articular cartilage layers and articular disc were included as finite-element (FE) models. As this combination of rigid-body and FE modeling had not been applied to musculoskeletal systems yet, its benefits and limitations were assessed by simulating both unloaded and loaded jaw movements. It was demonstrated that joint loads increase with muscle activation, irrespective of the external loads. With increasing joint load, the size of the stressed area of the articular surfaces was enlarged, whereas the peak stresses were much less affected. The results suggest that the articular disc enables distribution of local contact stresses over a much wider area of the very incongruent articular surfaces by transforming compressive principal stress into shear stress.     

Finite element simulations of a focal knee resurfacing implant applied to localized cartilage defects in a sheep model

Articular resurfacing metal implants have recently been tested in animal models to treat full thickness localized articular cartilage defects, showing promising results. However, the mechanical behavior of cartilage surrounding the metal implant has not been studied yet as it is technically challenging to measure in vivo contact areas, pressures, stresses and deformations from the metal implant. Therefore, we implemented a detailed numerical finite element model by approximating one of the condyles of the sheep tibiofemoral joint and created a defect of specific size to accommodate the implant. Using this model, the mechanical behavior of the surrounding of metal implant was studied. The model showed that the metal implant plays a significant role in the force transmission. Two types of profiles were investigated for metal implant. An implant with a double-curved profile, i.e., a profile fully congruent with the articular surfaces in the knee, gives lower contact pressures and stresses at the rim of the defect than the implant with unicurved spherical profile. The implant should be placed at a certain distance into the cartilage to avoid damage to opposing biological surface. Too deep positions, however, lead to high shear stresses in the cartilage edges around the implant. Mechanical sealing was achieved with a wedge shape of the implant, also useful for biochemical sealing of cartilage edges at the defect.     

Perspectives on chondrocyte mechanobiology and osteoarthritis

Osteoarthritis, the clinical syndrome of joint pain and dysfunction due to joint degeneration, is among the most frequent and symptomatic medical problems for middle aged and older people, and it is the most common cause of long term disability in most populations of people over 65. Currently there are no effective methods of preventing or curing osteoarthritis. Post-traumatic OA, the joint degeneration, pain and dysfunction that develop following joint injury, is the form of OA that is most directly related to elevated articular surface contact stress. However, mechanical stress that exceeds the tolerance of the articular surface can cause or accelerate the progression of joint degeneration in all individuals and in all synovial joints. In some patients, decreasing mechanical forces on degenerated joint surfaces stimulates formation of a new biologic articular surface. The advances in understanding of the effects of mechanical forces on chondrocytes and cartilage presented and discussed at the 4th Symposium on Mechanobiology: Cartilage and Chondrocyte will help in the efforts to develop new methods of preventing and treating osteoarthritis.     

Stresses in the local collagen network of articular cartilage: a poroviscoelastic fibril-reinforced finite element study

Osteoarthritis (OA) is a multifactorial disease, resulting in diarthrodial joint wear and eventually destruction. Swelling of cartilage, which is proportional to the amount of collagen damage, is an initial event of cartilage degeneration, so damage to the collagen fibril network is likely to be one of the earliest signs of OA cartilage degeneration. We propose that the local stresses and strains in the collagen fibrils, which cause the damage, cannot be determined dependably without taking the local arcade-like collagen-fibril structure into account. We investigate this using a poroviscoelastic fibril-reinforced FEA model. The constitutive fibril properties were determined by fitting numerical data to experimental results of unconfined compression and indentation tests on samples of bovine patellar articular cartilage. It was demonstrated that with this model the stresses and strains in the collagen fibrils can be calculated. It was also exhibited that fibrils with different orientations at the same location can be loaded differently, depending on the local architecture of the collagen network. To the best of our knowledge, the present model is the first that can account for these features. We conclude that the local stresses and strains in the articular cartilage are highly influenced by the local morphology of the collagen-fibril network.     

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.     

Dependence of nanoscale friction and adhesion properties of articular cartilage on contact load

Boundary lubrication of articular cartilage by conformal, molecularly thin films reduces friction and adhesion between asperities at the cartilage-cartilage contact interface when the contact conditions are not conducive to fluid film lubrication. In this study, the nanoscale friction and adhesion properties of articular cartilage from typical load-bearing and non-load-bearing joint regions were studied in the boundary lubrication regime under a range of physiological contact pressures using an atomic force microscope (AFM). Adhesion of load-bearing cartilage was found to be much lower than that of non-load-bearing cartilage. In addition, load-bearing cartilage demonstrated steady and low friction coefficient through the entire load range examined, whereas non-load-bearing cartilage showed higher friction coefficient that decreased nonlinearly with increasing normal load. AFM imaging and roughness calculations indicated that the above trends in the nanotribological properties of cartilage are not due to topographical (roughness) differences. However, immunohistochemistry revealed consistently higher surface concentration of boundary lubricant at load-bearing joint regions. The results of this study suggest that under contact conditions leading to joint starvation from fluid lubrication, the higher content of boundary lubricant at load-bearing cartilage sites preserves synovial joint function by minimizing adhesion and wear at asperity microcontacts, which are precursors for tissue degeneration.     

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.     

Hemiarthroplasty of hip joint: An experimental validation using porcine acetabulum

Biphasic properties of articular cartilage allow it to be an excellent bearing material and have been studied through several simplified experiments as well as finite element modelling. However, three-dimensional biphasic finite element (FE) models of the whole joint are rare. The current study was carried out to experimentally validate FE methodology for modelling hemiarthroplasty. Material properties such as equilibrium elastic modulus and permeability of porcine acetabular cartilage were initially derived by curve-fitting an experimental deformation curve with that obtained using FE. These properties were then used in the hemiarthroplasty hip joint modelling. Each porcine acetabular cup was loaded with 400N using a 34mm diameter CoCr femoral head. A specimen-specific FE model of each acetabular cup was created using μCT and a series of software processes. Each model was analysed under conditions similar to those tested experimentally. Contact stresses and contact areas predicted by the model, immediately after loading, were then compared with the corresponding experimentally measured values. Very high peak contact stresses (maximum experimental: 14.09MPa) were recorded. A maximum difference of 12.42% was found in peak contact stresses. The corresponding error for contact area was 20.69%. Due to a fairly good agreement in predicted and measured values of contact stresses and contact areas, the integrated methodology developed in this study can be used as a basis for future work. In addition, FE predicted total fluid load support was around 80% immediately after loading. This was lower than that observed in conforming contact problems involving biphasic cartilage and was due to a smaller local contact area and variable clearance making fluid exudation easier.     

Physical validation of a patient-specific contact finite element model of the ankle

A validation study was conducted to determine the extent to which computational ankle contact finite element (FE) results agreed with experimentally measured tibio-talar contact stress. Two cadaver ankles were loaded in separate test sessions, during which ankle contact stresses were measured with a high-resolution (Tekscan) pressure sensor. Corresponding contact FE analyses were subsequently performed for comparison. The agreement was good between FE-computed and experimentally measured mean (3.2% discrepancy for one ankle, 19.3% for the other) and maximum (1.5% and 6.2%) contact stress, as well as for contact area (1.7% and 14.9%). There was also excellent agreement between histograms of fractional areas of cartilage experiencing specific ranges of contact stress. Finally, point-by-point comparisons between the computed and measured contact stress distributions over the articular surface showed substantial agreement, with correlation coefficients of 90% for one ankle and 86% for the other. In the past, general qualitative, but little direct quantitative agreement has been demonstrated with articular joint contact FE models. The methods used for this validation enable formal comparison of computational and experimental results, and open the way for objective statistical measures of regional correlation between FE-computed contact stress distributions from comparison articular joint surfaces (e.g., those from an intact versus those with residual intra-articular fracture incongruity).     

Lubricin: a novel potential biotherapeutic approaches for the treatment of osteoarthritis

Osteoarthritis (OA) is a multi-factor disorder of sinovial joints, which characterized by escalated degeneration and loss of articular cartilage. Treatment of OA is a critical unmet need in medicine for regeneration of damaged articular cartilage in elderly. On the other hand, lubricin, a glycoprotein specifically synthesized by chondrocytes located at the surface of articular cartilage, has been shown to provide boundary lubrication of congruent articular surfaces under conditions of high contact pressure and near zero sliding speed. Lubrication of these surfaces is critical to normal joint function, while different gene expressions of lubricin had been found in the synovium of rheumatoid arthritis (RA) and OA. Moreover, mutations or lacking of lubricin gene have been shown to link to the joint disease such as camptodactyly-arthropathy-coxa vara-pericarditis syndrome (CACP), synovial hyperplasia and failure of joint function, suggesting an important role of lubricin in the pathogenesis of these joint disease. Recent studies demonstrate that administration with recombinant lubricin in the joint cavity would be effective in the prevention of cartilage degeneration in animal OA models. Therefore, a treatment with lubricin which would protect cartilage in vivo would be desirable. This article reviews recent findings with regard to the possible role of lubricin in the progression of OA, and further discusses lubricin as a novel potential biotherapeutic approaches for the treatment of OA.     

Post-traumatic osteoarthritis: the role of stress induced chondrocyte damage

Post-traumatic osteoarthritis is the form of osteoarthritis (OA) that develops following joint injury. Although its end-stage is indistinguishable from idiopathic OA, many patients with post-traumatic OA are younger than those with idiopathic OA, and they have a well-defined precipitating insult. Clinical and experimental studies suggest that excessive acute impact energy or chronic mechanical overload cause the degeneration of the articular surface responsible for post-traumatic OA. Yet, the mechanisms by which excessive mechanical force causes OA remain unknown. For these reasons it has not been possible to develop effective methods of preventing or decreasing the risk of post-traumatic OA. We hypothesized that mechanical loading that exceeds the tolerance of the articular surface causes chondrocyte damage due to oxidative stress. Our in vitro tests of human articular cartilage samples showed that shear stress causes chondrocyte death and that anti-oxidants decrease the shear stress induced cell death. These observations suggest that specific patterns of loading are particularly damaging to articular surfaces and that improved treatments of joint injuries may include mechanical methods of minimizing shear stresses and biologic methods of minimizing oxidative damage.     

Cartilage thickness distribution affects computational model predictions of cervical spine facet contact parameters

With motion-sparing disk replacement implants gaining popularity as an alternative to anterior cervical discectomy and fusion (ACDF) for the treatment of certain spinal degenerative disorders, recent laboratory investigations have studied the effects of disk replacement and implant design on spinal kinematics and kinetics. Particularly relevant to cervical disk replacement implant design are any postoperative changes in solid stresses or contact conditions in the articular cartilage of the posterior facets, which are hypothesized to lead to adjacent-level degeneration. Such changes are commonly investigated using finite element methods, but significant simplification of the articular geometry is generally employed. The impact of such geometric representations has not been thoroughly investigated. In order to assess the effects of different models of cartilage geometry on load transfer and contact pressures in the lower cervical spine, a finite element model was generated using cadaver-based computed tomography imagery. Mesh resolution was varied in order to establish model convergence, and cadaveric testing was undertaken to validate model predictions. The validated model was altered to include four different geometric representations of the articular cartilage. Model predictions indicate that the two most common representations of articular cartilage geometry result in significant reductions in the predictive accuracy of the models. The two anatomically based geometric models exhibited less computational artifact, and relatively minor differences between them indicate that contact condition predictions of spatially varying thickness models are robust to anatomic variations in cartilage thickness and articular curvature. The results of this work indicate that finite element modeling efforts in the lower cervical spine should include anatomically based and spatially varying articular cartilage thickness models. Failure to do so may result in loss of fidelity of model predictions relevant to investigations of physiological import.     

Contact analysis of biphasic transversely isotropic cartilage layers and correlations with tissue failure.

Failure of articular cartilage has been investigated experimentally and theoretically, but there is only partial agreement between observed failure and predicted regions of peak stresses. Since trauma and repetitive stress are implicated in the etiopathogenesis of osteoarthritis, it is important to develop cartilage models which correctly predict sites of high stresses. Cartilage is anisotropic and inhomogeneous, though it has been difficult to incorporate these complexities into engineering analyses. The objectives of this study are to demonstrate that a transversely isotropic, biphasic model of cartilage can provide agreement between predicted regions of high stresses and observed regions of cartilage failure and that with transverse isotropy cartilage stresses are more sensitive to convexity and concavity of the surfaces than with isotropy. These objectives are achieved by solving problems of diarthrodial joint contact by the finite-element method. Results demonstrate that transversely isotropic models predict peak stresses at the cartilage surface and the cartilage-bone interface, in agreement with sites of fissures following impact loading; isotropic models predict peak stresses only at the cartilage-bone interface. Also, when convex cartilage layers contacted concave layers in this study, the highest tensile stresses occur in the convex layer for transversely isotropic models; no such differences are found with isotropic models. The significance of this study is that it establishes a threshold of modeling complexity for articular cartilage that provides good agreement with experimental observations under impact loading and that surface curvatures significantly affect stress and strain within cartilage when using a biphasic transversely isotropic model.     

An approach for the stress analysis of transversely isotropic biphasic cartilage under impact load.

Stress analysis of contact models for isotropic articular cartilage under impacting loads shows high shear stresses at the interface with the subchondral bone and normal compressive stresses near the surface of the cartilage. These stress distributions are not consistent, with lesions observed on the cartilage surface of rabbit patellae from blunt impact, for example, to the patello-femoral joint. The purpose of the present study was to analyze, using the elastic capabilities of a finite element code, the stress distribution in more morphologically realistic transversely isotropic biphasic contact models of cartilage. The elastic properties of an incompressible material, equivalent to those of the transversely isotropic biphasic material at time zero, were derived algebraically using stress-strain relations. Results of the stress analysis showed the highest shear stresses on the surface of the solid skeleton of the cartilage and tensile stresses in the zone of contact. These results can help explain the mechanisms responsible for surface injuries observed during blunt insult experiments.     

Three-dimensional mathematical model analysis of the patellofemoral joint

This paper is concerned with a mathematical model analysis of the patellofemoral joint in the human knee, taking into account the articular surface geometry and mechanical properties of the ligament. It was made by the application of a computer-aided design theory (previously studied) and it was possible to express the articular surface geometries in a mathematical formulation and hence elucidate the joint movement mechanics. This method was then applied to a three-dimensional geometrical model of the patellofemoral joint. For the modelling of tendofemoral contact at large angles of knee flexion, the geodestic line theory was adopted. Applying the Newton-Raphson method and the Runge-Kutta Gil method to the model, variables such as patellar attitudes, patellofemoral contact force and tensile force of the patellar ligament for various knee flexion angles were computed. Applying the Hertzian elastic theory, contact stress was also computed. These results showed good agreement with the previously reported experimental results. As an application for the model, some parameter analyses were performed in terms of the contact stress variations and compared with those of the normal knee. The simulation results indicated that both the Q-angle increase and decrease increased contact stress, the patella alta showed undulating variations of stress while the patella infera showed little change of stress, and the tibial tuberositas elevation showed 20-30% reduction of stress.     

Evaluation of a subject-specific finite-element model of the equine metacarpophalangeal joint under physiological load

The equine metacarpophalangeal (MCP) joint is frequently injured, especially by racehorses in training. Most injuries result from repetitive loading of the subchondral bone and articular cartilage rather than from acute events. The likelihood of injury is multi-factorial but the magnitude of mechanical loading and the number of loading cycles are believed to play an important role. Therefore, an important step in understanding injury is to determine the distribution of load across the articular surface during normal locomotion. A subject-specific finite-element model of the MCP joint was developed (including deformable cartilage, elastic ligaments, muscle forces and rigid representations of bone), evaluated against measurements obtained from cadaver experiments, and then loaded using data from gait experiments. The sensitivity of the model to force inputs, cartilage stiffness, and cartilage geometry was studied. The FE model predicted MCP joint torque and sesamoid bone flexion angles within 5% of experimental measurements. Muscle–tendon forces, joint loads and cartilage stresses all increased as locomotion speed increased from walking to trotting and finally cantering. Perturbations to muscle–tendon forces resulted in small changes in articular cartilage stresses, whereas variations in joint torque, cartilage geometry and stiffness produced much larger effects. Non-subject-specific cartilage geometry changed the magnitude and distribution of pressure and the von Mises stress markedly. The mean and peak cartilage stresses generally increased with an increase in cartilage stiffness. Areas of peak stress correlated qualitatively with sites of common injury, suggesting that further modelling work may elucidate the types of loading that precede joint injury and may assist in the development of techniques for injury mitigation.     

Analysis of acute mechanical insult in an animal model of post-traumatic osteoarthrosis.

Chronic degeneration of articular cartilage and bone in a rabbit model of post-traumatic osteoarthrosis has been hypothesized to occur due to acute stresses that exceed a threshold for injury. In this study, we impacted the rabbit patellofemoral joint at low and high intensities. High-intensity impacts produced degenerative changes in the joint, such as softening of retropatellar cartilage, as measured by indentation, an increase in histopathology of the cartilage, and an increase in thickness of subchondral bone underlying the cartilage. Low-intensity impacts did not cause these progressive changes. These data suggest that low-intensity impacts produced acute tissue stresses below the injury threshold, while high-intensity impacts produced stresses that exceeded the threshold for disease pathogenesis. This study begins to identify "safe" and "unsafe" ranges of acute tissue stress, using the rabbit patella, which may have future utility in the design of injury prevention devices for the human.     

Glenohumeral joint cartilage contact in the healthy adult during scapular plane elevation depression with external humeral rotation

The shoulder (glenohumeral) joint has the greatest range of motion of all human joints; as a result, it is particularly vulnerable to dislocation and injury. The ability to non-invasively quantify in-vivo articular cartilage contact patterns of joints has been and remains a difficult biomechanics problem. As a result, little is known about normal in-vivo glenohumeral joint contact patterns or the consequences that surgery has on altering them. In addition, the effect of quantifying glenohumeral joint contact patterns by means of proximity mapping, both with and without cartilage data, is unknown. Therefore, the objectives of this study are to (1) describe a technique for quantifying in-vivo glenohumeral joint contact patterns during dynamic shoulder motion, (2) quantify normal glenohumeral joint contact patterns in the young healthy adult during scapular plane elevation depression with external humeral rotation, and (3) compare glenohumeral joint contact patterns determined both with and without articular cartilage data. Our results show that the inclusion of articular cartilage data when quantifying in-vivo glenohumeral joint contact patterns has significant effects on the anterior–posterior contact centroid location, the superior–inferior contact centroid range of travel, and the total contact path length. As a result, our technique offers an advantage over glenohumeral joint contact pattern measurement techniques that neglect articular cartilage data. Likewise, this technique may be more sensitive than traditional 6-Degree-of-Freedom (6-DOF) joint kinematics for the assessment of overall glenohumeral joint health. Lastly, for the shoulder motion tested, we found that glenohumeral joint contact was located on the anterior–inferior glenoid surface.     

The effect of collagen reinforcement in the behaviour of the temporomandibular joint disc

In this paper, the influence of collagen fibres in the behaviour of the temporomandibular joint disc is studied. A three-dimensional finite element model of the joint is developed from a set of medical images. The model comprises the mandible, part of the cranium and both temporomandibular joints. Joints have been considered to be composed of the articular discs and the temporomandibular ligaments. A fibre-reinforced porohyperelastic model was used to study the response under clenching of the fibrocartilage that composes the articular disc. This was divided in an intermediate zone, and two bands, an anterior and other posterior, in order to define the orientation of collagen fibres. The study demonstrates that the introduction of collagen fibres in the biphasic behaviour of the articular disc implies for a prescribed displacement not only an increase of the pressurization in the tissue, but also higher stresses in the anterior and posterior bands, as well as in the lateral zone of the disc. Thus, modelling the disc as an isotropic solid matrix leads in this case to an overestimation of the stresses in the intermediate zone, an underestimation of the pore pressure in this area, and an underestimation of the stresses in the rest of the disc.