小湾水电站蚀变岩蠕变力学特性及本构模型研究
|
摘要
小湾水电站大坝为混凝土双曲拱坝,其坝肩抗力体内发育不同性状的蚀变岩带,蚀变岩带的性状及其工程特性,特别是其蠕变特性对坝肩岩体变形及稳定性起控制作用,甚至影响大坝的安全。本论文就是结合这一工程问题,重点对蚀变岩的蠕变力学特性进行了深入分析、研究。
小湾蚀变岩属于热液蚀变成因,断裂构造是热液上升的通道,蚀变岩的规模仅仅局限于断裂破碎带及其两侧,岩体结构较破碎,故在工程中被视为软岩。
在研究中通过现场调研、勘察,分析了蚀变岩带在坝肩的空间分布规律,并且充分考虑了试验的代表性,在坝肩抗力体的不同部位,采集了大量蚀变岩,进行了室内试验,获得了蚀变岩的物理性质指标、长期强度和粘滞系数等蠕变力学特性参数。为处理坝肩岩体和坝体的设计提供了可靠的试验参数。在深入分析蚀变岩蠕变特性的基础上,归纳出了适合小湾电站坝肩蚀变岩的蠕变本构模型。利用岩土分析软件FLAC~(3D)对蚀变岩的单轴及三轴蠕变试验进行了数值模拟研究,并将单轴压缩蠕变试验模拟结果与室内试验结果进行了比较分析,得到满意结果。通过改变受力条件,对常规三轴压缩蠕变试验进行了数值模拟研究,分析了三向应力作用下,蚀变岩的蠕变特性。为进一步研究蚀变岩的蠕变特性进行了有益的探讨。
The Xiaowan Hydropower station has the concrete's dome dam, there is various Rotten-rock belt in the its dam abutment. the Rotten-rock belt's character and it's project characteristic, specially its creep property plays the control action to the distortion of dam abutment, even affects the security of the dam. This paper emphasizing analyze the Rotten-rock's creep mechanical characteristic unifing this project question.
The Rotten-rock of Xiaowan belongs to the hydrothermal alteration, and the rupture structure becomes the ascending channels of the heat liquid. the distribution of rotten-rock just is located in the rupture zone and its adjacent zone. the structure of its periphery rock mass is relatively fracture. So it be regarded as the soft rock in the project.
Analyzed the Rotten-rock belt's spatial distribution rule which is in the dam abutment, through the investigation and study in the locale. and Selected the different Rotten-rock to adopt the test specimen in the dam abutment, and tested in the indoor experiment. Obtained the data of physical property and the creep mechanics characteristic, such as long-term strength and viscosity coefficient of the Rotten-rock, and provided the reliable mechanics parameter for the dam abutment's design. Summarized creep constitutive-model which is suited the Xiaowan Rotten-rock, through analysis creep property of Rotten-rock. Utilized FLAC~(3D) which is a numerical simulation software to analyze the single and three axes compressive creep tests, and carried on the comparative analysis the indoor creep test with the test result of numerical simulation. Finally, operated the numerical simulation with the three axes compressive creep tests,and analyzed the creep property. This article carried on the beneficial discussion for further studying the creep property of the Rotten-rock.
小湾蚀变岩属于热液蚀变成因,断裂构造是热液上升的通道,蚀变岩的规模仅仅局限于断裂破碎带及其两侧,岩体结构较破碎,故在工程中被视为软岩。
在研究中通过现场调研、勘察,分析了蚀变岩带在坝肩的空间分布规律,并且充分考虑了试验的代表性,在坝肩抗力体的不同部位,采集了大量蚀变岩,进行了室内试验,获得了蚀变岩的物理性质指标、长期强度和粘滞系数等蠕变力学特性参数。为处理坝肩岩体和坝体的设计提供了可靠的试验参数。在深入分析蚀变岩蠕变特性的基础上,归纳出了适合小湾电站坝肩蚀变岩的蠕变本构模型。利用岩土分析软件FLAC~(3D)对蚀变岩的单轴及三轴蠕变试验进行了数值模拟研究,并将单轴压缩蠕变试验模拟结果与室内试验结果进行了比较分析,得到满意结果。通过改变受力条件,对常规三轴压缩蠕变试验进行了数值模拟研究,分析了三向应力作用下,蚀变岩的蠕变特性。为进一步研究蚀变岩的蠕变特性进行了有益的探讨。
The Xiaowan Hydropower station has the concrete's dome dam, there is various Rotten-rock belt in the its dam abutment. the Rotten-rock belt's character and it's project characteristic, specially its creep property plays the control action to the distortion of dam abutment, even affects the security of the dam. This paper emphasizing analyze the Rotten-rock's creep mechanical characteristic unifing this project question.
The Rotten-rock of Xiaowan belongs to the hydrothermal alteration, and the rupture structure becomes the ascending channels of the heat liquid. the distribution of rotten-rock just is located in the rupture zone and its adjacent zone. the structure of its periphery rock mass is relatively fracture. So it be regarded as the soft rock in the project.
Analyzed the Rotten-rock belt's spatial distribution rule which is in the dam abutment, through the investigation and study in the locale. and Selected the different Rotten-rock to adopt the test specimen in the dam abutment, and tested in the indoor experiment. Obtained the data of physical property and the creep mechanics characteristic, such as long-term strength and viscosity coefficient of the Rotten-rock, and provided the reliable mechanics parameter for the dam abutment's design. Summarized creep constitutive-model which is suited the Xiaowan Rotten-rock, through analysis creep property of Rotten-rock. Utilized FLAC~(3D) which is a numerical simulation software to analyze the single and three axes compressive creep tests, and carried on the comparative analysis the indoor creep test with the test result of numerical simulation. Finally, operated the numerical simulation with the three axes compressive creep tests,and analyzed the creep property. This article carried on the beneficial discussion for further studying the creep property of the Rotten-rock.
引文
[1]刘雄.岩石流变学概论[M].北京:地质出版社,1994.
[2]周德培.流变力学原理及其在岩土工程中的应用[M].成都:西南交通大学出版社,1995.
[3]中华人民共和国水利部.中华人民共和国国家标准《工程岩体分级标准》GB50218-94.北京:中国计划出版社,1995.
[4]中华人民共和国建设部.中华人民共和国国家标准《岩土工程勘探规范》GBS0021-94.北京:中国建筑工业出版社,1995.
[5]周维垣.高等岩石力学[M].北京:水利水电出版社,1993.
[6]蔡美峰,何满潮,刘东燕.岩石力学与工程[M].北京:科学出版社,2002.
[7]方安平,叶卫平等.Origin7.5科技绘图及数据分析[M].北京:机械工业出版社,2006.
[8]孙钧.岩土材料流变及其工程应用[M].北京:中国建筑工业出版社,1999.
[9]王思敬.西部大开发中的岩石力学与工程问题[A].见:第八次全国岩石力学与工程学术会议论文集[C].北京:科学出版社,2004,1-9.
[10]陈建叶.锦屏一级高拱坝坝肩稳定三维地质力学模型破坏试验研究[D].四川大学硕士学位论文,2004.
[11]中国岩石力学与工程学会主编.西部大开发中的岩石力学与工程问题[A]:见:第八次全国岩石力学与工程学术会议论文集[C].北京:科学出版社,2004,1-9.
[12]Griggs D.Creep of rocks[J].Geol,1939,47.
[13]ZischinskyU.On the deformation of slopes[A].In:Proc.1~(st)Int.Cong.RockMech.[C].Lissobon:[s.n.],1966.
[14]Langer M.Rheological behavior of rock masses[A].In:Proc.of4th Cong.ofLSRM[C].[s.1]:[s.n.],1979.
[15]郑哲敏.21世纪初的力学发展趋势[J].力学进展,1995,25(4):433-441.
[16]范庆忠,王素华,高延法.岩石流变试验与本构模型研究进展[J].山东农业大学学报(自然科学版),2006,37(1):136-140.
[17]胡洲,刁心宏,王继飞.浅谈软岩流变性的研究[J].江西科学,2006,24(3):309-313.
[18]李永盛.单轴压缩条件下四种岩石的蠕变和松弛试验研究[J].岩石力学与工程学报,1995,16(1):39-47.
[19]Li Yongsheng Xia Caichu Time- dependent tests on intact rocks in uniaxial compression[J].Int J Rock Mech Mine Sci and Geomech Abstr 2000,37:467-475.
[20]张向东,郑雨天,肖裕性.第三系软弱岩体蠕变理论[J].东北大学学报,1997,18(1):31-35.
[21]金丰年.岩石拉压特性的相似性[J].岩土工程学报,1998,20(3):31-33.
[22]吴刚.红砂岩卸荷破坏特性的试验研究[A].岩土力学与工程[C].第二届全国岩土力学与工程育年工作者学术讨论会论文集,大连:大连理工大学出版社,1995,228-236.
[23]Sun Jun,Hu Y Y.Time-dependent effects on the tensile strength of saturated granite at Three Gorges Project in China[J].Int J Rock Mech Mine Sci 1997,34:381-384.
[24]陈有亮,刘涛.岩石流变断裂扩展的力学分析[J].上海大学学报,2000,6(6):491-496.
[25]张忠亭,罗居剑.分级加载下岩石蠕变特性研究[J].岩石力学与工程学报,2004,23(2):218-222.
[26]曹树刚,边金,李鹏.岩石蠕变本构关系及改进的西原正夫模型[J].岩石力学与工程学报.2002,21(5):632-634.
[27]邓荣贵,周德培,张倬元,付小敏.一种新的岩石流变模型[J].岩石力学与工程学报,2001,20(6):780-784.
[28]陈沅江,潘长良,曹平等.软岩流变的一种新力学模型[J].岩土力学,2003,24(2):209-214.
[29]韦立德,徐卫亚,朱珍德等.岩石粘弹塑性模型的研究[J].岩土力学,2002,23(5):583-586.
[30]余启华.岩石的流变破坏过程及有限元分析[J].水利学报,1985,(1):55-61.
[31]范庆忠,高延法.软岩蠕变特性及非线性模型研究[J].岩石力学与工程学报.2007,26(2):391-396.
[32]余寿文,冯西桥.损伤力学[M].北京:清华大学出版社,1997,237-242.
[33]封志军,周德培,周应华.origin软件在岩土流变试验数据整理中的应用[J].路基工程,2005,(6):4-6.
[34]王来贵,何峰,刘向峰等.岩石试件非线性蠕变模型及其稳定性分析[J].岩石力学与工程学报,2004,23(10):1640-1642.
[35]范庆忠,高延法.分级加载条件下岩石流变特性的试验研究[J].岩土工程学报,2005,27(11):1273-1276.
[36]肖庆华,徐卫亚,巫德斌.岩石流变力学试验数据挖掘研究[J].工程地质学报,2004,12(02):215-220.
[37]李化敏,李振华,苏承东.大理岩蠕变特性试验研究[J].岩石力学与工程学报,2004,23(22):3745-3749.
[38]崔希海,付志亮.岩石流变特性及长期强度的试验研究[J].岩石力学与工程学报,2006,25(5):1021-1024.
[39]杨彩红,李剑光.非均匀软岩蠕变机理分析[J].采矿与安全工程学报.2006,23(4): 476-479.
[40]陈宗基.地下巷道长期稳定性的力学研究[J].岩石力学与工程学报,2001,1(1).
[41]夏熙伦,徐平,丁秀丽.岩石流变特性及高边破稳定性流变分析[J].岩石力学与工程学报,1996,15(4):312-322.
[42]陈有亮,孙钧.岩石的蠕变断裂特性分析[J].同济大学学报,1996,24(5):504-508.
[43]陈有亮.岩石蠕变断裂特性的试验研究[J].力学学报,2003,33(4):480-484.
[44]刘保国,孙钧.岩体流变本构模型的辨识及应用[J].北方交通大学学报,1998,22(4):10-14.
[45]张向东,李永靖,张树光等.软岩蠕变理论及其工程应用[J].岩石力学与工程学报,2004,23(10):1635-1639.
[46]陈炳瑞,冯夏庭等.基于模式搜索的岩石流变模型参数识别[J].岩石力学与工程学报,2005,24(2):207-211.
[47]何满潮,晏玉书,王同良等.软岩的概念及其分类.世纪之交软岩工程技术现状与展望.煤炭工业出版社,1999:37-47.
[48]徐平,杨挺青.岩石流变试验与本构模型辨识[J].岩石力学与工程学报,2001,20(增):1739-1744.
[49]杨根兰.蚀变岩特性及其工程响应研究—以澜沧江小湾水电站为例[D]成都理工大学博士毕业论文,2007.
[50]葛修润,丰定祥.岩体工程中流变问题的有限元分析,岩土力学,1980年12月,第3期:1-14
[51]黄润秋,张倬元,王士天.高边坡稳定性的系统工程地质研究,成都:成都科技大学出版社,1991
[52]丁秀丽,徐平,夏熙伦.三峡船闸高边坡岩体开挖卸荷变形及流变分[J].1995,12(4):37-63.
[53]夏熙伦,徐平,丁秀丽.岩石流变特性及高边坡稳定性流变分析[J].岩石力学与工程学报.1996,15(4):312-322.
[54]章根德,剡公瑞.岩体高边坡流变学性状有限元分析[J].岩土工程学报,1999,21(2),166-170.
[55]曹文贵,速宝玉.岩体蠕变大变形有限元分析及其在金川矿的应用[J].矿冶工程,1999,19(3):22-24.
[56]徐平,甘军,丁秀丽.三峡工程船闸高边坡岩体长期变形及稳定性有限元分析[J].长江科学院院报,1999,16(2):31-34.
[57]徐平,杨挺青,徐春敏等.三峡船闸高边坡岩体时效特性及长期稳定性分析[J].岩石力学与工程学报,2002,21(2):163-168.
[58]刘建华,朱维申,李术才小浪底水利枢纽地下厂房岩体流变与稳定性FLAC3D数值分析 [J].岩石力学与工程学报,2005,24(14):2 484-2 489.
[59]丁秀丽,付敬,刘建等.软硬互层边坡岩体的蠕变特性研究及稳定性分析[J],岩石力学与工程学报,2005,24(19)P:3 410-3 418.
[60]尹雪英,杨春和,陈剑文.金坛盐矿老腔储气库长期稳定性分析数值模拟[J].,岩土力学,2006,27(6),869-874.
[61]陈卫忠,伍国军,戴永浩等.废弃盐穴地下储气库稳定性研究[J].岩石力学与工程学报.2006,25(4),848-854.
[62]刘建华.岩体力学行为拉格朗日分析方法研究与工程应用[D].山东大学博士论文,2006.国外
[63]胡受奚,叶瑛,方长泉.交代蚀变岩岩石学及其找矿意义[M].北京:地质出版社,2004.
[64]李强,杨庆.尺寸效应对掩体破坏的影响[J].隧道、地下工程及岩石破碎理论与应用,2007.
[65]刘宝琛、张家生等.岩石抗压强度的尺寸效应[J].岩石力学与工程学报,1998,17(6):611-614.
[66]张学言等.岩土塑性力学基础[M].天津:天津大学出版社,2004.9.
[67]李同林,殷绥域.弹塑性力学[M].武汉:中国地质大学出版社,2006.9.
[68]唐辉明,晏鄂川,胡新丽.工程地质数值模拟的理论与方法[M].武汉:中国地质大学出版社,2001.10.
[69]Tan Tjong-kie(陈宗基).Kang Wen-fa.Lock in stress,dreep and dilatancy of rocks and constitutive equation.Rock Mechanics.1980.13:5-22.
[70]褚卫江,徐卫亚等.基于FLAC~(3D)岩石黏弹塑性流变模型的二次开发研究[J].岩土力学,2006,27(11):2005-2010.
[71]刘波,韩彦辉(美国).FLAC原理、实例与应用指南[M].北京:人民交通出版社,2005,9.
[2]周德培.流变力学原理及其在岩土工程中的应用[M].成都:西南交通大学出版社,1995.
[3]中华人民共和国水利部.中华人民共和国国家标准《工程岩体分级标准》GB50218-94.北京:中国计划出版社,1995.
[4]中华人民共和国建设部.中华人民共和国国家标准《岩土工程勘探规范》GBS0021-94.北京:中国建筑工业出版社,1995.
[5]周维垣.高等岩石力学[M].北京:水利水电出版社,1993.
[6]蔡美峰,何满潮,刘东燕.岩石力学与工程[M].北京:科学出版社,2002.
[7]方安平,叶卫平等.Origin7.5科技绘图及数据分析[M].北京:机械工业出版社,2006.
[8]孙钧.岩土材料流变及其工程应用[M].北京:中国建筑工业出版社,1999.
[9]王思敬.西部大开发中的岩石力学与工程问题[A].见:第八次全国岩石力学与工程学术会议论文集[C].北京:科学出版社,2004,1-9.
[10]陈建叶.锦屏一级高拱坝坝肩稳定三维地质力学模型破坏试验研究[D].四川大学硕士学位论文,2004.
[11]中国岩石力学与工程学会主编.西部大开发中的岩石力学与工程问题[A]:见:第八次全国岩石力学与工程学术会议论文集[C].北京:科学出版社,2004,1-9.
[12]Griggs D.Creep of rocks[J].Geol,1939,47.
[13]ZischinskyU.On the deformation of slopes[A].In:Proc.1~(st)Int.Cong.RockMech.[C].Lissobon:[s.n.],1966.
[14]Langer M.Rheological behavior of rock masses[A].In:Proc.of4th Cong.ofLSRM[C].[s.1]:[s.n.],1979.
[15]郑哲敏.21世纪初的力学发展趋势[J].力学进展,1995,25(4):433-441.
[16]范庆忠,王素华,高延法.岩石流变试验与本构模型研究进展[J].山东农业大学学报(自然科学版),2006,37(1):136-140.
[17]胡洲,刁心宏,王继飞.浅谈软岩流变性的研究[J].江西科学,2006,24(3):309-313.
[18]李永盛.单轴压缩条件下四种岩石的蠕变和松弛试验研究[J].岩石力学与工程学报,1995,16(1):39-47.
[19]Li Yongsheng Xia Caichu Time- dependent tests on intact rocks in uniaxial compression[J].Int J Rock Mech Mine Sci and Geomech Abstr 2000,37:467-475.
[20]张向东,郑雨天,肖裕性.第三系软弱岩体蠕变理论[J].东北大学学报,1997,18(1):31-35.
[21]金丰年.岩石拉压特性的相似性[J].岩土工程学报,1998,20(3):31-33.
[22]吴刚.红砂岩卸荷破坏特性的试验研究[A].岩土力学与工程[C].第二届全国岩土力学与工程育年工作者学术讨论会论文集,大连:大连理工大学出版社,1995,228-236.
[23]Sun Jun,Hu Y Y.Time-dependent effects on the tensile strength of saturated granite at Three Gorges Project in China[J].Int J Rock Mech Mine Sci 1997,34:381-384.
[24]陈有亮,刘涛.岩石流变断裂扩展的力学分析[J].上海大学学报,2000,6(6):491-496.
[25]张忠亭,罗居剑.分级加载下岩石蠕变特性研究[J].岩石力学与工程学报,2004,23(2):218-222.
[26]曹树刚,边金,李鹏.岩石蠕变本构关系及改进的西原正夫模型[J].岩石力学与工程学报.2002,21(5):632-634.
[27]邓荣贵,周德培,张倬元,付小敏.一种新的岩石流变模型[J].岩石力学与工程学报,2001,20(6):780-784.
[28]陈沅江,潘长良,曹平等.软岩流变的一种新力学模型[J].岩土力学,2003,24(2):209-214.
[29]韦立德,徐卫亚,朱珍德等.岩石粘弹塑性模型的研究[J].岩土力学,2002,23(5):583-586.
[30]余启华.岩石的流变破坏过程及有限元分析[J].水利学报,1985,(1):55-61.
[31]范庆忠,高延法.软岩蠕变特性及非线性模型研究[J].岩石力学与工程学报.2007,26(2):391-396.
[32]余寿文,冯西桥.损伤力学[M].北京:清华大学出版社,1997,237-242.
[33]封志军,周德培,周应华.origin软件在岩土流变试验数据整理中的应用[J].路基工程,2005,(6):4-6.
[34]王来贵,何峰,刘向峰等.岩石试件非线性蠕变模型及其稳定性分析[J].岩石力学与工程学报,2004,23(10):1640-1642.
[35]范庆忠,高延法.分级加载条件下岩石流变特性的试验研究[J].岩土工程学报,2005,27(11):1273-1276.
[36]肖庆华,徐卫亚,巫德斌.岩石流变力学试验数据挖掘研究[J].工程地质学报,2004,12(02):215-220.
[37]李化敏,李振华,苏承东.大理岩蠕变特性试验研究[J].岩石力学与工程学报,2004,23(22):3745-3749.
[38]崔希海,付志亮.岩石流变特性及长期强度的试验研究[J].岩石力学与工程学报,2006,25(5):1021-1024.
[39]杨彩红,李剑光.非均匀软岩蠕变机理分析[J].采矿与安全工程学报.2006,23(4): 476-479.
[40]陈宗基.地下巷道长期稳定性的力学研究[J].岩石力学与工程学报,2001,1(1).
[41]夏熙伦,徐平,丁秀丽.岩石流变特性及高边破稳定性流变分析[J].岩石力学与工程学报,1996,15(4):312-322.
[42]陈有亮,孙钧.岩石的蠕变断裂特性分析[J].同济大学学报,1996,24(5):504-508.
[43]陈有亮.岩石蠕变断裂特性的试验研究[J].力学学报,2003,33(4):480-484.
[44]刘保国,孙钧.岩体流变本构模型的辨识及应用[J].北方交通大学学报,1998,22(4):10-14.
[45]张向东,李永靖,张树光等.软岩蠕变理论及其工程应用[J].岩石力学与工程学报,2004,23(10):1635-1639.
[46]陈炳瑞,冯夏庭等.基于模式搜索的岩石流变模型参数识别[J].岩石力学与工程学报,2005,24(2):207-211.
[47]何满潮,晏玉书,王同良等.软岩的概念及其分类.世纪之交软岩工程技术现状与展望.煤炭工业出版社,1999:37-47.
[48]徐平,杨挺青.岩石流变试验与本构模型辨识[J].岩石力学与工程学报,2001,20(增):1739-1744.
[49]杨根兰.蚀变岩特性及其工程响应研究—以澜沧江小湾水电站为例[D]成都理工大学博士毕业论文,2007.
[50]葛修润,丰定祥.岩体工程中流变问题的有限元分析,岩土力学,1980年12月,第3期:1-14
[51]黄润秋,张倬元,王士天.高边坡稳定性的系统工程地质研究,成都:成都科技大学出版社,1991
[52]丁秀丽,徐平,夏熙伦.三峡船闸高边坡岩体开挖卸荷变形及流变分[J].1995,12(4):37-63.
[53]夏熙伦,徐平,丁秀丽.岩石流变特性及高边坡稳定性流变分析[J].岩石力学与工程学报.1996,15(4):312-322.
[54]章根德,剡公瑞.岩体高边坡流变学性状有限元分析[J].岩土工程学报,1999,21(2),166-170.
[55]曹文贵,速宝玉.岩体蠕变大变形有限元分析及其在金川矿的应用[J].矿冶工程,1999,19(3):22-24.
[56]徐平,甘军,丁秀丽.三峡工程船闸高边坡岩体长期变形及稳定性有限元分析[J].长江科学院院报,1999,16(2):31-34.
[57]徐平,杨挺青,徐春敏等.三峡船闸高边坡岩体时效特性及长期稳定性分析[J].岩石力学与工程学报,2002,21(2):163-168.
[58]刘建华,朱维申,李术才小浪底水利枢纽地下厂房岩体流变与稳定性FLAC3D数值分析 [J].岩石力学与工程学报,2005,24(14):2 484-2 489.
[59]丁秀丽,付敬,刘建等.软硬互层边坡岩体的蠕变特性研究及稳定性分析[J],岩石力学与工程学报,2005,24(19)P:3 410-3 418.
[60]尹雪英,杨春和,陈剑文.金坛盐矿老腔储气库长期稳定性分析数值模拟[J].,岩土力学,2006,27(6),869-874.
[61]陈卫忠,伍国军,戴永浩等.废弃盐穴地下储气库稳定性研究[J].岩石力学与工程学报.2006,25(4),848-854.
[62]刘建华.岩体力学行为拉格朗日分析方法研究与工程应用[D].山东大学博士论文,2006.国外
[63]胡受奚,叶瑛,方长泉.交代蚀变岩岩石学及其找矿意义[M].北京:地质出版社,2004.
[64]李强,杨庆.尺寸效应对掩体破坏的影响[J].隧道、地下工程及岩石破碎理论与应用,2007.
[65]刘宝琛、张家生等.岩石抗压强度的尺寸效应[J].岩石力学与工程学报,1998,17(6):611-614.
[66]张学言等.岩土塑性力学基础[M].天津:天津大学出版社,2004.9.
[67]李同林,殷绥域.弹塑性力学[M].武汉:中国地质大学出版社,2006.9.
[68]唐辉明,晏鄂川,胡新丽.工程地质数值模拟的理论与方法[M].武汉:中国地质大学出版社,2001.10.
[69]Tan Tjong-kie(陈宗基).Kang Wen-fa.Lock in stress,dreep and dilatancy of rocks and constitutive equation.Rock Mechanics.1980.13:5-22.
[70]褚卫江,徐卫亚等.基于FLAC~(3D)岩石黏弹塑性流变模型的二次开发研究[J].岩土力学,2006,27(11):2005-2010.
[71]刘波,韩彦辉(美国).FLAC原理、实例与应用指南[M].北京:人民交通出版社,2005,9.