循环荷载下饱与软黏土累积变形遗传分数阶导数开尔文模型.doc

循环荷载下饱和软黏土累积变形遗传分数阶导数开尔文模型 　　摘要：研究循环荷载作用下饱和软黏土的累积变形对计算交通荷载引起的地基沉降具有重要意义。对饱和软黏土分别进行围压为100kPa、150kPa、200kPa的等向固结动应力比为0.1、0.2和偏压固结动应力比为0.06、0.1、0.2的循环加载试验，得到轴向循环塑性应变与时间的关系，分析得到围压、动应力比及固结形式对轴向循环塑性应变的影响规律；将Abel黏壶引入开尔文模型，得到分数阶开尔文模型；利用遗传算法优化轴向循环塑性累积应变的开尔文模型和分数阶开尔文模型的参数，通过分析两组模型的计算值与试验值的对比曲线，得到分数阶开尔文模型更适合模拟计算循环荷再下饱和软黏土的累积变形 关键词：饱和软黏土；开尔文模型；分数阶开尔文模型；遗传算法 中图分类号：TU473 文献标志码：A 文章编号：1672-1098（2016）04-0052-08 在路基沉降计算中，长期交通荷载引起的软土地基累积沉降占主要部分，因此研究循环荷载作用下饱和软黏土的累积变形具有非常重要的意义 国内外学者通过大量的室内试验和现场试验，经拟合分析得到描述不同循环累积应变和循环次数之间关系的经验模型。文献[1]在考虑累积塑性应变速率和循环加载次数关系、循环加载动应力大小及应力历史等的影响因素条件下提出了指数模型。文献[2]提出了主要考虑应变与循环次数的指数关系的指数模型，而在此后提出的指数模型多是在此基础上提出的。如文献[3-4]提出的考虑初始静偏应力、静破坏偏应力和动偏应力影响的模型，文献[5-7]等提出的相对偏应力模型。这些由试验拟合得到的经验模型参数没有明确的物理意义，且难以确定。而由整数阶元件组合而成的整数阶模型具有理论性强，物理概念清晰的特点。但从以往的试验和研究得到，在运用整数阶微积分得到的整数阶微积分型本构方程描述岩土的蠕变曲线时，在蠕变或松弛的初始阶段不能与实验数据很好吻合。在整数阶微积分模型基础上得到的分数阶微积分模型不仅能体现整数阶微积分模型的优点，而且能改善其不足之处 本文是将Abel黏壶引入开尔文模型，得到分-数阶开尔文模型，用遗传算法优化参数来计算循环三轴轴向循环塑性累积应变。遗传算法是采用二进制遗传编码的一种自适应全局优化搜索算法。通过分析广义开尔文模型和分数阶开尔文模型的计算值与试验值，得到分数阶开尔文模型更能真实的反映循环荷载下饱和软黏土的累积应变特征 1.循环荷载下饱和软黏土三轴试验 1.1试样制备及方法 试验所用的土样深度为地表以下10-12m的上海第④层淤泥质软黏土。为了尽可能减少对土样的扰动，用挖土机挖至取样深度，将挖土面修平后，把取样筒垂直插入土中，然后把周围土挖掉取出土样 试样尺寸为39.1mm（直径）×79mm（高度），其物理指标如表1所示。试验所用的仪器是英国GDS多功能三轴仪。分别进行等向固结条件下动应力比为0.1、0.2和偏压固结条件下动应力比为O.06、0.1、0.2的循环加载试验，试验所采取的围压σ分别为100kPa、150kPa、200kPa。动应力比ηd。定义为ηd=qd/σ3c，其中qd为动应力峰值。试验采用的加载波形为半正弦波形，加载频率是1Hz，加载时间均为6000s 1.2试验结果及分析 偏压固结条件下循环三轴轴向循环塑性累积应变与循环时间的关系如图1所示，等向固结条件下循环三轴轴向循环塑性累积应变与循环时间的关系如图2所示 从图1可以得出，随着循环时间的增加，轴向循环塑性累积应变逐渐增加，最后逐渐趋于一个稳定值。对比图1中的曲线得到，当围压相同时，动应力比大的，对应的达到的稳定累积应变值也大，可见动应力比对轴向循环塑性累积应变影响显著 从图2中得出在等向固结条件下，轴向累积应变也是随循环时间的增加逐渐增加最后趋向一个稳定值。但对比图1、图2可以得到，当围压和动应力比相同时，等向固结条件下的轴向累积应变要大于偏压固结条件下的轴向累积应变，可见固结形式对轴向循环累积应变的影响也较显著 2.分数阶导数模型 2.1广义Kelvin模型 广义Kelvin模型元件是由一个开尔文元件和一个牛顿体串联而成（见图3） 2.2分数阶导数开尔文模型 采用分数阶导数的流变模型理论实质上就是用Abel黏壶取代经典模型理论中的Newton黏壶，分数阶导数开尔文模型组成如图4所示 2.3遗传算法优化模型参数 遗传算法（Genetic Algorithm，GA）是一种基于自然群体遗传演化原理而建立的一种最优化的高效探索算法，它是美国Michigan大学Holland教授于1975年首先提出来的。遗传算法是用于解决最佳化的搜索算法，是进化算法的一种。在遗传算法中，它将问题