所谓岩土工程位移反

1、所谓岩土工程位移反分析，即以现场测量到的位移为基础，通过数学物理反分析模型，得到岩土介质的本构模型及等效力学参数（如初始地应力、变形参数、强度参数等）的方法。最终目的是建立一个输出位移更接近现场实测位移的理论模型，以便较正确地反映或预测岩土结构的某些力学行为。 20 世纪 70 年代初人们开始岩土工程位移反分析的研究，随着岩土工程的发展，国内外众多学者对位移反分析的理论与应用进行了大量广泛而深入的研究。岩土工程位移反分析涉及的研究内容非常广泛，下面就从位移解析解、位移反分析的唯一性、位移测量点的优化布置、本构模型、数值计算方法、优化方法这六个方面对其进行综合地考察。 1.3.1 位移解析解 1

2、898 年，Kirsch92最早发表了弹性平板中圆孔周围的二维应力分布解，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求解了不同地

3、应力条件下任意形状巷道的位移解析解。1.3.2 位移反分析的唯一性 反分析的唯一性是位移反分析中最重要却研究得最不充分的理论问题之一。迄今为止，国外尚未有相关论文发表，国内的论文也是凤毛麟角。吕爱钟（1988）103推导了参数可辨识条件，论证了地下洞室弹性位移反分析的多种唯一性问题，并指出某些问题无论安装多少个位移测点其反分析的结果都不是唯一的。张路青（2001）101进一步研究了考虑剪应力时位移反分析的唯一性问题。杨志法104则用几何作图法证明了图谱反分析的唯一性。以上文献中均假设岩体为各向同性材料，国内外尚未有报道对各向异性岩体位移反分析的唯一性问题进行研究。为了更大限度地利用位移反分析方

4、法及更有效地指导施工与决策，本文将在以上研究的基础上，对横观各向同性岩体位移反分析的唯一性问题进行研究。 1.3.3 位移测量点的优化布置 测量点的布置对位移反分析的唯一性和反分析精度都有很大的影响，因此，众多研究人员对位移测量点的优化布置进行了研究。 测量点的布置包括测量点的数目和测量点的空间位置两方面的工作。 由于工程费用和测量现场条件的制约，测量点的数目应该控制在一定范围之内。关于测量点数目下限的确定，大家公认的原则是测量点数目至少要大于等于待反演参数的数目，否则就会因为信息量不足而导致反分析失败。但是在满足下限的情况下，测量点数目越多越好，还是越少越好，目前并没有形成统一的认识。Kem

5、evez（1978）107发现并不是测点数目越多越好，而是测量点的空间位置更加重要。Cividini（1981）109的大量计算结果表明多测点的反分析结果不一定比少测点的反分析结果好，应该综合考虑测点的数量和空间位置。沈新普（1995）105认为测量误差会导致反分析结果偏离真值较远，应该尽可能多地布置测点，从而消除测量误差的影响。孙钧（1996）108从工程计算中观察到太多的测点并不能显著地改进位移反分析结果。Jim（2000）106认为较多的测点数目能够提高反分析的效果，但是随着测点数目增加得越多，提高的效果就越不明显。对测量点的空间位置应该遵循的原则，目前主要有以下几种174： 最大位移原

6、则110。该原则认为位移绝对值比较大的测点，测量的相对误差就比较小，测量精度就比较高，所得数据的实用价值就越大。不过，由于现代测量技术的迅猛发展，对量值比较小的位移，测量精度已经有了较大的提高。 最大灵敏度原则111, 112。灵敏度反映了位移测量值相对于待反演参数的变化。该原则考虑到灵敏度越大就越有利于参数的反演，所以就以灵敏度最大作为测点布置的原则。但是，当待反演的参数不仅一个时，依照这种原则很难得到一个适用的综合判断指标。 最小方差原则113, 114。按照该原则布置测点时遵循的原则就是使参数估计误差的方差最小。这与 Fisher 信息矩阵115有关，许多学者进行了大量测量点布置的研究1

7、14-117。在实际工程中，提出的大部分测点布置原则都与D矩阵的性质有关。主要的测点布置原则包括： D最优：使D矩阵的行列式值最小化； A最优：使D矩阵主对角元素的和最小化； E 最优：使D矩阵的最小特征值最大化。对比三个原则表明118，A最优原则倾向选择灵敏度最大的测点，E 最优原则倾向选择灵敏度最小的测点，而按照D最优原则选择的测点的灵敏度介于A最优和E 最优两者之间。以上原则分别从不同的侧面对测点的优化布置提供了很有价值的参考。但是以上的研究大多是针对各向同性岩体，并且到目前位置，关于位移反分析的测点布置准则，并没有得到统一的认识。本文试图基于最大位移原则对横观各向同性岩体位移反分析中的

8、测点优化布置问题进行初步研究，探寻合适的位移测点布置准则。 1.3.4 本构模型 纵观岩体位移反分析的发展过程，反分析中使用的有代表性的本构模型如表1.1表 1.5 所示。 需要指出的是，在岩土工程的反分析中，还有另外一种选取物理模型的思路，即放弃传统的各种力学本构关系，而将岩土体作为一个黑箱系统，直接从该系统的激励和响应的样本数据中总结提炼规律。其中最具代表性的是基于人工神经网络的智能方法，这种方法具有很强的自学习能力和对环境的适应能力，因此它特别适用于模拟岩土体这类复杂的系统。目前已有很多学者在这方面取得了一些成果119-122，但由于供给神经网络进行学习的大量样本不容易获得，从而影响了这

9、种模型的实际应用。 1.3.5 数值计算方法 正分析的正确解决是反分析的基础，而合适的数值计算方法可以改进正分析的速度、精度和通用性，所以在反分析的正计算中要采用一种合适的数值计算方法。 近年来，伴随着计算机技术的飞速发展，岩体稳定性分析的数值计算方法日臻成熟。当前应用于岩体工程问题的主要数值分析方法有：有限单元法、边界元法、有限差分法、离散单元法、无限元法、界面单元法、无单元法、非连续变形分析、流形元法以及由以上各种方法相组合而得到的混合数值计算方法。当岩体被裂隙切割成块体集合时，非连续的数值方法如离散单元法、非连续变形分析等可以更逼真地反映岩体的内部结构，但块体的拓扑分析过于繁杂，所以目前

10、在岩土工程的数值计算中，应用较广的还是基于连续介质力学的数值计算方法。 纵观岩体位移反分析方法的发展过程，数值计算方法大多选用有限单元法和边界元法。这是因为，有限单元法是岩石力学数值计算方法中最为广泛应用的一种。自20 世纪50 年代发展至今，有限元已成功求解了许多复杂的岩石力学与工程问题。有限元法的突出优点是适于处理非线性、非均质和复杂边界等问题，而岩体应力变形分析就恰恰存在这些困难问题。边界元法在20 世纪70 年代得到迅速发展，有限元法是全区域离散化，而边界元法仅对边界离散化。这样使三维问题降为二维问题求解，使二维问题降为一维问题求解，当物体的表面积和体积之比比较小时，边界元的划分单元数

11、要比有限元少数倍或十几倍，这样也使待解的方程数目、处理和存储的数据量降低同样的倍数，大大节省了机时。边界元法比较适合求解无限区域和半无限区域问题，如深埋巷道是一个典型的例子。 有限差分法是从一般的物理现象出发建立相应的微分方程，经离散后得到差分方程，再进行求解的方法。差分方程在计算机出现以前用一般的手摇计算器也可以求解。20 世纪 60 年代以后，由于有限单元法和边界元法的异军突起，使差分法在岩土工程中的应用暂时趋于停滞，有限差分法曾一度受到冷遇。但 20 世纪80 年代末由美国ITASCA 公司开发的FLAC（Fast lagrangian analysis of continua）程序采用

12、差分方法进行求解，在岩土工程数值计算中得到了广泛的应用，使差分法重新焕发出了活力。岩土工程反分析中已有部分学者177-179开始采用 FLAC 程序用于位移反分析中的正计算。 与其它程序相比，FLAC 程序有如下特点： 完全动态运动方程使得FLAC 在模拟物理上的不稳定过程时不存在数值上的障碍； 采用了显式有限差分求解，与有限元计算相比，FLAC 具有较快的非线性求解速度； 因为不需要形成刚度矩阵，故占用微机内存小，便于求解大型工程问题。因此，在岩土工程的位移反分析中采用FLAC 程序做正计算是非常合适的。 1.3.6 优化方法 位移反分析，可归结为一个极值问题的优化求解，合适的优化方法可以提

13、高反分析的精度和速度。优化方法在岩土工程反分析中的应用，国内外很多学者进行了研究。1980年Gioda提出采用单纯形等优化方法求解岩体的弹性及弹塑性力学参数，并讨论了不同优化方法（单纯形法、Powell 法、Rosenbrok 法）在岩土工程反分析中的适用性180-184；Sakurai（1983）采用最小二乘法反算隧洞围岩地应力及岩体弹性模量126；1984 年Arai 采用二次梯度法求解弹性模量和泊松比132；Gioda（1987）提出了一个根据现场测量挡土结构位移来计算作用于墙体上土压力分布估计值的最小二乘法185；Gioda（1987）等总结了适用于岩土工程反分析的四种优化方法，即单纯

14、形法、拟梯度法、Rosenbrok 法和 Powell 法，验证表明，这四种方法计算量大、解的稳定性差、收敛速度慢184；冯紫良（1989）提出了多种位移反分析的计 算加速方法186；王芝银（1990）利用复合形法进行粘弹塑性增量位移反分析162；胡维俊（1991）等利用高斯-牛顿法和阻尼最小二乘法反分析坝体的多个弹性模量和坝基的多个变形模量187；李素华（1993）就不同的优化方法（单纯形加速法、复合形加速法、混合罚函数法和新鲍威尔法）在弹性横观各向同性以及弹塑性围岩位移反分析中的应用作了比较，并结合算例进行分析170；刘维倩（1995）等结合实例利用乘子梯度法一次反演初始地应力和材料参数并

15、分析了算法的可行性及计算精度188；Masumoto（1995）应用牛顿法反算三维渗透率的分布189；吕爱钟（1996）结合6 种最优化方法在巷道位移弹性、弹塑性反分析中的应用，从初始参数初始点的选择、收敛速度、收敛精度和可靠性方面评价了这几种方法的优劣190；Yang（1997）191利用 Powell 法研究地铁结构引起的地面沉降参数的反分析；沈振中（1997）提出了三维粘弹塑性位移反分析的可变容差法，并应用到三峡大坝的安全监测和反馈施工设计163；陈国荣等（1996，1998）利用阻尼最小二乘法进行三元件模型的粘弹性反分析，并应用到高速公路路基反分析及沉降预测154, 155；Ohkam

16、i（1997）利用牛顿法进行粘弹性参数辩识156；李仲奎（1997）215利用Pattern Search优化方法进行了二滩水电站地下厂房洞室群反馈分析；朱合华（1998）利用单纯形法反分析成层土体的弹性模量，进而进行深基坑的变形预测136。Ohkami（1999）利用非线性最小二乘法进行粘弹性材料的参数辩识192。以上为常规的优化方法，对高度非线性问题，搜索的最终结果为目标函数的极值点，并不能保证收敛到全局最优点，对反分析的结果影响较大。而以遗传算法（Genetic Algorithm，GA）为代表的现代优化方法，所要解决的是克服传统优化方法的不足，寻求问题的全局最优解，为岩土工程的优化反分

17、析提供了新的方法。遗传算法具有较强的鲁棒性和收敛到全局最优的能力，并且可以处理非解析式的目标函数，在岩体力学位移反分析中得到了广泛的应用158, 175。但是，随着反分析参数的增加求解空间也急剧增加，而且岩体工程计算模型的规模和复杂度也在不断增加，这常常会带来反分析计算量大、速度慢的问题，不能满足工程上对于反分析的及时性需求。因此提高遗传算法的运行速度便显得尤为突出，采用并行遗传算法（Parallel Genetic Algorithm，PGA）是提高搜索效率的方法之一。 105 沈新普, 徐秉业, 岑章志, 等. 计算岩体力学. 北京: 中国环境科学出版社, 1995. 106 Jim Y

18、T C, Liu S. Hydraulic tomography: Development of a new aquifer test method. Water Resources Research, 2000, 36(8): 2095-2105. 107 Kernevez J P, Knopf-Lenoir C, Touzot G, et al. An identification method applied to an orthotropic plate bending experiment. International Journal for Numerical Methods in

19、 Engineering, 1978, 12(1): 129-139. 108 孙钧, 蒋树屏, 袁勇, 等. 岩土力学反演问题的随机理论与方法. 汕头: 汕头大学出版社, 1996. 109 Cividini A, Jurina L, Gioda G. Some aspects of characterization problems in geomechnics. International Journal of Rock Mechanics and Mining Sciences, 1981, 18(6): 487-503. 110 杨林德, 彭敏. 地下洞室围岩位移量测的优化布置. 同

20、济大学学报：自然科学版, 1995, 23(2): 129-133. 111 李守巨, 马龙彪. 混凝土重力坝弹性参数识别的灵敏度分析及实例. 岩土工程学报, 2000, 22(3): 381-383. 112 Wardle L J, Gerrard C M. The equivalent anisotropic properties of layered rock and soil masses. Rock Mechanics, 1972, 4(3): 155-175. 113 Yeh W W G. Review of parameter identification procedures

21、in groundwater hydrology: the inverse problem. Water Resources Research, 1986, 22(2): 95-108. 114 Mejias M M, Orlande H R B, Ozisik M N. Effects of the heating process and body dimensions on the estimation of the thermal conductivity components of orthotropic solids. Inverse Problems in Engineering,

22、 2003, 11(1): 75-89. 115 Cunha J, Piranda J. Application of model updating techniques in dynamics for the identification of elastic constants of composite materials. Composites Part B (Engineering), 1999, 30B(1): 79-85. 116 Kavanagh K T. Experiment versus analysis: computational techniques for the d

23、escription of static material response. International Journal for Numerical Methods in Engineering, 1973, 5(4): 503-515. 117 Cunha J, Piranda J. Identification of stiffness properties of composite tubes from dynamic tests. Experimental Mechanics, 2000, 40(2): 211-218. 118 Haftka R T, Scott E P, Cruz

24、 J R. Optimization and experiments: A survey. Applied Mechanics Reviews, 1998, 51(7): 435-448. 119 邓建辉, 李焯芬, 葛修润. BP 神经网络和遗传算法在岩石边坡位移反分析中的应用. 岩石力学与工程学报, 2001, 20(1): 1-5. 120 冯夏庭, 张治强, 杨成祥. 位移反分析的进化神经网络方法研究. 岩石力学与工程学报, 1999, 18(5): 529-533. 121 刁心宏, 冯夏庭, 张士林, 等. 人工神经网络方法辨识岩体力学参数的可辨识性及其 177 稳定性探讨. 矿冶

25、, 2001, 10(3): 11-14. 122 冯夏庭. 智能岩石力学导论. 北京: 科学出版社, 2000. 123 Kirsten H A D. Determination of rock mass elastic moduli by back analysis of deformation measurement./ Proc. Symp. Exploration in Rock Engineering. Johannesburg, 1976: 1154-1160. 124 Kovari K, Amstad C, Fritzc P. Integrated measuring tech

26、nique for rock pressure determination./Proc. Int. Conf. on Field Measurements in Rock Mechanics. Zurich, 1977: 533-538. 125 Sakurai S, Abe S. A design approach to dimensioning underground opening. / Proc 3rd Int Conf Numerical Methods in Geomechanics. Aachen, 1979: 649-661. 126 Sakurai S, Takeuchi K

27、. Back analysis of measured displacement of tunnel. Rock Mechanics and Rock Engineering, 1983, 16(3): 173-180. 127 Gioda G. Indirect identification of the average elastic characterization of rock masses./ Proc. Int Conf on Structural Foundation on Rock. Sydney, 1980: 65-73. 128 杨志法, 刘竹华. 位移反分析法在地下工程

28、设计中的初步应用. 地下工程, 1981, (2): 20-24. 129 杨志法. 有限元法图谱. 北京: 科学出版社, 1988. 130 Jing F, Yang Z F. Tupu-Displacement back analysis method and its application./ Proc Int Cong on Progress and Innovation in Tunneling. Toronto, 1989: 2123-2126. 131 Arai K, Ohta H, Yasui T. Simple optimizaition techniques for eva

29、luating deformation moduli from field observations. Soils and Foundations, 1983, 23(1): 107-113. 132 Arai K. An inverse problems approach to the prediction of Multi-dimension consolidation behavior. Soil and Foundation, 1984, 24(1): 95-108. 133 Arai K, Ohta H, Kojima K. Estimation of soil parameters

30、 based on monitored movement of subsoil under consolidation. Soils and Foundations, 1984, 24(4): 95-108. 134 杨林德, 黄伟, 王聿. 初始地应力位移反分析的有限单元法. 同济大学学报(自然科学版), 1985, (4): 15-20. 135 Hisatake M. Three dimensional back analysis for tunnels./Proceedings of the International Symposium on Engineering in Compl

31、ex Rock Formations. Beijing: 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. Dire

32、ct search solution of an inverse problem in elastic-plasticity, identification of cohesion, friction angle and in-situ stress by pressure tunnel tests. International Journal for Numerical Methods in Engineering, 1980, 15: 1823-1848. 140 Sakurai S, Akutagawa S, Tokudome O. Characterization of yield f

33、unction and plastic potential function by back analysis./ Ichikawa Y, Kawamoto T, eds. Proceedings of the 8th International Conference on Computer Methods and Advances in Geomechanics. Part 3 (of 4). Rotterdam: Balkema, 1994: 2011-2016. 141 Sakurai S, Shimizu N, Matsumuro K. Evaluation of plastic zo

34、ne around underground openings by means of displacement measurements./Siriwardane H J, Zaman, eds. Proceedings of the Fifth International Conference. Rotterdam: Balkema, 1985: 111-118. 142 Zheng Y R, Wang C, Zhang D H. Back analysis from measured displacement based on elastoplastic theory in strain

35、space./ Int Symp on Geom Bridges Strc. Lanzhou, 1987: 505-508. 143 Anandarajah A, Agarwal D. Computer-aided calibration of a soil plasticity model. International Journal for Numerical and Analytical Methods in Geomechanics, 1991, 15(12): 835-856. 144 孙钧, 黄伟. 岩石力学参数弹塑性反演问题的优化问题. 岩石力学与工程学报, 1992, 11(3

36、): 221-229. 145 尹蓉蓉, 朱合华. 岩土介质弹塑性本构关系位移反分析. 华东船舶工业学院学报(自然科学版), 2004, 18(3): 21-25. 146 李宁, 段小强, 陈方方, 等. 围岩松动圈的弹塑性位移反分析方法探索. 岩石力学与工程学报, 2006, 25(7): 1304-1308. 147 Hisatake M, Ito T. Back analysis for tunnels by optimization method./ Ichikawa Y, Kawamoto T, eds. Proceedings of 5th International Confe

37、rence on Numerical Methods in Geomechanics. Rotterdam: Balkema, 1985: 1301-1307. 148 Wang Z Y, Liu H H. Back analysis of measured rheologic displacements of underground openings./ Proc 6th Conf on Num Methin Geom. Austria, 1988: 2291-2297. 149 Wang Z Y, Li Y P. Back analysis of viscoparameters and s

38、trata stress in underground openings./ Proc. Int. Symp. on Underground Eng. New Delhi, 1988: 181-186. 150 王芝银, 李云鹏. 地下工程位移反分析法及程序. 西安: 陕西科学技术出版社, 1993. 151 Li Y P, Wang Z Y, Liu H H. Three dimensional back analysis of viscoelastic creep displacements./ Proc 3rd Int Conf on Underground Space and Eart

39、h Sheltered Buildings. Shanghai: Tongji University Press, 1988: 383-387. 152 李云鹏, 王芝银. 粘弹性位移反分析的边界元法. 西安矿业学院学报, 1989, 9(1): 17-23. 175 92 Kirsch G. Die Theorie der Elastizität und die Bedürfnisse der Festigkeitslehre. Zeitschrift des Vereines Deutscher Ingenieure, 1898, (42): 797-807. 93 J

40、aeger J C, Cook N. Fundamentals of Rock Mechanics, Methuen & Co. London: Methuen, 1969. 94 Poulos H G, Davis E H. Elastic Solutions for soil and rock mechanics. New York: John Wiley & Sons, 1974. 95 Pender M J. Elastic solutions for a deep circular tunnel. Geotechnique, 1980, 30(2): 216-222.

41、 96 Carter J R. Elastic consolidation around a deep circular tunnel. International Journal of Solids and Structures, 1982, 18(12): 1059-1074. 97 Verruijt A. A complex variable solution for a deforming circular tunnel in an elastic half-plane. International Journal for Numerical and Analytical Method

42、s in Geomechanics, 1999, 21(2): 77-89. 98 Exadaktylos G E, Stavropoulou M C. A closed-form elastic solution for stresses and displacements around tunnels. International Journal of Rock Mechanics and Mining Sciences, 2002, 39(7): 905-916. 99 Muskhelishvili N I. Some basic problems of the mathematical

43、 theory of elasticity. Groningen: Noordhoff, 1953. 100 蔡晓鸿, 蔡勇斌, 蔡勇平, 等. 二向不等围压和内压作用下椭圆形洞室的计算. 地下空间与工程学报, 2008, 4(3): 453-459. 101 张路青, 贾正雪. 弹性位移反分析对地应力、弹模的反演唯一性. 岩土工程学报, 2001, 23(2): 172-177. 102 Lekhnitskii S G. Theory of elasticity of an anisotropic body. Moscow: Mir Publishers, 1981. 103 吕爱钟. 巷道

44、开挖围岩参数及地应力可辨性的探讨. 岩石力学与工程学报, 1988, 7(2): 155-164. 104 Yang Z F, Lee C F, Wang S J. Three-dimensional back-analysis of displacements in exploration edits - principles and application. International Journal of Rock Mechanics and Mining Sciences , 2000, 37(3): 525-533. 179 153 杨林德, 朱合华. 地层三维粘弹性反演分析. 岩土

45、工程学报, 1991, 13(6): 18-26. 154 陈国荣, 池永斌. 三维粘弹性参数反分析. 河海大学学报：自然科学版, 1996, 24(6): 25-28. 155 陈国荣, 高谦. 高速公路路基性态反分析及沉降预报. 工程地质学报, 1998, 6(4): 340-343. 156 Ohkami T, Ichikawa Y. A parameter identification procedure for viscoelastic materials. Computers and Geotechnics, 1997, 21(4): 255-275. 157 Ohkami T,

46、Murai J I. Identification procedure for viscoelastic materials using correspondence principle. Communications in Numerical Methods in Engineering, 1998, 14(6): 497-504. 158 高强, 郭杏林, 杨海天. 遗传算法求解粘弹性反问题. 大连理工大学学报, 2000, 40(6): 664-668. 159 王永岩, 李剑光, 魏佳, 等. 黏弹性有限元反分析方法及其在软岩流变问题中的应用. 煤炭学报, 2007, 32(11):

47、1162-1165. 160 Cividini A, Barla G. Calibration of rheological material model on the basis of field measurements./ Ichikawa Y, Kawamoto T, eds. Proceedings of 5th International Conference on Numerical Methods in Geomechanics. Rotterdam: Balkema, 1985: 1621-1628. 161 王芝银, 刘怀恒. 粘弹塑性有限元分析及其在岩石力学与工程中的应用

48、. 西安矿业学院学报, 1985, (1): 62-73. 162 王芝银, 李云鹏. 地下工程围岩粘弹塑性参数反分析. 水利学报, 1990, (9): 11-16. 163 沈振中. 三维粘弹塑性位移反分析的可变容差法. 水利学报, 1997, (9): 66-70. 164 Dziadziuszko P. Inverse analysis procedures for nonlinear and time-dependent constitutive models in geotechnicsDoctoral Dissertation. Nagoya: Nagoya University

49、, 2000. 165 赵新铭, 刘宁. 岩体弹-粘塑性参数的 Bayes 随机优化反分析. 南京航空航天大学学报, 2008, 40(1): 125-128. 166 Sakurai S. Numerical analysis for the interpretation of field measurements in geomechanics./ Desai C S, Gioda G, eds. Numerical Methods and constitutive Modelling in Geomechanics. New York: Springer Verlag, 1990: 35

50、1-407. 167 Sakurai S. Field measurement and back analysis./Beer G, Booker J R, Carter J P, eds. Computer Methods and Advances in Geomechanics. Rotterdam: Balkema, 1991: 1693-1701. 168 Sakurai S. Assessment of cut slope stability by means of back analysis of measured displacements./ Pasamehmetoglu, G

51、eunhan A, eds. Assessment and Prevention of failure Phenomena in rock Engineering. Rotterdam: Balkema, 1993: 3-9. 169 Ohkami T, Ichikawa Y, Kawamoto T. A boundary element method for identifying 180 orthotropic material parameters. International Journal for Numerical and Analytical Methods inGeomecha

52、nics, 1991, 15(9): 609-625. 170 李素华, 朱维申. 优化方法在弹性、横观各向同性以及弹塑性变形观测反分析中的应用. 岩石力学与工程学报, 1993, 12(2): 105-114. 171 Dong Q X. Parameter identification and its application to tunnellingDoctoral Dissertation. Innsbruck: Innsbruck University, 1997. 172 Swoboda G, Ichikawa Y, Dong Q, et al. Back analysis of

53、large geotechnical models. International Journal for Numerical and Analytical Methods in Geomechanics, 1999, 23(13): 1455-1472. 173 刘学增, 朱合华. 考虑动态施工过程的岩土介质横观各向同性粘弹性反分析及其工程应用. 岩土工程学报, 2002, 24(1): 89-92. 174 黄立新. 平面正交各向异性体材料参数识别的边界元法博士学位论文. 北京: 清华大学, 2005. 175 邝宏柱, 刘学增. 层状地层横观各向同性粘弹性优化反分析. 地下空间与工程学报,

54、 2006, 2(1): 112-114. 176 黄光明, 李云, 顾冲时, 等. 碾压混凝土坝横观各向同性粘弹性参数反演. 水利水运工程学报, 2006, (4): 15-20. 177 姜谙男, 茹忠亮, 张娇. 基于粒子群算法和 FLAC 的洞室围岩参数反分析. 矿业研究与开发, 2007, 27(5): 33-35. 178 罗润林, 阮怀宁, 黄亚哲, 等. 岩体初始地应力场的粒子群优化反演及在 FLAC（-3D）中的实现. 长江科学院院报, 2008, 25(4): 73-76. 179 张志增. 白云东矿山高陡边坡的稳定性分析和加固技术研究硕士学位论文. 北京: 北京科技大学大学, 2006. 180 Gioda G. Indirect identification of the average elastic characterization of rock masses./ Proc Int