삼차원 이중상태 연골 접촉 문제의 폐쇄형 해법

초록

본 논문은 강체인 피질골 위에 부착된 연골을 고체‑유체 이중상태(biphasic) 모델로 가정하고, 타원형 포물면 형태의 피질골을 경계조건으로 삼아 삼차원 단일 접촉 문제에 대한 폐쇄형 해석식을 도출한다. 장시간 하중 조건과 점증적 하중 증가 상황에서도 적용 가능한 해를 제공함으로써 연골의 변형·압력 분포를 정량적으로 예측한다.

상세 요약

연골은 고체 매트릭스와 관통성 체액으로 구성된 이중상태 물질이며, 하중이 가해질 때 체액이 압축을 완화하고 고체 매트릭스가 변형을 담당한다는 점에서 기존의 단일상태 탄성 모델보다 복잡한 거동을 보인다. 저자들은 이러한 특성을 반영하기 위해 마스코위츠(Murphy)와 마이클스(Mich) 등에서 제안된 2‑phase 이론을 채택하고, 피질골을 변형이 없는 강체로 가정함으로써 연골‑피질골 복합체의 경계조건을 단순화하였다. 특히 피질골 표면을 타원형 포물면(elliptic paraboloid) 형태로 모델링함으로써 실제 관절면의 비대칭 곡률을 수학적으로 구현하였다.

수학적 전개는 라플라스 변환을 이용한 시간 영역 해석과, 구면조화함수와 타원좌표계에서의 고유함수 전개를 결합한 혼합 해법을 사용한다. 이 과정에서 접촉 영역을 타원형으로 가정하고, 접촉 압력은 다항식 형태로 전개되며, 경계조건(무마찰, 비침투성)과 연골의 체적 보존 법칙을 동시에 만족하도록 계수를 결정한다. 핵심은 장시간(steady‑state) 하중 상황에서도 체액 흐름이 포화 상태에 도달하면서 압력-변형 관계가 일정한 형태로 수렴한다는 점이다. 이를 통해 저주파 하중이나 지속적인 압축 하중에 대한 연골의 응답을 폐쇄형 식으로 표현할 수 있었다.

또한, 저자들은 하중이 점진적으로 증가하는 경우를 고려하여, 하중 함수가 시간에 따라 선형적으로 증가하는 상황에서도 해가 연속적으로 변하도록 파라미터화하였다. 이때, 접촉 반경과 압력 중심의 이동을 동시에 추정함으로써, 실제 관절 움직임에서 발생하는 비대칭 접촉 현상을 정량화할 수 있다. 결과적으로, 해는 연골 두께, 고체 탄성계수, 체액 압축성, 그리고 피질골 곡률 반경 등 물성 파라미터와 직접적인 연관성을 갖는다.

이러한 해법은 기존의 수치해석(FEM) 대비 계산 비용이 현저히 낮으며, 파라미터 스터디나 임상 데이터와의 빠른 비교에 유리하다. 다만, 강체 피질골 가정과 접촉 면을 타원형으로 제한한 점은 실제 관절에서 복잡한 불규칙 형상이나 연골‑연골 간의 상대 움직임을 완전히 포착하지 못한다는 한계가 있다. 향후 연구에서는 비강체 피질골 모델링, 비선형 체액 흐름, 그리고 마찰 효과를 포함한 확장이 필요할 것으로 보인다.