[논문 리뷰] A closed-form solution of the three-dimensional contact problem for biphasic cartilage layers

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

Osteoarthritis Treatment and Mechanisms참고 문헌 12인용 수 3

한 줄 요약

이 논문은 장기 하중 조건 하에서 강성 타원형 포물면 뼈에 부착된 이상성 연골 층 간의 3차원 접촉 문제에 대한 정확한 폐쇄형 해석적 해를 제시한다. 일반화된 이상성 모델과 점점 가까워지는 변위 관계를 사용하여, 접촉 압력, 접촉 영역 변화, 압력 분포에 대한 명시적 표현을 유도하며, 이는 이전의 축대칭 해를 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.

연구 동기 및 목표

• 이상성 연골의 기존 축대칭 접촉 모델을 3차원 타원형 기하학으로 확장하기 위해.
• 장기 하중 조건 하에서 접촉 압력과 접촉 영역 변화에 대한 해석적 해를 개발하기 위해.
• 이상성 연골 층 내 시간에 따라 변화하는 유체 흐름과 고체 매트릭스 변형을 고려하기 위해.
• Wu 등 (1997)의 이전 해를 축대칭이 아닌 3차원 타원형 접촉 구형으로 일반화하기 위해.
• 폐쇄형으로 표현된 접촉 압력, 접촉 영역 형태, 변위 매개변수의 정확한 표현을 제공하기 위해.

제안 방법

• 강성 타원형 포물면 뼈에 부착된 두 개의 이상성 연골 층 간의 3차원 단일 접촉 조건을 사용하여 접촉 문제를 수립한다.
• Ateshian 등 (1994)의 점점 가까워지는 변위 관계를 적용하여, 유체 흐름과 고체 상의 강성에 의해 유도되는 시간에 따라 변화하는 수직 변위를 모델링한다.
• 유체 투과도와 고체 전단 모odulus를 포함하여, 접촉 압력, 변위, 시간 간의 연관을 나타내는 적분미분방정식 (식 1.7)을 유도한다.
• 문제를 해결 가능한 형태로 줄이기 위해 연산자 K를 사용한 변환을 도입하여, 압력 분포에 대한 폐쇄형 해를 도출한다.
• 그리고 그 결과 방정식을 그린 함수 방법과 헤비사이드 단위계단 함수를 사용하여 해석하여, 접촉 영역의 동적 확장을 모델링한다.
• 접촉 압력 분포, 접촉 영역의 반경 a(t)와 b(t), 그리고 압력 매개변수 p₀(t)에 대한 명시적 표현을 재료 및 기하학적 매개변수로 유도한다.

실험 결과

연구 질문

• RQ1 장기 하중 조건 하에서 이방성 연골 층 간의 3차원 접촉 문제는 어떻게 해석적으로 해결할 수 있는가?
• RQ2 이상성 연골에 대해 시간에 따라 변화하는 타원형 접촉 영역에서 접촉 압력 분포의 정확한 형태는 무엇인가?
• RQ3 하중 증가에 따라 접촉 영역은 어떻게 시간에 따라 확장되는가? 그리고 각 공간 점에서 접촉이 시작되는 조건은 무엇인가?
• RQ4 유체 흐름과 고체 매트릭스의 강성이 압력 분포와 접촉 영역 변화에 미치는 영향은 무엇인가?
• RQ5 유도된 3차원 해는 이전의 축대칭 해를 어떻게 일반화하는가?

주요 결과

• 접촉 압력 분포는 폐쇄형으로 주어지며, 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차원 해로 축소된다.
• 유도된 접촉 압력 및 영역 변화에 대한 표현식은 정확하며 장기간 동안 유효하여, 지속 하중 하에서 이방성 연골 반응을 정확하게 모델링할 수 있다.