Search

A novel 3D micromechanical Finite Element Model (FEM) has been developed to depict the Poynting effect in bi-phasic Representative volume element (RVE) with axons embedded in surrounding extra-cellular matrix (ECM) for simulating the brain white matter response under simple and pure shear. In the proposed 3D FEM, nonlinear Ogden hyper-elastic material model describes axons and ECM materials. The modeled bi-phasic RVEs have axons tied to surrounding matrix. In this proof-of-concept (POC) FEM, three simple shear loading configurations and a pure shear scenario were simulated. Root mean square deviation (RMSD) were computed for stress and deformation response plots to depict role of axon-ECM orientations & loading condition on the Poynting effect. Variations in normal stresses (S11, S22, or S33) perpendicular to the shear plane emphasized role of fiber-matrix interactions. At high strains, the stress-strain% plots also indicated modest strain stiffening effects and bending stresses in purely sheared axons.

Traumatic axonal injury occurs when loads experienced on the tissue-scale are transferred to the individual axons. Mechanical characterization of axon deformation especially under dynamic loads however is extremely difficult owing to their viscoelastic properties. The viscoelastic characterization of axon properties that are based on interpretation of results from in-vivo brain Magnetic Resonance Elastography (MRE) are dependent on the specific frequencies used to generate shear waves with which measurements are made. In this study, we aim to develop a fractional viscoelastic model to characterize the time dependent behavior of the properties of the axons in a composite white matter (WM) model. The viscoelastic powerlaw behavior observed at the tissue level is assumed to exist across scales, from the continuum macroscopic level to that of the microstructural realm of the axons. The material parameters of the axons and glia are fitted to a springpot model. The 3D fractional viscoelastic springpot model is implemented within a finite element framework. The constitutive equations defining the fractional model are coded using a vectorized user defined material (VUMAT) subroutine in ABAQUS finite element software. Using this material characterization, representative volume elements (RVE) of axons embedded in glia with periodic boundary conditions are developed and subjected to a relaxation displacement boundary condition. The homogenized orthotropic fractional material properties of the axon-matrix system as a function of the volume fraction of axons in the ECM are extracted by solving the inverse problem., Comment: Accepted for publication at the 12th International Conference on Mathematical Modeling in Physical Sciences

A novel finite element method (FEM) is developed to study mechanical response of axons embedded in extra cellular matrix (ECM) when subjected to harmonic uniaxial stretch under purely non-affine kinematic boundary conditions. The proposed modeling approach combines hyper-elastic (such as Ogden model) and time/frequency domain viscoelastic constitutive models to evaluate the effect of parametrically varying oligodendrocyte-axon tethering under harmonic stretch at 50Hz. A hybrid hyper-viscoelastic material (HVE) model enabled the analysis of repeated uniaxial load on stress propagation and damage accumulation in white matter. In the proposed FEM, oligodendrocyte connections to axons are depicted via a spring-dashpot model. This tethering technique facilitates contact definition at various locations, parameterizes connection points and varies stiffness of connection hubs. Results from a home-grown FE submodel configuration of a single oligodendrocyte tethered to axons at various locations are presented. Root mean square deviation (RMSD) are computed between stress-strain plots to depict trends in mechanical response. Steady-state dynamic (SSD) simulations show stress relaxation in axons. Gradual axonal softening under repetitive loads is illustrated employing Prony series - HVE models. Representative von-Mises stress plots indicate that undulated axons experience bending stresses along their tortuous path, suggesting greater susceptibility to damage accumulation and fatigue failure due to repeated strains.