创刊于1987年, 双月刊
主管:

江西理工大学

主办:

江西理工大学
江西省有色金属学会

ISSN:1674-9669
CN:36-1311/TF
CODEN YJKYA9

岩体变形及二次应力下矿体开采效应研究

向欣欣, 张耀平, 雷大星

向欣欣, 张耀平, 雷大星. 岩体变形及二次应力下矿体开采效应研究[J]. 有色金属科学与工程, 2024, 15(1): 87-96. DOI: 10.13264/j.cnki.ysjskx.2024.01.011
引用本文: 向欣欣, 张耀平, 雷大星. 岩体变形及二次应力下矿体开采效应研究[J]. 有色金属科学与工程, 2024, 15(1): 87-96. DOI: 10.13264/j.cnki.ysjskx.2024.01.011
XIANG Xinxin, ZHANG Yaoping, LEI Daxing. Study on ore-body mining effect under rock mass deformation and secondary stress[J]. Nonferrous Metals Science and Engineering, 2024, 15(1): 87-96. DOI: 10.13264/j.cnki.ysjskx.2024.01.011
Citation: XIANG Xinxin, ZHANG Yaoping, LEI Daxing. Study on ore-body mining effect under rock mass deformation and secondary stress[J]. Nonferrous Metals Science and Engineering, 2024, 15(1): 87-96. DOI: 10.13264/j.cnki.ysjskx.2024.01.011

岩体变形及二次应力下矿体开采效应研究

基金项目: 

江西省教育厅科学技术研究项目 GJJ218511

详细信息
    通讯作者:

    张耀平(1974— ),博士,教授,主要从事岩石力学及金属矿山地下开采工艺研究。E-mail:zyp2005csu@126.com

  • 中图分类号: TD803

Study on ore-body mining effect under rock mass deformation and secondary stress

  • 摘要: 为了解围岩在变形及二次应力下矿体开采效应的影响,本文选取-410 m中段采用现场试验与数值模拟相结合的方法对崩落法采矿进行矿体效应研究。结果表明:围岩应力发生震荡变化,在矿体上下盘围岩附近存在应力明显降低区域,下盘大于上盘,顶板最大主应力达0.44 MPa,局部呈现应力集中,且应力集中部位呈现类似水平零散分布的现象。同时开采边界呈现裂隙发育的状态。水平矿体有轻微下沉趋势,且在该中段的水平巷道边界和顶板裂隙附近易出现断层破碎带。竖直方向变形明显,且倾斜部位呈椭圆分布,最大下沉位移值达8.00 mm,最大上升方向变形值为3.03 mm,交界线附近易受到应力挤压现象。该研究对及时准确地掌握岩体变形规律和矿山安全信息使矿山安全有效合理地进行开采具有一定的指导意义。
    Abstract: In order to understand the influence of the mining effect of the surrounding rock under deformation and secondary stress,this paper selected the -410 m section to study the ore body effect of caving mining by combining field test and numerical simulation.The results show that the surrounding rock stress changes in shock,there is a stress reduction area near the upper and lower wall surrounding rock of the orebody,the foote is larger than the hanging wall,the maximum principal stress of the roof reaches 0.44 MPa,and local stress concentration.In addition,the stress concentration parts showed similar horizontal scattered distribution phenomenon.At the same time,the mining boundary presents the state of crack development.The horizontal ore body has a slight subsidence trend,and the fault fracture zone is easy to appear near the boundary of the horizontal roadway and the roof crack in the middle section.The deformation in the vertical direction is obvious,and the inclined part is elliptical distribution.The maximum subsidence displacement is 8.00 mm,and the maximum upward deformation is 3.03 mm.The stress extrusion phenomenon is easy to occur near the intersection line.The study has certain guiding significance to grasp the deformation law of rock mass and mine safety information in time and accurately so as to carry out mining safely,effectively and reasonably.
  • 7075铝合金是典型的航空用铝合金,具有强度高、质量轻、耐腐蚀等多重特性,对组织和性能要求较高[1-2]。近年来,7075铝合金不断实现民用化,已经应用于制造汽车零部件[3],在代替部分汽车用高强钢方面应用潜力巨大。7075铝合金的塑性变形工艺以挤压为主,目前有关7075铝合金轧制工艺的研究还较少,有关热处理工艺[2,4]和应力腐蚀性能[5]的研究较多。因此,有必要对7075铝合金的轧制工艺开展研究,找到实现批量化生产7075铝板的较优轧制工艺,进一步提高7075铝板的生产效率,优化生产成本。异步轧制是通过加入非对称元素,在轧制变形区内形成“搓轧区”,产生较大的剪切变形量[6]。在异步轧制基础上,又发展出了水平错位等速轧制和龙形轧制[7]等非对称轧制方法。异步轧制在极薄带材轧制方面表现出较大潜力,东北大学在20世纪80年代就采用异步轧制技术用直径90 mm的工作辊轧制出0.005 mm厚度的铜箔[8]。在铝板异步轧制方面,有实验表明异步轧制对AA 5454铝板力学性能改善不明显[9],但KAZEMI-NAVAEE等[10]的实验证实异步轧制对冷轧7075铝板的屈服强度和抗拉强度有提升作用,可以作为7075铝板轧制工艺改善的探索方向。以上实验的共同特点是轧制累积变形率均不足70%,因此有必要进行更高轧制累积变形率的异步轧制实验,并与常规轧制的铝板的性能对比,为进一步开发更高效的7075铝板轧制工艺提供借鉴。

    本文采用异速轧制和常规轧制2种方式,以不同轧制累积变形率(最高85%)对7075铝板进行轧制对比实验。考察了异速比对单道次异速轧制轧板弯曲程度的影响;采用相同道次分配和异速比,在不同轧制累积变形率的情况下,对比研究了异速轧制和常规轧制对7075铝板变形情况、晶粒组织、断口形貌和力学性能的影响。

    本文所用轧制板坯为7075铝合金T6态成品板,经锯切和铣面后获得板坯尺寸为300 mm × 100 mm × 20 mm(长 × 宽 × 厚)。实验所用轧机为二辊轧机,可在线调整异速比。坯料加热采用台车式电阻炉,加热工艺为440 ℃保温6 h。

    实验分2步进行:第1步进行异速比对轧板弯曲情况影响的轧制实验;第2步进行异速轧制和常规轧制的多道次轧制实验。第1步实验采用的异速比分别为1.02、1.06、1.12和1.16,每个异速比单道次轧制1块铝板,道次压下量设为4 mm、变形率20%,轧后立即水淬,再进行轧板弯曲程度测量。第2步进行多道次轧制实验,分2组进行,第1组为常规轧制,第2组为异速轧制,道次分配相同,如表1所列,异速轧制下轧辊转速快于上轧辊。每组实验轧制4块板,分别轧3、4、5、6道次,设计累积变形率分别为60%、75%、85%、95%。轧后铝板立即水淬,再进行T6热处理,工艺为460 ℃固溶1 h,再在120 ℃进行24 h时效处理。

    表  1  道次分配参数和设计累积变形率
    Table  1.  Pass allocation parameters and design cumulative deformation rates
    轧制道次设计道次变形量/mm设计累积变形率/%
    1420
    2440
    3460
    4375
    5285
    6295
    下载: 导出CSV 
    | 显示表格

    金相试样经磨抛后,采用4% HF + 6% HCl + 10% HNO3 + 80% H2O溶液侵蚀,在Observer 3 m ZEISS光学显微镜下观察金相组织。断口扫描电镜观察采用日立 SU5000热场式场发射扫描电镜, 力学性能测试采用 MTS 810 力学实验机。屈服强度和抗拉强度测量执行GB/T 228—2010《金属材料拉伸试验室温试验方法》,因为异速轧制后轧板呈现一定程度的弯曲,尽管只在较平直的位置截取拉伸试样,仍无法实现试样完全平直,因此,本文中异速轧制轧板的屈服强度和抗拉强度仅供对比参照。

    因为轧辊转速不同,轧板与轧辊接触的表面金属流动速度也不同,造成异速轧制轧板向慢速辊方向弯曲,所以对异速轧制轧板弯曲程度的表征和分析十分重要,也是改进工艺、克服弯曲所建立的数学模型的重要基础。为尝试建立轧板弯曲程度与异速比的数学关系,将轧板水平放置,使其开口朝向水平底面,轧板弯曲内面形成一个圆弧,测量该圆弧最高点距离水平底面距离y,以及圆弧即轧板弯曲内面与水平底面的两接触点距离x,用y/x值描述轧板弯曲程度[11],测量方法如图1(a)所示。

    图  1  轧板y/x值测量示意图和变化情况:(a)轧板y/x值测量示意;(b)y/x值随异速比变化的情况
    Figure  1.  Measurement and change of y/x value of rolled plates:(a)measurement of y/x value of rolled plates;(b)change of y/x value with different speed ratio

    所得y/x值随异速比变化的情况见图1(b),异速轧制轧板y/x值随着异速比增大而先升高后降低,在异速比1.12时达到最高0.222 18,之后在异速比1.16时降至0.191 14。JIANG等[12]对异步轧制厚钢板进行了有限元模拟(压下率12.8%),当异速比达到1.15时,轧板弯曲程度最大。由于轧机能力所限,本研究未能进行更大异速比的轧制,因此无法直接判断随着异速比进一步增大,轧板弯曲程度是否会持续减弱甚至消失。但是,采用其他工艺方法能够减弱甚至消除异速轧制轧板弯曲已被证实。例如龙形轧制[13],亦称蛇形轧制,是一种非对称轧制,在异步轧制的基础上,将转速不同的两辊在水平方向上前后偏移,对板材增加一个弯矩以抵消可能的弯曲,通过调节异速比、错位量、压下量等参数,可以实现无弯曲的轧制[14]。郑细昭等[15]采用非线性回归得到了龙形轧制中轧板弯曲控制的最佳错位量的数学模型。

    目前,轧板弯曲程度多采用弯曲曲率表征。FU等[16]建立了龙形轧制轧板曲率预测的解析模型,针对轧板曲率变化进行了大量计算。相比于曲率法,采用y/x值表征轧板弯曲程度的操作相对简便,可以较快地建立弯曲程度与异速比的数学关系。与用轧板的弯曲曲率表征轧板弯曲程度相同,采用y/x值表征轧板的弯曲程度同样是将轧板弯曲的弧线设想为一个标准圆形的其中一段,使得计算简化。但是,对于轧板实际弯曲的弧线形状以及适用的数学表达,目前还无具体研究结果。

    结合第1步实验的结果,异速比1.02对轧板的影响较小,不足以反映异速轧制的特征;异速比1.12导致轧板弯曲太大,进行多道次轧制存在咬入困难,并且可能损坏设备;异速比1.16已经达到设备可调的上限。因此,为顺利进行实验和保护设备安全,选择异速比1.06进行第2步实验。第2步实验中异速轧制铝板的表面形貌如图2所示。

    图  2  异速轧制铝板的表面形貌
    Figure  2.  Surface morphology of asymmetrical rolling plates

    为对比异速轧制和常规轧制后轧板的实际累积变形率是否存在差异,异速轧制和常规轧制采用了相同的道次分配方案。由于轧机存在一定的弹跳现象以及轧制能力有限,导致轧板的实际累积变形率达不到设计累积变形率。因此,本研究中每道次轧制的设计累积变形率偏大,以求达到预想的实际累积变形率,每道次轧制的设计累积变形率和实际累积变形率如图3所示。

    图  3  轧板实际累积变形率的变化情况
    Figure  3.  Changes of actual accumulative deformation ratio of rolled plates

    图3可见,常规轧制轧板的实际累积变形率与设计累积变形率相差8%~10%,这个数值在第3道次至第6道次都比较稳定。异速轧制轧板的实际累积变形率在第3道次至第5道次与设计累积变形率相差6%~8%,第6道次时相差10%。异速轧制和常规轧制轧板每道次的实际累积变形率相差不大,在1%~2%之间。

    在相同设计累积变形率的情况下,异速轧制与常规轧制所得的实际累积变形率不同的结果此前还未见报道。在同样存在异速轧制的情况下,雷军义等[17]在对龙形轧制的研究中提出辊缝补偿量的概念,即龙形轧制与常规轧制的实际变形量存在差异,并认为该辊缝补偿量与异速轧制无关。本研究中异速轧制与常规轧制所得的实际累积变形率相差不大,所以后续实验仍将同道次的两种轧板进行比较。

    异速轧制和常规轧制的轧板经T6热处理后,沿厚度方向在轧板中心部切取平行于轧向的侧面观察晶粒组织。如图4所示,异速轧制和常规轧制在第3道次时的T6板中心部以变形组织为主,此时实际累积压下率约为51%,原始板坯中的晶粒组织被破坏,只能看到垂直于轧制方向的条带状的变形组织。对比图4(a)和图4(e)可知,3道次异速轧制的T6板中心部变形更严重,组织变形更明显,变形条带更清晰、密集,但未观察到文献[10]中描述的剪切变形带。因为累积的变形储能有限,异速轧制和常规轧制3道次的轧板都无法在T6热处理过程中实现充分完全的再结晶,在T6热处理后仍存在明显的变形组织。这种以变形组织为主的特征在4道次轧制的T6板中没有明显改变。

    图  4  T6板晶粒组织(图中箭头所示为轧制方向):(a)常规轧制3道次;(b)常规轧制4道次;(c)常规轧制5道次;(d)常规轧制6道次;(e)异速轧制3道次;(f)异速轧制4道次;(g)异速轧制5道次;(h)异速轧制6道次
    (b)symmetrical rolling 4 passes; (c)symmetrical rolling 5 passes; (d)symmetrical rolling 6 passes; (e)asymmetrical rolling
    Figure  4.  Microstructure of T6 plates(the arrows in the figures show the rolling direction):(a)symmetrical rolling 3 passes;
    3 passes; (f)asymmetrical rolling 4 passes; (g)asymmetrical rolling 5 passes; (h)asymmetrical rolling 6 passes

    TAJALLY等[18]的研究表明,当冷轧变形量达到71%时,在265 ℃保温条件下仅需70 min即可完成再结晶。在图4(c)和图4(g)中,异速轧制和常规轧制在第5道次的T6板中心部已经可以看到变形组织减少,再结晶晶粒开始增加。经过6道次异速轧制和常规轧制的轧板,实际累积变形率分别达到83%和85%,轧制变形严重,累积了较充分的变形储能,而变形储能是再结晶的动力,因此在T6热处理过程中发生了较充分的静态再结晶。在图4(d)和图4(h)中,异速轧制和常规轧制的轧板经过6道次轧制和T6热处理后,中心部变形组织已经基本消除,转变为长条状的晶粒组织,尽管仍有垂直于轧制方向伸长的特征,但已可见比较清晰的晶粒和晶界。7075是Al-Zn-Mg-Cu合金,其主要强化模型为沉淀强化,细晶强化即晶界强化的效果较弱[19]

    图5所示为异速轧制和常规轧制第3道次和第6道次轧板T6热处理后的拉伸断口的SEM形貌。在7075铝合金固溶过程中,部分第二相(如Al7Cu2Fe)不会溶解,随着保温时间延长,Mg2Si会逐渐粗化长大[20]。结合图4可知,异速轧制和常规轧制的轧板第3道次的T6板中因为主要是变形组织,所以图5(a)和图5(d)中都有垂直于轧制方向的大量撕裂棱,局部可见由于大尺寸第二相颗粒造成的孔洞。由图5(b)和图5(e)可见,垂直于轧制方向的撕裂棱明显减少,以大小不一的韧窝断口和沿晶断口为主,是未溶相粒子和时效过程中弥散析出的第二相粒子共同存在的结果;虽然也可见第二相颗粒造成的孔洞,但相比图5(a)和图5(d)中都明显减少,是变形破碎大尺寸第二相颗粒的结果。对比可见,图5(e)中韧窝更多,大小韧窝相套,韧窝尺寸和深度的分布比图5(b)中更均匀。

    图  5  T6板断口SEM形貌:(a)常规轧制3道次;(b)常规轧制6道次;(c)常规轧制6道次,局部高倍形貌, 箭头所示为第二相粒子;(d)异速轧制3道次;(e)异速轧制6道次;(f)异速轧制6道次,局部高倍形貌, 箭头所示为第二相粒子
    Figure  5.  Fracture SEM morphology of T6 plates: (a)symmetrical rolling 3 passes; (b)symmetrical rolling 6 passes; (c)symmetrical rolling 6 passes, partial high magnification, the arrows show the second phase particles;(d)asymmetrical rolling 3 passes;(e)asymmetrical rolling 6 passes;(f)asymmetrical rolling 6 passes, partial high magnification, the arrows show the second phase particles

    图5(c)为图5(b)局部的高倍形貌,图5(f)为图5(e)局部的高倍形貌,图5(c)和图5(f)中均可见大量尺寸小于 1 μm 的韧窝,并且在韧窝底部出现弥散析出的尺寸小于 100 nm 的微小第二相粒子,图5(f)中的韧窝尺寸较图5(c)中更小。

    异速轧制和常规轧制的轧板经T6热处理后的屈服强度和抗拉强度如图6所示,异速轧制轧板T6热处理后的屈服强度和抗拉强度均高于常规轧制的轧板,但提高幅度不大。此结果与KRANER等[21]的铝板异步轧制实验结果相似,但对纯铝板进行的异步轧制实验发现,异步轧制轧板的强度明显高于对称轧制[22]

    图  6  异速轧制和常规轧制T6板强度变化
    Figure  6.  Strength changes of asymmetrical and symmetrical rolling T6 plates

    同时,常规轧制第4道次的轧板T6热处理后的强度最高,第5道次和第6道次轧制的轧板T6热处理后的强度不再上升,而异速轧制的轧板T6热处理后的强度随轧制道次增加,轧制累积变形率持续升高。据文献[23]报道,7075铝合金轧板在常规轧制累积变形率达到80%时,强度不再升高反而降低。本研究中常规轧制第4道次时累积变形率为67%,此后随累积变形率升高,T6热处理后的轧板强度不再明显变化。结合文献[23]与本研究中常规轧制的结果可知,7075铝合金轧板在常规轧制达到一定变形量时会出现强度峰值,出现峰值的变形量可能因板坯状态、均热条件和道次分配等的差异而在不同的实验中存在差异。但是,异速轧制的7075铝板T6热处理后的强度随轧制累积变形率升高而持续增大,未出现强度达到峰值后不再升高的情况(如图6所示)。

    本文采用4个不同的轧制累积变形率,最高轧制累积变形率达到80%以上,对比研究了异速轧制和常规轧制对7075铝板变形情况、晶粒组织、断口形貌和力学性能的影响,得到以下结论:

    1)当轧制异速比从1.02升至1.12时,y/x值达到最高0.222 18,此时轧板弯曲程度最大,之后在异速比1.16时y/x值降至0.191 14,即轧板弯曲程度较异速比1.12时减弱。

    2)异速轧制和常规轧制采用了相同的设计累计变形率,但轧板实际累计变形率相差1%~2%。

    3)异速轧制和常规轧制的7075铝板在累积实际变形率分别达到67%和69%时,T6热处理后的晶粒组织以变形组织为主;在累积实际变形率达到约76%时,随着累积实际变形率增加,再结晶更加充分,T6热处理后的7075铝板中以伸长的条状晶粒为主。

    4)经过6道次异速轧制和常规轧制的7075铝板,T6热处理后,在断口中可见大量尺寸小于1 μm的韧窝,并且在韧窝底部出现弥散析出的尺寸小于100 nm的微小第二相粒子,异速轧制7075铝板中韧窝尺寸较常规轧制7075铝板中更细小。

    5)与常规轧制相比,采用异速轧制的7075铝板的强度高,但提高幅度不大;7075铝合金轧板在常规轧制达到一定变形量时会出现强度峰值,而异速轧制的7075铝板的强度随轧制累积变形率升高而持续提高。

    赵中波
  • 图  1   光应力计受力示意

    Fig  1.   Schematic diagram of the force of the photostress meter

    图  2   -410 m中段监测点布设示意

    Fig  2.   Schematic diagram of the layout of mid-level monitoring points at -410 m

    图  3   开挖过程应力变化

    Fig  3.   Stress variation diagram during excavation

    图  4   各测点的应力-时间变化曲线

    Fig  4.   Stress-time curve of each measuring point

    图  5   -410 m水平位移变形曲线

    Fig  5.   -410 m horizontal displacement deformation curve

    图  6   -410 m沉降变形曲线

    Fig  6.   -410 m settlement deformation curve

    图  7   崩落法水平采矿进度及地压显现规律

    Fig  7.   Progress of horizontal mining in caving method and ground pressure manifestation law

    图  8   FLAC3D三维数值模型

    Fig  8.   FLAC3D numerical model diagram

    图  9   z方向位移云图(a)和顶板沉降变形(b)

    Fig  9.   z direction displacement nephogram(a) and roof settlement deformation(b)

    图  10   x方向位移云图:(a) 原有开采活动位移云图;(b) 开采过程位移云图

    Fig  10.   Displacement cloud map in the x direction:(a) original displacement cloud map of mining activities; (b) displacement cloud map during mining proces

    图  11   水平最小主应力等值线云图

    Fig  11.   Horizontal minimum principal stress contour nephogram

    图  12   水平最大主应力等值线云图和巷道岩体变形:(a) 水平最大主应力等值线云图;(b) 巷道岩体变形

    Fig  12.   Contour nephogram of maximum horizontal principal stress and deformation of roadway rock mass:(a) cloud map of horizontal maximum main stress; (b) rock mass deformation of roadway

    表  1   各种介质的力学参数

    Table  1   Mechanical parameters of various media

    介质名称密度/(kg/m3弹性模量/GPa泊松比内聚力/MPa单轴抗拉强度/MPa内摩擦角/(°)
    矿岩3 30031.600.283.803.547.0
    闪长岩2 73926.600.292.502.839.0
    大理岩2 70718.000.222.101.646.6
    尾胶充填体1 8910.360.250.170.538.0
    尾砂充填体1 7500.040.100.050.217.5
    下载: 导出CSV
  • [1] 李小瑞,侯公羽,吕文涛,等.巷道围岩二次应力状态弹性阶段的试验研究[J]. 采矿与安全工程学报,2021,38(2):269-275.
    [2] 刘允秋,李光明,徐放艳,等.深部高应力不稳固矿体开采方法研究[J]. 矿业研究与开发,2019,39(7):15-18.
    [3] 杨丹,刘洋.矩形巷道围岩弹性变形能积聚特征分析[J].力学与实践,2020,42(3):300-305.
    [4] 赵呈星,李英明.基于围岩变形的预应力锚杆受力特征分析及锚固机理[J]. 煤矿安全,2020,51(7):234-238.
    [5] 赵国彦,马驰,郭琪.深部硬岩二次应力状态下破裂的断裂力学分析[J]. 科技导报,2012,30(31):47-52.
    [6] 颜杜民.全断面隧道施工围岩变形多重分形特征研究[J].铁道建筑,2020,60(11):56-58.
    [7] 龚林金,任锐,王亚琼,等.隧道斜穿不同倾角断层破碎带围岩变形特征分析[J]. 公路,2021,66(7):313-319.
    [8] 王亚琼,杨强,潘红伟,等.基于3DEC模拟的高地应力水平层状隧道围岩变形破坏特征分析[J]. 现代隧道技术,2022,59(4):127-136.
    [9] 邓兆睿,孙德全,任仲久,等.近距离煤层群半煤岩上保护层开采卸压增透效应研究[J]. 采矿与岩层控制工程学报,2022,4(2):35-43.
    [10] 王学滨,岑子豪,薛承宇,等.正断层下盘开采应力演化及断层阻隔效应连续-非连续方法模拟[J]. 地球物理学进展,2022,37(6):2612-2621.
    [11] 刘雨涛,信连凯,谌志勇.高山矿近距离低透煤层上保护层开采底板卸压增透效应研究[J]. 煤矿安全,2020,51(3):32-37.
    [12] 王方田,李岗,班建光,等.深部开采充填体与煤柱协同承载效应研究[J]. 采矿与安全工程学报,2020,37(2):311-318.
    [13] 樊克松,申宝宏,张风达,等.厚煤层开采地表移动变形规律的深厚效应研究[J]. 煤炭科学技术,2018,46(3):194-199.
    [14]

    JIANG F F, ZHOU H,SHENG J,et al.Temperature effect on shear behavior of ore-backfill coupling specimens at various shear directions[J]. Journal of Central South University,2021,28(10):3173-3189.

    [15] 易璐,郑明贵. 中国稀土开采总量控制政策效应评估[J].有色金属科学与工程,2021,12(2):120-126.
    [16] 邹国良,刘娜娜. 基于组合赋权-云模型的离子型稀土矿开采工艺评价[J]. 有色金属科学与工程,2021,12(4):88-95.
    [17] 吴春成. 川藏铁路巴玉隧道开挖二次应力分布特征与岩爆风险分析[J]. 铁道标准设计,2021,65(7):131-135.
    [18] 郭延辉, 侯克鹏, 李克钢, 等. 不同原岩应力下厚度对急倾斜矿体开采岩移规律的影响[J]. 昆明理工大学学报(自然科学版),2020,45(2):31-37.
    [19] 张磊, 许琳, 李朝, 等. 煤层群开采垂直应力分布的层间岩性效应研究[J]. 煤炭技术,2021,40(11):31-34.
    [20] 李小双, 罗浪, 王运敏,等. 露天转地下开采覆岩采动响应的坡角影响效应研究[J]. 矿冶工程,2021,41(4):33-37.
    [21]

    CHEN M,ZANG C W,DING Z W,et al.Effects of confining pressure on deformation failure behavior of jointed rock[J].Journal of Central South University,2022,29(4):1305-1319.

    [22] 唐维军, 孙希奎, 王恒, 等. 膏体充填开采条带煤柱覆岩稳定效应研究[J]. 煤炭科学技术,2017,45(09):109-115.
    [23] 朱广安, 刘博文, 窦林名, 等. 基于采空区压实效应的工作面开采全过程模拟[J]. 中国矿业大学学报,2019,48(4):775-783.
    [24] 廖宝泉,柯愈贤,方立发,等.膏体充填开采覆岩移动变形规律研究[J].有色金属科学与工程,2022,13(3):99-105.
    [25]

    XUE Y C,XU T,WASANTHA P L P,et al.Control of dynamic disaster by backfill mining under thick magmatic rock in one side goaf[J].Journal of Central South University,2020,27(10):3103-3117.

  • 期刊类型引用(1)

    1. 刘志芳,杨舜,王志,唐德馨. 降雨入渗与采动卸荷耦合下裂隙岩体力学特性及裂隙演化规律研究. 有色金属科学与工程. 2024(04): 588-597 . 本站查看

    其他类型引用(0)

图(12)  /  表(1)
计量
  • 文章访问数:  65
  • HTML全文浏览量:  10
  • PDF下载量:  9
  • 被引次数: 1
出版历程
  • 收稿日期:  2022-11-23
  • 修回日期:  2023-03-24
  • 刊出日期:  2024-02-28

目录

/

返回文章
返回
x 关闭 永久关闭

敬  告

为进一步规范期刊查阅、引用、统计等,现对我部编辑出版的《有色金属科学与工程》英文刊名全称和英文刊名缩写公告如下:

英文刊名全称为:Nonferrous Metals Science and Engineering

英文刊名缩写为:Nonferr.Met.Sci.Eng.