[论文解读] A closed-form solution of the three-dimensional contact problem for biphasic cartilage layers

Ivan Argatov, Gennady Mishuris|arXiv (Cornell University)|Sep 22, 2010

Osteoarthritis Treatment and Mechanisms参考文献 12被引用 3

一句话总结

本文针对长期载荷下粘弹性软骨层与刚性椭圆抛物面骨之间三维接触问题，提出了一种精确的闭式解析解。基于广义双相模型和渐近位移关系，推导出接触压力、接触面积演化及压力分布的显式表达式，将先前的轴对称解推广至三维椭圆几何构型，其时间依赖行为由流体流动与固体基质力学共同决定。

ABSTRACT

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.

研究动机与目标

• 将现有的双相软骨轴对称接触模型扩展至三维椭圆几何构型。
• 为长期载荷下的接触压力与接触面积演化发展解析解。
• 在双相软骨层中考虑随时间变化的流体流动与固体基质变形。
• 将Wu等人（1997）的先前解推广至非轴对称、三维椭圆接触构型。
• 以闭式表达式提供接触压力、接触域形状及位移参数的精确表达式。

提出的方法

• 采用三维单面接触条件，建立附着于刚性椭圆抛物面骨的双相软骨层之间的接触问题模型。
• 应用Ateshian等人（1994）的渐近位移关系，模拟由流体流动与固体相刚度引起的时变垂直位移。
• 推导出一个积分微分方程（式1.7），将接触压力、位移与时间关联，引入流体渗透率与固体剪切模量参数。
• 通过引入算子K的变换，将问题简化为可解形式，进而获得压力分布的闭式解。
• 利用格林函数方法与海维赛德阶跃函数，建模接触区的动态扩展过程，求解所得方程。
• 推导出接触压力分布、接触面积半轴a(t)与b(t)，以及压力参数p₀(t)的显式表达式，其形式由材料与几何参数决定。

实验结果

研究问题

• RQ1 在长期载荷条件下，如何对双相软骨层的三维接触问题进行解析求解？
• RQ2 双相软骨在时变椭圆接触面积上的接触压力分布的确切形式是什么？
• RQ3 接触区如何随时间扩展？各空间点接触的起始时刻由什么决定？
• RQ4 流体流动与固体基质刚度在塑造压力分布与接触面积演化中起何种作用？
• RQ5 所推导的三维解如何推广先前的轴对称解？

主要发现

• 接触压力分布以闭式表达为P(x₁,x₂,t) = (1 - x₁²/a²(t) - x₂²/b²(t))² × p₀(t) 减去一个考虑流体流动历史的时间积分项。
• 接触面积随时间扩展，某点(x₁,x₂)的接触起始时间t*由接触轮廓到达该位置的条件确定。
• 对于阶跃加载，接触起始时间t*(x₁,x₂)由包含初始接触半径与材料参数的闭式表达式明确给出（式5.51）。
• 压力参数p₀(t)与位移δ₀(t)通过由力平衡与材料特性导出的时间相关关系相联系。
• 该解推广了先前的轴对称结果，并在几何对称时（s=1）退化为已知的二维情形。
• 所推导的接触压力与面积演化表达式为精确解，适用于长时间范围，可实现对双相软骨在持续载荷下响应的精确建模。