；
[0063]步骤10、将步骤8中所获得的边坡的当前破裂面积除以步骤9中所获得的边坡的临界破裂面积，获得台阶边坡在此次爆破载荷下的破裂度D，D = AJAc;
[0064]步骤11、判断破裂度所属范围，若所获破裂度小于等于0.5，则边坡处于稳定状态；若破裂度大于0.5且小于1，则边坡处于欠稳定状态；若破裂度等于1，则边坡处于稳定与不稳定的临界点；若破裂度大于1，则边坡处于失稳状态；
[0065]本发明实施例中，输入反演后的弹性参数及室内实验获得的强度参数，分析此次爆破作用引起的台阶边坡4至6的无量纲破裂面积分别为6 %、3 %及I %。遍历参数获得台阶边坡4至6在临界状态下的无量纲临界破裂面积，分别为42%、65%及54%。计算出台阶边坡4至6的破裂度分别为14%、4.6%及1.9%。从台阶边坡4至6的破裂度可以看出，3个台阶边坡在此次爆破作用下均处于稳定状态；
[0066]步骤12、反复执行步骤4至步骤11，获得多次爆破作用下，台阶边坡的破裂度及稳定性的演化情况。
[0067]本方法集成了现场调查、室内实验、振动监测、数值分析及破裂度评价等多种手段，是一种边坡爆破稳定性的综合分析法。该方法通过现场调查及室内实验获取爆区台阶边坡的基本几何形态、结构面空间展布特征、爆破参数、岩体及炸药的物理力学特性等，通过振动监测捕捉历次爆破作业时边坡典型位置的振动时程曲线，通过数值模拟与振动监测的结合反演修正数值计算参数，并利用数值模拟获得当前爆破参数下的破裂度，最后利用破裂度指标评价边坡的稳定性状态。
【主权项】
1.一种分析爆破作用下台阶边坡稳定性的方法，其特征在于，包括以下步骤: 步骤1、对现场台阶边坡进行数据采集，并对现场的台阶边坡进行取样； 步骤2、在现场台阶边坡的斜坡表面或平台中部安装振动监测设备，捕捉每次爆破过程中传感器安装位置的振动时程曲线； 步骤3、对现场取样的岩块进行实验，获得岩块的性能参数； 步骤4、根据取样岩石的性能参数和所采集的现场台阶边坡数据进行几何建模，并进行单元划分，获得边坡爆破稳定性分析的三维数值模型； 步骤5、对上述三维数值模型中的炸药的本构模型、岩块的本构模型、结构面的本构模型、孔隙渗流及裂隙渗流的本构模型进行模型选取，确定所选择模型的输入参数； 步骤6、设置三维数值模型的初边值条件； 步骤7、采用有限元与离散元耦合的数值模拟方法对三维数值模型进行分析计算，获得振动传感器埋设位置的振动时程曲线，通过调整三维数值模型中的弹性参数，使计算所获振动时程曲线与现场实际振动时程曲线的振动波形、振动幅值和纵波波速一致； 步骤8、采用有限元与离散元耦合的数值模拟方法对调整后的三维数值模型进行分析计算，获得当前状态下边坡的破裂面积； 步骤9、不断调整岩块及结构面的强度参数，以数值的收敛性为判别指标，找到爆破载荷作用下台阶边坡发生失稳滑移的临界状态，获得临界破裂面积； 步骤10、将步骤8中所获得的边坡的当前破裂面积除以步骤9中所获得的边坡的临界破裂面积，获得台阶边坡在此次爆破载荷下的破裂度； 步骤11、判断破裂度所属范围，若所获破裂度小于等于0.5，则边坡处于稳定状态；若破裂度大于0.5且小于I，则边坡处于欠稳定状态；若破裂度等于I，则边坡处于稳定与不稳定的临界点；若破裂度大于1，则边坡处于失稳状态； 步骤12、反复执行步骤4至步骤11，获得多次爆破作用下，台阶边坡的破裂度及稳定性的演化情况。2.根据权利要求1所述的分析爆破作用下台阶边坡稳定性的方法，其特征在于，步骤I所述的数据采集，包括台阶边坡的几何形态、地层特性、结构面分布情况、地下水情况和爆破参数； 所述的几何形态包括:台阶边坡高度、台阶边坡角度、安全平台宽度和台阶数量； 所述的地层特性包括地层数量、地层的岩石种类、各地层的倾向及倾角； 所述的结构面分布情况包括结构面组数、结构面产状、结构面间距和结构面迹长； 所述的地下水情况包括地下水位高度和水液面的空间形态； 所述的爆破参数包括炸药类型、耦合系数、炮孔直径、炮孔深度、堵塞长度、超深、单孔装药量、底盘抵抗线、炮孔数量、间排距、起爆顺序和延时。3.根据权利要求1所述的分析爆破作用下台阶边坡稳定性的方法，其特征在于，步骤3所述的对现场取样的岩块进行实验，获得岩块的性能参数，其中，所述的实验包括室内声波实验、单轴压缩实验、三轴实验、直剪实验、巴西劈裂实验、锤击波动实验及渗透实验，获得岩块的性能参数包括岩块的弹性模量、泊松比、粘聚力、内摩擦角、抗拉强度、孔隙率及渗透系数，获得岩块中结构面的性能参数包括结构面的法向接触刚度、切向接触刚度、粘聚力、内摩擦角、抗拉强度、拉伸断裂应变及剪切断裂应变。4.根据权利要求1所述的分析爆破作用下台阶边坡稳定性的方法，其特征在于，步骤4所述的几何建模采用点-线-面-体自下而上的方式进行建模，并采用Delaunay方法进行单元划分，获得边坡爆破稳定性分析的三维数值模型。5.根据权利要求1所述的分析爆破作用下台阶边坡稳定性的方法，其特征在于，步骤5所述的对上述三维数值模型中的炸药的本构模型、岩块的本构模型、结构面的本构模型、孔隙渗流及裂隙渗流的本构模型进行模型选取，确定所选择模型的输入参数，具体如下: 所述的炸药的本构模型选取朗道模型或JWL模型；当选择朗道模型时，输入参数包括:装药密度、爆速和爆热；当选择JWL模型时，输入参数包括:装药密度、爆炸产物的比内能及实验拟合参数； 所述的岩块的本构模型选取Mohr-Coulomb理想弹塑性模型或Drucker-Prager理想弹塑性模型，输入参数包括:岩体密度、弹性模量、泊松比、粘聚力、内摩擦角和抗拉强度； 所述的结构面的本构模型选取Mohr-Coulomb应变软化模型，输入参数包括:结构面法向接触刚度、切向接触刚度、粘聚力、内摩擦角、抗拉强度、拉伸断裂应变及剪切断裂应变； 所述的孔隙渗流及裂隙渗流的本构模型选取为非饱和达西渗流模型，输入参数包括:孔隙率、渗透系数和裂隙初始开度。6.根据权利要求1所述的分析爆破作用下台阶边坡稳定性的方法，其特征在于，步骤6所述的设置三维数值模型的初边值条件，包括:静力计算时在模型底部及四周设置法向约束条件，在整个模型上施加重力条件；爆炸计算时在模型底部及四周设置无反射边界条件。7.根据权利要求1所述的分析爆破作用下台阶边坡稳定性的方法，其特征在于，步骤7所述的通过调整三维数值模型中的弹性参数，所调整的弹性参数包括:岩体的弹性模量和结构面的接触刚度。8.根据权利要求1所述的分析爆破作用下台阶边坡稳定性的方法，其特征在于，步骤9所述的不断调整岩块及结构面的强度参数，包括粘聚力、内摩擦角、抗拉强度、拉伸断裂应变和剪切断裂应变；所述的以数值的收敛性为判别指标，即当数值发散时，即为临界状态。
【专利摘要】本发明一种分析爆破作用下台阶边坡稳定性的方法，属于露天矿开采技术领域，本发明与传统的方法相比，此方法将数值模拟与振动监测相结合，利用坡体的振动信息修正坡体的材料参数，可以为数值模拟提供更为准确的输入参数；采用有限元与离散元耦合的数值模拟方法进行力学分析，可以完整刻画边坡在爆破载荷作用下的渐进破坏过程，计算结果更为准确；采用破裂度进行边坡稳定性的评价，可直接反映坡体内部的损伤破裂情况，评价指标更为合理。
【IPC分类】G06F17/50
【公开号】CN105224742
【申请号】CN201510632042
【发明人】潘鹏飞, 孙厚广, 冯春, 韩忠和, 范晓明, 乔继延, 郭汝坤, 李世海
【申请人】鞍钢集团矿业公司
【公开日】2016年1月6日
【申请日】2015年9月29日