CN105224742A - 一种分析***作用下台阶边坡稳定性的方法 - Google Patents
一种分析***作用下台阶边坡稳定性的方法 Download PDFInfo
- Publication number
- CN105224742A CN105224742A CN201510632042.1A CN201510632042A CN105224742A CN 105224742 A CN105224742 A CN 105224742A CN 201510632042 A CN201510632042 A CN 201510632042A CN 105224742 A CN105224742 A CN 105224742A
- Authority
- CN
- China
- Prior art keywords
- model
- slope
- sillar
- parameter
- structural plane
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明一种分析***作用下台阶边坡稳定性的方法,属于露天矿开采技术领域,本发明与传统的方法相比,此方法将数值模拟与振动监测相结合,利用坡体的振动信息修正坡体的材料参数,可以为数值模拟提供更为准确的输入参数;采用有限元与离散元耦合的数值模拟方法进行力学分析,可以完整刻画边坡在***载荷作用下的渐进破坏过程,计算结果更为准确;采用破裂度进行边坡稳定性的评价,可直接反映坡体内部的损伤破裂情况,评价指标更为合理。
Description
技术领域
本发明属于露天矿开采技术领域,具体涉及一种分析***作用下台阶边坡稳定性的方法。
背景技术
露天矿的台阶边坡在***振动的反复作用下将出现材料性能的劣化,进而引起边坡的失稳破坏。目前,分析边坡稳定性的方法主要包括极限平衡法及数值模拟法等两大类。
极限平衡法基于刚塑性假设及一致性运动的假设,采用单一自由度来表征滑坡体的整体稳定性状态,因此无法刻画边坡的非一致性运动对其稳定性的影响。此外,极限平衡法是一种静态稳定性分析方法,动态问题需要采用拟静力法等效后方能进行分析,因此无法准确反映***振动对边坡材料性能的劣化效应。
数值模拟法包括基于连续介质的方法及基于非连续介质的方法等两类,与极限平衡法相比,数值模拟法可以实现***对边坡稳定性影响的全过程动态分析。基于连续介质的数值模拟方法包括有限元法、有限体积法、有限差分法等,此类方法通过在微元上引入损伤或塑性本构,实现了边坡损伤演化及塑性滑移过程的精确模拟,但却无法准确刻画边坡失稳过程中裂纹的扩展贯通过程及滑移过程中的碰撞解体过程。基于非连续介质的数值模拟方法包括块体离散元法、颗粒离散元法及DDA等,此类方法在模拟离散块体间的接触碰撞方面具有明显优势,但却无法准确刻画材料的连续变形及塑性、损伤状态。
目前,无论极限平衡法或数值模拟法,均采用安全系数作为边坡稳定性的评价指标。在极限平衡法中,安全系数的实质是抗滑力与下滑力的比值;在数值模拟法中,安全系数的实质是当前材料参数与弱化后材料参数的比值。上述两种方法中的安全系数均是一种***宏观稳定性的等效描述,无法准确反映坡体内部的损伤破裂程度。
由于地质结构的复杂性及演化性,***作用下地质体的材料参数将会不断演化。然而,无论极限平衡法或数值模拟法,均采用固定的材料参数进行分析计算,没有借助监测结果进行材料参数的动态反演分析及修正。
发明内容
针对现有技术的不足,本发明提出一种分析***作用下台阶边坡稳定性的方法,以达到将监测与数值模拟相结合,准确分析露天矿***开采对周边台阶边坡稳定性影响程度的目的。
一种分析***作用下台阶边坡稳定性的方法,包括以下步骤:
步骤1、对现场台阶边坡进行数据采集,并对现场的台阶边坡进行取样;
步骤2、在现场台阶边坡的斜坡表面或平台中部安装振动监测设备,捕捉每次***过程中传感器安装位置的振动时程曲线;
步骤3、对现场取样的岩块进行实验,获得岩块的性能参数;
步骤4、根据取样岩石的性能参数和所采集的现场台阶边坡数据进行几何建模,并进行单元划分,获得边坡***稳定性分析的三维数值模型;
步骤5、对上述三维数值模型中的***的本构模型、岩块的本构模型、结构面的本构模型、孔隙渗流及裂隙渗流的本构模型进行模型选取,确定所选择模型的输入参数;
步骤6、设置三维数值模型的初边值条件;
步骤7、采用有限元与离散元耦合的数值模拟方法对三维数值模型进行分析计算,获得振动传感器埋设位置的振动时程曲线,通过调整三维数值模型中的弹性参数,使计算所获振动时程曲线与现场实际振动时程曲线的振动波形、振动幅值和纵波波速一致;
步骤8、采用有限元与离散元耦合的数值模拟方法对调整后的三维数值模型进行分析计算,获得当前状态下边坡的破裂面积;
步骤9、不断调整岩块及结构面的强度参数,以数值的收敛性为判别指标,找到***载荷作用下台阶边坡发生失稳滑移的临界状态,获得临界破裂面积;
步骤10、将步骤8中所获得的边坡的当前破裂面积除以步骤9中所获得的边坡的临界破裂面积,获得台阶边坡在此次***载荷下的破裂度;
步骤11、判断破裂度所属范围,若所获破裂度小于等于0.5,则边坡处于稳定状态;若破裂度大于0.5且小于1,则边坡处于欠稳定状态;若破裂度等于1,则边坡处于稳定与不稳定的临界点;若破裂度大于1,则边坡处于失稳状态;
步骤12、反复执行步骤4至步骤11,获得多次***作用下,台阶边坡的破裂度及稳定性的演化情况。
步骤1所述的数据采集,包括台阶边坡的几何形态、地层特性、结构面分布情况、地下水情况和***参数;
所述的几何形态包括:台阶边坡高度、台阶边坡角度、安全平台宽度和台阶数量;
所述的地层特性包括地层数量、地层的岩石种类、各地层的倾向及倾角;
所述的结构面分布情况包括结构面组数、结构面产状、结构面间距和结构面迹长;
所述的地下水情况包括地下水位高度和水液面的空间形态;
所述的***参数包括***类型、耦合系数、炮孔直径、炮孔深度、堵塞长度、超深、单孔装药量、底盘抵抗线、炮孔数量、间排距、起爆顺序和延时。
步骤3所述的对现场取样的岩块进行实验,获得岩块的性能参数,其中,所述的实验包括室内声波实验、单轴压缩实验、三轴实验、直剪实验、巴西劈裂实验、锤击波动实验及渗透实验,获得岩块的性能参数包括岩块的弹性模量、泊松比、粘聚力、内摩擦角、抗拉强度、孔隙率及渗透系数,获得岩块中结构面的性能参数包括结构面的法向接触刚度、切向接触刚度、粘聚力、内摩擦角、抗拉强度、拉伸断裂应变及剪切断裂应变。
步骤4所述的几何建模采用点-线-面-体自下而上的方式进行建模,并采用Delaunay方法进行单元划分,获得边坡***稳定性分析的三维数值模型。
步骤5所述的对上述三维数值模型中的***的本构模型、岩块的本构模型、结构面的本构模型、孔隙渗流及裂隙渗流的本构模型进行模型选取,确定所选择模型的输入参数,具体如下:
所述的***的本构模型选取朗道模型或JWL模型;当选择朗道模型时,输入参数包括:装药密度、爆速和爆热;当选择JWL模型时,输入参数包括:装药密度、***产物的比内能及实验拟合参数;
所述的岩块的本构模型选取Mohr-Coulomb理想弹塑性模型或Drucker-Prager理想弹塑性模型,输入参数包括:岩体密度、弹性模量、泊松比、粘聚力、内摩擦角和抗拉强度;
所述的结构面的本构模型选取Mohr-Coulomb应***化模型,输入参数包括:结构面法向接触刚度、切向接触刚度、粘聚力、内摩擦角、抗拉强度、拉伸断裂应变及剪切断裂应变;
所述的孔隙渗流及裂隙渗流的本构模型选取为非饱和达西渗流模型,输入参数包括:孔隙率、渗透系数和裂隙初始开度。
步骤6所述的设置三维数值模型的初边值条件,包括:静力计算时在模型底部及四周设置法向约束条件,在整个模型上施加重力条件;***计算时在模型底部及四周设置无反射边界条件。
步骤7所述的通过调整三维数值模型中的弹性参数,所调整的弹性参数包括:岩体的弹性模量和结构面的接触刚度。
步骤9所述的不断调整岩块及结构面的强度参数,包括粘聚力、内摩擦角、抗拉强度、拉伸断裂应变和剪切断裂应变;所述的以数值的收敛性为判别指标,即当数值发散时,即为临界状态。
本发明优点:
本发明提出一种分析***作用下台阶边坡稳定性的方法,与传统的方法相比,此方法将数值模拟与振动监测相结合,利用坡体的振动信息修正坡体的材料参数,可以为数值模拟提供更为准确的输入参数;采用有限元与离散元耦合的数值模拟方法进行力学分析,可以完整刻画边坡在***载荷作用下的渐进破坏过程,计算结果更为准确;采用破裂度进行边坡稳定性的评价,可直接反映坡体内部的损伤破裂情况,评价指标更为合理。
附图说明
图1为本发明一种实施例的分析***作用下台阶边坡稳定性的方法流程图;
图2为本发明一种实施例的***作用下台阶边坡稳定性分析示意图,其中,1表示炮孔,2表示***,3表示堵塞,4表示第一台阶边坡,5表示第二台阶边坡,6表示第三台阶边坡,7表示顺倾节理,8表示第一速度传感器,9表示第二速度传感器,10表示第三速度传感器,11表示第四速度传感器,12表示地下水位。
具体实施方式
下面结合附图对本发明一种实施例做进一步说明。
本发明实施例中,对某露天铁矿的台阶边坡进行稳定性分析,具体方法步骤如图1所示,包括以下步骤:
步骤1、对现场台阶边坡进行数据采集,并对现场的台阶边坡的平台处进行取样;
所述的数据采集,包括台阶边坡的几何形态、地层特性、结构面分布情况、地下水情况和***参数;所述的几何形态包括:台阶边坡高度、台阶边坡角度、安全平台宽度和台阶数量;所述的地层特性包括地层数量、地层的岩石种类、各地层的倾向及倾角;所述的结构面分布情况包括结构面组数、结构面产状、结构面间距和结构面迹长;所述的地下水情况包括地下水位高度和水液面的空间形态;所述的***参数包括***类型、耦合系数、炮孔直径、炮孔深度、堵塞长度、超深、单孔装药量、底盘抵抗线、炮孔数量、间排距、起爆顺序和延时;
本发明实施例中,对某露天铁矿的台阶边坡进行稳定性分析,如图2所示,台阶边坡4至6的高度均为12m、坡面角65°、平台宽度8m;台阶边坡4至6被一组顺倾节理7切割,节理倾角45°,间距3m;炮孔1的直径为25cm、孔深为15m;炮孔1共3排,每排包含7个炮孔;炮孔1的底盘抵抗线为5m,间排距均为6.5m,采用等边三角形方式布孔,最后一排炮孔1到台阶边坡4坡脚的距离为12m;地下水位12位于***作业平台以下10m处;***2为乳化***,装药密度为1150kg/m3,爆轰速度为4250m/s,爆热为3.4MJ/kg;堵塞3的长度为7m,起爆顺序为逐孔起爆,孔间延时25ms,排间延时42ms;本发明实施例中,在台阶边坡4的底部平台进行取样;
步骤2、在现场台阶边坡的斜坡表面或平台中部安装振动监测设备,捕捉每次***过程中传感器安装位置的振动时程曲线;
本发明实施例中,所述的振动监测设备采用振动速度传感设备,在台阶边坡4至6对应的平台中部分别安装4支速度传感器8至11,用于测试***引起的边坡振动情况;基于时程曲线拾取各传感器的起振时间,计算出边坡不同位置的纵波波速;
步骤3、对现场取样的岩块进行实验,获得岩块的性能参数;
本发明实施例中,所述的对现场取样的岩块进行实验,获得岩块的性能参数,其中,所述的实验包括室内声波实验、单轴压缩实验、三轴实验、直剪实验、巴西劈裂实验、锤击波动实验及渗透实验,获得岩块的性能参数包括岩块的弹性模量、泊松比、粘聚力、内摩擦角、抗拉强度、孔隙率及渗透系数,获得岩块中结构面的性能参数包括结构面的法向接触刚度、切向接触刚度、粘聚力、内摩擦角、抗拉强度、拉伸断裂应变及剪切断裂应变;
本发明实施例中,测试结果表明,台阶边坡4至6处铁矿的密度为3200kg/m3,弹性模量60GPa,泊松比0.25、粘聚力36MPa、抗拉强度12MPa、内摩擦角40°;节理7的法向接触刚度为50GPa/m,切向接触刚度为35GPa/m,粘聚力为8MPa,抗拉强度为3MPa,内摩擦角为30°,拉伸断裂应变为0.1%,剪切断裂应变为0.3%;
步骤4、根据取样岩石的性能参数和所采集的现场台阶边坡数据进行几何建模,并进行单元划分,获得边坡***稳定性分析的三维数值模型;
本发明实施例中,所述的几何建模采用点-线-面-体自下而上的方式进行建模(采用几何建模及网格划分软件,例如:GiD软件、Patran软件、Gmsh软件),并采用Delaunay方法进行单元划分,获得边坡***稳定性分析的三维数值模型;
步骤5、对上述三维数值模型中的***的本构模型、岩块的本构模型、结构面的本构模型、孔隙渗流及裂隙渗流的本构模型进行模型选取,确定所选择模型的输入参数;
本发明实施例中,所述的对上述三维数值模型中的***的本构模型、岩块的本构模型、结构面的本构模型、孔隙渗流及裂隙渗流的本构模型进行模型选取,确定所选择模型的输入参数,具体如下:
所述的***的本构模型选取朗道模型或JWL模型;当选择朗道模型时,输入参数包括:装药密度、爆速和爆热;当选择JWL模型时,输入参数包括:装药密度、***产物的比内能及实验拟合参数,所述实验拟合参数包括压力A、压力B、无量纲数R1、无量纲数R2及无量纲数w;所述的岩块的本构模型选取Mohr-Coulomb理想弹塑性模型或Drucker-Prager理想弹塑性模型,输入参数包括:岩体密度、弹性模量、泊松比、粘聚力、内摩擦角和抗拉强度;所述的结构面的本构模型选取Mohr-Coulomb应***化模型,输入参数包括:结构面法向接触刚度、切向接触刚度、粘聚力、内摩擦角、抗拉强度、拉伸断裂应变及剪切断裂应变;所述的孔隙渗流及裂隙渗流的本构模型选取为非饱和达西渗流模型,输入参数包括:孔隙率、渗透系数和裂隙开度;
本发明实施例中,利用通用前后处理软件GiD建立三维数值模型,共剖分四面体网格51.2万。采用有限元与离散元耦合的数值方法进行力学分析,爆源模型选择为朗道模型,岩块模型选择为Mohr-Coulomb理想弹塑性模型,结构面模型选择为Mohr-Coulomb应***化模型。
步骤6、设置三维数值模型的初边值条件,包括:静力计算时在模型底部及四周设置法向约束条件,在整个模型上施加重力条件;***计算时在模型底部及四周设置无反射边界条件。
步骤7、采用有限元与离散元耦合的数值模拟方法对三维数值模型进行分析计算,获得振动传感器埋设位置的振动时程曲线,通过调整三维数值模型中的弹性参数,包括:岩体的弹性模量和结构面的接触刚度,使计算所获振动时程曲线与现场实际振动时程曲线的振动波形、振动幅值和纵波波速一致;
本发明实施例中,采用有限元与离散元耦合的数值模拟方法及软件(如CDEM方法、FEM/DEM方法、Rockfield软件、UDEC/3DEC软件);
本发明实施例中,借助有限元与离散元耦合的数值方法及振动速度传感器8至11的振动时程曲线,对铁矿及节理面的弹性参数进行反分析,反分析获得的铁矿弹性模量为52GPa、泊松比为0.30,节理7的法向接触刚度为60GPa/m、切向接触刚度为40GPa/m。
步骤8、采用有限元与离散元耦合的数值模拟方法对调整后的三维数值模型进行分析计算,获得当前状态下边坡的破裂面积An;
步骤9、不断调整岩块及结构面的强度参数,包括粘聚力、内摩擦角、抗拉强度、拉伸断裂应变和剪切断裂应变,以数值的收敛性为判别指标,即当数值发散时,找到***载荷作用下台阶边坡发生失稳滑移的临界状态,获得临界破裂面积Ac;
步骤10、将步骤8中所获得的边坡的当前破裂面积除以步骤9中所获得的边坡的临界破裂面积,获得台阶边坡在此次***载荷下的破裂度D,D=An/Ac;
步骤11、判断破裂度所属范围,若所获破裂度小于等于0.5,则边坡处于稳定状态;若破裂度大于0.5且小于1,则边坡处于欠稳定状态;若破裂度等于1,则边坡处于稳定与不稳定的临界点;若破裂度大于1,则边坡处于失稳状态;
本发明实施例中,输入反演后的弹性参数及室内实验获得的强度参数,分析此次***作用引起的台阶边坡4至6的无量纲破裂面积分别为6%、3%及1%。遍历参数获得台阶边坡4至6在临界状态下的无量纲临界破裂面积,分别为42%、65%及54%。计算出台阶边坡4至6的破裂度分别为14%、4.6%及1.9%。从台阶边坡4至6的破裂度可以看出,3个台阶边坡在此次***作用下均处于稳定状态;
步骤12、反复执行步骤4至步骤11,获得多次***作用下,台阶边坡的破裂度及稳定性的演化情况。
本方法集成了现场调查、室内实验、振动监测、数值分析及破裂度评价等多种手段,是一种边坡***稳定性的综合分析法。该方法通过现场调查及室内实验获取爆区台阶边坡的基本几何形态、结构面空间展布特征、***参数、岩体及***的物理力学特性等,通过振动监测捕捉历次***作业时边坡典型位置的振动时程曲线,通过数值模拟与振动监测的结合反演修正数值计算参数,并利用数值模拟获得当前***参数下的破裂度,最后利用破裂度指标评价边坡的稳定性状态。
Claims (8)
1.一种分析***作用下台阶边坡稳定性的方法,其特征在于,包括以下步骤:
步骤1、对现场台阶边坡进行数据采集,并对现场的台阶边坡进行取样;
步骤2、在现场台阶边坡的斜坡表面或平台中部安装振动监测设备,捕捉每次***过程中传感器安装位置的振动时程曲线;
步骤3、对现场取样的岩块进行实验,获得岩块的性能参数;
步骤4、根据取样岩石的性能参数和所采集的现场台阶边坡数据进行几何建模,并进行单元划分,获得边坡***稳定性分析的三维数值模型;
步骤5、对上述三维数值模型中的***的本构模型、岩块的本构模型、结构面的本构模型、孔隙渗流及裂隙渗流的本构模型进行模型选取,确定所选择模型的输入参数;
步骤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方法进行单元划分,获得边坡***稳定性分析的三维数值模型。
5.根据权利要求1所述的分析***作用下台阶边坡稳定性的方法,其特征在于,步骤5所述的对上述三维数值模型中的***的本构模型、岩块的本构模型、结构面的本构模型、孔隙渗流及裂隙渗流的本构模型进行模型选取,确定所选择模型的输入参数,具体如下:
所述的***的本构模型选取朗道模型或JWL模型;当选择朗道模型时,输入参数包括:装药密度、爆速和爆热;当选择JWL模型时,输入参数包括:装药密度、***产物的比内能及实验拟合参数;
所述的岩块的本构模型选取Mohr-Coulomb理想弹塑性模型或Drucker-Prager理想弹塑性模型,输入参数包括:岩体密度、弹性模量、泊松比、粘聚力、内摩擦角和抗拉强度;
所述的结构面的本构模型选取Mohr-Coulomb应***化模型,输入参数包括:结构面法向接触刚度、切向接触刚度、粘聚力、内摩擦角、抗拉强度、拉伸断裂应变及剪切断裂应变;
所述的孔隙渗流及裂隙渗流的本构模型选取为非饱和达西渗流模型,输入参数包括:孔隙率、渗透系数和裂隙初始开度。
6.根据权利要求1所述的分析***作用下台阶边坡稳定性的方法,其特征在于,步骤6所述的设置三维数值模型的初边值条件,包括:静力计算时在模型底部及四周设置法向约束条件,在整个模型上施加重力条件;***计算时在模型底部及四周设置无反射边界条件。
7.根据权利要求1所述的分析***作用下台阶边坡稳定性的方法,其特征在于,步骤7所述的通过调整三维数值模型中的弹性参数,所调整的弹性参数包括:岩体的弹性模量和结构面的接触刚度。
8.根据权利要求1所述的分析***作用下台阶边坡稳定性的方法,其特征在于,步骤9所述的不断调整岩块及结构面的强度参数,包括粘聚力、内摩擦角、抗拉强度、拉伸断裂应变和剪切断裂应变;所述的以数值的收敛性为判别指标,即当数值发散时,即为临界状态。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510632042.1A CN105224742B (zh) | 2015-09-29 | 2015-09-29 | 一种分析***作用下台阶边坡稳定性的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510632042.1A CN105224742B (zh) | 2015-09-29 | 2015-09-29 | 一种分析***作用下台阶边坡稳定性的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105224742A true CN105224742A (zh) | 2016-01-06 |
CN105224742B CN105224742B (zh) | 2018-01-23 |
Family
ID=54993708
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510632042.1A Active CN105224742B (zh) | 2015-09-29 | 2015-09-29 | 一种分析***作用下台阶边坡稳定性的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105224742B (zh) |
Cited By (33)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105758509A (zh) * | 2016-04-25 | 2016-07-13 | 宏大矿业有限公司 | 一种露天矿山岩体声速的现场测量方法 |
CN105866464A (zh) * | 2016-04-21 | 2016-08-17 | 宏大矿业有限公司 | 一种现场快速测量深孔***时***爆速的***及方法 |
CN105928601A (zh) * | 2016-04-25 | 2016-09-07 | 宏大矿业有限公司 | 一种测量露天矿山爆区下方台阶岩体声速的方法 |
CN105954357A (zh) * | 2016-04-25 | 2016-09-21 | 宏大矿业有限公司 | 一种同时测量多台阶水平岩体声速的方法 |
CN105973456A (zh) * | 2016-04-25 | 2016-09-28 | 宏大矿业有限公司 | 一种现场测量深孔******与岩石匹配关系的方法 |
CN106339547A (zh) * | 2016-08-26 | 2017-01-18 | 中国水利水电科学研究院 | 一种***的数值模拟方法 |
CN106372314A (zh) * | 2016-08-31 | 2017-02-01 | 广东省交通规划设计研究院股份有限公司 | 一种线性工程岩质边坡稳定性半定量自动评价*** |
CN106596294A (zh) * | 2016-12-30 | 2017-04-26 | 中国科学院地质与地球物理研究所 | 一种模拟顺向边坡动力反应的试验装置及方法 |
CN107169236A (zh) * | 2017-06-15 | 2017-09-15 | 中国水利水电科学研究院 | 一种基于有限元与离散元耦合的虚拟三轴试验仿真方法 |
CN107328919A (zh) * | 2017-08-15 | 2017-11-07 | 绍兴文理学院 | 矿山边坡岩体工程稳定性分级分析方法 |
CN107328920A (zh) * | 2017-08-15 | 2017-11-07 | 绍兴文理学院 | 矿山边坡岩体工程稳定性精准评价方法 |
CN107655691A (zh) * | 2017-08-28 | 2018-02-02 | 西北工业大学 | 一种飞行器进气道喘振锤击波试验装置及方法 |
CN108254262A (zh) * | 2016-12-28 | 2018-07-06 | 中国石油天然气股份有限公司 | 岩石层理裂缝剪切参数预测方法及装置 |
CN108287112A (zh) * | 2018-01-31 | 2018-07-17 | 成都理工大学 | 一种基于三轴压缩试验测定岩石损伤参数的方法 |
CN108331042A (zh) * | 2018-03-05 | 2018-07-27 | 武汉科技大学 | 一种模拟***开挖对边坡影响的试验装置及试验方法 |
CN108920754A (zh) * | 2018-05-25 | 2018-11-30 | 重庆交通大学 | 基于应***化剪切本构模型和位移变化的结构面控制边坡稳定性动态评价方法 |
CN109142671A (zh) * | 2018-07-17 | 2019-01-04 | 中国地质大学(武汉) | 研究***振动下饱水软弱结构面强度渐变劣化规律的方法 |
CN109270243A (zh) * | 2018-07-26 | 2019-01-25 | 重庆交通大学 | 冻土***模型试验***及其试验方法 |
CN109613607A (zh) * | 2018-12-04 | 2019-04-12 | 保利新联***工程集团有限公司 | 一种不同高差***地震波传播测试方法 |
CN109992841A (zh) * | 2019-03-11 | 2019-07-09 | 长江水利委员会长江科学院 | 一种***块度空间全方位分区耦合数值仿真方法 |
CN110360907A (zh) * | 2019-07-31 | 2019-10-22 | 中冶集团资源开发有限公司 | 一种遍布节理破碎岩体的边坡台阶扩帮***方法 |
CN110502825A (zh) * | 2019-08-19 | 2019-11-26 | 青岛理工大学 | 一种提取三维破裂面的方法 |
CN111024588A (zh) * | 2019-12-31 | 2020-04-17 | 山东大学 | 一种体现渗流对岩土体强度弱化的dem接触模型构建方法 |
CN111199108A (zh) * | 2020-01-07 | 2020-05-26 | 包钢集团矿山研究院(有限责任公司) | 一种网格化切分数值模拟爆堆模型的***块度预测方法 |
CN111324942A (zh) * | 2019-12-27 | 2020-06-23 | 昆明理工大学 | 一种考虑滑面动力渐进破坏的地震边坡稳定性分析方法 |
CN112182731A (zh) * | 2020-10-27 | 2021-01-05 | 上海工程技术大学 | 非极限状态二维边坡稳定性评价方法 |
CN112213767A (zh) * | 2020-09-11 | 2021-01-12 | 中国矿业大学 | 一种用于顶板超前预裂***效果的评价方法 |
CN113324831A (zh) * | 2021-06-16 | 2021-08-31 | 山东黄金矿业科技有限公司深井开采实验室分公司 | 一种露天矿边坡内采空区动力失稳破坏机理试验方法 |
CN113357987A (zh) * | 2021-04-15 | 2021-09-07 | 西安科技大学 | 一种自动***设计方法 |
CN113654419A (zh) * | 2020-05-12 | 2021-11-16 | 中国水利水电第三工程局有限公司 | 面板堆石坝过渡料***参数与最优级配曲线的分析方法 |
CN114062127A (zh) * | 2021-10-29 | 2022-02-18 | 煤炭科学研究总院 | 岩体稳定性的检测方法及其装置 |
CN115079248A (zh) * | 2022-05-26 | 2022-09-20 | 中国科学院武汉岩土力学研究所 | 基于裂纹密度及声波阈值反演***损伤范围的方法及设备 |
CN116796453A (zh) * | 2023-04-18 | 2023-09-22 | 江汉大学 | 一种减振孔设计方法、装置、设备及可读存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130206027A1 (en) * | 2010-02-03 | 2013-08-15 | Jan-Åke Bengtsson | Rock cracker cartridge and ignition capsule |
CN104406983A (zh) * | 2014-12-09 | 2015-03-11 | 鞍钢集团矿业公司 | 一种矿石***块度分布的***评价方法 |
CN104850695A (zh) * | 2015-05-14 | 2015-08-19 | 辽宁工程技术大学 | 一种确定有效采空区放顶***方案的方法 |
-
2015
- 2015-09-29 CN CN201510632042.1A patent/CN105224742B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130206027A1 (en) * | 2010-02-03 | 2013-08-15 | Jan-Åke Bengtsson | Rock cracker cartridge and ignition capsule |
CN104406983A (zh) * | 2014-12-09 | 2015-03-11 | 鞍钢集团矿业公司 | 一种矿石***块度分布的***评价方法 |
CN104850695A (zh) * | 2015-05-14 | 2015-08-19 | 辽宁工程技术大学 | 一种确定有效采空区放顶***方案的方法 |
Non-Patent Citations (2)
Title |
---|
付波等: "岩石高边坡***振动局部放大效应分析", 《***》 * |
陈占军等: "***荷载作用下岩石边坡动态响应的FLAC模拟研究", 《***》 * |
Cited By (52)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105866464A (zh) * | 2016-04-21 | 2016-08-17 | 宏大矿业有限公司 | 一种现场快速测量深孔***时***爆速的***及方法 |
CN105866464B (zh) * | 2016-04-21 | 2019-07-30 | 宏大矿业有限公司 | 一种现场快速测量深孔***时***爆速的***及方法 |
CN105758509A (zh) * | 2016-04-25 | 2016-07-13 | 宏大矿业有限公司 | 一种露天矿山岩体声速的现场测量方法 |
CN105928601A (zh) * | 2016-04-25 | 2016-09-07 | 宏大矿业有限公司 | 一种测量露天矿山爆区下方台阶岩体声速的方法 |
CN105954357A (zh) * | 2016-04-25 | 2016-09-21 | 宏大矿业有限公司 | 一种同时测量多台阶水平岩体声速的方法 |
CN105973456A (zh) * | 2016-04-25 | 2016-09-28 | 宏大矿业有限公司 | 一种现场测量深孔******与岩石匹配关系的方法 |
CN105973456B (zh) * | 2016-04-25 | 2019-07-30 | 宏大矿业有限公司 | 一种现场测量深孔******与岩石匹配关系的方法 |
CN105758509B (zh) * | 2016-04-25 | 2019-05-10 | 宏大***有限公司 | 一种露天矿山岩体声速的现场测量方法 |
CN105928601B (zh) * | 2016-04-25 | 2019-07-30 | 宏大矿业有限公司 | 一种测量露天矿山爆区下方台阶岩体声速的方法 |
CN105954357B (zh) * | 2016-04-25 | 2019-10-11 | 宏大矿业有限公司 | 一种同时测量多台阶水平岩体声速的方法 |
CN106339547A (zh) * | 2016-08-26 | 2017-01-18 | 中国水利水电科学研究院 | 一种***的数值模拟方法 |
CN106339547B (zh) * | 2016-08-26 | 2019-04-23 | 中国水利水电科学研究院 | 一种***的数值模拟方法 |
CN106372314A (zh) * | 2016-08-31 | 2017-02-01 | 广东省交通规划设计研究院股份有限公司 | 一种线性工程岩质边坡稳定性半定量自动评价*** |
CN108254262A (zh) * | 2016-12-28 | 2018-07-06 | 中国石油天然气股份有限公司 | 岩石层理裂缝剪切参数预测方法及装置 |
CN108254262B (zh) * | 2016-12-28 | 2021-03-30 | 中国石油天然气股份有限公司 | 岩石层理裂缝剪切参数预测方法及装置 |
CN106596294A (zh) * | 2016-12-30 | 2017-04-26 | 中国科学院地质与地球物理研究所 | 一种模拟顺向边坡动力反应的试验装置及方法 |
CN107169236B (zh) * | 2017-06-15 | 2020-09-22 | 中国水利水电科学研究院 | 一种基于有限元与离散元耦合的虚拟三轴试验仿真方法 |
CN107169236A (zh) * | 2017-06-15 | 2017-09-15 | 中国水利水电科学研究院 | 一种基于有限元与离散元耦合的虚拟三轴试验仿真方法 |
CN107328920A (zh) * | 2017-08-15 | 2017-11-07 | 绍兴文理学院 | 矿山边坡岩体工程稳定性精准评价方法 |
CN107328919A (zh) * | 2017-08-15 | 2017-11-07 | 绍兴文理学院 | 矿山边坡岩体工程稳定性分级分析方法 |
CN107328920B (zh) * | 2017-08-15 | 2019-07-12 | 绍兴文理学院 | 矿山边坡岩体工程稳定性精准评价方法 |
CN107328919B (zh) * | 2017-08-15 | 2019-07-12 | 绍兴文理学院 | 矿山边坡岩体工程稳定性分级分析方法 |
CN107655691B (zh) * | 2017-08-28 | 2020-07-10 | 西北工业大学 | 一种飞行器进气道喘振锤击波试验装置及方法 |
CN107655691A (zh) * | 2017-08-28 | 2018-02-02 | 西北工业大学 | 一种飞行器进气道喘振锤击波试验装置及方法 |
CN108287112A (zh) * | 2018-01-31 | 2018-07-17 | 成都理工大学 | 一种基于三轴压缩试验测定岩石损伤参数的方法 |
CN108331042A (zh) * | 2018-03-05 | 2018-07-27 | 武汉科技大学 | 一种模拟***开挖对边坡影响的试验装置及试验方法 |
CN108331042B (zh) * | 2018-03-05 | 2023-10-13 | 武汉科技大学 | 一种模拟***开挖对边坡影响的试验装置及试验方法 |
CN108920754A (zh) * | 2018-05-25 | 2018-11-30 | 重庆交通大学 | 基于应***化剪切本构模型和位移变化的结构面控制边坡稳定性动态评价方法 |
CN109142671A (zh) * | 2018-07-17 | 2019-01-04 | 中国地质大学(武汉) | 研究***振动下饱水软弱结构面强度渐变劣化规律的方法 |
CN109142671B (zh) * | 2018-07-17 | 2021-04-30 | 中国地质大学(武汉) | 研究***振动下饱水软弱结构面强度渐变劣化规律的方法 |
CN109270243A (zh) * | 2018-07-26 | 2019-01-25 | 重庆交通大学 | 冻土***模型试验***及其试验方法 |
CN109270243B (zh) * | 2018-07-26 | 2021-11-30 | 重庆交通大学 | 冻土***模型试验***及其试验方法 |
CN109613607A (zh) * | 2018-12-04 | 2019-04-12 | 保利新联***工程集团有限公司 | 一种不同高差***地震波传播测试方法 |
CN109992841A (zh) * | 2019-03-11 | 2019-07-09 | 长江水利委员会长江科学院 | 一种***块度空间全方位分区耦合数值仿真方法 |
CN109992841B (zh) * | 2019-03-11 | 2022-12-06 | 长江水利委员会长江科学院 | 一种***块度空间全方位分区耦合数值仿真方法 |
CN110360907A (zh) * | 2019-07-31 | 2019-10-22 | 中冶集团资源开发有限公司 | 一种遍布节理破碎岩体的边坡台阶扩帮***方法 |
CN110502825A (zh) * | 2019-08-19 | 2019-11-26 | 青岛理工大学 | 一种提取三维破裂面的方法 |
CN111324942A (zh) * | 2019-12-27 | 2020-06-23 | 昆明理工大学 | 一种考虑滑面动力渐进破坏的地震边坡稳定性分析方法 |
CN111324942B (zh) * | 2019-12-27 | 2023-04-07 | 昆明理工大学 | 一种考虑滑面动力渐进破坏的地震边坡稳定性分析方法 |
CN111024588A (zh) * | 2019-12-31 | 2020-04-17 | 山东大学 | 一种体现渗流对岩土体强度弱化的dem接触模型构建方法 |
CN111199108A (zh) * | 2020-01-07 | 2020-05-26 | 包钢集团矿山研究院(有限责任公司) | 一种网格化切分数值模拟爆堆模型的***块度预测方法 |
CN113654419A (zh) * | 2020-05-12 | 2021-11-16 | 中国水利水电第三工程局有限公司 | 面板堆石坝过渡料***参数与最优级配曲线的分析方法 |
CN112213767A (zh) * | 2020-09-11 | 2021-01-12 | 中国矿业大学 | 一种用于顶板超前预裂***效果的评价方法 |
CN112182731A (zh) * | 2020-10-27 | 2021-01-05 | 上海工程技术大学 | 非极限状态二维边坡稳定性评价方法 |
CN113357987A (zh) * | 2021-04-15 | 2021-09-07 | 西安科技大学 | 一种自动***设计方法 |
CN113324831A (zh) * | 2021-06-16 | 2021-08-31 | 山东黄金矿业科技有限公司深井开采实验室分公司 | 一种露天矿边坡内采空区动力失稳破坏机理试验方法 |
CN114062127A (zh) * | 2021-10-29 | 2022-02-18 | 煤炭科学研究总院 | 岩体稳定性的检测方法及其装置 |
CN114062127B (zh) * | 2021-10-29 | 2022-06-28 | 煤炭科学研究总院有限公司 | 岩体稳定性的检测方法及其装置 |
CN115079248A (zh) * | 2022-05-26 | 2022-09-20 | 中国科学院武汉岩土力学研究所 | 基于裂纹密度及声波阈值反演***损伤范围的方法及设备 |
CN115079248B (zh) * | 2022-05-26 | 2024-06-07 | 中国科学院武汉岩土力学研究所 | 基于裂纹密度及声波阈值反演***损伤范围的方法及设备 |
CN116796453A (zh) * | 2023-04-18 | 2023-09-22 | 江汉大学 | 一种减振孔设计方法、装置、设备及可读存储介质 |
CN116796453B (zh) * | 2023-04-18 | 2024-03-08 | 江汉大学 | 一种减振孔设计方法、装置、设备及可读存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN105224742B (zh) | 2018-01-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105224742A (zh) | 一种分析***作用下台阶边坡稳定性的方法 | |
Silva et al. | Practical assessment of rock damage due to blasting | |
Banadaki | Stress-wave induced fracture in rock due to explosive action | |
Hao et al. | Numerical analysis of blast-induced stress waves in a rock mass with anisotropic continuum damage models part 1: equivalent material property approach | |
Hu et al. | Numerical simulation of the complete rock blasting response by SPH–DAM–FEM approach | |
Yu et al. | Experimental study of the mechanical behavior of sandstone affected by blasting | |
CN108331042A (zh) | 一种模拟***开挖对边坡影响的试验装置及试验方法 | |
Aliabadian et al. | Wave and fracture propagation in continuum and faulted rock masses: distinct element modeling | |
Siamaki et al. | Degradation of a discrete infilled joint shear strength subjected to repeated blast-induced vibrations | |
Sainoki et al. | Dynamic modelling of fault slip induced by stress waves due to stope production blasts | |
US11966672B2 (en) | Method and system for simulating contact and interaction between support member and chamber surrounding rock mass | |
Aliabadian et al. | Numerical modeling of presplitting controlled method in continuum rock masses | |
Hu et al. | A new safety control method of blasting excavation in high rock slope with joints | |
Yang et al. | Numerical and experimental study of the presplit blasting failure characteristics under compressive stress | |
CN116894291A (zh) | 基于fem-dem数值计算的不同地质隧洞围岩***振动安全判据构建方法 | |
Haghnejad et al. | Numerical investigation of the impact of geological discontinuities on the propagation of ground vibrations | |
Parra | Study of blast-induced damage in rock with potential application to open pit and underground mines | |
Jianhua et al. | Discussion on blasting vibration monitoring for rock damage control in rock slope excavation | |
Simioni et al. | Field measurements of snowpack response to explosive loading | |
CN207944494U (zh) | 一种模拟***开挖对边坡影响的试验装置 | |
Bakhshandeh Amnieh et al. | Numerical Analysis of the Primer Location Effect on Ground Vibration Caused by Blasting | |
Aydan et al. | Effect of blasting on the bench stability at the Demirbilek open-pit mine | |
Yi et al. | Numerical modelling of dynamic response of underground openings under blasting based on field tests | |
CN110427688A (zh) | 一种基于实测振动的地应力大小预测方法 | |
Tomori et al. | An integrated study on the response of unsupported underground cavity to the nearby construction of piles of Gushikawa By-Pass Bridge |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information |
Address after: 114001 Anshan District, Liaoning, No. 219 Road, No. 39, Tiedong Applicant after: Anshan Iron and Steel Group Mining Co., Ltd. Address before: 114001 Anshan District, Liaoning, No. 219 Road, No. 39, Tiedong Applicant before: Angang Group Mine Company |
|
COR | Change of bibliographic data | ||
GR01 | Patent grant | ||
GR01 | Patent grant |