On the Application of Fractional Order Derivatives for Characterizing Brain White Matter Viscoelasticity

TL;DR

This paper introduces a novel 3D fractional viscoelastic finite element model for brain white matter, revealing nonlinear, directional, and microstructural influences on tissue mechanics, with faster computational implementation.

Contribution

It is the first to develop and implement a 3D fractional viscoelastic FE model of brain white matter, demonstrating microstructural effects and improving computational efficiency.

Findings

Material parameters vary nonlinearly with axon volume fraction.

Directional dependence of viscoelastic properties is observed.

Model predicts two stiffening stages related to microstructure.

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.