[논문 리뷰] On the Application of Fractional Order Derivatives for Characterizing Brain White Matter Viscoelasticity

P. Pasupathy, J. G. Georgiadis외 1명|arXiv (Cornell University)|2026. 01. 26.

Automotive and Human Injury Biomechanics인용 수 0

한 줄 요약

이 논문은 축삭(axon) + ECM으로 구성된 뇌 백색질의 3D 부분적 점탄성 유한요소 모델을 스프링-포트 접근법으로 개발하고, VUMAT 구현을 검증하며, 축삭 체적분율의 함수로 균질화 매개변수를 추출합니다.

ABSTRACT

Conventional viscoelastic characterization of brain white matter (BWM), typically described using Prony series models, remains a largely empirical representation that is difficult to interpret physically. Growing evidence suggests that BWMviscoelasticity follows power-law behavior. Under the assumptions of linear viscoelasticity and causality, a power-law model in the frequency domain yields a fractional viscoelastic model in the time domain. A fractional viscoelastic constitutive model for the axon and extracellular matrix (ECM) is implemented via a Fortran VUMAT subroutine. A biphasic periodic finite element (FE) model of hexagonally packed representative volume elements (RVEs) of axons embedded in an ECM is constructed in Abaqus under quasi-static loading. The inverse problem of extracting homogenized material properties is solved using an optimization workflow. The model predicts that the springpot coefficient, which determines the solid-fluid behavior and, the power-law exponent, which encodes information about the underlying tissue architecture, follows a bi-logistic function along the transverse normal and shear directions. The nonlinear variation of the parameters reveals two distinct stiffening stages: a lower rate at low axon volume fractions, followed by a higher rate as increased axonal content reinforces the RVE. To our knowledge, this study is the first to propose and implement a 3D fractional viscoelastic FE model of the corpus callosum of BWM in the time domain. The thread-safe implementation of the VUMAT achieves significantly faster performance than existing approaches. The results reveal nonlinear variation in material parameters, directional dependence of BWM mechanics, and the complex interplay among microstructural elements.

연구 동기 및 목표

• Prony-시리즈 모델을 넘어 뇌 백색질의 물리적으로 해석 가능한 점탄성 특성화를 동기화합니다.
• 축삭과 ECM에 대한 멱법칙(분수) 스프링-포트 표현을 도입합니다.
• BVW 미구조를 시뮬레이션하기 위한 VUMAT 서브루틴을 가진 3D 유한요소 프레임워크를 개발합니다.
• 로딩 방향 전반에서 축삭 체적분율에 따른 균질화된 분수점탄성 매개변수를 보정합니다.

제안 방법

• 주파수 영역 데이터로 축삭과 ECM의 멱함수 파라미터를 로지스틱 회귀 방법으로 적합합니다.
• Riemann–Liouville 도함수(스프링-포트)를 이용한 분수점탄성 구성모형을 형식화하고 이를 3D FEM VUMAT에 구현합니다.
• 축삭이 ECM에 의해 육각적으로 포장된 RVEs를 주기경계조건으로 구성하고 6개의 로딩 방향을 적용합니다.
• FEM과 표적 응력 사이의 RMSD를 최소화하여 균질화된 파라미터를 회복하기 위한 최적화 워크플로(Nelder–Mead)를 사용합니다.
• Explicit 시간적 적분에서 과거 메모리 의존 계산을 줄이기 위한 단기 기억 잘라내기(truncated short-memory)를 도입합니다.

실험 결과

연구 질문

• RQ1 분수(멱법칙) 점탄성 모형이 Prony-시리즈 모델보다 BWM 미구조의 효과를 더 잘 포착할 수 있는가?
• RQ2 축삭 체적분율 및 로딩 방향에 따라 균질화된 스프링-포트 파라미터가 어떻게 달라지는가?
• RQ3 Abaqus에서 BVW 미구조에 대해 3D 분수점탄성 VUMAT를 구현하고 검증하는 것이 가능한가?
• RQ4 Explicit 시뮬레이션에서 메모리 절단(짧은 기억)으로 인한 정확도 및 계산 효율의 영향은 어떠한가?

주요 결과

• 스프링-포트 계수(E_beta)는 섬유 방향으로의 축삭 체적분율에 따라 선형적으로 변하고, 횡방향 및 전단 방향에서는 비선형적으로 변합니다.
• 멱함수 지수 beta는 섬유 방향으로는 포화하는 지수적 경향을 보이고, 횡방향 및 전단 방향에서는 이중 로지스틱 변화로 나타납니다.
• 균질화된 E_beta와 beta는 축삭 체적분율이 높아질수록 증가하며, 방향 의존성과 미세구조가 BWM 역학에 미치는 영향을 반영합니다.
• 두 개의 경직 증가 기전이 관찰되는데, 낮은 축삭 함량에서의 느린 증가와 축삭 함량이 증가함에 따라 더 빠르게 증가하는 경직 증가가 있습니다.
• 짧은 기억 접근법은 준정적(quasi-static) 케이스에서 정확도 손실이 허용 가능한 수준(약 2–3% 편차까지)으로 계산 시간을 줄입니다.
• 연구는 뇌의 백질 다발 부위인 교뇌의 3D 분수점탄성 FE 모델에서의 새로움을 주장하며, 기존 방법에 비해 thread-safe하고 더 빠른 구현을 제공합니다.