Budday S, Sommer G, Holzapfel G A, Steinmann P, Kuhl E
Department of Mechanical Engineering, University of Erlangen-Nuremberg, 91058 Erlangen, Germany.
Institute of Biomechanics, Graz University of Technology, 8010 Graz, Austria.
J Mech Behav Biomed Mater. 2017 Oct;74:463-476. doi: 10.1016/j.jmbbm.2017.07.014. Epub 2017 Jul 11.
Understanding the constitutive behavior of the human brain is critical to interpret the physical environment during neurodevelopment, neurosurgery, and neurodegeneration. A wide variety of constitutive models has been proposed to characterize the brain at different temporal and spatial scales. Yet, their model parameters are typically calibrated with a single loading mode and fail to predict the behavior under arbitrary loading conditions. Here we used a finite viscoelastic Ogden model with six material parameters-an elastic stiffness, two viscoelastic stiffnesses, a nonlinearity parameter, and two viscous time constants-to model the characteristic nonlinearity, conditioning, hysteresis and tension-compression asymmetry of the human brain. We calibrated the model under shear, shear relaxation, compression, compression relaxation, and tension for four different regions of the human brain, the cortex, basal ganglia, corona radiata, and corpus callosum. Strikingly, unconditioned gray matter with 0.36kPa and white matter with 0.35kPa were equally stiff, whereas conditioned gray matter with 0.52kPa was three times stiffer than white matter with 0.18kPa. While both unconditioned viscous time constants were larger in gray than in white matter, both conditioned constants were smaller. These rheological differences suggest a different porosity between both tissues and explain-at least in part-the ongoing controversy between reported stiffness differences in gray and white matter. Our unconditioned and conditioned parameter sets are readily available for finite element simulations with commercial software packages that feature Ogden type models at finite deformations. As such, our results have direct implications on improving the accuracy of human brain simulations in health and disease.
了解人类大脑的本构行为对于解释神经发育、神经外科手术和神经退行性变过程中的物理环境至关重要。人们已经提出了各种各样的本构模型来描述不同时空尺度下的大脑。然而,它们的模型参数通常是在单一加载模式下校准的,无法预测任意加载条件下的行为。在这里,我们使用了一个具有六个材料参数的有限粘弹性奥格登模型——一个弹性刚度、两个粘弹性刚度、一个非线性参数和两个粘性时间常数——来模拟人类大脑的特征非线性、调节、滞后和拉伸-压缩不对称性。我们在剪切、剪切松弛、压缩、压缩松弛和拉伸条件下,对人类大脑的四个不同区域——皮质、基底神经节、放射冠和胼胝体——进行了模型校准。令人惊讶的是,无调节的灰质(0.36kPa)和白质(0.35kPa)同样坚硬,而有调节的灰质(0.52kPa)比白质(0.18kPa)硬三倍。虽然无调节的粘性时间常数在灰质中都比白质中更大,但有调节的常数都更小。这些流变学差异表明两种组织之间存在不同的孔隙率,并至少部分解释了关于灰质和白质刚度差异报道之间持续存在的争议。我们的无调节和有调节参数集可直接用于使用具有有限变形奥格登类型模型的商业软件包进行有限元模拟。因此,我们的结果对提高健康和疾病状态下人类大脑模拟的准确性具有直接影响。