摘要：本文回顾了常用的几种岩石粘弹性本构模型，依据蠕变柔量的概念，用位移和蠕变柔量两步反分析法从粘弹性本构模型的一般表达式中辨识出岩石的本构模型.先由相应的位移实测值用解析法反演岩石的蠕变柔量，再由岩石的蠕变柔量运用非线性优化技术辨识出岩石本构模型表达式，最后进一步由模型参数反算出岩体的粘弹性参数.文中最后给出了工程实例.
关键词：粘弹性本构模型 蠕变柔量 两步反分析法 反问题 岩石
岩体属于各向异性流变介质，由于其客观复杂性，在岩体理论分析和数值模拟方面，参数和模型的正确给定是岩石力学研究中的两大难题，它们的合理性将大大增强岩体工程分析设计的可靠性.新奥法施工及目前广泛应用的现场监控法或信息反馈施工法是把岩石的变形观测作为对原设计及施工方案进行修正的依据[1，2].这些方法紧密结合工程实际指导修改原设计，收到了良好的效果，已广泛应用于岩土工程实践领域.参数取值，传统的方法大多凭经验、工程类比来进行，常常不能作出科学的判断.参数反演作为参数辨识的一种方法，是基于实测位移反求系统某种参数的一种逆问题，对确定认识系统和进一步正演分析相当重要.参数的辨识，包括模型参数的辨识，是在模型的结构式给定的情况下进行的，模型的给定与实际相符合的程度显然相当重要.工程岩体是复杂的不确定系统，岩体的流变力学模型较多，模型识别的问题有着重要的理论意义和实用价值.在岩体本构模型的辨识，尤其是与时间有关的粘弹性本构关系辨识方面当前还有许多问题值得研究.
本文回顾了工程应用和科学研究中几种常用的岩石粘弹性本构模型，然后依据蠕变柔量概念，用位移和蠕变柔量两步反分析法从粘弹性本构模型的一般表达式中辨识岩石本构模型，第一步由相应的位移实测值用解析法反演岩石的蠕变柔量，第二步由岩石的蠕变柔量运用非线性优化技术辨识出岩石本构模型的表达式，最后进一步由模型参数与粘弹性参数的关系反算出岩石的粘弹性参数.
1 粘弹性本构模型回顾
工程岩体的长期稳定性是当今岩土工程领域的一个十分重大的前沿问题[3].工程岩体的变形破坏是有时间过程的，具有时间效应，岩体流变是岩体的重要基本力学特性和行为，对于某些实际岩体工程，研究变形的时间效应有着极其重要的意义[8,9].真实岩体是一种非连续、非均质、各向异性的流变介质，目前很难用数学手段作出适当精确描述.在一定情况下，可将岩体视为似连续、均匀化、类各向同性介质，采用经验与理论相结合，定量分析定性使用的原则，对工程设计与施工起重要的定性指导作用[10,11].在这一前提下，在岩土工程应用和科学研究中提出的粘弹性本构模型主要包括经验模型和组合模型[5](弹性和粘性元件的不同组合而成).表1列出了常用的几种微分型组合模型的公式及其特征(h代表弹性元件，n代表粘性元件).在一般状态下，这些组合模型的本构方程微分形式的一般式可表示为[5]
表1 常用岩石粘弹性微分型组合本构模型
式中：
式(1)即线粘弹性微分型本构关系的一般表达式.
2 粘弹性模型识别
2.1岩体介质蠕变柔量的确定在弹性岩体中开挖任意形状断面的隧道，应用平面复变影射、保角变换方法可得围岩内任一点在曲线坐标中的位移为[7]
由弹性理论确定.其中：ω(ζ)为垂直隧道轴线的物理平面z上的非圆形洞室外域到数学平面ζ内单位圆外域的映射函数；φ(ζ)和ψ(ζ)为满足应力边界条件的复势函数.
应用弹性—粘弹性对应原理，对照式(2)，得出任意断面的隧洞内任一点因开挖引起的粘弹性位移在曲线坐标中的表达式为[8]
d(x)为荷载逐步释放系数，若计算断面远离开挖面，则可作为平面应变问题处理；若接近开挖面，则严格说来是空间问题.为使问题简化，考虑了开挖面空间效应后，仍作为平面问题来处理.开挖面空间效应通过洞周释放荷载逐步释放代替瞬间完全释放来体现.据研究资料[4]表明，其取如下形式：d(x)=1-0.7exp(-3.15x/2a).其中，a为隧洞半径；p0为垂直向初始地应力；g为剪切模量；r为所考虑点到隧道轴线的距离；λ为侧压力系数；k为体积模量，d(x)为应力释放系数，x为计算断面距开挖面距离.
在流变岩体中开挖任意形状的隧道，在开挖之前布置量测仪器，可测到由于开挖引起的全部位移，在ti时刻可测到n个测点由于开挖引起的相对径向全部位移，记为uk(ti),k=1,2,，…，n,i=1,2,,…，l(假定量测了l个时段)，代入式(4)，则有蠕变柔量jl(ti)和广义蠕变柔量j2(ti)的线性方程组
2.2粘弹性应力应变关系的确定线粘弹性微分型本构关系的一般表达式如式(1)所示，对其进行拉氏变换，考虑光滑化假定，得[7]
对于
(s)=0有重根或虚根的情况，j(t)中含有时间t的阶次项和正余弦项，且未知参量个数增加，给反演优化大大增加了难度，收敛性差，不利于程序的统一编制，这里作了一些简化，在工程实际中，复杂的高阶次模型也不必要.
式(10)中，在模型阶次确定时，参数ai(i=1～2m+1)与模型式(1)中的模型参数p0,p1,p2,…,pm,q0,q1,q2,…,qn有一定的对应关系，由式(9)和式(10)可推导得出，模型参数均可用ai简单运算表达.
表2 常用岩石粘弹性模型参数表
利用实测位移由式(6)可求出众多不同时段ti的蠕变柔量ji(ti)的最佳估计值，又由式(10)可得相对应的j(ti)(i=1～l)，因而构造如下的非线性优化目标函数：
（11）
式(10)是ai的非线性函数，该问题是一非线性最小二乘优化问题，带约束的隐式非线性优化问题，这里采用了可变容差优化方法，它是在单纯形法和复合形法的基础上变化而来，具体见文献[6].
上述模型参数估计是在模型已经确定的情况下进行的，即偏微分方程的阶次m,n为已知的.这里模型的最佳阶次m,n的确定不是通过数学推导求出，而是通过试验来确定的.令模型阶次m,n的取值从1开始，分别以步长1递增，比较m,n取不同值时最优估计式(11)的值，即拟合误差.通常随着阶次增大，拟合误差下降.但当阶次为最佳阶次时，拟合误差达到最小；然后随着m,n的增大，拟合误差趋向增大.
考虑河岸冲刷的弯曲河道水流及河床变形的数值模拟
黄土高原沟壑坝系建设有关问题探讨