所谓岩土工程位移反分析资料

所谓岩土工程位移反分析，即以现场测量到的位移为基础，通过数学物理反分析模型，得到
岩土介质的本构模型及等效力学参数（如初始地应力、变形参数、强度参数等）的方法。

最终目的是建立一个输出位移更接近现场实测位移的理论模型，以便较正确地反映或预测岩土
结构的某些力学行为。

20世纪70年代初人们开始岩土工程位移反分析的研究，随着岩土工程的发展，国内外众多学者对位移反分析的理论与应用进行了大量广泛而深入的研究。

岩土工程位移反分析
涉及的研究内容非常广泛，下面就从位移解析解、位移反分析的唯一性、位移测量点的优化布置、本构模型、数值计算方法、优化方法这六个方面对其进行综合地考察。

1.3.1 位移解析解
1898年，Kirsch［92］最早发表了弹性平板中圆孔周围的二维应力分布解，Jaeger和Cook （1969）［93］对Kirsch 方程进行了详细的推导。

此后，Poulos 和Davi（s1974）［94］、Pender
（1980）［95］、Carter（1982）［96］和Verruijt（1999）［97］分别在不同的边界条件下给出了圆形巷道的位移解析解。

Exadaktylos（2002）给出了半圆形巷道的位移解析解［98］。

Muskhelishvili （1953）［99］和蔡晓鸿（2008）［100］分别在不同的边界条件下给出了椭圆形巷道的位移解析解。

吕爱钟（1998）［10］、张路青（2001）［101］求解了不同地应力条件下任意形状巷道的位移解析解。

1.3.2 位移反分析的唯一性
反分析的唯一性是位移反分析中最重要却研究得最不充分的理论问题之一。

迄今为止，国外尚未有相关论文发表，国内的论文也是凤毛麟角。

吕爱钟（1988）［103］推导了参数可
辨识条件，论证了地下洞室弹性位移反分析的多种唯一性问题，并指出某些问题无论安装多
少个位移测点其反分析的结果都不是唯一的。

张路青（2001）［101］进一步研究了考虑剪应
力时位移反分析的唯一性问题。

杨志法［104］则用几何作图法证明了图谱反分析的唯一性。

以上文献中均假设岩体为各向同性材料，国内外尚未有报道对各向异性岩体位移反分析的唯
一性问题进行研究。

为了更大限度地利用位移反分析方法及更有效地指导施工与决策，本文
将在以上研究的基础上，对横观各向同性岩体位移反分析的唯一性问题进行研究。

1.3.3 位移测量点的优化布置
测量点的布置对位移反分析的唯一性和反分析精度都有很大的影响，因此，众多研究人
员对位移测量点的优化布置进行了研究。

测量点的布置包括测量点的数目和测量点的空间
位置两方面的工作。

由于工程费用和测量现场条件的制约，测量点的数目应该控制在一定
范围之内。

关于测量点数目下限的确定，大家公认的原则是测量点数目至少要大于等于待反
演参数的数目，否则就会因为信息量不足而导致反分析失败。

但是在满足下限的情况下，测
量点数目越多越好，还是越少越好，目前并没有形成统一的认识。

Kemevez（1978）［107］发现并不是测点数目越多越好，而是测量点的空间位置更加重要。

Cividini （1981）［109］的大量计算结果表明多测点的反分析结果不一定比少测点的反分析结果好，应该综合考虑测点
的数量和空间位置。

沈新普（1995）［105］认为测量误差会导致反分析结果偏离真值较远，应该尽可能多地布置测点，从而消除测量误差的影响。

孙钧（1996）［108］从工程计算中观察到太多的测点并不能显著地改进位移反分析结果。

Jim（2000）［106］认为较多的测点数目
能够提高反分析的效果，但是随着测点数目增加得越多，提高的效果就越不明显。

对测量点的空间位置应该遵循的原则，目前主要有以下几种［174］:
最大位移原则［110］。

该原则认为位移绝对值比较大的测点，测量的相对误差就比较小，测量精度就比较高，所得数据的实用价值就越大。

不过，由于现代测量技术的迅猛发展，对量值比较小的位移，测量精度已经有了较大的提高。

最大灵敏度原则［111, 112］。

灵敏度反映了位移测量值相对于待反演参数的变化。

该原
则考虑到灵敏度越大就越有利于参数的反演，所以就以灵敏度最大作为测点布置的原则。

但是，当待反演的参数不仅一个时，依照这种原则很难得到一个适用的综合判断指标。

最小方差原则［113, 114］。

按照该原则布置测点时遵循的原则就是使参数估计误差的方
差最小。

这与Fisher信息矩阵［115］有关，许多学者进行了大量测量点布置的研究［114-117］。

在实际工程中，提出的大部分测点布置原则都与D矩阵的性质有关。

主要的测点布置原则
包括：① D最优：使D矩阵的行列式值最小化；② A最优：使D矩阵主对角元素的和最小化；③E 最优：使D矩阵的最小特征值最大化。

对比三个原则表明［118］，A最优原
则倾向选择灵敏度最大的测点，E最优原则倾向选择灵敏度最小的测点，而按照D最优原
则选择的测点的灵敏度介于A最优和E最优两者之间。

以上原则分别从不同的侧面对测点
的优化布置提供了很有价值的参考。

但是以上的研究大多是针对各向同性岩体，并且到目前
位置，关于位移反分析的测点布置准则，并没有得到统一的认识。

本文试图基于最大位移原
则对横观各向同性岩体位移反分析中的测点优化布置问题进行初步研究，探寻合适的位移测
点布置准则。

1.3.4 本构模型
纵观岩体位移反分析的发展过程，反分析中使用的有代表性的本构模型如表 1.1~表1.5
所示。

需要指出的是，在岩土工程的反分析中，还有另外一种选取物理模型的思路，即放弃传统的各种力学本构关系，而将岩土体作为一个黑箱系统，直接从该系统的激励和响应的
样本数据中总结提炼规律。

其中最具代表性的是基于人工神经网络的智能方法，这种方法具
有很强的自学习能力和对环境的适应能力，因此它特别适用于模拟岩土体这类复杂的系统。

目前已有很多学者在这方面取得了一些成果［119-122］，但由于供给神经网络进行学习的大量
样本不容易获得，从而影响了这种模型的实际应用。

农1.1反分折中使用的并向同性弹性模型
研究者时间研究内無
Kirsten^1976辨识岩体的弹性模量
Kovari1^1977捉岀反舁地层压力参数的方仏
Sakurai^ ]26J1979. 1983捉出了卩面应变弹性问题位移反分析.硏究了隧洞嗣
岩地应力及岩休弹性模量的逆反分析法Goida1^1980辨识君体的蝉性力学参数
利用爭先建a的图谱反分析乘直r洞轴的小坍地应力畅志法吟闻1981分录及弹性犠就，在国内首次将该法成功地用于甘肃
'TI大人秦”工程的试验洞
Arai[l3W33]1983. 1984同时辨识弹性参数和灌透系数：指出应力和净水压现测虽无助于弹性参数的反分析
畅林德凹1985用冇限沁法辨识地应力和弹性模吐
1986垮虑隧道衬砌的二维弹性反分析•才人彳十川玮.
口
1998. 2000研究动态施匚过程屮的弹H模址反分析石安池32006将模拟担火BP网络用丁边坡弹杵模量反分析
表1・3反分析中使用的备向同性粘弹性模型
研究者时间研究内容
Hisatakc114711985辨识地应力和粘弹性参数
王芝银X8W]1988. 1993辨识地应力和粘弹性参数
李云鹏^*5211988. 1989用竹限元和边界元法辨识粘弹性参数杨林徳皿1991用边界元法辨识地应力和粘弹性参数
陈国荣直⑧1996、1998利用阻尼域小二乘法进行三元件模型的粘弹件反分析Ohkamt^15711997、1998Laplace变换及对偶边界控制法辨识粘弹性参数高强晦2000应用遗传算法进彳f粘弹性位移反分析
王永岩晦2007粘弹性有限元反分析在软岩流变问题中的应用
表1.2反分析中使用的各向同性弹型性模熨
研究者时间研究内容
Goida皿1980辨识岩体的粘结力和内原擦角
Sakurai^,4,)1985. 1994辨识岩体的塑性区和各种教性参数
在应力空IHJ及应变空间中川边界元法进行弹塑性位移郑颍人凹1987反分析，根据围岩位移反算初始地应力和弹性模量，可用
于解决二维和三维问题
Anandarajah^11991用变尺度法反演粘土塑性参数，用半解析法计算梯度孙钧皿1992提出多◎次参数反分析的概念以寻找全局最优解
尹蓉蓉国2004对采用莫尔■库仑准则、徳曾克-普拉格准则两种弹塑性反分
析的结果进行了对比
李宁【146]2006围岩松动圈的弹塑性位移反分析方法探索表1. 4反分析中使用的各向同性粘呦性模型
研究者时间研究内容
Cividini1國1985辨识岩土的粘塑性参数
王芝银皿问1985、1990辨识地应力和粘塑性参数
沈振中血1997利用可变容差法进行三维粘弹测性位移反分析Dziadziuszko1-^-12000研究了粘塑性模型的参数辨识问题赵新铭皿2008岩体弹粘塑性参数的Bayes随机优化反分析
1.5反分析川使用的再向异性模型
研究者时间研究内容
Ohkami11^11990. 1993用调节参数的方法來尚动辨识岩体的力学模型1991用边界元法辨识平面正父齐向异件材料参数
李素华皿1993横观各向同性互层岩体的优化反分析
Dm严网IQ97辨识横观齐向冋件材料参数和节理单元参数
刘学增皿2002研究动态施1,过程的横观各向同性粘弹性反分析
黄立新皿2005研究平面川交齐向并性体参数识别的边界汇仏
邙宏柱皿2006壊用遗传算法进押橫观各1-0同性粘弹性反分析
黄光明皿2006进行碾压混凝十坝横观齐向同性粘弹性参数的反分析
1.3.5 数值计算方法
正分析的正确解决是反分析的基础，而合适的数值计算方法可以改进正分析的速度、精度和通用性，所以在反分析的正计算中要采用一种合适的数值计算方法。

近年来，伴随着计算机技术的飞速发展，岩体稳定性分析的数值计算方法日臻成熟。

当前应用于岩体工程问
题的主要数值分析方法有：有限单元法、边界元法、有限差分法、离散单元法、无限元法、界面单元法、无单元法、非连续变形分析、流形元法以及由以上各种方法相组合而得到的混合数值计算方法。

当岩体被裂隙切割成块体集合时，非连续的数值方法如离散单元法、非连续变形分析等可以更逼真地反映岩体的内部结构，但块体的拓扑分析过于繁杂，所以目前在岩土工程的数值计算中，应用较广的还是基于连续介质力学的数值计算方法。

纵观岩体位移反分析方法的发展过程，数值计算方法大多选用有限单元法和边界元法。

这是因为，有限单元法是岩石力学数值计算方法中最为广泛应用的一种。

自20世纪50年代发展至今，有
限元已成功求解了许多复杂的岩石力学与工程问题。

有限元法的突出优点是适于处理非线性、非均质和复杂边界等问题，而岩体应力变形分析就恰恰存在这些困难问题。

边界元法在20世纪70年代得到迅速发展，有限元法是全区域离散化，而边界元法仅对边界离散化。

这样使三维问题降为二维问题求解，使二维问题降为一维问题求解，当物体的表面积和体积
之比比较小时，边界元的划分单元数要比有限元少数倍或十几倍，这样也使待解的方程数目、
处理和存储的数据量降低同样的倍数，大大节省了机时。

边界元法比较适合求解无限区域和
半无限区域问题，如深埋巷道是一个典型的例子。

有限差分法是从一般的物理现象出发建立相应的微分方程，经离散后得到差分方程，再进行求解的方法。

差分方程在计算机出现以前用一般的手摇计算器也可以求解。

20世纪60年代以后，由于有限单元法和边界元法的异军突起，使差分法在岩土工程中的应用暂时趋于
停滞，有限差分法曾一度受到冷遇。

但20世纪80年代末由美国ITASCA公司开发的FLAC
(Fast lagrangian analysis of continua )程序采用差分方法进行求解，在岩土工程数值计算中得到了广泛的应用，使差分法重新焕发出了活力。

岩土工程反分析中已有部分学者[177-179]
开始采用FLAC程序用于位移反分析中的正计算。

与其它程序相比，FLAC程序有如下
特点：① 完全动态运动方程使得FLAC在模拟物理上的不稳定过程时不存在数值上的障
碍；② 采用了显式有限差分求解，与有限元计算相比，FLAC具有较快的非线性求解速
度；③ 因为不需要形成刚度矩阵，故占用微机内存小，便于求解大型工程问题。

因此，在岩土工程的位移反分析中采用FLAC程序做正计算是非常合适的。

1.3.6 优化方法
位移反分析，可归结为一个极值问题的优化求解，合适的优化方法可以提高反分析的精
度和速度。

优化方法在岩土工程反分析中的应用，国内外很多学者进行了研究。

1980年Gioda
提出采用单纯形等优化方法求解岩体的弹性及弹塑性力学参数，并讨论了不同优化方法（单
纯形法、Powell法、Rosenbrok法）在岩土工程反分析中的适用性[180-184] ；Sakurai （1983）
采用最小二乘法反算隧洞围岩地应力及岩体弹性模量[126] ；1984年Arai采用二次梯度法
求解弹性模量和泊松比[132] ；Gioda（1987）提出了一个根据现场测量挡土结构位移来计算作用于墙体上土压力分布估计值的最小二乘法[185] ；Gioda（1987）等总结了适用于岩土工程反分析的四种优化方法，即单纯形法、拟梯度法、Rosenbrok法和Powell法，验证表明，
这四种方法计算量大、解的稳定性差、收敛速度慢[184];冯紫良（1989）提出了多种位移
反分析的计算加速方法[186];王芝银（1990）利用复合形法进行粘弹塑性增量位移反分析[162];胡维俊（1991）等利用高斯-牛顿法和阻尼最小二乘法反分析坝体的多个弹性模量和坝基的多个变形模量[187];李素华（1993）就不同的优化方法（单纯形加速法、复合形加速法、混合罚函数法和新鲍威尔法）在弹性横观各向同性以及弹塑性围岩位移反分析中的应用作了比较，并结合算例进行分析[170];刘维倩（1995）等结合实例利用乘子梯度法一次
反演初始地应力和材料参数并分析了算法的可行性及计算精度[188] ; Masumoto（1995）应用牛顿法反算三维渗透率的分布[189];吕爱钟（1996）结合6种最优化方法在巷道位移弹
性、弹塑性反分析中的应用，从初始参数初始点的选择、收敛速度、收敛精度和可靠性方面
评价了这几种方法的优劣[190] ; Ya ng（1997）[191]利用Powell法研究地铁结构引起的地
面沉降参数的反分析；沈振中（1997）提出了三维粘弹塑性位移反分析的可变容差法，并应用到三峡大坝的安全监测和反馈施工设计[163];陈国荣等（1996，1998）利用阻尼最小二
乘法进行三元件模型的粘弹性反分析，并应用到高速公路路基反分析及沉降预测[154, 155]; Ohkami （1997）利用牛顿法进行粘弹性参数辩识[156];李仲奎（1997）[215]利用Pattern Search 优化方法进行了二滩水电站地下厂房洞室群反馈分析；朱合华（1998）利用单纯形法反分析
成层土体的弹性模量，进而进行深基坑的变形预测[136]。

Ohkami（1999）利用非线性最小
二乘法进行粘弹性材料的参数辩识[192]。

以上为常规的优化方法，对高度非线性问题，搜索的最终结果为目标函数的极值点，并不能保证收敛到全局最优点，对反分析的结果影响较大。

而以遗传算法（Genetic Algorithm，GA ）为代表的现代优化方法，所要解决的是克服传统优化方法的不足，寻求问题的全局最优解，为岩土工程的优化反分析提供了新的方法。

遗传算法具有较强的鲁棒性和收敛到全局
最优的能力，并且可以处理非解析式的目标函数，在岩体力学位移反分析中得到了广泛的应
用[158, 175]。

但是，随着反分析参数的增加求解空间也急剧增加，而且岩体工程计算模型
的规模和复杂度也在不断增加，这常常会带来反分析计算量大、速度慢的问题，不能满足工
程上对于反分析的及时性需求。

因此提高遗传算法的运行速度便显得尤为突出，采用并行遗传算法（Parallel Genetic Algorithm，PGA ）是提高搜索效率的方法之一。

[105] 沈新普，徐秉业，岑章志，等.计算岩体力学•北京：中国环境科学出版社,1995.
[106] Jim Y T C, Liu S. Hydraulic tomography: Developme nt of a new aquifer test method.
Water Resources Research, 2000, 36(8): 2095-2105.
[107] Kernevez J P, Kn opf-Le noir C, Touzot G, et al. An ide ntificati on
method applied to an
orthotropic plate bending experime nt. Intern ati onal Jour nal for Numerical
Methods in
En gi neeri ng, 1978, 12(1): 129-139.
[108] 孙钧，蒋树屏，袁勇，等.岩土力学反演问题的随机理论与方法.汕头:
汕头大学出
版社,1996.
[109] Cividini A, Jurina L, Gioda G. Some aspects of characterization problems in geomech nics.
In ternatio nal Journal of Rock Mecha nics and Mining Scie nces, 1981, 18(6): 487-503.
[110] 杨林德，彭敏.地下洞室围岩位移量测的优化布置.同济大学学报：自然科学版,1995,
23(2): 129-133.
[111] 李守巨，马龙彪.混凝土重力坝弹性参数识别的灵敏度分析及实例.岩土工
程学报，
2000, 22(3): 381-383.
[112] Wardle L J, Gerrard C M. The equivale nt ani sotropic properties of layered rock and soil masses. Rock Mecha nics, 1972, 4(3): 155-175.
[113] Yeh W W G. Review of parameter identification procedures in groundwater hydrology: the
in verse problem. Water Resources Research, 1986, 22(2): 95-108.
[114] Mejias M M, Orla nde HRB, Ozisik M N. Effects of the
heating process and body
dime nsions on the estimati on of the thermal con ductivity comp onents of orthotropic solids.
In verse Problems in En gi neeri ng, 2003, 11(1): 75-89.
[115] Cunha J, Pira nda J. Applicati on of model updat ing tech niq ues in dyn amics for the identification of elastic constants of composite materials. Composites Part B (E ngi neeri ng), 1999, 30B(1): 79-85.
[116] Kavanagh K T. Experiment versus analysis: computational techniques for the descripti on
of static material resp on se. Intern ati onal Jour nal for Numerical Methods
in
En gi neeri ng, 1973, 5(4): 503-515.
[117] Cunha J, Piranda J. Identification of stiffness properties of composite tubes from dyn amic
tests. Experime ntal Mecha nics, 2000, 40(2): 211-218.
[118] Haftka R T, Scott E P, Cruz J R. Optimizatio n and experime nts:
A survey. Applied
Mecha nics Reviews, 1998, 51(7): 435-448.
[119] 邓建辉，李焯芬，葛修润.BP神经网络和遗传算法在岩石边坡位移反分析中的应用.
岩石力学与工程学报，2001,20(1): 1-5.
[120] 冯夏庭，张治强，杨成祥.位移反分析的进化神经网络方法研究.岩石力学与工程
学
报,1999, 18(5): 529-533.
[121] 刁心宏，冯夏庭，张士林，等.人工神经网络方法辨识岩体力学参数的可辨识性及其177
稳定性探讨•矿冶，2001, 10(3): 11-14.
[122] 冯夏庭.智能岩石力学导论.北京：科学出版社，2000.
[123] Kirste n HAD. Determ in ati on of rock mass elastic moduli by back an alysis of
deformati on measureme nt.〃Proc. Symp. Explorati on in Rock Engin eeri ng. Joha nn esburg, 1976: 1154-1160.
[124] Kovari K, Amstad C, Fritzc P. In tegrated measuri ng tech nique for rock pressure
determ in ati on. //Proc. Int. Conf. on Field Measureme nts in Rock Mecha nics. Zurich, 1977: 533-538.
[125] Sakurai S, Abe S. A desig n approach to dime nsioning un dergr ound ope ning. // Proc 3rd Int
Conf Numerical Methods in Geomecha nics. Aache n, 1979: 649-661.
[126] Sakurai S, Takeuchi K. Back an alysis of measured displaceme nt of
tunn el. Rock
Mecha nics and Rock En gi neeri ng, 1983, 16(3): 173-180.
[127] Gioda G. In direct ide ntificati on of the average elastic characterizati on
of rock masses.//
Proc. I nt Conf on Structural Fou ndation on Rock. Syd ney, 1980: 65-73.
[128] 杨志法，刘竹华.位移反分析法在地下工程设计中的初步应用.地下工程，1981, (2):
20-24.
[129] 杨志法.有限元法图谱.北京：科学出版社，1988.
[130] Jing F, Yang Z F. Tupu-Displaceme nt back an alysis method and its applicati on.// Proc Int
Co ng on Progress and Inno vation in Tunn eli ng. Toro nto, 1989: 2123-2126.
[131] Arai K, Ohta H, Yasui T. Simple optimizaitio n tech ni ques for evaluat ing deformati on
moduli from field observatio ns. Soils and Fou ndatio ns, 1983, 23(1): 107-113.
[132] Arai K. An in verse problems approach to the predictio n of Multi-dime nsion con solidati on
behavior. Soil and Fou ndati on, 1984, 24(1): 95-108.
[133] Arai K, Ohta H, Kojima K. Estimation of soil parameters based on monitored moveme nt of subsoil under con solidatio n. Soils and Fou ndatio ns, 1984, 24(4): 95-108.
[134] 杨林德，黄伟，王聿.初始地应力位移反分析的有限单元法.同济大学学报(自然科
学版),1985, (4): 15-20.
[135] Hisatake M. Three dime nsional back an alysis for tunn els.//Proceedi ngs
of the
Intern ati onal Symposium on Engin eeri ng in Complex Rock Formati ons. Beiji ng: Science Press, 1986: 791-797.
[136] 朱合华，桥本正.深基坑工程动态施工反演分析与变形预报.岩土工程学报，1998,
20(4): 30-35.
[137] 朱合华，傅德明，等.软土深基坑粘弹性动态增量反演分析与变形预测.岩
土力学，
2000, 21(4): 381-384.
[138] 石安池，徐卫亚，周家文，等.边坡弹性模量反分析的模拟退火BP网络方法.河海大
178
学学报：自然科学版，2006, 34(1): 69-73.
[139] Gioda G Maier G. Direct search soluti on of an in verse problem in elastic-plasticity,
ide ntificati on of cohesi on, friction an gle and in-situ stress by pressure tunnel tests.
In ternatio nal Journal for Numerical Methods in En gi neeri ng, 1980, 15: 1823-1848.
[140] Sakurai S, Akutagawa S, Tokudome O. Characterization of yield function and plastic
pote ntial function by back an alysis.// Ichikawa Y, Kawamoto T, eds. Proceedi ngs of the
8th Intern atio nal Conference on Computer Methods and Adva nces in Geomecha nics. Part
3 (of 4). Rotterdam: Balkema, 1994: 2011-2016.
[141] Sakurai S, Shimizu N, Matsumuro K. Evaluatio n of plastic zone around un dergro und ope nings by means of displaceme nt measureme nts.//Siriwarda ne H J, Zama n, eds.
Proceedi ngs of the Fifth In ternatio nal Conferen ce. Rotterdam: Balkema, 1985: 111-118.
[142] Zhe ng Y R, Wang C, Zha ng D H. Back an alysis from measured displaceme nt based on
elastoplastic theory in strain space.// Int Symp on Geom Bridges Strc. Lan zhou, 1987:
505-508.
[143] Anan darajah A, Agarwal D. Computer-aided calibrati on of a soil plasticity model.
Intern ati onal Journal for Numerical and An alytical Methods in Geomecha nics, 1991,
15(12): 835-856.
[144] 孙钧，黄伟.岩石力学参数弹塑性反演问题的优化问题.岩石力学与工程学
报，1992,
11(3): 221-229.
[145] 尹蓉蓉，朱合华.岩土介质弹塑性本构关系位移反分析.华东船舶工业学院
学报(自
然科学版),2004, 18(3): 21-25.
[146] 李宁，段小强，陈方方，等.围岩松动圈的弹塑性位移反分析方法探索 .岩
石力学与
工程学报，2006, 25(7): 1304-1308. [147] Hisatake M, Ito T. Back an alysis for tunn els by optimizati on
method.// Ichikawa Y,
Kawamoto T, eds. Proceed ings of 5th Intern ati onal Conference on Numerical Methods in Geomecha nics. Rotterdam: Balkema, 1985: 1301-1307.
[148] Wang Z Y, Liu H H. Back analysis of measured rheologic displacements of un dergro und
ope nin gs.〃 Proc 6th Conf on Num Meth in Geom. Austria, 1988: 2291-2297.
[149] Wang Z Y, Li Y P. Back an alysis of viscoparameters and strata stress in underground ope nin gs.// Proc. I nt. Symp. on Un dergrou nd Eng. New Delhi, 1988: 181-186.
[150] 王芝银，李云鹏.地下工程位移反分析法及程序.西安：陕西科学技术出版 社,1993. [151] Li Y P, Wang Z Y, Liu H H. Three dime nsio nal back an alysis
of viscoelastic creep
displaceme nts.// Proc 3rd Int Conf on Un dergro und Space and Earth Sheltered Buildi ngs. Sha nghai: To ngji Uni versity Press, 1988: 383-387.
[152] 李云鹏，王芝银.粘弹性位移反分析的边界元法.西安矿业学院学报，1989, 9(1):
17-23. 175 [92] Kirsch G Die Theorie der Elastizit?t und die Bed u rfnisse der Festigkeitslehre. Zeitschrift des Verei nes Deutscher In ge ni eure, 1898, (42): 797-807.
[93] Jaeger J C, Cook N. Fun dame ntals of Rock Mecha nics, Methuen & Co. London:
Methue n,
1969. Wiley & Sons, 1974. 216-222. Solids and Structures, 1982, 18(12): 1059-1074. tunnel in an elastic
half-pla ne. Intern atio nal Jour nal for Numerical and An alytical Methods in Geomecha nics, 1999, 21(2): 77-89.
[98] Exadaktylos G E, Stavropoulou M C. A closed-form elastic soluti on for stresses and
[94] Poulos H G, Davis E H. Elastic Solutio ns for soil and rock
mecha nics. New York:
Joh n [95] Pen der
Geotech ni que, 1980, J. Elastic
soluti ons for deep circular tu nnel.
30(2):
[96] Carter J
R. Elastic con solidati on
around
deep
circular
tu nnel.
Intern ati onal Jour nal of
[97] Verruijt A. A complex variable soluti on for a deform ing circular
displaceme nts around tunn els. Intern ati onal Jour nal of Rock Mecha nics and Mining Scie nces, 2002, 39(7): 905-916. [99]
Muskhelishvili N I. Some basic problems of the mathematical theory of elasticity.
Groningen: Noordhoff, 1953.
[100] 蔡晓鸿，蔡勇斌，蔡勇平，等.二向不等围压和内压作用下椭圆形洞室的计 算.地下 空间与工程学报,2008, 4(3): 453-459.
[101] 张路青，贾正雪.弹性位移反分析对地应力、弹模的反演唯一性.岩土工程学 报,2001, 23(2): 172-177. [102] Lekh nitskii S
G.
Theory of elasticity of an
ani sotropic
body.
Moscow: Mir Publishers,
1981.
[103] 吕爱钟.巷道开挖围岩参数及地应力可辨性的探讨 .岩石力学与工程学报
，
1988, 7(2):
155-164.
Computers and Geotech nics, 1997, 21(4): 255-275.
[157]
Ohkami T, Murai J I.
Ide ntificati on
procedure for viscoelastic
materials using
corresp ondence prin ciple. Communi cati ons in Numerical Methods in Engin eeri ng, 1998, 14(6): 497-504. [158] 高强，郭杏林,杨海天.遗传算法求解粘弹性反问题.大连理工大学学报,
2000, 40(6):
664-668.
[159] 王永岩,李剑光，魏佳，等.黏弹性有限元反分析方法及其在软岩流变问题 中的应用.
煤炭学报,2007, 32(11): 1162-1165. [160] Cividi ni A, Barla G Calibratio n of rheological material model on the basis of field
measureme nts.//
Ichikawa Y, Kawamoto T, eds. Proceed ings of 5th
[104] Ya ng Z F, of displaceme nts in
explorati on edits - Mecha nics
and Mining Scie nces 179 [153] 杨林德， 18-26.
[154]
陈国荣，
1996, 24(6):
25-28. [155] 陈国荣，i
6(4):
340-343.
[156] Ohkami
Lee C F, Wang S J.
prin ciples and applicati on.
,2000, 37(3): 525-533. 朱合华.地层三维粘弹性反演分析 池永斌.三维粘弹性参数反分析
Three-dime nsio nal back-a nalysis
Intern ati onal Journal of Rock
.岩土工程学报，1991, 13(6):
河海大学学报：自然科学版
高谦.高速公路路基性态反分析及沉降预报.工程地质学报，1998,
T, Ichikawa Y. A parameter
ide ntificati on
procedure for
Intern ati onal
Conference on Numerical Methods in Geomecha nics. Rotterdam: Balkema, 1985:
1621-1628.
[161] 王芝银，刘怀恒.粘弹塑性有限元分析及其在岩石力学与工程中的应用.西
安矿业学
院学报，1985, (1): 62-73.
[162] 王芝银，李云鹏.地下工程围岩粘弹塑性参数反分析.水利学报，1990, (9):
11-16.
[163] 沈振中.三维粘弹塑性位移反分析的可变容差法.水利学报，1997, (9): 66-70.
[164] Dziadziuszko P. In verse an alysis procedures for non li near and time-depe ndent con stitutive
models in geotech ni cs[Doctoral Dissertati on]. Nagoya: Nagoya Uni versity, 2000.
[165] 赵新铭，刘宁.岩体弹-粘塑性参数的Bayes随机优化反分析.南京航空航
天大学学报，
2008, 40(1): 125-128.
[166] Sakurai S. Numerical an alysis for the in terpretati on of field measurements in
geomecha nics.// Desai C S, Gioda G, eds. Numerical Methods and con stitutive Modelli ng in Geomecha nics. New York: Sprin ger Verlag, 1990: 351-407.
[167] Sakurai S. Field measurement and back analysis.//Beer G, Booker J R, Carter J P, eds.
Computer Methods and Adva nces in Geomecha nics. Rotterdam: Balkema, 1991:
1693-1701.
[168] Sakurai S. Assessme nt of cut slope stability by means of back an alysis of measured displaceme nts.// Pasamehmetoglu, Geunhan A, eds. Assessme nt and Preven tio n of failure Phe nomena in rock Engin eeri ng. Rotterdam: Balkema, 1993: 3-9.
[169] Ohkami T, Ichikawa Y, Kawamoto T. A boun dary eleme nt method for ide ntify ing
180
orthotropic material parameters. Intern atio nal Journal for Numerical and
An alytical
Methods in Geomecha nics, 1991, 15(9): 609-625.
[170] 李素华，朱维申.优化方法在弹性、横观各向同性以及弹塑性变形观测反分析中的应用.岩石力学与工程学报，1993, 12(2): 105-114.
[171] Dong Q X. Parameter ide ntificati on and its applicati on to tunn elli ng[Doctoral
Dissertati on ]. I nn sbruck: Inn sbruck Uni versity, 1997.
[172] Swoboda G, Ichikawa Y, Do ng Q, et al. Back an alysis of large
geotech ni cal models.
Intern ati onal Journal for Numerical and An alytical Methods in Geomecha nics, 1999,
23(13): 1455-1472.
[173] 刘学增,朱合华.考虑动态施工过程的岩土介质横观各向同性粘弹性反分析及其工
程应用.岩土工程学报,2002, 24(1): 89-92.
[174] 黄立新.平面正交各向异性体材料参数识别的边界元法[博士学位论文].北京:。