Loading [MathJax]/jax/output/HTML-CSS/jax.js

搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

水中微小波纹圆柱体声散射低频共振调控

周彦玲 王斌 范军

SUN Kangsheng, HAN Chao, QIN Haifeng, GU Tao, LI Wei, YU Cheng. Fast generation algorithm of high-quality holograms based on attention convolutional neural network. Acta Phys. Sin.. doi: 10.7498/aps.74.20241713
Citation: SUN Kangsheng, HAN Chao, QIN Haifeng, GU Tao, LI Wei, YU Cheng. Fast generation algorithm of high-quality holograms based on attention convolutional neural network. Acta Phys. Sin.. doi: 10.7498/aps.74.20241713

水中微小波纹圆柱体声散射低频共振调控

周彦玲, 王斌, 范军

Manipulation of low-frequency resonance scattering from submerged slightly deformed polymer cylinders

Zhou Yan-Ling, Wang Bin, Fan Jun
Article Text (iFLYTEK Translation)
PDF
HTML
导出引用
  • 塑料类高分子材料甲基丙烯酸甲酯-亚克力(PMMA)圆柱中亚音速Rayleigh波低频隧穿共振可引起反向散射增强, 在低频标准散射体设计等领域具有重要应用价值. 提出一种微弱形变的规则波纹表面结构, 可实现水中PMMA圆柱反向散射低频共振频率的无源调控. 利用微扰法推导了水中微弱形变规则波纹圆柱反向散射低频共振频率偏移的近似解, 讨论了波纹微扰系数、周期对规则波纹圆柱共振频率偏移的影响规律. 基于Rayleigh波相位匹配方法分析了低频共振频率偏移的机理. 研究表明: 微弱形变规则波纹圆柱中亚音速Rayleigh波沿微弱形变波纹表面传播, 与光滑圆柱体相比, 传播路径的改变引起Rayleigh波传播相位变化, 导致了Rayleigh波低频共振频率发生偏移. 最后开展了微弱形变规则波纹圆柱体声散射特性水池实验, 获取了其反向散射共振频率, 明显观察到了规则波纹圆柱共振频率偏移现象, 与理论预报结果吻合较好.
    Backscattering enhancement associated with subsonic Rayleigh wave resonance of a polymethlmethacrylate (PMMA) cylinder is observed at low frequencies in water, which suggests that the PMMA cylinders may have essential applications in the low-frequency standard scatterer design. A slightly deformed surface structure with regular corrugation is presented to manipulate the low-frequency backscattering resonance frequency of PMMA cylinder passively. Using the perturbation method, the approximate resonance frequency shift is derived for an infinite slightly deformed cylinder with regular corrugations. Based on the phase matching of Rayleigh waves, the mechanism of low-frequency resonance frequency shift is revealed. Studies show that compared with a bare cylinder, a small boundary deformation can cause the propagation path of Rayleigh waves to change, namely, the Rayleigh waves propagate along the slightly deformed surface with regular corrugations. The modified propagation path can cause the propagation phase to change, which brings about the low-frequency backscattering resonance frequency shift of a PMMA cylinder. Furthermore, how the resonance frequency shifts with the deformation coefficient and period of the corrugation is discussed in detail. The backscattering resonance frequency of the regular corrugated cylinder shifts to low frequency with the increase of ξ under the condition of the deformation coefficient ξ > 0, but shifts to high frequency with the decrease of ξ at ξ < 0, and the resonance frequency shift increases with the increase of absolute value of deformation coefficient. When corrugation period m < 10, the phase variation with m is too small, so the corrugation period has little effect on the resonance frequency shift. Finally, acoustic scattering experiment of the regular corrugated cylinder is conducted in the tank. The resonance frequency shift is obviously observed in the experiment, which is in good agreement with the theoretical prediction. Hence, the characteristics of backscattering enhancement associated with subsonic Rayleigh wave and the resonance frequency shifts make the PMMA deformed cylinder have potential applications such as in standard scatter design and identification using “AcoustiCode”.
      通信作者: 王斌, bin_wang@sjtu.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 11774229)资助的课题
      Corresponding author: Wang Bin, bin_wang@sjtu.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant No. 11774229)

    在光学领域, Ge等[1,2]利用微扰理论研究了边界微弱形变光学微腔的共振特性, 研究表明微腔边界形变导致不同共振模态之间耦合, 使得其共振频率发生偏移, 并在很大程度改变腔内射线动力学特性和输出方向, 且光波的输出方向对边界变形极其灵敏, 在光束控制和高分辨率检测方面具有重要应用前景. 在声学领域, Fawcett[3]用数值方法计算了外径不变内径微弱形变水中球壳声散射频域特征, 结果显示边界的微弱形变可以使得由于球壳中传播的a0波与刚性背景耦合形成的中频吻合共振[4]得到抑制, 但未对抑制中频吻合共振的机理进行研究; 提出了一种时域Kirchhoff/衍射混合算法[5], 此算法只适用于微弱形变刚性目标高频声散射的计算. Antonio等[6]提出一种数值曲面积分方法, 可以用来计算不规则表面或随机粗糙表面几何目标的声散射问题, 但此方法只适用于刚性目标. 本文基于微扰法分析了微弱形变弹性圆柱体低频声散射特征, 提出了一种利用微弱形变规则波纹表面结构调控塑料类高分子材料甲基丙烯酸甲酯-亚克力(PMMA)圆柱反向散射低频共振频率的方法.

    塑料类高分子聚合物材料PMMA是3D打印领域的一类主要材料, 其密度与水接近, 剪切波速度小于水中声速[7], 近年来在声学领域逐渐受到关注. Hefner和Marston[8]通过理论和实验研究了水中PMMA球体低频反向散射增强现象, 研究表明由于亚音速Rayleigh波共振使得PMMA球体目标低频反向散射大幅度增强, 并指出这类PMMA目标在水中低频标准散射体设计等方面具有重要应用价值. PMMA目标低频共振频率的人工调控, 是设计满足不同频率需求的声散射标准体, 是水中被动声学标记体的基础[9,10]. 本文提出了一种应用于水中PMMA圆柱体低频共振频率调控的规则波纹表面, 第2节基于微扰法推导了无限长微弱形变规则波纹圆柱体声散射共振频率偏移的近似解, 分析了其低频共振频率偏移的机理, 即规则波纹表面导致圆柱中亚音速Rayleigh波传播路径改变, 从而引起传播相位变化, 使得低频共振频率发生偏移. 第3节讨论波纹微扰系数、周期等波纹参数对水中PMMA规则波纹圆柱低频共振散射频率偏移的影响规律. 第4节进行了光滑PMMA圆柱体和规则波纹PMMA圆柱体(相同波纹周期、不同波纹微扰系数)声散射特性水池试验, 获取了其反向散射共振特性, 观察到共振频率偏移特征, 同时验证了第2节规则波纹圆柱共振频率偏移近似公式的有效性和正确性.

    无限长微弱形变规则波纹表面圆柱体示意如图1所示.

    图 1 无限长微弱形变规则波纹表面圆柱体\r\nFig. 1. Infinite regular-corrugated cylinder.
    图 1  无限长微弱形变规则波纹表面圆柱体
    Fig. 1.  Infinite regular-corrugated cylinder.

    若圆柱表面存在微小变形[3], 即

    R(a,φ)=a+εf(φ), (1)

    其中, a为未变形圆柱半径, ε为波纹微扰系数, f(φ)=sin2(mφ),m=0,1,2,3,, m[0,π]的波纹周期. 则波纹圆柱圆周长度S

    S(ε,m)=2π0[a+εsin2(mφ)]2+[2mεsin(mφ)cos(mφ)]2dφ. (2)

    由(2)式可以看出, 微弱形变规则波纹圆柱圆周长度与未变形圆柱体半径a、波纹微扰系数ε和波纹周期m有关.

    假设|εf(φ)|a, 入射平面波沿着x轴正向垂直入射到无限长波纹圆柱, 设散射声场具有如下形式[11,12]:

    pscat(r,φ)=n=0εnin[bnH(1)n(kr)cos(nφ)+εpnαbpH(1)p(kr)cos(pφ)], (3)

    这里忽略时间因子e–iωt, 其中k = ω/c是入射平面波波数, ω为角频率, c为水中声速, r为接收点到目标中心的距离, φ为入射平面波与x轴正向的夹角, εn={1,n=0,2,else, bn是散射系数, H(1)n(kr)为第一类汉克尔函数. 则微弱形变规则波纹圆柱外部总声场为

    pt(r,φ)=n=0εnin{[Jn(kr)+bnH(1)n(kr)]cos(nφ)+εpnαbpH(1)p(kr)cos(pφ)}. (4)

    微弱形变规则波纹圆柱内部声场:

    Φ(r,φ)=n=0εnin[cnJn(kdr)cos(nφ)+εpnαcpJp(kdr)cos(pφ)],Ψ(r,φ)=n=0εnin[enJn(ksr)sin(nφ)+εpnαepJp(ksr)sin(pφ)], (5)

    其中, Jnn阶Bessel函数, αbp, αcpαep是待定的耦合系数, kdks分别为纵波波数和横波波数, p为耦合阶次.

    根据边界条件:

    τrr(a+εf(φ),φ)=pt(a+εf(φ),φ),
    ur(2)(a+εf(φ),φ)=ur(1)(a+εf(φ),φ),
    τrφ(a+εf(φ),φ)=0,

    其中, τrr是结构中的法向应力, u(1)ru(2)r分别为结构和水中的法向位移, τrφ是结构中的切向应力. 基于微扰理论将(6)式按照ε的一阶泰勒级数展开, 则有

    τrr(a,φ)+pt(a,φ)=εf(φ)[τrrr+ptr](a,φ),
    u(2)r(a,φ)u(1)r(a,φ)=εf(φ)[u(2)rru(1)rr](a,φ),
    τrφ(a,φ)=εf(φ)τrφr(a,φ),

    由(7)式, 得

    [Dd11Dd12Dd13000Dd21Dd22Dd230000Dd32Dd33000Dd41Dd42Dd43Dd44Dd45Dd46Dd51Dd52Dd53Dd54Dd55Dd560Dd62Dd630Dd65Dd66]×(bncnenαbpαcpαep)=(A1A20A4A50), (8)

    根据克莱姆法则得到:

    bn=BdnDdn,αbp=AbpnDdn, (9)

    其中Bdn, AbpnDdn为6 × 6的行列式, 行列式具体元素见附录A.

    无限长规则波纹圆柱远场反向散射形态函数[13]表示为

    f(x)=2iπxn=0εn[Bdn(x)Ddn(x)(1)n+εpnAbpn(x)Ddn(x)(1)p]. (10)

    x=kR为无限长微弱形变规则波纹圆柱特征方程Ddn(x)=0的根. 当ε=0时, (10)式为无限长光滑圆柱远场反向散射形态函数, 此时通过文献[14]的方法求解特征方程得到无限长光滑圆柱特征值的实部, 记为x0. 根据共振散射理论[15], 无限长圆柱特征值的实部对应共振峰的频率, 则无限长圆柱共振频率f0=x0c2πa, 对应的相速度[16]:

    cph=2πf0a/n,n=0,1,2,3,. (11)

    ε0时, 无限长微弱形变规则波纹圆柱共振频率记为f. 将xx0附近按照ε的一阶展开为

    x=x0+εx1. (12)

    将(6)式在x0附近按照ε的一阶展开, 具体推导过程见附录A, 经过一系列数学化简得

    x1=x0Ann, (13)

    其中

    Ann=εnππ0f(φ)acos(nφ)cos(nφ)dφ.

    则微弱形变规则波纹圆柱特征值偏移量为

    Δx=xx0=εx0Ann. (14)

    对应的微弱形变规则波纹圆柱共振频率偏移量可表示为

    Δf=ff0=εf0Ann. (15)

    无限长微弱形变规则波纹圆柱共振频率近似为

    f=f0εf0Ann. (16)

    定义微弱形变规则波纹圆柱归一化波纹微扰系数ξ = ε/a, 利用(10)式计算半径a = 0.05 m材料为PMMA的无限长光滑圆柱和波纹参数m = 6, ξ = –3%的无限长微弱形变规则波纹圆柱反向散射形态函数幅频特性, 所用材料参数见表1. 计算频率范围为1—18 kHz, 频率间隔为5 Hz, 所对应无因次频率ka的范围为0.20—4.19, ka较小, 文中所研究频段为低频.

    表 1  计算所用材料参数
    Table 1.  Material parameters used in the calculations.
    材料密度/(kg·m–3)纵波波速/(m·s–1)剪切波波速/(m·s–1)
    PMMA119026901340
    10001500
    下载: 导出CSV 
    | 显示表格

    图2(a)对比了无限长光滑圆柱(黑色实线)和微弱形变规则波纹圆柱(红色点划线)反向散射形态函数幅频特性, 研究发现边界微弱变形使得PMMA圆柱共振频率发生了偏移. 图2(b)图2(a)蓝色虚线框中共振峰的局部放大图, 无限长光滑圆柱对应的共振峰频率f0 = 6745 Hz, 微弱形变规则波纹圆柱的共振峰频率f1 = 6845 Hz与利用(16)式计算的微弱形规则波纹圆柱共振频率f2 = 6846.2 Hz基本一致. 因此, 可以用(16)式近似估计无限长微弱形变规则波纹圆柱的共振频率.

    图 2 基于微扰法无限长规则波纹圆柱形态函数幅频特性 (a)反向散射; (b)局部放大\r\nFig. 2. Form function of the infinite regular-corrugated cylinder based on the perturbation method: (a) Backscattering; (b) local enlargement.
    图 2  基于微扰法无限长规则波纹圆柱形态函数幅频特性 (a)反向散射; (b)局部放大
    Fig. 2.  Form function of the infinite regular-corrugated cylinder based on the perturbation method: (a) Backscattering; (b) local enlargement.

    根据(11)式计算得到无限长PMMA圆柱Rayleigh波相速度频散曲线, 如图3所示. 从图3可以看出, Rayleigh波相速度(黑色实线)随着频率的增大逐渐趋于常数cR1060 m/s, 小于水中声速, 为亚音速波. 文献[8,17]利用射线近似理论详细描述了PMMA球体中亚音速Rayleigh波对低频反向散射增强的贡献及传播路径. 对于平面波垂直入射无限长PMMA圆柱, Rayleigh波沿圆柱表面周向传播, 如图4(a)所示. 亚音速Rayleigh波在临界角θl=π/2B1点被激发, 在B1, B2点(箭头处)耦合沿圆柱表面周向传播, D1, D2点处耦合从D1点出射进入散射声场. 在焦散线上, Rayleigh波的投影切向速度等于水中的声速. 声场从渐消到传播的转变发生在半径b=ac/cR附近. 在a<r<b区域为渐消波区域, Rayleigh波通过耦合隧穿此渐消波区域再辐射形成反向散射声场. 因为入射平面波对称激发, 所以环绕波成对出现, 它们的绕行方向相反, 从D1B1方向传播的波这里没有显示.

    图 3 PMMA圆柱中Rayleigh波相速度频散曲线(黑色线), 曲线$2{\text{π}}a{f_0}/M$(蓝色点划线)和曲线$Sf/M$(红色虚线)\r\nFig. 3. Dispersion curve of subsonic Rayleigh waves (black line), curve $2{\text{π}}a{f_0}/M$ (blue dashed-dotted line) and curve $Sf/M$(red dotted line).
    图 3  PMMA圆柱中Rayleigh波相速度频散曲线(黑色线), 曲线2πaf0/M(蓝色点划线)和曲线Sf/M(红色虚线)
    Fig. 3.  Dispersion curve of subsonic Rayleigh waves (black line), curve 2πaf0/M (blue dashed-dotted line) and curve Sf/M(red dotted line).
    图 4 亚音速Rayleigh波传播路径 (a)光滑圆柱; (b)规则波纹圆柱\r\nFig. 4. Ray diagram for subsonic Rayleigh waves propagating around: (a) Bare cylinder; (b) regular-corrugated cylinder.
    图 4  亚音速Rayleigh波传播路径 (a)光滑圆柱; (b)规则波纹圆柱
    Fig. 4.  Ray diagram for subsonic Rayleigh waves propagating around: (a) Bare cylinder; (b) regular-corrugated cylinder.

    当圆柱表面具有微弱形变时, Rayleigh波传播过程与光滑圆柱相同, 但传播路径发生改变, 如图4(b)所示. Rayleigh波沿规则波纹表面传播, 由微弱形变表面所确定的焦散半径bd=[a+εsin2(mφ)]c/cR发生变化. Rayleigh波传播路径的改变引起传播相位变化, 由环绕波分析理论[18]可知, Rayleigh波沿规则波纹圆柱表面周向传播的圆周数v=0,1,, 在接收点形成程差为S的等间隔波列, 这些波列的相位差为

    Δφ(f,ξ,m)=ωS(ξ,m)cR=2πfS(ξ,m)cR. (17)

    若这些波列的相位差是2πM倍(M为整数), 同相叠加形成共振, 由此得到共振条件

    Δφ=2Mπ,M=1,2,3,, (18)

    由(17)式、(18)式得

    cR=S(ξ,m)Mf. (19)

    ε=0时, cR=2πaf0/M, 无限长PMMA圆柱中Rayleigh波频散曲线与曲线2πaf0/M(图3中蓝色点划线)的交点所对应的横坐标为f0, 无限长PMMA圆柱的共振频率为f0, 见表2. 当ε0时, 利用(16)式求得的规则波纹圆柱共振频率为fξ, 无限长PMMA圆柱中Rayleigh波频散曲线与曲线Sf/M(图3中红色虚线)的交点所对应的横坐标为fξ, 这里m = 6, ξ=5%. 定义由共振条件确定的共振频率相对误差为δf=|ff|/f×100%, 由图5可知对于光滑圆柱和规则波纹圆柱, 由共振条件确定的共振频率相对误差均在0.3%以下. 因此, 利用Rayleigh波共振条件可以近似估计无限长PMMA圆柱反向散射低频共振频率. 同时, 这也证明了规则波纹PMMA圆柱中Rayleigh波速度不变, 但Rayleigh波传播路径发生改变引起传播相位的变化, 使得PMMA圆柱低频共振频率发生偏移, 从而实现了对PMMA圆柱反向散射低频共振频率的调控.

    表 2  无限长PMMA圆柱和规则波纹圆柱共振频率
    Table 2.  Resonance frequency of PMMA bare and corrugated cylinder
    lf0/Hzf0/Hzfξ /Hzfξ/Hz
    16419642062596240
    29953995597049700
    313383133851304813050
    416790167901637016375
    下载: 导出CSV 
    | 显示表格
    图 5 不同方法获取共振频率误差\r\nFig. 5. The relative error of resonance frequencies obtained by different methods.
    图 5  不同方法获取共振频率误差
    Fig. 5.  The relative error of resonance frequencies obtained by different methods.

    以上研究表明, 边界微弱形变可以使得PMMA圆柱反向散射低频共振频率发生偏移. 由(16)式很直观地看出无限长规则波纹圆柱反向散射低频共振频率与波纹微扰系数、波纹周期有关, 基于微扰法分析边界变形结构参数如波纹微扰系数、波纹周期对无限长PMMA圆柱反向散射低频共振频率的影响. 图6(a)为规则波纹圆柱反向散射波纹微扰系数-频率谱, 颜色表示形态函数幅值, 横轴是频率, 单位为kHz, 纵轴是归一化波纹微扰系数ξ范围–10%—10%, 间隔1%. 图中黑色虚线为根据(16)式获取的各个共振峰频率随波纹微扰系数变化曲线, 当归一化波纹微扰系数ξ > 0时, 波纹在[a,a+ε]区域内起伏, 规则波纹圆柱反向散射共振频率随着ξ的增大向低频偏移, 且频率偏移量随着ξ的增大而增大; 当归一化波纹微扰系数ξ < 0时, 波纹在[a+ε,a]区域内起伏, 规则波纹圆柱反向散射共振频率随着ξ的减小向高频偏移, 且频率偏移量随着ξ的减小而增大. 因此可以通过改变波纹微扰系数对水中PMMA圆柱反向散射低频共振频率进行无源调控. 图6(b)为规则波纹圆柱反向散射频率-波纹周期谱, 颜色表示形态函数幅值, 横轴是频率, 单位kHz, 纵轴是波纹周期m范围2—6, 这里ξ = –3%. 图中黑色虚线为根据(16)式获取的共振峰频率随波纹周期变化曲线, 可见当波纹周期比较小时, 波纹周期对无限长规则波纹圆柱反向散射共振频率偏移影响较小. 不同波纹周期所对应相位变化量为ΔΦm=Δφ(f0,ξ,m)Δφ(f0,ξ,1), 图7给出共振频率f0 = 6419 Hz, ξ = –3%, m取值0—20所对应的相位变化量ΔΦm, 可见当m < 10时, m变化引起的相位变化较小, 不足以对共振频率偏移产生重要作用.

    图 6 (a)频率-波纹微扰系数谱; (b)频率-波纹周期谱\r\nFig. 6. (a) Frequency-height spectra; (b) frequency-period spectra.
    图 6  (a)频率-波纹微扰系数谱; (b)频率-波纹周期谱
    Fig. 6.  (a) Frequency-height spectra; (b) frequency-period spectra.
    图 7 不同波纹周期对应相位变化量\r\nFig. 7. Phase varying with corrugated period.
    图 7  不同波纹周期对应相位变化量
    Fig. 7.  Phase varying with corrugated period.

    针对PMMA规则波纹圆柱开展声散射特性水池实验, 此次实验是在中国科学院声学研究所北海站消声水池完成. 3个实验模型(一个光滑圆柱, 两个规则波纹圆柱, 波纹周期m = 6, 归一化波纹微扰系数分别为ξ = –3%和ξ = –10%)均为亚克力材料, 密度ρ = 1150 kg/m3, 半径a = 0.05 m, 长度L = 0.5 m. 为了方便吊放, 模型两端各打4个孔, 实验模型垂直置于水池中绕中轴线转动. 发射换能器采用刚性连接吊放位置距离目标3.65 m, 水听器型号为B&K8103, 吊放位置距离目标2.11 m. 发射换能器、水听器和目标三者处于同一深度, 距离水面3.87 m, 实验模型及具体布放如图8图9所示. 为排除市电干扰, 测量过程中所有设备全程处于UPS直流供电状态.

    图 8 实验模型\r\nFig. 8. Experimental objects.
    图 8  实验模型
    Fig. 8.  Experimental objects.
    图 9 实验布放\r\nFig. 9. Diagram of experimental system setup.
    图 9  实验布放
    Fig. 9.  Diagram of experimental system setup.

    发射信号为线性调频信号, 频率为5—15 kHz, 脉宽0.8 ms, 发射周期为500 ms. 将水听器接收到的时域信号进行宽带信号处理[17], 截取直达波(即发射阵直接到达水听器的信号)时域信号pi(t)和回波时域信号ps(t)进行傅里叶变换得到其频谱Pi(f)Ps(f), 再将回波信号频谱除以直达波信号频谱并对距离进行修正获得3个实验模型反向散射形态函数幅频特性曲线, 如图10所示. 实验中在所研究频段明显观察到圆柱反向散射形态函数3个共振峰(黑色线), 微弱形变规则波纹圆柱反向散射共振峰明显的向高频偏移, 且ξ = –10%的变形圆柱(蓝色虚线)共振峰频率偏移量大于ξ = –3%的变形圆柱(红色点划线)共振峰偏移量, 与第3节理论分析结果一致.

    图 10 光滑圆柱和规则波纹表面圆柱反向散射形态函数实验结果\r\nFig. 10. Backscattering form function of regular-corrugated cylinders in the experiment.
    图 10  光滑圆柱和规则波纹表面圆柱反向散射形态函数实验结果
    Fig. 10.  Backscattering form function of regular-corrugated cylinders in the experiment.

    图10中3个实验模型的反向散射形态函数中提取共振峰频率, 光滑圆柱共振峰频率记为f0, 规则波纹圆柱共振峰频率表示为fξ, 利用(16)式获取的共振峰频率为fξ, 详细数值见表3. 这里定义实验获取的微弱形变规则波纹圆柱共振峰频率相对误差为

    表 3  实验模型共振峰频率
    Table 3.  Resonance frequencies of experimental objects.
    lf0/Hzfξ1/Hzfξ1 /Hzfξ2/Hzfξ2/Hz
    166496748.768526981.46966
    286438772.690989075.19294
    31241012596.01299613031.013192
    下载: 导出CSV 
    | 显示表格
    δξ=|fξlfξl|fξl×100%. (20)

    为了更为直观地观察各个共振频率, 以图11(a), (b)柱状图形式呈现, 分别对应归一化波纹微扰系数ξ = –3%和ξ = –10%规则波纹的共振频率. 实验获取微弱形变规则波纹圆柱共振峰频率与近似(16)式所估算的结果吻合较好, 相对误差均在5%以内(见图11(c)). 因此, 可以用(16)式快速估计微弱形变规则波纹圆柱的共振频率, 从而实现对PMMA圆柱反向散射低频共振频率的调控.

    图 11 实验和(16)式获取得的共振频率对比 (a) ξ = –3%; (b) ξ = –10%; (c)相对误差\r\nFig. 11. Resonance frequencies obtained by experiment and Eq. (16): (a) ξ = –3%; (b) ξ = –10%; (c) relative error.
    图 11  实验和(16)式获取得的共振频率对比 (a) ξ = –3%; (b) ξ = –10%; (c)相对误差
    Fig. 11.  Resonance frequencies obtained by experiment and Eq. (16): (a) ξ = –3%; (b) ξ = –10%; (c) relative error.

    图12为实验获取的微弱形变规则波纹圆柱的频率-角度谱, 横轴表示入射角度, 纵轴表示频率, 单位为kHz, 颜色表示形态函数幅值, 黑色虚线表示共振峰频率. 从图12可以看出: 微弱形变规则波纹圆柱低频共振频率在各个角度都发生了偏移, 且在各个方位角偏移量基本相同. 主要的原因是归一化微扰系数ξ = –3%和ξ = –10%所对应的微扰系数分别为ε = 0.0015 m和ε = 0.005 m, 在所研究频段内Rayleigh波波长λ > 0.07 m (所对应频率为f = 15 kHz), 此时Rayleigh波波长远大于波纹微扰系数, 对散射起主要作用的仍然是低频Rayleigh波共振, 波纹多镜面反射的影响尚未清晰显示出来[19], 因此微弱形变规则波纹PMMA圆柱低频共振频率偏移量与入射方向无关.

    图 12 规则波纹圆柱频率-角度谱 (a) ξ = 0; (b) ξ = –3%; (c) ξ = –10%\r\nFig. 12. Frequency-angle spectra: (a) ξ = 0; (b) ξ = –3%; (c) ξ = –10%.
    图 12  规则波纹圆柱频率-角度谱 (a) ξ = 0; (b) ξ = –3%; (c) ξ = –10%
    Fig. 12.  Frequency-angle spectra: (a) ξ = 0; (b) ξ = –3%; (c) ξ = –10%.

    本文提出了一种微弱形变的规则波纹表面结构, 可实现PMMA圆柱体反向散射低频共振频率的无源调控. 利用微扰法推导了无限长微弱形变规则波纹圆柱共振频率的近似解, 研究发现由于边界微弱形变Rayleigh波传播路径发生改变, 沿着规则波纹表面传播, 传播路径的改变引起传播相位的变化, 从而导致共振峰频率发生偏移. 讨论了波纹微扰系数、周期对规则波纹圆柱共振频率偏移的影响规律: 当归一化波纹微扰系数ξ > 0时, 规则波纹圆柱反向散射共振频率随着ξ的增大向低频偏移, 当ξ < 0时, 规则波纹圆柱反向散射共振频率随着ξ的减小而向高频偏移, 共振峰频率偏移量随着波纹微扰系数绝对值的增大而增大; 当波纹周期m < 10时, 波纹周期对共振频率偏移影响较小. 最后开展了规则波纹圆柱体声散射特性水池实验, 获取了不同波纹微扰系数规则波纹圆柱体的反向散射共振频率, 明显观察到了规则波纹圆柱共振峰频率偏移现象, 与理论预报结果吻合较好. PMMA圆柱低频亚音速Rayleigh波隧穿共振引起反向散射幅度增强和规则波纹结构引起低频共振频率偏移的特征, 使得规则波纹PMMA圆柱在标准散射体设计、利用“声条形码”进行身份识别方面具有潜在的应用价值.

    此外, 文中基于微扰法推导的用于计算无限长规则波纹圆柱声散射的解, 只适用于微弱形变, 即波纹微扰系数和波纹周期较小的情况, 此时忽略了由于法向变化引起的微扰, 只考虑了函数值引起的微扰, 是一种近似解. 若要对本问题或相关问题进行深入探讨, 获取更为精确的解, 可参考文献[20]中规则形腔中简正模式的微扰近似方法.

    对于正横方向入射无限长规则波纹圆柱体, 行列式Ddn, BdnAbpn分别为

    Ddn=|Dd11Dd12Dd13000Dd21Dd22Dd230000Dd32Dd33000Dd41Dd42Dd43Dd44Dd45Dd46Dd51Dd52Dd53Dd54Dd55Dd560Dd62Dd630Dd65Dd66|,
    Bdn=|A1Dd12Dd13000A2Dd22Dd230000Dd32Dd33000A4Dd42Dd43Dd44Dd45Dd46A5Dd52Dd53Dd54Dd55Dd560Dd62Dd630Dd65Dd66|,
    Abpn=|Dd11Dd12Dd13A100Dd21Dd22Dd23A2000Dd32Dd33000Dd41Dd42Dd43A4Dd45Dd46Dd51Dd52Dd53A5Dd55Dd560Dd62Dd630Dd65Dd66|,
    A1=[A1+ε(A12A1)Ann],A2=[A2+ε(A22A2)Ann],A4=(A12A1)Anp,A5=(A22A2)Anp,
    Dd11=d11+ε(d112d11)Ann,Dd12=d12+ε(d122d12)Ann,Dd13=d13+ε(d132d13)Ann,Dd21=d21+ε(d21d21)Ann,Dd22=d22+ε(d22d22)Ann,Dd23=d23+ε(d23d23)Ann,Dd32=d32+ε(d322d32)A1nn,Dd33=d33+ε(d332d33)A1nn,Dd41=(d112d11)Anp,Dd42=(d122d12)Anp,Dd43=(d132d13)Anp,Dd44=d11p+ε(d11p2d11p)App,Dd45=d12p+ε(d12p2d12p)App,Dd46=d13p+ε(d13p2d13p)App,Dd51=(d21d21)Anp,Dd52=(d22d22)Anp,Dd53=(d23d23)Anp,Dd54=d21p+ε(d21pd21p)App,Dd55=d22p+ε(d22pd22p)App,Dd55=d22p+ε(d22pd22p)App,Dd62=(d322d32)A1np,Dd63=(d332d33)A1np,Dd65=d32p+ε(d32p2d32p)A1pp,Dd66=d33p+ε(d33p2d33p)A1pp,
    A1=ρ1ρ21ρ1ω2k2sa2Jn(ka),   A2=kaρ1ω2Jn(ka),A1=ρ1ρ21ρ1ω2[2k2sa2Jn(ka)+kak2sa2Jn(ka)],A2=kaρ1ω2Jn(ka)k2a2ρω2Jn(ka),
    d11=ρ1ρ21ρ1ω2k2sa2H(1)n(ka),d12=(2n2k2sa2)Jn(kda)2kdaJn(kda),d13=2n[ksaJn(ksa)Jn(ksa)],d11=ρ1ρ21ρ1ω2[2k2sa2H(1)n(ka)+kak2sa2H(1)n(ka)],d12=2k2sa2Jn(kda)+kda(2n2k2sa22)Jn(kda)2k2da2Jn(kda),d13=2nk2sa2Jn(ksa),
    d21=kaρ1ω2H(1)n(ka),d22=kdaJn(kda),   d23=nJn(ksa),d21=kaρ1ω2H(1)n(ka)k2a2ρ1ω2H(1)n(ka),d22=kdaJn(kda)+k2da2Jn(kda),  d23=nksaJn(ksa),
    d32=2n[Jn(kda)kdaJn(kda)],d33=2ksaJn(ksa)+(k2sa22n2)Jn(ksa),d32=2nk2da2Jn(kda),d33=2k2sa2Jn(ksa)+2k2sa2Jn(ksa)+ksa[(ksa)22n2+2]Jn(ksa),
    A1nn=εnππ0f(φ)asin(nφ)sin(nφ)dφ,A1pp=εpππ0f(φ)asin(pφ)sin(pφ)dφ,App=εpππ0f(φ)acos(pφ)cos(pφ)dφ,

    其中, d11p, d12p, d13p, d21p, d22p, d23p, d32p, d33pd11p, d12p, d13p, d21p, d22p, d23p, d32p, d33p分别等于d11, d12, d13, d21, d22, d23, d32, d33d11, d12, d13, d21, d22, d23, d32, d33中将阶次n替换成p的值. JnJn分别代表贝塞尔函数的一阶和二阶导数.

    kdR=c1cdx, ksR=c1csx代入(6a)式, 化简得

    τrrn(x)+ptn(x)=μk2({ρ1ρ21ρ1ω2c21c2sH(1)n(x)Bn[(2n2x2c21c2s)Jn(c1cdx)2c1cd1xJn(c1cdx)]Cn+2n[c1cs1xJn(c1csx)1x2Jn(c1csx)]En+ρ1ρ21ρ1ω2c21c2sJn(x)}cos(nφ)+εpn{ρ1ρ21ρ1ω2c21c2sH(1)p(x)αbp[(2p2x2c21c2s)Jp(c1cdx)2c1cd1xJp(c1cdx)]αcp+2p[c1cs1xJp(c1csx)1x2Jp(c1csx)]αep}cos(pφ)),

    将(A1)式中xx0处按照ε的一阶展开得

    τrrn(x0)+ptn(x0)=μk2[εx1{ρ1ρ21ρ1ω2c21c2sH(1)n(x0)Bn+[c1cd(2n2x02c21c2s+2x02)Jn(c1cdx0)4n2x30Jn(c1cdx0)2c12cd21x0Jn(c1cdx0)]Cn+2n[c12cs21x0Jn(c1csx0)c1cs2x02Jn(c1csx0)+2x03Jn(c1csx0)]En+ρ1ρ21ρ1ω2c21c2sJn(x0)}cos(nφ)+ε(pn{ρ1ρ21ρ1ω2c21c2sH(1)p(x0)αbp[(2p2x20c21c2s)Jp(c1cdx0)2c1cd1x0Jp(c1cdx0)]αcp+2p[c1cs1x0Jp(c1csx0)1x20Jp(c1csx0)]αep})cos(pφ)],
    τrrnr+ptnr=μk3[(ρ1ρ21ρ1ω2c21c2sH(1)n(x)Bn+[c1cd(2n21x2c21c2s+21x2)Jn(c1cdx)2c21c2d1xJn(c1cdx)4n21x3Jn(c1cdx)]Cn+{2nc21c2s1xJn(c1csx)4n[c1cs1x2Jn(c1csx)1x3Jn(c1csx)]}En+ρ1ρ21ρ1ω2c21c2sJn(x))cos(nφ)+εpn(ρ1ρ21ρ1ω2c21c2sH(1)p(x)αbp+{c1cd(2p21x2c21c2s+21x2)Jp(c1cdx)2c21c2d1xJp(c1cdx)4p21x3Jp(c1cdx)}αcp+{2pc21c2s1xJp(c1csx)4p[c1cs1x2Jp(c1csx)1x3Jp(c1csx)]}αep)cos(pφ)],
    τrrn(x0)+ptn(x0)=εf(φ)[τrrnr(x0)+ptnr(x0)].

    将(A2)式和(A3)式代入(A4)式得到(13)式.

    [1]

    Ge L, Song Q H, Redding B, Eberspacher A, Wiersig J, Cao H 2013 Phys. Rev. A. 88 043801Google Scholar

    [2]

    Ge L, Song Q H, Redding B, Cao H 2013 H Phys. Rev. A. 87 023833.1Google Scholar

    [3]

    Fawcett J A 2016 IEEE J. Ocean. Eng. 41 682Google Scholar

    [4]

    Zhang L G, Sun N H, Marston P L 1992 J. Acoust. Soc. Am. 91 1862Google Scholar

    [5]

    Fawcett J A 2001 J. Acoust. Soc. Am. 109 1312Google Scholar

    [6]

    Guel-Tapia J A, Villa-Villa F, Mendoza-Suarez A, Perez-Aguilar H 2016 Archives of Acoustics 41 461Google Scholar

    [7]

    Hartmann B, Jarzynski B 1974 J. Acoust. Soc. Am. 56 1469Google Scholar

    [8]

    Hefner B T, Marston P L 2000 J. Acoust. Soc. Am. 107 1930Google Scholar

    [9]

    Satish A, Trivett D, Sabra K G 2020 J. Acoust. Soc. Am. 147 EL517Google Scholar

    [10]

    Srivastava P, Nichols B, Sabra K G 2017 J. Acoust. Soc. Am. 142 EL573Google Scholar

    [11]

    Doolittle R D, Überall H, Uginčius P 1968 J. Acoust. Soc. Am. 43 1Google Scholar

    [12]

    Dubertrand R, Bogomolny E, Djellali N, Lebental M, Schmit C 2008 Phys. Rev. A. 77 013804.1Google Scholar

    [13]

    Donald B, Gaunaurd G C 1983 J. Acoust. Soc. Am. 73 1448Google Scholar

    [14]

    龙云亮, 文希理, 谢处方 1994 数值计算与计算机应用 2 88Google Scholar

    You Y L, Wen X L, Xie C F 1994 Journal on Numerical Methods and Computer Applications 2 88Google Scholar

    [15]

    汤渭霖 1995 声学学报 06 456Google Scholar

    Tang W L 1995 ACTA. ACUSTICA. 06 456Google Scholar

    [16]

    Marston P L, Sun N H 1995 J. Acoust. Soc. Am. 97 777Google Scholar

    [17]

    周彦玲, 范军, 王斌 2019 物理学报 68 214301Google Scholar

    Zhou Y L, Fan J, Wang B 2019 Acta Phys. Sin. 68 214301Google Scholar

    [18]

    汤渭霖, 范军, 马忠诚 著 2018 水中目标声散射 (北京: 科学出版社) 第99—104页

    Tang W L, Fan J, Ma Z C 2018 Acoustic Scattering of Underwater Targets (Beijing: Science Press) pp99–104 (in Chinese)

    [19]

    彭茜蕤, 周彦玲, 范军 2018 声学技术 37 528

    Peng X R, Zhou Y L, Fan J 2018 Technical Acoustics 37 528

    [20]

    程建春 著 2019 声学原理 (北京: 科学出版社) 第615—619页

    Cheng J C 2019 Acoustical Principle (Beijing: Science Press) pp615–619 (in Chinese)

  • 图 1  无限长微弱形变规则波纹表面圆柱体

    Fig. 1.  Infinite regular-corrugated cylinder.

    图 2  基于微扰法无限长规则波纹圆柱形态函数幅频特性 (a)反向散射; (b)局部放大

    Fig. 2.  Form function of the infinite regular-corrugated cylinder based on the perturbation method: (a) Backscattering; (b) local enlargement.

    图 3  PMMA圆柱中Rayleigh波相速度频散曲线(黑色线), 曲线2πaf0/M(蓝色点划线)和曲线Sf/M(红色虚线)

    Fig. 3.  Dispersion curve of subsonic Rayleigh waves (black line), curve 2πaf0/M (blue dashed-dotted line) and curve Sf/M(red dotted line).

    图 4  亚音速Rayleigh波传播路径 (a)光滑圆柱; (b)规则波纹圆柱

    Fig. 4.  Ray diagram for subsonic Rayleigh waves propagating around: (a) Bare cylinder; (b) regular-corrugated cylinder.

    图 5  不同方法获取共振频率误差

    Fig. 5.  The relative error of resonance frequencies obtained by different methods.

    图 6  (a)频率-波纹微扰系数谱; (b)频率-波纹周期谱

    Fig. 6.  (a) Frequency-height spectra; (b) frequency-period spectra.

    图 7  不同波纹周期对应相位变化量

    Fig. 7.  Phase varying with corrugated period.

    图 8  实验模型

    Fig. 8.  Experimental objects.

    图 9  实验布放

    Fig. 9.  Diagram of experimental system setup.

    图 10  光滑圆柱和规则波纹表面圆柱反向散射形态函数实验结果

    Fig. 10.  Backscattering form function of regular-corrugated cylinders in the experiment.

    图 11  实验和(16)式获取得的共振频率对比 (a) ξ = –3%; (b) ξ = –10%; (c)相对误差

    Fig. 11.  Resonance frequencies obtained by experiment and Eq. (16): (a) ξ = –3%; (b) ξ = –10%; (c) relative error.

    图 12  规则波纹圆柱频率-角度谱 (a) ξ = 0; (b) ξ = –3%; (c) ξ = –10%

    Fig. 12.  Frequency-angle spectra: (a) ξ = 0; (b) ξ = –3%; (c) ξ = –10%.

    表 1  计算所用材料参数

    Table 1.  Material parameters used in the calculations.

    材料密度/(kg·m–3)纵波波速/(m·s–1)剪切波波速/(m·s–1)
    PMMA119026901340
    10001500
    下载: 导出CSV

    表 2  无限长PMMA圆柱和规则波纹圆柱共振频率

    Table 2.  Resonance frequency of PMMA bare and corrugated cylinder

    lf0/Hzf0/Hzfξ /Hzfξ/Hz
    16419642062596240
    29953995597049700
    313383133851304813050
    416790167901637016375
    下载: 导出CSV

    表 3  实验模型共振峰频率

    Table 3.  Resonance frequencies of experimental objects.

    lf0/Hzfξ1/Hzfξ1 /Hzfξ2/Hzfξ2/Hz
    166496748.768526981.46966
    286438772.690989075.19294
    31241012596.01299613031.013192
    下载: 导出CSV
  • [1]

    Ge L, Song Q H, Redding B, Eberspacher A, Wiersig J, Cao H 2013 Phys. Rev. A. 88 043801Google Scholar

    [2]

    Ge L, Song Q H, Redding B, Cao H 2013 H Phys. Rev. A. 87 023833.1Google Scholar

    [3]

    Fawcett J A 2016 IEEE J. Ocean. Eng. 41 682Google Scholar

    [4]

    Zhang L G, Sun N H, Marston P L 1992 J. Acoust. Soc. Am. 91 1862Google Scholar

    [5]

    Fawcett J A 2001 J. Acoust. Soc. Am. 109 1312Google Scholar

    [6]

    Guel-Tapia J A, Villa-Villa F, Mendoza-Suarez A, Perez-Aguilar H 2016 Archives of Acoustics 41 461Google Scholar

    [7]

    Hartmann B, Jarzynski B 1974 J. Acoust. Soc. Am. 56 1469Google Scholar

    [8]

    Hefner B T, Marston P L 2000 J. Acoust. Soc. Am. 107 1930Google Scholar

    [9]

    Satish A, Trivett D, Sabra K G 2020 J. Acoust. Soc. Am. 147 EL517Google Scholar

    [10]

    Srivastava P, Nichols B, Sabra K G 2017 J. Acoust. Soc. Am. 142 EL573Google Scholar

    [11]

    Doolittle R D, Überall H, Uginčius P 1968 J. Acoust. Soc. Am. 43 1Google Scholar

    [12]

    Dubertrand R, Bogomolny E, Djellali N, Lebental M, Schmit C 2008 Phys. Rev. A. 77 013804.1Google Scholar

    [13]

    Donald B, Gaunaurd G C 1983 J. Acoust. Soc. Am. 73 1448Google Scholar

    [14]

    龙云亮, 文希理, 谢处方 1994 数值计算与计算机应用 2 88Google Scholar

    You Y L, Wen X L, Xie C F 1994 Journal on Numerical Methods and Computer Applications 2 88Google Scholar

    [15]

    汤渭霖 1995 声学学报 06 456Google Scholar

    Tang W L 1995 ACTA. ACUSTICA. 06 456Google Scholar

    [16]

    Marston P L, Sun N H 1995 J. Acoust. Soc. Am. 97 777Google Scholar

    [17]

    周彦玲, 范军, 王斌 2019 物理学报 68 214301Google Scholar

    Zhou Y L, Fan J, Wang B 2019 Acta Phys. Sin. 68 214301Google Scholar

    [18]

    汤渭霖, 范军, 马忠诚 著 2018 水中目标声散射 (北京: 科学出版社) 第99—104页

    Tang W L, Fan J, Ma Z C 2018 Acoustic Scattering of Underwater Targets (Beijing: Science Press) pp99–104 (in Chinese)

    [19]

    彭茜蕤, 周彦玲, 范军 2018 声学技术 37 528

    Peng X R, Zhou Y L, Fan J 2018 Technical Acoustics 37 528

    [20]

    程建春 著 2019 声学原理 (北京: 科学出版社) 第615—619页

    Cheng J C 2019 Acoustical Principle (Beijing: Science Press) pp615–619 (in Chinese)

  • [1] 田志富, 吴迪, 胡涛. 圆柱曲面单光子量子雷达散射截面的理论研究. 物理学报, 2022, 71(3): 034204. doi: 10.7498/aps.71.20211295
    [2] 周彦玲, 范军, 王斌, 李兵. 水下环形凹槽圆柱体散射声场空间指向性调控. 物理学报, 2021, 70(17): 174301. doi: 10.7498/aps.70.20210111
    [3] 田志富, 吴迪, 胡涛. 圆柱曲面单光子量子雷达散射截面的理论研究. 物理学报, 2021, (): . doi: 10.7498/aps.70.20211295
    [4] 张兴坊, 刘凤收, 闫昕, 梁兰菊, 韦德全. 同心椭圆柱-纳米管结构的双重Fano共振研究. 物理学报, 2019, 68(6): 067301. doi: 10.7498/aps.68.20182249
    [5] 金国梁, 尹剑飞, 温激鸿, 温熙森. 基于等效参数反演的敷设声学覆盖层的水下圆柱壳体声散射研究. 物理学报, 2016, 65(1): 014305. doi: 10.7498/aps.65.014305
    [6] 潘安, 范军, 王斌, 陈志刚, 郑国垠. 双层周期加肋有限长圆柱壳声散射精细特征研究. 物理学报, 2014, 63(21): 214301. doi: 10.7498/aps.63.214301
    [7] 刘建晓, 张郡亮, 苏明敏. 基于时域有限差分法的各向异性铁氧体圆柱电磁散射分析. 物理学报, 2014, 63(13): 137501. doi: 10.7498/aps.63.137501
    [8] 潘安, 范军, 卓琳凯. 准周期加隔板有限长圆柱壳声散射. 物理学报, 2013, 62(2): 024301. doi: 10.7498/aps.62.024301
    [9] 刘敏, 侯志林, 傅秀军. 二维正方排列圆柱状亥姆赫兹共振腔阵列局域共振声带隙的研究. 物理学报, 2012, 61(10): 104302. doi: 10.7498/aps.61.104302
    [10] 潘安, 范军, 卓琳凯. 周期性加隔板有限长圆柱壳声散射. 物理学报, 2012, 61(21): 214301. doi: 10.7498/aps.61.214301
    [11] 王运华, 张彦敏, 郭立新. 两相邻有限长圆柱的复合电磁散射研究. 物理学报, 2011, 60(2): 021102. doi: 10.7498/aps.60.021102
    [12] 王泽锋, 胡永明, 罗洪, 孟洲, 倪明, 熊水东. 腔壁弹性对水下小型圆柱形亥姆霍兹共振器共振频率的影响. 物理学报, 2009, 58(4): 2507-2512. doi: 10.7498/aps.58.2507
    [13] 王泽锋, 胡永明, 孟 洲, 倪 明. 水下圆柱形Helmholtz共振器的声学特性分析. 物理学报, 2008, 57(11): 7022-7029. doi: 10.7498/aps.57.7022
    [14] 伍瑞新, 王耀俊. 各向同性媒质中球体声散射的共振模式缺损. 物理学报, 2003, 52(8): 1948-1953. doi: 10.7498/aps.52.1948
    [15] 周宇峰, 王耀俊, 马力, 高天赋. 流体饱和多孔圆柱体的声波散射. 物理学报, 2000, 49(3): 480-486. doi: 10.7498/aps.49.480
    [16] 董正超. 多层金属圆柱线中量子输运的表面和界面散射效应. 物理学报, 1999, 48(1): 127-133. doi: 10.7498/aps.48.127
    [17] 宫野, 温晓军, 张鹏云, 邓新绿. 圆柱模型下电子回旋共振微波等离子体离子输运过程的数值研究. 物理学报, 1997, 46(12): 2376-2383. doi: 10.7498/aps.46.2376
    [18] 叶芃生. 重子共振态的选择规则. 物理学报, 1965, 21(11): 1924-1926. doi: 10.7498/aps.21.1924
    [19] 竺迺宜. 低能π-N散射p波共振中交叉过程的贡献. 物理学报, 1964, 20(5): 411-417. doi: 10.7498/aps.20.411
    [20] 蒲富恪, 郑庆祺. 自旋波之间的散射对铁磁共振曲线的影响. 物理学报, 1962, 18(2): 81-90. doi: 10.7498/aps.18.81
计量
  • 文章访问数:  5531
  • PDF下载量:  61
出版历程
  • 收稿日期:  2020-09-15
  • 修回日期:  2020-10-28
  • 上网日期:  2021-02-26
  • 刊出日期:  2021-03-05

/

返回文章
返回