-
When a hypersonic vehicle flies, it will have friction with the atmosphere, ionizing the surrounding air, and producing a plasma sheath containing a large number of free electrons. The plasma sheath will cause the electromagnetic wave to seriously attenuate, which will result in communication interruption and other problems. With the gradual realization of terahertz wave technology, its high penetrability and anti-interference performance provides an important way to solve the blackout problem. Thus, the using of the terahertz wave to solve the blackout problem encountered during vehicle reentry is of great significance to studying the transmission performance of terahertz wave in the plasma sheath. This article refers to the public data of the plasma sheath during the reentry of the RAM vehicle. Considering the asymmetry of the sheath density distribution, a double Gaussian distribution is used to simulate the longitudinal electron density distribution. Based on the SMM algorithm, the article uses the magnetization direction, electron density, external magnetic field strength, collision frequency of the non-uniformly magnetized plasma as variables, and their effects on left-hand and right-hand polarized terahertz wave under normal incidence are studied. The results show that these parameters have obvious effects on the transmission performance of terahertz wave in high-density plasma sheath. The right-hand polarized terahertz wave will produce a power absorption peak near the cyclotron frequency due to cyclotron resonance. Changing the magnetization angle in a certain direction will bring an opposite effect on the transmission rate to left-hand polarized and right-hand polarized terahertz wave. Reducing the magnetization intensity can avoid the absorption peak of the right-hand polarized wave by the plasma to a certain extent. Increasing the magnetization can increase the transmission power of the left-hand polarized wave to a certain extent. Moreover, reducing the collision frequency can narrow the absorption band of the right-hand polarized wave in the plasma and increase the transmission power of left-hand polarized wave. In general, the transmission performance of left-hand polarized terahertz wave in non-uniformly magnetized plasma is better than that of right-hand polarized terahertz wave. These results provide a theoretical basis for investigating the blackout phenomenon. The adjusting of these parameters studied in the article is expected to be able to alleviate the blackout problem to a certain extent.
1. 引 言
高超音速飞行器再入大气过程中会与空气发生摩擦, 在产生大量热量的同时飞行器表面的空气会发生电离, 生成一层薄等离子体鞘套[1]. 飞行器的天线窗口被等离子体覆盖后, 其发射和接受信号的能力大大下降, 严重的时候通信会完全中断, 产生“黑障”现象[2].
为了解决“黑障”问题, 大量的科研人员对等离子体特性[3]和电磁波在等离子体中的传播行为进行了研究. 研究人员发现: 电磁波频率低于等离子体频率是导致“黑障”问题的重要原因之一, 因而就连高频段的毫米波也无法完全穿透飞行器表面形成的高电子密度等离子体鞘套. 近年来THz波技术不断进步, 其优秀的穿透性和抗干扰性使得“黑障”问题有了解决的可能, 因此, 研究THz波在等离子体鞘套中的传播有着重要意义. 另外, 研究人员开发了多种算法来分析电磁波与等离子体的相互作用, 如时域有限差分法(FDTD)[4-7]、解析法[8]、温泽尔-克莱默-布里渊法(WKB)[9,10]以及散射矩阵法(SMM)[11,12]等. 本文采用SMM方法进行研究分析, SMM方法较好地运用了分层思想来实现对等离子体不均匀性的体现, 能直接高效地计算出总的反射、透射和吸收率, 相较于FDTD方法也不需要进行大量的网格化操作, 所以对于飞行器表面等离子体鞘套的分析, SMM方法具有较高的效率和较低的算法复杂度. 正因如此, Chen等[13]用SMM方法研究了时变等离子体鞘层对雷达回波信号的影响, Guo和Guo[14]用SMM方法分析了THz波在非均匀磁化热等离子体中的传输特性, Zhang等[15]用SMM方法研究了THz波在右极化并且外加磁场角度与入射角度平行时, 在非均匀磁化等离子体中的传输特性. 我们发现, 这些研究大多只是基于固定的外加磁场角度对传输特性进行研究, 没有考虑外加磁场方向的变化对THz波传输的影响. 一些论文也只是基于单一极化方向的分析, 然而不同的极化方向也是影响THz波传输性能的因素. 所以, 为了更全面地探究磁化角度的改变对左极化和右极化THz波在非均匀磁化等离子体中传输的影响, 本文加入磁化角度作为分析传输特性的变量之一, 使用SMM方法分别研究了左极化和右极化THz波与等离子体之间不同的作用特点.
2. SMM理论
对于分布不均匀的等离子体可以设定为图1所示的分层模型, 将等离子体平板平均分成n层等离子体薄片, 每一薄片中的电子密度都视为均匀分布, 总的等离子体层电子密度按给定的密度函数分布. 电磁波垂直入射等离子板表面, 入射区(0)和透射区(p)都为自由空间.
在此模型中每一层薄片等离子体中都有其传播常数
k(m) :k(m)=ωc√ε(m)r, (1) 其中,
ε(m)r 为第m层薄片中的相对介电常数,c 为真空中的光速,ω 为入射电磁波角频率. 由于飞行器表面形成的等离子体在不同的高度或状态下有不同的形式[16], 十分复杂, 而冷等离子体模型在一些情况下能够近似描述等离子体鞘套, 这也是当前研究鞘套的主要模型之一[17,18]. 在这种情况下, 鞘套的相对介电常数用Appleton公式[19]表示:ε(m)r=1−ZX−Ysin2θ2(X−Z)±[Y2sin4θ4(X−Z)2+Ycos2θ]1/2, (2) 其中, 正负号代表左极化或者右极化,
θ 是波传播方向与磁化方向的夹角,X=1−iv/ω ,v 是碰撞角频率,Y=ωce2/ω2 ,ωce 是电子回旋角频率,Z= ω2p,m/ω2 ,ωp,m 是第m层的等离子体角频率并由(3)式表示[20]:ωp,m=√Nee2meε0, (3) 其中
Ne 为自由电子密度,e 和me 分别为电子电荷量和电子质量,ε0 为真空中的介电常数.将图1中入射区域(0)和薄层(1)的电磁场, 用入射(透射)电场
Eiz 、反射电场Esz 、入射(透射)磁场Hiy 和反射磁场Hsy 表示, 并分别写为以下形式:{Ei(0)z=E0e−ik(0)x,Es(0)z=AE0eik(0)x,Hi(0)y=−E0η0e−ik(0)x,Hi(0)y=AE0η0eik(0)x, (4) {Ei(1)z=B1E0e−ik(1)x,Es(1)z=C1E0eik(1)x,Hi(1)y=−B1E0η1e−ik(1)x,Hs(1)y=C1E0η1eik(1)x, (5) 其中
A 是射入到第一个等离子层表面的反射系数,B1 和C1 为电磁波从区域(0)射入等离子体薄层(1)的透射系数和薄层(2)反射回薄层(1)的反射系数,k(0) 是自由空间的波数. 这样可以确定入射区域(0)和薄层(1)的总电磁场关系, 通过匹配x=0 处分界面的边界条件能得到矩阵关系:(B1C1)=S1(A1), (6) 其中
S1=12k(1)x(k(1)x−k(0)xk(1)x+k(0)xk(1)x+k(0)xk(1)x−k(0)x). (7) 和上面相同, 继续匹配第2到n个分界面能得到矩阵关系:
(BmCm)=Sm(Bm−1Cm−1), (8) 其中
Sm=(e−ik(m)xdmeik(m)xdmk(m)xe−ik(m)xdm−k(m)xeik(m)xdm)−1×(e−ik(m−1)xdmeik(m−1)xdmk(m−1)xe−ik(m−1)xdm−k(m−1)xeik(m−1)xdm). (9) 在最后一个分界面的左右是薄层(n)和自由空间(p), 在区域(p)中只有透射电磁波, 没有反射电磁波, 所以匹配
x=d 处的边界条件能得到矩阵关系:(BnCn)=Vp⋅D, (10) 其中
Vp=12k(n)x((k(n)x+k(p)x)ei(k(n)x−k(p)x)dp(k(n)x−k(p)x)e−i(k(n)x+k(p)x)dp). (11) Sg(A1)=Vp⋅D (12) 其中
Sg=Sn⋅Sn−1⋅Sn−2⋯⋅S1=(2∏m=nSm)S1. (13) 接着将
Sg 表示为Sg=(Sg1,Sg2) , 其中Sg1 和Sg2 表示Sg 的第一列和第二列向量, 所以(12)式可以变换为(Sg1−Vp)(AD)=−Sg2. (14) 最后解得系数A和D的矩阵:
(AD)=−(Sg1−Vp)−1⋅Sg2. (15) 在接下来的分析中将电磁波功率归一化处理, 归一化反射功率表示为
|A|2 , 归一化透射功率表示为|D|2 , 那么归一化吸收功率为1−|A|2−|D|2 .3. 结果与分析
根据NASA钝头RAM飞行器的公开数据[16], 在不同高度下有不同的电子密度分布, 并且通常来说可以用指数分布、高斯分布、抛物分布等模型来模拟等离子体中的电子密度分布[21]. 但单一的分布模型不能模拟实际分布的非对称性, 所以本文采用双高斯分布模型来模拟计算以更好地贴近实际分布情况, 分布函数表示如(16)式, 函数图像如图2所示.
Ne={N0exp[−100(x−x0)2],0⩽ (16) 3.1 外加磁场方向对左极化和右极化THz波传输功率的影响
根据现有飞行器表面形成的等离子体鞘套的公开数据和流场计算[22], 在此设置等离子体平板的厚度为0.1 m, 为保证计算收敛性将其分成50层薄片, 碰撞频率
{f_v} = 0.1\;{\rm{ THz}} , 最大电子密度{N_0} = 1 \times {10^{18}}\;{{\rm{m}}^{ - 3}} , 外加磁场强度B = 5\;{\rm{ T}} . 根据(2)式, 可以对外加磁场的角度进行分析.由图3可以看出, 随着磁场角度由0°增加到90°时, 右极化THz波传输功率基本逐渐升高, 同时左极化THz波的传输功率基本逐渐降低. 由此可以看出, 当磁场方向与波的传播方向垂直时, 最有利于右极化太赫兹波信号的接收; 当磁场方向与波的传播方向平行时, 最有利于左极化太赫兹波信号的接收. 另外, 当电磁波在外加静态磁场的等离子体中传播时, 等离子体中的电子将受到磁场作用做回旋运动, 而右极化波会与回旋运动的电子同步旋转形成回旋共振, 导致能量吸收. 所以在相同磁化角度情况下, 左极化波和右极化波的传输功率曲线产生了不同的变化趋势.
那么以磁化角度分别为90°和0°来分析右极化太赫兹波和左极化太赫兹波在非均匀磁化等离子体中的传输特性. 当磁化角度
\theta = {90^ \circ } 且为右极化时相对介电常数如(17)式; 当磁化角度\theta = {0^ \circ } 且为左极化时相对介电常数如(18)式.\varepsilon _{\rm{r}}^{(m)} = 1 - \dfrac{{\omega _{{\rm{p}},m}^2/{\omega ^2}}}{{1 - {\rm{i}}\frac{v}{\omega } - \dfrac{{\omega _{{\rm{ce}}}^2/{\omega ^2}}}{{1 - \omega _{{\rm{p}},m}^2/{\omega ^2} - {\rm{i}}\,{v}/{\omega }}}}}, (17) \varepsilon _{\rm{r}}^{(m)} = 1 - \frac{{\omega _{{\rm{p}},m}^2}}{{{\omega ^2} - \omega ({\rm{i}}v + {\omega _{{\rm{ce}}}})}}. (18) 在这两种条件下就可以进一步分析在不同电子密度、碰撞频率和外加磁场强度条件下右极化和左极化THz波在等离子体中的传输特性.
3.2 电子密度对右极化和左极化THz波反射、传输、吸收功率的影响
图4(a)—图4(c)分别给出了右极化THz波在不同最大电子密度下的反射、透射和吸收功率. 其中等离子体平板的总厚度
d = 10\;{\rm{ cm}} , 碰撞频率{f_v} = 0.1\;{\rm{ THz}} , 外加磁场强度为B = 5\;{\rm{ T}} , 最大电子密度的值分别为{N_0} = {\rm{1}} \times {10^{1{\rm{6}}}}\;{{\rm{m}}^{ - 3}} ,{N_0} = {\rm{1}} \times {10^{1{\rm{7}}}}\;{{\rm{m}}^{ - 3}} ,{N_0} = {\rm{1}} \times {10^{1{\rm{8}}}}\;{{\rm{m}}^{ - 3}} 和{N_0} = {\rm{1}} \times {10^{1{\rm{9}}}}\;{{\rm{m}}^{ - 3}} . 可以看出, 较高电子密度的等离子体对太赫兹波的反射能力非常弱. 随着电子密度的增加, THz波传输功率相应降低, 且吸收功率增加, 在电子密度为{N_0} = 1 \times {10^{19}}\;{{\rm{m}}^{ - 3}} 时出现了传输阻带, 阻带内几乎所有的电磁波功率都被等离子体层吸收. 另外, 当电子密度{N_0} < 1 \times {10^{19}}\;{{\rm{m}}^{ - 3}} 时, 右极化电磁波的吸收功率在0—0.6\;{\rm{ THz}} 内先增加后减少, 在0.1—0.2\;{\rm{ THz}} 间出现一个突起峰值, 此时外加磁场条件下的电子回旋频率{f_{{\rm{ce}}}} = 0.14\;{\rm{ THz}} , 正好在峰值区域, 所以可以推断由于右旋波的偏振方向与电子回旋方向相同, 引起了电子回旋共振, 导致能量吸收增强, 并形成了突起峰值. 而电子密度的增加也引起了吸收峰值的升高, 这是因为电子密度的增加意味着等离子体中电子数量的增加, 从而导致电子从电磁波中吸收的能量增加. 图4(d)—图4(f) 分别给出了左极化THz波在不同最大电子密度下的反射、透射和吸收功率, 参数与右极化THz波保持一致. 可以看到的是, 由于左旋波的偏振方向与电子回旋方向不同, 左极化THz波相对于右极化THz波最大的差异在于没有了电子回旋共振的影响, 并且在相同的电子密度下左极化THz波对等离子体有更好的穿透性能.由于THz波在等离子体中的反射功率非常低, 传输功率和吸收功率的变化规律相反, 但传输功率在一定程度上能更好地体现信号传输的能力, 所以接下来主要对传输功率进行分析.
3.3 外加磁场强度对右极化和左极化THz波传输功率的影响
图5(a)和图5(b)以及图5(c)和图5(d)分别给出了外加磁场强度对右极化和左极化THz波传输功率的影响, 此时等离子体平板的总厚度
d = 10\;{\rm{ cm}} , 碰撞频率{f_v} = {\rm{0}}{\rm{.1 \;THz}} , 最大电子密度的值{N_0} = {\rm{1}} \times {10^{1{\rm{8}}}}\;{\rm{ }}{{\rm{m}}^{ - 3}} . 从图5(a)和图5(c)可以看出, 左极化波传输功率在整个频带上随外加磁场强度的增加而增加. 而右极化波传输功率随着外加磁场强度的增加在低频端和高频端有着不同的变化趋势: 在低频端右极化波传输功率随外加磁场强度的增加而增加, 但由于磁场强度的增加引起回旋共振的频带向高频端偏移, 进而在高频端产生更多的能量吸收, 所以在高频端右极化波传输功率会随外加磁场强度的增加而降低.3.4 碰撞频率对右极化和左极化THz波传输功率的影响
图6(a)和图6(b)以及图6(c)和图6(d)分别给出了碰撞频率对右极化和左极化THz波传输功率的影响, 此时等离子体平板的总厚度
d = 10\;{\rm{ cm}} , 外加磁场强度为B = {\rm{5 T}} , 最大电子密度的值{N_0} = {\rm{1}} \times {10^{1{\rm{8}}}}\;{{\rm{m}}^{ - 3}} . 可以看出, 右极化波传输功率的最小值随碰撞频率的增加而减小, 但吸收频带增大. 左极化波传输功率随着碰撞频率的增大而减少. 总体来说左极化波更能穿透磁化冷等离子体, 有助于更好地改善飞行器再入“黑障”问题.4. 结 论
本文以NASA钝头RAM实验数据为基础, 研究了左极化和右极化太赫兹波在非均匀磁化冷等离子体中的传输特性. 结果表明, 外加磁场的角度对左右极化太赫兹波的影响有所不同, 在磁场方向与波传输方向垂直时最有利于右极化太赫兹波的传输; 在磁场方向与波传播方向平行时最有利于左极化太赫兹波的传输. 另外, 随着等离子体最大电子密度的增加, 等离子体对太赫兹波的能量吸收增强, 电磁波在右极化时会发生回旋共振, 在回旋频率附近产生吸收峰, 并且外加磁场的强度会改变吸收峰在频谱上的左右位置, 此外碰撞频率大小也会影响太赫兹波的传输. 总的来说左极化太赫兹波在非均匀磁化等离子体中的传输性能要优于右极化太赫兹波. 这些结果可为改善和解决再入飞行“黑障”问题提供理论参考.
[1] Gupta R N, Yos J M, Thompson R A, Lee K P 1990 A Review of Reaction Rates and Thermodynamic and Transport Properties for an 11-species Air Model for Chemical and Thermal Nonequilibrium Calculations to 30000 K (Hampton: Langley Research Center) NASA-RP-1232
[2] 姚博 2019 博士学位论文 (西安: 西安电子科技大学)
Yao B 2019 Ph. D. Dissertation (Xi’an: Xidian University) (in Chinese)
[3] 杨楠, 杜海伟 2014 红外与毫米波学报 33 237
Google Scholar
Yang N, Du H W 2014 J. Infrared Millimeter Waves 33 237
Google Scholar
[4] Huang S J, Li F 2004 Int. J. Infrared Millimeter Waves 25 815
Google Scholar
[5] 陈伟, 郭立新, 李江挺, 淡荔 2017 物理学报 66 084102
Google Scholar
Chen W, Guo L X, Li J T, Dan L 2017 Acta Phys. Sin. 66 084102
Google Scholar
[6] Chen W, Guo L X, Li J T, Liu S H 2016 IEEE Trans. Plasma Sci. 44 3235
Google Scholar
[7] Cheng G X, Liu L 2010 IEEE Trans. Plasma Sci. 38 3109
Google Scholar
[8] Jazi B, Rahmani Z, Shokri B 2013 IEEE Trans. Plasma Sci. 41 290
Google Scholar
[9] 曹建章, 李景镇, 陈国瑞 2002 电波科学学报 17 125
Google Scholar
Cao J Z, Li J Z, Chen G R 2002 Chin. J. Radio Sci. 17 125
Google Scholar
[10] 林敏, 徐浩军, 魏小龙, 梁华, 张艳华 2015 物理学报 64 055201
Google Scholar
Lin M, Xu H J, Wei X L, Liang H, Zhang Y H 2015 Acta Phys. Sin. 64 055201
Google Scholar
[11] Helaly A, Soliman E A, Megahed A A 1997 IEE Proc. MIicrow. Antennas Propag. 144 61
Google Scholar
[12] Hu B J, Wei G, Lai S L 1999 IEEE Trans. Plasma Sci. 27 1131
Google Scholar
[13] Chen X, Li K, Liu Y, Zhou Y, Li X, Liu Y 2017 IEEE Trans. Plasma Sci. 45 3166
Google Scholar
[14] Guo L X, Guo L J 2017 Phys. Plasmas 24 112119
Google Scholar
[15] Zhang Y Y, Xu G J, Zheng Z Q 2019 Optik 182 618
Google Scholar
[16] Dunn M G, Kang S W 1973 Theoretical and Experimental Studies of Reentry Plasmas (Washington: National Aeronautics and Space Administration) NASA-CR-2232
[17] Liu J F, Xi X L, Wang L L 2011 IEEE Trans. Plasma Sci. 39 852
Google Scholar
[18] Tian Y X, Yan W Z, Gu X L, Jin X L, Li J Q, Li B 2017 AIP Adv. 7 125325
Google Scholar
[19] Heald M A, Wharton C B, Furth H P 1965 Phys. Today 18 72
Google Scholar
[20] Yeh C, Rusch W V T 1965 J. Appl. Phys. 36 2302
Google Scholar
[21] 孙朋飞 2019 硕士学位论文 (西安: 西安电子科技大学)
Sun P F 2019 M. S. Thesis (Xi’an: Xidian University) (in Chinese)
[22] 薄勇, 赵青, 罗先刚, 刘颖, 陈禹旭, 刘建卫 2016 物理学报 65 035201
Google Scholar
Bo Y, Zhao Q, Luo X G, Liu Y, Chen Y X, Liu J W 2016 Acta Phys. Sin. 65 035201
Google Scholar
-
-
[1] Gupta R N, Yos J M, Thompson R A, Lee K P 1990 A Review of Reaction Rates and Thermodynamic and Transport Properties for an 11-species Air Model for Chemical and Thermal Nonequilibrium Calculations to 30000 K (Hampton: Langley Research Center) NASA-RP-1232
[2] 姚博 2019 博士学位论文 (西安: 西安电子科技大学)
Yao B 2019 Ph. D. Dissertation (Xi’an: Xidian University) (in Chinese)
[3] 杨楠, 杜海伟 2014 红外与毫米波学报 33 237
Google Scholar
Yang N, Du H W 2014 J. Infrared Millimeter Waves 33 237
Google Scholar
[4] Huang S J, Li F 2004 Int. J. Infrared Millimeter Waves 25 815
Google Scholar
[5] 陈伟, 郭立新, 李江挺, 淡荔 2017 物理学报 66 084102
Google Scholar
Chen W, Guo L X, Li J T, Dan L 2017 Acta Phys. Sin. 66 084102
Google Scholar
[6] Chen W, Guo L X, Li J T, Liu S H 2016 IEEE Trans. Plasma Sci. 44 3235
Google Scholar
[7] Cheng G X, Liu L 2010 IEEE Trans. Plasma Sci. 38 3109
Google Scholar
[8] Jazi B, Rahmani Z, Shokri B 2013 IEEE Trans. Plasma Sci. 41 290
Google Scholar
[9] 曹建章, 李景镇, 陈国瑞 2002 电波科学学报 17 125
Google Scholar
Cao J Z, Li J Z, Chen G R 2002 Chin. J. Radio Sci. 17 125
Google Scholar
[10] 林敏, 徐浩军, 魏小龙, 梁华, 张艳华 2015 物理学报 64 055201
Google Scholar
Lin M, Xu H J, Wei X L, Liang H, Zhang Y H 2015 Acta Phys. Sin. 64 055201
Google Scholar
[11] Helaly A, Soliman E A, Megahed A A 1997 IEE Proc. MIicrow. Antennas Propag. 144 61
Google Scholar
[12] Hu B J, Wei G, Lai S L 1999 IEEE Trans. Plasma Sci. 27 1131
Google Scholar
[13] Chen X, Li K, Liu Y, Zhou Y, Li X, Liu Y 2017 IEEE Trans. Plasma Sci. 45 3166
Google Scholar
[14] Guo L X, Guo L J 2017 Phys. Plasmas 24 112119
Google Scholar
[15] Zhang Y Y, Xu G J, Zheng Z Q 2019 Optik 182 618
Google Scholar
[16] Dunn M G, Kang S W 1973 Theoretical and Experimental Studies of Reentry Plasmas (Washington: National Aeronautics and Space Administration) NASA-CR-2232
[17] Liu J F, Xi X L, Wang L L 2011 IEEE Trans. Plasma Sci. 39 852
Google Scholar
[18] Tian Y X, Yan W Z, Gu X L, Jin X L, Li J Q, Li B 2017 AIP Adv. 7 125325
Google Scholar
[19] Heald M A, Wharton C B, Furth H P 1965 Phys. Today 18 72
Google Scholar
[20] Yeh C, Rusch W V T 1965 J. Appl. Phys. 36 2302
Google Scholar
[21] 孙朋飞 2019 硕士学位论文 (西安: 西安电子科技大学)
Sun P F 2019 M. S. Thesis (Xi’an: Xidian University) (in Chinese)
[22] 薄勇, 赵青, 罗先刚, 刘颖, 陈禹旭, 刘建卫 2016 物理学报 65 035201
Google Scholar
Bo Y, Zhao Q, Luo X G, Liu Y, Chen Y X, Liu J W 2016 Acta Phys. Sin. 65 035201
Google Scholar
期刊类型引用(3)
1. 石玉仁,万晓欢,赵月星,张娟,樊小贝. 太赫兹波斜入射高温磁化等离子体的解析法研究. 西北师范大学学报(自然科学版). 2022(01): 44-50 . 百度学术
2. 蔡丽君,袁英豪. 电子密度对太赫兹波在均匀非磁化等离子体中的传输影响. 光学仪器. 2022(01): 80-86 . 百度学术
3. 张洁,杨光华. 等离子体随机分布对THz波传输特性的影响. 西北师范大学学报(自然科学版). 2022(04): 58-65 . 百度学术
其他类型引用(6)
计量
- 文章访问数: 6029
- PDF下载量: 81
- 被引次数: 9