地震耦合作用下顺层岩质边坡动力响应研究
摘要
我国是一个多山的国家,山地、高原和丘陵占国土面积的69%,这就从客观上决定了我国有大量的自然边坡。同时我国地处世界上两条最活跃的地震带之间,强震分布广,其中西部地区地震活动性强、频度高,近年来的地震活动表明我国地震活动正趋于活跃。地震作用下边坡的动力响应研究是边坡动力稳定分析的基础,近半个世纪以来,人们对边坡动力反应问题的研究,倾注了极大的热情,并取得了一些成绩,但研究对象多局限于均质概化模型边坡,研究程度较低。
本文首先归纳总结了边坡动力响应方面的研究现状,接着针对顺层岩质边坡,使用FLAC3D软件模拟了二维顺层岩质边坡模型在双向耦合地震情况下的动力响应特征、三维顺层岩质边坡模型在双向耦合地震情况下的动力响应特征和三维顺层岩质边坡模型在三向耦合地震情况下的动力响应特征,分析了坡面的位移时程曲线,最大位移分布,PGA放大系数分布,取得了如下一些主要认识：
1.二维顺层岩质边坡在双向耦合地震情况下的坡面位移的变化主要由X方向加速度控制,坡面上的位移时程变化和X方向加速度的峰值变化成对应关系；耦合地震作用下X方向最大位移增大了39.3%,说明耦合地震作用模拟比单纯的水平地震作用模拟更加接近实际情况,对岩土体的破坏效应更大；双向耦合地震情况下,随着av/aH的增大,坡面各点横向最大位移基本呈增大趋势,证明地震波产生的水平向与竖向拉裂耦合作用是触发边坡产生初期崩滑破坏的主控因素,竖向地震作用起了重要的破坏作用,而av/aH越大,破坏越剧烈,这也印证了许多地震工作者发现极震区竖向地震力起破坏作用的事实。
2.三维边坡模型的动力响应特征有其作为三维体而独有的特点,特别是坡面长度(Z向)上的动力响应特征分布是二维边坡模型所不具有的,所以对于边坡的数值模型计算应该根据边坡的实际长度来建立三维模型,这样可以使模拟结果更加接近实际情况。
3.三向耦合地震是最接近现实地震情况的模拟方式,在三向耦合地震作用下三维顺层岩质边坡模型的模拟结果最接近实际情况,不仅可以作为定性研究的结论,在一定程度上也可以作为定量分析的依据；坡面方向(Z向)动力响应特征值的不同分布也符合了现实中边坡的滑移存在先后顺序的情况。
China is a mountainous country and mountain, plateau and hills accounts for69%of land area, so it objectively determines that there are plenty of natural slopes in our country. Meanwhile, our country is located between two of the most active earthquake zones, which makes earthquake distribution is wide and seismic activity is strong and frequency is high in the Middle Western region. And the seismic activity in recent years show that seismic activity in China is tend to be more active. Study on slope dynamic response under earthquake force is the basis of stability analysis of slope. Nearly half a century, reaearchers have poured nenormous enthusiasm into the dynamic response of slope and have made some achivment. But most of the research object modek is the relatively simple homogeneous slope, and the research degree of this problem at home and abroad is relatively low.
This paper summarized the research status quo of the slope dynamic response, then the bedding rock slope is selected for the study direction. By using FLAC3D, author simulates the dynamic response characteristics of2D model of bedding rock slope under coupling vibration in two cases, the dynamic response characteristics of3D model of bedding rock slope under coupling vibration in two cases and the dynamic response characteristics of3D model of bedding rock slope under the condition of three-way coupling vibration. By analyzing of the slope surface displacement, time history curve of maximum displacement, and change of PGA amplification coefficient, author has made some major conclutions as follows:
The change of slope surface displacement in the2D model of bedding rock slope under coupling vibration in two cases are mainly controlled by the level of the acceleration in X direction, and the slope displacement time history on the surface has a corresponding relationship with the peak value of acceleration of X direction. X direction maximum displacement under the coupling vibration increases by39.3%, which shows that coupling vibration effect is closer to actual situation and to the destruction of the rock mass effect is bigger than simple horizontal vibration simulation. Under coupling vibration in two case, all point slope lateral maximum displacement showed a trend of increase with the increase of the aV/aH, which proves slope body collapse is trigger by the effect of the coupling of horizontal and vertical tensile splitting seismic wave and the destruction of the vertical seismic action plays an important role. Meanwhile, the greater the aV/aH, the more severe damage, which confirms many workers'founds that the vertical seismic force play damage effect in magistoseismic area.
The slope dynamic response characteristics of3D model has its unique characteristics as3D body, especially the distribution of the dynamic response characteristics along the slope length (Z-direction) dosen't exist in2D slope model. Therefore, the numerical model calculation for slope should build3D model according to the actual length of the slope. As a result, we can make the simulation results closer to actual situation.
Three-way coupling vibration is the simulation methods which is most close to the real situation of the earthquake and simulation results of3D bedding rock slope model under three-way coupling vibration is quite close to the actual situation, which not only can be used as a qualitative conclusion, to a certain extent, can also be used as the basis for quantitative. Different distribution of dynamic response characteristics along slope direction (Z-direction) also conforms to the reality of the slope slip points successively.
本文首先归纳总结了边坡动力响应方面的研究现状,接着针对顺层岩质边坡,使用FLAC3D软件模拟了二维顺层岩质边坡模型在双向耦合地震情况下的动力响应特征、三维顺层岩质边坡模型在双向耦合地震情况下的动力响应特征和三维顺层岩质边坡模型在三向耦合地震情况下的动力响应特征,分析了坡面的位移时程曲线,最大位移分布,PGA放大系数分布,取得了如下一些主要认识：
1.二维顺层岩质边坡在双向耦合地震情况下的坡面位移的变化主要由X方向加速度控制,坡面上的位移时程变化和X方向加速度的峰值变化成对应关系；耦合地震作用下X方向最大位移增大了39.3%,说明耦合地震作用模拟比单纯的水平地震作用模拟更加接近实际情况,对岩土体的破坏效应更大；双向耦合地震情况下,随着av/aH的增大,坡面各点横向最大位移基本呈增大趋势,证明地震波产生的水平向与竖向拉裂耦合作用是触发边坡产生初期崩滑破坏的主控因素,竖向地震作用起了重要的破坏作用,而av/aH越大,破坏越剧烈,这也印证了许多地震工作者发现极震区竖向地震力起破坏作用的事实。
2.三维边坡模型的动力响应特征有其作为三维体而独有的特点,特别是坡面长度(Z向)上的动力响应特征分布是二维边坡模型所不具有的,所以对于边坡的数值模型计算应该根据边坡的实际长度来建立三维模型,这样可以使模拟结果更加接近实际情况。
3.三向耦合地震是最接近现实地震情况的模拟方式,在三向耦合地震作用下三维顺层岩质边坡模型的模拟结果最接近实际情况,不仅可以作为定性研究的结论,在一定程度上也可以作为定量分析的依据；坡面方向(Z向)动力响应特征值的不同分布也符合了现实中边坡的滑移存在先后顺序的情况。
China is a mountainous country and mountain, plateau and hills accounts for69%of land area, so it objectively determines that there are plenty of natural slopes in our country. Meanwhile, our country is located between two of the most active earthquake zones, which makes earthquake distribution is wide and seismic activity is strong and frequency is high in the Middle Western region. And the seismic activity in recent years show that seismic activity in China is tend to be more active. Study on slope dynamic response under earthquake force is the basis of stability analysis of slope. Nearly half a century, reaearchers have poured nenormous enthusiasm into the dynamic response of slope and have made some achivment. But most of the research object modek is the relatively simple homogeneous slope, and the research degree of this problem at home and abroad is relatively low.
This paper summarized the research status quo of the slope dynamic response, then the bedding rock slope is selected for the study direction. By using FLAC3D, author simulates the dynamic response characteristics of2D model of bedding rock slope under coupling vibration in two cases, the dynamic response characteristics of3D model of bedding rock slope under coupling vibration in two cases and the dynamic response characteristics of3D model of bedding rock slope under the condition of three-way coupling vibration. By analyzing of the slope surface displacement, time history curve of maximum displacement, and change of PGA amplification coefficient, author has made some major conclutions as follows:
The change of slope surface displacement in the2D model of bedding rock slope under coupling vibration in two cases are mainly controlled by the level of the acceleration in X direction, and the slope displacement time history on the surface has a corresponding relationship with the peak value of acceleration of X direction. X direction maximum displacement under the coupling vibration increases by39.3%, which shows that coupling vibration effect is closer to actual situation and to the destruction of the rock mass effect is bigger than simple horizontal vibration simulation. Under coupling vibration in two case, all point slope lateral maximum displacement showed a trend of increase with the increase of the aV/aH, which proves slope body collapse is trigger by the effect of the coupling of horizontal and vertical tensile splitting seismic wave and the destruction of the vertical seismic action plays an important role. Meanwhile, the greater the aV/aH, the more severe damage, which confirms many workers'founds that the vertical seismic force play damage effect in magistoseismic area.
The slope dynamic response characteristics of3D model has its unique characteristics as3D body, especially the distribution of the dynamic response characteristics along the slope length (Z-direction) dosen't exist in2D slope model. Therefore, the numerical model calculation for slope should build3D model according to the actual length of the slope. As a result, we can make the simulation results closer to actual situation.
Three-way coupling vibration is the simulation methods which is most close to the real situation of the earthquake and simulation results of3D bedding rock slope model under three-way coupling vibration is quite close to the actual situation, which not only can be used as a qualitative conclusion, to a certain extent, can also be used as the basis for quantitative. Different distribution of dynamic response characteristics along slope direction (Z-direction) also conforms to the reality of the slope slip points successively.
引文
[1]张咸恭,王思敬,张倬元,等.中国工程地质学[M].北京,科学出版社,2000.
[2]陈建君.复杂山区斜坡的地震动力响应分析[N].成都：成都理工大学.2009.
[3]杜晓丽.地震荷载作用下岩质边坡稳定性研究[D].西安科技大学.2008.
[4]祁生文.岩质边坡动力反应分析[M].北京：科学出版社.2007
[5]Keh-Jian Shou, Cheng-Fung Wang. Analysis of the Chiufengershan landslide triggered by the 1999 Chi-Chi earthquake in Taiwan[J]. Engineering Geology,2003,68:237-250.
[6]张勇,许力生,陈运泰.2010年青海玉树地震震源过程[J].中国科学：地球科学：2010.7(40)：819-821.
[7]陈建君.复杂山区斜坡的地震动力响应分析[D].导师：许强.成都理工大学,2009.
[8]祁生文.边坡动力响应分析及应用研究[D].博士论文,中国科学研究生院,2002.
[9]刘红帅,薄景山,刘德东.岩土边坡地震稳定性分析研究评述[J].地震工程与工程振动,2005,8(1)：165-167.
[10]Leshchinsky D.,San K.Ching.Pseudo-static seismic stability of slopes[J].Design charts. J.of Geotech.Engrg,1994,120(9):1514-1532.
[11]Baker R.,Garber M.Theoretical analysis of the stability of slopes[J].Geotechnique,1978, 28(4):395-411.
[12]Taylor D.W.Stability of earth slopes[J].J.Boston Society of Civ.Engrgs,1937,24(3):197-246.
[13]姚爱军,苏永华.复杂岩质边坡锚固工程地震敏感性分析[J].土木工程学报,2003,(11)：34-37.
[14]Seed H.B.Stability of earth and rock fill dams during earthquakes in Embankment-Dam Engrg[M].Casagrande,Vol.(Eds Hirschfeld and Poulos),JohnWiley,1973.
[15]Newmark N.M.Effects of earthquakes on dams and embankments[J].Geotechnique,1965, 15(2):139-160.
[16]RobertoR.Seismieallyindueedlandslidedisplaeements:apredietivemode[J].Eng Geol,2000, 55(3):337-351.
[17]刘汉东.考虑地震历时影响的岩质边坡楔体稳定性分析与计算[J].华北水利水电学院院报,1991,(4)：35-40.
[18]薛守义,王思敬,刘建中.块状岩体边坡地震滑动位移分析[J].工程地质学报,1997,5(2)：131-136.
[19]黄建梁,王威中,薛宏交.坡体地震稳定性的动态分析[J].地震工程与工程振动,1997,(04)：114-123.
[20]王思敬.岩石边坡动态稳定性的初步探讨[J].地质科学,1977,(04)：372-376.
[21]王思敬,张菊明.边坡岩体滑动稳定性的动力分析[J].地质科学.1982,(2)：162-170.
[22]王思敬,薛守义.岩体边坡楔形体动力学分析[J].地质科学.1992,(2)：177-182.
[23]王勇智,戚炜,门玉明,等.大柳树坝址岩体松动机制物理模拟研究[J].西安科技大学学报,2004,24(3)：305—308.
[24]王勇智,戚炜,门玉明,等.强震区岩体动力破坏机制研究[J].岩土力学,2005,26(11)：1841—1844.
[25]郝建斌,门玉明,彭建兵,等.层状岩体边坡动力稳定性模型试验研究[J].公路交通科技,2005,22(6)：72—75.
[26]门玉明,彭建兵,李寻昌,等.层状结构岩质边坡动力稳定性试验研究[J].世界地震工程,2004,20(4)：131—136.
[27]许强,陈建君,冯文凯.斜坡地震响应的物理模拟试验研究[J].四川大学学报(工学版),2009,41(3)：266—272.
[28]徐光兴,姚令侃,高召宁,等.斜坡动力特性与动力响应的大型振动台模型试验研究[J].岩石力学与工程学报,2008,27(3)：624—632.
[29]许强,刘汉香,邹威,等.斜坡加速度动力响应特性的大型振动台试验研究[J].岩石力学与工程学报,2010,29(12)：2420—2428.
[30]艾畅,冯春,李世海,等.地震作用下顺层岩质边坡动力响应的试验研究[J].岩石力学与工程学报.2010,29(9)：1825-1832.
[31]董金玉,杨国香,伍法权,等.地震作用下顺层岩质边坡动力响应和破坏模式大型振动台试验研究[J].岩土力学.2011,32(10)：2977-2983.
[32]董金玉,杨继红,杨国香,等.基于正交设计的模型试验相似材料的配比试验研究[J].煤炭学报.2012,37(1)：44-49.
[33]杨国香,伍法权,董金玉,等.地震作用下岩质边坡动力响应特性及变形破坏机制研究[J].岩石力学与工程学报.2012,31(4)：696-702.
[34]杨国香,叶海林,伍法权,等.反倾层状结构岩质边坡动力响应特性及破坏机制振动台模型试验研究[J].岩石力学与工程学报.2012,31(11)：2214-2221.
[35]Hayashi H, N ishikawa J.Taniguch i K. Seismic behavior of road embankment[C]//Proc Intl Conf Centrifuge 98.Rotterdam:A. A. Balkema,l 998:243-248.
[36]M adabhushi S P QHaigh S K.Subedi B R.Seismic behavior of steep slopes[C]//Proc Physical Modeling in Geotechnics. Rotterdam:A. A.Balkema,2002:489-494.
[37]王年香,章为民.混凝土面板堆石坝地震反应离心模型试验[J].水利水运工程学报,2003(1)：18-22.
[38]Ng C.W.W.,Li Xiang-Song,Van Laak P.A.,et al.Centrifuge modeling of loose fill embankment subjected to uni-axial and bi-axial earthquakes[J].Soil Dynamics& Earthquake Engineering,2004,24(4):305-318.
[39]邓丽军,于玉贞,陈正发.应用于边坡动力离心试验的模型箱设计与试验检测[C].第24届全国土工测试会议论文集.郑州：黄河水利出版社,2005：246-252.
[40]Ishizaki and Hatekeyama.Consideration on the dynamical behavior of earth dams[R].Bul. No.52.Disaster Prevention Research Inst,Kyoto Unvi.,1963.
[41]陈育民,徐鼎平FLAC/FLAC3D基础与工程实例[M].北京：中国水利水电出版社,2008.
[42]Clough R.W.and Chopra A.K.Earthquake stress analysis in earth dams[J].J.Engrg.Mech., ASCE 1966,92,EM2.
[43]何蕴龙,陆述远.岩石边坡地震作用近似计算方法[J].岩土工程学报,1998,20(2)：66-68
[44]钱胜国,陆秋蓉.长江三峡船闸高边坡地震稳定性分析[R].武汉：长江科学院科研报告,1991.
[45]彭定,陈玲玲.三峡船闸高边坡在地震作用下的稳定性分析[J].长江科学院院报,2001,18(6)：47-49.
[46]徐卫亚,谢守益,徐瑞春,等.水布垭马崖高陡边坡三维数值模拟研究[J].武汉水利电力大学 (宜昌)学报,1998,20(3)：11-15.
[47]蔡德所.地震输入波形对三峡工程左岸坝段的动力影响研究[J].岩石力学与工程学报,1998,17(2)：183-187.
[48]刘立平,杨实君,李英民.软夹层参数对边坡动力特性的影响[J].重庆大学学报,2007,30(5)：31-34.
[49]狄仁东,石玉成,付长华.高边坡在水平动荷载下的动力响应规律研究[J].2007,23(2)：131-138.
[50]张学东,言志信,张森.ANSYS在岩质边坡动力响应分析中的应用[J].西北地震学报,2010,32(2)：117-121.
[51]徐志英,沈珠江.地震液化有效应力二维动力分析方法[J].华东水利学院学报,1981,No.3：1-14.
[52]徐志英,沈珠江.土坝地震孔隙水压力产生、扩散和消散的有限单元法动力分析[J].华东水利学院学报,1981,No.4：1-16.
[53]徐志英,沈珠江.尾矿高堆坝地震反应的综合分析与液化计算[J].水利学报,1983,No.5：30-39.
[54]徐志英,周健.土坝地震孔隙水压力产生、扩散和消散的三维动力分析[J].地震工程与工程振动,1985,5(4)：57-72.
[55]徐志英,周健.奥罗维尔土坝三维排水有效应动力分析[J].水利学报,1991,No.6：19-27.
[56]徐志英.奥罗维尔土坝三维简化动力分析[J].岩土工程学报,1996,18(2)82-87.
[57]周健,徐志英.土尾矿坝的三维有效应力的动力反应分析[J].地震工程与工程振动,1984,4(3)：60-70.
[58]刘波,韩彦辉,等FLAC原理-实例与应用指南[M].北京：人民交通出版社,2005.
[59]陈志坚.层状岩质边坡工程安全监控建模理论及关键技术研究[D]河海大学,2001,11
[60]许强,黄润秋.5.12汶川大地震诱发大型崩滑灾害动力特征初探[J].工程地质学报,2008,16(6)：721-729.
[60]许强,陈建君,冯文凯,肖锐铧,等.斜坡地震响应的物理模拟试验研究[J].四川大学学报(工程科学版).2009,41(3).
[61]肖克强,李海波,刘亚群,等.地震荷载作用下顺层岩体边坡变形特征分析[J].岩土力学.2007,28(8)：1557-1564.
[62]言志信,张森,张学东,等.顺层岩质边坡地震动力响应及地震动参数影响研究[J].岩石力学与工程学报.2011,30(增2)：3522-3528.
[63]孙博,彭宁波,王逢睿.云冈石窟第19窟西耳窟地震动力响应[J].西南交通大学学报.2012,47(4)：573-579.
[64]于翔,钱七虎,赵跃堂,郭志昆.地铁工程结构破坏的竖向地震力影响分析[J].解放军工大学学报(自然科学版),2001,2(3)：75-77.
[65]尹紫红.地震作用下的高速公路滑坡稳定性分析[J].重庆工学院学报(自然科学版),2007,21(9)：92-95.
[66]黄润秋.汶川8.0级地震触发崩滑灾害机制及其地质力学模式[J].岩石力学与工程学报,2009,28(6)：1239-1249.
[67]许强,黄润秋.5.12汶川大地震诱发大型崩滑灾害动力特征初探[J].工程地质学报,2008,16(6)：721-729.
[68]殷跃平.汶川八级地震地质灾害研究[J].工程地质学报,2008,6(24)：433-444.
[69]非明伦,樊跃新,崔建文,包一峰,陈坤华,卢吉高.施甸5.9级地震震害与竖向地震力作用[J].地震研究,2003,26(3)：292-298.
[70]胡晓燕,周健,胡晓虎.地震引起的竖向力对地铁隧道的影响[J].工程抗震,2006,2：32-35.
[71]梁庆国,韩文峰,赵士耀.垂直向地震作用对节理岩体失稳破坏的影响[J],西北地震学报,2007,26(4)：307-313.
[72]侯超群,王晓谋,石恒俊.地震对边坡稳定性影响分析[J].路基工程,2008,6：187-188.
[73]何蕴龙,陆述远.岩石边坡地震作用近似计算方法[J].岩土工程学报,1998,20(2)：66-68.
[74]崔芳鹏,胡瑞林,殷跃平,许强,张明.地震纵横波时差耦合作用的斜坡崩滑效应研究[J].工程地质学报,2009,(04)：455-462.
[75]崔芳鹏,殷跃平,许强,胡瑞林,张明.地震纵横波时差耦合作用的斜坡崩滑响应研究[J].岩土工程学报,2010,(08)：1266-1273.
[76]崔芳鹏,许强,谭儒蛟,殷跃平.地震动力作用触发的斜坡崩滑效应模拟[J].同济大学学报(自然科学版),2011,(03)：445-450.
[77]崔芳鹏,殷跃平,许强,胡瑞林.地震动力作用触发的斜坡崩滑响应研究——以青川东河口滑坡为例[J].三峡大学学报(自然科学版),2011,(03)：23-28.
[78]崔芳鹏,胡瑞林,殷跃平,许强,张明.纵横波时差耦合作用的斜坡崩滑效应离散元分析——以北川唐家山滑坡为例[J].岩石力学与工程学报,2010,(02)：319-327.
[79]方建瑞,许志雄,庄晓莹.三维边坡稳定弹塑性有限元分析与评价[J].岩土力学,2008,(10)：2667-2672.
[80]张向东,张玉,陈慧刚.黄土边坡二维和三维变形破坏模型的对比分析[J].中国地质灾害与防治学报,2012,(03)：18-20+38.
[81]袁伟,黄茂松,刘怡林,王浩然.含软弱夹层土坡稳定性的三维弹塑性数值模拟[J].郑州大学学报(工学版),2012,(05)：1-4.
[82]胡聿贤.地震工程学[M].北京：地震出版社,2006.
[83]Biggs J. M. Introduction to Structural Dynamics[M].New York:McGraw-Hill,1964
[84]万洪.岩质高边坡动力响应分析及其应用研究[D].武汉：武汉理工大学,2009.
[85]叶海林,郑颖人,黄润秋,等.地震边坡破坏机制及其破裂面的分析探讨[J].岩石力学与工程学报,2009,28(8)：1714-1723.
[86]言志信,蔡汉成,王群敏,等.岩土体在地震作用下的破坏研究[J].煤炭学报,2010,35(10)：1621-1626.
[87]言志信,张刘平,曹小红,等.地震作用下锚固浅层红粘土边坡动力分析[J].岩土工程学报,2011,33(s1)：54-58.
[88]吴世明.土介质中的波[M].北京：科学出版社,1997.
[89]林淋.竖向地震动特征分析[D].哈尔滨：中国地震局工程力学研究所,2005.
[90]翟芳鹏.地震P-S波时差耦合作用的斜坡崩滑效应研究[D].北京：中国科学院地质与地球物理研究所,2009.
[91]刘晓.汶川地震区斜坡动力反应研究[D].成都：成都理工大学,2010.
[92]蔡汉成,言志信,王群敏,等.交通荷载作用下边坡动力响应数值分析[J].西北地震学报.2010,32(2)：220-224.
[93]Itasca Consulting Group. FLAC3D Version 4.0 (fast lagrangian analysis of continua in 3 dimensions) User's manual[M]. USA:Itasca Consulting Group Inc,2009.
[94]张鲁渝.用有限元强度折减法进行边坡稳定分析[J].中国工程科学,2002,4(10)：57—62.
[95]言志信,蔡汉成,王群敏,等.岩土体在地震作用下的破坏研究[J].煤炭学报,2010,35(10)：1621-1626.
[96]刘汉龙,费康,高玉峰.边坡地震稳定时程分析方法[J].岩土力学,2003,24(4)：553-556.
[97]Lin Meei-ling,Wang Kuo-lung.Seismic slope behavior in a large-scale shaking table modeltest[J].Engineering Geology,2006,86(2-3):118-133.
[98]唐洪祥,邵龙潭.地震动力作用下有限元土石坝边坡稳定分析[J].岩石力学与工程学报,2004,23(8)：1318-1324.
[99]周永江,王开云,符文熹,等.高地震烈度区堆积体边坡动力响应时程特征分析[J].山地学报,2007,25(1)：93-98.
[100]徐光兴,姚令侃,高召宁,等.边坡动力特性与动力响应的大型振动台模型试验研究[J].岩力学与工程学报,2008,27(3)：624-632.
[101]韩文峰.黄河黑山峡河段开发重大工程地质问题研究[M].北京：科学出版社,2004.
[2]陈建君.复杂山区斜坡的地震动力响应分析[N].成都：成都理工大学.2009.
[3]杜晓丽.地震荷载作用下岩质边坡稳定性研究[D].西安科技大学.2008.
[4]祁生文.岩质边坡动力反应分析[M].北京：科学出版社.2007
[5]Keh-Jian Shou, Cheng-Fung Wang. Analysis of the Chiufengershan landslide triggered by the 1999 Chi-Chi earthquake in Taiwan[J]. Engineering Geology,2003,68:237-250.
[6]张勇,许力生,陈运泰.2010年青海玉树地震震源过程[J].中国科学：地球科学：2010.7(40)：819-821.
[7]陈建君.复杂山区斜坡的地震动力响应分析[D].导师：许强.成都理工大学,2009.
[8]祁生文.边坡动力响应分析及应用研究[D].博士论文,中国科学研究生院,2002.
[9]刘红帅,薄景山,刘德东.岩土边坡地震稳定性分析研究评述[J].地震工程与工程振动,2005,8(1)：165-167.
[10]Leshchinsky D.,San K.Ching.Pseudo-static seismic stability of slopes[J].Design charts. J.of Geotech.Engrg,1994,120(9):1514-1532.
[11]Baker R.,Garber M.Theoretical analysis of the stability of slopes[J].Geotechnique,1978, 28(4):395-411.
[12]Taylor D.W.Stability of earth slopes[J].J.Boston Society of Civ.Engrgs,1937,24(3):197-246.
[13]姚爱军,苏永华.复杂岩质边坡锚固工程地震敏感性分析[J].土木工程学报,2003,(11)：34-37.
[14]Seed H.B.Stability of earth and rock fill dams during earthquakes in Embankment-Dam Engrg[M].Casagrande,Vol.(Eds Hirschfeld and Poulos),JohnWiley,1973.
[15]Newmark N.M.Effects of earthquakes on dams and embankments[J].Geotechnique,1965, 15(2):139-160.
[16]RobertoR.Seismieallyindueedlandslidedisplaeements:apredietivemode[J].Eng Geol,2000, 55(3):337-351.
[17]刘汉东.考虑地震历时影响的岩质边坡楔体稳定性分析与计算[J].华北水利水电学院院报,1991,(4)：35-40.
[18]薛守义,王思敬,刘建中.块状岩体边坡地震滑动位移分析[J].工程地质学报,1997,5(2)：131-136.
[19]黄建梁,王威中,薛宏交.坡体地震稳定性的动态分析[J].地震工程与工程振动,1997,(04)：114-123.
[20]王思敬.岩石边坡动态稳定性的初步探讨[J].地质科学,1977,(04)：372-376.
[21]王思敬,张菊明.边坡岩体滑动稳定性的动力分析[J].地质科学.1982,(2)：162-170.
[22]王思敬,薛守义.岩体边坡楔形体动力学分析[J].地质科学.1992,(2)：177-182.
[23]王勇智,戚炜,门玉明,等.大柳树坝址岩体松动机制物理模拟研究[J].西安科技大学学报,2004,24(3)：305—308.
[24]王勇智,戚炜,门玉明,等.强震区岩体动力破坏机制研究[J].岩土力学,2005,26(11)：1841—1844.
[25]郝建斌,门玉明,彭建兵,等.层状岩体边坡动力稳定性模型试验研究[J].公路交通科技,2005,22(6)：72—75.
[26]门玉明,彭建兵,李寻昌,等.层状结构岩质边坡动力稳定性试验研究[J].世界地震工程,2004,20(4)：131—136.
[27]许强,陈建君,冯文凯.斜坡地震响应的物理模拟试验研究[J].四川大学学报(工学版),2009,41(3)：266—272.
[28]徐光兴,姚令侃,高召宁,等.斜坡动力特性与动力响应的大型振动台模型试验研究[J].岩石力学与工程学报,2008,27(3)：624—632.
[29]许强,刘汉香,邹威,等.斜坡加速度动力响应特性的大型振动台试验研究[J].岩石力学与工程学报,2010,29(12)：2420—2428.
[30]艾畅,冯春,李世海,等.地震作用下顺层岩质边坡动力响应的试验研究[J].岩石力学与工程学报.2010,29(9)：1825-1832.
[31]董金玉,杨国香,伍法权,等.地震作用下顺层岩质边坡动力响应和破坏模式大型振动台试验研究[J].岩土力学.2011,32(10)：2977-2983.
[32]董金玉,杨继红,杨国香,等.基于正交设计的模型试验相似材料的配比试验研究[J].煤炭学报.2012,37(1)：44-49.
[33]杨国香,伍法权,董金玉,等.地震作用下岩质边坡动力响应特性及变形破坏机制研究[J].岩石力学与工程学报.2012,31(4)：696-702.
[34]杨国香,叶海林,伍法权,等.反倾层状结构岩质边坡动力响应特性及破坏机制振动台模型试验研究[J].岩石力学与工程学报.2012,31(11)：2214-2221.
[35]Hayashi H, N ishikawa J.Taniguch i K. Seismic behavior of road embankment[C]//Proc Intl Conf Centrifuge 98.Rotterdam:A. A. Balkema,l 998:243-248.
[36]M adabhushi S P QHaigh S K.Subedi B R.Seismic behavior of steep slopes[C]//Proc Physical Modeling in Geotechnics. Rotterdam:A. A.Balkema,2002:489-494.
[37]王年香,章为民.混凝土面板堆石坝地震反应离心模型试验[J].水利水运工程学报,2003(1)：18-22.
[38]Ng C.W.W.,Li Xiang-Song,Van Laak P.A.,et al.Centrifuge modeling of loose fill embankment subjected to uni-axial and bi-axial earthquakes[J].Soil Dynamics& Earthquake Engineering,2004,24(4):305-318.
[39]邓丽军,于玉贞,陈正发.应用于边坡动力离心试验的模型箱设计与试验检测[C].第24届全国土工测试会议论文集.郑州：黄河水利出版社,2005：246-252.
[40]Ishizaki and Hatekeyama.Consideration on the dynamical behavior of earth dams[R].Bul. No.52.Disaster Prevention Research Inst,Kyoto Unvi.,1963.
[41]陈育民,徐鼎平FLAC/FLAC3D基础与工程实例[M].北京：中国水利水电出版社,2008.
[42]Clough R.W.and Chopra A.K.Earthquake stress analysis in earth dams[J].J.Engrg.Mech., ASCE 1966,92,EM2.
[43]何蕴龙,陆述远.岩石边坡地震作用近似计算方法[J].岩土工程学报,1998,20(2)：66-68
[44]钱胜国,陆秋蓉.长江三峡船闸高边坡地震稳定性分析[R].武汉：长江科学院科研报告,1991.
[45]彭定,陈玲玲.三峡船闸高边坡在地震作用下的稳定性分析[J].长江科学院院报,2001,18(6)：47-49.
[46]徐卫亚,谢守益,徐瑞春,等.水布垭马崖高陡边坡三维数值模拟研究[J].武汉水利电力大学 (宜昌)学报,1998,20(3)：11-15.
[47]蔡德所.地震输入波形对三峡工程左岸坝段的动力影响研究[J].岩石力学与工程学报,1998,17(2)：183-187.
[48]刘立平,杨实君,李英民.软夹层参数对边坡动力特性的影响[J].重庆大学学报,2007,30(5)：31-34.
[49]狄仁东,石玉成,付长华.高边坡在水平动荷载下的动力响应规律研究[J].2007,23(2)：131-138.
[50]张学东,言志信,张森.ANSYS在岩质边坡动力响应分析中的应用[J].西北地震学报,2010,32(2)：117-121.
[51]徐志英,沈珠江.地震液化有效应力二维动力分析方法[J].华东水利学院学报,1981,No.3：1-14.
[52]徐志英,沈珠江.土坝地震孔隙水压力产生、扩散和消散的有限单元法动力分析[J].华东水利学院学报,1981,No.4：1-16.
[53]徐志英,沈珠江.尾矿高堆坝地震反应的综合分析与液化计算[J].水利学报,1983,No.5：30-39.
[54]徐志英,周健.土坝地震孔隙水压力产生、扩散和消散的三维动力分析[J].地震工程与工程振动,1985,5(4)：57-72.
[55]徐志英,周健.奥罗维尔土坝三维排水有效应动力分析[J].水利学报,1991,No.6：19-27.
[56]徐志英.奥罗维尔土坝三维简化动力分析[J].岩土工程学报,1996,18(2)82-87.
[57]周健,徐志英.土尾矿坝的三维有效应力的动力反应分析[J].地震工程与工程振动,1984,4(3)：60-70.
[58]刘波,韩彦辉,等FLAC原理-实例与应用指南[M].北京：人民交通出版社,2005.
[59]陈志坚.层状岩质边坡工程安全监控建模理论及关键技术研究[D]河海大学,2001,11
[60]许强,黄润秋.5.12汶川大地震诱发大型崩滑灾害动力特征初探[J].工程地质学报,2008,16(6)：721-729.
[60]许强,陈建君,冯文凯,肖锐铧,等.斜坡地震响应的物理模拟试验研究[J].四川大学学报(工程科学版).2009,41(3).
[61]肖克强,李海波,刘亚群,等.地震荷载作用下顺层岩体边坡变形特征分析[J].岩土力学.2007,28(8)：1557-1564.
[62]言志信,张森,张学东,等.顺层岩质边坡地震动力响应及地震动参数影响研究[J].岩石力学与工程学报.2011,30(增2)：3522-3528.
[63]孙博,彭宁波,王逢睿.云冈石窟第19窟西耳窟地震动力响应[J].西南交通大学学报.2012,47(4)：573-579.
[64]于翔,钱七虎,赵跃堂,郭志昆.地铁工程结构破坏的竖向地震力影响分析[J].解放军工大学学报(自然科学版),2001,2(3)：75-77.
[65]尹紫红.地震作用下的高速公路滑坡稳定性分析[J].重庆工学院学报(自然科学版),2007,21(9)：92-95.
[66]黄润秋.汶川8.0级地震触发崩滑灾害机制及其地质力学模式[J].岩石力学与工程学报,2009,28(6)：1239-1249.
[67]许强,黄润秋.5.12汶川大地震诱发大型崩滑灾害动力特征初探[J].工程地质学报,2008,16(6)：721-729.
[68]殷跃平.汶川八级地震地质灾害研究[J].工程地质学报,2008,6(24)：433-444.
[69]非明伦,樊跃新,崔建文,包一峰,陈坤华,卢吉高.施甸5.9级地震震害与竖向地震力作用[J].地震研究,2003,26(3)：292-298.
[70]胡晓燕,周健,胡晓虎.地震引起的竖向力对地铁隧道的影响[J].工程抗震,2006,2：32-35.
[71]梁庆国,韩文峰,赵士耀.垂直向地震作用对节理岩体失稳破坏的影响[J],西北地震学报,2007,26(4)：307-313.
[72]侯超群,王晓谋,石恒俊.地震对边坡稳定性影响分析[J].路基工程,2008,6：187-188.
[73]何蕴龙,陆述远.岩石边坡地震作用近似计算方法[J].岩土工程学报,1998,20(2)：66-68.
[74]崔芳鹏,胡瑞林,殷跃平,许强,张明.地震纵横波时差耦合作用的斜坡崩滑效应研究[J].工程地质学报,2009,(04)：455-462.
[75]崔芳鹏,殷跃平,许强,胡瑞林,张明.地震纵横波时差耦合作用的斜坡崩滑响应研究[J].岩土工程学报,2010,(08)：1266-1273.
[76]崔芳鹏,许强,谭儒蛟,殷跃平.地震动力作用触发的斜坡崩滑效应模拟[J].同济大学学报(自然科学版),2011,(03)：445-450.
[77]崔芳鹏,殷跃平,许强,胡瑞林.地震动力作用触发的斜坡崩滑响应研究——以青川东河口滑坡为例[J].三峡大学学报(自然科学版),2011,(03)：23-28.
[78]崔芳鹏,胡瑞林,殷跃平,许强,张明.纵横波时差耦合作用的斜坡崩滑效应离散元分析——以北川唐家山滑坡为例[J].岩石力学与工程学报,2010,(02)：319-327.
[79]方建瑞,许志雄,庄晓莹.三维边坡稳定弹塑性有限元分析与评价[J].岩土力学,2008,(10)：2667-2672.
[80]张向东,张玉,陈慧刚.黄土边坡二维和三维变形破坏模型的对比分析[J].中国地质灾害与防治学报,2012,(03)：18-20+38.
[81]袁伟,黄茂松,刘怡林,王浩然.含软弱夹层土坡稳定性的三维弹塑性数值模拟[J].郑州大学学报(工学版),2012,(05)：1-4.
[82]胡聿贤.地震工程学[M].北京：地震出版社,2006.
[83]Biggs J. M. Introduction to Structural Dynamics[M].New York:McGraw-Hill,1964
[84]万洪.岩质高边坡动力响应分析及其应用研究[D].武汉：武汉理工大学,2009.
[85]叶海林,郑颖人,黄润秋,等.地震边坡破坏机制及其破裂面的分析探讨[J].岩石力学与工程学报,2009,28(8)：1714-1723.
[86]言志信,蔡汉成,王群敏,等.岩土体在地震作用下的破坏研究[J].煤炭学报,2010,35(10)：1621-1626.
[87]言志信,张刘平,曹小红,等.地震作用下锚固浅层红粘土边坡动力分析[J].岩土工程学报,2011,33(s1)：54-58.
[88]吴世明.土介质中的波[M].北京：科学出版社,1997.
[89]林淋.竖向地震动特征分析[D].哈尔滨：中国地震局工程力学研究所,2005.
[90]翟芳鹏.地震P-S波时差耦合作用的斜坡崩滑效应研究[D].北京：中国科学院地质与地球物理研究所,2009.
[91]刘晓.汶川地震区斜坡动力反应研究[D].成都：成都理工大学,2010.
[92]蔡汉成,言志信,王群敏,等.交通荷载作用下边坡动力响应数值分析[J].西北地震学报.2010,32(2)：220-224.
[93]Itasca Consulting Group. FLAC3D Version 4.0 (fast lagrangian analysis of continua in 3 dimensions) User's manual[M]. USA:Itasca Consulting Group Inc,2009.
[94]张鲁渝.用有限元强度折减法进行边坡稳定分析[J].中国工程科学,2002,4(10)：57—62.
[95]言志信,蔡汉成,王群敏,等.岩土体在地震作用下的破坏研究[J].煤炭学报,2010,35(10)：1621-1626.
[96]刘汉龙,费康,高玉峰.边坡地震稳定时程分析方法[J].岩土力学,2003,24(4)：553-556.
[97]Lin Meei-ling,Wang Kuo-lung.Seismic slope behavior in a large-scale shaking table modeltest[J].Engineering Geology,2006,86(2-3):118-133.
[98]唐洪祥,邵龙潭.地震动力作用下有限元土石坝边坡稳定分析[J].岩石力学与工程学报,2004,23(8)：1318-1324.
[99]周永江,王开云,符文熹,等.高地震烈度区堆积体边坡动力响应时程特征分析[J].山地学报,2007,25(1)：93-98.
[100]徐光兴,姚令侃,高召宁,等.边坡动力特性与动力响应的大型振动台模型试验研究[J].岩力学与工程学报,2008,27(3)：624-632.
[101]韩文峰.黄河黑山峡河段开发重大工程地质问题研究[M].北京：科学出版社,2004.