[发明专利]一种分析爆破作用下台阶边坡稳定性的方法

权利要求书

1.一种分析爆破作用下台阶边坡稳定性的方法，其特征在于，包括以下步骤：

步骤1、对现场台阶边坡进行数据采集，并对现场的台阶边坡进行取样；

步骤2、在现场台阶边坡的斜坡表面或平台中部安装振动监测设备，捕捉每次爆破过程中传感器安装位置的振动时程曲线；

步骤3、对现场取样的岩块进行实验，获得岩块的性能参数；

步骤4、根据取样岩石的性能参数和所采集的现场台阶边坡数据进行几何建模，并进行单元划分，获得边坡爆破稳定性分析的三维数值模型；

步骤5、对上述三维数值模型中的炸药的本构模型、岩块的本构模型、结构面的本构模型、孔隙渗流及裂隙渗流的本构模型进行模型选取，确定所选择模型的输入参数；

具体如下：

所述的炸药的本构模型选取朗道模型或JWL模型；当选择朗道模型时，输入参数包括：装药密度、爆速和爆热；当选择JWL模型时，输入参数包括：装药密度、爆炸产物的比内能及实验拟合参数；

所述的岩块的本构模型选取Mohr-Coulomb理想弹塑性模型或Drucker-Prager理想弹塑性模型，输入参数包括：岩体密度、弹性模量、泊松比、粘聚力、内摩擦角和抗拉强度；

所述的结构面的本构模型选取Mohr-Coulomb应变软化模型，输入参数包括：结构面法向接触刚度、切向接触刚度、粘聚力、内摩擦角、抗拉强度、拉伸断裂应变及剪切断裂应变；

所述的孔隙渗流及裂隙渗流的本构模型选取为非饱和达西渗流模型，输入参数包括：孔隙率、渗透系数和裂隙初始开度；

步骤6、设置三维数值模型的初边值条件；

步骤7、采用有限元与离散元耦合的数值模拟方法对三维数值模型进行分析计算，获得振动传感器埋设位置的振动时程曲线，通过调整三维数值模型中的弹性参数，使计算所获振动时程曲线与现场实际振动时程曲线的振动波形、振动幅值和纵波波速一致；

步骤8、采用有限元与离散元耦合的数值模拟方法对调整后的三维数值模型进行分析计算，获得当前状态下边坡的破裂面积；

步骤9、不断调整岩块及结构面的强度参数，以数值的收敛性为判别指标，找到爆破载荷作用下台阶边坡发生失稳滑移的临界状态，获得临界破裂面积；

步骤10、将步骤8中所获得的边坡的当前破裂面积除以步骤9中所获得的边坡的临界破裂面积，获得台阶边坡在此次爆破载荷下的破裂度；

步骤11、判断破裂度所属范围，若所获破裂度小于等于0.5，则边坡处于稳定状态；若破裂度大于0.5且小于1，则边坡处于欠稳定状态；若破裂度等于1，则边坡处于稳定与不稳定的临界点；若破裂度大于1，则边坡处于失稳状态；

步骤12、反复执行步骤4至步骤11，获得多次爆破作用下，台阶边坡的破裂度及稳定性的演化情况。

2.根据权利要求1所述的分析爆破作用下台阶边坡稳定性的方法，其特征在于，步骤1所述的数据采集，包括台阶边坡的几何形态、地层特性、结构面分布情况、地下水情况和爆破参数；

所述的几何形态包括：台阶边坡高度、台阶边坡角度、安全平台宽度和台阶数量；

所述的地层特性包括地层数量、地层的岩石种类、各地层的倾向及倾角；

所述的结构面分布情况包括结构面组数、结构面产状、结构面间距和结构面迹长；

所述的地下水情况包括地下水位高度和水液面的空间形态；

所述的爆破参数包括炸药类型、耦合系数、炮孔直径、炮孔深度、堵塞长度、超深、单孔装药量、底盘抵抗线、炮孔数量、间排距、起爆顺序和延时。

3.根据权利要求1所述的分析爆破作用下台阶边坡稳定性的方法，其特征在于，步骤3所述的对现场取样的岩块进行实验，获得岩块的性能参数，其中，所述的实验包括室内声波实验、单轴压缩实验、三轴实验、直剪实验、巴西劈裂实验、锤击波动实验及渗透实验，获得岩块的性能参数包括岩块的弹性模量、泊松比、粘聚力、内摩擦角、抗拉强度、孔隙率及渗透系数，获得岩块中结构面的性能参数包括结构面的法向接触刚度、切向接触刚度、粘聚力、内摩擦角、抗拉强度、拉伸断裂应变及剪切断裂应变。

4.根据权利要求1所述的分析爆破作用下台阶边坡稳定性的方法，其特征在于，步骤4所述的几何建模采用点-线-面-体自下而上的方式进行建模，并采用Delaunay方法进行单元划分，获得边坡爆破稳定性分析的三维数值模型。