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

摘要

本发明一种分析爆破作用下台阶边坡稳定性的方法，属于露天矿开采技术领域，本发明与传统的方法相比，此方法将数值模拟与振动监测相结合，利用坡体的振动信息修正坡体的材料参数，可以为数值模拟提供更为准确的输入参数；采用有限元与离散元耦合的数值模拟方法进行力学分析，可以完整刻画边坡在爆破载荷作用下的渐进破坏过程，计算结果更为准确；采用破裂度进行边坡稳定性的评价，可直接反映坡体内部的损伤破裂情况，评价指标更为合理。

主权项

一种分析爆破作用下台阶边坡稳定性的方法，其特征在于，包括以下步骤：步骤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，获得多次爆破作用下，台阶边坡的破裂度及稳定性的演化情况。