锦屏水电站地应力场回归分析与坝肩槽开挖高边坡稳定性三维非线性有限元分析
|
摘要
锦屏一级水电站位于四川省盐源县和木里县境内，是雅砻江干流上的重要梯级电站。电站设计装机容量3300MW，采用混凝土抛物型双曲拱坝，最大坝高305.0m，是拟建的国际上最高拱坝。无论是对大坝还是对地下厂房洞室群的稳定而言，地应力场的研究都具有非常重要的意义，同时，大坝基础开挖后形成的300m级高边坡稳定问题也是一个有待深入研究的课题。
本文根据锦屏一级水电站坝址区的地形地质特性，运用三维非线性有限元仿真计算方法，对坝址区初始地应力场和坝肩槽开挖高边坡的稳定性及工程加固措施开展了系统研究。
首先结合锦屏一级水电站坝址区的地形地质条件与地应力实测资料，运用三维有限元回归分析方法，建立工程区域的三维初始地应力场且真实反映坝址区地质构造的有限元模型，并用逐步回归方法探讨了地应力场的待回归因素与构造应力场的模拟方式，对比分析了两种回归方案的特点，为坝肩槽开挖高边坡和地下厂房洞室群围岩稳定性分析提供合理的初始地应力场。
在确定的坝址区初始地应力场条件下，采用三维非线性有限元的仿真技术，模拟左右岸坝肩槽分级开挖过程，研究了岩质边坡的卸荷效应，以及应力、位移和点安全系数的分布规律，揭示了边坡可能出现的破坏部位与发育程度。在此基础上提出了相应的加固处理措施，并研究了加固措施的工程效应，论证了工程加固措施的合理性与可靠性。
本文的特色在于：(1)针对锦屏水电站工程坝区地形地质条件，建模时充分考虑了坝址区内的断层、深部裂缝以及节理裂隙的分布特性，客观地反映了工程的实际情况；(2)在地应力场回归分析计算时，对构造应力场的模拟方
摘要
式作了探讨与新的尝试，引入剪切位移模拟地质构造作用对主应力方向的调
整作用，采用逐步回归的方法来选择待回归因素，讨论了岩体自重应力场在
初始地应力场中的权重问题，并考虑了两种地应力场回归方案，综合评价了
各自的特点;(3)在三维非线性有限元仿真计算中，考虑了岩体的Drucke二
Prager弹塑性模型和遍布节理的定向破坏模型，研究分级开挖过程中坝肩槽
边坡岩体的应力场、位移场及破坏信息分布特征;(4)通过综合分析论证，提
出了坝肩槽高边坡的加固处理方案并进行了作用效果研究;(5)首次采用隐含
锚杆有限元子结构形式，系统研究锚杆对边坡岩体的加固作用以及锚杆在随
加固区变形过程中存在的反向“锁固力”对岩体变形的抑制作用。
Jinping hydropower station is located on Yalong River in Sichuan province with installed capacity 3300 MW. The maximum height of the concrete parabolic double-curvature arch dam is 305.2 meters and it is the highest of the plans in the world. The study of initial stress field is of far-reach significance to the stability of both dam and the underground houses, at the same time, the stability of 300 meters high slope caused by the excavation of the foundation need further research.
In combination with the geological data and the measured stress data of Jinping, using three-dimension (3-D) nonlinear finite element method (FEM), initial stress field, the stability of abutment high slopes and the reinforcement measurements of the project are systematically studied.
First of all, integrated with the all data what was mentioned, built the 3D-model which can factually image the geological structure of Jinping, the regression estimation and the simulation technique are searched by regression step by step, and one of the two regression precepts which gains by contrast is to provide reasonable initial stress field for the following, the stability analysis of excavation of the high slop and surrounding rock of underground houses.
As the initial stress field is determined, the excavation process of abutment is simulated with 3D-FEM, the load-off effect of rock slop and the distribution of stress, distance and the safe coefficient are studied, and then the conceivable failure parts and develop degree of rock are opened out. Based on these, corresponding reinforcement measurements are put forward and their reasonability and reliability are demonstrated.
The features of this paper are: (1) according to the geological conditions of Jinping project, model is established fully considering the main faults, deep cracks and joint fissures that accord with actual conditions. (2) In the calculation of initial stress field, the maiden attempt that shear distance is introduced to adjust the direction of principal stress is entered for the simulation of tectonism, and the regression estimate is selected step by step, and then each character of two regression precepts is appraised by synthesis after the weight of the self-weight stress field in the initial geological stress field is discussed. (3) when the distribution of stress field, displacement field and failure information during each excavation process are studied, the nonlinear elastoplastic model named Drucker-Prager and the direct failure model in rock masses are considered in 3D-FME. (4) Through comprehensive analysis and demonstration, reinforcement measurements are put forward for the high slop and their working effectiveness is studied. (5) The sub-structure concerning implicit bolts is firstly used to research the bolts' strengthening function and their restrain for deformation.
本文根据锦屏一级水电站坝址区的地形地质特性，运用三维非线性有限元仿真计算方法，对坝址区初始地应力场和坝肩槽开挖高边坡的稳定性及工程加固措施开展了系统研究。
首先结合锦屏一级水电站坝址区的地形地质条件与地应力实测资料，运用三维有限元回归分析方法，建立工程区域的三维初始地应力场且真实反映坝址区地质构造的有限元模型，并用逐步回归方法探讨了地应力场的待回归因素与构造应力场的模拟方式，对比分析了两种回归方案的特点，为坝肩槽开挖高边坡和地下厂房洞室群围岩稳定性分析提供合理的初始地应力场。
在确定的坝址区初始地应力场条件下，采用三维非线性有限元的仿真技术，模拟左右岸坝肩槽分级开挖过程，研究了岩质边坡的卸荷效应，以及应力、位移和点安全系数的分布规律，揭示了边坡可能出现的破坏部位与发育程度。在此基础上提出了相应的加固处理措施，并研究了加固措施的工程效应，论证了工程加固措施的合理性与可靠性。
本文的特色在于：(1)针对锦屏水电站工程坝区地形地质条件，建模时充分考虑了坝址区内的断层、深部裂缝以及节理裂隙的分布特性，客观地反映了工程的实际情况；(2)在地应力场回归分析计算时，对构造应力场的模拟方
摘要
式作了探讨与新的尝试，引入剪切位移模拟地质构造作用对主应力方向的调
整作用，采用逐步回归的方法来选择待回归因素，讨论了岩体自重应力场在
初始地应力场中的权重问题，并考虑了两种地应力场回归方案，综合评价了
各自的特点;(3)在三维非线性有限元仿真计算中，考虑了岩体的Drucke二
Prager弹塑性模型和遍布节理的定向破坏模型，研究分级开挖过程中坝肩槽
边坡岩体的应力场、位移场及破坏信息分布特征;(4)通过综合分析论证，提
出了坝肩槽高边坡的加固处理方案并进行了作用效果研究;(5)首次采用隐含
锚杆有限元子结构形式，系统研究锚杆对边坡岩体的加固作用以及锚杆在随
加固区变形过程中存在的反向“锁固力”对岩体变形的抑制作用。
Jinping hydropower station is located on Yalong River in Sichuan province with installed capacity 3300 MW. The maximum height of the concrete parabolic double-curvature arch dam is 305.2 meters and it is the highest of the plans in the world. The study of initial stress field is of far-reach significance to the stability of both dam and the underground houses, at the same time, the stability of 300 meters high slope caused by the excavation of the foundation need further research.
In combination with the geological data and the measured stress data of Jinping, using three-dimension (3-D) nonlinear finite element method (FEM), initial stress field, the stability of abutment high slopes and the reinforcement measurements of the project are systematically studied.
First of all, integrated with the all data what was mentioned, built the 3D-model which can factually image the geological structure of Jinping, the regression estimation and the simulation technique are searched by regression step by step, and one of the two regression precepts which gains by contrast is to provide reasonable initial stress field for the following, the stability analysis of excavation of the high slop and surrounding rock of underground houses.
As the initial stress field is determined, the excavation process of abutment is simulated with 3D-FEM, the load-off effect of rock slop and the distribution of stress, distance and the safe coefficient are studied, and then the conceivable failure parts and develop degree of rock are opened out. Based on these, corresponding reinforcement measurements are put forward and their reasonability and reliability are demonstrated.
The features of this paper are: (1) according to the geological conditions of Jinping project, model is established fully considering the main faults, deep cracks and joint fissures that accord with actual conditions. (2) In the calculation of initial stress field, the maiden attempt that shear distance is introduced to adjust the direction of principal stress is entered for the simulation of tectonism, and the regression estimate is selected step by step, and then each character of two regression precepts is appraised by synthesis after the weight of the self-weight stress field in the initial geological stress field is discussed. (3) when the distribution of stress field, displacement field and failure information during each excavation process are studied, the nonlinear elastoplastic model named Drucker-Prager and the direct failure model in rock masses are considered in 3D-FME. (4) Through comprehensive analysis and demonstration, reinforcement measurements are put forward for the high slop and their working effectiveness is studied. (5) The sub-structure concerning implicit bolts is firstly used to research the bolts' strengthening function and their restrain for deformation.
引文
[1]沈明荣主编．岩体力学[M]．上海：同济大学出版社，1999．3
[2]周维垣主编．高等岩石力学[M]．北京：中国水利电力出版社，1989．3
[3]许文年主编．水电工程若干技术问题的理论与实践[M]．北京：中国水利水电出版社，1999．11
[4]崔政权，李宁编．边坡工程—理论与实践的最新发展[M]．北京：中国水利水电出版社，1999．12
[5]杨林德等编．岩土工程问题的反演理论与工程实践[M]．北京：科学出版社，1995．4
[6]杨林德．初始地应力位移反分析的有限元法[J]．同济大学学报，1985，(4)：69～77
[7]Feng Z. L. &lewis, R. W, Optimal Estimation of In-situ Ground Stresses from Displacement Measurements, Int. J. Numer. Ana. Methods in Geomech[J], 1987, (11): 391～408
[8]郑宏，葛修润等．关于岩土工程有限元分析中的若干问题[J]．岩土力学，1995(3)：7～12
[9]莫海鸿等．竖向线性分布初始地应力场的反演分析．第二届全国岩石力学数值计算与模型实验学术研讨会论文集[M]．上海同济大学出版社，1990，669～674
[10]肖明，刘志明．锦屏二级水电站三维地应力场反演回归分析[J]．人民长江，2000(9)，42～44
[11]杨法志．位移反分析及其应用[M]．中国岩石力学与工程学会编，岩石力学新进展，186～19．东北工学院出版社，1989
[12]邓建辉，廖勇龙．边坡初始地应力场的位移反分析初步研究[J]．贵州工业大学学报，1997(10)，102～107
[13]范景伟，何江达等．高土石坝关键技术研究勘测关键技术及其综合应用研究．地应力测试和分析新技术研究．四川大学“八五”国学科技攻关报告，1995．1
[14]韦敏才．某地下厂房计算中的岩石力学问题[J]．岩土工程学报，1986(4)
[15]张有天，胡惠昌．地应力场的趋势分析[J]．水利学报，1984(4)：31～38
[16]邵国建．岩体初始地应力场的反演回归分析[J]．水利水电科技进展，2000(10)：36～38
[17]刑建民．浅谈高边坡稳定性研究发展[J]．长沙大学学报，2000，14(2)：65～67．
[18]汪小刚，夏万仁，郭映龙等．水电工程高边坡稳定问题研究现状和发展方向[J]．水力发电，1994(5)15～18
[19]夏元友，朱瑞赓，李新平．边坡稳定性研究的综述与展望[J]．金属矿山，1995(12)
[20]胡柳青，李夕兵，温世游．边坡稳定性研究及其发展趋势[J]．矿业研究与开发，2000，20(5)：7～8
[21]邹成杰，庞声宽等．典型层状岩体边坡稳定分析与工程治理[M]．北京：水利电力出版社，1995．5
[22]张有天主编．岩石边坡的变形与稳定[M]．北京：中国水利水电出版社，1999
[23]水利部水利水电规划设计总院编．预应力锚固技术[M]．北京：中国水利水电出版社，2001．12
[24]夏元友，李梅．边坡稳定性评价方法研究及发展趋势[J]．岩石力学与工程学报，2002，21(7)：1087～1091
[25]陈清远．数值模拟法在边坡稳定分析中的应用．99力学与工程学术讨论会论文集．1999
[26]杨静，何江达等．土质边坡稳定分析中条分法与有限元法的比较[J]．四川水力发电，2003(1)：27～29
[27]周维垣，杨若琼等．清华大学“九五攻关”报告．拱坝与地基整体失稳机理及分析方法研究．1999，12
[28]油新华，李晓．国外离心试验技术在边坡工程中的应用现状与展望[J]．工程地质学报，2000年第04期
[29]狄京．数采系统在岩土工程离心模型试验中的应用[J]．矿山压力与顶板管理，1994年第1期
[30]潘亨永，胡小明，胡定．小湾左岸坝前高边坡稳定性离心模型试验研究[J]．红水河，2000，19(3)：29～33
[31]G．Gioda．有限元法在岩土力学中的应用[M]．张清，张弥译．北京：中国铁道
出版社．1983
[32]汪益敏．有限元法在边坡岩体稳定分析中的应用[J]．西安公路学院学报，1994，14(2)：13～18
[33]朱伯芳．有限单元法原理与应用(第二版)[M]．北京：中国水利水电出版社，1998
[34]Zienkiewicz O C. The Finite Element Method.3rd ed.London: McGraw-Hill, 1997
[35]Bathe K J. Finite Element Procedures in Engineering Analysis. Englewood Cliffs: Prentice-Hall, 1982
[36]T.J.R.Hughes,W.K.Liu. Implicit-explicit finite element in Transient Analysis. Stability Theory. J.Appl.Mechs.Vol.45, 1978
[37]I.M.Smith and D.V.Griffiths. Programming the Finite Element method, 2nd Edition
[38]Probert E J, Hasson O & Morgan K, Peraire J. An adaptive finite element for transient compressible flows with moving boundary. Int.J.Num.Meth. Eng.1991,32(4): 751～765
[39]Goodman,R.E., Taylor, R.L., Brekke, T.L: A Model for the Mechanics of Jointed Rock.J. ASCE. Soil. Mechanics and Foundation Division, 1968
[40]陈胜宏．三峡工程船闸边坡的弹粘塑性自适应有限元分析[J]．岩土力学，1998，19(1)：13～19
[41]高大钊主编．岩土工程的回顾与前瞻[M]．人民交通出版社，2001．6
[42]焦玉勇，葛修润．基于静态松驰法求解的三维离散单元法[J]．岩石力学与工程学报，2001，19(4)：453～458
[43]张建海，范景伟，胡定．刚体弹簧元理论及应用[M]．成都：成都科技大学出版社，1997．11
[44]张社贵，贾世军，郭怀志．岩石边坡稳定的可靠性分析[J]．岩土力学，1999，20(2)：57～61
[45]刘汉东．水电工程岩质高边坡安全度标准研究[J]．华北水利水电学院学报，1998，19(1)：40～44
[46]裴觉民．DDA在裂隙岩体边坡工程中的应用[J]．西北水电，1997年第1期
[47]王如路，孙钧．DDA理论及在三峡船闸岩提高边坡工程中的应用研究．第三届
全国岩石力学数值计算与模型实验学术研讨会论文集[M]．2000
[48]石根华．块体系统不连续变形数值分析新方法[M]．北京：科学出版社，1993
[49]冠晓东，周维垣．三维快速拉格朗日法进行小湾拱坝稳定分析[J]．水利学报，2000(9)：8～14
[50]寇晓东，周维垣，杨若琼．FLAC-3D进行三峡船闸高边坡稳定分析[J]．岩石力学与工程学报，2001，20(1)：6～10
[51]刘锦华，吕祖珩．块体理论在工程岩体稳定分析中的应用[M]．1988
[52]卓家寿，章青．不连续介质力学问题的界面元法[M]．北京：科学出版社，2000
[53]周维垣，杨若琼．流形元法及其在工程中的应用[J]．岩石力学与工程学报，1996，15(3)：211～218
[54]陈新民，罗国煜．基予经验的边坡稳定性灰色系统分析与评价[J]．岩石力学与工程学报，1999，21(5)：638～641．
[55]王亮．岩石边坡稳定性的可拓学分析[J]．河北冶金，1999年第一期
[56]刘瑞伶等．岩石边坡稳定性和Fuzzy综合评判法[J]．岩石力学与工程学报，1999，18(2)：170～175
[57]徐卫亚，邵建富．边坡工程反馈设计研究的人工神经网络方法[J]．武汉水利电力大学(宜昌)学报，1997，1 9f2)：9～11
[58]冯夏庭．智能岩石力学导论[M]．北京：科学出版社，2000．8
[59]樱井春辅．关于岩石力学领域的研究．世界隧道[J]，1999(2)：1～7
[60]国家电力公司成都勘测设计研究院编．雅砻江锦屏一级水电站坝址选择研究报告之四．工程地质条件(初稿)，2001．3
[61]何江达等．锦屏一级水电站坝址区地应力场回归及坝后左右岸边坡稳定性研究(初稿)，四川大学水电学院，2003．1
[62]张宁．岩体初始地应力场发育规律研究．浙江大学硕士学位论文[D]，2002．4
[63]陈庆宣，王维襄等．岩石力学与构造应力场分析[M]．北京：地质出版社，1998
[64]庞作会，陈文胜等．复杂仞始地应力场的反分析[J]．岩土工程学报，1998，20(4)：44～47
[65]朱伯芳．岩体初始地应力反分析[J]．水利学报，1994(10)：30～34
[66]朱焕春，陶振宇．河谷走向与河谷地应力分布[J]．岩石力学与工程学报，14(1)．17～24
[67]朱焕春，李浩．论岩体构造应力[J]．水利学报，2001(9)：81～84
[68]张奇华，钟作武等．施加边界位移产生纯剪应力及反分析应用[J]．长江科学院院报，2000，17(2)：34～36
[69]刘永川等．乌江索风营水电站施工阶段右岸地下厂房洞室群围岩稳定三维有限元分析(中间报告)．国家电力公司贵州勘测设计院等，2003．2
[70]张林、何江达等．含诱导缝碾压混凝土拱坝开裂和破坏机制研究(“九五”国家重点科技攻关项目)．四川大学水电学院，2000．1
[71]雷晓燕．三维锚杆单元理论及其应用[J]．工程力学，1996，13(2)：50～60
[72]ZHAO De-an. A Consideration of Plasticity in a Rock-Bolt Element[J]. JOURNAL OF LANZHOU RAILWAY UNIVERSITY(Natural Sciences), 20(1): 12～17
[2]周维垣主编．高等岩石力学[M]．北京：中国水利电力出版社，1989．3
[3]许文年主编．水电工程若干技术问题的理论与实践[M]．北京：中国水利水电出版社，1999．11
[4]崔政权，李宁编．边坡工程—理论与实践的最新发展[M]．北京：中国水利水电出版社，1999．12
[5]杨林德等编．岩土工程问题的反演理论与工程实践[M]．北京：科学出版社，1995．4
[6]杨林德．初始地应力位移反分析的有限元法[J]．同济大学学报，1985，(4)：69～77
[7]Feng Z. L. &lewis, R. W, Optimal Estimation of In-situ Ground Stresses from Displacement Measurements, Int. J. Numer. Ana. Methods in Geomech[J], 1987, (11): 391～408
[8]郑宏，葛修润等．关于岩土工程有限元分析中的若干问题[J]．岩土力学，1995(3)：7～12
[9]莫海鸿等．竖向线性分布初始地应力场的反演分析．第二届全国岩石力学数值计算与模型实验学术研讨会论文集[M]．上海同济大学出版社，1990，669～674
[10]肖明，刘志明．锦屏二级水电站三维地应力场反演回归分析[J]．人民长江，2000(9)，42～44
[11]杨法志．位移反分析及其应用[M]．中国岩石力学与工程学会编，岩石力学新进展，186～19．东北工学院出版社，1989
[12]邓建辉，廖勇龙．边坡初始地应力场的位移反分析初步研究[J]．贵州工业大学学报，1997(10)，102～107
[13]范景伟，何江达等．高土石坝关键技术研究勘测关键技术及其综合应用研究．地应力测试和分析新技术研究．四川大学“八五”国学科技攻关报告，1995．1
[14]韦敏才．某地下厂房计算中的岩石力学问题[J]．岩土工程学报，1986(4)
[15]张有天，胡惠昌．地应力场的趋势分析[J]．水利学报，1984(4)：31～38
[16]邵国建．岩体初始地应力场的反演回归分析[J]．水利水电科技进展，2000(10)：36～38
[17]刑建民．浅谈高边坡稳定性研究发展[J]．长沙大学学报，2000，14(2)：65～67．
[18]汪小刚，夏万仁，郭映龙等．水电工程高边坡稳定问题研究现状和发展方向[J]．水力发电，1994(5)15～18
[19]夏元友，朱瑞赓，李新平．边坡稳定性研究的综述与展望[J]．金属矿山，1995(12)
[20]胡柳青，李夕兵，温世游．边坡稳定性研究及其发展趋势[J]．矿业研究与开发，2000，20(5)：7～8
[21]邹成杰，庞声宽等．典型层状岩体边坡稳定分析与工程治理[M]．北京：水利电力出版社，1995．5
[22]张有天主编．岩石边坡的变形与稳定[M]．北京：中国水利水电出版社，1999
[23]水利部水利水电规划设计总院编．预应力锚固技术[M]．北京：中国水利水电出版社，2001．12
[24]夏元友，李梅．边坡稳定性评价方法研究及发展趋势[J]．岩石力学与工程学报，2002，21(7)：1087～1091
[25]陈清远．数值模拟法在边坡稳定分析中的应用．99力学与工程学术讨论会论文集．1999
[26]杨静，何江达等．土质边坡稳定分析中条分法与有限元法的比较[J]．四川水力发电，2003(1)：27～29
[27]周维垣，杨若琼等．清华大学“九五攻关”报告．拱坝与地基整体失稳机理及分析方法研究．1999，12
[28]油新华，李晓．国外离心试验技术在边坡工程中的应用现状与展望[J]．工程地质学报，2000年第04期
[29]狄京．数采系统在岩土工程离心模型试验中的应用[J]．矿山压力与顶板管理，1994年第1期
[30]潘亨永，胡小明，胡定．小湾左岸坝前高边坡稳定性离心模型试验研究[J]．红水河，2000，19(3)：29～33
[31]G．Gioda．有限元法在岩土力学中的应用[M]．张清，张弥译．北京：中国铁道
出版社．1983
[32]汪益敏．有限元法在边坡岩体稳定分析中的应用[J]．西安公路学院学报，1994，14(2)：13～18
[33]朱伯芳．有限单元法原理与应用(第二版)[M]．北京：中国水利水电出版社，1998
[34]Zienkiewicz O C. The Finite Element Method.3rd ed.London: McGraw-Hill, 1997
[35]Bathe K J. Finite Element Procedures in Engineering Analysis. Englewood Cliffs: Prentice-Hall, 1982
[36]T.J.R.Hughes,W.K.Liu. Implicit-explicit finite element in Transient Analysis. Stability Theory. J.Appl.Mechs.Vol.45, 1978
[37]I.M.Smith and D.V.Griffiths. Programming the Finite Element method, 2nd Edition
[38]Probert E J, Hasson O & Morgan K, Peraire J. An adaptive finite element for transient compressible flows with moving boundary. Int.J.Num.Meth. Eng.1991,32(4): 751～765
[39]Goodman,R.E., Taylor, R.L., Brekke, T.L: A Model for the Mechanics of Jointed Rock.J. ASCE. Soil. Mechanics and Foundation Division, 1968
[40]陈胜宏．三峡工程船闸边坡的弹粘塑性自适应有限元分析[J]．岩土力学，1998，19(1)：13～19
[41]高大钊主编．岩土工程的回顾与前瞻[M]．人民交通出版社，2001．6
[42]焦玉勇，葛修润．基于静态松驰法求解的三维离散单元法[J]．岩石力学与工程学报，2001，19(4)：453～458
[43]张建海，范景伟，胡定．刚体弹簧元理论及应用[M]．成都：成都科技大学出版社，1997．11
[44]张社贵，贾世军，郭怀志．岩石边坡稳定的可靠性分析[J]．岩土力学，1999，20(2)：57～61
[45]刘汉东．水电工程岩质高边坡安全度标准研究[J]．华北水利水电学院学报，1998，19(1)：40～44
[46]裴觉民．DDA在裂隙岩体边坡工程中的应用[J]．西北水电，1997年第1期
[47]王如路，孙钧．DDA理论及在三峡船闸岩提高边坡工程中的应用研究．第三届
全国岩石力学数值计算与模型实验学术研讨会论文集[M]．2000
[48]石根华．块体系统不连续变形数值分析新方法[M]．北京：科学出版社，1993
[49]冠晓东，周维垣．三维快速拉格朗日法进行小湾拱坝稳定分析[J]．水利学报，2000(9)：8～14
[50]寇晓东，周维垣，杨若琼．FLAC-3D进行三峡船闸高边坡稳定分析[J]．岩石力学与工程学报，2001，20(1)：6～10
[51]刘锦华，吕祖珩．块体理论在工程岩体稳定分析中的应用[M]．1988
[52]卓家寿，章青．不连续介质力学问题的界面元法[M]．北京：科学出版社，2000
[53]周维垣，杨若琼．流形元法及其在工程中的应用[J]．岩石力学与工程学报，1996，15(3)：211～218
[54]陈新民，罗国煜．基予经验的边坡稳定性灰色系统分析与评价[J]．岩石力学与工程学报，1999，21(5)：638～641．
[55]王亮．岩石边坡稳定性的可拓学分析[J]．河北冶金，1999年第一期
[56]刘瑞伶等．岩石边坡稳定性和Fuzzy综合评判法[J]．岩石力学与工程学报，1999，18(2)：170～175
[57]徐卫亚，邵建富．边坡工程反馈设计研究的人工神经网络方法[J]．武汉水利电力大学(宜昌)学报，1997，1 9f2)：9～11
[58]冯夏庭．智能岩石力学导论[M]．北京：科学出版社，2000．8
[59]樱井春辅．关于岩石力学领域的研究．世界隧道[J]，1999(2)：1～7
[60]国家电力公司成都勘测设计研究院编．雅砻江锦屏一级水电站坝址选择研究报告之四．工程地质条件(初稿)，2001．3
[61]何江达等．锦屏一级水电站坝址区地应力场回归及坝后左右岸边坡稳定性研究(初稿)，四川大学水电学院，2003．1
[62]张宁．岩体初始地应力场发育规律研究．浙江大学硕士学位论文[D]，2002．4
[63]陈庆宣，王维襄等．岩石力学与构造应力场分析[M]．北京：地质出版社，1998
[64]庞作会，陈文胜等．复杂仞始地应力场的反分析[J]．岩土工程学报，1998，20(4)：44～47
[65]朱伯芳．岩体初始地应力反分析[J]．水利学报，1994(10)：30～34
[66]朱焕春，陶振宇．河谷走向与河谷地应力分布[J]．岩石力学与工程学报，14(1)．17～24
[67]朱焕春，李浩．论岩体构造应力[J]．水利学报，2001(9)：81～84
[68]张奇华，钟作武等．施加边界位移产生纯剪应力及反分析应用[J]．长江科学院院报，2000，17(2)：34～36
[69]刘永川等．乌江索风营水电站施工阶段右岸地下厂房洞室群围岩稳定三维有限元分析(中间报告)．国家电力公司贵州勘测设计院等，2003．2
[70]张林、何江达等．含诱导缝碾压混凝土拱坝开裂和破坏机制研究(“九五”国家重点科技攻关项目)．四川大学水电学院，2000．1
[71]雷晓燕．三维锚杆单元理论及其应用[J]．工程力学，1996，13(2)：50～60
[72]ZHAO De-an. A Consideration of Plasticity in a Rock-Bolt Element[J]. JOURNAL OF LANZHOU RAILWAY UNIVERSITY(Natural Sciences), 20(1): 12～17