罗茨真空泵新型偏心渐开线型不对称转子的构建与性能研究

上一篇

下一篇

张轩宇, 刘俊材, 王君, 杨舒然. 罗茨真空泵新型偏心渐开线型不对称转子的构建与性能研究[J]. 真空科学与技术学报, 2024, 44(1): 78-86. doi: 10.13922/j.cnki.cjvst.202309001
引用本文: 张轩宇, 刘俊材, 王君, 杨舒然. 罗茨真空泵新型偏心渐开线型不对称转子的构建与性能研究[J]. 真空科学与技术学报, 2024, 44(1): 78-86. doi: 10.13922/j.cnki.cjvst.202309001
Xuanyu ZHANG, Juncai LIU, Jun WANG, Shuran YANG. Construction and Performance Study of a Novel Eccentric Involute Asymmetric Rotor of Roots Vacuum Pumps[J]. zkkxyjsxb, 2024, 44(1): 78-86. doi: 10.13922/j.cnki.cjvst.202309001
Citation: Xuanyu ZHANG, Juncai LIU, Jun WANG, Shuran YANG. Construction and Performance Study of a Novel Eccentric Involute Asymmetric Rotor of Roots Vacuum Pumps[J]. zkkxyjsxb, 2024, 44(1): 78-86. doi: 10.13922/j.cnki.cjvst.202309001

罗茨真空泵新型偏心渐开线型不对称转子的构建与性能研究

    通讯作者: Tel: 13589269049; E-mail: wangjun@upc.edu.cn
  • 中图分类号: TH36

Construction and Performance Study of a Novel Eccentric Involute Asymmetric Rotor of Roots Vacuum Pumps

    Corresponding author: Jun WANG, wangjun@upc.edu.cn
  • MSC: TH36

  • 摘要: 罗茨转子作为罗茨真空泵的重要组成部件,对罗茨真空泵性能影响极大。为了解决传统渐开线型转子在渐开线与圆弧连接处存在不光滑连接点,导致泄漏加剧,运行效率下降的问题,文章采用偏心渐开线与齿顶、齿根圆弧平滑相连,提出了一种偏心渐开线型不对称转子。建立了偏心渐开线型不对称转子的几何模型,分析了转子的几何参数对面积利用率等性能的影响。结果表明:随着圆心角增大及节圆半径减小,所提出的偏心渐开线型不对称转子的面积利用率提高;与传统渐开线型转子相比,所提出的偏心渐开线型不对称转子每转排量提升了17.2%,出口流量脉动降低9.8%,出口压力脉动降低7.9%。研究内容对罗茨真空泵的发展具有重要的意义。
  • 加载中
  • 图 1  转子坐标系

    Figure 1.  Coordinate systems for rotor

    图 2  带尖点的渐开线凸轮转子

    Figure 2.  Involute type Roots rotor with sharp points

    图 3  左转子四分之一型线

    Figure 3.  One quarter of the left rotor profile

    图 4  偏心圆光滑连接两圆弧

    Figure 4.  Left rotor profile connected by arc

    图 5  渐开线沿顺时针方向展开

    Figure 5.  The involute unfolds clockwise

    图 6  右转子四分之一型线:BP<CP

    Figure 6.  One quarter of the right rotor profile: BP<CP

    图 7  右转子四分之一型线:BP=CP

    Figure 7.  One quarter of the right rotor profile: BP=CP

    图 8  右转子四分之一型线:BP>CP

    Figure 8.  One quarter of the right rotor profile: BP>CP

    图 9  渐开线的发生角。(a) BP<CP,(b) BP>CP

    Figure 9.  The involute initial angle. (a) BP<CP, (b) BP>CP

    图 10  不同半圆心角θ下的转子

    Figure 10.  Design results of new involute type rotor profile

    图 11  转子运动轨迹图

    Figure 11.  Locus of rotor profile

    图 12  圆心角θ与面积利用系数的关系

    Figure 12.  Relationship between central angle θ and area utilization ratio

    图 13  三维网格。(a)不对称渐开线型转子,(b)带尖点的渐开线型转子

    Figure 13.  Three dimensional grid. (a) Asymmetric involute rotor, (b) involute rotor with sharp points

    图 14  两种转子的每转平均质量流量

    Figure 14.  The average mass flow at outlet of the two kinds of rotors

    图 15  两种转子出口处的流量脉动

    Figure 15.  Flow pulsation at the outlet of two kinds of rotors

    图 16  两种转子出口处的压力脉动

    Figure 16.  Pressure pulsation at the outlet of two kinds of rotors

    表 1  验证网格无关性

    Table 1.  The grid independence verification

    组别网格数量平均质量流量(kg/s)
    A1538270.0156
    B2052650.0170
    C2703850.0177
    D3155730.0179
    下载: 导出CSV
  • [1] 徐文兵, 于振华, 胡焕林. 偏心圆弧罗茨转子型线的研究[J]. 真空,2006,43(1):6−8(in Chinese) Xu W B, Yu Z H, Hu H L. On the profile of eccentric-arc rotor of Roots pumps[J]. Vacuum,2006,43(1):6−8
    [2] Liu L L, Chu J P, Hu J Z. Study on profile of roots vacuum pump rotor[J]. Journal of Machine Design, 2007(03):64-67(刘林林, 初嘉鹏, 胡建中. 罗茨真空泵转子型线的研究[J]. 机械设计, 2007(03):64-67(in Chinese)
    [3] Hsieh C F, Hwang Y W. Tooth profile of a Roots rotor with a variable trochoid ratio[J]. Mathematical and Computer Modelling,2008,48(1-2):19−33 doi: 10.1016/j.mcm.2007.08.008
    [4] Hwang Y W, Hsieh C F. Study on high volumetric efficiency of the roots rotor profile with variable trochoid ratio[J]. Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, 2006, 220(9):1375-1384
    [5] Hsieh C F, Hwang Y W. Study on the high-sealing of Roots rotor with variable trochoid ratio[J]. Journal of Mechanical Design,2007,129(12):1278−1284 doi: 10.1115/1.2779897
    [6] 秦丽秋, 刘玉岱. 罗茨泵圆弧转子型线研究[J]. 真空, 1990(01):32-39(in Chinese) Qin L Q, Liu Y D, Study on the profile of circular arc rotor of roots pump[J]. Vacuum, 1990(01):32-39
    [7] Wang J, Liu R, Yang S, et al. Geometric study and simulation of an elliptical rotor profile for roots vacuum pumps[J]. Vacuum,2018,153:168−175 doi: 10.1016/j.vacuum.2018.04.014
    [8] 刘厚根, 李皓, 周鹏. 基于齿条法的罗茨机械增压器转子型线设计[J]. 机械设计与研究, 2018, 34(03):86−88+93(in Chinese) Liu H G, Li H, Zhou P. Rotor profile design of roots supercharger based upon rack generation method[J]. Machine Design and Research, 2018, 34(03):86−88+93
    [9] 朱晓东, 赵厚继. 罗茨鼓风机渐开线型转子的改进分析[J]. 风机技术, 2009(5):19-21(in Chinese) Liu H G, Zhu X D, Zhao H J. The improvement analysis of involute profile type rotor in roots blower[J]. Chinese Journal of Turbomachinery, 2009(5):19-21
    [10] 李海洋, 赵玉刚, 胡柳, 等. 渐开线型罗茨真空泵转子型线的改进研究[J]. 机床与液压, 2011, 39(22):37-39(in Chinese) Li H Y, Zhao Y G, Hu L, et al. The improvement study on involute profile type rotor profile in roots vacuum pump[J]. Machine Tool & Hydraulics, 2011, 39(22):37−39
    [11] 沈浩, 贾琨, 黎义斌, 等. 凸轮泵转子型线优化设计及模拟计算[J]. 兰州理工大学学报,2014,40(3):44−48(in Chinese) Shen H, Jia K, Li Y B, et al. Optimization design of cam pump rotor profile and simulation of pump performance[J]. Journal of Lanzhou University of Technology,2014,40(3):44−48
    [12] 朱超颖, 林景殿, 苏中地. 一种新罗茨转子型线的构成方法[J]. 真空科学与技术学报, 2015(12):1449−1452(in Chinese) Zhu C Y, LIN J D, Su Z D. A novel type of roots rotor profile[J]. Chinese Journal of Vacuum Science and Technology, 2015(12):1449−1452
    [13] 刘坤, 巴德纯, 张振厚, 等. 圆弧包络线在罗茨转子型线设计中的应用[J]. 真空, 2007(01):14−17(in Chinese) Liu K, Ba C D, Zhang Z H, et al. Application of enveloping curve of circular arcs in design of roots rotor profile[J]. Vacuum, 2007(01):14−17
    [14] Hsieh C F. A new curve for application to the rotor profile of rotary lobe pumps[J]. Mechanism and Machine Theory,2015,87:70−81 doi: 10.1016/j.mechmachtheory.2014.12.018
    [15] 李玉龙, 冉光泽, 张宸赫. 渐开线转子型线的参数化与最优化研究[J]. 成都大学学报(自然科学版), 2017, 36(04):398−401(in Chinese) Li Y L, Ran G Z, Zhang C H. Parameterization and optimization researchon involute profile of rotor[J]. Journal of Chengdu University(Natural Science Edition), 2017, 36(04):398−401
    [16] 肖芝, 张小萍, 翟旭军, 等. 扭叶罗茨鼓风机转子型线改进设计与内流数值模拟分析[J]. 制造业自动化, 2012, 34(21):47−51(in Chinese) Xiao Z, Zhang X P, Zhai X J, et al. The modification design of rotor profiles and numerical simulation analysis of interior field for twisted roots blower[J]. Manufacturing Automation, 2012, 34(21):47-51
    [17] 戴映红. 气冷式罗茨真空泵的转子型线设计及流场分析[D]. 浙江工业大学, 2010(in Chinese) Dai Y H. Profile design and numerical analysis on air-cooled roots vacuum pump[D]. Zhejiang University of Technology, 2010
    [18] 张宇. 罗茨真空泵内部流场的研究[D]. 东北大学, 2010(in Chinese) Zhang Y. Flow Field investigation on roots pump[D]. Northeastern University, 2010
    [19] 岳向吉, 巴德纯, 刘坤, 等. 干式罗茨真空泵吸气级内流动的瞬态模拟[J]. 真空科学与技术学报, 2012, 32(09):850−855(in Chinese) Yue X J, Ba D C, Liu K, et al. Transient numerical simulation of gas flow in suction stage of dry roots-type vacuum pump[J]. Chinese Journal of Vacuum Science and Technology, 2012, 32(09):850−855
    [20] 陈长琦, 王旭迪. 罗茨泵径向间隙泄漏量的研究与建模[J]. 流体机械, 2002(11):29−31(in Chinese) Chen C Q, Wang X D. Research and modeling of the leakage through the radial clearance of the roots vacuum pump[J]. Fluid Machinery, 2002(11):29−31
  • 加载中
图( 16) 表( 1)
计量
  • 文章访问数:  162
  • HTML全文浏览数:  162
  • PDF下载数:  1
  • 施引文献:  0
出版历程
  • 收稿日期:  2023-09-04
  • 刊出日期:  2024-01-31

罗茨真空泵新型偏心渐开线型不对称转子的构建与性能研究

    通讯作者: Tel: 13589269049; E-mail: wangjun@upc.edu.cn
  • 中国石油大学(华东)新能源学院 青岛 266580

摘要: 罗茨转子作为罗茨真空泵的重要组成部件,对罗茨真空泵性能影响极大。为了解决传统渐开线型转子在渐开线与圆弧连接处存在不光滑连接点,导致泄漏加剧,运行效率下降的问题,文章采用偏心渐开线与齿顶、齿根圆弧平滑相连,提出了一种偏心渐开线型不对称转子。建立了偏心渐开线型不对称转子的几何模型,分析了转子的几何参数对面积利用率等性能的影响。结果表明:随着圆心角增大及节圆半径减小,所提出的偏心渐开线型不对称转子的面积利用率提高;与传统渐开线型转子相比,所提出的偏心渐开线型不对称转子每转排量提升了17.2%,出口流量脉动降低9.8%,出口压力脉动降低7.9%。研究内容对罗茨真空泵的发展具有重要的意义。

English Abstract

  • 罗茨真空泵由一对相互啮合的罗茨转子、泵腔、吸气口和排气口组成。罗茨真空泵工作时,两个转子由齿轮驱动,在相反方向上进行同步旋转运动。

    罗茨真空泵结构简单、可靠性高、抽速大、效率高,在冶金、化工,以及航空航天、食品工业等领域的各环节中应用广泛。

    罗茨转子是罗茨真空泵的关键部件,转子型线设计是提高泵性能的关键因素。徐文兵等[1]提出了一种包含偏心圆弧的转子型线,有面积利用率高、啮合严密的优点。刘林林等[2]提出了一种转子型线,由渐开线、摆线、圆弧组成,提高了转子运行的稳定性。Hsieh C F和Hwang Y W[3-5]将变摆线比的摆线应用于罗茨转子的设计中,并介绍了其无根切型线构建方法,相较传统摆线具有设计多样性,研究表明所提出的变摆线型转子有较高的面积利用系数和密封性。秦丽秋等[6]使用圆弧及其共轭曲线设计了两种罗茨泵圆弧转子型线。王君等[7]提出一种相较圆弧转子有更大面积利用系数的椭圆弧型转子,研究了消除该型线的根切现象和余隙容积的方法。刘厚根等[8]将齿条法应用于罗茨机械的型线设计之中,采用Auto CAD模拟加工以得到转子型线。刘厚根等[9]通过改进传统渐开线转子齿顶圆弧,使渐开线长度增加,理论效率有所提高,改善了啮合重合度。

    在工作过程中,转子型线的不光滑连接点易于磨损,导致转子的力学性能降低,严重限制了罗茨转子的工作性能和可靠性。李海洋等[10]优化了传统渐开线型罗茨泵的转子型线,避免了叶片间的干涉,提高了转子的径距比和面积利用率。沈浩等[11]采用高次曲线,消除了转子型线的不光滑点,改进后的罗茨真空泵出口处的脉动现象得到了改善。朱超颖等[12]基于传统渐开线型转子,提出了一种由渐开线和多段圆弧组成的新型罗茨转子,新型线光滑性较好,面积利用系数较高。

    为研究罗茨转子的工作特性,得到不同型线样式对罗茨转子工作性能的影响规律,刘坤等[13]对圆弧型转子的面积利用率与形状参数等的关系进行了研究;Hsieh C F[14]提出了一种椭圆轮盘转子,研究了不同椭圆轴比对效率、瞬时质量流量、流速、压力脉动的影响。李玉龙等[15]通过改变节圆半径等参数对渐开线转子型线进行了优化。肖芝[16]等利用Fluent软件对改进后的扭叶转子罗茨风机的内部气流流动进行了研究。戴映红[17]设计了用于气冷式罗茨真空泵的罗茨转子,并对其进行了数值模拟。张宇[18]采用了结构化动网格,对空冷罗茨泵进行了数值模拟,研究了不同位置的速度及压力分布,分析了其脉动特性。岳向吉等[19]通过数值模拟研究了罗茨真空内部流场,分析了罗茨真空泵内旋涡的尺寸、方位与强度。陈长琦等[20]建立了罗茨泵间隙泄漏模型,对不同间隙下的部分工作参数进行了理论计算。

    渐开线-圆弧型罗茨转子在渐开线与圆弧连接处存在不光滑连接点,长期运行转子易发生磨损,转子的力学性能变差。此外,不光滑连接点的磨损,导致转子间气体泄漏严重,真空泵的性能降低。针对以上问题,本文采用偏心渐开线与齿顶、齿根两段圆弧平滑相连。采用数值模拟对转子几何参数与面积利用率之间的关系进行了研究。对比了罗茨真空泵新型不对称转子与传统带尖点的罗茨转子的内部流动特性。研究内容丰富了罗茨转子型线,对提高罗茨真空泵的性能具有重要意义和价值。

  • 转子动静坐标系如图1所示。左转子型线方程由x1O1y1变换至x2O2y2的方程为:

    式中:r1为左转子型线在x1O1y1中确定的方程,r2为左转子型线在x2O2y2中确定的方程。

    传统带尖点的渐开线型罗茨转子截面型线如图2所示,为消除尖点构建新型渐开线不对称转子模型。

  • 左转子曲线组成如图3所示,转子呈中心对称结构。转子由销齿圆弧ABCD,偏心渐开线BC组成。圆弧与偏心渐开线连接处可平滑过渡。

    各段曲线的方程如下所示。

    圆弧AB的方程为:

    圆弧CD的方程为:

    对于偏心渐开线BC的方程可分为以下三种情况:

    延长圆弧AB的半径BM和圆弧CD的半径NC交与点P,得到线段BPCP

    BP<CP时,渐开线BC基圆圆心Ob在点P之下,此时偏心渐开线沿逆时针方向展开,如图3所示,其方程为:

    BP=CP时,渐开线BC变作以点P为圆心的圆弧,如图4所示,其方程为:

    BP>CP时,基圆圆心Ob在点P之上,此时偏心渐开线沿顺时针方向展开,如图5所示,其方程为:

  • 右转子由圆弧和偏心渐开线的共轭曲线组成,右转子齿顶圆弧与左转子齿底圆弧相同,右转子齿底圆弧与左转子齿顶圆弧相同,故求出偏心渐开线的共轭曲线方程即可得到右转子型线。

    由啮合原理可知,互为共轭曲线的两曲线一定满足如下啮合方程:

    式中,N1为过两曲线啮合点的公法矢量:

    $ {\boldsymbol{V}}_1^{(12)} $为啮合点处的相对速度矢量:

    式中:

    为了便于计算,取dφ/dt值为1。由此,可化简式(9)得到:

    将式(8)和式(13)代入式(7)可得:

    将式(4)-(6)分别代入式(14),并进一步化简得到曲线BC及其共轭曲线的啮合条件。

    BP<CP时:

    BP=CP时:

    BP>CP时:

    将式(4)-(6)分别代入式(1)可以得到三种情形下的共轭曲线方程:

    BP<CP时:

    图6所示,为BP<CP时的1/4右转子型线。

    BP=CP时:

    图7所示,为BP=CP时1/4右转子型线。

    BP>CP时:

    图8所示,为BP>CP时1/4右转子型线。

  • 根据渐开线的基本性质,无论BPCP的长度关系如何,仅当BPCPBC的基圆相切时,偏心渐开线可以与圆弧ABCD平滑连接,且有以下几何关系:

    根据图4中的几何关系可知,$ BP \bot CP $。因此,基圆半径为:

    因存在图5BP>CP的情况,基圆半径为:

    根据图3中的几何关系可知,点B、点C的坐标分别为:

    直线BP的方程可表示为:

    直线CP的方程可表示为:

    将式(27)、(28)联立解得交点P的坐标为:

    线段BPCP的长度分别为:

  • 由于BC的基圆和BPCP都是相切的,所以可将BPCP平移来求出基圆圆心的坐标。

    BP<CP时,基圆的圆心在P点之下。把直线BP往下平移Rb/sinθ可得:

    把直线CP往下平移Rb/cosθ可得:

    联立式(31)、(32)可得基圆圆心点Ob的坐标为:

    BP=CP时,点P即为圆弧BC的圆心,其坐标如式(29)所示。

    BP>CP时,基圆的圆心在P点之方。把直线BP往上平移Rb/sinθ可得:

    把直线CP往上平移Rb/cosθ可得:

    联立式(34)、(35)可得基圆圆心点Ob的坐标为:

  • 把发生角为0的渐开线围绕基圆的圆心Ob转动可求出偏心渐开线的发生角。

    BP<CP时,如图9(a)所示,在RtCJOb中,根据勾股定理可知:

    根据渐开线的基本性质,代入各点的坐标后可得:

    联立式(38)、(39)可得:

    则角度参数t可取区间为:t∈[t1t1+π/2]。

    因此,C点在渐开线旋转之前的相应点C´坐标:

    则可得偏心渐开线的发生角α

    同理,当BP>CP时,如图9(b)所示,可得:

    则角度参数t可取区间为:t∈[t1t1+π/2]。

    B点在渐开线旋转之前的相应点B´坐标:

    则可得偏心渐开线的发生角α

  • 面积利用系数是指罗茨真空泵每转的排气面积与容腔面积的比值。面积利用系数越高,真空泵的理论吸气量越大。面积利用系数η的表达式如下:

    其中,Al表示左转子的截面积,mm2Ar表示右转子的截面积,mm2

    为将圆弧圆心角度θ取不同大小时的转子面积利用率进行比较,设定R1为60 mm和R2为42.5 mm作为参考。圆心角θ取不同数值时的转子截面构建结果如图10所示。转子的运行轨迹如图11所示,两转子可以实现正确啮合。面积利用系数如图12所示,从图中可以看出:当转子外圆半径R1为一固定值时,加大圆心角2θ、减小节圆半径R2,可以使转子面积利用率η得到提高。

  • 为便于对比分析带尖点的渐开线型和新型不对称渐开线型罗茨真空泵的工作性能,两种转子采用相同的模型参数并设置相同的边界条件。转子的长度l为170 mm,转子的外圆半径R1为60 mm,节圆半径R2为42.5 mm,圆心角θ为25º,转子顶部轮廓与气缸内壁之间的距离β1为0.15 mm,两转子边缘间的距离β2为0.3 mm。工作介质选择理想空气,进口压力为2×104 Pa,出口压力为3×104 Pa,设定进口温度为295 K,采用标准k-ε模型,转速为3000 r/min,设定计算的一个周期为0.06 s。

    采用四种网格数量对偏心渐开线型罗茨泵开展数值模拟,其中流体域为结构化网格,进出口为非结构化网格,如表1图13所示:

    表1中计算结果可知,C组求得质量流量与D组相比低1.12%,B组求得质量流量与C组相比低3.95%。可知网格数量高于270385时,网格数量对数值模拟影响较小,为提高计算效率,采用网格数量270385的计算结果进行后续数值分析。

  • 转子区域的面积利用率可以通过出口的平均质量流量来体现。从如图14所示的两类罗茨转子真空泵转子的平均质量流量与时间关系曲线中可以看出,在旋转的第二周(0.02 s),两种转子的出口流量脉动均逐渐变得稳定,至第三周时(0.04 s)便已经稳定,带尖点的渐开线型罗茨转子真空泵和不对称渐开线型罗茨转子真空泵每转平均质量流量分别为:0.0151 kg/s、0.0177 kg/s,不对称渐开线型罗茨转子真空泵每转排量比带尖点的渐开线型罗茨转子真空泵高出约17.2%。

    图15所示,为两种转子在第三转时出口处流量脉动图。从图中可以看出,转子每转动一圈,脉动周期的个数为4。

    现定义流量脉动系数F,其计算公式为:

    式中:Fmax为质量流量的瞬时最大值,kg/s;Fmin为质量流量的瞬时最小值,kg/s;Fave为质量流量的平均值(每转),kg/s。

    带尖点的渐开线型罗茨转子真空泵和不对称渐开线型罗茨转子真空泵出口处的流量脉动系数分别为:4.2089、3.7980,与带尖点的渐开线型罗茨转子真空泵相比,不对称渐开线型罗茨转子真空泵在出口处的流量脉动降低约9.8%。

    图16所示,为两种转子在第三转时出口处压力脉动图。从图中可以看出,转子每转动一圈,脉动周期的个数为4。

    定义压力脉动系数p,其计算公式为:

    式中:pmax为出口压力最大值,Pa;pmin为出口压力最小值,Pa;pave为出口压力均值,Pa。

    带尖点的渐开线型罗茨转子真空泵和不对称渐开线型罗茨转子真空泵出口流量脉动系数分别为:0.0777、0.0716,与带尖点的渐开线型罗茨转子真空泵相比,不对称渐开线型罗茨转子真空泵在出口处的压力脉动降低约7.9%。

  • (1)为优化罗茨转子型线,提高转子面积利用率,进而提高转子的工作性能,将偏心渐开线与齿顶、齿根圆弧相连实现平滑过渡,减少了转子的不光滑连接点,改善了转子的力学性能。此外转子的型线种类的减少,简化了转子的加工过程,降低了加工成本。运用啮合原理求解了渐开线平移、旋转变换后的共轭曲线,构建了不对称渐开线型转子的数学模型,推导了转子截面型线方程。

    (2)得到BP<CPBP>CP时渐开线角度参数的取值范围及偏心渐开线的发生角。研究了圆弧的圆心角θ与转子面积利用系数之间的关系,发现当外圆半径R1为一固定值时,通过加大圆心角2θ及减小节圆半径R2,可以使转子面积利用率得到提高。

    (3)通过数值模拟对比了所提出的不对称渐开线型罗茨转子真空泵和传统渐开线型罗茨空泵的内部流场,发现相较于传统罗茨真空泵:每转排量提升17.2%,出口流量脉动降低9.8%,出口压力脉动降低7.9%。

参考文献 (20)

目录

/

返回文章
返回