Analysis on stability of ionic rare earth slope under heavy rainfall
-
摘要: 降雨是诱发离子型稀土矿山滑坡灾害的主要因素之一。以赣南某离子型稀土矿山边坡为研究对象,采用数值模拟方法得到强降雨下边坡孔隙水压力的变化情况,分别分析降雨类型、停雨间隔时间和降雨强度对边坡孔隙水压力的影响,并应用有限元强度折减法计算拟定工况下的边坡稳定性系数。结果表明:①不同降雨类型下边坡孔隙水压力变化情况不同,增强型降雨对边坡稳定性系数影响最大。②边坡表面孔隙水压力回落幅度随停雨间隔时间的增加而增大。总降雨量相同时,停雨间隔时间缩短对边坡稳定性不利。③随着降雨强度的增大,边坡孔隙水压力的变化速率增大,边坡稳定性系数下降速率增大并且下降幅度增大。Abstract: Rainfall is one of the main factors inducing the landslide disasters in ionic rare earth mines. In this paper, with the slope of an ionic rare earth mine in southern Jiangxi as the research object, the changes in pore water pressure of the slope under heavy rainfall were obtained by a numerical simulation method. The influence of rainfall typ, the interval time between rain stops and rainfall intensity on the slope pore water pressure were analyzed respectively, and the coefficient of the slope stability was calculated by applying finite element strength subtraction under the proposed working conditions. The results show that pore water pressures of the slope vary under different rainfall types, and the enhanced rainfall has the greatest influence on the slope stability coefficient. The drop amplitude of the pore water pressure on the slope surface increases with the increase of rain stop interval. When the total rainfall is the same, the shortened rain stop interval is detrimental to slope stability. With the increase of rainfall intensity, the change rate of slope pore water pressure increases, the decline rate of slope stability coefficient increases, and the decline range increases.
-
Keywords:
- rainfall infiltration /
- ionic rare earth /
- numerical simulation /
- slope stability
-
钨基高比重合金(WHAs)因具有高密度、高强度、高硬度、吸收射线能力强、导电导热性好、耐腐蚀性强以及良好的延展性等一系列优异的物理力学性能,使得其作为配重块、动能穿甲弹弹芯和屏蔽材料等,在国防工业、航空航天和民用工业等领域中得到广泛的应用,是一种非常重要的军民两用材料[1-3].
传统钨基高比重合金的制备是通过将合金粉末混合,经模压成型后再进行液相烧结而成. 液相烧结温度较高(通常在1 480~1 600℃之间),此时得到的钨颗粒尺寸粗大,一般在30~60 μm,因此合金性能低下[4]. 随着科学技术的迅猛发展,对钨合金的使用性能提出许多更高更新的要求,传统钨基高比重合金的性能已经无法满足现代材料向小型化、智能化、高集成、高密度存储与快速传输化的发展需求. 研究表明,钨晶粒的细化对于钨基高比重合金的综合性能,尤其是其穿甲“自锐”性能的提高效果显著[5-6]. 因此,高强韧、结构纳米化高比重钨合金的研究开发成为其最重要的发展方向[7-9]. 许多材料工作者尝试采用添加合金元素或者变形加工的方法来提高合金的某些性能和细化晶粒,但在全面改善其综合性能方面仍无突破性进展,大大限制其广泛应用[10-15]. 放电等离子体烧结(Spark Plasma Sintering,SPS)是通过在粉末颗粒之间施加脉冲电流进行加热和加压,从而实现材料烧结致密化,是一种新型的粉末热压烧结新技术. 具有加热均匀、升温速度快、烧结温度低、烧结时间短、冷却迅速、生产效率高、节能环保等独特优势[16-18],为获得细小均匀的组织、改善钨基高比重合金的性能提供一条全新的思路. 目前,已有部分学者尝试采用SPS技术制备高强韧、超细晶钨基高比重合金,但绝大部分研究仅仅局限于合金显微组织结构和力学性能方面[19-20],而对于合金在特殊用途过程中所必须具备的耐磨性能方面的研究则鲜有报道.基于此,论文作者重点研究球磨工艺对SPS W-4.9Ni-2.1Fe合金摩擦磨损特性的影响,为高强韧钨基高比重合金的研究开发提供翔实的基础数据和可靠的理论指导.
1 实 验
实验采用微米级W、Ni、Fe元素粉末作为原料粉末,按W-4.9Ni-2.1Fe(质量分数,%)的成分配比进行机械预混合后,于QM-3SP04型球磨机中进行湿磨(无水乙醇环境),球磨时间为2~48 h,球磨罐材质为304不锈钢,球磨介质为硬质合金球,球料比为5∶1,转速为400 r/min,球磨过程中通入氩气以防止合金粉末氧化. 球磨后的合金粉末放入真空干燥箱中(70 ℃)进行干燥.
不同球磨工艺得到的合金粉末的致密化过程在SPS-10T-5型真空放电等离子体烧结系统(上海晨华电炉有限公司)上进行,具体烧结工艺如下:烧结温度为1 250 ℃,保温时间为5 min,烧结压力为35 MPa,温升率为100 ℃/min,每次烧结用粉量约20 g,烧结气氛为真空,真空度<10 Pa,烧结完成后随炉冷却.
采用阿基米德排水法精确测量试样的密度;于HR-150A型洛氏硬度计(济南华银试验仪器有限公司)上测量烧结合金的洛氏硬度;利用荷兰帕纳科Empyrean型衍射仪对不同球磨时间的合金粉末进行物相分析,管电流为30 mA,管电压为40 kV;使用HSR-2M型高速往复摩擦磨损试验机测量试样的摩擦磨损性能,加载载荷为700 g,往复长度为5 mm,加载时间为60 min.对偶摩擦副采用直径5mm的Si3N4球. 实验前,用超声波和酒精溶液清洗试样和摩擦副. 摩擦磨损试验后,采用NanoMap-500LS型扫描三维轮廓仪测量试验的磨痕轮廓(扫描距离1 000 μm,扫描速度100 μm/s),计算体积磨损量,并采用TM3030型扫描电子显微镜观察烧结后合金显微组织结构以及磨损后的表面组织和磨痕形貌.
2 结果与讨论
2.1 物相分析
图 1所示为经过不同球磨时间得到的合金粉末的XRD图谱. 从图 1中可以发现,当球磨时间较短(2 h)时,图谱中不仅存在主元素W的衍射峰,而且出现Ni、Fe元素的衍射峰,表明此时Ni、Fe元素仍以单质形式存在.当球磨时间为24 h时,Ni、Fe峰消失,仅存在W元素的衍射峰.同时可以发现,随着球磨时间的逐渐延长,W衍射峰强度逐渐降低并且产生明显的宽化,究其原因是由于Ni(Fe)溶入W晶格中形成W的过饱和固溶体,同时长时间的球磨过程使晶粒细化,晶粒内部产生严重的晶格畸变,最终导致Ni、Fe峰的消失和W峰衍射线条的宽化.
2.2 显微组织观察
钨基高比重合金的微观组织是由黏结相Ni(Fe)固溶体包围着硬质相W颗粒所构成. W颗粒与W颗粒之间直接接触的地方是断裂最容易发生处. 因此,使黏结相均匀分布在W颗粒周围,减少W颗粒之间的直接接触,降低W-W颗粒连接度,是制备高强韧细晶钨合金的有效途径[21]. 图 2所示为经过不同球磨时间的W-4.9Ni-2.1Fe合金粉末SPS烧结态合金SEM形貌. 从图 2中可知,合金呈现出明显的固相烧结态组织,为黏结相和硬质相双相组织结构. 当球磨时间为2 h时,合金中的Ni(Fe)黏结相量相对较少且分布不均匀(图 2(a)). 原因主要有以下2点:一方面是由于球磨时间短,Ni、Fe元素仍以单质形式存在,未溶入W晶格中形成W的过饱和固溶体,这一点从XRD分析中可以得到印证;另一方面,SPS烧结时间和保温时间均很短,黏结相无法得到充分扩散,致使合金中的黏结相仅局部形成,而不能有效地将W颗粒包裹起来,最终导致W颗粒与黏结相分布极不均匀. 进一步观察可以看出,随着球磨时间的逐渐延长,黏结相分布趋于均匀化(图 2(b)~(d)),有效地避免了钨颗粒之间的直接接触.
密度、硬度分析
图 3所示为SPS烧结态合金试样的相对密度随球磨时间的变化曲线. 从图 3中可以看出,当球磨时间为2 h时,合金试样的相对密度最高,达98 %. 随着球磨时间的延长,球磨的破碎效应占据主导作用,粉末晶粒尺寸逐渐变细小,相应的比表面积显著增大,从而提高球磨粉末与空气之间的接触面积,导致粉末颗粒表面吸附大量的气体杂质.球磨时间越长,吸附效果越明显,吸附的气体杂质量越大. 吸附的气体杂质在随后的SPS烧结过程中无法完全脱除而残留在合金内部,形成闭孔隙,最终导致合金试样相对密度的降低. 同时,长时间的球磨还会使球磨罐的成分元素进入合金粉末当中,改变合金成分结构,从而影响合金粉末的性能. 因此,当球磨时间从2 h增加到48 h时,合金试样的相对密度从98 %下降至85 %,下降幅度约15 %.这充分表明适宜的球磨时间对获得高性能的细晶钨基高比重合金是至关重要的.
图 4所示为合金试样的洛氏硬度值随球磨时间的变化关系曲线. 从图 4中可以发现,经过高能球磨和SPS烧结后,试样硬度值较高,均明显大于70 HRA. 当球磨时间为2 h时,合金试样的硬度值最大,可达78.6 HRA,这种异常现象的产生可能与硬质合金研磨球中高硬度脆性相WC的脱落有关,这一点在Xiang D.P.的研究过程中也有类似发现[22]. 当球磨时间从24 h延长至48 h时,合金试样硬度值从71.4 HRA逐渐提高到77.3 HRA,增长幅度约7 %. 硬度值的这种变化规律与合金试样中W晶粒的细化程度是息息相关的. 随着球磨时间的延长,晶粒细化明显,硬度值相应增大. 但长时间的研磨,又会使合金粉末中的杂质元素和氧含量增大,阻碍烧结过程中烧结颈的形成和长大,使合金粉末致密化过程受阻,合金相对密度下降.
2.4 摩擦磨损行为分析
图 5显示的是不同球磨时间的SPS烧结态合金摩擦磨损系数随磨损时间的变化关系曲线. 从图 5中可以看出,随着球磨时间的延长,合金试样的摩擦系数呈现出先下降后升高的趋势,且当球磨时间为24 h时,合金试样的摩擦系数最低,其平均值仅为0.48. 此外,通过观察磨损曲线随时间的变化趋势可以发现,球磨2 h的SPS烧结态合金其摩擦系数随着磨损时间的延长出现较大波动,36 h和48 h的波动次之,而24 h的烧结态合金摩擦系数曲线波动最小,摩擦系数值比较稳定.
图 6所示为不同球磨时间的SPS烧结合金摩擦磨损试验后的磨损形貌SEM像. 由图 6可以明显看出球磨时间为24 h时,合金试样的磨痕宽度最窄,仅为293 μm,磨痕表面因长时间的摩擦,受磨损区域在接触应力作用下,表面产生裂纹,开始形成片状脱落,但只有少数区域脱落. 而相较于球磨24 h的合金试样磨痕形貌,球磨时间分别为2 h、36 h以及48 h的烧结合金,其磨痕宽度较大,分别为581 μm、585 μm和659 μm,明显大于球磨24 h的烧结合金,且其磨痕形貌存在明显犁皱而形成的沟槽.同时,通过磨痕的二维轮廓图和三维轮廓图(图 7、图 8所示)也可以看出,球磨2 h、36 h以及48 h的烧结合金在磨损1 h后具有较大的磨痕深度与宽度,表明其磨痕表面存在大量的脱落物.这些成片状脱落的表面材料,在压应力的作用下,不断被碾碎和沿运动方向被横推出摩擦区域,最终导致合金试样磨痕形貌表面出现较多的沟槽.通过对合金试样磨损1 h后的磨损率和磨损量进行精确测定和计算(图 9所示),可以看出不同球磨时间的合金试样的磨损率和磨损量的大小顺序为:48 h>36 h>2 h>24 h.综上,当球磨时间为24 h,合金试样具有最优异的耐磨性能.
分析认为,合金试样摩擦磨损行为的这种变化趋势主要和试样的相对密度、硬度、黏结相的量以及分布等因素有关.当球磨时间为2 h时,合金试样的相对密度较高,但黏结相的量相对较少,仅在局部生成,且分布极不均匀如图 2(a)所示,这大大增加W颗粒与W颗粒直接接触的几率,而W-W之间的界面结合强度很低,导致磨损曲线出现大的波动,致使合金试样的磨痕宽度和深度变大,耐磨性能下降. 当球磨时间延长至36 h和48 h时,合金粉末中的杂质元素和氧含量显著升高,SPS烧结试样无法有效致密化,合金材料相对密度降低,孔隙度较大如图 2(c)、2(d)所示. 大量闭孔隙的存在,一方面影响晶界和位错的滑移,使材料的塑性变形能力下降;另一方面,孔隙处一般为应力集中区,初始裂纹源易于在该处形成,在随后的摩擦磨损试验过程中,这些微裂纹会扩展成大的裂纹,导致磨损曲线产生较大的波动,并加速合金材料的摩擦磨损[23-24]. 当球磨时间为24 h时,既保证合金试样中粘结相的量和均匀分布,又不至于引入大量的杂质元素,引起合金成分的改变. 因此,合金试样显示出优异的耐磨性能.
3 结 论
采用高能球磨和放电等离子体烧结技术制备W-4.9Ni-2.1Fe高比重合金,研究球磨工艺参数对合金摩擦行为的影响,得出如下结论:
1)球磨时间为2 h时,合金粉末中Ni、Fe元素仍以单质的形式存在;当球磨时间大于24 h时,Ni、Fe峰消失,仅存在W元素的衍射峰,Ni(Fe)溶入W晶格中形成W的过饱和固溶体,同时衍射峰出现明显宽化现象;
2)当球磨时间为2 h时,合金相对密度和硬度值较高,但黏结相量相对较少且分布不均匀;当球磨时间为36 h、48 h时,黏结相量较多且分布较为均匀,但过量杂质的引入容易造成合金成分的变化,使合金试样难以烧结致密化;
3)适量的球磨时间(24 h),既可以保证合金中粘结相的含量和均匀分布,又可以使合金成分不产生突变,因此合金拥有最优的耐摩擦磨损性能.
-
表 1 数值模型相关参数
Table 1 Relevant parameters of each layer of soil
表 2 4 d时刻不同停雨间隔时间下的边坡稳定性系数
Table 2 Slope stability coefficient under different persistence on the 4th day
-
[1] 郭钟群, 赵奎, 金解放, 等. 离子型稀土开发面临的问题与绿色提取研究进展[J]. 化工进展, 2019, 38(7): 3425-3433. https://www.cnki.com.cn/Article/CJFDTOTAL-HGJZ201907045.htm [2] 罗嗣海, 袁磊, 王观石, 等. 浸矿对离子型稀土矿强度影响的试验研究[J]. 有色金属科学与工程, 2013, 4(3): 58-61. doi: 10.13264/j.cnki.ysjskx.2013.03.006 [3] 王晓军, 李永欣, 黄广黎, 等. 浸矿过程离子型稀土矿孔隙结构演化规律研究[J]. 中国稀土学报, 2017(4): 528-536. https://www.cnki.com.cn/Article/CJFDTOTAL-XTXB201704013.htm [4] 吴丹瑞, 吴安琪, 何新玥, 等. 赣南地区近58年来极端气候变化趋势分析[J]. 长江科学院院报, 2017, 34(1): 24-29. https://www.cnki.com.cn/Article/CJFDTOTAL-CJKB201701007.htm [5] 邓文建, 饶运章, 王丹, 等. 降雨入渗对离子型稀土边坡的稳定性影响研究[J]. 矿业研究与开发, 2019, 39(11): 6-12. https://www.cnki.com.cn/Article/CJFDTOTAL-KYYK201911002.htm [6] TANG D, JIANG Z M, YUAN T, et al. Stability analysis of soil slope subjected to perched water condition[J]. KSCE Journal of Civil Engineering, 2020, 24(9): 2581-2590.
[7] 王力, 李高, 陈勇, 等. 赣南地区人工切坡降雨致灾机制现场模型试验[J]. 岩土力学, 2021, 42(3): 846-854. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202103027.htm [8] 邱祥, 蒋煌斌, 欧健, 等. 降雨条件下边坡暂态饱和区形成条件与演化特征数值分析[J]. 水利学报, 2020, 51(12): 1525-1535. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB202012010.htm [9] 钟文, 欧阳健, 黄震, 等. 不同降雨条件下稀土尾矿堆入渗与破坏规律试验研究[J]. 中南大学学报(自然科学版), 2021, 52(11): 4099-4109. https://www.cnki.com.cn/Article/CJFDTOTAL-ZNGD202111028.htm [10] ZHANG Y, LU H F. Unsaturated seepage and stability of double-layer slope under rainfall infiltration[J]. E3S Web of Conferences, 2021, 248(7) : 03024.
[11] 张蔚榛. 地下水与土壤水动力学[M]. 北京: 中国水利水电出版社, 1996: 1-5. [12] 戚蓝, 王守甲, 傅长锋. 短历时强降雨下边坡滑移特征时空分布规律的研究[J]. 水资源与水工程学报, 2018, 29(6): 190-195. https://www.cnki.com.cn/Article/CJFDTOTAL-XBSZ201806029.htm [13] VAN GENUCHTEN M T. A closed form equation for prediction the hydraulic conductivity of unsaturated soils[J]. Soil Science Society of America Journal, 1980, 44(5): 892-898.
[14] 万思豪, 张继勋, 张景风. 不同降雨条件下裂隙性土质边坡渗流特性分析[J]. 三峡大学学报(自然科学版), 2019, 41(1): 33-37. https://www.cnki.com.cn/Article/CJFDTOTAL-WHYC201901009.htm [15] 彭珂, 彭红霞, 梁峰, 等. 赣州市地质灾害分布特征及孕灾环境分析[J]. 安全与环境工程, 2017, 24(1): 33-39. https://www.cnki.com.cn/Article/CJFDTOTAL-KTAQ201701006.htm [16] 赵奎, 赵康, 徐贻赣, 等. 赣州境内地质灾害实例分析[J]. 采矿技术, 2005, 5(2): 36-38. https://www.cnki.com.cn/Article/CJFDTOTAL-SJCK200502019.htm [17] 荣冠, 张伟, 周创兵. 降雨入渗条件下边坡岩体饱和非饱和渗流计算[J]. 岩土力学, 2005(10): 24-29. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200510004.htm [18] LI S H, JIANG Z L, QUE Y, et al. Water field distribution characteristics under slope runoff and seepage coupled effect based on the finite element method[J]. Water, 2021, 13(24): 3569-3569.
[19] 年庚乾, 陈忠辉, 张凌凡, 等. 边坡降雨入渗问题中两种边界条件的处理及应用[J]. 岩土力学, 2020, 41(12): 4105-4115. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202012031.htm [20] 邬长福, 姚贵佳, 陈亮, 等. 基于Geo-Studio的离子型稀土矿边坡稳定性分析[J]. 中国地质灾害与防治学报, 2016, 27(2): 72-77. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGDH201602011.htm [21] 何书, 陈飞, 刘强, 等. 稀土矿层厚度对浸矿边坡渗流及稳定性的影响[J]. 防灾减灾工程学报, 2018, 38(3): 528-534. https://www.cnki.com.cn/Article/CJFDTOTAL-DZXK201803017.htm [22] 张昊, 余巍伟, 林杭, 等. 不同安全系数对应的边坡滑动面位置分析[J]. 岩土力学, 2012, 33(2): 449-452. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201202023.htm [23] 高冯, 李小军, 迟明杰. 基于有限元强度折减法的单双面边坡稳定性分析[J]. 工程地质学报, 2020, 28(3): 650-657. https://www.cnki.com.cn/Article/CJFDTOTAL-GCDZ202003022.htm [24] 栾茂田, 武亚军, 年廷凯. 强度折减有限元法中边坡失稳的塑性区判据及其应用[J]. 防灾减灾工程学报, 2003, 23(3): 1-8. https://www.cnki.com.cn/Article/CJFDTOTAL-DZXK200303001.htm