碎石土路堑高边坡地震动力响应过程分析

2020-03-05吴红刚郭春香

中国地质灾害与防治学报 2020年1期

王飞，吴红刚，郭春香

(1.兰州交通大学土木工程学院，甘肃兰州 730070；2.西部环境岩土及场地修复技术工程实验室，甘肃兰州 730000；3.中铁西北科学研究院有限公司，甘肃兰州 730000；4.中国中铁滑坡工程实验室，甘肃兰州 730000)

0 引言

地震诱发的边坡滑动是主要的地震地质灾害类型之一，在山区和丘陵地带，地震诱发的滑坡往往具有分布广、数量多、危害大的特点[1]。不同的开挖、支挡顺序对边坡岩土体变形有很大的影响[2-3]，在施工中及时支护有利于开挖边坡的稳定性[4]，目前公共交通建设向施工难度更大的地区拓展，在高烈度地区高填深挖已成为常态，对高大边坡开挖一般都采用分级开挖及时支护的工序，因此有必要在考虑此施工顺序对其地震作用下的坡体变形及动力响应规律进行研究。

目前，开挖对边坡的影响主要集中在静力状态下开挖。夏熙伦等[5]根据试验观测资料，结合适当的有限元计算分析，研究探讨了三峡船闸高边坡岩体因开挖卸荷及应力调整而产生的松动区及其性状，进而为研究边坡的稳定性及加固方案提供了重要依据。姚裕春[4]用坡体开挖松弛区来描述坡体影响范围，认为在坡体开挖影响范围内，不同的稳定系数对应不同的松弛区，开挖后及时支护，能有效控制边坡的变形和限制下级边坡开挖后的应力松弛程度，对边坡的稳定性有利。赵尚毅等[6]利用FLAC3D研究了支护顺序对边坡变形及锚索受力的影响，表明及时支护(开挖一级支护一级)，能减小边坡的变形，但对锚索受力不利。肖春龙[7]用数值模拟得出结论，及时支护能有效保证边坡稳定，避免经济损失。王章琼等[8]分析某工地地质条件及开挖卸荷扰动的基础上,结合UDEC,总结了该边坡的变形破坏模式，并与实际监测进行对比进行了稳定性预测。曾芮等[9]以鄂西地区典型赵家岩崩塌为研究对象,通过现场调研和数值模拟分析等手段,对其变形破坏机理进行研究,认为开挖造成的应力重分布加速了坡体破坏。

在地震高度烈度地区，高边坡开挖及时支护施工完毕后，可能会受到强烈地震作用，如云南功东高速大营盘深路堑高边坡施工完毕后(2018年4月)即遇到通海5.0级地震(8月13日)、墨江5.9级地震(9月5日)，但地震作用下考虑开挖及时支护工序的边坡变形规律和动力响应还无人研究。故本文以功东高速大营盘深路堑高边坡为原型，在静力开挖及时支护工序的基础上，对自然和降雨工况下的地震动力响应进行对比分析，以期得到高边坡工程的变形规律及动力响应特征。

1 工程概况

1.1 研究区坡体形态

功东高速大营盘路堑区该坡体地势西南高东北低，地貌类型属中低山侵蚀剥蚀堆积地貌，地形起伏较大，山体坡向63°，自然坡度为20°～45°。坡型主要受构造侵蚀、风化堆积作用控制，图1为边坡地质横断面图。

图1 边坡地质横断面图Fig.1 Geological cross-section of slope

1.2 地层岩性特征

根据钻探揭露，路堑区揭露地层为第四系人工成因种植土、第四系中更新统冲洪积粉质黏土和碎石土。各个岩层特征自上而下描述如下：

(1)种植土：结构松散，潮湿，主要以黏性土为主，局部有少量碎石，含量约10%。分布在坡体表层，厚度为0.7～1.3 m。

(2)粉质黏土：该层土体呈硬塑状态，土质不均，切面粗糙，局部含砾石，粒径在10～20 mm，含量约占20%。层厚在3.0 m以内。

(3)碎石土：土体密实，呈稍湿状态，粒径大于20 mm颗粒质量约占总质量的75%，最大粒径约100 mm，颗粒磨圆度较差，呈棱角状-次棱角状，级配较好，母岩成份以砂岩和灰岩为主，隙间由泥砂质充填，局部含块石，层顶埋深0.7～24.2 m，厚度巨大，勘察未揭穿。

1.3 开挖支护方式

此路堑边坡开挖后左侧形成边坡最大高度78.9 m，边坡坡向约51°，属于高边坡，在施工时采用六级放坡，及时支护的施工顺序，即在开挖后立即对所开坡面进行支护，按照设计方案，坡体中部锚索较长，坡顶和坡脚锚杆较短：第一级坡和第六级坡锚杆分别为20 m和12 m；第二、五级坡锚索全长25 m，其中锚固段长10 m；第三、四级坡单根锚索全长30 m，锚固段长10 m。图2为此深路堑高边坡的支护方案示意图，坡体基本特征、锚索(杆)位置和长度示意图如图2所示。

图2 研究断面简化图(单位：m)Fig.2 The studied cross-section

1.4 研究区地震动参数

根据《中国地震动参数区划图》GB l8306—2015，项目起点功山至东川段地震动峰值加速度≥0.40g(对应的地震基本烈度为9度)，地震动反应谱特征周期为0.40 s，设计地震分组为第二组。

2 建立计算模型

2.1 地质模型简化

按照FLAC3D建立模型和计算的需要，保留边坡的真实尺寸分布，对其他附属结构如排水沟、表层土等不影响结构计算的部分忽略，化简后的边坡断面如图2所示。

根据实际的地形和结构物的布置位置，为了防止边界效应对研究区产生影响，对模型计算范围进行放大，以坡体为中心，将坡顶和坡脚各延伸一倍的坡高，取模型的长为180 m，高为98 m，宽7.2 m，遵循实际的尺寸为原型，进行相应的简化，得到三维地质模型。

锚索和框架梁采用结构单元cable和beam模拟。横纵梁的beamsel之间的link定义为固接(6个自由度全部设置为rigid)，同样的cablesel 与beamsel之间也为固接。开挖后设置支护结构的数值计算模型示意图如图3所示。

图3 开挖完成后支护结构示意图Fig.3 Schematic diagram of support structure after excavation

2.2 网格精度要求

根据Lysmer和Kuhlemeyer(1969)的研究，动力计算中要想精确的展示模型中波的传播，建模时网格尺寸必须小于输入波最高频率对应的波长的1/8到1/10[10]。

本次滤波后最大频率为10 Hz，由根据钻探及波速测试成果中钻孔波速测试结果可知，场地等效剪切波速为269.0 m/s，则最小网格尺寸应小于2.69 m；再根据《公路工程抗震设计规范》，硬土的剪切波速为250～500 m/s，故网格尺寸最小应当小于2.5 m。

本模型网格最小尺寸为2.5 m，满足要求。

2.3 参数取值

通过原位直剪试验设备，对大营盘碎石土进行水平推剪实验根据实验所测的推剪应力和剪切位移，计算得到大营盘强度参数(C、φ)如表1所示，斜杠后为饱和度70%降雨后的强度参数。

表1 岩土体强度参数取值Table 1 Strength parameters of rock and soil mass

本次采用新型玄武岩纤维复合(BFRP)筋是以玄武岩纤维为增强材料，以合成树脂为基体材料，并掺入适量辅助剂，经特殊的工艺处理形成的一种新型非金属复合材料，其轻质高强、耐腐蚀、抗疲劳、环保等优良的工程特性使之在支护工程中崭露头角，可行性和优越性已经被很多学者所证明[11-12]。锚固只考虑锚索的轴向作用，在锚索的锚头部位设置大的黏结参数以模拟锚头的作用，锚索材料参数选择BFRP厂家提供的数据，框架梁的截面为0.4 m×0.4 m的正方形，具体取值见表2。

表 2 锚索和混凝土框架梁参数取值Table 2 Parameters of the anchor cable and concrete frame beam

2.4 监测点及工况设置

在开挖时，每次开挖两级坡，每开挖一次，就将开挖的两极坡面进行锚索框架梁支护。

为研究开挖及时支护对边坡变形的影响，本次模拟分两个阶段进行，第一阶段为静力阶段：分级开挖及时支护具体工况设置如表3所示。第二阶段为地震作用下的动力计算：将静力计算的最终结果(工序3和4)，分别作为动力计算的自然和降雨工况的初始状态，详细的计算工况如表4所示。

在每级坡面中点设置坡面监测点，在内部设置一纵、二纵测点(测点为红色圆圈)(图2)。

表3 施工顺序(静力工况)Table 3 Construction sequence (static condition)

表4 地震工况设置Table 4 Setting of seismic conditions

2.5 地震波输入

动力计算时取地震波为实测鲁甸加速度波，该波持续时间为50 s，加速度峰值为1.2 m/s2。截取其中能量最强的12 s作为输入波，且只研究破坏力强的剪切波。高频部分振幅很小，故滤去10 Hz以上的地震波，以减小网格尺寸。处理后的输入波时程曲线如图4所示。

图4 位移和输入波时程曲线Fig.4 Time-history curve of displacement and input wave

2.6 动力条件设置

在动力计算时，模型边界条件和阻尼的选取对计算结果影响很大。FLAC3D中提供了安静边界和自由场边界来减少模型边界上波的反射[13]。此次模拟四周采用自由场边界(图5)，底面设置安静边界使边坡基础部分置于一个半空间无限的自然地质体内，以更接近实际土体状态。

图5 设置自由边界的模型Fig.5 Model with free boundaries

FLAC3D动力计算中提供了三种阻尼形式，瑞利阻尼、局部阻尼和滞后阻尼[14]。本文采取在模拟中较符合实际的瑞利阻尼进行计算。确定瑞利阻尼需要两个参数：最小临界阻尼比和最小中心频率，对于岩土材料来说，临界阻尼比一般为0.02～0.05[13]。最小中心频率取地震波的卓越频率3.3 Hz，两个参数分别为0.05、3.3。

3 结果分析

3.1 静力结果分析

通过对各个工况的结果(位移、塑性区和剪切应变增量云图(篇幅原因不在列出))分析可得，开挖主要是对开挖面一定范围内有影响，且主要是由于回弹变形所引起的位移，对稳定性影响不大。

图6为坡面测点按施工顺序的水平向位移时步曲线，可以看出，每次开挖后，基本上都有向着坡内弹性回弹的负位移产生，随着开挖深度增加，坡体负位移增加量逐渐降低，尤其是下部(第一、二和三级坡)坡面负位移增加量降低幅度较大，在第二次开挖时，增加量约有坡体上部(第四、五和六级坡)的1/2，第三次开挖时，回弹的负位移增加更少，甚至在坡脚有向下坡方向位移的趋势。在降雨时，坡面各测点均开始向下坡方向位移，且越向坡脚此趋势越明显。

图6 各测点水平向位移的时步曲线Fig.6 Time step curve of horizontal displacement

表5为每个工况下相对于上个工况的水平负位移增加率，则负值表示向下坡方向的正位移增加量。可以看出每次开挖坡体中部第三、四级坡为负位移增加率即弹性回弹最大的部位。在开挖支护完成(第三次开挖支护)后，下滑的趋势从坡脚最先开始，在降雨后陡增并逐渐向上部迁移，增加率逐渐降低。

表5 负位移增加率Table 5 Negative displacement increase rate

3.2 地震作用下结果分析

3.2.1剪应变分析

在坡体中剪应变增量较大的部位，是土体受剪切破坏最严重的部位，在不利条件下可能沿着这一剪切带发生失稳破坏。

图7为地震作用下的剪应变增量云图，可以看出在天然工况下，剪应变增量最大值为1.915 8×10-2，剪切带在顶部没有贯通，形成坡中和坡脚两个剪出口，而在降雨工况下，剪应变增量最大值增大一倍，剪切带贯通，从一级坡平台剪出。可以认为在自然工况下，坡体整体强度较高，在地震作用下会在坡中部产生局部剪出口，倾泄地震动能量，有利于坡体保持整体稳定性。而在降雨工况下，土中剪应变在剪切带迅速增加，继而形成贯通滑带，产生类似于滑溃型的地震失稳模式[15]。

图7 地震作用下剪应变增量云图Fig.7 Incremental shear strain nephogram under earthquake

3.2.2位移响应

在地震作用下，主要分析可以表征坡体稳定的水平XX向位移，图8为两种工况下的水平向位移云图，在两种工况下，位移变化规律相似：在坡顶位移量较小且方向平行于坡面，到坡体中下部坡体位移逐渐增加至最大且呈抛出状，到坡脚时又开始收敛，方向逐渐平行于坡面。主要不同体现在降雨工况下的位移迅速增加，最大值为自然工况下的2倍，达到42 cm。

图8 地震作用下水平向位移云图Fig.8 Horizontal displacement cloud map under earthquake

图4为一级坡降雨工况下的位移和输入波时程曲线，其他点在地震作用下变化规律相似，只是大小有所差异，故取次点加以说明。可以看出，位移是由加速度幅值较大的地震波引起的，且最大位移波动滞后最大加速度幅值T=1.387 s，在此之后，在较小振幅下引起了比最大振幅所引起的更大的位移。

图9为各测点在动力计算过程中的最大值，图中百分比为各坡面测点水平位移峰值在降雨工况下较自然工况下的增加率，可以看出在两种工况下，均为中部坡体位移最大，在降雨工况下，越往坡体上部，位移增加率越大：从一级坡的102.6%到六级坡的154.8%，最大位移也从自然工况的二级坡最大到降雨工况上移至三级坡。

图9 地震作用下坡面最大位移Fig.9 Maximum displacement at slope surface under earthquake

3.2.3加速度响应

图10为坡面(P)、一纵(Z1)和二纵(Z2)测点在两种工况下沿高程的加速度峰值曲线。P测点在两个工况下，除了坡顶均沿高程整体呈减小趋势。Z1测点在高程为55 m的中部开始降低，总体呈“V”字型。在Z2测点，沿高程增加方向虽有小幅波动，但整体呈沿高程增加的趋势，且两种工况的不同在坡内不明显，到坡顶逐渐变得清晰。分析可得在地震中，离开挖面越近，两种工况在加速度相应中差别越大，即受降雨影响越大。与位移图8对比，发现位移较大部分的加速度相应会降低，越往位移小的坡体内部(Z2)沿高程放大越明显，降雨的影响越小。其原因可能是土体内部裂缝增加，土体刚度降低，阻尼增加，从而导致其加速度响应降低。