Demirci Nagehan, Tönük Ergin
Department of Mechanical Engineering, Middle East Technical University, Ankara, Turkey.
Department of Mechanical Engineering, Graduate Program of Biomedical Engineering, Biomaten, Center of Exellence in Biomaterials and Tissue Engineering, Middle East Technical University, 06800 Ankara, Turkey.
Acta Bioeng Biomech. 2014;16(4):13-21.
During the last decades, derivatives and integrals of non-integer orders are being more commonly used for the description of constitutive behavior of various viscoelastic materials including soft biological tissues. Compared to integer order constitutive relations, non-integer order viscoelastic material models of soft biological tissues are capable of capturing a wider range of viscoelastic behavior obtained from experiments. Although integer order models may yield comparably accurate results, non-integer order material models have less number of parameters to be identified in addition to description of an intermediate material that can monotonically and continuously be adjusted in between an ideal elastic solid and an ideal viscous fluid.
In this work, starting with some preliminaries on non-integer (fractional) calculus, the "spring-pot", (intermediate mechanical element between a solid and a fluid), non-integer order three element (Zener) solid model, finally a user-defined large strain non-integer order viscoelastic constitutive model was constructed to be used in finite element simulations. Using the constitutive equation developed, by utilizing inverse finite element method and in vivo indentation experiments, soft tissue material identification was performed.
The results indicate that material coefficients obtained from relaxation experiments, when optimized with creep experimental data could simulate relaxation, creep and cyclic loading and unloading experiments accurately.
Non-integer calculus viscoelastic constitutive models, having physical interpretation and modeling experimental data accurately is a good alternative to classical phenomenological viscoelastic constitutive equations.
在过去几十年中,非整数阶导数和积分越来越多地用于描述包括软生物组织在内的各种粘弹性材料的本构行为。与整数阶本构关系相比,软生物组织的非整数阶粘弹性材料模型能够捕捉从实验中获得的更广泛的粘弹性行为。尽管整数阶模型可能会产生相当准确的结果,但非整数阶材料模型除了描述一种可以在理想弹性固体和理想粘性流体之间单调连续调整的中间材料外,需要识别的参数数量更少。
在这项工作中,从关于非整数(分数)微积分的一些预备知识开始,构建了“弹簧罐”(固体和流体之间的中间力学元件)、非整数阶三元件(齐纳)固体模型,最后构建了一个用户定义的大应变非整数阶粘弹性本构模型,用于有限元模拟。利用所建立的本构方程,通过逆有限元方法和体内压痕实验进行软组织材料识别。
结果表明,从松弛实验中获得的材料系数,当用蠕变实验数据进行优化时,可以准确模拟松弛、蠕变和循环加载卸载实验。
具有物理解释且能准确模拟实验数据的非整数微积分粘弹性本构模型是经典唯象粘弹性本构方程的一个很好的替代方案。