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

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.