1、第 42 卷 第 10 期2023 年 10 月中国材料进展MATERIALS CHINAVol.42 No.10Oct.2023收稿日期:2021-08-18 修回日期:2021-12-06第一作者:刘 权,男,1997 年生,硕士 通讯作者:刘百宽,男,1961 年生,教授级高工,硕士生导师,Email:liubaikuan DOI:10.7502/j.issn.1674-3962.202108018基于正交试验的电工级氧化镁改性粉配比优化刘 权1,刘百宽1,2,田晓利2,李志勋2,许 雨1(1.西安建筑科技大学材料科学与工程学院,陕西 西安 710055)(2.濮阳濮耐高温材料(集团)股
2、份有限公司,河南 濮阳 457100)摘 要:为改善电工级氧化镁普粉电性能的不足,使改性后电工级氧化镁粉更加适配电热管制造行业,通过加入白炭黑、电工级氧化镁用添加剂和硅油的方法对普粉进行改性,采用正交试验法设计 25 组配比组合进行试验。通过极差法分析单一指标(冷态绝缘耐压强度、泄漏电流和热态绝缘耐压强度)下的最优配比,并利用矩阵分析法实现多目标综合分析优化,得到优化后的组合配比,运用指标值预估法预测最终配比的电工级氧化镁改性粉的电性能。结果表明:以冷态和热态绝缘耐压强度为指标时,白炭黑掺量为主要影响因素;以泄漏电流为指标时,硅油掺量为主要影响因素,其次是电工级氧化镁用添加剂掺量。多目标优化选
3、择的改性粉配比为白炭黑掺量0.25%(质量分数,下同),电工级氧化镁用添加剂掺量2%,硅油掺量 0.25%。最终配比电工级氧化镁改性粉电性能预估值为冷态绝缘耐压强度 3365 V,泄漏电流 0.112 mA,热态绝缘耐压强度 2399 V。关键词:电工级氧化镁改性;正交试验;极差法;矩阵分析法;指标值预估;最优配比中图分类号:TQ132.2 文献标识码:A 文章编号:1674-3962(2023)10-0826-07引用格式:刘权,刘百宽,田晓利,等.基于正交试验的电工级氧化镁改性粉配比优化J.中国材料进展,2023,42(10):826-832.LIU Q,LIU B K,TIAN X L,
4、et al.Optimization Research of Modified Electrical Grade Magnesium Oxide Powder Ratio Based on Orthogo-nal TestJ.Materials China,2023,42(10):826-832.Optimization Research of Modified Electrical GradeMagnesium Oxide Powder Ratio Based on Orthogonal TestLIU Quan1,LIU Baikuan1,2,TIAN Xiaoli2,LI Zhixun2
5、,XU Yu1(1.College of Materials Science and Engineering,Xian University of Architecture and Technology,Xian 710055,China)(2.Puyang Refractories Group Co.,Ltd.,Puyang 457100,China)Abstract:To improve its insufficient electrical property and make it more suitable for the electric heating tube manufactu
6、r-ing industry after modification,the plain electrical grade magnesium oxide powder was modified by adding silica,additive used in electrical grade magnesium oxide powder and silicone oil.Therefore,25 groups of ratios were designed to carry out combined experiments.On the basis of the orthogonal exp
7、eriment method.Analyze the optimal ratio under a single indicator(normal withstand voltage,leakage current and heating withstand voltage)with the range method.Also the matrix analysis method was used to analysis and optimize to obtain the optimal ratio,and the index value prediction method was used
8、to pre-dict the electrical property of the final ratio of modified electrical grade magnesium oxide powder.Finally,the results showed that the content of silica is the main influencing factor when the normal and heating withstand voltage is used as the indica-tor.However,when using leakage current a
9、s an indicator,the content of silicone oil is the main influencing factor,followed by the content of additive used in electrical grade magnesium oxide powder.In addition,the final parameters of the optimally selected backfill are as follows:0.25wt%of silica,2wt%of additive used in electrical grade m
10、agnesium oxide powder,0.25wt%of silicone oil.Its estimated values of electrical properties of the final ratio of modified electrical grade magnesium oxide powder are normal withstand voltage of 3365 V,leakage current of 0.112 mA,heating withstand voltage of 2399 V.Key words:electrical grade magnesiu
11、m oxide powder modification;orthogonal test;range method;matrix analysis method;index value prediction method;optimal mix proportion1 前 言菱镁矿是一种以碳酸镁为主要成分的天然矿物1,第 10 期刘 权等:基于正交试验的电工级氧化镁改性粉配比优化根据结晶状态不同,一般分为晶质和非晶质 2 种2。我国川藏地区的菱镁矿矿石致密,细腻有润滑感,坚硬但质脆,具有均匀的陶土状宏观结构,同时具有典型的晶格常数和完整的结晶形貌,被称为微晶菱镁矿3。普通电熔镁砂可以采用菱镁矿作
12、为原料,在电熔炉中经熔融和缓慢冷却而制得4-6,具有纯度高、熔点高(2800)、结晶晶粒大、结构致密、耐高温、化学性能稳定、耐压强度大和绝缘性能好等特点,基于这些优良性能,该材料广泛应用于高温电气绝缘材料7。电工级氧化镁是由电熔氧化镁制备而来的,具有高电阻、高导热性、耐高压和低泄漏电流等优良电热特点,广泛应用于家用及工业用电热电器,是理想的绝缘材料8。一直以来,国内电工级氧化镁行业的发展相对于电热管行业存在着滞后。随着国内对电加热元件品质要求的提升,许多企业生产的电工级氧化镁质量并不达标,主要原因有原矿开采品位无法保证、选用浮选法优化低品位菱镁矿、电熔氧化镁冶炼时间变短以及以电熔氧化镁皮砂为原
13、料9。为了提高电工级氧化镁的品质,相关研究主要从两方面进行。一方面,提升电熔氧化镁的品质,通过提升电熔冶炼技术以及优化配制电熔时的菱镁矿矿石配比,从而合理控制电熔氧化镁的各项化学指标。另一方面,电工级氧化镁改性,通过有机硅化物对电工级氧化镁进行表面处理,达到防潮效果,主要研究包括有机硅化物的种类、用量、涂覆方式以及新的耐高温高分子材料的研制等;通过添加无机物添加剂提高电工级氧化镁电性能,主要研究包括无机物添加剂的种类、用量及制备方法等10,11。目前以川藏地区微晶菱镁矿为原料制备电熔氧化镁的相关研究已经开展12,以此电熔氧化镁为原料制备电工级氧化镁的相关研究也开始起步。本试验以提升电工级氧化镁
14、的电性能为目标,在电工级氧化镁中掺加白炭黑、电工级氧化镁用添加剂及硅油,对它进行改性,通过正交试验设计配比试验方案。在通过极差法和方差法分析每种指标(冷态绝缘耐压强度、泄漏电流和热态绝缘耐压强度)影响因素的主次顺序、影响程度以及最优配比的基础上,利用多指标单一化的矩阵分析法确定最终配比,并运用指标值预估法预估最终配比改性粉的电性能,从而为川藏地区微晶菱镁矿的后期大规模工业生产电工级氧化镁提供理论依据和指导。2 实 验2.1 原材料电熔氧化镁购自青海濮耐高新材料有限公司,体积密度为 3.47 g/cm3,主要化学成分和物相组成见表 1 和图 1;白炭黑、电工级氧化镁用添加剂和硅油分别购自沈阳市华
15、铬新化工产品有限公司、营口松辽镁业公司和上海三荣化工科技有限公司;其中电工级氧化镁用添加剂的主要成分是含镁硅酸盐13。图 1 电熔氧化镁原料的物相组成Fig.1 The XRD pattern of the fused magnesium oxide表 1 电熔氧化镁原料的化学成分Table 1 Chemical compositions of the fused magnesium oxide(w/%)CompositionMgOCaOSiO2Fe2O3Al2O3TiO2K2ONa2OPercentage97.561.580.580.070.030.030.0302.2 试样制备试验选取白炭
16、黑、电工级氧化镁用添加剂、硅油3 种改性材料的掺量为考察因子,分别以 A、B、C 表示,每个因子设置 5 个水平,选择 L25(56)正交表前 4 列设计正交试验,D 为空白组,如表 2 所示。电工级氧化镁改性粉试样制备流程如下:将电熔氧化镁原料破碎筛分至 40325 目,制得电工级氧化镁普粉。根据正交试验设定的掺量添加改性材料,混合均匀,即可得电工级氧化镁改性粉,按照正交表分为 25 组,依次标记为 A1A25。表 2 正交试验因素水平表Table 2 The factor level table of orthogonal experimentLevelfactorABCDSilicaco
17、ntent/(wt%)Additive used inelectrical grademagnesium oxide powdercontent/(wt%)Silicone oilcontent/(wt%)Blank10.050.50.05020.110.1030.151.50.15040.220.2050.252.50.250728中国材料进展第 42 卷2.3 电性能测试将制得的 25 组试样装管、缩管、平头、1050 过炉20 min 后制得相应的电热管,每组 3 支,共 75 支。对每支电热管进行相关电性能检测分析(管表负荷为 9 W/cm3),利用 ZHZ8A 型耐电压测试仪检测冷态
18、和热态绝缘耐压强度,利用 PA30 型数字泄漏电流测试仪检测泄漏电流。以 3 次检测的平均值作为结果,试验结果中超出平均值10%的,应剔除后再取平均值作为试验结果。3 结果与讨论3.1 正交试验结果冷态绝缘耐压强度、泄漏电流和热态绝缘耐压强度均为电工级氧化镁粉的电性能指标,能够表征其绝缘性能,绝缘耐压强度越高越好,泄漏电流越小越好。通过上述试验测量并计算获得了冷态和热态绝缘耐压强度和泄漏电流参数,如表 3。表 3 正交试验结果Table 3 The results of orthogonal experimentGroupABCDNormalwithstandvoltage/VLeakagec
19、urrent/mAHeatingwithstandvoltage/VA1111123900.3492065A2122227930.2012270A3133328470.1372260A4144426170.1241960A5155527200.1081980A6212330700.2402130A7223432130.1402205A8234524230.1132255A9245129700.0982105A10251228400.2001995A11313528050.2192240A12324126870.1612075A13335232430.1352070A14341328330.21
20、12160A15352429330.1752340A16414227950.2052380A17425333400.1582300A18431431400.2222260A19442531170.1622340A20453130030.1402313A21515428300.1882565A22521530500.2512400A23532134830.1872423A24543234150.1382520A25554336800.12123803.2 极差分析极差分析简称 R 法,根据 R 的大小可以判断因素重要性的主次,R 越大,该因素对试验指标影响越大。极差的计算如式(1)和式(2)所示
21、:kmn=1NNi=1Pi(1)Rm=max(km1,km2,kmn)-min(km1,km2,kmn)(2)式中:kmn为第 m 个因素第 n 水平对应指标值的平均值;Pi为指标值;Rm为第 m 个因素的极差值。根据式(1)和式(2)对表 3 正交试验结果进行极差分析,结果见表 4,D 代表未考虑的因素,如果该因素的极差大于已知因素,则说明还有更重要的因素未考虑到,反之亦然。根据表 4 做出指标趋势图,见图 2。从表 4 可以看出,以冷态绝缘耐压强度为指标时,影响因素的主次顺序为:ADBC,其中未考虑因素影响不是最大的,根据 k 值大小选取最优解为 A5B5C2。以泄漏电流为指标时,影响因素
22、的主次顺序为:CBAD,其中未考虑因素影响不是最大的,根据 k 值大小选取最优解为A2B4C5。以热态绝缘耐压强度为指标时,影响因素的主次顺序为:ACBD,其中未考虑因素影响不是最大的,根据 k 值大小选取最优解为 A5B1C3。从图 2a 和 2c表 4 极差分析结果Table 4 The results of range analysisFactorLevelABCDNormalwithstandvoltage/Vk12673277828512907k22903301730793017k32900302730573154k43079299028402947k5329230353021282
23、3R619257239331Leakagecurrent/mAk10.1840.2400.2470.187k20.1580.1820.1930.176k30.1800.1590.1560.173k40.1780.1470.1450.170k50.1770.1490.1380.171R0.0260.0930.1090.017Heatingwithstandvoltage/Vk12107227621762196k22138225023012247k32177225423082246k42319221722102266k52458220222042243R3517413270828 第 10 期刘
24、权等:基于正交试验的电工级氧化镁改性粉配比优化图 2 因素与性能指标的关系:(a)冷态绝缘耐压强度指标,(b)泄漏电流指标,(c)热态绝缘耐压指标Fig.2 Relationship between factors and performance indicators:(a)normal withstand voltage,(b)leakage current,(c)heating withstand voltage可以看出,随着白炭黑掺量的增多,绝缘耐压强度不断增大。白炭黑基于其小粒径和高表面能来吸附在粉体表面并形成一个表层,提高了粉体分散性,从而使电工级氧化镁获得较好的流动性,提高其电性能
25、14。从图 2b 中可以看出,随着电工级氧化镁用添加剂掺量的增多,泄漏电流先减小后增大,在 B4 出现拐点,而电工级氧化镁用添加剂通常利用硅、镁化合物进行合成,可通过调节电工级氧化镁化学成分尤其是硅和钙的物质的量的比来提高其电性能15,说明 B4 对应的硅和钙的物质的量的比最合适。从图 2b 中还可以看出,随着硅油掺量的增多,泄漏电流不断减小,但减小速率有所下降,硅油是通过吸附作用包裹于颗粒表面,显著提高粉体表面的亲油性,从而达到改善氧化镁的防潮效果,提高电工级氧化镁电性能16。3.3 方差分析方差分析是一种检验多个总体均值是否相等的统计方法,对表 3 正交试验结果用 Minitab 软件进行
26、方差分析,结果见表 5。由表 5 中计算得到的 F 值与给定显著水平的标准 F 值相比较可知,A 对改性粉的冷态绝缘耐压强度影响显著,对热态绝缘耐压强度影响明显,即在工作状态下,A 对绝缘耐压强度影响更显著;对泄漏电流影响明显的是 C,其次是 B。表 5 方差分析结果Table 5 The results of variance analysisIndexFactorSquaresumDOFFCriticalvalueNormalwithstandvoltage/VLeakagecurrent/mAHeatingwithstandvoltage/VA1 063 06043.96B234 866
27、40.88C265 17940.99A0.001 9842.54B0.0303438.88 C0.0409452.53 A427 945410.24 B17 91940.43C72 66741.74F0.10=2.48,F0.05=3.26,F0.01=5.41 Notes:F represents this factor influences index significantly;F represents this factor influences index modestly;F represents this factor influences index small.3.4 矩阵分
28、析矩阵分析法主要用于解决多指标正交试验设计的方案优化问题,能够将多个指标单一化,使多指标问题简化。本工作将冷态绝缘耐压强度、泄漏电流和热态绝缘耐压强度 3 个指标单一化,获取 3 个指标下的最优配比928中国材料进展第 42 卷以及该配比下的因素影响重要性顺序。该方法需建立一个 3 层结构模型:指标层、因素层和水平层,根据各个层次的数据,矩阵定义如下17:首先建立指标层,矩阵 M:M=K1100K1200K1n000K2100K2200K2n000Km100Km200Kmn(3)式中,当指标越大越好时,Kmn=kmn;当指标越小越好时,Kmn=1/kmn。令 tj=1ni=1Ki,再建立因素层
29、矩阵 T:T=t1000t20000tm(4)rj=Rjmi=1Ri,最后建立水平层矩阵 S:S=r1r2rm(5)在 p 个指标中第 i 个指标的矩阵如下:i=MiTiSi(6)最终结果为:=1+2+pp(7)根据以上公式可以得出矩阵分析结果,如表 6 所列。从表 6 中可以看出每个指标及平均值下的结果总和均等于 1,说明每个指标的重要性一致,在各指标重要性一致的条件下,各因素各水平下指标均值最大值分别为 A5(0.0796)、B4(0.0489)、C5(0.0611)和 D3(0.0284)。根据各数值大小可以判断出影响因素重要顺序为 ACBD,D 影响较小,因此最优解为 A5B4C5。表
30、 6 矩阵分析结果Table 6 The results of matrix analysisGroupMatrix analysis resultsNormalwithstandvoltageLeakage currentHeating withstandvoltageIndex averageA10.07710.02020.10530.0675A20.08370.02350.10690.0714A30.08360.02060.10880.0710A40.08880.02090.11590.0752A50.09490.02100.12290.0796B10.03330.05370.02400
31、.0370B20.03610.07080.02370.0435B30.03620.08100.02380.0470B40.03580.08760.02340.0489B50.03630.08650.02320.0487C10.03170.06050.04090.0444C20.03430.07740.04330.0517C30.03400.09580.04340.0577C40.03160.10300.04150.0587C50.03360.10820.04140.0611D10.04480.01300.02190.0266D20.04650.01380.02240.0276D30.04860
32、.01410.02240.0284D40.04540.01430.02260.0274D50.04350.01420.02240.02673.5 指标值预估指标值预估是指按照已经确定的方案通过特定的公式进行各指标值预估,其计算公式如下18:=x-+ai+bj+ck+(8)式中,为指标值预估值;x-为试验数据值的总平均值;ai、bj、ck为各因素在最优水平下的效应值。其中各效应值计算如下:ai=Ai-x-bj=Bj-x-ck=Ck-x-(9)式中,Ai、Bj、Ck为各因素在最优水平下的指标平均值。将式(9)代入式(8)可得如下公式:=mi=1Vi-(m-1)x-(10)038 第 10 期刘 权
33、等:基于正交试验的电工级氧化镁改性粉配比优化式中,m 为因素数;Vi为各因素所选水平对应的指标平均值。根据以上公式可以求出各最优解的预估值,其中以矩阵分析法得出的最优解为例,计算结果如表 7,各最优解预估值结果如表 8。从表 8 中可以看出,根据冷态绝缘耐压强度得出的最优解的预估值中泄漏电流较大;根据泄漏电流得出的最优解的预估值中冷态和热态绝缘耐压强度都是最小的,与试验值相比基本上一致;根据热态绝缘耐压强度得出的最优解的预估值中泄漏电流是最大的。综合所有因素来看,矩阵分析法得出的最优解的预估值与单一指标得出的最优解的预估值相比,各项指标值均接近最佳指标值,因此,矩阵分析法得到的最优解的预估值最
34、好。选取矩阵分析法得到的最优解 A5B4C5 进行各项性能测试,测试结果为冷态绝缘耐压为 3349 V,泄漏电流为 0.115 mA,热态绝缘耐压强度为 2420 V。该测试结果与指标值预估法得到的预估值基本一致,表明指标值预估法合理。表 7 矩阵分析法最优解的预估值过程及结果Table 7 Process and results of optimal solution estimation by matrix analysis methodFactorABCVx-Estimated valueNormal withstand voltage/V329229903021930329693365
35、Leakage current/mA0.1770.1470.1380.4620.1750.112Heating withstand voltage/V245822172204687922402399表 8 各最优解的预估值Table 8 Estimated values of each optimal solutionFactorEstimated valueOptimal solutionNormal withstandvoltage/VLeakage current/mAHeating withstandvoltage/VNormal withstand voltageA5B5C23469
36、0.1692481Leakage currentA2B4C529760.0932079Heating withstand voltageA5B1C331890.2232562Matrix methodA5B4C533650.11223994 结 论(1)极差分析试验结果表明:各因素对电工级氧化镁改性粉冷态绝缘耐压强度影响程度的顺序为:白炭黑掺量电工级氧化镁用添加剂掺量硅油掺量;对改性粉泄漏电流影响程度的顺序为:硅油掺量电工级氧化镁用添加剂掺量白炭黑掺量;对改性粉热态绝缘耐压强度影响程度的顺序为:白炭黑掺量硅油掺量电工级氧化镁用添加剂掺量。(2)方差分析试验结果表明:白炭黑掺量对冷态和热态绝缘耐
37、压强度的影响最为显著,并且在工作状态下,白炭黑掺量对绝缘耐压强度影响更大;泄漏电流的主要影响因素是硅油掺量,其次是电工级氧化镁用添加剂掺量。(3)利用矩阵分析法优化得到的改性粉最终配比为:白炭黑掺量 0.25%,电工级氧化镁用添加剂掺量 2%,硅油掺量 0.25%(皆为质量分数)。(4)运用指标值预估法计算得到最终配比下的指标值分别是冷态绝缘耐压强度 3365 V,泄漏电流 0.112 mA,热态绝缘耐压强度 2399 V。参考文献 References1 高陟,马北越,任鑫明,等.耐火材料J,2020,54(1):88-92.GAO Z,MA B Y,REN X M,et al.Refrac
38、toriesJ,2020,54(1):88-92.2 胡庆福.镁化合物生产与应用M.北京:化学工业出版社,2004:11-77.HU Q F.Production and Application of Magnesium CompoundsM.Beijing:Chemical Industry Press,2004:11-77.3 郭杰,刘百宽,田晓利,等.非金属矿J,2019,42(5):45-48.GUO J,LIU B K,TIAN X L,et al.Non-Metallic Mines J,2019,42(5):45-48.4 蒋小蒙.菱镁矿熔炼电熔镁砂节能新工艺及方镁石结晶过程研究
39、138中国材料进展第 42 卷D.沈阳:东北大学,2011.JIANG X M.Study on Advanced Energy Conservation Technology for Fused Magnesia Smelt from Magnesite and the Crystallization Process of PericlaseD.Shenyang:Northeastern University,2011.5 FU Y,WANG N H,WANG Z,et al.Journal of Materials Process-ing TechnologyJ,2017,244:231-
40、239.6 CHAI T Y,WU Z W,WANG H.IFAC PapersOnLineJ,2017,50(1):14992-14999.7 罗仙平,李建康,徐徽,等.盐业与化工J,2016,45(8):8-13.LUO X P,LI J K,XU H,et al.Journal of Salt and Chemical In-dustryJ,2016,45(8):8-13.8 胡庆福,宋丽英,刘宝树,等.无机盐工业J,2009,41(4):8-10.HU Q F,SONG L Y,LIU B S,et al.Inoaganic Chemicals IndustryJ,2009,41(4)
41、:8-10.9 宗俊.特种镁化物(氧化镁)的现状与发展趋势C/2020 年镁化合物行业年会暨专家工作会议论文集.河北:中国无机盐工业协会镁化合物分会,2020:11-22.ZONG J.The Current Situation and Development Trend of Special Magnesium Compounds(Magnesium Oxide)C/2020 Magnesium Compound Industry Annual Meeting and Expert Working Conference Proceedings.Hebei:Magnesium Compound
42、 Branch of China Inorgan-ic Salt Industry Association,2020:11-22.10 胡庆福,胡晓湘,宋丽英,等.化工进展J,2005,24(1):28-31.HU Q F,HU X X,SONG L Y,et al.Chemical Industry and Engi-neering ProgressJ,2005,24(1):28-31.11 胡庆福,刘景泽,宋丽英,等.无机盐工业J,2006,38(9):13-16.HU Q F,LIU J Z,SONG L Y,et al.Inoaganic Chemicals IndustryJ,200
43、6,38(9):13-16.12 杨树,刘百宽,田晓利,等.硅酸盐通报J,2019,38(9):2777-2781.YANG S,LIU B K,TIAN X L,et al.Bulletin of the Chinese Ce-ramic SocietyJ,2019,38(9):2777-2781.13 季克捷,杜永平,王洪元.一种电工级氧化镁粉添加剂的制造方法:CN107200591BP.2020-04-07.JI K J,DU Y P,WANG H Y.A Preparation Method for Electrical Grade Magnesium Oxide Powder Add
44、itive:CN107200591BP.2020-04-07.14 倪忠伟,王继奇,史生川,等.一种防潮改性电工级氧化镁的制备方法以及应用:CN109135339A P.2019-01-04.NI Z W,WANG J Q,SHI S C,et al.A Preparation Method and Application for Moisture-Proof Modified Electrical Grade Magnesium Oxide Powder:CN109135339AP.2019-01-04.15 宗俊.电工级氧化镁发展现状和趋势C/2006 年中国镁盐行业年会暨新技术、新产品、
45、新设备推介会论文集.河北:中国无机盐工业协会钙镁盐分会,2006:71-76.ZONG J.Development Status and Trend of Electrical Grade Magnesi-um Oxide PowderC/2006 China Magnesium Salt Industry Annual Conference and New Technology,New Products,New Equipment Pro-motion Conference Proceedings.Hebei:Calcium and Magnesium Salt Branch of Chin
46、a Inorganic Salt Industry Association,2006:71-76.16 郑水林,王彩丽,李春全.粉体表面改性M.北京:中国建材工业出版社,2019:12.ZHENG S L,WANG C L,LI C Q.Powder Surface ModificationM.Beijing:China Building Materials Industry Press,2019:12.17 周玉珠.数学的实践与认识J,2009,39(2):202-207.ZHOU Y Z.Mathematics in Practice and TheoryJ,2009,39(2):202-207.18 赵风文,胡建华,曾平平,等.中国有色金属学报J,2021,31(4):1096-1105.ZHAO F W,HU J H,ZENG P P,et al.The Chinese Journal of Nonferrous MetalsJ,2021,31(4):1096-1105.(编辑 费蒙飞)238