瀑布沟水电站泄洪洞出口段开挖变形及稳定性非线性有限元分析
|
摘要
我国大型水电工程多建在高山峡谷中，具有水头高、流量大，泄洪功率巨大等显著特点。泄洪洞下泄流量大、洞内流速高，能否安全运行，将直接影响土石坝的安全，也是整个水利枢纽能否正常安全运行的关键。瀑布沟水电站位于四川省汉源县境内大渡河中游，装机总容量3300MW，总库容50.6亿m~3。工程以发电为主，兼有防洪、拦沙等综合效益，拦河大坝为直心墙土石坝，最大坝高186.0m。泄洪洞布置于左岸地下厂房的花岗岩体内，洞线顺直，隧洞为无压明流洞，断面形式为圆拱直墙型。
瀑布沟泄洪洞出口段右侧壁最小厚度仅约10.0m，左侧山体高达100.0m，隧洞偏压问题十分突出。隧洞偏压对边坡稳定性、隧洞出口段衬砌的影响以及偏压所致变形不对称性是工程非常关心的问题。考虑到岩体的不连续性、非均匀性、应力—应变关系的非线性特性以及力学性质的各向异性，本文采用非线性有限元法，计算围岩的变形和应力，定量评价岩体的稳定性，并追踪岩体破坏状态的发展演变过程。
非线性有限元法在处理复杂结构、复杂边界及荷载条件、非线性问题方面具有独特的效能，并能获得更为详尽、准确的结果，现已成为隧洞计算最有效的方法之一。本文根据泄洪洞总体施工方案和隧洞出口的地质条件，建模时严格模拟岩层界面、地形、F_1断层结构面的产状，以及泄洪洞体形和各种工程处理措施，充分反映了斜坡地形地质条件对隧洞的影响。在有限元计算中将泄洪洞分三级开挖用以模拟施工开挖过程，采用遍在节理包络强度模型模拟岩体材料，运用点强度储备安全系数法研究隧洞开挖过程中边坡的稳定性和隧洞的偏压应力。本文的主要工作：(1)基于工程的地质CAD剖面图
四川大学硕士学位论文
形，建立二维有限元网格图，利用AutoLIsP和Fortran语言开发二维有限元
前处理优化程序，实现二维有限元网格节点的自动编码、点差优化排序，和
网格单元编码和材料信息的自动生成，使复杂工程中有限元计算前处理的功
效得到较大的提高。(2)分析计算了施工开挖前围岩的应力情况和自然边坡稳
定性;(3)在边坡进行支护和不支护、分级开挖无衬砌和有衬砌等各种工况下，
分别计算围岩、边坡稳定性和结构安全性，并对支护方案作出评价;(4)在运
行、检修和地震工况下，对比分析衬砌厚为1.5m、2.0m时围岩稳定性、结
构安全性，推荐衬砌厚度方案。(5)将W类围岩指标改为V类围岩指标作力学
参数敏感性分析。
通过综合分析，建议在泄洪洞开挖前，先期进行岸坡加固，同时，在泄
洪洞分级(步)开挖过程中，及时作好一次喷锚支护，在上半洞喷锚支护完
成条件下进行下半洞开挖，有条件时尽可能及早进行二次硅衬砌。本文研究
成果已得到工程应用。
Many large scale hydropower stations of our country are built in gorges among high mountains. High hydraulic head, large-discharge, huge power of flood discharge are notable characteristics of these water powers. The operational safety of earth-rock dam, or even the whole hydraulic project, is directly influenced by discharge tunnel. Pubugou hydropower station, located in Hanyuan county, Sichuan province, has 3300MW power capacity, and 5060Mm3 total storage capacity. The project is built for electricity generation, flood control, sand barrage, etc. The maximal height of the vertical core-walled rockfill dam is 186m. The free-flow discharge tunnel, whose form of section is circular arch and straight wall, is located in granite rock below underground workshop of left hillside.
Because the minimal thickness of the right side of tunnel is only 10m and the height of the left side of the tunnel is over 100m, the surrounding pressure distribution at the discharge tunnel of Pubugou Hydropower Station is much unbalanced. The unbalanced pressure of the tunnel and the stability of the slope during the course of excavation are studied using nonlinear finite element method (FEM). The nonlinear model strictly simulates the overall geological condition, such as fault and fissures, the tunnel's section form and various schemes of construction. The stability of the slope, under deformation of excavation and the unbalanced stress of the tunnel, is evaluated and appraised by point safety factor distribution.
The major works of this thesis are as follows: (1) Using AutoLISP and
Fortran, a 2D mesh generation program is developed to achieve automatic node coding and obtaining the information of element and corresponding material.(2) The unbalanced stress state and stability of the slope before excavation are calculated and analyzed. (3)Stability analysis of surrounding rock and slope, and appraisal of structure safety are carried out, under the conditions of support or non-support, with lining or without lining. (4)The stress state of 1.5m thick lining is compared with that of 2.0m thick lining. (5)By degrading rank IV surrounding rock to rank V rock, the sensitivity of slope stability to mechanical parameter is also checked.
It's recommended to reinforce the slope before excavation, to bolt and to shotcrete in time during excavation, and to excavate the lower-half of tunnel after bolting and shotcreting of the upper-half, and to make the second shotcrete lining as early as possible. These suggestions have been adopted in practice.
瀑布沟泄洪洞出口段右侧壁最小厚度仅约10.0m，左侧山体高达100.0m，隧洞偏压问题十分突出。隧洞偏压对边坡稳定性、隧洞出口段衬砌的影响以及偏压所致变形不对称性是工程非常关心的问题。考虑到岩体的不连续性、非均匀性、应力—应变关系的非线性特性以及力学性质的各向异性，本文采用非线性有限元法，计算围岩的变形和应力，定量评价岩体的稳定性，并追踪岩体破坏状态的发展演变过程。
非线性有限元法在处理复杂结构、复杂边界及荷载条件、非线性问题方面具有独特的效能，并能获得更为详尽、准确的结果，现已成为隧洞计算最有效的方法之一。本文根据泄洪洞总体施工方案和隧洞出口的地质条件，建模时严格模拟岩层界面、地形、F_1断层结构面的产状，以及泄洪洞体形和各种工程处理措施，充分反映了斜坡地形地质条件对隧洞的影响。在有限元计算中将泄洪洞分三级开挖用以模拟施工开挖过程，采用遍在节理包络强度模型模拟岩体材料，运用点强度储备安全系数法研究隧洞开挖过程中边坡的稳定性和隧洞的偏压应力。本文的主要工作：(1)基于工程的地质CAD剖面图
四川大学硕士学位论文
形，建立二维有限元网格图，利用AutoLIsP和Fortran语言开发二维有限元
前处理优化程序，实现二维有限元网格节点的自动编码、点差优化排序，和
网格单元编码和材料信息的自动生成，使复杂工程中有限元计算前处理的功
效得到较大的提高。(2)分析计算了施工开挖前围岩的应力情况和自然边坡稳
定性;(3)在边坡进行支护和不支护、分级开挖无衬砌和有衬砌等各种工况下，
分别计算围岩、边坡稳定性和结构安全性，并对支护方案作出评价;(4)在运
行、检修和地震工况下，对比分析衬砌厚为1.5m、2.0m时围岩稳定性、结
构安全性，推荐衬砌厚度方案。(5)将W类围岩指标改为V类围岩指标作力学
参数敏感性分析。
通过综合分析，建议在泄洪洞开挖前，先期进行岸坡加固，同时，在泄
洪洞分级(步)开挖过程中，及时作好一次喷锚支护，在上半洞喷锚支护完
成条件下进行下半洞开挖，有条件时尽可能及早进行二次硅衬砌。本文研究
成果已得到工程应用。
Many large scale hydropower stations of our country are built in gorges among high mountains. High hydraulic head, large-discharge, huge power of flood discharge are notable characteristics of these water powers. The operational safety of earth-rock dam, or even the whole hydraulic project, is directly influenced by discharge tunnel. Pubugou hydropower station, located in Hanyuan county, Sichuan province, has 3300MW power capacity, and 5060Mm3 total storage capacity. The project is built for electricity generation, flood control, sand barrage, etc. The maximal height of the vertical core-walled rockfill dam is 186m. The free-flow discharge tunnel, whose form of section is circular arch and straight wall, is located in granite rock below underground workshop of left hillside.
Because the minimal thickness of the right side of tunnel is only 10m and the height of the left side of the tunnel is over 100m, the surrounding pressure distribution at the discharge tunnel of Pubugou Hydropower Station is much unbalanced. The unbalanced pressure of the tunnel and the stability of the slope during the course of excavation are studied using nonlinear finite element method (FEM). The nonlinear model strictly simulates the overall geological condition, such as fault and fissures, the tunnel's section form and various schemes of construction. The stability of the slope, under deformation of excavation and the unbalanced stress of the tunnel, is evaluated and appraised by point safety factor distribution.
The major works of this thesis are as follows: (1) Using AutoLISP and
Fortran, a 2D mesh generation program is developed to achieve automatic node coding and obtaining the information of element and corresponding material.(2) The unbalanced stress state and stability of the slope before excavation are calculated and analyzed. (3)Stability analysis of surrounding rock and slope, and appraisal of structure safety are carried out, under the conditions of support or non-support, with lining or without lining. (4)The stress state of 1.5m thick lining is compared with that of 2.0m thick lining. (5)By degrading rank IV surrounding rock to rank V rock, the sensitivity of slope stability to mechanical parameter is also checked.
It's recommended to reinforce the slope before excavation, to bolt and to shotcrete in time during excavation, and to excavate the lower-half of tunnel after bolting and shotcreting of the upper-half, and to make the second shotcrete lining as early as possible. These suggestions have been adopted in practice.
引文
[1] 潘家铮．水工隧洞和调压室(水工隧洞部分) 水工建筑物设计丛书．水利电力出版社．1990．6
[2] 汪胡桢．水工隧洞的设计理论和计算．北京：水利电力出版社，1977，9
[3] 潘家铮，何景主编．中国大坝50年．北京：中国水利水电出版社，2000．8
[4] 潘家铮，何璟主编．中国水力发电工程(工程地质卷)．北京：中国电力出版社，2000，8
[5] 潘家铮，何璟主编．中国水力发电工程(水工卷)中国电力出版社．北京，．2000．8
[6] 周维恒，高等岩石力学．北京：水利水电出版社，1989
[7] 郑颍人，龚晓南．岩土塑性力学基础．北京：建筑工程出版社，1988
[8] 姜弘道等．水工结构工程与岩土工程的现代计算方法及程序．南京：河海大学出版社，1992．
[9] 岩石力学进展与工程应用译文集．水利水电岩石力学情报网编译．科学出版社．1987．12
[10] 张忠亭，陶振宇．洞室围岩收敛变形的时间效应特性研究．武汉水利电力学院学报．1992，8
[11] 邓聚龙．灰色系统理论教程．华中理工大学出版社．1990
[12] 张忠庭．洞室困岩收敛变形灰色理论预测．四川水利发电．1995，12．
[13] 朱伯芳．有限单元法原理与应用(第二版)．北京：中国水利水电出版社，1998
[14] 殷有泉．固体力学非线性有限元引论．北京：北京大学出版社，1987
[15] G．哥德赫．有限元法在岩土力学中的应用．张清，张弥译．北京：中国铁道出版社，1983
[16] Zienkiewicz O C. The Fininte Element Method.3rd ed. London: McGraw-Hill,1997
[17] 张学言编著．岩土塑性力学．北京：人民交通出版社，1993．9
[18] 冯夏庭．智能岩石力学导论．[M]，北京：科学出版社，2000，8
[19] 王泳嘉．离散单元法—一种适用节理岩石力学分析的数值方法．第一届全国岩石力学数值计算机模型试验讨论文集．
[20]符文喜．地应力环境场下岩体的变形特性及预测研究．成都理工学院博士论文．2000．5
[21]全国青年岩上力学与工程会议论文集．南京：1998
[22]陆培炎，史永胜主编．岩土力学数值分析与解析方法—第六届全国岩土力学数值分析方法讨论会论文集．广州：广东科技出版社，1998，7
[23]张建海，范景伟，胡定．刚体弹簧元理论及应用．成都：成都科技大学出版社，1997
[24]张建海，何江达．小湾高拱坝三维有限元弹塑性分析．云南水力发电．2000．16(1)：61～64．
[25]范景伟，何江达．含定向闭合断续节理岩体的强度特性．岩土力学与工程学报．1992．11(2)．
[26]迟世春，胡志强，林皋．丰满大坝静动力特性分析．世界地震工程．2001(9)：91～97
[27]许百立、刘世煌．地下工程设计理论与实践．地下工程技术．1999．3：1-12
[28]李宁、张平．地下洞室设计计算中的几个误区．地下工程技术．1999．4：
[29] Zhang Jianhai, et. al: Static and dynamic stability assessment of rock slopes and dam foundations using rigid body-spring element method. Int. J. Rock Mechs& Mining Scis. 2001. 38 (8): 1081-1090.
[30]李瑞遐编著：有限元法与边界元法，上海：科技教育出版社，1991
[31]安跃洲．地基与上部结构相互作用的边界元理论分析．福建建筑．2002．4
[32]朱济祥，薛玺成．岩体工程问题的数值分析方法．水利水电工程设计．1995．2
[33]钱康．天生桥一级水电站导流隧洞进口边坡稳定性分析．红水河．Vol．14，No．4．1995
[34]祁文刚．二滩水电站泄洪洞典型围岩稳定性分析．Vol．18 No．6．2002
[35]张大成，张跃民，陈晴，高永辉，宋新峰．大朝山水电站地下洞室群开挖支护．水力发电．2001(12)
[36]甘三才，王颂，王雪波．江口水电站地下厂房围岩稳定性研究．人民长江．Vol．32，No．3，Mar．2001
[37]许文年主编．水电工程若干技术问题的理论与实践．北京：中国水利水电出版社，1999．11
[38]范景伟．岩石力学论文集．成都：成都科技大学出版社，1991
[39]陶振宇．试论岩石力学的最新进展．力学进展．1999(2)
[40]李世辉．隧道围岩稳定系统分析．1991
[41]冯玉国．灰色理论在地下工程位移稳定性分类中的应用．岩土工程学报，1996．3
[42]殷有泉．考虑损伤的节理本构模型．工程地质学报．1994．2(4)
[43]朱维申，张强勇，李术才．三维弹塑性断裂损伤模型在裂隙岩体工程中的应用[J]．同体力学学报，1999，20(2)：164-170．
[44]水利电力部水利水电规划设计院主编．水利水电工程地质手册．北京：水利电力出版社，1985
[45]水利水电规划设计总院，水利水电地下建筑物情报网编．水利水电工程地下建筑物设计手册．成都：四川科学技术出版社，1993
[46]99’岩上力学与工程学术讨论会议文集．重庆：《地下空间》编辑部，1999．
[47]李文纲．鲁布革水电站地下厂房位置选择及围岩稳定性评价．四川水力发电．1986(2)
[48]周维垣，杨若琼，杨强．岩石力学发展的最新趋势．岩石力学与工程学报．2000．19(1)
[49] Kumsar H., Aydan O., Ulusay. R. : D ynamic and static stability assessment of rock slope against wedge failures. Rock Mechanics and Rock Engineering. 2000. 9. 33 (1), 31-51.
[50]吴刚，孙钧．卸荷应力状态下裂隙岩体的变形和强度特性．岩石力学与工程学报．1998，17(16)
[51]陈卫念，朱维中，李术才．节理岩体断裂损伤耦合的流变模型及其应用[J]．水利学报，1999，(12)：33—37．
[52]电力工业部东北勘测设计院主编．地下洞室的锚喷支护．北京：中国水利水电出版社
[53]中国岩石力学与工程学会，新世纪岩石力学与工程的开拓和发展，北京：中国科学技术出版社，2000．10
[54]滕丽，张建海等．利用AutoLISP语言实现有限元网格节点自动编码与坐标拾取．CAD／CAM与制造业信息化．2002．11
[55]陈道洁，王德庆．AutoCAD及其应用。成都：成都科技大学出版社．1988
[56] 王璋奇，安利强．基于SuperSAP的有限元前处理接口程序设计，华北电力大学学报，1999．Vol．26 NO．1：74-77
[57] 熊铁虹，顾伯达．ADINA程序前处理系统AED程序设计及实现，武汉工学院学报，1990，Vol．12 NO．4：24-28
[58] 何发详，刘浩吾．有限元法前处理技术的发展，水利水电科技进步，2000，12(6)：15-18
[59] 谭也平，有限元前处理系统的图形集成实现，计算机辅助工程，2000，Mar．NO．1：53-58
[60] 罗特军，罗季军，汪榴．有限元网格优化方法，四川联合大学学报(工程科学版)，1999，Vol．3 NO．3：65-72
[61] 袁明道，肖明．有限元网格半自动剖分及显示，武汉水利电力大学学报，1999，Vol．32 NO．3：72-75
[62] 张广建，黄由玲．用CAGD进行有限元前处理，工程力学，1995，Vol．12 NO．4：124-130
[63] 安晓卫，邵伟平，用于二维和三维问题的有限元前处理程序机械设计与制造，1995，(1)，P21-23
[64] 梁雪春，崔洪斌，吴义忠，曹康．AutoLISP实用教程．北京：人民邮电出版社．1999
[65] 李水乡，刘剑飞，袁明武．Auto CAD平台上的有限元前处理软件Auto MESH，计算力学学报，1999，Vol．16 NO．3：375-378
[66] K.Holc, Finite element mesh generation methods, a review and classification,Computer-aided Design, 1998, 20(1),27-32
[67] 水工建筑物抗震设计规范．DL5073-2000．北京：中国电力出版社，2001．
[68] 沈聚敏，周锡元，高小旺等．抗震工程学．北京：中国建筑工业出版社，2000．
[69] 后荣昌，韩渭宾主编．四川活动断裂与地震．北京：地震出版社，1993
[70] 李善邦．中国地震．北京：地震出版社，1981
[71] 陈厚群等．水电工程抗震设防概率水准和地震作用的概率模型．自然灾害学报．1993(2)
[72] 李桂青．结构动力可靠性理论及其应用[M]．北京：地震出版社，1993．
[73] 严松宏．结构地震响应分析的广义脉冲函数法．兰州铁道学院学报(自然科学版)．2002，Vol．21，No．6 Dec．
[74] 王进廷，杜修力，张楚汉．瑞利阻尼介质有限元离散模型．动力分析的数值稳定性．2002，Dec．Vol．22 No．6
[75] 李小军，廖振鹏，杜修力．有阻尼体系动力问题的一种显式解法[J]．地震工程与工程振动，1992，12(4)：74-79
[76] 王进廷，杜修力．有阻尼体系动力分析的一种显式差分法[J]．工程力学．2002，19(3)：109～112．
[77] 杜修力，王进廷．阻尼弹性结构动力计算的显式差分法[J]．工程力学，2002，17(5)：37～43．
[78] 寥振鹏，刘晶波．离散网格中的弹性波动(Ⅰ)[J]．地震工程与工程振动，1986，6(2)：1～15．
[79] Zhang Chuhan、Q. A. Pekau and Jin Fene. Application of FE—BE—IBE coupling to dynamic inteaction between alluvial soil and rock canyons[J]. Earthquake Engineering and Structural Dynamics, 1992, 21 (5): 367—385.
[80] 刘纪陆．求解材料非线性系统动力响应的新方法．建筑结构学报．2002，Vol．23，No．4，Aug．
[81] 王岚，陈中林．反应谱的数值计算及其应用中的误差问题．地壳形变与地震，1993，13(3)：85—89
[82] Ju JW. On energy based coupled elastoplastic damage theories: constitutive modeling and computational aspects[J] . Int. J. Solids Structure, 1989, 25(7): 803-833
[83] 苏克忠等．大坝强震安全监测．北京：中国水利水电出版社，1996
[84] 潘吕实，张弥，吴鸿庆．隧道力学数值方法[M]．北京：中国铁道出版社，1995
[2] 汪胡桢．水工隧洞的设计理论和计算．北京：水利电力出版社，1977，9
[3] 潘家铮，何景主编．中国大坝50年．北京：中国水利水电出版社，2000．8
[4] 潘家铮，何璟主编．中国水力发电工程(工程地质卷)．北京：中国电力出版社，2000，8
[5] 潘家铮，何璟主编．中国水力发电工程(水工卷)中国电力出版社．北京，．2000．8
[6] 周维恒，高等岩石力学．北京：水利水电出版社，1989
[7] 郑颍人，龚晓南．岩土塑性力学基础．北京：建筑工程出版社，1988
[8] 姜弘道等．水工结构工程与岩土工程的现代计算方法及程序．南京：河海大学出版社，1992．
[9] 岩石力学进展与工程应用译文集．水利水电岩石力学情报网编译．科学出版社．1987．12
[10] 张忠亭，陶振宇．洞室围岩收敛变形的时间效应特性研究．武汉水利电力学院学报．1992，8
[11] 邓聚龙．灰色系统理论教程．华中理工大学出版社．1990
[12] 张忠庭．洞室困岩收敛变形灰色理论预测．四川水利发电．1995，12．
[13] 朱伯芳．有限单元法原理与应用(第二版)．北京：中国水利水电出版社，1998
[14] 殷有泉．固体力学非线性有限元引论．北京：北京大学出版社，1987
[15] G．哥德赫．有限元法在岩土力学中的应用．张清，张弥译．北京：中国铁道出版社，1983
[16] Zienkiewicz O C. The Fininte Element Method.3rd ed. London: McGraw-Hill,1997
[17] 张学言编著．岩土塑性力学．北京：人民交通出版社，1993．9
[18] 冯夏庭．智能岩石力学导论．[M]，北京：科学出版社，2000，8
[19] 王泳嘉．离散单元法—一种适用节理岩石力学分析的数值方法．第一届全国岩石力学数值计算机模型试验讨论文集．
[20]符文喜．地应力环境场下岩体的变形特性及预测研究．成都理工学院博士论文．2000．5
[21]全国青年岩上力学与工程会议论文集．南京：1998
[22]陆培炎，史永胜主编．岩土力学数值分析与解析方法—第六届全国岩土力学数值分析方法讨论会论文集．广州：广东科技出版社，1998，7
[23]张建海，范景伟，胡定．刚体弹簧元理论及应用．成都：成都科技大学出版社，1997
[24]张建海，何江达．小湾高拱坝三维有限元弹塑性分析．云南水力发电．2000．16(1)：61～64．
[25]范景伟，何江达．含定向闭合断续节理岩体的强度特性．岩土力学与工程学报．1992．11(2)．
[26]迟世春，胡志强，林皋．丰满大坝静动力特性分析．世界地震工程．2001(9)：91～97
[27]许百立、刘世煌．地下工程设计理论与实践．地下工程技术．1999．3：1-12
[28]李宁、张平．地下洞室设计计算中的几个误区．地下工程技术．1999．4：
[29] Zhang Jianhai, et. al: Static and dynamic stability assessment of rock slopes and dam foundations using rigid body-spring element method. Int. J. Rock Mechs& Mining Scis. 2001. 38 (8): 1081-1090.
[30]李瑞遐编著：有限元法与边界元法，上海：科技教育出版社，1991
[31]安跃洲．地基与上部结构相互作用的边界元理论分析．福建建筑．2002．4
[32]朱济祥，薛玺成．岩体工程问题的数值分析方法．水利水电工程设计．1995．2
[33]钱康．天生桥一级水电站导流隧洞进口边坡稳定性分析．红水河．Vol．14，No．4．1995
[34]祁文刚．二滩水电站泄洪洞典型围岩稳定性分析．Vol．18 No．6．2002
[35]张大成，张跃民，陈晴，高永辉，宋新峰．大朝山水电站地下洞室群开挖支护．水力发电．2001(12)
[36]甘三才，王颂，王雪波．江口水电站地下厂房围岩稳定性研究．人民长江．Vol．32，No．3，Mar．2001
[37]许文年主编．水电工程若干技术问题的理论与实践．北京：中国水利水电出版社，1999．11
[38]范景伟．岩石力学论文集．成都：成都科技大学出版社，1991
[39]陶振宇．试论岩石力学的最新进展．力学进展．1999(2)
[40]李世辉．隧道围岩稳定系统分析．1991
[41]冯玉国．灰色理论在地下工程位移稳定性分类中的应用．岩土工程学报，1996．3
[42]殷有泉．考虑损伤的节理本构模型．工程地质学报．1994．2(4)
[43]朱维申，张强勇，李术才．三维弹塑性断裂损伤模型在裂隙岩体工程中的应用[J]．同体力学学报，1999，20(2)：164-170．
[44]水利电力部水利水电规划设计院主编．水利水电工程地质手册．北京：水利电力出版社，1985
[45]水利水电规划设计总院，水利水电地下建筑物情报网编．水利水电工程地下建筑物设计手册．成都：四川科学技术出版社，1993
[46]99’岩上力学与工程学术讨论会议文集．重庆：《地下空间》编辑部，1999．
[47]李文纲．鲁布革水电站地下厂房位置选择及围岩稳定性评价．四川水力发电．1986(2)
[48]周维垣，杨若琼，杨强．岩石力学发展的最新趋势．岩石力学与工程学报．2000．19(1)
[49] Kumsar H., Aydan O., Ulusay. R. : D ynamic and static stability assessment of rock slope against wedge failures. Rock Mechanics and Rock Engineering. 2000. 9. 33 (1), 31-51.
[50]吴刚，孙钧．卸荷应力状态下裂隙岩体的变形和强度特性．岩石力学与工程学报．1998，17(16)
[51]陈卫念，朱维中，李术才．节理岩体断裂损伤耦合的流变模型及其应用[J]．水利学报，1999，(12)：33—37．
[52]电力工业部东北勘测设计院主编．地下洞室的锚喷支护．北京：中国水利水电出版社
[53]中国岩石力学与工程学会，新世纪岩石力学与工程的开拓和发展，北京：中国科学技术出版社，2000．10
[54]滕丽，张建海等．利用AutoLISP语言实现有限元网格节点自动编码与坐标拾取．CAD／CAM与制造业信息化．2002．11
[55]陈道洁，王德庆．AutoCAD及其应用。成都：成都科技大学出版社．1988
[56] 王璋奇，安利强．基于SuperSAP的有限元前处理接口程序设计，华北电力大学学报，1999．Vol．26 NO．1：74-77
[57] 熊铁虹，顾伯达．ADINA程序前处理系统AED程序设计及实现，武汉工学院学报，1990，Vol．12 NO．4：24-28
[58] 何发详，刘浩吾．有限元法前处理技术的发展，水利水电科技进步，2000，12(6)：15-18
[59] 谭也平，有限元前处理系统的图形集成实现，计算机辅助工程，2000，Mar．NO．1：53-58
[60] 罗特军，罗季军，汪榴．有限元网格优化方法，四川联合大学学报(工程科学版)，1999，Vol．3 NO．3：65-72
[61] 袁明道，肖明．有限元网格半自动剖分及显示，武汉水利电力大学学报，1999，Vol．32 NO．3：72-75
[62] 张广建，黄由玲．用CAGD进行有限元前处理，工程力学，1995，Vol．12 NO．4：124-130
[63] 安晓卫，邵伟平，用于二维和三维问题的有限元前处理程序机械设计与制造，1995，(1)，P21-23
[64] 梁雪春，崔洪斌，吴义忠，曹康．AutoLISP实用教程．北京：人民邮电出版社．1999
[65] 李水乡，刘剑飞，袁明武．Auto CAD平台上的有限元前处理软件Auto MESH，计算力学学报，1999，Vol．16 NO．3：375-378
[66] K.Holc, Finite element mesh generation methods, a review and classification,Computer-aided Design, 1998, 20(1),27-32
[67] 水工建筑物抗震设计规范．DL5073-2000．北京：中国电力出版社，2001．
[68] 沈聚敏，周锡元，高小旺等．抗震工程学．北京：中国建筑工业出版社，2000．
[69] 后荣昌，韩渭宾主编．四川活动断裂与地震．北京：地震出版社，1993
[70] 李善邦．中国地震．北京：地震出版社，1981
[71] 陈厚群等．水电工程抗震设防概率水准和地震作用的概率模型．自然灾害学报．1993(2)
[72] 李桂青．结构动力可靠性理论及其应用[M]．北京：地震出版社，1993．
[73] 严松宏．结构地震响应分析的广义脉冲函数法．兰州铁道学院学报(自然科学版)．2002，Vol．21，No．6 Dec．
[74] 王进廷，杜修力，张楚汉．瑞利阻尼介质有限元离散模型．动力分析的数值稳定性．2002，Dec．Vol．22 No．6
[75] 李小军，廖振鹏，杜修力．有阻尼体系动力问题的一种显式解法[J]．地震工程与工程振动，1992，12(4)：74-79
[76] 王进廷，杜修力．有阻尼体系动力分析的一种显式差分法[J]．工程力学．2002，19(3)：109～112．
[77] 杜修力，王进廷．阻尼弹性结构动力计算的显式差分法[J]．工程力学，2002，17(5)：37～43．
[78] 寥振鹏，刘晶波．离散网格中的弹性波动(Ⅰ)[J]．地震工程与工程振动，1986，6(2)：1～15．
[79] Zhang Chuhan、Q. A. Pekau and Jin Fene. Application of FE—BE—IBE coupling to dynamic inteaction between alluvial soil and rock canyons[J]. Earthquake Engineering and Structural Dynamics, 1992, 21 (5): 367—385.
[80] 刘纪陆．求解材料非线性系统动力响应的新方法．建筑结构学报．2002，Vol．23，No．4，Aug．
[81] 王岚，陈中林．反应谱的数值计算及其应用中的误差问题．地壳形变与地震，1993，13(3)：85—89
[82] Ju JW. On energy based coupled elastoplastic damage theories: constitutive modeling and computational aspects[J] . Int. J. Solids Structure, 1989, 25(7): 803-833
[83] 苏克忠等．大坝强震安全监测．北京：中国水利水电出版社，1996
[84] 潘吕实，张弥，吴鸿庆．隧道力学数值方法[M]．北京：中国铁道出版社，1995