[发明专利]一种岩石变形破坏模式预测方法

摘要

本发明公开了一种岩石变形破坏模式预测方法,包括以下步骤：S1：将三轴应力作用下的岩石爆破抽象为一个含球形空腔的无限弹性介质模型的动静加载过程：S2：测量地下原岩应力场，获得地应力场的三个主应力分量σx0，σy0，σz0；S3：测定爆破对象岩石的密度ρ、泊松比μ以及弹性模量E；S4：根据岩石的三轴试验结果，获得岩石在不同围压下的屈服应力；解决了现有技术中设备费用高、有一定危险性、试验不方便、非三轴应力状态、软件价格高、使用门槛高和模拟步骤繁琐的问题。

主权项

1.一种岩石变形破坏模式预测方法,其特征在于，包括以下步骤：S1：将三轴应力作用下的岩石爆破抽象为一个含球形空腔的无限弹性介质模型的动静加载过程；S2：测量地下原岩应力场，获得地应力场的三个主应力分量σ_x0，σ_y0，σ_z0；S3：测定爆破对象岩石的密度ρ、泊松比μ以及弹性模量E；S4：根据岩石的三轴试验结果，获得岩石在不同围压下的屈服应力；S5：计算对应的应力张量第一不变量I₁以及应力偏张量第二不变量J₂，以‑I₁值为横坐标，对应的为纵坐标，在平面直角坐标系中标记不同围压下的屈服应力状态对以上散点数据进行拟合，结合岩石单轴拉伸试验或巴西圆盘试验测定的抗拉强度σ_t，获得岩石的完整屈服面；S6：根据地应力场的三个主应力分量σ_x0，σ_y0，σ_z0，首先计算在σ_z0作用下的静态应力场，根据球坐标系内弹性介质中的点的应力状态计算方法，计算出S7：在σ_y0，σ_x0分别作用下的静态应力场，根据球坐标系内弹性介质中的点的应力状态计算方法，计算在σ_y0，σ_x0分别作用下含球形空腔的弹性介质中的应力状态S8：将分别由σ_x0，σ_y0，σ_z0产生的应力场叠加起来，得到在远场三轴地应力作用下的静态应力分布S9：对于在球腔壁上的爆炸荷载作用下的岩石动态响应，根据胡克定律，写出各应力分量与径向位移之间的关系式：S10：根据达朗贝尔原理，写出球面波传播的控制方程、边界条件以及初始条件：式中，φ(r,t)是位移势函数，t是弹性波传播的时间，c₁是纵波波速，P(t)是爆炸压力函数，其表达式为P(t)＝P_VN(eγ/n)ⁿtⁿe^‑γt，其中，P_VN是爆炸峰值压力，一般n＝3，γ＝0.7，求解以上控制方程可得：式中，λ₁＝‑C₄S11：根据位移势力函数的表达式，写出径向位移的表达式为：S12：根据径向位移表达式及关于r的偏导代入到各应力分量与径向位移之间的关系式中，可求得在爆炸荷载作用下介质的瞬时应力状态S13：将静态应力状态与动态应力响应叠加，可得到在真三轴地应力作用下的岩石爆炸动力响应S14：将σ_ij投影至应力不变量空间(‑I₁ vs.J₂)中，绘制点P的应力轨迹，根据应力轨迹与屈服面的相互关系，分析岩石的爆炸响应。