高级检索

环缝气流冲击雾化制粉喷嘴流场特性的仿真研究

马海英, 张鹏, 郭志军

马海英, 张鹏, 郭志军. 环缝气流冲击雾化制粉喷嘴流场特性的仿真研究[J]. 粉末冶金技术, 2018, 36(2): 148-152. DOI: 10.19591/j.cnki.cn11-1974/tf.2018.02.012
引用本文: 马海英, 张鹏, 郭志军. 环缝气流冲击雾化制粉喷嘴流场特性的仿真研究[J]. 粉末冶金技术, 2018, 36(2): 148-152. DOI: 10.19591/j.cnki.cn11-1974/tf.2018.02.012
MA Hai-ying, ZHANG Peng, GUO Zhi-jun. Simulation research on flow field characteristics of atomizing nozzle by circular seam air impact[J]. Powder Metallurgy Technology, 2018, 36(2): 148-152. DOI: 10.19591/j.cnki.cn11-1974/tf.2018.02.012
Citation: MA Hai-ying, ZHANG Peng, GUO Zhi-jun. Simulation research on flow field characteristics of atomizing nozzle by circular seam air impact[J]. Powder Metallurgy Technology, 2018, 36(2): 148-152. DOI: 10.19591/j.cnki.cn11-1974/tf.2018.02.012

环缝气流冲击雾化制粉喷嘴流场特性的仿真研究

基金项目: 

国家自然科学基金资助项目 51175150

详细信息
    通讯作者:

    马海英, E-mail:aotemen2012@126.com

  • 中图分类号: TD403

Simulation research on flow field characteristics of atomizing nozzle by circular seam air impact

More Information
  • 摘要: 分析了一种环缝气流冲击雾化制粉喷嘴的结构和工作原理, 并对该雾化喷嘴模型进行了流场分布数值仿真, 研究了不同进气压力和窄缝倾角对流场速度分布的影响。研究结果表明, 从进气管入口到气流冲击雾化喷嘴出口, 加速气流速度保持最大值不变, 在距进气管入口中心线上80 mm位置之后逐渐减小; 随着进气压力的增加, 气流冲击雾化喷嘴出口速度表现出明显的增加趋势, 且仿真结果与理论计算结果误差在10%范围以内。随着窄缝倾角的增加, 负压涡流逐步向进粉管靠近, 回流区、分离区和混合区的最大速度值表现出单调增加的变化规律; 综合考虑气流冲击效果和喷射速度, 确定最优的窄缝倾角50°~60°。
    Abstract: The structure and working principle of atomizing nozzle by circular seam air impact were analyzed, the numerical simulation model on flow field distribution of atomizing nozzle was studied, and the effects of inlet pressure and circular seam inclination on velocity distribution in flow field were investigated. The results show that, the acceleration of air velocity remains unchanged from the inlet of intake pipe to the exit of air impact atomizing nozzle, and the maximum is achieved at 80 mm away from the center line of intake pipe inlet. With the increasing of inlet pressure, the exit velocity of air impact atomizing nozzle shows an increasing trend obviously, and the error between simulation results and theoretical calculation result is within 10%. With the increasing of circular seam inclination, the negative pressure eddy current is gradually closed to the import of powder tube, the maximum velocity in backflow zone, separation zone, and mixed zone increases monotonously. Considering the air impact effect and jet velocity, the optimal circular seam inclination is determined as 50°~60°.
  • 在切削工具领域, 自从上世纪50年代WC潜在替代者TiC出现以来, Ti(C, N)基金属陶瓷已引起广大科研工作者的高度关注[1]。Ti(C, N)基金属陶瓷具有较高的红硬性和耐磨性、优良的抗月牙洼磨损能力和化学稳定性、极低的摩擦系数, 是较为理想的切削刀具材料[2-5]。不过, 由于Ti(C, N)与Ni或Co的润湿性较差, 直接由Ti(C, N)与Ni或Co制备的Ti(C, N)基金属陶瓷的韧性和强度较低, 其抗弯强度仅为相应WC–Co硬质合金的70%左右[6-8], 这不仅不利于Ti(C, N)基金属陶瓷特有优势的发挥, 也大大限制了Ti(C, N)基金属陶瓷在机加工领域的广泛应用。

    为改善Ti(C, N)对Ni或Co的润湿性, 进而提高Ti(C, N)基金属陶瓷的力学性能和使用性能, 通常加入WC或Mo2C等第二类碳化物。大量研究表明, 在Ti(C, N)基金属陶瓷制备过程中添加适量的WC, 不仅可有效改善其润湿性和烧结性能、提高密度、减少金属陶瓷的孔隙率, 而且可提高韧性、减小晶粒长大速度。不过, 含WC的Ti(C, N)基金属陶瓷其微观组织结构非常复杂, 虽然一般均呈现出“芯/环” (Core/Rim)结构, 但WC粒度和添加量不同以及Ti(C, N)种类和金属陶瓷基体成分不同, 其微观组织结构也非常不同, 例如, 芯部未熔Ti(C, N)硬质相的大小和形状及环部固溶体的成份、厚度和形貌等, 从而导致其性能变化也不尽相同[9-11]

    本文拟探讨WC添加量(质量分数)对以Ti(C0.7N0.3)为基体、Ni和Co为黏结相、Mo2C及TaC为增强相的多组元复杂成分体系Ti(C0.7N0.3)基金属陶瓷微观组织结构及性能的影响, 以期为这类高性能Ti(C, N)基金属陶瓷的研究与应用提供技术参考。

    以Ti(C0.7N0.3)为基体, WC、Mo2C、TaC为增强相, Ni和Co为黏结相, 实验用Ti(C0.7N0.3)基金属陶瓷组元成分如表 1所示, 其中主要原料的化学成分及平均粒径如表 2所示。

    表  1  Ti(C0.7N0.3)基金属陶瓷组元成分(质量分数)
    Table  1.  Component of Ti(C0.7N0.3) -based cermets  %
    原料 Ti(C0.7N0.3) WC Mo2C TaC Co+Ni
    a 65 0 10 5 20
    b 60 5 10 5 20
    c 55 10 10 5 20
    d 50 15 10 5 20
    e 45 20 10 5 20
    下载: 导出CSV 
    | 显示表格
    表  2  金属陶瓷主要原料的化学成分和粒度
    Table  2.  Chemical component and particle size of raw powders
    原料 化学成分(质量分数)/ % 粒度/ μm
    Ti Ta Mo W C N O Ni Co
    Ti(C0.7N0.3) 77.87 13.56 6.73 0.61 1.82
    WC 93.23 6.16 0.16 0.42
    Mo2C 93.65 5.76 0.05 0.58 1.62
    TaC 93.41 6.23 0.05 0.12 1.14
    Ni 0.011 0.49 99.95 1.96
    Co 0.015 0.47 99.97 0.86
    下载: 导出CSV 
    | 显示表格

    将原料粉按组分配比称量装入球磨罐, 加入适量无水乙醇做为球磨介质, 采用硬质合金球(球料质量比7:1)在罐磨球磨机内进行湿磨混匀, 转速45 r·min-1, 球磨时间36 h, 之后在80℃下干燥后掺蜡以提高粉末成形性, 过筛后在100 MPa压力下压制成形。压坯经6 MPa氩气低压烧结制得Ti(C, N)基金属陶瓷, 烧结温度为1420℃, 保温时间60 min。将烧结后的压块经线切割并用金刚石砂轮磨削制得直径3 mm×4 mm×35 mm的标准试样条进行力学、磁学性能测试以及微观结构观察。

    利用排水法测量密度; 用三点弯曲法在力学性能试验机上测量抗弯强度, 跨距30 mm, 加载速度为0.5 mm·min-1, 测试结果为5个样品的平均值; 用HR-150A洛氏硬度计测量洛氏硬度(HRA), 取5个点的平均值; 用钴磁计和矫顽磁力计分别测定试样的钴磁和矫顽磁力; 将试样抛光后用扫描电子显微镜(scanning electron microscope, SEM)观察Ti(C0.7N0.3)基金属陶瓷的显微组织, 并对组织中的不同相区进行能谱(energy dispersive spectrum, EDS)分析。

    在扫描电子显微镜背散射电子模式下(scanning electron microscope-backscattered electron, SEM-BSE), 高原子序数的元素含量越多, SEM-BSE照片中显示的颜色就越白[12]。WC含量(质量分数)对Ti(C0.7N0.3)基金属陶瓷显微组织的影响如图 1所示, 从图 1可以看出, WC质量分数从0增大到20%的过程中, 显微组织都有明显的芯/壳结构。对不同相区进行能谱分析, 结果如图 2所示, 黑色芯部主要是由Ti、C和N组成的Ti(C, N), 灰色环部组织是由Mo2C、WC、TaC和Ti(C0.7N0.3)固溶形成的固溶体(Ti, Mo, Ta, W) (C, N), 白色相区主要由黏结剂Co和Ni及溶解在黏结相中起固溶强化作用的Ti、Mo等组成。

    图  1  添加不同WC质量分数Ti(C0.7N0.3)基金属陶瓷的背散射电子像: (a) 0; (b) 5%; (c) 10%; (d) 15%; (e) 20%
    Figure  1.  Backscattered electron images of Ti(C0.7N0.3) -based cermets added by WC in different mass fraction: (a) 0; (b) 5%; (c) 10%; (d) 15%; (e) 20%
    图  2  添加质量分数为15%WC制备的Ti(C0.7N0.3)基金属陶瓷不同相区的能谱分析: (a)扫描电子显微形貌; (b)位置a处能谱分析; (c)位置b处能谱分析; (d)位置c处能谱分析
    Figure  2.  EDS patterns of different regions in Ti(C0.7N0.3) -based cermets added by 15%WC by mass: (a) SEM image; (b) EDS in region a; (c) EDS in region b; (d) EDS in region c

    当WC质量分数为0时, 图 1 (a)显示完整的黑芯-灰环结构, 这是由于在烧结过程中, Mo2C溶解并在未溶解的Ti(C0.7N0.3)颗粒表面以(Ti, Mo) (C, N)的形式析出并长大, 形成灰色的环形相。随着WC质量分数增大, 黑色芯部的数量及体积比减小, 环形结构和黏结相的体积比增大, 表明WC促进了Ti(C0.7N0.3)在黏结相中溶解[13]。在液相烧结过程中, 环形相组成的变化主要受形成温度和系统中碳化物溶解速率影响, 体现在最终的微观结构中。在烧结过程中Ti(C0.7N0.3)部分溶解, 然后在降温过程中在未溶解的Ti(C0.7N0.3)颗粒表面以固溶形式析出。Ti(C0.7N0.3)在黏结相中溶解饱和时开始形成环形相, WC在黏结相中的溶解速率是Ti(C0.7N0.3)的2倍以上, 在液相烧结过程中, 相对不稳定的WC先溶解并沉淀在硬质相周围, 形成稳定的、高W质量分数的内环相[4, 14]。故随着WC质量分数增加, 颗粒粗化的趋势得到了抑制, 且环形结构的外环相变薄。当WC质量分数超过15%时, 由于WC对硬质相Ti(C0.7N0.3)溶解的促进作用且Ti(C0.7N0.3)的粒径较小, 在固相烧结阶段通过扩散形成(Ti, Mo) (C, N)固溶体, 很多部分无原始的硬质相Ti(C0.7N0.3)残留, 在SEM-BES模式下呈现出白芯结构, 如图 1 (e)所示。

    WC质量分数从0增大到20%的过程中, Ti(C0.7N0.3)基金属陶瓷的密度从6.0916 g·cm-3逐渐增大到7.1035 g·cm-3, 其抗弯强度及硬度随WC质量分数变化关系如图 3所示。WC质量分数从0增大到15%的过程中, 金属陶瓷的抗弯强度及硬度均提高。达到15%之后, 随着WC质量分数继续增加, 抗弯强度无明显变化, 硬度有所下降。

    图  3  WC质量分数对Ti(C0.7N0.3)基金属陶瓷硬度(a)和抗弯强度(b)的影响
    Figure  3.  Effect of WC mass fraction on the hardness (a) and bending strength (b) of Ti(C0.7N0.3) -based cermets

    随着WC质量分数增加, Ti(C0.7N0.3)基金属陶瓷的抗弯强度增加, 这是由于WC改善了Ti(C0.7N0.3)相的润湿性, 使得更多的Ti、Mo、W等元素溶解到黏结相中, 通过固溶强化增大了抗弯强度。同时, 随着WC质量分数增多, 金属陶瓷的密度也在增加, WC的溶解增大了内环相的体积比, 抑制了硬质相长大, 使得晶粒细化, 硬度提高。在WC质量分数达到15%以后, 溶解逐渐达到饱和, 环形相相对稳定, 因而抗弯强度不再随WC质量分数变化而变化。但是随着WC增多, 硬质相Ti(C0.7N0.3)的质量分数相对下降, 而Ti(C0.7N0.3)的显微硬度比WC高, 因而硬度有所下降。Mo2C同样具有改善润湿性, 细化晶粒的作用。综上, 当添加WC质量分数为15%、Mo2C质量分数为10%时, Ti(C0.7N0.3)基金属陶瓷的环形相厚度稳定, 晶粒细化和固溶强化效果最好, 因此此时综合力学性能最佳, HRA 91.6, 抗弯强度2160 MPa。

    Ti(C0.7N0.3)基金属陶瓷的钴磁及矫顽磁力随WC质量分数变化的关系如图 4所示, 随WC质量分数增加, 金属陶瓷的钴磁逐渐降低, 矫顽磁力逐渐增大。钴和镍可以形成无限互溶的固溶体, 硬质相中的W、Ti、Ta等原子在液相烧结过程中发生高温扩散和在液相中溶解析出, 合金烧结过程中, W、Ti、Ta、Mo等原子在钴镍固溶体中置换一定量的Co、Ni原子而形成置换固溶体。各原子在钴镍固溶体中的数量是不同的, 它们在钴镍固溶体中的数量、形态及分布决定了Ti(C0.7N0.3)基金属陶瓷的钴磁和矫顽磁力。

    图  4  WC质量分数对Ti(C0.7N0.3)基金属陶瓷钴磁(a)和矫顽磁力(b)的影响
    Figure  4.  Effect of WC mass fraction on the cobalt magnetization (a) and coercivity (b) of Ti(C0.7N0.3) -based cermets

    Ti(C0.7N0.3)基金属陶瓷中W、Ti、Ta、Mo等元素为顺磁或抗磁金属, 它们基本不显示磁性, 而Co和Ni为强磁性金属, 因此金属陶瓷的磁性来自具有磁性的钴镍黏结金属[15]。随着WC质量分数的增加, 黏结相中固溶的非磁性合金元素总质量分数增加, 其晶格膨胀程度也会增大, 使得黏结金属原子间的磁交换作用减弱, 从而导致金属陶瓷的钴磁降低。矫顽磁力表示对应于磁壁移动的阻抗量, 可以反映平均晶粒度, 但无法反映是否存在异常长大晶粒、是否存在游离碳等。在黏结相含量一定的情况下, 合金碳化物晶粒度越细, 则黏结相厚度越薄, 矫顽磁力就越大[16]。由于原料WC粉末费氏粒度最小, 而且WC的添加可防止晶界滑动和晶粒长大, 可有效抑制粗晶生长, 因此矫顽磁力随着WC质量分数增大而增大。综上所述, WC质量分数从0增大到20%的过程中, Ti(C0.7N0.3)基金属陶瓷的钴磁减小, 矫顽磁力增大。

    (1) 经1420℃氩气低压烧结1 h后, Ti(C0.7N0.3)基金属陶瓷具有典型的芯/壳结构。随着WC质量分数增大, 硬质相在黏结相中的溶解度增大, 其数量和体积减小, 环形相和黏结相体积变大, 晶粒度逐渐细化。

    (2) 随着WC质量分数从0增大到20%, Ti(C0.7N0.3)基金属陶瓷的密度从6.0916 g·cm-3逐渐增大到7.1035 g·cm-3; 硬度先增大后减小, 抗弯强度先增大后趋于稳定。WC质量分数为15%时, 力学性能最佳, HRA 91.6, 抗弯强度2160 MPa。

    (3) 随着WC质量分数增大, 饱和磁化强度逐渐降低, 矫顽磁力逐渐增大。

  • 图  1   环缝气流冲击雾化制粉喷嘴结构图

    Figure  1.   Structure of atomizing nozzle by circular seam air impact

    图  2   环缝气流冲击雾化喷嘴计算区域及边界类型

    Figure  2.   Calculation region and boundary of atomizing nozzle by circular seam air impact

    图  3   喷嘴内部流场速度分布(a)及速度矢量(b)

    Figure  3.   Velocity distribution in flow field of atomizing nozzle (a) and velocity vector (b)

    图  4   喷嘴中心线速度分布

    Figure  4.   Centerline velocity distribution of nozzle

    图  5   流场速度在不同窄缝倾角下的分布云图: (a) 30°; (b) 55°; (c) 80°

    Figure  5.   Velocity distribution in flow field of atomizing nozzle at different circular seam inclination: (a) 30°; (b) 55°; (c) 80°

    表  1   不同进气压力下喷嘴出口速度

    Table  1   Exit velocity of atomizing nozzle in different inlet pressure

    进气压力/ MPa 喷嘴出口速度/ (m·s-1)
    仿真结果 理论计算结果
    4 23.4 22.2
    6 28.6 26.8
    8 34.2 31.5
    10 36.8 33.6
    下载: 导出CSV

    表  2   雾化喷嘴不同区域最大速度随窄缝倾角的变化结果

    Table  2   Maximum velocity in different zone at different circular seam inclination of atomizing nozzle

    窄缝倾角/ (°) 最大速度/ (m·s-1)
    涡流区 回流区 分离区 混合区
    30 55.4 29.6 25.8 25.0
    35 55.2 31.2 27.9 26.9
    40 55.2 33.5 31.8 29.8
    45 55.7 34.9 33.1 30.1
    50 55.1 36.0 34.2 31.2
    55 55.3 37.1 35.0 33.3
    60 55.6 38.2 35.8 34.2
    65 54.8 41.8 36.6 35.8
    70 54.5 43.4 36.9 36.1
    75 55.1 43.8 37.1 36.6
    80 53.8 44.3 37.5 36.9
    85 54.6 45.2 37.9 37.1
    下载: 导出CSV
  • [1] 陈仕奇, 黄伯云.金属粉末气体雾化制备技术的研究现状与进展.粉末冶金技术, 2004, 22(5): 297 DOI: 10.3321/j.issn:1001-3784.2004.05.009

    Chen S Q, Huang B Y. The status and development of gas atomization for production of metal powders. Powder Metall Technol, 2004, 22(5): 297 DOI: 10.3321/j.issn:1001-3784.2004.05.009

    [2] 陈刚, 陈振华, 严红革, 等.一种新型的气体雾化制粉方法.中国有色金属学报, 2001, 11(增刊2): 33 https://www.cnki.com.cn/Article/CJFDTOTAL-ZYXZ2001S2007.htm

    Chen G, Chen Z H, Yan H G, et al. Novel gas atomization method for powder preparation. Chin J Nonferrous Met, 2001, 11(Suppl 2): 33 https://www.cnki.com.cn/Article/CJFDTOTAL-ZYXZ2001S2007.htm

    [3] 赵学华, 王俊文, 王晓斌, 等.纳米WC粉体的制备及其催化应用.粉末冶金技术, 2012, 30(4): 307 DOI: 10.3969/j.issn.1001-3784.2012.04.011

    Zhao X H, Wang J W, Wang X B, et al. The progress in the preparation and catalytic application of WC nano-powders. Powder Metall Technol, 2012, 30(4): 307 DOI: 10.3969/j.issn.1001-3784.2012.04.011

    [4] 熊建敏, 邹光明, 李萍, 等.扇形雾化喷嘴的射流传热实验研究.机床与液压, 2015, 15: 121 https://www.cnki.com.cn/Article/CJFDTOTAL-JCYY201515032.htm

    Xiong J M, Zou G M, Li P, et al. Experimental study on heat transfer of fan-shaped spray atomizer. Mach Tool Hydraul, 2015, 15: 121 https://www.cnki.com.cn/Article/CJFDTOTAL-JCYY201515032.htm

    [5] 王瑞祥, 张焕然, 刘建华.雾化干燥法制备球形超细氧化镍粉工艺研究.粉末冶金技术, 2013, 31(6): 413 DOI: 10.3969/j.issn.1001-3784.2013.06.003

    Wang R X, Zhang H R, Liu J H. Study on the technology of preparing ultrafine spherical Ni O by spray drying. Powder Metall Technol, 2013, 31(6): 413 DOI: 10.3969/j.issn.1001-3784.2013.06.003

    [6] 郑捷庆, 罗惕乾, 张军, 等.气力式雾化喷嘴最佳气耗率的试验.农业机械学报, 2007, 38(10): 196 https://www.cnki.com.cn/Article/CJFDTOTAL-NYJX200710051.htm

    Zheng J Q, Luo X Q, Zhang J, et al. Test on the optimum gas consumption rate of pneumatic atomizing nozzle. Trans Chin Soc Agric Mach, 2007, 38(10): 196 https://www.cnki.com.cn/Article/CJFDTOTAL-NYJX200710051.htm

    [7] 刘福平.环孔雾化喷嘴设计参数的研究.粉末冶金工业, 2011, 21(1): 40 DOI: 10.3969/j.issn.1006-6543.2011.01.007

    Liu F P. Analysis on design parameters of the atomizing nozzle with ring holes. Powder Metall Ind, 2011, 21(1): 40 DOI: 10.3969/j.issn.1006-6543.2011.01.007

    [8]

    Zeoli N, Gu S. Numerical modelling of droplet break-up for gas atomization. Comput Mater, 2006, 38: 282. DOI: 10.1016/j.commatsci.2006.02.012

    [9] 曹俊华, 饶锡新, 林立俊, 等.制粉雾化器喷嘴的流场仿真研究.机械设计与制造, 2013(7): 105 DOI: 10.3969/j.issn.1001-3997.2013.07.033

    Cao J H, Rao X X, Lin L J, et al. Research on fluid simulation of atomizer nozzle based on spray powder preparation. Mach Des Manuf, 2013(7): 105 DOI: 10.3969/j.issn.1001-3997.2013.07.033

    [10] 张维涛, 牛显明, 胡春莲, 等.双级耦合气雾化制备Ni60A粉末的技术研究.粉末冶金工业, 2012, 22(5): 1 DOI: 10.3969/j.issn.1006-6543.2012.05.001

    Zhang W T, Niu X M, Hu C L, et al. Study on the preparation Ni60A powder by gas atomization with two-stage coupling. Powder Metall Ind, 2012, 22(5): 1 DOI: 10.3969/j.issn.1006-6543.2012.05.001

    [11] 赵承庆, 姜毅.气体射流动力学.北京: 北京理工大学出版社, 1998

    Zhao C Q, Jiang Y. Gas Jet Dynamics. Beijing: Beijing Institute of Technology Press, 1998

    [12]

    Miller B, Sallam K A, Bingabr M, et al. Secondary breakup of aerated liquid jets in subsonic crossflow(Postprint). J Propul Power, 2008, 24(2): 253. DOI: 10.2514/1.30390

    [13]

    Choi J S, Mclaughlin D K, Thompson D E. Experimental on the unsteady flow field and noise generation in a centrifugal pump impeller. J Sound Vib, 2003, 263(3): 493. DOI: 10.1016/S0022-460X(02)01061-1

    [14] 李会平, 何伟.气体雾化与喷射成形中雾化过程的数学模型.材料导报, 2008, 22(3): 107 DOI: 10.3321/j.issn:1005-023X.2008.03.026

    Li H P, He W. A Mathematical model of metallic melts break-up process during gas atomization and spray forming. Material Review, 2008, 23(3): 107 DOI: 10.3321/j.issn:1005-023X.2008.03.026

    [15] 贾月梅, 赵秋霞, 赵广慧.流体力学.北京: 国防工业出版社, 2006

    Jia Y M, Zhao Q X, Zhao G H. Fluid Mechanics. Beijing: National Defend Industry Press, 2006

  • 期刊类型引用(3)

    1. 杨亚飞,姚草根,吕宏军,谢飞,孟烁. 粉末特性对W-Ti[w(Ti)=10%]合金致密度和组织的影响. 宇航材料工艺. 2024(05): 31-39 . 百度学术
    2. 杨亚飞,吕宏军,姚草根,李圣刚,王思伦. 制备工艺对W-10Ti合金组织的影响. 宇航材料工艺. 2023(05): 38-43 . 百度学术
    3. 莫蓉,李晖,孙洪涛,白英丽,陈红. 惰性熔融-红外吸收光谱法测定钨钛合金中氧含量. 广州化工. 2020(01): 95-96+146 . 百度学术

    其他类型引用(2)

图(5)  /  表(2)
计量
  • 文章访问数:  206
  • HTML全文浏览量:  143
  • PDF下载量:  38
  • 被引次数: 5
出版历程
  • 收稿日期:  2017-10-11
  • 刊出日期:  2018-04-26

目录

/

返回文章
返回