[論文レビュー] On the Application of Fractional Order Derivatives for Characterizing Brain White Matter Viscoelasticity

P. Pasupathy, J. G. Georgiadis+1 more|arXiv (Cornell University)|Jan 26, 2026

Automotive and Human Injury Biomechanics被引用数 0

ひとこと要約

この論文は、軸索と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-seriesモデルを超える物理的に解釈可能な脳白質の粘弾性特性の動機づけ。
• 軸索とECMのためのべき-law（分数）スプリング-ポット表現の導入。
• RVEのビパシック挙動を模擬するVUMATサブルーチンを備えた3D有限要素フレームワークの開発。
• 軸索体積分率に対する水平・垂直・せん断方向の均質化分数粘弾性パラメータを較正。

提案手法

• 周波数領域データからロジスティック回帰を用いて軸索とECMのべき法パラメータを適合。
• リーマン-リウラ derivative（スプリング-ポット）を用いた分数粘弾性構成モデルを定式化し、3D FEMのVUMATに実装。
• 周期境界条件付きの軸索-ECMの六方向荷重を適用した六角形充填RVEsを構築。
• FEMとターゲット応力間の RMSDを最小化して均質化パラメータを回復する最適化ワークフロー（Nelder–Mead）を適用。
• 明示的時間積分における履歴依存計算を減らすための短期記憶切断を組み込む。

実験結果

リサーチクエスチョン

• RQ1 分数（べき-law）粘弾性モデルはProny-seriesモデルよりBWMマイクロ構造の影響をうまく捉えられるか。
• RQ2 均質化スプリング-ポットパラメータは軸索体積分率と荷重方向でどのように変化するか。
• RQ3 AbaqusでBVWマイクロ構造の3D分数粘弾性VUMATを実装・検証することは可能か。
• RQ4 明示的シミュレーションにおけるメモリ切断（ショートメモリ）の精度と計算効率への影響はどの程度か。

主な発見

• スプリング-ポット係数(E_beta)は繊維方向に沿って軸索体積分率と直線的に変化し、横断・せん断方向では非線形。
• べき法指数betaは繊維方向に沿って飽和する指数関数的傾向を示し、横断・せん断方向では二重対数ロジスティック変化。
• 均質化したE_betaとbetaは軸索体積分率の増加に伴い増加し、方向依存性とマイクロ構造の影響を反映。
• 低軸索含有量での低速剛化と、軸索含有量が増加するにつれての高速度剛化という二つの硬化レジームが観察される。
• ショートメモリアプローチは、準静的ケースで計算時間を削減しつつ精度は許容範囲内（約2–3%の偏差まで）で維持。
• 本研究は大脳梁の3D分数粘弾性FEモデルの新規性を主張し、既存手法と比較して安全で高速な実装を提供。