盐岩流变本构的分数阶导数描述及其参数拟合分析.pdfVIP

盐岩流变本构的分数阶导数描述及参数拟合分析 王春萍1，周宏伟1，韩冰冰1，段志强1，侯正猛2 (1中国矿业大学(北京)煤炭资源与安全开采国家重点实验室，北京100083) (2德国Clausthai工业大学) 摘要：从分数阶导数出发．讨论了一种新的流变元件即Abcl粘壹．并通过AbeI粘壶代替 经典西原模型中Newton粘壶的方法，构建了基于分数阶导数的盐岩流变本构模型，给出了 流变本构模型的求解过程．论文进一步根据盐岩的流变实验数据，对分数阶导数盐岩流变 本构模型的参数进行了拟合分析，确定了模型中的参数，通过比较发现，分数阶流变本构 模型可更好地与实验数据相吻合． 关键词：盐岩；流变模型；分数阶导数；西原模型：参数拟合 一、引言 盐岩以其低孔隙率、低渗透性、损伤自恢复等良好性能而成为国际上公认的地下战 略能源储备、CO，封存、高放废物处置的理想介质，对盐岩流变特性的研究是一个十分 and 重要的课题。分数阶微积分在力学中的应用最早始于材料粘弹性本构的研究(Bagley 数出发，讨论了一种新的流变元件即Abel粘壶，并构建了基于分数阶导数的盐岩流变 本构模型(分数阶硒元模型)。 二、基于分数阶导数的盐岩流变本构模型 2．1Riemann．Liouville分数阶导数及流变本构模型的建立与求解 1 一 。 y阶分数阶积分闭：z)-7／(f)=熹【(t-善)v-z八f)蟛 l LyJ” Y阶分数阶导数：∥巾)=ZY[D-”叭，))】，m是大于厂的最小整数，且y=朋·yO- Abel粘壶的本构关系为：or(t)=r／7D716(0I(o≤y≤1)。 本文用Abel粘壶取代经典西原正夫模型中的Newton粘壶，构造了基于分数阶导数 的粘弹塑性模型，并得其本构方程为： 、 “ d力：云+吾喜页霎嘉乏芸晶何巧，，￡c。：薏+吾砉 +Tcr-正丽g何≥吒)。瞄r(1+力 2．2分数阶导数流变本构模型的参数拟合分析 5 本文根据德国Clausthal工业大学进行的单轴盐岩流变实验的实验数据，通过拟合分 析来确定口≥吒时分数阶微分流变模型的参数‘31。根据非线性最小二乘的方法，参数拟 合的结果由表l给出，实验数据及模型拟合所得的蠕变曲线由图l给出。 表l参数拟合的结果 毛 巨 彳 孵 厂 误差 Gpa GpaOpa·h7Opa·h7 分数阶模 22．989．88 1．23 13．76 0．3I4．3lxlo．． 型 M I仃 图1实验数据及分数阶流变本构模型 三、 结论 本文通过将西原模型中Newton粘壶替换为分数阶Abel粘壶的方法，构建了基于分 数阶导数的新的盐岩流变本构模型。经过分析我们发现Abel粘壶，与阻尼器元件比较， 能更好地反应流变问题的应力或应变的非线性渐变过程。分数阶流变模型中当导数阶次 等于1时，分数阶导数描述的流变模型即可退化为经典西原模型。另外，由本文的分数 阶流变本构模型拟合所得的蠕变曲线与实验数据的吻合效果还是不错的。 参考文献 l P．J。Fractionalcalculus-Adifferent tothe of