Uniaxial poroelastic tendon model with crimped fibre recruitment

기술 요약: 크림프 섬유 모집을 포함한 단축성 다공 탄성 건 모델

문제 정의
건(Tendon)은 주로 콜라겐 피브릴, 비콜라겐성 기질(NCM), 그리고 물(55–70%)로 구성된 복잡한 연조직이다. 피브릴이 펴지면서 하중을 모집하기 시작하는 '크림프(crimped)' 현상의 기계적 역할은 실험적으로 잘 입증되어 있으나, 기존의 수학적 모델들은 이러한 모집 행동을 다공 탄성 프레임워크 내에 통합하는 데 종종 실패해 왔다. 기존의 다공 탄성 모델들은 섬유를 완전히 무시하여 고체 골격을 균질한 네오-후크(neo-Hookean) 재료로 취급하거나, 크림프에 의한 점진적인 모집을 고려하지 않은 채 섬유를 포함하기도 했다. 본 연구는 비선형 다공 탄성 프레임워크 내에서 섬유 모집이 조직의 하중 및 언로딩(unloading) 반응, 특히 확산율, 응력-변형 비대칭성, 히스테리시스에 어떻게 영향을 미치는지에 대한 이해의 공백을 다룬다.

방법론
저자들은 근육과 뼈 사이에 위치한 건을 나타내는 1차원 비선형 다공 탄성 모델(FIB 모델)을 개발하였다. 고체 골격은 다음의 복합체로 모델링된다:

1. 크림프된 피브릴(Crimped Fibrils): 기질 내에 매립되어 초기에는 크림프 상태인 단단한 콜라겐 피브릴. 이들은 펴지기(uncrimped) 전까지는 하중을 받지 않는다. 모집은 피브릴이 점진적으로 펴짐에 따라 발생하는 피브릴 크림프 각도 $\theta(r)$의 분포에 의해 제어된다.
2. 비콜라겐성 기질(Non-Collagenous Matrix, NCM): 지속적으로 하중을 받는 부드러운 네오-후크 배경 기질.

모델은 불압축성 고체 입자와 유체를 가정하며, 유체 흐름은 Darcy의 법칙과 공극률 의존적 투과율(Kozeny-Carman 법칙)에 의해 지배된다. 지배 방정식은 공극률에 대한 비선형 확산 방정식이며, 여기서 확산 계수 $D$는 국소 강성($\partial s' / \partial \hat{\Phi}$)에 따라 달라진다.

FIB 모델은 등가 거시적 영률(Young's moduli)을 사용하는, 섬유 모집이 없는 순수 네오-후크 재료인 제어 모델(NH 모델)과 비교된다. 연구는 근육 인터페이스에서의 두 가지 경계 조건을 조사한다:

• 하중 인가(Loading): 이완된 건에 일정한 인장 응력을 즉각적으로 가함.
• 언로딩(Unloading): 사전 신장된 건으로부터 하중을 즉각적으로 제거함.

시뮬레이션은 공간에 대해서는 유한 체적법(finite volume method)을, 시간에 대해서는 암시적 가변 스텝 사이즈(implicit variable-step size) 방법을 사용하여 수행되었다.

주요 기여

• 다공 탄성성에 모집의 통합: 본 연구는 피브릴 크림프에 기반한 섬유 모집 규칙을 명시적으로 포함하도록 기존의 다공 탄성 프레임워크를 확장하여, 섬유 보강 조직과 균질한 네오-후크 조직 간의 직접적인 비교를 가능하게 하였다.
• 맞춤형 응력 법칙 유도: 저자들은 피브릴이 크림프 상태(하중 0)에서 팽팽한(linear load-bearing) 상태로 전환되는 과정을 고려하여, 피브릴과 NCM의 응력 기여를 구분하는 특정 응력 법칙을 유도하였다.
• 확산-강성 결합의 식별: 본 연구는 확산 계수가 상수가 아니라 섬유 모집에 따라 동적으로 변한다는 점을 강조한다. 피브릴이 모집됨에 따라 강성이 증가하여 확산율이 높아지는 반면, 피브리들이 다시 크림프되면 강성이 떨어져 확산율이 감소한다.

결과

• 하중 단계(Loading Phase):

  • 더 부드러운 반응: FIB 모델은 NH 모델보다 더 부드러운 응력-변형 반응을 보인다. 결과적으로 동일한 하중 하에서 FIB 건은 더 많이 늘어난다(높은 정상 상태 변형률/공극률).
  • 정상 상태 도달 속도 저하: 더 많이 늘어남에도 불구하고, FIB 모델은 NH 모델보다 정상 상태에 더 느리게 도달한다. 이는 초기 모집 단계가 낮은 강성(제로에 가까운 확산)을 수반하여, 조직 전체로 하중 반응의 전파를 늦추기 때문이다.
  • 공간적 불균질성: FIB 모델은 피브릴이 막 팽팽해지는 지점($Z_c$)에서 응력, 공극률 및 플럭스 프로파일의 급격한 굴곡점을 나타내며, 이는 NH 모델의 매끄러운 프로파일에서는 나타나지 않는 특징이다.

• 언로딩 단계(Unloading Phase):

  • 비대칭성 및 느린 이완: FIB 모델은 하중 인가와 언로딩 사이의 상당한 비대칭성을 보인다. 언로딩 시 조직이 이완됨에 따라 피브릴이 다시 크림프되며 조직이 부드러워진다. 이 연화 현상은 확산 계수를 감소시켜, 하중 반응의 전파를 극도로 늦추는 피드백 루프를 생성한다. 조직은 시뮬레이션 시간 내에 완전히 이완된 상태에 도달하지 못하는 경우가 많으며, 여전히 상당한 비율의 피브릴이 팽팽한 상태로 남아 있다.
  • 히스테리시스: 이러한 비대칭성은 평균 응력-변형 곡선에서 히스테리시스 루프로 나타난다. 히스테리시스는 언로딩 단계 초기에 피브릴이 조기에 재크림프되는 낮은 인가 하중에서 가장 두드드러진다. 높은 하중에서는 피브릴이 더 오래 팽팽하게 유지되어 초기 이완이 더 빠르고 히스테리시스가 감소하지만, 시스템은 여전히 주어진 시간 내에 완전히 이완되지 않는다.

• 보호 기제: 모델은 크림프된 피브릴의 존재가 건이 피브릴이나 NCM에 과도한 변형을 가하지 않으면서도 전체적인 신장을 더 크게 할 수 있게 하여, 구조적 구성 요소에 대한 자연스러운 보호 기제로 작용함을 보여준다.

의의 및 주장
본 논문은 크림프된 섬유 모집을 다공 탄성 프레임워크에 통합하는 것이 균질한 네오-후크 모델보다 건 역학에 대해 더 생물학적으로 현실적인 해석을 제공한다고 주장한다. 주요 의의는 히스테리시스와 하중/언로딩 비대칭성의 기원을 단순한 점탄성이 아니라, 유체 흐름의 확산적 특성과 강도 의존적 섬유 모집 사이의 상호작용을 통해 설명하는 데 있다.

저자들은 FIB 모델의 파라미터(크림프 각도, 피브릴 분율, 피브릴 탄성률)가 네오-후크 법칙에 필요한 현상론적 상수들보다 더 직접적으로 측정 가능하고 생물학적으로 해석 가능하다는 점을 언급하였다. 모델은 토(toe) 영역과 히스테리시스를 성공적으로 재현하였다. 그러나 저자들은 크림프된 피브릴이 하중을 전혀 받지 않는다는 가정(엄밀히 말하면 사실이 아닐 수 있음), 섬유 슬라이딩(fibre sliding)의 제외, 그리고 1차원 기하학적 단순화를 포함한 한계를 겸허히 인정하였다. 이 프레임워크는 향후 건 손상(예: 피브릴 파손) 및 고차원 모델링을 위한 기초로 제시된다.