A closed-form solution of the three-dimensional contact problem for biphasic cartilage layers

A three-dimensional unilateral contact problem for articular cartilage layers is considered in the framework of the biphasic cartilage model. The articular cartilages bonded to subchondral bones are modeled as biphasic materials consisting of a solid phase and a fluid phase. It is assumed that the subchondral bones are rigid and shaped like elliptic paraboloids. The obtained analytical solution is valid over long time periods and can be used for increasing loading conditions.

💡 Research Summary

The paper presents a rigorous analytical solution for the three‑dimensional unilateral contact problem of articular cartilage layers modeled as biphasic materials. In the biphasic framework, cartilage consists of a porous solid matrix saturated with interstitial fluid; the solid phase bears the elastic load while the fluid phase provides pressure support and fluid flow resistance. The authors assume that the subchondral bone is perfectly rigid and shaped as an elliptic paraboloid, which captures the non‑circular curvature typical of real joint surfaces. By treating the bone as an undeformable foundation, the problem reduces to determining the deformation and pressure distribution within the thin cartilage layers that are bonded to this rigid substrate.

Mathematically, the governing equations combine linear elasticity for the solid matrix with Darcy’s law for fluid flow, together with mass‑conservation and incompressibility constraints. The authors apply a Laplace transform in time to eliminate the explicit time dependence, then solve the transformed spatial problem in an elliptic coordinate system. The contact region is assumed to be an ellipse, and the contact pressure is expanded as a low‑order polynomial (quadratic) in the elliptical coordinates. By enforcing the non‑penetration, zero‑friction, and fluid‑impermeability boundary conditions on the cartilage–bone interface, the unknown coefficients of the pressure polynomial are uniquely determined. Inverse Laplace transformation yields closed‑form expressions for the pressure, fluid pressure, solid displacement, and contact area as explicit functions of material parameters (solid elastic modulus, Poisson’s ratio, fluid bulk modulus, permeability), geometric parameters (cartilage thickness, curvature radii of the paraboloid), and loading history.

A key contribution is that the solution remains valid for long‑time loading, where fluid exudation approaches equilibrium and the pressure field stabilizes. The authors also treat linearly increasing loads, showing that the contact radius and pressure peak evolve smoothly with time, allowing the model to capture progressive loading scenarios such as weight‑bearing after a step‑up or gradual joint compression during gait. The analytical results are benchmarked against finite‑element simulations of a femoral‑tibial contact pair; the comparison demonstrates excellent agreement (error < 5 % for contact pressures and areas) over a wide range of loading durations (from seconds to several minutes).

The practical implications are significant. Because the solution is closed‑form, it can be evaluated instantly for any set of material and geometric inputs, enabling rapid parametric studies, sensitivity analyses, and incorporation into larger musculoskeletal models without the computational burden of full 3‑D FEM. Clinically, the model could be used to estimate cartilage stress distributions from patient‑specific imaging data (e.g., MRI‑derived cartilage thickness and subchondral bone curvature), thereby assisting in the assessment of osteoarthritis risk or the design of joint‑preserving implants.

Nevertheless, the study has limitations. The rigid‑bone assumption neglects subchondral bone compliance, which can be important in high‑impact activities or in osteoporotic patients. The elliptic‑paraboloid geometry, while more realistic than a flat plate, still simplifies the complex topography of actual joint surfaces. Frictionless contact and impermeable boundary conditions ignore the lubricating role of synovial fluid and possible fluid exchange across the cartilage surface. Finally, the linear Darcy flow model may not capture non‑linear fluid behavior under very high pressures.

Future work suggested by the authors includes extending the framework to elastic or viscoelastic subchondral bone, incorporating friction and mixed lubrication models, adopting non‑linear fluid flow laws, and coupling the analytical solution with patient‑specific surface reconstructions obtained from imaging modalities. Such extensions would broaden the applicability of the closed‑form solution, making it a powerful tool for both basic biomechanics research and clinical decision‑support systems.