滑塌式危岩破坏断裂力学探析.doc

滑塌式危岩破坏断裂力学探析 　　摘 要：危岩是一种全球性高频率的地质灾害，其中滑塌式危岩破坏是危岩主要的破坏类型之一。该文针对滑塌式岩石破坏类型，通过考虑自重、裂隙水压力和地震力3种组合荷载构建滑塌式危岩力学模型，采用断裂力学的方法推导其Ⅰ-Ⅱ型复合应力强度因子。根据最大周向应力理论，使用三角形万能公式计算变形得出其理论断裂角 关键词：断裂力学 Ⅰ-Ⅱ型复合应力强度因子 断裂角 滑塌式危岩 中图分类号：TU45 文献标识码：A 文章编号：1674-098X（2017）01（c）-0009-0 岩石是由一种或多种矿物在地质作用下天然产生的复杂结构体。大多岩体工程中的岩石属于压剪状态，其岩石裂缝的应力场应为压剪应力场。滑塌式危岩破坏的本质是压剪破坏，因此对于滑塌式危岩破坏的研究具有必要性和迫切性。众多国内外学者都致力于对岩石破坏的机理进行了相应的研究。Nara Y等[1]以处于亚临界裂纹扩展状态的花岗岩为研究对象，他发现裂隙水将加速花岗岩裂纹扩展的速度并影响花岗岩的结构强度。Zygouri V等[2]对Skolis山和Acrocorinthos的岩石进行了研究，研究表明浅层地震会引起大范围的岩石崩塌。Chen H K等[3]提出了处于激励效应下的不稳定岩石的破坏准则并建立了对其安全性的评估方法。Johari A等[4]采用联合分布随机变量的方法对处于临界状态的岩石的稳定性进行了可靠的评估。Li Y等[5]使用FLAC3D软件模拟了节理裂隙水压力作用下岩体的裂隙发育，结果表明裂隙水对节理岩体的强度和稳定性有明?@的降低作用。Liang L等[6]在龙马溪地区对页岩进行取样，研究表明水性液体对页岩形成的裂纹扩展具有显著的正向影响。迄今为止，学者多以试验和数值分析的方法对岩石裂缝进行研究，对于岩石裂缝的理论分析略有欠缺。该文采用断裂力学的方法对滑塌式危岩进行理论推导，其结果对于防灾减灾和工程安全评估等方面具有一定的理论意义和经济价值 1 应力分量的坐标变换 如图1所示，建立滑塌式危岩模型。其中为危岩体重心，为主控结构面的孔隙水压力，为单位长度水平地震力，为单位长度竖向地震力，为单位长度的岩体重力 在图1中的裂纹尖端选取单元体，在岩石自身重力作用下形成压剪应力（如图2） 在单元体平行于轴的直线上，设其外法线与坐标轴的夹角分别为、，则有： ，，， 将已知参数带入边界条件方程： （1） （2） 可得： ， 在单元体平行于轴直线上，同理可得：，。综上可得： （3） 建立新、旧坐标系如图3所示，旧坐标系的应力张量为，新坐标系的应力张量为。新坐标系轴在旧坐标系2个坐标轴上的投影（即与2个坐标轴的夹角余弦）分别为、；新坐标系轴在旧坐标系2个坐标轴上的投影（即与2个坐标轴的夹角余弦）分别为、 推导已知： （4） （5） （6） 将（3）（4）（5）代入（6），则有： （7） 由上式应力张量可知： （8） 2 推导Ⅰ型应力强度因子 设，，根据Westergaard提出的应力函数[7]为： （9） 因为，，所以： （10） 同理可得： （11） （12） 将Westergaard应力函数偏微分，可得： （13） 同理可得： （14） （15） 将（13）（14）（15）整理，可知： （16） 如图4所示，在无限体内有一个长为的中心贯穿裂纹，在其无穷远处受双向均匀应力作用。其边界条件为：当，a时，；当， 由（16）式可知，当，时，；当，时，。可将其应力函数假设为： 由于其裂纹关于轴对称，因此假设其应力函数为： 由于上述应力函数是在，的特殊情况下推导的。对于大多数情况下，因此可将上式中用代替，则有： （17） 将原坐标点移新坐标点，设新坐标系中任意一点的复数坐标为，则两坐标系的换算关系如下： （18） 即： （19） 将（19）式带入（17）式： （20） 当时，令 （21） 由于当→0时，为一常数。因此，可令常数为： （22） 联解（20）（21）（22），可得： （23） 3 推导Ⅱ型应力强度因子 设，，根据Westergaard提出的应力函数[7]为 （24） 因为，，所以： （25） 同理可得： （26） （27） 将Westergaard应力函数偏微分，可得： （28） 同理可得： （29） （30） 将（28）（29）（30）整理，可知： （