2.吉林大学 机械科学与工程学院, 吉林 长春 130025
2.School of Mechanical Science and Engineering, Jilin University, Changchun 130025, China
破碎过程模型的建立是矿物处理工业的核心.随着对高精度、高复杂度模型的需求,一些计算机软件(如EDEM)被广泛采纳并应用于各工业领域,取得了很好的效果[1-4].由于破碎过程模型的建立主要依赖岩石本构参数,因此如何快捷、准确地获得岩石微观力学参数是建立破碎过程模型的首要任务.目前研究者多从黏结颗粒模型(BPM)着手[5-8],采用夹逼处理的方法,以岩石宏观力学特性为指标进行微观参数标定[9],但该方法工作量大、耗时长.因此,合适的BPM建模方法对提高离散元破碎仿真的效率具有重大意义.
本文从破碎原理出发,基于EDEM对岩石进行单轴压缩数值实验,找到影响岩石宏观力学性能的关键因素及影响规律;进而,通过全析因实验设计得到影响岩石宏观力学性能的关键主效应和交互效应,建立BPM宏观力学性能参数理论方程;最后,通过回归方法对该方程的系数进行求解,并进行离散元仿真和物理实验验证.
1 BPM模型的预测 1.1 影响岩石宏观力学性能关键参数分析经虚拟实验参数匹配,并保证岩石弹性变形阶段显著,选择五莲花柱状花岗岩作为标准试样,直径50 mm,高径比为2:1,构建BPM,微观颗粒半径为1.5 mm,平行键法向刚度和强度分别为1.1×1012 Pa和1.8×108 Pa·m-1,微观颗粒泊松比、剪切模量和密度分别为0.2,43 GPa,2 780 kg·m-3,颗粒之间摩擦系数为0.3.根据Refahi的研究[10],岩石经N个时间步后,当外部累加载荷和大于岩石强度极限时,岩石发生破碎.又因为岩石内部受力情况与摩擦系数和杨氏模量相关,因此以岩石材料内摩擦系数、剪切模量、平行键法向强度和刚度以及微观颗粒数量为因子进行多次单因素数值实验,结果如图 1所示,可看出这些微观参数对岩石抗压强度和宏观刚度有显著影响,且存在一定规律.
在EDEM环境下进行五因子两水平全析因实验[11],对BPM岩石宏观性能进行预测,并基于此计算出平行键参数.设置4个中心点,36组实验数据如表 1所示.所得响应结果与实验运行序关系如图 2a和2b所示,其中0代表中心数据点,1代表其他实验数据点.从散点图可看出,趋势正常,并未出现连续上升或下降趋势,且岩石宏观刚度数据没有明显偏倚,因此可视为正态分布;而抗压强度响应指标明显向低值偏倚,但样本量大于30,可认为样本均值趋于正态.从散点图中还可看出,中心点所对响应值稳定,说明实验环境趋于一致.
为考虑各因素效应对响应影响,并采用序贯实验对结果进行分析,对响应可采用式(1)进行预估:
(1) |
其中:y为响应变量;xi为主效应;xi·xj为交互效应.对各效应对响应的显著性进行判断,并提取关键主效应和交互效应.设置T检验显著性水平α=0.05,构建帕累托效应图,如图 3所示.根据各效应对响应影响的显著性,忽略量纲影响,可建立岩石宏观刚度和抗压强度的解析方程,如式(2)和式(3)所示.通过计算,可得解析式的各项系数,如表 2所示.
(2) |
(3) |
对宏观力学参数模型的有效性、失拟性和曲率进行假设检验,其中主效应、弯曲和失拟所对应的P值如表 3所示.
从表 3中可以看出,该线性回归方程总的来说是有效的.虽忽略了很多交互效应,但模型并无发生明显失拟和弯曲.
为判断回归方程优劣,还需对其拟合优度和残差进行研究.回归分析中的离差总平方和SST由回归平方和SSR和残差平方和SSE两部分构成,如式(4)所示;其拟合优度R2及扣除回归方程中包含项数影响的拟合优度Radj2分别如式(5)和式 (6)所示.其中n为观测值总个数,将调整前后拟合优度汇总成表 4.
(4) |
(5) |
(6) |
从表 4可看出,模型性态良好.模型残差与拟合值的关系如图 4所示.从图中可看出,去除异常点后,抗压强度和宏观刚度的误差率均小于10%.
在EDEM环境中,构造3个验证岩石试样,进行单轴压缩数值实验,按照表 5设置变量值,岩石抗压强度和宏观刚度预测结果和仿真结果如表 6所示,从表中可见模型有很好的预测精度.
为进一步验证模型有效性,还需进行物理实验验证,选取五莲花花岗岩标准试样进行单轴压缩,如图 5a所示.以岩石抗压强度及宏观刚度为验证指标,在EDEM环境下进行仿真,如图 5b所示.按照预测模型设置平行键刚度和强度,比较物理实验和数值实验的应力应变图,如图 6所示.从图中可看出,物理实验和数值实验曲线图接近,但略有差异,符合统计学规律.
1) 从岩石破碎原理出发,进行大量岩石压缩数值分析.结果表明对岩石宏观力学性能的影响因素为岩石内摩擦系数、剪切模量、平行键刚度、平行键强度和微观颗粒个数.
2) 基于全析因实验设计方法,采用多项式回归模型建立岩石宏观刚度和强度与关键BPM参数函数关系,并对该函数的方差、拟合优度和残差进行分析.
3) 以岩石宏观强度和刚度为指标,采用离散元方法和物理实验对模型进行验证,其误差率均小于10%.因此该预测模型可根据物理实验所得岩石宏观力学性能准确、快捷地得出平行键强度和刚度,提高了破碎仿真效率.
[1] | Weerasekara N S, Powell M S, Cleary P W, et al. The contribution of DEM to the science of comminution[J]. Powder Technology , 2013, 248 : 3–24. DOI:10.1016/j.powtec.2013.05.032 (0) |
[2] | Lichter J, Lim K, Potapov A, et al. New development in cone crusher performance optimization[J]. Minerals Engineering , 2009, 22 (7) : 613–617. (0) |
[3] | Ohannes J.Cone crusher modeling and simulation-development of a virtual rock crushing environment based on the discrete element method with industrial scale experiments for validation [D].Nybro:Chalmers University of Technology ,2012. (0) |
[4] |
蒋明镜, 陈贺, 刘芳.
岩石微观胶结模型及离散元数值仿真方法初探[J]. 岩石力学与工程学报 , 2013, 32 (1) : 15–23.
( Jiang Ming-jing, Chen He, Liu Fang. A microscopic bond model for rock and preliminary study of numerical simulation method by distinct element method[J]. Chinese Journal of Rock Mechanics and Engineering , 2013, 32 (1) : 15–23. ) (0) |
[5] | Potyondya D O, Cundall P A. A bonded-particle model for rock[J]. International Journal of Rock Mechanics & Mining Sciences , 2004, 41 (8) : 1329–1364. (0) |
[6] | Ding X B, Zhang L Y. A new contact model to improve the simulated ratio of unconfined compressive strength to tensile strength in bonded particle models[J]. International Journal of Rock Mechanics & Mining Sciences , 2014, 69 (3) : 111–119. (0) |
[7] | Manouchehrian A, Sharifzadeh M, Marji M F, et al. A bonded particle model for analysis of the flaw orientation effect on crack propagation mechanism in brittle materials under compression[J]. Archives of Civil and Mechanical Engineering , 2014, 14 (1) : 40–52. DOI:10.1016/j.acme.2013.05.008 (0) |
[8] | Obermayr M, Dressler K, Vrettos C, et al. A bonded-particle model for cemented sand[J]. Computers and Geotechnics , 2013, 49 : 299–313. DOI:10.1016/j.compgeo.2012.09.001 (0) |
[9] | Quist J.Cone crusher modelling and simulation[D].Nybro:Chalmers University of Technology,2012. (0) |
[10] | Refahi A, Aghazadeh M J, Rezaib B. Discrete element modeling for predicting breakage behavior and fracture energy of a single particle in a jaw crusher[J]. International Journal of Mineral Processing , 2010, 94 (2) : 83–91. (0) |
[11] | Hassan J, Mohd H I, Amirreza S. Evaluation of significant manufacturing parameters in lost foam casting of thin-wall Al-Si-Cu alloy using full factorial design of experiment[J]. Ransactions of Nonferrous Metals Society of China , 2013, 23 (10) : 2843–2851. DOI:10.1016/S1003-6326(13)62805-8 (0) |