论文部分内容阅读
地震荷载作用下块状结构岩体边坡的稳定性研究是我国西部强地震区大型水利水电工程建设、山区高速公路工程、铁路客运专线以及核电工程建设中需要解决的难点和热点问题。对边坡在地震作用下的动力响应规律、失稳机制、稳定性评价方法等科学问题进行系统深入的研究,不仅能推动地震边坡稳定性理论的发展和创新,而且能够为岩土工程设施的规划、设计和施工提供必要的技术支持,具有重要的理论意义和现实意义。岩体是由结构面网络及其所围限的岩石块体所组成。结构面对地震波在岩体中的传播过程及规律有着很重要的影响。研究地震波在节理岩体中的传播特征是分析岩体边坡动力响应和动力稳定性的基础。边坡的地震动力响应包括加速度、速度、位移、动应力和动应变响应等等多方面的内容。广义上,除了上述这些标准参量外,某些抽象的、人为定义的变量,例如动态的稳定系数,也可以视为边坡对地震动力的响应。在所有的这些能够表征边坡动力响应的物理量、信号量中,加速度无疑是居于基础性地位。这主要有如下两个方面的原因:第一,速度和位移可分别由加速度对时间积分和二次积分获得,因此不具有独立性。第二,边坡体内任意微元的动应力状态与其加速度、质量之间的关系符合牛顿第二定律,并且动应力与动应变通过模量建立联系,从这个意义上来说,动应力和动应变对加速度也具有依存关系。事实上,边坡的各种动力响应皆由地震加速度引发。因此,研究加速度的动力响应规律(时间和空间分布)已成为这一领域的基本问题之一,同时对于评价边坡的动力稳定性也有着重要的价值。在岩体边坡动力稳定性计算中,永久位移分析和动力稳定性系数分析方法是两种主要的动力稳定性评价方法。永久位移标准较稳定性系数标准能够更为准确地反映边坡在地震工况下的动力响应,在理论上永久位移标准比单一的稳定性系数更为合理,但是永久位移法却存在诸多不足,首先从本质上看永久位移法是一个较为粗略的评价方法,在计算边坡地震永久位移时不能考虑边坡材料自身的变形、动力条件下岩体的强度以及变形特性;其次,永久位移法计算的地震永久位移其实是平均位移,不能反映具有复杂节理岩体在地震作用下真实的位移分布情况;另外,永久位移法无法得出地震荷载作用下边坡的稳定性系数,这在实际工程中进行边坡稳定性评价是极为不便的。在此背景下,本文将尝试通过将永久位移法与动力强度折减法相结合,以永久位移作为边坡的动力失稳评价指标,并通过强度折减法来计算得出边坡的稳定性系数。由于强度折减法物理意义明确,计算过程简洁,因此国内外很多学者采用强度折减法对边坡的稳定性系数数进行了分析。目前,主要存在以下两种边坡的动力失稳判据:(1)从关键点永久位移和折减系数的关系曲线来判断。(2)以地震荷载结束后关键点永久位移的时程曲线是否发散来判断。然而就目前国内外的研究进展来看,上述两种失稳判据均有它们的不足之处,本文将在上述动力失稳判据的基础上分别提出新的计算方法。通常情况下,岩体质量较好、岩体结构完整、倾向坡外的结构面未发育或较少发育的块状岩体边坡整体稳定性较好。但在强震作用下,岩体震裂松动、结构面加剧扩展、结构面抗剪强度显著降低,该类岩体边坡的整体稳定性也大大降低。目前,关于该类岩体边坡地震动力响应以及地震条件下的失稳破坏方式方面的研究尚少。鉴于地震的强大破坏力,正确分析强震作用下块状岩体边坡的动力响应及失稳破坏方式则显得尤为迫切和重要。本文主要针对块状结构岩体边坡,采用三维离散单元法(3DEC)对其在地震荷载作用下的动力响应规律、稳定性分析以及破坏方式进行了深入研究。目的是揭示地震荷载对块状结构岩体边坡稳定性的影响机制,为更准确地分析块状结构岩体边坡的动力稳定性提供有效的理论依据。论文完成的主要工作包括以下几个方面:(1)基于三维离散元软件3DEC,建模分析了地震波在含单个节理岩体中的传播特征。首先,建立简谐地震波在均质岩体中传播的数值模型,通过参数研究选取了合理的建模参数和动力分析参数,包括边界条件、网格单元尺寸、动力时间步等,为选取合理的参数应用于动力学问题方面的研究提供了依据。接着建立地震波在单个节理岩体中的传播模型,并分别考虑了地震波垂直入射结构面和地震波倾斜入射结构面时结构面刚度对波传播规律的影响,并结合前人的理论研究成果,探讨数值模拟方法的可行性,结果表明采用离散元软件3DEC来分析节理岩体的动力问题是非常可靠的。(2)通过密集的数值模拟计算,对地震作用下各结构面因素对岩体边坡动力响应的影响规律进行了研究。岩体边坡地震动力响应分析是一个十分复杂的研究课题。除了地震荷载的复杂性之外,岩石材料的性质、结构面的物理力学特性及其在岩体内的分布和规模等因素,都将对岩体边坡的地震响应产生一定的影响。如果要全面地考虑上述因素的影响,将使得这一问题变得相当复杂以至于难以求解。因此为了突出重点,本文主要研究了结构面物理力学特性以及其在岩体边坡内的分布特征对岩体边坡动力响应的影响。鉴于结构面的刚度对地震波的透射系数和反射系数有着重大的影响,在结构面物理力学特性方面主要考虑结构面刚度对岩体边坡动力响应的影响。此外在结构面分布特征方面,主要从结构面的产状、结构面的起始位置以及结构面的密度这几个角度来分别进行研究。研究结果表明,①、结构面产状的变化会使得结构面对地震波的反射方向和折射方向发生变化,进而影响地震波场能量在坡体中的分布情况,最终体现为边坡PGA放大系数等值线分布的空间变化。对于顺倾结构面边坡,当结构面倾角较小时,结构面对地震波的传播起着一定的衰减作用,而当结构面倾角较大时,则起着相反的作用。边坡的动力稳定性随着结构面倾角的增大而不断变差。逆倾结构面的存在会对地震波在岩质边坡中的传播产生衰减作用,并且随着结构面倾角的增大而不断弱化。相比于顺倾结构面边坡,逆倾结构面边坡的动力稳定性大体上要强于前者。②、结构面的起始位置对于岩体边坡的地震动力响应有着比较明显的影响,其位置越高,岩体边坡的地震动力响应就越强,反之则越弱。③、随着结构面刚度的增大,岩体边坡的地震动力响应的强度呈现出明显的递增趋势。结构面刚度对各反射波和透射波的影响作用主要体现在它们的能量分配关系上,而对于它们在岩体边坡中的传播路径则没有什么影响。当入射波为剪切波时,岩体边坡地震动力响应对结构面剪切刚度的敏感性要明显地高于其法向刚度。④、结构面发育愈密集,边坡岩体的完整性愈差,岩体边坡的地震动力响应就越强烈,其动力稳定性也就越差。(3)运用正交试验设计法,通过数值计算与数理统计分析评价了岩体边坡PGA放大系数极大值影响因素的敏感性。虽然PGA放大系数极大值并不能全面地反映岩体边坡地震动力响应的整体情况,但其能在一定程度上诠释岩体边坡地震动力响应的强度问题。在分析各种结构面因素单独作用时PGA放大系数极大值的变化规律的基础上,综合考虑各种结构面因素之间的相关性和随机联动性,分别采用极差分析法和方差分析法对岩体边坡PGA放大系数极大值影响因素的敏感性进行了排序。参与分析的5个结构面因素中,敏感性由大到小的顺序依次为:结构面倾角、结构面起始位置、结构面剪切刚度、结构面法向刚度、结构面间距。因此,当岩体边坡中发育有倾角较陡的结构面且结构面发育深度较浅时,其动力响应的强度往往就会越强,在地震荷载的作用下,这类岩体边坡的动力稳定性问题也将显得更为突出。(4)针对动力强度折减法中第一种失稳判据在实际应用中的局限性,提出了“基于永久位移比的岩体边坡动力稳定性计算方法”。地震荷载作用下,地震永久位移随强度折减系数的增加而增加,且在强度折减系数较小时增幅不明显,当超过某一定值时,地震永久位移的增幅迅速增大,呈明显的指数函数形式变化。关于如何在地震永久位移与折减系数的关系曲线中标定出边坡的动力稳定性系数,目前的解决方法主要是以该曲线“拐点”即永久位移值突变点所对应的折减系数来确定,然而对于“拐点”的位置该如何确定则存在较大的人为主观因素,并且缺乏合理的数学理论公式做支撑,使得该方法在使用上存在着一定的局限性。在此背景下,本文提出永久位移比理论新概念,即当岩体边坡沿滑动面的地震永久位移达到滑动面长度的一定比值时认为该边坡失稳破坏。将此方法运用于汶川地震区块状岩体边坡实例研究中,计算得出该实例边坡在汶川地震荷载作用下的动力稳定性系数为1.04,稳定性计算结果与实际地质调查情况相符,验证了本方法的工程实用性。(5)针对动力强度折减法中第二种失稳判据的不足之处,提出了“考虑结构面退化的岩体边坡动力稳定性计算方法”。当采用第二种失稳判据来分析岩体边坡的动力稳定性时,通过结构面静态强度参数计算所得的岩体边坡动力稳定性系数与其静力稳定性系数进是相等的,地震荷载的振幅对边坡动力稳定性系数的影响得不到任何体现,在此背景下,基于前人有关结构面摩擦系数随累积位移和相对运动速率不断衰减的研究成果,提出了考虑结构面退化的岩体边坡动力稳定性计算方法,给出了相应的实施流程,并用汶川地震区的工程实例进行了说明,计算结果表明,使用该方法来计算岩体边坡的动力稳定性是可行的。(6)老虎嘴岩体边坡是汶川地震区典型的受优势结构面控制的块状岩体边坡,5·12汶川大地震使得该边坡内部的结构面贯通,形成楔形体在地震作用下向临空面方向运动,以至产生严重的崩滑现象,并导致岷江堵塞和交通封闭。以老虎嘴块状结构岩体边坡的地质原型为基础,应用三维离散元数值模拟软件3DEC,建立了老虎嘴块状结构岩体边坡的三维离散元数值计算模型,并通过动力计算还原了该边坡在汶川地震荷载作用下变形破坏的演化过程,对其动力失稳机制进行了详细分析,为研究具有相同结构类型的岩体边坡在强震作用下的破坏方式提供了参考依据。