高级检索

高温液态熔渣飞行及碰撞过程研究

王凯, 仪垂杰, 胡凤超, 战胜

王凯, 仪垂杰, 胡凤超, 战胜. 高温液态熔渣飞行及碰撞过程研究[J]. 粉末冶金技术, 2022, 40(6): 535-540. DOI: 10.19591/j.cnki.cn11-1974/tf.2020050018
引用本文: 王凯, 仪垂杰, 胡凤超, 战胜. 高温液态熔渣飞行及碰撞过程研究[J]. 粉末冶金技术, 2022, 40(6): 535-540. DOI: 10.19591/j.cnki.cn11-1974/tf.2020050018
WANG Kai, YI Chui-jie, HU Feng-chao, ZHAN Sheng. Study on flight and collision process of molten blast furnace slag[J]. Powder Metallurgy Technology, 2022, 40(6): 535-540. DOI: 10.19591/j.cnki.cn11-1974/tf.2020050018
Citation: WANG Kai, YI Chui-jie, HU Feng-chao, ZHAN Sheng. Study on flight and collision process of molten blast furnace slag[J]. Powder Metallurgy Technology, 2022, 40(6): 535-540. DOI: 10.19591/j.cnki.cn11-1974/tf.2020050018

高温液态熔渣飞行及碰撞过程研究

基金项目: 国家重点研发计划资助项目(2017YFB0603602-03)
详细信息
    通讯作者:

    仪垂杰: E-mail: chuijieyi@vip.163.com

  • 中图分类号: TF09

Study on flight and collision process of molten blast furnace slag

More Information
  • 摘要:

    对离心粒化后高温熔渣的飞行过程建立数学模型,通过Runge-Kutta方法对建立的数学模型进行离散求解。结果表明,熔渣液滴沿x方向飞行距离与液滴直径和初始速度成正比;由于空气绕流阻力和重力作用,熔渣液滴速度先降低后增加。对熔渣液滴撞壁后过剩反弹能进行分析,获得了熔渣液滴的临界撞击速度。结果表明,临界撞击速度为区间,存在上界和下界,且上界和下界同时随直径增加而降低。对初始速度为10、12和14 m·s‒1三种粒化工况进行实验,结果表明,由于熔渣液滴从粒化盘抛出时速度小于粒化盘边缘线速度,熔渣液滴实际下降距离大于其理论值;三种工况下熔渣液滴撞击速度在临界撞击速度区间内,均未产生粘结。

    Abstract:

    The mathematical model was established for the flight process of the molten blast furnace slag after the centrifugal graining, and the model was discretely solved by the Runge-Kutta method. The results show that, the flight distance of the slag droplets along the x direction is proportional to the diameter and the initial velocity of droplets. Due to the air flow resistance and gravity, the velocity of droplets decreases with time, and then increases slightly. Critical impact velocity of the droplets was proposed by analyzing the excess rebound energy of the slag droplets after collision with wall. The results indicate that, the critical impact velocity is an interval, and the larger the droplets, the lower the bounds. The experiments were performed with the initial droplet velocities of 10, 12, and 14 m·s‒1. The results indicate that, the actual falling distance of the slag droplets is greater than the theoretical values, because the initial velocity of the slag droplets is less than the linear velocity of the granulation plate. Meanwhile, the impact velocity of the slag droplets is between the upper and lower bound, thus no adhesion occurs.

  • 稀有金属钼(Mo)是重要的高熔点金属,其熔点为2610 ℃,仅次于碳、钨、铼、钽和锇。金属Mo呈银白色,外形近似钢铁,具有高的硬度和弹性模量,低的蒸气压和蒸发速度,低的线膨胀系数,高的抗腐蚀能力等一系列优异特性,在现代国防、原子能工业、电真空、电光源等工程应用领域占有重要地位,在一些特殊高温应用领域甚至具有不可取代的作用[14]

    研究表明,在金属Mo基体中引入稀土氧化物粒子(可称作“稀土氧化物–Mo基材料”)可进一步提高材料的性能,拓展材料的应用。例如,通过引入氧化镧(La2O3)、氧化钇(Y2O3)等粒子对材料弥散强化,不仅可以大大提高金属Mo的室温强度和硬度,而且可以提高材料的再结晶温度,增强高温力学性能,显著延长作为高温发热体材料的使用寿命[58]。此外,在金属Mo基体中引入氧化钪(Sc2O3)、Y2O3等稀土氧化物粒子还可以提高材料的电子发射能力,用作优秀的阴极材料[9]

    作为改善金属Mo性能的稀土氧化物粒子,其尺寸大小及在Mo基体中的分布直接影响所制材料的性能。通常认为,粒子越细小,在Mo基体中分布越均匀,越有利于材料性能的提高[412],因此,设法获得粒度细小的稀土氧化物粒子、并使其均匀分布在基体中,是制备高性能稀土氧化物–Mo基材料的基础。由于熔点较高,目前难熔金属主要采用粉末冶金方法制备,而在粉末冶金工艺中,原料粉末是决定材料性能和制造成本的关键一环,要获得高性能的稀土氧化物–Mo基材料,需要首先制备出高纯度、细粒度、稀土氧化物粒子细小且掺杂分布均匀的Mo基粉末原料。与传统制备稀土氧化物–Mo基粉末的机械合金化法相比,溶液燃烧法具有掺杂少、合成效率高、能耗低等优点。特别是溶液燃烧法的合成原料均为水溶性物质,目标金属在水溶液中以离子形态存在,能够很容易实现各组分在原子或分子水平上的均匀分散和混合,这为最终得到Mo基材料中稀土氧化物弥散相的粒径细化和均匀分布提供了有利条件。

    为了增加溶液燃烧合成法的应用范围,同时为La2O3掺杂Mo合金的制备提供新思路,本文以七钼氨酸((NH4)6Mo7O24·4H2O)作为金属源,甘氨酸(C2H5O2N)为燃料,硝酸铵(NH4NO3)为氧化剂,采用溶液燃烧法合成不同质量分数La2O3掺杂的Mo前驱体粉末,并对前驱体粉末进行还原、烧结,研究La2O3掺杂量(质量分数)对粉体性能及对烧结后Mo合金各项性能的影响。

    以高可溶性的七钼氨酸((NH4)6Mo7O24·4H2O)为金属源,硝酸铵(NH4NO3)(≥99.0%)为氧化剂,甘氨酸(C2H5O2N)为燃料及添加剂,添加不同质量分数La(NO3)3·6H2O(以La2O3含量占最终合金粉末质量的比例为计算标准,分别为0、0.3%、0.7%、1.0%),通过溶液燃烧反应合成前驱体。在700 ℃下氢气氛围中还原,制备出La2O3掺杂Mo粉。对制备的粉末进行放电等离子体烧结(spark plasma sintering,SPS),烧结温度1600 ℃。

    采用X射线衍射仪(X-ray diffraction,XRD;PANalytical X-Pert PRO MPD)对未添加La2O3的氧化钼前驱体及Mo–La2O3前驱体的物相组成进行表征。采用场发射扫描电子显微镜(field emission scanning electron microscope,FESEM;Hitachi SU8020)和透射电子显微镜(transmission electron microscope,TEM)对产物的显微组织进行观察。采用能谱仪(energy disperse spectroscope,EDS)对试样中Mo和La的元素分布进行测定。

    图1为不同La2O3掺杂量的前驱体粉末微观形貌,可以清楚地发现,当不掺杂La2O3时,获得的前躯体粉末为片状结构,厚度为200 nm,片的尺寸约为0.5~2.0 μm。随着La2O3掺杂量的增加,其形貌开始变为细长颗粒状,且颗粒尺寸逐渐变小。当La2O3掺杂含量达到1.0%(质量分数)时,粉末晶粒尺寸以小于200 nm为主,且出现严重团聚现象。

    图  1  La2O3掺杂量对前驱体粉末显微形貌的影响:(a)0;(b)0.3%;(c)0.7%;(d)1.0%
    Figure  1.  Effect of La2O3 doping content (mass fraction) on the microstructure of the precursor powders: (a) 0; (b) 0.3%; (c) 0.7%; (d) 1.0%

    对不同La2O3掺杂量的前驱体粉末在700 ℃下进行还原,图2为还原产物扫描电子显微形貌。由图可以看出,制备出的La2O3掺杂Mo粉尺寸在纳米级别,随着La2O3添加量的增加,Mo粉的晶粒尺寸逐渐减小,其中掺杂质量分数为0、0.3%、0.7%和1.0%La2O3的Mo粉晶粒尺寸分别为220、180、150以及100 nm,这是由于添加La2O3抑制了Mo晶粒长大。另外,由于纳米粉末尤其是难熔金属的纳米粉末的表面积非常大,为了降低体系能量,还原后的粉末颗粒自发的聚集在一起,从而出现了不均匀的团聚现象。

    图  2  掺杂不同质量分数La2O3的Mo粉700 ℃还原产物显微形貌:(a)0;(b)0.3%;(c)0.7%;(d)1.0%
    Figure  2.  SEM images of the reduction products of the Mo powders doped by La2O3 in different mass fraction: (a) 0; (b) 0.3%; (c) 0.7%; (d) 1.0%

    图3为掺杂不同质量分数La2O3的Mo粉在700 ℃还原产物的X射线衍射图谱,由图可知,氧化钼前驱体均被还原成了Mo粉,这说明通过溶液燃烧法可以获得高纯度的La2O3掺杂Mo粉。此外,虽然在Mo粉中掺杂了不同含量的La2O3第二相粒子,但是在图中并未发现La的峰,可能是加入的La2O3所占比例非常小,在X射线衍射检测中未能发现。为了验证La2O3粒子的掺杂,实验对还原后的粉末进行了能谱分析,结果如图4所示,在掺杂质量分数为1.0%La2O3的Mo粉中发现了La特征峰,证明了La元素的存在。

    图  3  掺杂不同质量分数La2O3的Mo粉700 ℃还原产物X射线衍射图谱
    Figure  3.  XRD patterns of the Mo powders doped by La2O3 in different mass fraction after reduction at 700 ℃
    图  4  掺杂质量分数1.0%La2O3的Mo粉在700 ℃还原产物的扫描电子显微形貌(a)和对应的能谱分析(b)
    Figure  4.  SEM image (a) and the corresponding EDS analysis (b) of the Mo powders doped by 1.0%La2O3 after reduction at 700 ℃

    对还原后的粉末做进一步分析,通过透射电子显微镜对掺杂质量分数0.7%La2O3的Mo粉进行表征,结果见图5。从图中可以清楚地观察到,还原后的粉末粒径大约为150~200 nm,而且分散性较好。这主要是因为溶液燃烧法在反应过程中产生的前驱体晶粒细小,团聚体中存在大量的孔隙(如图1所示),因此在较低温度还原后,合金粉末的晶粒能够保持在纳米尺寸且分散性较好[13]

    图  5  Mo–0.7La2O3前驱体粉末透射电子显微镜照片:(a)低倍;(b)高倍
    Figure  5.  TEM images of the Mo–0.7La2O3 precursor powders: (a) low magnification; (b) high magnification

    图6为经1600 ℃烧结后La2O3掺杂Mo合金的断口形貌。和纯Mo相比,La2O3掺杂Mo合金材料的晶粒更为细小,并且随La2O3质量分数的提高,细化作用逐渐明显。可以看出,在La2O3质量分数为0.7%时,Mo晶粒尺寸为500 nm左右,继续增加La2O3质量分数至1.0%,其晶粒尺寸降至300 nm。随着La2O3掺杂量的增加,Mo–La2O3烧结体中空隙数量增加,La2O3质量分数为1.0%时,其断口形貌中孔隙数量最多。

    图  6  经1600 ℃烧结后不同质量分数La2O3掺杂Mo合金的断口形貌:(a)0;(b)0.3%;(c)0.7%;(d)1.0%
    Figure  6.  Fracture morphology of the Mo alloys doped by La2O3 in different mass fraction sintered at 1600 ℃: (a) 0; (b) 0.3%; (c) 0.7%; (d) 1.0%

    图7所示为不同La2O3掺杂量对Mo–La2O3合金相对密度的影响。可以明显看出,随着La2O3质量分数的提高,Mo合金的相对密度逐渐减小。这一方面是因为La2O3的实际密度低于纯Mo,随着掺杂量的提高,其相对密度必然会下降;另一方面,La2O3的加入会阻碍晶粒与烧结颈长大,同时阻碍晶界的迁移,使得材料的致密化行为变得困难,降低其相对密度[14]。这也与图6(d)中大量空隙相对应。

    图  7  1600 ℃烧结Mo–La2O3合金相对密度随La2O3质量分数变化
    Figure  7.  Relative density of the Mo–La2O3 alloys doped by La2O3 in different mass fraction sintered at 1600 ℃

    图8所示为Mo–La2O3合金材料的显微硬度随着La2O3掺杂量的变化。从图中可以看出,合金材料的显微硬度呈现先增加后减小的趋势,在La2O3质量分数为0.7%时,显微硬度达到最高,为HV0.2546。这是由于La2O3的加入会阻碍晶粒生长,细化晶粒,提高材料的力学性能[15]。同时,第二相粒子La2O3可以起到钉扎作用,阻碍位错的迁移,使得材料硬度提高。但是,当La2O3掺杂量过多时,样品密度降低,孔隙数量增加,从而引起硬度降低[1516]。因此当La2O3掺杂量超过0.7%时,硬度值又出现下降的趋势。

    图  8  1600 ℃烧结Mo–La2O3合金显微硬度随La2O3质量分数变化
    Figure  8.  Microhardness of the Mo–La2O3 alloys doped by La2O3 in different mass fraction sintered at 1600 ℃

    (1)将溶液燃烧法应用于纳米稀土氧化物掺杂Mo基材料的制备,成功制备出La2O3掺杂Mo合金粉,并经烧结获得合金样品,所制备合金样品具有优异的力学性能。

    (2)随着La2O3掺杂量(质量分数)的增加,溶液燃烧合成制备的前驱体粉末逐渐由片状大颗粒变成细小的不规则颗粒。在掺杂量为1.0%时,前驱体粉末晶粒尺寸在200 nm左右。经还原后得到的Mo–La2O3粉末晶粒尺寸随着La2O3掺杂量的增加而减小,在掺杂量为1.0%时,晶粒尺寸为100 nm左右。

    (3)所制得的La2O3掺杂Mo粉经1600 ℃烧结后产物相对密度在均在95%以上,随着La2O3掺杂量的增加(La2O3质量分数在0~1.0%范围内),相对密度逐渐降低,而显微硬度呈现先上升后下降的趋势。在La2O3掺杂量为0.7%时,Mo–La2O3合金显微硬度呈现出最大值,此时晶粒尺寸为500 nm左右,显微硬度达到HV0.2564。

  • 图  1   粒化过程示意图(a)和熔渣液滴受力分析(b)

    Figure  1.   Schematic diagram of the granulating (a) and the force analysis for the slag droplets (b)

    图  2   熔渣液滴飞行轨迹(d=2 mm,t=1 s)

    Figure  2.   Flight path of the slag droplets (d=2 mm, t=1 s)

    图  3   熔渣液滴沿x轴飞行距离(t=1 s)

    Figure  3.   Flight distance of the slag droplets along x direction (t=1 s)

    图  4   熔渣液滴速度随时间变化图(d=2 mm)

    Figure  4.   Velocity of the slag droplets with the different times (d=2 mm)

    图  5   熔渣液滴撞击过剩反弹能与撞击速度关系(T=1400 ℃,α=60°)

    Figure  5.   Dependence between the excess rebound energy and the impact velocity (T=1400 ℃, α=60°)

    图  6   渣粒直径分布

    Figure  6.   Distribution of the slag diameter

    图  7   熔渣液滴沿y轴下降距离(x=1 m)

    Figure  7.   Falling distance of the slag droplets along y direction (x=1 m)

    图  8   撞击壁面时的速度(x=1 m)

    Figure  8.   Impact velocity of the slag droplets (x=1 m)

    表  1   熔渣化学成分(质量分数)

    Table  1   Chemical composition of the blast slag %

    SiO2 CaO Al2O3 MgO FeO MnO TiO2 S 其他
    33.92 40.15 14.50 8.17 0.33 0.30 0.50 1.03 1.10
    下载: 导出CSV

    表  2   熔渣液滴直径与临界撞击速度关系

    Table  2   Relationship between the slag droplet diameter and the critical impact velocity

    d / mm v* / (m·s‒1)
    0.5 0.45~227.32
    1.5 0.33~75.75
    2.5 0.27~45.44
    3.5 0.24~32.45
    4.5 0.22~25.23
    下载: 导出CSV

    表  3   实验工况

    Table  3   Experimental condition

    序号 Dm / m D / m N / (r·min‒1) v0 / (m·s‒1)
    1 2 0.25 764 10
    2 2 0.25 917 12
    3 2 0.25 1069 14
    下载: 导出CSV
  • [1]

    Wu J J, Wang H, Zhu X, et al. Centrifugal granulation performance of liquid with various viscosities for heat recovery of blast furnace slag. Appl Therm Eng, 2015, 89: 494 DOI: 10.1016/j.applthermaleng.2015.06.031

    [2]

    Worldsteel Association. Steel statistical yearbook [J/OL]. Worldsteel Association [2020-1]. Steel-Statistical-Yearbook-2019-concise-version. pdf (worldsteel. org)

    [3]

    Purwanto H, Akiyama T. Hydrogen production from biogas using hot slag. Int J Hydrogen Energy, 2006, 31(4): 491 DOI: 10.1016/j.ijhydene.2005.04.021

    [4]

    Li P, Yu Q B, Qin Q, et al. Kinetics of CO2/coal gasification in molten blast furnace slag. Ind Eng Chem Res, 2012, 51(49): 15872 DOI: 10.1021/ie301678s

    [5]

    Pickering S J, Hay N, Roylance T F, et al. New process for dry granulation and heat recovery from molten blast-furnace slag. Ironmaking Steelmaking, 1985, 12(1): 14

    [6]

    Kashiwaya Y, In-nami Y, Akiyama T. Mechanism of the formation of slag particles by the rotary cylinder atomization. ISIJ Int, 2010, 50(9): 1252 DOI: 10.2355/isijinternational.50.1252

    [7]

    Akiyama T, Toshio M, Jun-Ichiro Y, et al. Feasibility study of hydrogen generator with molten slag granulation. Steel Res Int, 2004, 75(2): 122 DOI: 10.1002/srin.200405937

    [8]

    Xie D, Washington B M, Norgate T, et al. Dry granulation of slags–turning waste into valuable products. CAMP-ISIJ, 2005, 18(4): 1088

    [9]

    Xie D, Pan Y, Flann R, et al. Heat recovery from slag through dry granulation // 1st CSRP Annual Conference. Melbourne, 2007: 1

    [10] 于庆波, 刘军祥, 窦晨曦, 等. 转杯法高炉渣粒化实验研究. 东北大学学报(自然科学版), 2009, 30(8): 1163 DOI: 10.3321/j.issn:1005-3026.2009.08.025

    Yu Q B, Liu J X, Dou C X, et al. Dry granulation experiment of blast furnace slag by rotary cup atomizer. J Northeast Univ Nat Sci, 2009, 30(8): 1163 DOI: 10.3321/j.issn:1005-3026.2009.08.025

    [11] 刘军祥, 于庆波, 李朋, 等. 高炉渣干法粒化试验研究. 钢铁, 2010, 45(2): 95 DOI: 10.13228/j.boyuan.issn0449-749x.2010.02.032

    Liu J X, Yu Q B, Li P, et al. Experimental study on dry-granulation of molten blast furnace slag. Iron Steel, 2010, 45(2): 95 DOI: 10.13228/j.boyuan.issn0449-749x.2010.02.032

    [12]

    Lin B, Wang H, Zhu X, et al. Crystallization properties of molten blast furnace slag at different cooling rates. Appl Therm Eng, 2016, 96: 432 DOI: 10.1016/j.applthermaleng.2015.11.075

    [13]

    Zhu X, Ding B, Wang H, et al. Phase evolution of blast furnace slags with variation in the binary basicity in a variable cooling process. Fuel, 2018, 219: 132 DOI: 10.1016/j.fuel.2018.01.075

    [14]

    Zhu X, Ding B, Wang H, et al. Numerical study on solidification behaviors of a molten slag droplet in the centrifugal granulation and heat recovery system. Appl Therm Eng, 2018, 130: 1033 DOI: 10.1016/j.applthermaleng.2017.11.080

    [15]

    Luo S Y, Fu J, Zhou Y M, et al. The production of hydrogen-rich gas by catalytic pyrolysis of biomass using waste heat from blast-furnace slag. Renew Energy, 2017, 101: 1030 DOI: 10.1016/j.renene.2016.09.072

    [16] 米沙, 谢锴, 孙岱, 等. 冶金渣颗粒与空气间的换热和阻力特性. 中国有色金属学报, 2015, 25(7): 1993 DOI: 10.19476/j.ysxb.1004.0609.2015.07.032

    Mi S, Xie K, Sun D, et al. Heat transfer and resistance characteristics between metallurgical slag particles and air. Chin J Nonferrous Met, 2015, 25(7): 1993 DOI: 10.19476/j.ysxb.1004.0609.2015.07.032

    [17] 王崇琳. 粉末飞行之研究I粉末在静止气体场中的飞行轨迹. 粉末冶金技术, 2008, 26(4): 243

    Wang C L. Investigation on the flying of powder particles Ⅰ Trajectories of flying powder particles in static atmosphere. Powder Metall Technol, 2008, 26(4): 243

    [18]

    Richter A, Nikrityuk P A. Drag forces and heat transfer coefficients for spherical, cuboidal and ellipsoidal particles in cross flow at sub-critical Reynolds numbers. Int J Heat Mass Trans, 2012, 55(4): 1343 DOI: 10.1016/j.ijheatmasstransfer.2011.09.005

    [19]

    Ni J J, Yu G S, Guo Q H, et al. Submodel for predicting slag deposition formation in slagging gasification systems. Energy Fuels, 2011, 25(3): 1004 DOI: 10.1021/ef101696a

    [20]

    Liu J X, Yu Q B, Duan W J, et al. Experimental investigation on ligament formation for molten slag granulation. Appl Therm Eng, 2014, 73(1): 888 DOI: 10.1016/j.applthermaleng.2014.08.042

    [21]

    Sun Y Q, Shen H W, Wang H, et al. Experimental investigation and modeling of cooling processes of high temperature slags. Energy, 2014, 76: 761 DOI: 10.1016/j.energy.2014.08.073

    [22]

    Ding B, Wang H, Zhu X, et al. Crystallization behaviors of blast furnace (BF) slag in a phase-change cooling process. Energy Fuels, 2016, 30(4): 3331 DOI: 10.1021/acs.energyfuels.5b03000

    [23]

    Han C, Chen M, Zhang W D, et al. Viscosity model for iron blast furnace slags in SiO2–Al2O3–CaO–MgO system. Steel Res Int, 2015, 86(6): 678 DOI: 10.1002/srin.201400340

    [24] 王利明, 袁意林, 邵毅敏, 等. 二十辊轧机轧辊磨床砂轮动不平衡对磨削颤振的影响. 工程科学学报, 2015, 37(suppl1): 78

    Wang L M, Yuan Y L, Shao Y M, et al. Chatter analysis about roll grinder of twenty-high rolling mill in grinding process with grinding wheel dynamic imbalance fault. Chin J Eng, 2015, 37(Suppl1): 78

    [25]

    Dhirhi R, Prasad K, Shukla A, et al. Experimental study of rotating dry slag granulation unit: Operating regimes, particle size analysis and scale up. Appl Therm Eng, 2016, 107: 898 DOI: 10.1016/j.applthermaleng.2016.07.049

    [26]

    Wang L Y, Sun W Q, Li X L, et al. Flight dynamics and sensible heat recovery of granulated blast furnace slag. Open Fuels Energy Sci J, 2015, 8: 356 DOI: 10.2174/1876973X01508010356

  • 期刊类型引用(3)

    1. 沙逢源,荣川,杨劼人,刘颖,邢辉,李志强. Ti-48Al-3Nb-2V-(Nd, Pr)合金显微组织演化规律与结构调控. 内蒙古工业大学学报(自然科学版). 2024(05): 443-452 . 百度学术
    2. 马静,王铁凝,姜秋月,冯志浩,张欣,李建辉. Al_2O_3-Ce_2O_3复合薄膜对304不锈钢抗高温氧化性能的影响. 材料保护. 2024(10): 19-26 . 百度学术
    3. 徐强林,吴小成,张玲. 建筑铝合金塑性变形件中裂纹源位置对力学性能和抗氧化性能的影响. 铸造. 2024(12): 1714-1720 . 百度学术

    其他类型引用(0)

图(8)  /  表(3)
计量
  • 文章访问数:  196
  • HTML全文浏览量:  46
  • PDF下载量:  48
  • 被引次数: 3
出版历程
  • 收稿日期:  2020-12-13
  • 录用日期:  2020-12-13
  • 刊出日期:  2022-12-27

目录

/

返回文章
返回