不规则铜钼矿颗粒断裂强度分布实验研究.pdf

不规则铜钼矿颗粒断裂强度分布实验研究周强1,2,3，汪轶凡1，肖庆飞1，刘向阳1，邵云丰1，黄守向1，王庆凯2,3，邹海11.昆明理工大学 国土资源工程学院，云南 昆明 650093；2.矿冶过程自动控制技术国家重点实验室，北京 102628；3.矿冶过程自动控制技术北京重点实验室，北京 102628中图分类号：TD921+.2；TD453文献标识码：A文章编号：10010076（2023）04003310DOI：10.13779/ki.issn1001-0076.2023.04.004摘要不规则矿石颗粒的断裂是矿物加工过程中常见的现象，其强度分布决定了矿石破碎特性。为了定量分析不规则颗粒的强度分布，通过对 5 种不同粒级的铜钼矿颗粒进行准静态单轴压缩实验，确定破碎过程中的最大破碎力和断裂能。选取 3 种常见的统计学模型分别对不同定义下的（最大破碎力、断裂应力、断裂能和断裂比能）颗粒强度进行拟合，并研究了其与颗粒尺寸和材料特性之间的定量关系。试验结果表明：Weibull 模型较 Lognormal 和 Lmogintic 两种模型更适合描述铜钼矿颗粒的强度分布，其模型中强度分布的离散程度 D 只与材料特性有关，与颗粒尺寸呈弱函数关系；F63.20、E63.20与颗粒尺寸呈正比关系，而 63.20和 Em63.20随着颗粒尺寸的增大呈幂函数规律减小；不同定义下的颗粒强度之间（最大破碎力断裂能和断裂应力断裂比能）的关系都只与材料特性有关，而与颗粒尺寸无关，在双对数坐标系下斜率分别为 1.49 与 0.67。关键词不规则颗粒；颗粒强度分布；选择函数；Weibull 模型 引言不规则矿石颗粒的断裂是碎矿与磨矿过程中的常见现象，具有重要的理论研究与工业应用价值。随着时代快速发展，碎磨设备趋于大型化和复杂化，但破碎过程中的能量利用效率依然很低。许多学者对单颗粒破碎、颗粒床破碎、冲击破碎和球磨 4 种破碎方式的能量利用效率进行研究1，实验结果表明：以单晶断裂的比表面能作为基准来评估不同破碎方式的能量利用效率，其中不规则单颗粒破碎 1.0%，颗粒床破碎 0.7%，冲击破碎 0.5%，球磨 0.15%。此外，影响颗粒破碎性的因素有很多，主要取决于颗粒的物理性质：粒径、级配、形状、孔隙比、断裂强度等，导致颗粒的破碎过程是复杂难测的，目前还未有数学模型能描述所有破碎过程。因此如何提高破碎过程的能量利用效率与构建统一的破碎数学模型是目前亟待解决的问题1-2。构建描述矿石断裂过程的破碎数学模型需要确定断裂概率和断裂函数3。近年来，国内外学者对断裂概率与断裂函数开展了大量的实验研究，取得了重要的研究成果4-7：根据实验设备的不同，可以选择对应的实验工况来确定断裂概率和断裂函数。例如，对于颚式破碎机、圆锥破碎机和高压辊磨机等准静态加载设备，可以通过压力实验机来确定颗粒的断裂概率与断裂函数8；对于冲击式破碎机、球磨机等加载速率高的设备，可以通过落重冲击实验、双摆锤实验或UFLC 实验系统来确定断裂概率与断裂函数。大量实验结果表明：即使矿石颗粒尺寸一致、性质相同，其断裂强度与破碎产品粒度分布也会呈现出显著的离散性，主要原因是颗粒内部的原生缺陷在大小、分布、成核、拓展和汇聚上存在很大的随机性9-12。颗粒强度的定义方法有很多种，最直观的方式是以最大破碎力来定义颗粒强度13，相关数据可以直接通过实验测量得到，但通过不同的处理方式，我们也可以通过断裂应力10、断裂能14和断裂比能15来定义强度。利用统计学方法对不规则颗粒强度进行分析时，为保证所得数值的可信度，必须有足够多的样本数据。国内学者通常选用双参数或三参数的 Weibull分布应用于描述颗粒的强度分布以及破碎产品粒度 收稿日期：2023 06 01基金项目：国家自然科学基金地区科学基金项目（51964044）；矿冶过程自动控制技术国家重点实验室开放基金项目（BGRIMMKZSKL20221）；云南省科技厅基础研究项目（202301AT070392）；云南省教育厅科学研究基金项目（2023J0125）作者简介：周强（1991），男，山东枣庄人，博士，讲师，硕士生导师，主要从事碎磨理论及碎磨过程数值模拟研究，E-mail：。通信作者：肖庆飞（1980），男，云南昆明人，博士，教授，博士生导师，主要从事碎磨理论与工艺的研究，E-mail：。第 4 期矿产保护与利用No.42023 年 8 月Conservation and Utilization of Mineral ResourcesAug.2023分布7,9-10。因此分析不规则颗粒的断裂强度时，可以利用该统计函数。除此之外，本文还采用了其他统计函数去拟合颗粒的强度分布，如Logistic 分布11、Lognormal分布12等，通过比较选择最优的统计函数来描述不规则铜钼矿的强度分布。选用最优的统计函数对不规则颗粒的强度分布进行拟合，并对选择函数中的参数进行分析，建立了不规则铜钼矿的选择函数与颗粒尺寸关系的定量关系式。对于颗粒强度的定义有很多种，但很少有学者深入研究对它们之间的关系，因此本文还对不同的强度定义进行对比，分析并建立了它们之间的函数关系式。通过研究分析不规则颗粒断裂概率与颗粒尺寸以及不同强度定义之间的关系，确定了选择函数，为破碎模型的建立提供依据。1实验材料与实验方案1.1实验材料实验研究对象为华泰龙不规则铜钼矿石颗粒，铜钼矿在常温常压下的密度为 2.93103 kg/m3，平均抗压强度为 66.32 MPa，弹性模量平均为 4.35104 MPa，泊松比为 0.16，整体呈灰白色，含铜黄色金属光泽颗粒。本次实验选用五个粒级进行强度分布的研究，分别为0.9 2 mm、2 5 mm、5 8 mm、8 10 mm、1012 mm，每个粒级选择 3035 个不规则颗粒进行实验。1.2实验方案本次实验采用 YCDW5E 型电子伺服实验系统，该系统主要测量颗粒的最大破碎力以及断裂能，可提供最大轴向压力 5 000 N，实验力测量范围 1005 000 N，实验力示值精度为1%，丝杠最大移动速度 500 mm/min，最大压缩空间 800 mm，实验过程中砧板压缩速度1 mm/min。实验装置 YCDW5E 型电子伺服实验系统如图 1 所示，实验材料在加载前被放在上下压头之间，为减少摩擦，可在压头与试样之间滴几滴润滑油。所有实验通过位移控制，加载速度低，惯性效应可不计，固定上压头，下压头缓缓上升压缩试样。当实验材料被完全破碎，或者检测到轴向力发生急速下降现象时，装置停止加载。该实验系统的加载装置承载能力较好，因此可忽略装置自身变形的影响。电脑在加载过程中记录实验数据，后续通过软件绘出力位移曲线，该曲线将用来整理试样的相关数据，例如最大破碎力以及对应的断裂能等。本次实验每个粒级的实验次数皆在 30 次以上，共计 150 余组实验。1.3颗粒强度计算方法不同定义下固体颗粒强度可由最大破碎力、断裂应力、断裂能、断裂比能等来表示，其计算方法如下。1.3.1最大破碎力最大破碎力是指物料发生灾难性破坏时所对应的力，反映在接触力位移曲线上，如图 2 所示。Fc为最大破碎力，c为最大破碎力对应的位移。由于不规则矿石颗粒的不规则性，导致加载点附近产生应力集中，产生局部破碎，因此需要满足以下两个条件则实验结束：（1）不规则颗粒与上下压头的接触点之间出现贯穿性裂纹；（2）需要满足图 2（a）中出现的“悬崖”（即急剧下降线），以保证不规则颗粒为宏观的张拉性破坏。在实验中发现，由于矿石颗粒的不规则性，其接触力位移曲线的类别也存在很大的差异。本文分析了 150 条不规则铜钼矿单颗粒破碎过程的接触力与位移曲线特征，通过总结归纳，可将其特征定义成两大类，如图 2（a）与（b）所示。第一类曲线：在加载初期会出现一段平坦区，此时处于重心稳定期，实验力变化不大，在这个阶段主要是对不规则颗粒位置的调整。随后进入弹性变形期，实验装置不断加大对颗粒的压缩，并在实验过程中会发现颗粒表面与压头接触部位的棱角发生断裂，为局部断裂，如图 2（a）箭头 1 所指位置，此时处于局部棱角剥离期。随着接触力继续增大，当实验力达到Fc时，颗粒发生完全断裂，反映在力位移曲线上则是载荷的突然跌落，如图 2（a）箭头 2 所指位置，此时处于完全断裂期。箭头 3 所指位置反映的是加载装置对颗粒的二次破碎，试样处于二次破碎期，此时试样破碎所需的力要小于 Fc。这种曲线是较为标准的，我们可以很容易定义其最大破碎力。图 1试验装置Fig.1 Testing device 34 矿产保护与利用2023 年第二类曲线：整个破碎实验过程会形成多个荷载峰值，且随着轴向位移的增加，荷载峰值会逐渐水平，无减弱的趋势。出现此类形态的接触力位移曲线时，需辅助高速摄影来确定不规则颗粒的最大破碎力，根据接触点之间出现宏观裂纹来确定。实验中发现，初始形状扁平，形状极不规则的颗粒会出现该类别的接触力位移曲线。1.3.2断裂应力断裂应力被广泛用来表征颗粒的强度。当试样是标准试件时，其与外界作用力接触的面积容易测量，进而确定断裂应力；但在实际应用中，尤其对于不规则颗粒，接触面积很难准确测量，因此，有学者提出很多经验公式10，通过经验修正系数减小误差来计算断裂应力。本文采用式（1）计算不规则颗粒断裂应力 c16：c=4Fcda2（1）式中：Fc为颗粒最大破碎力，N；da为颗粒的几何平均径，mm；取 3.141 592 6。1.3.3断裂能如图 3 所示，在接触力位移曲线上，颗粒的形变程度与其承受的破碎力一一对应，而断裂点前的力位移曲线的线下面积在数值上等同于断裂能 Ec，可以通过式（2）得到14：Ec=wc0Fd（2）式中：F 为破碎力，N；c为最大破碎力对应的位移，m。1.3.4断裂比能断裂比能 Em为破坏单位质量或体积所需要的能量，假定颗粒破坏前后尺寸近似为常数，可以根据式（3）计算15。Em=6Ecda3（3）式中：Ec为颗粒断裂能，J；为颗粒密度，kg/m3；da为颗粒的几何平均径，m；取 3.141 592 6。1.3.5断裂概率计算方法当讨论不规则颗粒的断裂强度分布时，必须对足够数量的样本进行实验。对于有限数量的样本，颗粒的断裂概率可以通过概率估算因子 P 进行计算17。P=in+1（4）式中：n 为实验样品的总数；i 为对所有试样的强度数据进行升序排列后某一个颗粒强度在序列中的排名。2结果与讨论2.1不同定义下颗粒的强度分布通过公式（4）计算出的断裂概率，对不同定义下的强度数据进行整理分析，得到不同定义下颗粒的强度分布，选取 58 mm 铜钼矿颗粒的强度分布图进行举例，如图 4 所示。由图 4 所示，强度数据的分布集中中间及较低的部分，较大的强度数据稀疏。以图 4（a）为例，0600 N之间就有 33 个数据点，而 600800 N 之间只有 2 个 0.00.51.01.52.02.50100200300400500600试验力/N位移/mm断裂点(Fc,c)Fcc123平坦区(a)0.00.51.01.52.0020406080100120试验力/N位移/mm断裂点(Fc,c)cFc1连续峰值区平坦区(b)（a）第一类的接触力位移曲线（b）第二类的接触力位移曲线图 2两类铜钼矿单颗粒压缩破碎接触力位移曲线Fig.2 Two kinds of forcedisplacement curves of single coppermolybdenum particle crushing 0.00.51.01.52.00100200300400500600700试验力/N位移/mmcFc断裂点(Fc,c)图 3断裂能定义图Fig.3 Fracture energy defines the diagram第 4 期周强，等：不规则铜钼矿颗粒断裂强度分布实验研究 35 数据点。2.2颗粒强度分布函数分别选择 Logistic 模型（公式（5）、Weibull 模型（公式（6）、Lognormal 模型（公式（7）去拟合不同定义下的强度数据与断裂概率之间的函数关系，为了研究结果的准确性，先从中找出不规则颗粒的最优统计模型。P=111+(Fca)b（5）P=1exp(Fca)b（6）P=121+erf(ln(Fc)a2b)（7）式中：P 为颗粒的断裂概率；a，b 为模型参数。当用其他方式来定义颗粒的强度时，只需替换 Fc即可。以最大破碎力为例，比较 3 种统计模型与实验数据的拟合精度，拟合结果如图 5 所示。由图 5 可知，3 种统计模型均与实验数据吻合度较好，而确定最优模型最直接的方法就是通过相关系数 R2比较拟合精度。同样，利用统计模型对其他 3 种不同定义强度的数据进行拟合，最后得到不同统计模型下的 R2，见表 1。由表 1 可知，3 个模型均可以用来描述不规则颗粒的强度分布，Logistic 模型和 Lognormal 模型的拟合相关系数接近，Logistic 模型的精度较 Lognormal 模型来说要略大，但两者均小于 Weibull 模型拟合的相关系数。Weibull 模型的数学形式简单，且模型中的参数有明确的物理意义，其中参数 a 为断裂概率为 63.20%对应的强度，参数 b 为强度分布的离散程度。再将公式（6）转化为更简单的形式，得到的函数分别以最大破碎力、断裂应力、断裂能和断裂比能表示，如公式（8）（11）：Fc=F63.20ln(11PF)DF（8）0200400600800100012000.00.20.40.60.81.0 铜钼矿颗粒58 mm断裂概率最大破碎力/N(a)05101520253035400.00.20.40.60.81.0 铜钼矿颗粒58 mm断裂概率断裂应力/MPa(b)0.00.20.40.60.81.00.00.20.40.60.81.0断裂概率断裂能/J 铜钼矿颗粒58 mm(c)0501001502002503000.00.20.40.60.81.0 铜钼矿颗粒58 mm断裂概率断裂比能/(Jkg-1)(d)（a）最大破碎力断裂概率（b）断裂应力断裂概率（c）断裂能断裂概率（d）断裂比能断裂概率图 4不同定义下颗粒的强度分布Fig.4 Strength distribution of particles under different definitions 01002003004005006007008009000.00.20.40.60.81.0 断裂概率 Logistic 模型 Lognormal 模型 Weibull 模型断裂概率最大破碎力/N图 5试验数据与 3 个统计模型拟合的对比Fig.5 Comparison of experimental data with three statisticalmodel fits 36 矿产保护与利用2023 年c=63.20ln(11P)D（9）Ec=E63.20ln(11PE)DE（10）Em=Em63.20ln(11PEm)DEm（11）综上所述，Fc为最大破碎力；c为断裂应力；Ec为断裂能；Em为断裂比能；F63.20为断裂概率 63.20%对应的破碎力；63.20为断裂概率 63.20%对应的断裂应力；E63.20为断裂概率 63.20%对应的断裂能；Em63.20为断裂概率 63.20%对应的断裂比能；PF、P、PE和 PEm分别为以最大破碎力、断裂应力、断裂能和断裂比能定义的断裂概率；DF、D、DE和 DEm分别为以最大破碎力、断裂应力、断裂能和断裂比能定义的分布的离散程度，该值越大，离散程度越小即数据较集中。对三个模型再进行整体的比较，比较三个模型对不同尺寸下不同强度分布定义的拟合精度。图 6 是用 Weibull 模型拟合铜钼矿颗粒在不同尺寸下不同强度分布定义的拟合结果，可以看出拟合曲线与实验数据吻合度均较好。图 7 是用 Logistic 模型拟合铜钼矿颗粒在不同尺寸下不同强度分布定义的拟合结果，可以看出拟合精度低于 Weibull 模型。图 8 是用 Lognormal 模型拟合铜钼矿颗粒在不同尺寸下不同强度分布定义的拟合结果，可以看出拟合精度低于 Weibull 模型和 Logistic 模型。将总实验数据拟合的相关系数绘制相关系数误差分析图，如图 9。通过综合比较，最后选择 Weibull 模型来描述不规则颗粒的强度分布。0500100015002000250030000.00.20.40.60.81.0 0.92 mm 25 mm 58 mm 810 mm 1012 mm Weibull 模型 R2=0.9748断裂概率最大破碎力/N(a)0501001500.00.20.40.60.81.0 0.92 mm 25 mm 58 mm 810 mm 1012 mm Weibull 模型 R2=0.9748断裂概率断裂应力/MPa(b)01230.00.20.40.60.81.0 0.92 mm 25 mm 58 mm 810 mm 1012 mm Weibull 模型 R2=0.9564断裂概率断裂能/J(c)05001000150020000.00.20.40.60.81.0 0.92 mm 25 mm 58 mm 810 mm 1012 mm Weibull 模型 R2=0.9564断裂概率断裂比能/(Jkg-1)(d)（a）最大破碎力断裂概率（b）断裂应力断裂概率（c）断裂能断裂概率（d）断裂比能断裂概率图 6Weibull 模型对不同尺寸下的试验数据的拟合Fig.6 Weibull model fitting test data of materials with different sizes 表 1 部分试验数据拟合的相关系数Table 1 Correlation coefficient fitted to some experimental data类型模型R2最大破碎力Logistic0.984Lognormal0.981Weibull0.989断裂应力Logistic0.984Lognormal0.981Weibull0.989断裂能Logistic0.979Lognormal0.981Weibull0.990断裂比能Logistic0.979Lognormal0.981Weibull0.990第 4 期周强，等：不规则铜钼矿颗粒断裂强度分布实验研究 37 0500100015002000250030000.00.20.40.60.81.0 0.92 mm 25 mm 58 mm 810 mm 1012 mm Logistic 模型 R2=0.9672断裂概率最大破碎力/N(a)0501001500.00.20.40.60.81.0 0.92 mm 25 mm 58 mm 810 mm 1012 mm Logistic 模型 R2=0.9672断裂概率断裂应力/MPa(b)012340.00.20.40.60.81.0 0.92 mm 25 mm 58 mm 810 mm 1012 mm Logistic 模型 R2=0.9514断裂概率断裂能/J(c)05001000150020000.00.20.40.60.81.0 0.92 mm 25 mm 58 mm 810 mm 1012 mm Logistic 模型 R2=0.9514断裂概率断裂比能/(Jkg-1)(d)（a）最大破碎力断裂概率（b）断裂应力断裂概率（c）断裂能断裂概率（d）断裂比能断裂概率图 7Logistic 模型对不同尺寸下的材料的试验数据的拟合Fig.7 Logistic model fitting test data of materials with different sizes 0500100015002000250030000.00.20.40.60.81.0 0.92 mm 25 mm 58 mm 810 mm 1012 mm Lognormal 模型 R2=0.9662断裂概率最大破碎力/N(a)0501001500.00.20.40.60.81.0 0.92 mm 25 mm 58 mm 810 mm 1012 mm Lognormal 模型 R2=0.9662断裂概率断裂应力/MPa(b)012340.00.20.40.60.81.0 0.92 mm 25 mm 58 mm 810 mm 1012 mm Lognormal 模型 R2=0.9496断裂概率断裂能/J(c)05001000150020000.00.20.40.60.81.0 0.92 mm 25 mm 58 mm 810 mm 1012 mm Lognormal 模型 R2=0.9496断裂概率断裂比能/(Jkg-1)(d)（a）最大破碎力断裂概率（b）断裂应力断裂概率（c）断裂能断裂概率（d）断裂比能断裂概率图 8Lognormal 模型对不同尺寸下的材料的试验数据的拟合Fig.8 Lognormal model fitting test data of materials with different sizes 38 矿产保护与利用2023 年2.3颗粒强度分布函数的扩展利用式（6）计算出来的不规则铜钼矿的 Weibull分布模型中的参数 b 展示在表 2 中。DF、D、DE和DEm分别为以破碎力、断裂强度、断裂能和断裂比能定义下的分布的离散程度，该值越大，离散程度越小即数据较集中。由表 2 可知，DF、D、DE和 DEm与颗粒的尺寸没有表现出显著的关系，这也进一步证实了之前学者提出 Weibull 分布模量是个常数的假设17-18。取各尺寸的平均值 DF、DE，作为函数参数初始值并固定，重新拟合不同定义下的强度分布数据，计算误差范围，若误差范围允许，则可简化 Weibull 模型，见表 3。拟合结果表明，当 DF、DE不变时，与之前的拟合结果比较，其相关系数较低，但误差仍在可接受的范围内，验证了 DF、DE只与材料特性有关，与颗粒尺寸存在弱函数关系。因此可固定 D 值，从而简化 Weibull模型。表 3 DF、DE与颗粒尺寸的关系Table 3 The relationship of DF、DE and particle size粒级/mm最大破碎力/N断裂能/JDFF50R2DFF50R2DEE50R2DEE50R20.921.66109.140.9451.72108.350.9441.5127.790.8811.2330.080.873251.62349.060.9901.72347.360.9891.25145.300.9891.23145.760.989581.95557.380.9641.72559.270.9601.18430.100.9891.23429.210.9888101.52902.900.9901.72893.580.9850.95911.010.9431.23928.010.92610121.871 317.510.9851.721 306.470.9831.301 366.010.9801.231 373.010.980 本文也对 Weibull 模型中的 F63.20、63.20、E63.20和Em63.20进行研究，发现 F63.20、E63.20与颗粒的尺寸呈线性关系，如图 10 所示。F63.20、E63.20与颗粒尺寸呈正比关系，颗粒的尺寸越大，F63.20、E63.20也越大，并由此得到式（12）以及（13）。F63.20=AF+BFd（12）E63.20=AE+BEd（13）式中：AF、BF、AE、BE均为只与材料特性有关的特性参数。同时对 63.20、Em63.20进行研究，对其进行 Tavares模型的拟合，发现两者的相关系数分别达到 0.993，0.983，拟合精度均较高，如图 11、12，并由此得到公式（14）以及（15）：63.20=A1+(dB)C（14）Em63.20=AEm1+(dBEm)CEm（15）式中：A、B、C、AEm、BEm、CEm均为只与材料特性有 0246810120200400600800100012001400 F63.20=119.56d-104.85 R2=0.964 E63.20=136.40d-307.94 R2=0.942粒度/mmF63.20/N02004006008001000120014001600E63.20/J图 10铜钼矿颗粒的 F63.20和 E63.20与颗粒尺寸的关系Fig.10 Relations between F63.20 and E63.20 of coppermolybdenumore particles and particle size 断裂比能断裂能最大破碎力 断裂应力0.900.920.940.960.981.00R2颗粒不同定义下的强度分布 Logistic Lognormal Weibull图 9相关系数误差分析图Fig.9 Correlation coefficient error analysis plot 表 2 不同粒级的铜钼矿颗粒强度分布拟合离散程度汇总Table 2 Summary of the degree to which the intensitydistribution fits discretely for coppermolybdenum ore particlesof different sizes粒级/mmDFDDEDEm0.921.661.661.511.51251.621.621.251.25581.951.951.181.188101.521.520.950.9510121.871.871.301.30平均值1.721.721.241.24第 4 期周强，等：不规则铜钼矿颗粒断裂强度分布实验研究 39 关的特性参数。2.4最大破碎力与断裂能的关系在 4 种定义颗粒强度的方式中，最大破碎力可以通过实验系统直接测量，而断裂能则可以经过力位移曲线积分可得。将测量的最大破碎力和与其对应的断裂能绘制在双对数坐标系中，发现不同尺寸的实验数据可用一条直线拟合，如图 13 所示。通过拟合可得铜钼矿颗粒的斜率为 1.49。由图 13 可得断裂能与最大破碎力之间的关系只与材料特性有关，而与颗粒尺寸无关。2.5断裂应力与断裂比能的关系TAVARES10研究发现，对于脆性材料的荷载与位移的曲线可以近似用 Hertzian 方程拟合，只在加载的初始阶段由于颗粒的错动会引起一些误差，而在弹性阶段几乎重合。并通过理论推导断裂应力与断裂比能存在定量关系,即c=Ke2/5Em3/5（16）式中，Ke为物料特性参数；为颗粒的形状系数。本文也对断裂应力与断裂比能之间的关系进行研究，将断裂应力和与其对应的断裂比能绘制在双对数坐标系中，不同尺寸的实验数据同样均可用一条直线来拟合，进一步说明不同定义强度分布之间的关系仅与材料特性有关，而与颗粒尺寸无关。如图 14 所示，直线的斜率为 0.67。1101001000110100断裂应力/MPa断裂比能/(Jkg-1)0.92 mm 25 mm 58 mm 810 mm 1012 mm y=0.67x-0.12 R2=0.964图 14铜钼矿颗粒的断裂应力与断裂比能的关系Fig.14 Relations between crushing stress and breakage specificenergy of Cupromolybdenum ore particles 3结论（1）对比 Logistic 模型、Lognormal 模型和 Weibull模型，前 2 种模型拟合精度均低于 Weibull 模型，再考虑到三种模型数学形式的复杂程度以及模型参数的物理意义，选择 Weibull 模型来描述颗粒的强度分布，其中参数 a 为断裂概率为 63.20%对应的强度，参数 b为强度分布的离散程度。（2）Weibull 模型中 D 与材料特性有关，与颗粒尺寸呈弱函数关系，因此在使用 Weibull 模型拟合时，可固定 D，从而简化模型。模型中的 F63.20、E63.20随颗粒尺寸增大呈线性规律增大，而 63.20和 Em63.20随着颗粒尺寸的增大而呈幂函数规律减小。（3）将最大破碎力断裂能以及断裂应力断裂比能绘制在双对数坐标系中，发现即使颗粒的尺寸不同，0246810120102030405060708063.20/MPa粒度/mm 63.20 y=0.0331+(d/23052.26)-0.79 R2=0.993图 11铜钼矿颗粒的 63.20与颗粒尺寸的关系Fig.11 Relations between 63.20 of coppermolybdenum oreparticles and particle size 0246810120100200300400500600700800900 Em63.20 y=37.141+(d/13.88)-1.304 R2=0.983Em63.20/(Jkg-1)粒度/mm图 12铜钼矿颗粒的 Em63.20与颗粒尺寸的关系Fig.12 Relations between Em63.20 of coppermolybdenum oreparticles and particle size 101001000101001000断裂能/J最大破碎力/N 0.92 mm 25 mm 58 mm 810 mm 1012 mm y=1.49x-1.53 R2=0.973图 13铜钼矿颗粒的断裂能与最大破碎力的关系Fig.13 Relations between breakage energy and maximumcrushing force of coppermolybdenum ore particles 40 矿产保护与利用2023 年其实验数据均可用一条直线来拟合，最大破碎力断裂能和断裂应力与断裂比能之间的关系只与材料特性有关，而与颗粒尺寸无关，因此可以用这种特定的关系区分不同属性的物料。参考文献：FUERSTENAU D W,ABOUZEID A Z M.The energy efficiency of ballmilling in comminutionJ.International Journal of Mineral Processing,2002,67(1):161185.1 TAVARES L M.Handbook of Powder TechnologyM.Oxford:Elsevier,2007:368.2 POWELL M S,MORRISON R D.The future of comminutionmodellingJ.International Journal of Mineral Processing,2007,84(1/2/3/4):228239.3 VOGEL L,PEUKERT W.Breakage behaviour of differentmaterialsconstruction of a mastercurve for the breakage probabilityJ.Powder Technology,2003,129(1/2/3):101110.4 SALMAN,ALEXANDER RUSSELL,SERGEJ AMAN,et al.Breakageprobability of granules during repeated loadingJ.Powder Technology,2015,269:541547.5 AMAN S,JRGEN T,KALMAN H.Breakage probability of irregularlyshaped particlesJ.Chemical Engineering Science,2010,65(5):15031512.6 刘建远.威布尔分布在颗粒碎裂描述和粉碎数学建模中的应用J.矿冶,2009(3):18.LIU J Y.Application of Weibull distribution in particle fragmentationdescription and comminution mathematical modelingJ.Mining andMetallurgy,2009(3):18.7 ROZENBLAT Y,PORTNIKOV D,LEVY A,et al.Strength distributionof particles under compressionJ.Powder Technology,2011,208(1):215224.8 HUANG J,XU S,YI H,HU S.Size effect on the compression breakagestrengths of glass particlesJ.Powder Technology,2014,268(1):8694.9 GHADIRI M,BRUNARD N,KOLENDA F.Weibull analysis ofquasistatic crushing strength of catalyst particlesJ.ChemicalEngineering Research&Design,2003,81(8):953962.10 PETUKHOV Y,KALMAN H.Empirical breakage ratio of particlesdue to impactJ.Powder Technology,2004,143(26):160169.11 JRGEN T,SCHREIER M,TORSTEN G.Impact crushing of concretefor liberation and recyclingJ.Powder Technology,1999,105(1):3951.12 LIN Y L,WANG D M,LU W M.Compression and deformation of softspherical particlesJ.Chemical Engineering Science,2008,63(1):195203.13 ANTONYUK S,TOMAS J,HEINRICH S.Breakage behaviour ofspherical granulates by compressionJ.Chemical Engineering Science,2005,60(14):40314044.14 WEICHERT R.Theoretical prediction of energy consumption andparticle size distribution in grinding and drilling of brittle materialsJ.Particle&Particle Systems Characterization,1991,8(1):5562.15 MA L,LI Z,WANG M.Effects of size and loading rate on themechanical properties of single coral particlesJ.Powder Technology,2019,342:961971.16 黄俊宇.冲击载荷下脆性颗粒材料多尺度变形破碎特性研究D.合肥:中国科学技术大学,2016.HUANG J Y.Research on multiscale deformation and breakagecharacteristics of brittle particle materials under impact loadD.Hefei:University of Science and Technology of China,2016.17 周强,潘永泰,朱长勇,等.准静态下脆性材料强度分布试验研究J.煤炭学报,2019,44(S2):708716.ZHOU Q,PA