地震作用下呷爬滑坡变形特征分析.docVIP

地震作用下呷爬滑坡变形特征分析 　　摘要：运用有限差分软件FLAC，建立了呷爬滑坡体三维模型，分别施加1倍和3倍El-Centro地震波。计算结果表明，1倍地震波作用下会发生一定的永久位移，但不会滑动失稳；3倍地震波作用下将发生滑动失稳。进一步发现边坡在1倍和3倍地震波激励下的变形滑移特征不同，第一种情况属于累积效应阶段，变形主要集中于滑体O，第二种情况属于触发效应阶段，滑体N的滑动变形同样剧烈。文中从滑坡体的组成物质和材料参数、滑动区域的地形因素和滑动区域的动力响应三方面进行了分析。运用Newmark有限滑动位移法，基于滑体N和滑体O的加速度响应时程，计算了两者的永久位移，结果与数值计算吻合良好。表明滑体O处的动力响应会产生更大的永久位移，从而验证了数值计算的结果。同时，基于结构动力响应而计算得到的永久位移也可以作为衡量地震作用对该处结构变形影响的指标，从而估计地震对局部结构的破坏作用。 　　关键词：呷爬滑坡；变形特征；Newmark 方法；永久位移；加速度响应 　　中图分类号：TU443 文献标识码：A 文章编号：1672-1683（2017）05-0110-06 　　对于人口分布稀疏的山区，地震触发的滑坡、山体崩塌以及泥石流等灾害所造成的损失往往比地震造成的直接损失更大。对于20世纪末到21世纪初发生的台湾集集地震、汶川地震、雅安芦山地震，均伴随多次大规模滑坡和其它次生灾害，其人员伤亡和财产损失不亚于一场地震。因此，地表运动激励下的滑坡体稳定性分析逐渐被岩土工程界和地震工程界所重视。近年来，专家学者对于地震作用下边坡的失稳机制和动力特性进行了大量研究。何蕴龙等通过有限元对滑坡体进行了动力分析，可以近似计算滑坡所受到的地震作用；祁升文等对各种情况下滑坡的动力响??进行了数值模拟与分析，总结出了高、低边坡的动力反应形式；言志信等对于顺层岩质滑坡体地震激励下的加速度响应规律、频谱特性、破裂面发生位置以及失稳判定条件进行了研究，并对耦合地震波作用、多级平台和黄土地质等情况进行了较为全面的研究；徐光兴等结合Newmark方法，提出了根据加速度预测边坡永久位移的数学模型；张雪东等运用FLAG3D有限差分商业软件分析了锦屏库区呷爬滑坡在静力作用下的变形特征和稳定性；殷跃平等结合大量的现场资料，分析了大光包巨型滑坡的失稳机制和动力响应。本文根据呷爬滑坡体变形监测图及地形图，建立了滑坡体模型，运用FLAC对其动力作用下的变形特征进行了分析，并结合New-mark数值计算法利用边坡地震响应计算了边坡的永久位移。 　　1模型建立 　　1.1滑坡体地质概况 　　呷爬滑坡位于锦屏水电站库区上游右岸，距离坝址约11.5 km。雅砻江流经该滑体处的河道较顺直，流向约为N12°E。该处河谷宽阔，左岸山体较为平缓，并有河漫滩、阶地发育；右岸山体呈反向坡，岸坡的坡度约为35°。 　　雅砻江岸坡岩层产状在呷爬滑坡处呈N 10～30°E/SE ∠85°，岩层主要由三叠系杂谷脑组深-灰黑色粉砂质板岩夹变质细砂岩组成。 　　呷爬滑坡体的平面投影呈长条状，具有完整的圈谷地貌。滑体靠近雅砻江一侧高程较低，为l 655 m，与枯水期雅砻江水位相同，靠近山体一侧高程较高，为2 120m，滑体前后缘高差为465 m。滑坡体的纵长约为880m，宽度在260到300 m之间，表面面积约为Q 28 km2，滑坡残体体积约为1 300万m3。 　　滑坡表面可见3级平台，3级平台高程分别为1775～1840 m、1880～1925m、2010～2100 m，坡度为17～18°，原为耕地。在第一和第二级平台中发育有两条纵向冲沟。其中，一条位于滑坡轴线的北侧（1号沟），另一条位于滑坡轴线的南侧（2号沟），南侧冲沟比北侧冲沟更为深长。 　　在滑坡体临江侧到中部，表层岩体主要由块碎石土组成，块碎石土的成分为粉砂质板岩、泥岩和砾岩；在滑坡中部到靠近山体一侧，表层岩体主要由孤石，块碎石土，解体不完全的粉砂质板岩、变质细砂岩和砾岩等组成；滑体物质分层明显，据钻孔揭露，孔深27.97～39.57 m以上，为褐黄、灰黄色块碎石土，以下至77.08～81.80 m为破碎的岩石碎块夹黏土，局部岩体保留有原始层面。 　　在滑坡临江侧和两条冲沟内，可见明显的滑带出露。出露的滑带在剪出口处反翘，根据地址勘察结果，其产状为N15°E/SE∠22°。滑带的主要组成物质为灰黑色泥夹碎石并带有大量石英颗粒，其中碎石成分为炭质板岩、变质细砂岩等，滑带在滑体范围内平均厚度为1～3 m，局部可达8 m，上述滑带组成物质在干燥时刚度较大，含有一定水分后刚度显著减小。 　　1.2建立数值模型 　　模型建立是FLAC3D的短板，本文采用先在CAD中生成面域，导入AN SYS划分网格，再利用清华大学郑文