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

所谓岩土工程位移反分析， 即以现场测量到的位移为基础， 通过数学物理反分析模型， 得到岩土介质的本构模型及等效力学参数（如初始地应力、变形参数、强度参数等）的方法。最终目的是建立一个输出位移更接近现场实测位移的理论模型， 以便较正确地反映或预测岩土结构的某些力学行为。20世纪 70年代初人们开始岩土工程位移反分析的研究，随着岩土工程的发展，国内外众多学者对位移反分析的理论与应用进行了大量广泛而深入的研究。 岩土工程位移反分析涉及的研究内容非常广泛， 下面就从位移解析解、 位移反分析的唯一性、 位移测量点的优化布置、本构模型、数值计算方法、优化方法这六个方面对其进行综合地考察。1.3.1 位移解析解1898年，Kirsch[92]最早发表了弹性平板中圆孔周围的二维应力分布解， Jaeger和Cook（1969）[93]对Kirsch 方程进行了详细的推导。 此后，Poulos和Davi（s1974）[94]、Pender1980）[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.3.5数值计算方法正分析的正确解决是反分析的基础，而合适的数值计算方法可以改进正分析的速度、精度和通用性，所以在反分析的正计算中要采用一种合适的数值计算方法。近年来，伴随着计算机技术的飞速发展，岩体稳定性分析的数值计算方法日臻成熟。当前应用于岩体工程问题的主要数值分析方法有：有限单元法、边界元法、有限差分法、离散单元法、无限元法、界面单元法、无单元法、非连续变形分析、流形元法以及由以上各种方法相组合而得到的混合数值计算方法。当岩体被裂隙切割成块体集合时，非连续的数值方法如离散单元法、非连续变形分析等可以更逼真地反映岩体的内部结构，但块体的拓扑分析过于繁杂，所以目前在岩土工程的数值计算中，应用较广的还是基于连续介质力学的数值计算方法。纵观岩体位移反分析方法的发展过程，数值计算方法大多选用有限单元法和边界元法。这是因为，有限单元法是岩石力学数值计算方法中最为广泛应用的一种。自20世纪50年代发展至今，有限元已成功求解了许多复杂的岩石力学与工程问题。有限元法的突出优点是适于处理非线性、非均质和复杂边界等问题，而岩体应力变形分析就恰恰存在这些困难问题。边界元法在20世纪70年代得到迅速发展，有限元法是全区域离散化，而边界元法仅对边界离散化。这样使三维问题降为二维问题求解，使二维问题降为一维问题求解，当物体的表面积和体积之比比较小时，边界元的划分单元数要比有限元少数倍或十几倍，这样也使待解的方程数目、处理和存储的数据量降低同样的倍数，大大节省了机时。边界元法比较适合求解无限区域和半无限区域问题，如深埋巷道是一个典型的例子。有限差分法是从一般的物理现象出发建立相应的微分方程，经离散后得到差分方程，再进行求解的方法。差分方程在计算机出现以前用一般的手摇计算器也可以求解。20世纪60年代以后，由于有限单元法和边界元法的异军突起，使差分法在岩土工程中的应用暂时趋于停滞，有限差分法曾一度受到冷遇。但20世纪80年代末由美国ITASCA公司开发的FLACFastlagrangiananalysisofcontinua）程序采用差分方法进行求解，在岩土工程数值计算中得到了广泛的应用，使差分法重新焕发出了活力。岩土工程反分析中已有部分学者[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）利用复合形法进行粘弹塑性增量位移反分析；胡维俊（1991）等利用高斯-牛顿法和阻尼最小二乘法反分析坝体的多个弹性模量和坝基的多个变形模量[187]；李素华（1993）就不同的优化方法（单纯形加速法、复合形加速法、混合罚函数法和新鲍威尔法）在弹性横观各向同性以及弹塑性围岩位移反分析中的应用作了比较，并结合算例进行分析[170]；刘维倩（1995）等结合实例利用乘子梯度法一次反演初始地应力和材料参数并分析了算法的可行性及计算精度[188]；Masumoto（1995）应用牛顿法反算三维渗透率的分布[189]；吕爱钟（1996）结合6种最优化方法在巷道位移弹性、弹塑性反分析中的应用，从初始参数初始点的选择、收敛速度、收敛精度和可靠性方面评价了这几种方法的优劣[190]；Yang（1997）[191]利用Powell法研究地铁结构引起的地面沉降参数的反分析；沈振中（1997）提出了三维粘弹塑性位移反分析的可变容差法，并应用到三峡大坝的安全监测和反馈施工设计[163]；陈国荣等（1996，1998）利用阻尼最小二乘法进行三元件模型的粘弹性反分析，并应用到高速公路路基反分析及沉降预测[154,155]；Ohkami（1997）利用牛顿法进行粘弹性参数辩识[156]；李仲奎（1997）[215]利用PatternSearch优化方法进行了二滩水电站地下厂房洞室群反馈分析；朱合华（1998）利用单纯形法反分析成层土体的弹性模量，进而进行深基坑的变形预测[136]。Ohkami（1999）利用非线性最小二乘法进行粘弹性材料的参数辩识[192]。以上为常规的优化方法，对高度非线性问题，搜索的最终结果为目标函数的极值点，并不能保证收敛到全局最优点，对反分析的结果影响较大。而以遗传算法（GeneticAlgorithm，GA）为代表的现代优化方法，所要解决的是克服传统优化方法的不足，寻求问题的全局最优解，为岩土工程的优化反分析提供了新的方法。遗传算法具有较强的鲁棒性和收敛到全局最优的能力，并且可以处理非解析式的目标函数，在岩体力学位移反分析中得到了广泛的应用[158,175]。但是，随着反分析参数的增加求解空间也急剧增加，而且岩体工程计算模型的规模和复杂度也在不断增加，这常常会带来反分析计算量大、速度慢的问题，不能满足工程上对于反分析的及时性需求。因此提高遗传算法的运行速度便显得尤为突出，采用并行遗传算法（ParallelGeneticAlgorithm，PGA）是提高搜索效率的方法之一。沈新普,徐秉业,岑章志,等.计算岩体力学.北京:中国环境科学出版社,1995.[106]JimYTC,LiuS.Hydraulictomography:Developmentofanewaquifertestmethod.WaterResourcesResearch,2000,36(8):2095-2105.[107]KernevezJP,Knopf-LenoirC,TouzotG,etal.Anidentificationmethodappliedtoanorthotropic plate bending experiment. International Journal for NumericalMethods inEngineering,1978,12(1):129-139.[108] 孙钧, 蒋树屏, 袁勇, 等. 岩土力学反演问题的随机理论与方法 . 汕头:汕头大学出版社,1996.CividiniA,JurinaL,GiodaG.Someaspectsofcharacterizationproblemsingeomechnics.InternationalJournalofRockMechanicsandMiningSciences,1981,18(6):487-503.[110] 杨林德, 彭敏. 地下洞室围岩位移量测的优化布置 . 同济大学学报：自然科学版,1995,23(2):129-133.[111] 李守巨, 马龙彪. 混凝土重力坝弹性参数识别的灵敏度分析及实例 . 岩土工程学报,2000,22(3):381-383.WardleLJ,GerrardCM.Theequivalentanisotropicpropertiesoflayeredrockandsoilmasses.RockMechanics,1972,4(3):155-175.YehWWG.Reviewofparameteridentificationproceduresingroundwaterhydrology:theinverseproblem.WaterResourcesResearch,1986,22(2):95-108.[114] Mejias M M, Orlande H R B, Ozisik M N. Effects of theheating process and bodydimensionsontheestimationofthethermalconductivitycomponentsoforthotropicsolids.InverseProblemsinEngineering,2003,11(1):75-89.CunhaJ,PirandaJ.Applicationofmodelupdatingtechniquesindynamicsfortheidentificationofelasticconstantsofcompositematerials.CompositesPartB(Engineering),1999,30B(1):79-85.KavanaghKT.Experimentversusanalysis:computationaltechniquesforthedescriptionof static material response. International Journal for Numerical MethodsinEngineering,1973,5(4):503-515.CunhaJ,PirandaJ.Identificationofstiffnesspropertiesofcompositetubesfromdynamictests.ExperimentalMechanics,2000,40(2):211-218.[118] Haftka R T, Scott E P, Cruz J R. Optimization and experiments:A survey. Applied邓建辉,李焯芬,葛修润.BP神经网络和遗传算法在岩石边坡位移反分析中的应用.岩石力学与工程学报

,2001,20(1):1-5.[120] 冯夏庭

,

张治强

,

杨成祥

.

位移反分析的进化神经网络方法研究

.

岩石力学与工程学报,1999,18(5):529-533.刁心宏,冯夏庭,张士林,等.人工神经网络方法辨识岩体力学参数的可辨识性及其177稳定性探讨. 矿冶,2001,10(3):11-14.冯夏庭.智能岩石力学导论.北京:科学出版社,2000.[123] Kirsten H A D. Determination of rock mass elastic moduli byback analysis ofdeformationmeasurement.//Proc.Symp.ExplorationinRockEngineering.Johannesburg,1976:1154-1160.[124] Kovari K, Amstad C, Fritzc P. Integrated measuring technique forrock pressuredetermination.//Proc.Int.Conf.onFieldMeasurementsinRockMechanics.Zurich,1977:533-538.SakuraiS,AbeS.Adesignapproachtodimensioningundergroundopening.//Proc3rdIntConfNumericalMethodsinGeomechanics.Aachen,1979:649-661.[126] Sakurai S, Takeuchi K. Back analysis of measured displacement oftunnel. RockMechanicsandRockEngineering,1983,16(3):173-180.[127] Gioda G. Indirect identification of the average elastic characterizationof rock masses.//Proc.IntConfonStructuralFoundationonRock.Sydney,1980:65-73.[128] 杨志法, 刘竹华. 位移反分析法在地下工程设计中的初步应用 . 地下工程,1981, (2):20-24.杨志法.有限元法图谱.北京:科学出版社,1988.JingF,YangZF.Tupu-Displacementbackanalysismethodanditsapplication.//ProcIntCongonProgressandInnovationinTunneling.Toronto,1989:2123-2126.[131] Arai K, Ohta H, Yasui T. Simple optimizaition techniques forevaluating deformationmodulifromfieldobservations.SoilsandFoundations,1983,23(1):107-113.AraiK.AninverseproblemsapproachtothepredictionofMulti-dimensionconsolidationbehavior.SoilandFoundation,1984,24(1):95-108.AraiK,OhtaH,KojimaK.Estimationofsoilparametersbasedonmonitoredmovementofsubsoilunderconsolidation.SoilsandFoundations,1984,24(4):95-108.[134]杨林德,黄伟,王聿.初始地应力位移反分析的有限单元法.同济大学学报(自然科学版),1985,(4):15-20.[135]HisatakeM.Threedimensionalbackanalysisfortunnels.//Proceedingsof theInternationalSymposiumonEngineeringinComplexRockFormations.Beijing:SciencePress,1986:791-797.朱合华,桥本正.深基坑工程动态施工反演分析与变形预报.岩土工程学报,1998,20(4):30-35.[137] 朱合华

,

傅德明

,

等.

软土深基坑粘弹性动态增量反演分析与变形预测

.

岩土力学

,2000,21(4):381-384.[138] 石安池, 徐卫亚

,

周家文

,

等.

边坡弹性模量反分析的模拟退火

BP

网络方法.

河海大178in

学学报：自然科学版[139] Gioda G,elastic-plasticity,identification of

,2006,34(1):69-73.Maier G. Directcohesion, friction

searchangle

solution ofand in-situ

an inversestress by

problempressuretunnel tests.InternationalJournalforNumericalMethodsinEngineering,1980,15:1823-1848.[140]SakuraiS,AkutagawaS,TokudomeO.Characterizationofyieldfunctionandplasticpotentialfunctionbybackanalysis.//IchikawaY,KawamotoT,eds.Proceedingsofthe8thInternationalConferenceonComputerMethodsandAdvancesinGeomechanics.Part3(of4).Rotterdam:Balkema,1994:2011-2016.[141]SakuraiS,ShimizuN,MatsumuroK.Evaluationofplasticzonearoundundergroundopenings by means of displacement measurements.//Siriwardane H J, Zaman,eds.ProceedingsoftheFifthInternationalConference.Rotterdam:Balkema,1985:111-118.ZhengYR,WangC,ZhangDH.Backanalysisfrommeasureddisplacementbasedonelastoplastictheoryinstrainspace.//IntSymponGeomBridgesStrc.Lanzhou,1987:505-508.[143]AnandarajahA,AgarwalD.Computer-aidedcalibrationofasoilplasticitymodel.InternationalJournalforNumericalandAnalyticalMethodsinGeomechanics,1991,15(12):835-856.[144]孙钧,黄伟.岩石力学参数弹塑性反演问题的优化问题.岩石力学与工程学报,1992,11(3):221-229.尹蓉蓉,朱合华.岩土介质弹塑性本构关系位移反分析.华东船舶工业学院学报(自然科学版),2004,18(3):21-25.[146] 李宁, 段小强, 陈方方, 等. 围岩松动圈的弹塑性位移反分析方法探索 . 岩石力学与工程学报,2006,25(7):1304-1308.[147]HisatakeM,ItoT.Backanalysisfortunnelsbyoptimizationmethod.//IchikawaY,KawamotoT,eds.Proceedingsof5thInternationalConferenceonNumericalMethodsinGeomechanics.Rotterdam:Balkema,1985:1301-1307.WangZY,LiuHH.Backanalysisofmeasuredrheologicdisplacementsofundergroundopenings.//Proc6thConfonNumMethinGeom.Austria,1988:2291-2297.[149] Wang Z Y, Li Y P. Back analysis of viscoparameters and stratastress in undergroundopenings.//Proc.Int.Symp.onUndergroundEng.NewDelhi,1988:181-186.王芝银,李云鹏.地下工程位移反分析法及程序.西安:陕西科学技术出版社,1993.[151] Li Y P, Wangof viscoelastic creep

Z Y, Liu H H. Three

dimensional back analysisdisplacements.//Proc3rdIntConfonUndergroundSpaceandEarthShelteredBuildings.Shanghai:TongjiUniversityPress,1988:383-387.[152] 李云鹏, 王芝银. 粘弹性位移反分析的边界元法 . 西安矿业学院学报 , 1989,9(1):17-23.175[92]KirschG.DieTheoriederElastizit?tunddieBedürfnissederFestigkeitslehre.ZeitschriftdesVereinesDeutscherIngenieure,1898,(42):797-807.JaegerJC,CookN.FundamentalsofRockMechanics,Methuen&Co.London:Methuen,1969.[94] Poulos H G,Davis E H. Elastic Solutions for soil and rockmechanics. New York: JohnWiley&Sons,1974.[95] Pender M J. Elastic solutions for a deep circular tunnel.Geotechnique, 1980, 30(2):216-222.CarterJR.Elasticconsolidationaroundadeepcirculartunnel.InternationalJournalofSolidsandStructures,1982,18(12):1059-1074.[97]VerruijtA.Acomplexvariablesolutionforadeformingcirculartunnelinanelastichalf-plane.InternationalJournalforNumericalandAnalyticalMethodsinGeomechanics,1999,21(2):77-89.[98]ExadaktylosGE,StavropoulouMC.Aclosed-formelasticsolutionfor stresses anddisplacements around tunnels. International Journal of Rock Mechanics andMiningSciences,2002,39(7):905-916.[99] Muskhelishvili N I. Some basic problems of the mathematical theoryof elasticity.Groningen:Noordhoff,1953.蔡晓鸿,蔡勇斌,蔡勇平,等.二向不等围压和内压作用下椭圆形洞室的计算.地下空间与工程学报 ,2008,4(3):453-459.张路青,贾正雪.弹性位移反分析对地应力、弹模的反演唯一性.岩土工程学报,2001,23(2):172-177.[102]LekhnitskiiSG.Theoryofelasticityofananisotropicbody.Moscow:MirPublishers,1981.[103]吕爱钟.巷道开挖围岩参数及地应力可辨性的探讨.岩石力学与工程学报,1988,7(2):155-164.[104]YangZF,LeeCF,WangSJ.Three-dimensionalback-analysisofdisplacementsinexplorationedits-principlesandapplication.InternationalJournalofRockMechanicsandMiningSciences,2000,37(3):525-533.179[153]杨林德,朱合华.地层三维粘弹性反演分析.岩土工程学报,1991,13(6):18-26.陈国荣,池永斌.三维粘弹性参数反分析.河海大学学报：自然科学版,1996,24(6):25-28.陈国荣,高谦.高速公路路基性态反分析及沉降预报.工程地质学报,1998,6(4):340-343.OhkamiT,IchikawaY.Aparameteridentificationprocedureforviscoelasticmaterials.ComputersandGeotechnics,1997,21(4):255-275.[157] Ohkami T, Murai J I. Identification procedure for viscoelasticmaterials usingcorrespondenceprinciple.CommunicationsinNumericalMethodsinEngineering,1998,14(6):497-504.高强,郭杏林,杨海天.遗传算法求解粘弹性反问题.大连理工大学学报,2000,40(6):664-668.王永岩,李剑光,魏佳,等.黏弹性有限元反分析方法及其在软岩流变问题中的应用.煤炭学报,2007,32(11):1162-1165.[160]CividiniA,BarlaG.Calibrationofrheologicalmaterialmodelonthebasisoffieldmeasurements.//IchikawaY,KawamotoT,eds.Proceedingsof5thInternationalConferenceonNumericalMethodsinGeomechanics.Rotterdam:Balkema,1985:1621-1628.[161]王芝银,刘怀恒.粘弹塑性有限元分析及其在岩石力学与工程中的应用.西安矿业学院学报,1985,(1):62-73.[162]王芝银,李云鹏.地下工程围岩粘弹塑性参数反分析.水利学报,1990,(9):11-16.[163]沈振中.三维粘弹塑性位移反分析的可变容差法.水利学报,1997,(9):66-70.DziadziuszkoP.Inverseanalysisproceduresfornonlinearandtime-dependentconstitutivemodelsingeotechnics[DoctoralDissertation].Nagoya:NagoyaUniversity,2000.赵新铭,刘宁.岩体弹-粘塑性参数的Bayes随机优化反分析.南京航空航天大学学报,2008,40(1):125-128.[166]SakuraiS.Numericalanalysisfortheinterpretationoffieldmeasurementsingeomechanics.//DesaiCS,GiodaG,eds.NumericalMethodsandconstitutiveModellinginGeomechanics.NewYork:SpringerVerlag,1990:351-407.[167]SakuraiS.Fieldmeasurementandbackanalysis.//BeerG,BookerJR,CarterJP,eds.ComputerMethodsandAdvancesinGeomechanics.Rotterdam:Balkema,1991:1693-1701.[168]SakuraiS.Assessmentofcutslopestabilitybymeansofbackanalysisofmeasureddisplacements.//Pasamehmetoglu,GeunhanA,eds.AssessmentandPreventionoffailurePhenomenainrockEngineering.Rotterdam:Balkema,1993:3-9.[169]OhkamiT,IchikawaY,KawamotoT.Aboundaryelementmethodforidentifying180orthotropicmaterialparameters.InternationalJournalforNumericalandAnalyticalMethodsinGeomechanics,1991,15(9):609-625.李素华,朱维申.优化方法在弹性、横观各向同性以及弹塑性变形观测反分析中的应用.

岩石力学与工程学报

,1993,12(2):105-114.[171] Dong

Q X.

Parameter

identification

and

its

application

totunnelling[DoctoralDissertation].Innsbruck:InnsbruckUniversity,1997.SwobodaG,IchikawaY,DongQ,etal.Backanalysisoflargegeotechnicalmodels.InternationalJournalforNumericalandAnalyticalMethodsinGeomechanics,1999,23(13):1455-1472.刘学增,朱合华.考虑动态施工过程的岩土介质横观各向同性粘弹性反分析及其工程应用.岩土工程学报,2002,24(1):89-92.[174] 黄立新. 平面正交各向异性体材料参数识别的边界元法 [博士学位论文 ]. 北京:清华大学,2005.[175] 邝宏柱, 刘学增. 层状地层横观各向同性粘弹性优化反分析 . 地下空间与工程学报,2006,2(1):112-114.[176] 黄光明, 李云, 顾冲时, 等. 碾压混凝土坝横观各向同性粘弹性参数反演 .水利水运工程学报,2006,(4):15-20.[177]姜谙男,茹忠亮,张娇.基于粒子群算法和FLAC的洞室围岩参数反分析.矿业研究与开发,2007,27(5):33-35.罗润林,阮怀宁,黄亚哲,等.岩体初始地应力场的粒子群优化反演及在FLAC（-3D）中的实现.长江科学院院报,2008,25(4):73-76.[179]张志增.白云东矿山高陡边坡的稳定性分析和加固技术研究[硕士学位论文].北京:北京科技大学大学,2006.[180]GiodaG.Indirectidentificationoftheaverageelasticcharacterizationofrockmasses.//ProcIntConfonStructuralFoundationonRock.Sydney,1980:65-73.[181]GiodaG.Someremarksonbackanalysisandcharacterizationproblemsingeomechanics.//ProceedingsoftheFifthInternationalConferenceforNumericalMethodsinGeomechanics.Nagoya,1985:47-61.[182]MaierG,GiodaG.Optimizationmethodsforparametricidentificationofgeotechnicalsystem.//MaitinsJ,Bed,eds.NumericalMethodsinGeomechanics.Boston,1980:431-436.[183]GiodaG,MaierG.Directsearchsolutionofaninverseprobleminelectroplasticity:Identificationofcohesion,frictionangleandinsitustressbypressuretunneltests.InternationalJournalforNumericalMethodsinEngineering,1980,15(12):1823-1848.[184] Gioda G, Pandolfi A,some back analysisalgorithms and their application

Cividinito in-situ

A. A comparativeload tests.// Proc

evaluation of2nd Int Sympon

field1821964:281-298.[204] 王开健, 刘西拉, 顾雷. 基于 MPI机群环境下的广义逆力法并行化初探 .岩石力学与工程学报,2005,24(1):57-65.刘晓平,安竹林,郑利平.基于MPI的主从式并行遗传算法框架.系统仿真学报,2004,16(9):1938-1940.HollandJH.AdaptationinNaturalandArtificialSystem.Michigan:UniversityofMichiganPress,1975.[207] 李邵军, 冯夏庭, 廖红建. 并行遗传算法对硅藻软岩本构关系的求解 . 岩土力学,2002,23(1):1-5.WilkinsonB.,AllenM.并行程序设计.陆鑫达译.北京:机械工业出版社,2005.李仲奎,戴荣,姜逸明.FLAC3D分析中的初始应力场生成及在大型地下洞室群计算中的应用. 岩石力学与工程学报 ,2002,21(A02):2387-2392.[210] 张学民. 岩石材料各向异性特征及其对隧道围岩稳定性影响研究

[博士学位论文].

长沙:

中南大学

,2007.李仲奎,周钟,汤雪峰,等.锦屏一级水电站地下厂房洞室群稳定性分析与思考.岩石力学与工程学报,2009,28(11):2167-2175.李仲奎,徐千军,张志增,等.雅砻江锦屏一级水电站地下厂房洞室群施工期快速监测与反馈分析科研成果报告 (一)~(九). 清华大学水利水电工程系 , 2007年 4月~2008年12月.周钟,汤雪峰,廖成刚,等.雅砻江锦屏一级水电站地下厂房洞室群围岩稳定与支护、厂房主要结构设计报告 . 中国水电顾问集团成都勘测设计研究院 ,2007年4月.PeiJM,LiZK.Criterionfor3Dand2Dfailuremodesandstrengthofrockmassunderconsiderationoftheeffectofweakplanes.JournalofGeotechnicalEngineering,1988,.10(2):3.[215] Li Z K. Modeling