[論文レビュー] 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）になると既知の2次元解に還元される。
• 導出された接触圧力および接触領域発展の式は正確であり、長時間にわたる持続的荷重下における二相性関節軟骨の応答を正確にモデル化可能である。