Processing math: 5%

过渡金属和磷共掺杂多孔碳作为氧气还原电催化剂

万泽远, 李桂林, 吴娇

万泽远, 李桂林, 吴娇. 过渡金属和磷共掺杂多孔碳作为氧气还原电催化剂[J]. 复合材料学报, 2023, 40(2): 844-851. DOI: 10.13801/j.cnki.fhclxb.20220228.002
引用本文: 万泽远, 李桂林, 吴娇. 过渡金属和磷共掺杂多孔碳作为氧气还原电催化剂[J]. 复合材料学报, 2023, 40(2): 844-851. DOI: 10.13801/j.cnki.fhclxb.20220228.002
WAN Zeyuan, LI Guilin, WU Jiao. Transition metal and phosphorus co-doped porous carbon as electrocatalyst for oxygen reduction[J]. Acta Materiae Compositae Sinica, 2023, 40(2): 844-851. DOI: 10.13801/j.cnki.fhclxb.20220228.002
Citation: WAN Zeyuan, LI Guilin, WU Jiao. Transition metal and phosphorus co-doped porous carbon as electrocatalyst for oxygen reduction[J]. Acta Materiae Compositae Sinica, 2023, 40(2): 844-851. DOI: 10.13801/j.cnki.fhclxb.20220228.002

过渡金属和磷共掺杂多孔碳作为氧气还原电催化剂

基金项目: 山西省基础研究计划项目(201901D211157);山西省高校科技创新项目(2020L0044)Shanxi Basic Research Program Project (201901D211157); Science and Technology Innovation Project of Colleges and Universities in Shanxi Province (2020L0044)
详细信息
    通讯作者:

    吴娇,博士,副教授,硕士生导师,主要研究方向为高性能电催化材料 E-mail:jiaow@sxu.edu.cn

  • 中图分类号: TM911.4

Transition metal and phosphorus co-doped porous carbon as electrocatalyst for oxygen reduction

  • 摘要: 碳基材料作为非贵金属催化剂具有导电性能高、稳定性能好、价格低廉、环境友好等优点,在燃料电池阴极催化剂领域中引起了广泛的关注,尤其是过渡金属和异原子共掺杂能够显著提高碳材料的氧气还原催化活性。本文采用聚醚(F127)作为软模版,苯酚、甲醛作为碳源,四苯基溴化膦作为磷源,硝酸盐作为过渡金属来源,通过挥发溶剂自组装及高温煅烧过程制备了过渡金属(Co、Fe、Ni、Mn)和磷(P)共掺杂多孔碳材料(TM-P-C)。通过旋转环盘电极研究了TM-P-C在0.1 mol/L KOH电解液中的氧气还原电催化性能。研究结果表明:TM-P-C催化剂具有较高的氧化还原反应(ORR)电催化性能,其ORR活性为P-Co-C>P-Ni-C>P-Fe-C>P-Mn-C,其中P-Co-C的ORR电催化性能可与商业20wt%Pt/C催化剂相媲美,其电流密度与20wt%Pt/C催化剂的电流密度相当,与20wt%Pt/C仅存在66 mV的半波电位差,表现为ORR的4e转移途径。制备的TM-P-C催化剂所具有的较高氧气还原电催化活性主要来自于过渡金属和P原子之间的协同作用。此外,TM-P-C催化剂表现出优异的长期稳定性和抗甲醇毒化性能,优于商业化20wt%Pt/C催化剂。
    Abstract: Carbon-based materials, as non-noble catalysts, have attracted extensive attention in the field of cathode catalysts for fuel cells due to their high conductivity, long-term stability, low-cost and environmental friendliness. Especially, the oxygen reduction reaction (ORR) activity of carbon materials can be significantly improved after co-doping with transition metal and heteroatoms. In this work, transition metal (Co, Fe, Ni, Mn) and phosphorus (P) co-doped porous carbon (TM-P-C) was prepared via self-assembly method combined with a carbonization process, in which polyether (F127) was introduced as soft template, phenol and formaldehyde as carbon precursor, tetraphenylphosphine bromide as phosphorus source, and nitrate as transition metal source. The electrocatalytic activity of TM-P-C for ORR in alkaline electrolyte was studied by using the rotating ring-disk electrode (RRDE) technique. The results reveal that TM-P-C exhibits high electrocatalytic performance for ORR in 0.1 mol/L KOH, and the activity follows P-Co-C>P-Ni-C>P-Fe-C>P-Mn-C. Moreover, the ORR performance of P-Co-C is compared to that of commercial 20wt%Pt/C catalyst. The diffusion limiting current density of P-Co-C reaches that of 20wt%Pt/C and a negative shift of about 66 mV exists in the half-wave potential of P-Co-C as compared to 20wt%Pt/C, indicating the four-electron pathway during the ORR. The enhancement in the activity for ORR is mainly attributed to the synergistic effect of P and transition metal doping in carbon of TM-P-C. Moreover, TM-P-C shows excellent long-term stability and methanol toxicity resistance, superior to that of commercial 20wt%Pt/C.
  • 树脂基纤维增强复合材料(简称“复材”)具有轻质、高强、耐腐蚀等突出优点,被广泛应用于建筑、交通、能源、航空航天及体育等行业[1]。全球复合材料市场规模预计将从2019年的906亿美元增长至2024年的1316亿美元[2]。随着复合材料应用日益广泛,制造期间产生的废料、测试样品以及退役产品等复合材料固体废弃物(复材固废)总量也逐年攀升[3-4]。以风电领域为例,根据我国现有风电装机容量和未来总体规划预测,2010至2024期间年我国风力发电叶片废弃物累计总量为1~1.3万吨,其后5年将增至19.7万吨,预计到2040年累计总量将达到240万吨左右[5]。这些复材废弃物不可降解,其数量的快速增长将对环境造成巨大压力[6]

    传统复材固废的处理方法包括直接填埋和焚烧。前者简单便捷、应用最广,但占用大量土地,并对填埋地的自然环境造成持久性污染[6]。各国日益收紧的固废填埋政策也导致固废填埋费用逐步攀升[7-8]。焚烧处理成本同样较高,且会产生有毒有害气体,而原固废体积一半左右的残余灰烬仍需填埋处理[9-11]

    常见复材固废回收再利用方法主要有化学回收法、热解回收法以及机械回收法[6, 11-13]。化学回收法采用不同溶剂使聚合物的化学键断裂,从而将树脂基体与纤维分离[14]。化学法成本高、耗时长,目前处于实验室研究阶段,尚未推广应用[6]。热解法一般指采用450℃以上的高温将复材固废进行裂解,从而将纤维从复材固废中分离出来[11]。该法适用于化学性能稳定且产品附加值较高的碳纤维复材固废(碳纤固废)的再生利用[15-16],已成功商业化,采用此法对碳纤固废进行再生利用的公司包括英国ELG Carbon Fiber Ltd.、意大利Karborek Spa、美国Material Innovation Technology和国内的复源新材料科技有限公司等[11, 17]。然而,热解法对回收利用占总体复材市场规模95%以上的玻璃纤维(玻纤)复材固废并不经济适用[6]。其原因为玻璃纤维附加值低,且力学性能在回收处理后显著降低[8, 11]

    机械回收法指通过破碎、粉碎、研磨等工序将复材固废处理成不同形状的回收产物,并作为填充材料[12, 18-20]或抗裂增强材料[21-23]用于制备水泥基产品或其他产品。有学者将复材固废处理成不同尺寸的块体、颗粒以及粉末,用于部分替代砂石骨料以制备混凝土[12, 18-20]。相比砂石骨料,玻纤固废粉碎料的弹性模量小、与包裹基体的粘结性能弱,从而导致混凝土强度相对较低[12, 18-20]。将复材固废处理成短条状或细纤维,作为抗裂增强材料制备纤维混凝土,所获混凝土的力学性能要优于基于粗、细骨料替代法制备的混凝土[21-24]。综上所述,机械回收法可能是解决玻纤固废回收利用难题的有效方法之一。

    本研究团队此前已经提出了一种新的玻纤复材固废回收利用方法:通过机械切割及铣削将玻纤复材固废加工成长度30~100 mm,宽度3~5 mm,厚度0.5~1.5 mm的粗纤维,并用于制备玻纤复材固废粗纤维混凝土(Macro fibre reinforced concrete,MFRC)[25]。目前,关于MFRC的研究还处于起步探索阶段,仅有少数文献对MFRC材料开展了试配试验研究和基本力学性能测试[25-30]。结果表明,玻纤复材固废粗纤维具有作为混凝土中抗裂增强材料的潜力,然而关于粗纤维对混凝土力学性能的影响规律尚未有系统的研究成果。鉴于此,本文通过一系列轴压、劈拉试验,研究了粗纤维体积掺量、厚度、长度对两种配比混凝土的力学性能的影响,评估了现有FRC劈拉强度预测公式对MFRC劈拉强度的预测效果,并基于试验结果提出了MFRC劈拉强度预测公式。

    本研究所用粗纤维由某废弃风机叶片通过机械切割及铣削处理得到,主要步骤包括:(1)将完整叶片解体为若干便于运输的粗切节段,并将其切割为长度为0.5~3 m的板材(图1(a));(2)将板材切割为长度30~100 mm的块体(图1(b));(3)采用自动铣削机床(图1(c)),顺着原始叶片的纤维铺层主方向将块体进一步切割成长度30~100 mm,宽度3~5 mm,厚度0.5~1.5 mm的粗纤维(图1(d))。经测试,本文所用粗纤维的密度为1.97 g/cm3

    图  1  风机叶片制备粗纤维的机械方法
    Figure  1.  Mechanical method for producing macro fibres from a wind turbine blade
    GFRP—Glass fibre reinforced polymer

    光学显微镜(OM)和扫描电镜(SEM)图像显示:粗纤维由树脂基体和纤维两相组成,大部分纤维丝嵌置于树脂基体中,但也有少量纤维丝裸露于粗纤维表面,粗纤维表面存在纤维丝损伤断裂和树脂基体剥落的现象(图2)。回收粗纤维的力学性能低于新出厂的风机叶片,主要原因为:(1)叶片在服役期间的性能退化;(2)运输退役叶片过程中造成的损伤;(3)切割及铣削处理造成的损伤[25]。单根粗纤维可能包含多层原始复合材料铺层,而层间在铣削加工剪切力作用下容易受到损伤,因而发生层间脱离现象,造成粗纤维力学性能的降低。对于本研究所采用的机械加工工艺,所制备的粗纤维厚度超过1.5 mm时层间脱离现象较为明显。

    图  2  粗纤维表面的粗糙纹理
    Figure  2.  Rough texture of macro fibre surface

    本研究采用5种不同尺寸的粗纤维,从中各随机挑选200根以对其实际几何尺寸进行了统计(表1)。对于名义尺寸为90 mm×4.5 mm×0.5 mm的粗纤维,其实际平均长度、宽度和厚度分别为91.20、4.49和0.54 mm。其余4种规格粗纤维的实际几何尺寸也接近对应的名义尺寸。

    表  1  回收粗纤维统计参数
    Table  1.  Statistical properties of recycled macro fibres (MF)
    Designation Nominal
    dimensions/
    mm
    Length/mm Width/mm Thickness/mm Tensile strength/MPa Elastic modulus/GPa
    Average Standard
    deviation
    Average Standard
    deviation
    Average Standard
    deviation
    Average Standard
    deviation
    Average Standard
    deviation
    MF-1 90×4.5×0.5 91.20 1.94 4.49 0.43 0.54 0.16 554 135 37.7 7.3
    MF-2 90×4.5×1.2 90.57 4.56 4.57 0.36 1.26 0.26
    MF-3 90×4.5×2.2 87.60 2.97 4.57 0.75 2.26 0.53
    MF-4 60×4.5×0.5 58.25 1.94 4.63 0.80 0.54 0.19
    MF-5 30×4.5×0.5 33.30 1.20 4.81 0.48 0.50 0.12
    下载: 导出CSV 
    | 显示表格

    根据ASTM 3039标准[31]对9个粗纤维测试试件的静力拉伸性能进行测试。拉伸试件名义尺寸为150 mm×4.5 mm×1 mm,名义标距为60 mm。为防止夹持部位发生破坏,采用环氧树脂将尺寸为45 mm×34 mm×1 mm的铝片粘贴于拉伸试件端部。拉伸试验破坏后的粗纤维试件如图3所示。粗纤维平均拉伸强度为554 MPa,标准差为135 MPa,其弹性模量平均值为37.7 GPa,标准差为7.3 GPa。可见粗纤维在力学性能上有显著的离散性,这与前期研究中的观察一致[25, 32]。其离散性来源于:(1)固废来源不同;即便来源于同一叶片,不同部位的纤维铺层方向、纤维含量也不同;(2)切割和铣削处理造成的损伤差异。

    图  3  拉伸破坏后的粗纤维试件
    Figure  3.  Macro fibre specimens after tensile failure

    本文设置两批次共11组试件,如表2所示。批次1的试件用于研究粗纤维体积分数对MFRC轴压和劈拉性能的影响,批次2的试件用于研究粗纤维长度、厚度以及混凝土配合比的影响。每组制备尺寸均为150 mm×300 mm(直径×高度)的圆柱抗压试件和劈拉试件各3个。在试件组别的命名中,“S”代表批次,“W”代表水灰比,“L”代表纤维长度,“V”代表纤维体积分数,“T”代表纤维厚度。例如,S1W0.32L90V-0.5%T0.5-1代表第一批,混凝土水灰比为0.32,粗纤维长度为90 mm,粗纤维体积分数为0.5vol%,粗纤维厚度为0.5 mm的同组3个试件中的第一个。此外,为便于后续讨论,每个组别也被赋予一个简化编号。例如,组别1-3代表批次1中的第3组。

    表3为本研究采用的两种依据普通混凝土配合比设计规程JGJ 55—2011[33] 确定的配合比(C60和C40);在本文中,对混凝土配比无特别说明时,所讨论对象均为配比一混凝土。在预试验中,掺入1.5vol%粗纤维的C40混凝土的坍落度由180 mm急剧下降至30 mm以下,且掺入2vol%粗纤维的混凝土搅拌困难,未能成型。原因在于粗纤维尺寸较大且表面粗糙,容易“结团”。因此,在正式试验中,除原C40配比(水灰比0.47,砂率0.4)外,新加入了一种C60配比(水灰比0.32,砂率0.45)。C60配比的水灰比较低,胶凝材料用量大,容易形成粘稠的浆体,因此提升砂率可以改善流动性,减少空隙率,提升密实度。因掺入粗纤维会显著影响混凝土的流动性,在试配试验中将不含粗纤维的两种配比混凝土的目标坍落度定为200 mm,以给后续因粗纤维而导致的坍落度下降提供空间。此外,通过试配确定了两种混凝土配比的减水剂用量。两种不同配比混凝土的实测坍落度分别为220 mm和195 mm,接近预期的200 mm。

    表  2  试件设计
    Table  2.  Specimen design
    No. Mix Water/cement ratio Fibre length/mm Fibre thickness/mm Fibre volume
    fraction/vol%
    No. of specimens
    Compression Splitting
    1-1 S1W0.32L0V0%T0 0.32 0 3 3
    1-2 S1W0.32L90V0.5%T0.5 0.32 90 0.5 0.5 3 3
    1-3 S1W0.32L90V1%T0.5 0.32 90 0.5 1 3 3
    1-4 S1W0.32L90V1.5%T0.5 0.32 90 0.5 1.5 3 3
    1-5 S1W0.32L90V2%T0.5 0.32 90 0.5 2 3 3
    2-1 S2W0.32L90V1.5%T1.2 0.32 90 1.2 1.5 3 3
    2-2 S2W0.32L90V1.5%T2.2 0.32 90 2.2 1.5 3 3
    2-3 S2W0.32L60V1.5%T0.5 0.32 60 0.5 1.5 3 3
    2-4 S2W0.32L30V1.5%T0.5 0.32 30 0.5 1.5 3 3
    2-5 S2W0.47L90V1.5%T0.5 0.47 90 0.5 1.5 3 3
    2-6 S2W0.47L0V0%T0 0.47 0 3 3
    下载: 导出CSV 
    | 显示表格
    表  3  混凝土基体配合比
    Table  3.  Mix proportions of concrete matrix
    Mix Strength grade Fresh water/
    cement ratio
    Fresh water/
    (kg·m−3)
    Cement/
    (kg·m−3)
    Fine aggregate/
    (kg·m−3)
    Coarse aggregate/
    (kg·m−3)
    Water reducer/%
    1 C60 0.32 160 500 795 970 0.45
    2 C40 0.47 185 395 715 1075 0.19
    下载: 导出CSV 
    | 显示表格

    粗骨料为5~20 mm碎石集料,细骨料为天然河砂,水泥采用强度等级为42.5R的石井牌普通硅酸盐水泥。在混凝土浇筑前对碎石集料进行清洗以去除其表面淤泥及其他杂质,河砂和碎石集料经过阳光晒干以去除其内部水分。在混凝土搅拌过程中,先将其他组份拌合均匀,然后继续搅拌1 min,在此期间匀速加入粗纤维。对坍落度进行测试后,将混凝土拌合物浇筑至模具内,24 h后对试件进行脱模,然后将其置于湿润室温环境下进行养护。需要说明的是,因疫情影响,至加载时试件的实际龄期约为6个月。所有试件的试验工作都在13 d内完成。

    采用意大利MATEST 4000 kN伺服控制压力试验机对MFRC圆柱体进行轴压强度测试(图4(a)),具体操作依照ASTM C39规范执行[34]。试件加载前,采用高强石膏对圆柱体试件上下两端表面进行找平,确保荷载能够均匀地施加于试件上下端。试验全程采用位移控制,加载速率设为0.18 mm/min,当轴压荷载降至峰值荷载的20%时,加载结束。

    图  4  试验装置
    Figure  4.  Test setups

    在圆柱试件半高处设置有两个标距为80 mm的环向应变片,以及两个标距为100 mm的纵向应变片,纵向应变片主要用于监测加载初始阶段试件的同心度。采用位移计测量试件的轴向变形,通过套箍将两个位移计对称设置于圆柱体试件的中心高度120 mm(即位移计的标距=120 mm)区域上,其平均读数除以标距为轴向应变。需要说明的是,当荷载达到峰值后,试件表面裂纹逐渐发展,套箍可能发生移位,致使位移传感器数值不能正确反映轴向应变。因此,上述位移计读数仅用于测量峰值荷载前的轴向应变。在加载平台上另外设置两个位移传感器以监测试件全高的轴向变形。峰值荷载后的实际轴向应变由以下方法求得:实际轴向应变=(全高位移计平均读数−峰值荷载处的全高位移计平均读数)/试件全高+峰值荷载处的套箍位移计平均读数/位移计标距。依据试验所得圆柱体轴压应力-应变全曲线(包含上升和下降段),依照ASTM C469 规程[35]获取弹性模量(E)、峰值应变(εco)及其他力学参数。

    依照ASTM C496 实验规程,采用前述相同压力试验机对MFRC圆柱体进行劈拉强度测试(图4(b))[36]。为保证劈拉力沿试件长度方向均匀分布,将尺寸为300 mm×2.5 mm×3 mm的木垫片粘贴于圆柱试件上下承压表面,并与中心线对齐;其后,将圆柱试件置于劈拉夹具内,加载条沿试件中心线对齐。试验全程采用劈拉应力控制,以1.2 MPa/min的速率加载至试件破坏。劈拉强度fsp由下式计算:

    fsp=2Pπld (1)

    其中:P为劈拉峰值荷载;ld分别为圆柱体试件的长度和直径。

    粗纤维掺入降低了混凝土拌合物的流动性,对采用配合比一的混凝土,坍落度由不含纤维时的220 mm降低至2vol%纤维掺量时的140 mm(图5(a))。原因在于:(1)各组分总表面积增加,需水量上升;(2)粗纤维的“棚架”效应。粗纤维的厚度和长度,以及不同配比对坍落度的影响如图5(b)所示。2-1组和2-2组的坍落度分别为165 mm和175 mm。这表明,在其他参数不变的情况下,随着纤维厚度的增加,坍落度略有提升,这是由于粗纤维的根数降低。粗纤维掺量为1.5vol%时,掺入30 mm和60 mm粗纤维的混凝土(2-4组和2-3组)的塌落度均在160 mm左右。

    图  5  粗纤维掺入对混凝土坍落度的影响
    Figure  5.  Effect of macro fibres on concrete slump value

    采用配合比一的情况下,1-1组和1-4组的坍落度分别为220 mm和170 mm。而采用配合比二的情况下,2-6组和2-5组的坍落度分别为195 mm和30 mm。同样掺入1.5vol%的粗纤维,配合比二混凝土的坍落度下降幅度显著大于配合比一混凝土。这表明不同基体配比的粗纤维混凝土具有不同的流动性。值得注意的是,纤维的掺量较高(≥1.5vol%)且长度较长(如90 mm)时,可能出现明显的粗纤维“棚架”效应。此时,振捣不够均匀可能引起浇筑质量控制不佳,而振捣过久易引起混凝土泌水。

    受压试件典型破坏模态如图6所示。轴压试验过程中,随着荷载上升,试件开始出现表面裂纹。达到峰值荷载后,未掺纤维的控制试件发生部分块体剥落(图6(a));MFRC试件无大块体剥落,可清晰看到起桥接作用的粗纤维,试件裂而不散,整体性保持较好(图6(b))。

    图  6  粗纤维混凝土(MFRC)轴压试件典型破坏模态
    Figure  6.  Typical failure modes of macro fibre reinforced concrete (MFRC) cylinders under compression

    轴压实验数据如表4所示。受压试件的轴压强度如图7所示。在批次1的试件中(均为配合比一混凝土),纤维掺量为0vol%、0.5vol%和1vol%的试件的平均轴压强度分别为51.7 MPa、46.0和49.1 MPa,表明粗纤维的掺入对轴压强度产生了轻微不利影响。然而,当纤维含量增加到1.5vol%和2vol%时,轴压强度和弹性模量均有明显下降,该两组试件在浇筑中有轻微泌水现象。如前文所述,粗纤维的“棚架”效应会引起混凝土的浇筑质量问题,在纤维的掺量和长度较高的情况下,可采取添加外加剂等措施以保证浇筑质量。

    表  4  MFRC轴压实验结果
    Table  4.  Key results of MFRC compression tests
    No. Group fc,m/MPa CoV/% εco,m/% CoV/% Em/GPa CoV/% vm CoV/%
    1-1 S1W0.32L0V0%T0 51.7 2.1 0.223 13.7 35.5 5.9 0.18 6.2
    1-2 S1W0.32L90V0.5%T0.5 46.0 6.4 0.175 10.6 36.2 0.7 0.19 0.0
    1-3 S1W0.32L90V1%T0.5 49.1 4.7 0.183 3.8 35.4 3.6 0.20 5.0
    1-4 S1W0.32L90V1.5%T0.5 37.4 5.3 0.167 9.0 32.7 6.8 0.19 11.2
    1-5 S1W0.32L90V2%T0.5 38.0 7.3 0.194 6.2 28.6 6.7 0.17 3.5
    2-1 S2W0.32L90V1.5%T1.2 74.5 1.6 0.264 4.7 39.7 1.2 0.17 0.0
    2-2 S2W0.32L90V1.5%T2.2 65.5 1.2 0.245 11.9 39.5 1.9 0.21 4.8
    2-3 S2W0.32L60V1.5%T0.5 57.4 1.0 0.236 16.2 37.1 5.4 0.19 6.0
    2-4 S2W0.32L30V1.5%T0.5 60.2 2.7 0.208 1.7 40.2 2.5 0.19 5.3
    2-5 S2W0.47L90V1.5%T0.5 50.9 2.9 0.272 0.7 32.4 4.6 0.19 16.6
    2-6 S2W0.47L0V0%T0 50.1 2.4 0.218 7.2 36.0 2.9 0.22 16.2
    Notes: fc,m—Average compressive strength; εco,m—Average value of strain at compressive strength; Em—Average value of modulus of elasticity in compression; vm—Average Poisson's ratio; CoV—Coefficient of variation.
    下载: 导出CSV 
    | 显示表格
    图  7  MFRC轴压强度
    Figure  7.  Compressive strength of MFRC

    在第2批次的试件中,纤维长度从30 mm提升至60 mm时,混凝土轴压强度从60.2 MPa降至57.4 MPa。纤维厚度从1.2 mm提升至2.2 mm时,混凝土轴压强度从74.5 MPa降至65.5 MPa。此外,在采用配合比二的情况下,粗纤维掺入对轴压强度没有显著影响。相关研究表明[25, 32],粗纤维掺入的主要影响有:(1)吸收部分自由水,致使水灰比轻微下降;(2)限制裂纹的发展;(3)产生粗纤维和水泥基体之间的弱界面。当前两种效应占主导地位时,混凝土强度受到有利影响。

    表4所示,批次2中C60配比的抗压强度和模量接近预期,而批次1中C60配比的抗压性能略低于预期,导致其与C40配比的数据相近。造成批次1中C60配比的抗压性能低于预期的原因有:(1)对于同批次试件的浇筑,所用粗细骨料均为同批次采购且同时进行清洗晾晒处理,而不同批次间所采用的原料可能有细微不同,进而造成了试件抗压性能的差异。因此,本文仅在同一批次试件中进行单因素影响分析。例如,采用批次2中的2-1组(粗纤维厚度为1.2 mm)和2-2组(粗纤维厚度为2.2 mm)进行粗纤维厚度的影响分析,而批次1中的1-4组(粗纤维厚度为0.5 mm)不包括在内;(2)如前文所述,粗纤维对混凝土强度有不利影响,尤其是纤维掺量较高(如≥1.5vol%)且长度较长(如90 mm)时可能出现明显的粗纤维“结团”效应,进而影响该组混凝土的浇筑质量(如1-4和1-5组)。

    圆柱体试件轴压应力-应变曲线如图8所示。粗纤维掺入对峰值前曲线的斜率(即弹性模量)影响较小,而对峰值后曲线形状产生较为明显的影响。如图8(a)所示,随着纤维掺量的上升(0vol%~1.5vol%),曲线峰值点有所下降,同时峰值后曲线形状逐渐变得平缓,应力下降速率降低。但纤维掺量提升至2vol%时,峰值后曲线的下降速率并未进一步降低。如前文所述,粗纤维掺入能有效桥接裂缝以限制其发展,因而试件能更持久地维持高轴压应力状态,表现为更平缓的峰值后曲线。

    图  8  MFRC轴压应力-应变曲线
    Figure  8.  Compressive stress-strain curves of MFRC

    粗纤维厚度对曲线的影响如图8(b)所示。纤维厚度从1.2 mm (2-1组)提升至2.2 mm (2-2组)时,曲线峰值和抗压延性均有所降低。本研究所用的机械加工工艺所能产生的粗纤维最薄为0.5 mm左右,厚度超过1.5 mm时粗纤维自身的层间脱离现象较为明显。因此,本研究认为粗纤维厚度的合理范围为0.5~1.5 mm。

    粗纤维长度对曲线的影响如图8(c)所示。粗纤维越长,曲线下降段越平缓,这与采用其他种类纤维的混凝土特性类似[37]。对于未掺入粗纤维的两种不同配比混凝土(图8(a)中S1W0.32L0V0T0组和图8(d)中S2W0.47L0V0T0组),峰值后应力下降至峰值的20%时其轴压应变均在0.5%左右。而掺入1.5vol%粗纤维后,该值均超过1%,增韧效果明显。

    图9为劈拉试件的两种典型破坏模态。未掺粗纤维的控制试件均于峰值荷载处沿中心线劈裂成两半。对于MFRC试件,达到峰值荷载后其加载面出现可视裂缝,荷载突然下降;由于有粗纤维横跨裂缝,MFRC试件此时并未被完全劈开,裂而不断。继续加载可使试件最终被劈裂成两半,在劈裂面可以看到许多横跨主裂缝分布的粗纤维。

    图  9  MFRC劈拉试件典型破坏模态
    Figure  9.  Typical failure modes of MFRC splitting tensile test specimens

    圆柱体劈拉试验结果见图10。粗纤维掺量从0vol%提升至2vol%时,劈拉强度从4.2 MPa提升至5.2 MPa。如图10(b)所示,粗纤维厚度从1.2 mm提升至2.2 mm时,劈拉强度的变化不明显。粗纤维长度从30 mm提升至60 mm时,劈拉强度从4.8 MPa提升至5.4 MPa。这与采用其他纤维的混凝土特性类似,即长细比越大,强度越高。此外,在采用配合比二的情况下,2-6组和2-5组的劈拉强度分别为3.5 MPa和4.9 MPa。这表明,同样掺入1.5vol%的粗纤维,采用配合比二的混凝土劈拉强度提升幅度(40%)显著大于配合比一(10%)。

    图  10  MFRC劈拉强度
    Figure  10.  Splitting tensile strength of MFRC

    纤维的掺入对混凝土拉伸性能影响显著。影响纤维混凝土(Fibre reinforced concrete,FRC)劈拉强度的因素包括基体强度,纤维体积掺量[38-39]、纤维直径和长度[40-47]。现有的众多纤维混凝土劈拉强度模型[38-47](表5)通常包括基体贡献项和纤维贡献项,两项结合方式或为乘积形式[40, 42, 44, 47]或为相加形式[38-39, 41, 43, 45-46]

    表  5  纤维混凝土(FRC)劈拉强度预测模型
    Table  5.  Models for predicting the splitting tensile strength of fibre reinforced concrete (FRC)
    Parameter Equation Note Ref.
    fc,Vf fspf=0.58(fc)0.5+3.02Vf fc [38]
    f_{\rm{c}}^{\prime},{V_{\rm{f}}},L,{d_{\rm{f}}} {f_{\rm{spf}}} = 0.83{(f_{\rm{c}}^{\prime})^{0.47}}\left(\dfrac{{ - 2}}{{1 + {{(0.125 R_{\mathrm{i}})}^{0.8}}}} + 3\right) \begin{gathered} R_{\mathrm{i}}{\text{ = }}{V_{\rm{f}}}L/{d_{\rm{f}}} \\ f_{\rm{c}}^{\prime} \leqslant 150\;{\text{MPa}} \\ \end{gathered} [40]
    f_{\rm{c}}^{},{V_{\rm{f}}} f_{\rm{spf}}=0.63(f_{\mathrm{c}})^{0.5}+3.01V_{\rm{f}}-0.02V_{\rm{f}}^2 {f_{\rm{c}}} \leqslant 100\;{\text{MPa}} [39]
    f_{\rm{c}}^{},{V_{\rm{f}}},L,{d_{\rm{f}}},{b_{\rm{f}} } {f_{\rm{spf}}} = 0.47{(f_{\rm{c}}^{})^{0.5}} + 4.2 F \begin{gathered} F{\text{ = }}{V_{\rm{f}}}{{\text{b}}_{\rm{f}}}L/{d_{\rm{f}}} \\ 41 \leqslant f_{\rm{c}}^{} \leqslant 115\;{\text{MPa}} \\ \end{gathered} [41]
    f_{\rm{cf} }^{},{V_{\rm{f}} },L,{d_{\rm{f}} } {f_{\rm{spf} }} = {({f_{\rm{cf} }})^{0.5}}\left(0.6 + 0.4{V_{\rm{f}}}\dfrac{L}{{{d_{\rm{f}} }}}\right) 4 \leqslant {f_{\rm{cf} }} \leqslant 120\;{\text{MPa}} [42]
    f_{\rm{cf} }^{},{V_{\rm{f}} },L,{d_{\rm{f}} },{b_{\rm{f}}} {f_{\rm{spf} }} = 0.095\left[{f_{\rm{cf} }} + 10\sqrt F \left(3 - \dfrac{{20}}{{{f_{\rm{cf} }}}}\right)\right] F{\text{ = }}{V_{\rm{f}} }{{\text{b}}_{\rm{f}}}L/{d_{\rm{f}}} [43]
    f_{\rm{cf} }^{},{V_{\rm{f}}},L,{d_{\rm{f}}},{b_{\rm{f}}} {f_{\rm{spf}}} = {({f_{\rm{cf} }})^{0.5}}\left(0.94{V_{\rm{f}} }\dfrac{L}{{{d_{\rm{f}}}}}{b_{\rm{f}} } + 0.67\right) 120 \leqslant {f_{\rm{cf} }} \leqslant 200\;{\text{MPa}} [44]
    Notes: fc'—Characteristic compressive strength of plain concrete; fc— Measured compressive strength of plain concrete; fcf—Measured compressive strength of FRC; Vf—Fibre volume fraction; L—Fibre length; df—Fibre diameter; bf—Bond factor; fspf—Splitting tensile strength of FRC; Ri—Reinforcing index; F—Fibre factor.
    下载: 导出CSV 
    | 显示表格
    {f}_{{\mathrm{sp}}}=U({f}'_{{\mathrm{c}}}) G(R_{\mathrm{i}}) (2)

    f_{\mathrm{sp}}=U(f_{\mathrm{c}}')+G(R_{\mathrm{i}}) (3)

    其中:fsp为混凝土劈拉强度;fc'为混凝土轴压强度标准值;U(fc')为基体贡献项,常采用以fc'为自变量的幂函数,JSCE和FIB规范采用的指数为0.67[48-49],而ACI规范采用的指数为0.5[50-51]G(Ri)为纤维贡献项,是以Ri为自变量的函数,多采用幂函数或多项式函数;参数Ri (Reinforcing index)为纤维影响系数,可表示为

    {R_{\text{i}}} = \frac{{V_{{\mathrm{f}}} }L} {d_{{\mathrm{f}}} } (4)

    其中:Vf为纤维体积掺量;L为纤维长度;df为纤维直径。在预测纤维混凝土强度、弹性模量、泊松比以及坍落度时,采用型如Vf(L/df)αndLα/dfβ (nd为单位体积纤维数量,aβ均为拟合系数)的纤维影响系数可以获得较好的拟合效果[52]

    采用预测值/实验值之比的统计参数和积分绝对误差(Integral absolute error,IAE)评估了现有FRC劈拉强度预测公式对MFRC劈拉强度的预测效果,其统计参数见表6。IAE由下式计算:

    表  6  FRC模型对MFRC劈拉强度的预测效果
    Table  6.  Evaluation of FRC models for predicting the splitting tensile strength of MFRC
    Ratio of predicted value and test value IAE/% Ref.
    Average SD CoV/%
    0.83 0.14 16 19.8 [38]
    1.32 0.14 10 30.5 [40]
    0.98 0.16 16 12.6 [39]
    0.92 0.13 15 13.5 [41]
    1.19 0.12 10 18.5 [42]
    1.28 0.16 13 27.5 [43]
    1.34 0.13 10 33.5 [44]
    1.00 0.08 8 6.0 Eq. (5)
    1.02 0.15 15 7.4 Eq. (6)
    Notes: SD—Standard deviation; IAE—Integral absolute error.
    下载: 导出CSV 
    | 显示表格
    \mathrm{IAE}=\sum_{ }^{ }\frac{\left| {{O_{\rm{i}}} - {P_{\rm{i}}}} \right|}{\displaystyle\sum_{ }^{ }O_{\mathrm{i}}} (5)

    其中:Oi为实验值;Pi为预测值。IAE小于10%,代表模型预测效果较好[40]。如表6所示,各模型预测效果不一,其预测值与实验值之比的平均值在0.83~1.34,IAE在12.6%~33.5%之间。

    参考现有FRC模型的特征和基本形式,经过试算及优选,建议以下MFRC劈拉强度预测模型:

    {f_{{\mathrm{sp}}}} = 0.53{(f_{\mathrm{c}}^{})^{0.5}} + 8.58{V_{\mathrm{f}}}{L^{0.51}}{\left(\frac{\psi }{A}\right)^{0.056}} (6)

    其中,ψA分别为粗纤维的截面周长及截面积,因此该模型适用于任意截面形状的粗纤维。考虑到式(5)较为复杂,经过试算也可以采用以下较为简单的表达式:

    f_{\mathrm{sp}}=(f_{\mathrm{c}})^{0.67}\left(\text{0}\text{.28}+0.021V_{\mathrm{f}}L\frac{\psi}{A}\right) (7)

    当粗纤维掺量为0vol%时,公式(6)退化为普通混凝土劈拉强度预测公式:

    {f_{{\mathrm{sp}}}} = {\text{0}}{\text{.28}}{(f_{\mathrm{c}}^{})^{0.67}} (8)

    其形式与JSCE[48]和FIB[49]规范中普通混凝土劈拉强度预测公式相同:

    f_{\mathrm{sp}}=\text{0}\text{.23}(f_{\mathrm{c}}^{ }\mathrm{'})^{0.67},\begin{array}{*{20}{c}} & 20\; \text{MPa}\leqslant f_{\mathrm{c}}^{ }\mathrm{'}\leqslant\end{array}80\; \text{MPa} (9)
    f_{\mathrm{sp}}=\text{0}\text{.3}(f_{\mathrm{c}}^{ }\mathrm{'})^{0.67},\begin{array}{*{20}{c}} & f_{\mathrm{c}}^{ }\mathrm{'}\leqslant\end{array}50\; \text{MPa} (10)

    式(6)中的拟合系数0.28 处于JSCE和FIB规范公式中的拟合系数0.23~0.3之间。

    本文所提式(5)和式(6)的预测效果见图11。式(5)的R2为0.75,IAE为5.96%,式(6)的R2为0.70,IAE为7.4%,总体预测效果较好。值得注意的是,由于砂石等原材料的性能差异,部分试件的抗压性能低于预期。然而,这些实验数据仍被纳入数据库和预测公式推导中。一方面,尽管抗压强度数据存在误差,劈裂强度数据仍然具有价值。另一方面,在实际应用中,砂石骨料来源不同所带来的性能差异也应在数据库中得到体现。在进一步研究中,建议扩大实验规模以研究更广泛的实验参数范围对数据结果的影响。同时,预测公式中应考虑除抗压强度外的更多指标。

    图  11  MFRC劈拉强度实验值和预测值的比较
    Figure  11.  Predicted vs. test values of MFRC splitting tensile strength
    ψ—Section circumference of macro fibre; A—Sectional area of macro fibre; R2—Coefficient of determination

    (1)添加粗纤维降低了混凝土拌合物的流动性,坍落度随粗纤维掺量提升而降低。例如,在配比一混凝土中掺入2vol%体积掺量的粗纤维后,混凝土坍落度从220 mm降至140 mm。纤维掺量较高且长度较长时,其“棚架”效应可能会影响混凝土的浇筑质量。

    (2)粗纤维的掺入对轴压强度有轻微不利影响,粗纤维过厚会造成混凝土强度明显降低。例如,在配比一混凝土中,纤维厚度从1.2 mm提升至2.2 mm时,混凝土轴压强度从74.5 MPa降至65.5 MPa。

    (3)粗纤维的掺入提升了混凝土的劈拉强度。在配比二混凝土中,粗纤维掺量从0vol%提升至1.5vol%时,劈拉强度由3.5 MPa提升至4.9 MPa,提升了40%。

    (4)对于本文的粗纤维混凝土试验数据,所提出的劈拉强度公式的积分绝对误差(IAE)为5.96%~7.40%,总体预测效果较好。

    致谢:感谢广东工业大学提供的场地和试验设备。

  • 图  1   P-Co-C (a)、P-Fe-C (b),P-Ni-C (c)、P-Mn-C (d)的SEM图像

    Figure  1.   SEM images of P-Co-C (a), P-Fe-C (b), P-Ni-C (c) and P-Mn-C (d)

    图  2   P-Co-C、P-Fe-C、P-Ni-C和P-Mn-C样品的XRD图谱(a)和Raman图谱(b)

    Figure  2.   XRD patterns (a) and Raman spectra (b) of P-Co-C, P-Fe-C, P-Ni-C and P-Mn-C

    图  3   P-Co-C、P-Ni-C、P-Fe-C和P-Mn-C样品的N2吸脱附曲线(a)和相对应孔径分布(b)

    Figure  3.   N2 adsorption and desorption isotherm (a) and the corresponding pore size distribution curves (b) of P-Co-C, P-Ni-C, P-Fe-C and P-Mn-C

    图  4   Co2p(a)、Fe2p(b)、Ni2p(c)和Mn2p(d)高分辨XPS图谱

    Sat.—Satellite peak

    Figure  4.   High-resolution XPS spectra of Co2p (a), Fe2p (b), Ni2p (c) and Mn2p (d)

    图  5   (a) P-Co-C、P-Ni-C、P-Fe-C和P-Mn-C在N2饱和0.1 mol/L KOH及扫速10 mV·s−1下的循环伏安曲线;(b) O2饱和0.1 mol/L KOH 及转速1600 r/min条件下测得的线性扫描伏安曲线;(c)转移电子数n;(d)中间产物

    Figure  5.   (a) CVs of P-Co-C, P-Ni-C, P-Fe-C and P-Mn-C in N2-saturated 0.1 mol/L KOH with a scan rate of 10 mV·s−1; (b) LSVs of P-Co-C, P-Ni-C, P-Fe-C and P-Mn-C in O2-saturated 0.1 mol/L KOH with a scan rate of 10 mV·s−1 at rotating speed of 1600 r/min; (c) Electron transfer number n; (d) Determined peroxide percentage

    图  6   P-Co-C、P-Fe-C、P-Ni-C、P-Mn-C、P-C 及Pt/C的塔菲尔曲线

    Figure  6.   Tafel plots of P-Co-C, P-Fe-C, P-Ni-C, P-Mn-C, P-C and Pt/C

    图  7   P-Co-C、P-Ni-C、P-Fe-C、P-Mn-C和商业化Pt/C催化剂氧气还原反应计时电流响应曲线(a)及抗甲醇毒化性能曲线(b)(O2饱和0.1 mol/L KOH、−0.3 V恒电压、扫速10 mV·s−1、转速为1600 r/min)

    Figure  7.   Current-time (i-t) chronoamperometric responses (a) for the ORR of P-Co-C, P-Ni-C, P-Fe-C, P-Mn-C and commercial Pt/C and chronoamperometric response (b) upon addition of 3 mol/L methanol into O2-saturated 0.1 mol/L KOH at −0.3 V with a scan rate of 10 mV·s−1 and a rotating speed of 1600 r/min

  • [1]

    CHIU Y S P, LIAN J H, CHIU V, et al. Mathematical modeling for a multiproduct manufacturing system featuring postponement, external suppliers, overtime, and scrap[J]. International Journal of Industrial Engineering Computations,2020,13:1-12.

    [2]

    MA M J, YANG X X, QIAO J S, et al. Progress and challenges of carbon-fueled solid oxide fuel cells anode[J]. Journal of Energy Chemistry,2021,56:209-222. DOI: 10.1016/j.jechem.2020.08.013

    [3]

    LYU D, YAO S X, BAHARI Y, et al. In situ molecular-level synthesis of N, S co-doped carbon as efficient metal-free oxygen redox electrocatalysts for rechargeable Zn-air batteries[J]. Applied Materials Today,2020,20:100737-100747. DOI: 10.1016/j.apmt.2020.100737

    [4]

    LIU T, ZHANG L M, TIAN Y. Earthworm-like N, S-doped carbon tube-encapsulated Co9S8 nanocomposites derived from nanoscaled metal-organic frameworks for highly efficient bifunctional oxygen catalysis[J]. Journal of Materials Chemistry A,2018,6:5935-5943. DOI: 10.1039/C7TA11122C

    [5]

    AMIJNU I S, PU Z H, LIU X B, et al. Multifunctional Mo-N/C@MoS2 electro-catalysts for HER, OER, ORR, and Zn-air batteries[J]. Advanced Functional Material,2017,27:1702300-1702311. DOI: 10.1002/adfm.201702300

    [6]

    WEI J, HU Y X, LIANG Y, et al. Nitrogen-doped nanoporous carbon/graphene nano-sandwiches: Synthesis and application for efficient oxygen reduction[J]. Advanced Functional Material,2015,25:5768-5777. DOI: 10.1002/adfm.201502311

    [7]

    LI J, QIN X P, HOU P X, et al. Identification of active sites in nitrogen and sulfur co-doped carbon-based oxygen reduction catalysts[J]. Carbon,2019,147:303-311. DOI: 10.1016/j.carbon.2019.01.018

    [8]

    ZHAO Y, KAMYIYA K, HASHIMOTO K, et al. Efficient bifunctional Fe/C/N electrocatalysts for oxygen reduction and evolution reaction[J]. Journal of Physical Chemistry C,2015,119:2583-2588. DOI: 10.1021/jp511515q

    [9]

    QIAN Y H, HU Z G, GE X M. A metal-free ORR/OER bifunctional electrocatalyst derived from metal-organic frameworks for rechargeable Zn-air batteries[J]. Carbon,2017,111:641-650. DOI: 10.1016/j.carbon.2016.10.046

    [10]

    ZENG K, SU J M, CAO X C, et al. B, N co-doped ordered mesoporous carbon with enhanced electrocatalytic activity for the oxygen reduction reaction[J]. Journal of Alloys and Compounds,2020,824:153908-153914. DOI: 10.1016/j.jallcom.2020.153908

    [11]

    ZHANG D, SUN P P, ZUO Z, et al. N, P-co-doped carbon nanotubes coupled with Co2p nanoparticles as bifunctional oxygen electrocatalyst[J]. Journal of Electroanalytical Chemistry,2020,871:114327-114335. DOI: 10.1016/j.jelechem.2020.114327

    [12]

    HUANG N, YANG L, ZHANG M Y, et al. Cobalt-embedded N-doped carbon arrays derived in situ as trifunctional catalyst toward hydrogen and oxygen evolution, and oxygen reduction[J]. ChemElectroChem,2019,6:4522-4532. DOI: 10.1002/celc.201901106

    [13]

    QUÍLEZ-BERMEJ J, GONZALEZ-GAITAN C, EMILIA M, et al. Effect of carbonization conditions of polyaniline on its catalytic activity towards ORR-Some insights about the nature of the active sites[J]. Carbon,2017,119:62-71. DOI: 10.1016/j.carbon.2017.04.015

    [14]

    SON S, LIM D, NAM D, et al. N, S-doped nanocarbon derived from ZIF-8 as a highly efficient and durable electro-catalyst for oxygen reduction reaction[J]. Journal of Solid State Chemistry,2019,274:237-242. DOI: 10.1016/j.jssc.2019.03.036

    [15]

    ZHANG W M, YAO X Y, ZHOU S N, et al. ZIF-8/ZIF-67-derived co-Nx-embedded 1D porous carbon nanofibers with graphitic carbon-encased Co nanoparticles as an efficient bifunctional electrocatalyst[J]. Small,2018,14:1800423-1800430. DOI: 10.1002/smll.201800423

    [16]

    YUAN H L, WANG Y Q, ZHOU S M. Low-temperature preparation of superparamagnetic CoFe2O4 microspheres with high saturation magnetization[J]. Nanoscale Research Letters,2010,5:1817-1821. DOI: 10.1007/s11671-010-9718-7

    [17]

    LI J S, LI S L, TANG Y J, et al. Nitrogen-doped Fe/Fe3C@graphitic layer/carbon nanotube hybrids derived from MOFs: Efficient bifunctional electrocatalysts for ORR and OER[J]. Chemical Communications,2015,51:2710-2713. DOI: 10.1039/C4CC09062D

    [18]

    SHARMA Y, SHARMA N, RAO G V S, et al. Studies on spinel cobaltites, FeCo2O4 and MgCo2O4 as anodes for Li-ion batteries[J]. Solid State Ionics,2008,179:587-597. DOI: 10.1016/j.ssi.2008.04.007

    [19]

    OKU M, HIROKAWA K. X-ray photoelectron spectra of inequivalent atoms in inorganic compounds[J]. Journal of Solid State Chemistry,1979,30:45-53. DOI: 10.1016/0022-4596(79)90128-2

    [20]

    XIAO J W, YANG S. Bio-inspired synthesis of NaCl-type CoxNi1-xO (0 ≤ x < 1) nanorods on reduced graphene oxide sheets and screening for asymmetric electrochemical capacitors[J]. Journal of Material Chemistry,2012,22:12253-12262. DOI: 10.1039/c2jm31057k

    [21]

    FU Y Y, XU L, ZHAO W K, et al. Spinel CoMn2O4 nanosheet arrays grown on nickel foam for high-performance supercapacitor electrode[J]. Applied Surface Science,2015,357:2013-2021. DOI: 10.1016/j.apsusc.2015.09.176

    [22]

    MA S C, SUN L Q, CONG L N. Multiporous MnCo2O4 microspheres as an efficient bifunctional catalyst for nonaqueous Li-O2 batteries[J]. Journal of Physical Chemistry C,2013,117:25890-25897. DOI: 10.1021/jp407576q

    [23]

    WU J, YANG Z R, SUN Q J, et al. Synthesis and electrocatalytic activity of phosphorus-doped carbon xerogel for oxygen reduction[J]. Electrochimica Acta,2014,127:53-60. DOI: 10.1016/j.electacta.2014.02.016

  • 期刊类型引用(3)

    1. 范浩博,豆书亮,李垚. 二氧化钒智能热控涂层光学结构原理及研究进展. 材料导报. 2025(01): 34-43 . 百度学术
    2. 张新宇,徐慧妍,杨立凯,王尚,杨帅军,蒋绪川. 智能窗用二氧化钒基复合结构薄膜的制备及研究进展. 复合材料学报. 2025(02): 679-691 . 本站查看
    3. 张丽萍,孟晓荣,宋锦峰,杜金晶. VO_2@KH550/570@PS复合薄膜的制备及其热致相变性能. 化工学报. 2024(09): 3348-3359 . 百度学术

    其他类型引用(7)

图(7)
计量
  • 文章访问数:  1013
  • HTML全文浏览量:  552
  • PDF下载量:  55
  • 被引次数: 10
出版历程
  • 收稿日期:  2022-01-05
  • 修回日期:  2022-02-07
  • 录用日期:  2022-02-11
  • 网络出版日期:  2022-02-28
  • 刊出日期:  2023-02-14

目录

/

返回文章
返回