摘 要：为了更好地描述人工冻土的蠕变特性，将广义开尔文模型中的经典Newton黏壶中的黏滞系数定义为与时间有关的非定常参数，通过推算得到非定常开尔文模型，用粒子群优化算法识别模型中的参数。非定常开尔文模型能模拟不同应力下的蠕变试验数据，效果都很好，充分说明了该模型的合理性。目前采用非定常开尔文模型来研究人工冻土的蠕变规律的很少，此方法的提出为人工冻土领域的计算开辟了新思路。
关键词：人工冻土；蠕变特性；非定常开尔文模型
岩石力学理论与工程实践中的重要课题之一是本构模型的参数辨识。目前蠕变模型参数辨识方法主要有：回归方法、最小二乘法及其改进方法、模式搜索法以及蠕变曲线分解法等[3]。以上的优化算法都对问题有一定的约束条件，如要求优化函数可微等。粒子群优化算法作为一种仿生算法是一种模拟生物智能行为的优化算法，该算法利用生物群体内个体的合作和竞争等复杂行为产生群体智能，由于其几乎不存在对问题的约束，因此粒子群优化算法在各种优化问题中得到广泛应用[4]。粒子群优化算法是一种新兴的演化群体智能算法，其采用速度-位置搜索模型，每个例子代表解空间的一个候选解，解的优劣程度由适应度函数来决定[5-6]。运用模拟退火粒子群优化方法，对不同加载应力下蠕变试验数据建立的蠕变模型参数进行识别，得到蠕变模型的参数。
最后通过对比非定常蠕变模型计算值与试验结果，发现非定常开尔文蠕变模型能很好地模拟人工冻土的蠕变规律。非定常开尔文模型是人工冻土领域计算的一种新方法。
1 模型的建立
1.1 广义开尔文模型
广义开尔文模型由弹簧Ek和阻尼器并联，然后串联弹簧E1而成（见图1）。
非定常开尔文模型就是用与时间有关的函数替代广义开尔文模型中的黏滞系数η和串联的的弹性模量E1。
1.3 粒子群算法的描述
粒子群优化算法是基于群体智能理论的优化算法，通过群体中粒子间的合作与竞争产生的智能指导优化搜索[8]，为工程优化问题提供高效的解决方法。
在粒子群优化算法中，每个优化问题的潜在解都是搜索空间中的一只鸟，称之为“粒子”。所有的粒子都有一个由被优化的函数决定的适应值，每个粒子还有一个速度决定他们飞翔的方向和距离。粒子们就追随当前的最优粒子在解空间中搜索。优化开始时先初始化为一群随机粒子（随机解）。然后通过迭代找到最优解。在每一次迭代中，粒子通过跟踪两个极值来更新自己。第一个极值就是整个种群目前找到的最优解。这个极值是全局极值。另外也可以不用整个种群而只是用其中一部分作为粒子的邻居，那么在所有邻居中的极值就是局部极值。邻近种群数是在最优粒子附近的种群个数。第二个极值是粒子本身所找到的最优解，称为个体极值。这是因为粒子仅仅通过跟踪全局极值或者局部极值来更新位置，不可能总是获得较好的解。这样在优化过程中，粒子在追随全局极值或局部极值的同时追随个体极值则圆满的解决了这个问题。这就是粒子群优化算法的原理[9-10]。
基本粒子群优化算法可描述如下：
将V限定一个范围，使粒子每一维的运动速度都被限制在[-Vmin，Vmax]之间，以防止粒子运动速度过快而错过最优解，这里的Vmax根据实际问题来确定。当粒子的飞行速度足够小或达到预设的迭代步数时，算法停止迭代，输出结果[11]。
2 应用实例
2.1 人工冻土单轴蠕变试验
2.2 试验验证
采用粒子群算法来优化分数阶导数蠕变模型参数，非定常开尔文蠕变模型参数如表4所示。
3 结论
用非定常开尔文模型来模拟人工冻土的蠕变过程，不仅具有明确物理意义参数，而且能与试验结果很好地拟合。 1） 考虑到某些参数会受到时间的影响，以广义开尔文模型为基础，建立了非定常开尔文蠕变模型，其公式推导严密，模型具有明确的物理意义且模型参数少；
2）用非定常开尔文模型结合粒子群优化方法对不同温度下的蠕变规律进行了模拟，得到了相应的参数。
3）通过与试验值对比，发现非定常开尔文模型能很好的模拟蠕变曲线，相关系数基本都在0.995以上，为人工冻土领域的计算提供了一个新思路。本文只研究了一种土性的模拟，研究模拟结果比较好，是否能很好的模拟其他土性和是否能很好的模拟其他温度的蠕变模型将是下一步研究的方向。
