搜索

x

留言板

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

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

水下环形凹槽圆柱体散射声场空间指向性调控

周彦玲 范军 王斌 李兵

引用本文:
Citation:

水下环形凹槽圆柱体散射声场空间指向性调控

周彦玲, 范军, 王斌, 李兵

Manipulating spatial directivity of acoustic scattering from a submerged cylinder by means of annular grooves

Zhou Yan-Ling, Fan Jun, Wang Bin, Li Bing
PDF
HTML
导出引用
  • 本文提出一种具有深度梯度的环形凹槽结构, 可用于调控水中有限长刚性圆柱体散射声场空间指向性. 基于声学相位阵列理论分析了环形凹槽圆柱声散射空间指向性改变的机理, 研究表明: 凹槽深度方向相位延迟和凹槽间Bragg散射的相互作用使得平面声波垂直于圆柱方向入射其正横方向散射声波发生偏转. 采用有限元方法讨论了凹槽结构参数如占空比、梯度等对圆柱散射声场空间分布特征的影响规律. 多个不同深度梯度环形凹槽单元组合圆柱体散射声场数值计算和实验结果显示: 具有环形凹槽结构圆柱体正横方向散射声波均匀偏转到预定的空间范围内, 使得圆柱体声散射场空间指向性均衡化, 改变了圆柱整体的散射特征, 这为水下目标声隐身设计和声波定向传播提供了新的方法.
    An annular groove (AG) structure with depth gradient is proposed which can manipulate the spatial distribution of the acoustic scattering field for a finite rigid cylinder in water. An analytical analysis is given for better understanding the underlying mechanism of the abnormal scattered wave, which can be accomplished by using the phased array theory. When the plane acoustic wave is normally incident, the scattering acoustic wave in the transverse direction of the cylinder deflects, which is due to the interaction between the phase delay modulated by the AG structure with varying groove depths and the Bragg scattering of adjacent grooves. The finite element method is used to calculate the acoustic scattering field of a finite rigid cylinder with annular grooves and obtain the frequency and spatial distribution characteristics. How the structural parameters such as depth, gradient, and duty ratio of the annular grooves affect the acoustic scattering field is discussed in detail. The results show that the target strength in the transverse direction decreases linearly with duty ratio increasing while the target strength in the deflection direction of the acoustic wave increases with the duty ratio until δ = 30%, after which it remains almost constant. When the incident acoustic wave is fixed, the acoustic scattering wave of the AG cylinder can be deflected by designing the gradient appropriately, and the deflection direction is independent of the frequency. Numerical and experimental results for a cylinder with multiple annular-groove units show that the spatial directivity of the scattering field of the grooved cylinder changes, and the target strength is enhanced at six pre-designed deflection angles. Meanwhile, the deflected acoustic wave has a certain width and the interference among periodic structures of the AG units exists, which makes the spatial directivity of the scattering field of the cylinder equalize and changes the scattering characteristics of the cylinder, thereby providing a theoretical basis for designing three-dimensional underwater objects each with an acoustic stealth.
      通信作者: 范军, fanjun@sjtu.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 11774229)资助的课题
      Corresponding author: Fan Jun, fanjun@sjtu.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant No. 11774229)
    [1]

    Li Y, Liang B, Gu Z M, Zou X Y, Cheng J C 2013 Sci. Rep. 3 2546Google Scholar

    [2]

    Zhao J J, Li B W, Chen Z N, Qiu C W 2013 Sci. Rep. 3 2537Google Scholar

    [3]

    Zhao J J, Li B W, Chen Z N, Qiu C W 2013 Appl. Phys. Lett. 103 151604Google Scholar

    [4]

    Li Y, Xue J, Li R Q, Liang B, Zou X Y, Yin L L, Cheng J C 2014 Phys. Rev. Appl. 2 064002Google Scholar

    [5]

    Mei J, Wu Y 2014 New J. Phys. 16 123007Google Scholar

    [6]

    Fan X D, Zhu Y F, Liang B, Yang J, Cheng J C 2016 Appl. Phys. Lett. 109 243501Google Scholar

    [7]

    Zhu Y F, Fan X D, Liang B, Yang J, Yang J, Yin L L, Cheng J C 2016 AIP Adv. 6 121702Google Scholar

    [8]

    Ye Y T, Ke M Z, Li Y X, Wang T, Liu Z Y 2013 J. Appl. Phys. 114 154504Google Scholar

    [9]

    Chen J, Sun Z Q, Fan Z 2019 Appl. Phys. Lett. 114 254102Google Scholar

    [10]

    Li J F, Wang W Q, Xie Y B, Popa B I, Cummer S A 2016 Appl. Phys. Lett. 109 091908Google Scholar

    [11]

    Long H Y, Cheng Y, Liu X J 2017 Appl. Phys. Lett. 111 143502Google Scholar

    [12]

    Long H Y, Gao S X, Cheng Y, Liu X J 2018 Appl. Phys. Lett. 112 033507Google Scholar

    [13]

    Li Y and M. Assouar B 2016 Appl. Phys. Lett. 108 063502Google Scholar

    [14]

    Melde K, Mark A G, Qiu T, Fischer P 2016 Nature 537 518Google Scholar

    [15]

    Tian Y, Wei Q, Cheng Y, Liu X J 2017 Appl. Phys. Lett. 110 191901Google Scholar

    [16]

    Xie Y B, Shen C, Wang W Q, Li J F, Suo D J, Popa B Jing Y, Cummer S A 2016 Sci. Rep. 6 35437Google Scholar

    [17]

    Zhu Y F, Assouar B 2019 Phys. Rev. Mater. 3 045201Google Scholar

    [18]

    Christensen J, Fernandez D A I, De L P F, Martin M L, Garcia V F J 2007 Nat. Phys. 3 851Google Scholar

    [19]

    Zhu J, Chen Y Y, Zhu X F, Garcia V F J, Yin X B, Zhang W L, Zhang X 2013 Sci. Rep. 3 1728Google Scholar

    [20]

    Jia H, Lu M H, Ni X, Bao M, Li X D 2014 J. Appl. Phys. 116 124504Google Scholar

    [21]

    Jia H, Lu M H, Wang Q C, Bao M, Li X D 2013 Appl. Phys. Lett. 103 103505Google Scholar

    [22]

    Zhu Y F, Zou X Y, Li R Q, Jiang X, Tu J, Liang B, Cheng J C 2015 Sci. Rep. 5 10966Google Scholar

    [23]

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

    [24]

    Wu X X, Xia X X, Tian J X, Liu Z Y, Wen W J 2016 Appl. Phys. Lett. 108 163502Google Scholar

    [25]

    Liu J F, Declercq N F 2016 Appl. Phys. Lett. 109 261603Google Scholar

    [26]

    Lee H K, Jung M, Kim M, Shin R, Kang S, Ohm W, Kim Y T 2018 J. Acoust. Soc. Am. 143 1534Google Scholar

    [27]

    程建春 2019 声学原理 (北京: 科学出版社) 第89−90页

    Cheng J C 2019 Acoustical Principle (Beijing: Science Press) pp89−90 (in Chinese)

    [28]

    朱一凡, 梁彬, 程建春 2018 应用声学 37 53Google Scholar

    Zhu Y F, Liang B, Cheng J C 2018 J. Appl. Acoust. 37 53Google Scholar

    [29]

    汤渭霖, 陈德智 1988 声学学报 1 29

    Tang W L, Chen D Z 1988 Acta. Acustica. 1 29

    [30]

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

    Tang W L, Fan J, Ma Z C 2018 Acoustic Scattering of Underwater Targets (Beijing: Science Press) pp55−66 (in Chinese)

  • 图 1  (a)具有环形凹槽结构圆柱示意图; (b)蓝色虚线框局部放大; (c)圆柱散射声场收发分置示意图

    Fig. 1.  (a) Schematic illustration of an annular groove cylinder; (b) details of the blue dotted box; (c) bistatic diagram of simulation

    图 2  频率-角度谱 (a)刚性圆柱; (b)环形凹槽结构圆柱

    Fig. 2.  Frequency-angle spectra of target strength of the finite: (a) Rigid cylinder; (b) annular groove cylinder

    图 3  凹槽圆柱归一化指向性函数频率-角度谱

    Fig. 3.  Frequency-angle spectra of the normalized directional factors for the annular groove cylinder by Eq. (5)

    图 4  频率f = 80 kHz占空比δ = 0和δ = 83.3%凹槽圆柱目标强度空间指向性

    Fig. 4.  Spatial directivity of target strength of the annular groove cylinder with δ = 0 and 83.3% at f = 80 kHz.

    图 5  f = 80 kHz, 不同占空比凹槽圆柱在45°和90°方位目标强度

    Fig. 5.  Target strength of the annular groove cylinder with varying δ in the 45° and 90° direction at f = 80 kHz.

    图 6  f = 80 kHz, 不同梯度环形凹槽结构圆柱目标强度空间指向性

    Fig. 6.  Spatial directivity of target strength of the annular groove cylinder with different g at f = 80 kHz.

    图 7  6个环形凹槽单元组合圆柱目标强度指向性, f = 80 kHz

    Fig. 7.  Spatial directivity of target strength of the cylinder with six annular groove units at f = 80 kHz

    图 8  实验模型

    Fig. 8.  Experimental models.

    图 9  实验装置布放图

    Fig. 9.  Diagram of experimental system setup.

    图 10  时间-角度谱 (a)刚性圆柱; (b)环形凹槽单元组合圆柱

    Fig. 10.  Time-angle spectra: (a) Rigid cylinder; (b) annular groove cylinder.

    图 11  凹槽单元组合圆柱目标强度频率-角度谱 (a)数值计算结果; (b)实验结果; (c) 6个凹槽单元圆柱二维图

    Fig. 11.  Frequency-angle spectra of target strength for annular groove cylinder: (a) Numerical result; (b) experimental result; (c) 2-D geometry of the annular groove cylinder.

    图 12  频率f = 80 kHz凹槽圆柱目标强度指向性

    Fig. 12.  The normalized directivity of target strength for cylinder with six annular groove units at f = 80 kHz

    图 13  频率f = 80 kHz圆柱和凹槽圆柱目标强度实验结果对比

    Fig. 13.  The comparison of target strength in experiment between cylinder and annular groove cylinder at f = 80 kHz

  • [1]

    Li Y, Liang B, Gu Z M, Zou X Y, Cheng J C 2013 Sci. Rep. 3 2546Google Scholar

    [2]

    Zhao J J, Li B W, Chen Z N, Qiu C W 2013 Sci. Rep. 3 2537Google Scholar

    [3]

    Zhao J J, Li B W, Chen Z N, Qiu C W 2013 Appl. Phys. Lett. 103 151604Google Scholar

    [4]

    Li Y, Xue J, Li R Q, Liang B, Zou X Y, Yin L L, Cheng J C 2014 Phys. Rev. Appl. 2 064002Google Scholar

    [5]

    Mei J, Wu Y 2014 New J. Phys. 16 123007Google Scholar

    [6]

    Fan X D, Zhu Y F, Liang B, Yang J, Cheng J C 2016 Appl. Phys. Lett. 109 243501Google Scholar

    [7]

    Zhu Y F, Fan X D, Liang B, Yang J, Yang J, Yin L L, Cheng J C 2016 AIP Adv. 6 121702Google Scholar

    [8]

    Ye Y T, Ke M Z, Li Y X, Wang T, Liu Z Y 2013 J. Appl. Phys. 114 154504Google Scholar

    [9]

    Chen J, Sun Z Q, Fan Z 2019 Appl. Phys. Lett. 114 254102Google Scholar

    [10]

    Li J F, Wang W Q, Xie Y B, Popa B I, Cummer S A 2016 Appl. Phys. Lett. 109 091908Google Scholar

    [11]

    Long H Y, Cheng Y, Liu X J 2017 Appl. Phys. Lett. 111 143502Google Scholar

    [12]

    Long H Y, Gao S X, Cheng Y, Liu X J 2018 Appl. Phys. Lett. 112 033507Google Scholar

    [13]

    Li Y and M. Assouar B 2016 Appl. Phys. Lett. 108 063502Google Scholar

    [14]

    Melde K, Mark A G, Qiu T, Fischer P 2016 Nature 537 518Google Scholar

    [15]

    Tian Y, Wei Q, Cheng Y, Liu X J 2017 Appl. Phys. Lett. 110 191901Google Scholar

    [16]

    Xie Y B, Shen C, Wang W Q, Li J F, Suo D J, Popa B Jing Y, Cummer S A 2016 Sci. Rep. 6 35437Google Scholar

    [17]

    Zhu Y F, Assouar B 2019 Phys. Rev. Mater. 3 045201Google Scholar

    [18]

    Christensen J, Fernandez D A I, De L P F, Martin M L, Garcia V F J 2007 Nat. Phys. 3 851Google Scholar

    [19]

    Zhu J, Chen Y Y, Zhu X F, Garcia V F J, Yin X B, Zhang W L, Zhang X 2013 Sci. Rep. 3 1728Google Scholar

    [20]

    Jia H, Lu M H, Ni X, Bao M, Li X D 2014 J. Appl. Phys. 116 124504Google Scholar

    [21]

    Jia H, Lu M H, Wang Q C, Bao M, Li X D 2013 Appl. Phys. Lett. 103 103505Google Scholar

    [22]

    Zhu Y F, Zou X Y, Li R Q, Jiang X, Tu J, Liang B, Cheng J C 2015 Sci. Rep. 5 10966Google Scholar

    [23]

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

    [24]

    Wu X X, Xia X X, Tian J X, Liu Z Y, Wen W J 2016 Appl. Phys. Lett. 108 163502Google Scholar

    [25]

    Liu J F, Declercq N F 2016 Appl. Phys. Lett. 109 261603Google Scholar

    [26]

    Lee H K, Jung M, Kim M, Shin R, Kang S, Ohm W, Kim Y T 2018 J. Acoust. Soc. Am. 143 1534Google Scholar

    [27]

    程建春 2019 声学原理 (北京: 科学出版社) 第89−90页

    Cheng J C 2019 Acoustical Principle (Beijing: Science Press) pp89−90 (in Chinese)

    [28]

    朱一凡, 梁彬, 程建春 2018 应用声学 37 53Google Scholar

    Zhu Y F, Liang B, Cheng J C 2018 J. Appl. Acoust. 37 53Google Scholar

    [29]

    汤渭霖, 陈德智 1988 声学学报 1 29

    Tang W L, Chen D Z 1988 Acta. Acustica. 1 29

    [30]

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

    Tang W L, Fan J, Ma Z C 2018 Acoustic Scattering of Underwater Targets (Beijing: Science Press) pp55−66 (in Chinese)

  • [1] 周彦玲, 王斌, 范军. 水中微小波纹圆柱体声散射低频共振调控. 物理学报, 2021, 70(5): 054301. doi: 10.7498/aps.70.20201535
    [2] 马瑞轩, 王益民, 张树海, 武从海, 王勋年. 旋涡声散射特性的尺度效应数值研究. 物理学报, 2021, 70(10): 104301. doi: 10.7498/aps.70.20202206
    [3] 王益民, 马瑞轩, 武从海, 罗勇, 张树海. 旋涡声散射的空间尺度特性数值研究. 物理学报, 2021, 70(19): 194302. doi: 10.7498/aps.70.20202232
    [4] 冯康艺, 王成会. 超声场中空化泡对弹性粒子微流的影响. 物理学报, 2019, 68(24): 244301. doi: 10.7498/aps.68.20191253
    [5] 周彦玲, 范军, 王斌. 塑料类高分子聚合物材料水中目标声学参数反演. 物理学报, 2019, 68(21): 214301. doi: 10.7498/aps.68.20190991
    [6] 范雨喆, 李海森, 徐超, 陈宝伟, 杜伟东. 基于声散射的水下气泡群空间关联性研究. 物理学报, 2017, 66(1): 014305. doi: 10.7498/aps.66.014305
    [7] 秦修培, 耿德路, 洪振宇, 魏炳波. 超声悬浮过程中圆柱体的旋转运动机理研究. 物理学报, 2017, 66(12): 124301. doi: 10.7498/aps.66.124301
    [8] 金国梁, 尹剑飞, 温激鸿, 温熙森. 基于等效参数反演的敷设声学覆盖层的水下圆柱壳体声散射研究. 物理学报, 2016, 65(1): 014305. doi: 10.7498/aps.65.014305
    [9] 潘安, 范军, 王斌, 陈志刚, 郑国垠. 双层周期加肋有限长圆柱壳声散射精细特征研究. 物理学报, 2014, 63(21): 214301. doi: 10.7498/aps.63.214301
    [10] 何春涛, 王聪, 何乾坤, 仇洋. 圆柱体低速入水空泡试验研究. 物理学报, 2012, 61(13): 134701. doi: 10.7498/aps.61.134701
    [11] 仝焕平, 章林溪. 受限于圆柱体内半刚性高分子链的螺旋结构转变. 物理学报, 2012, 61(5): 058701. doi: 10.7498/aps.61.058701
    [12] 吴海军, 蒋伟康, 鲁文波. 三维声学多层快速多极子边界元及其应用. 物理学报, 2012, 61(5): 054301. doi: 10.7498/aps.61.054301
    [13] 张超英, 李华兵, 谭惠丽, 刘慕仁, 孔令江. 椭圆柱体在牛顿流体中运动的格子Boltzmann方法模拟. 物理学报, 2005, 54(5): 1982-1987. doi: 10.7498/aps.54.1982
    [14] 董全林, 王 坤, 张春熹, 刘 彬. 圆柱体相对转动动力学方程的积分解. 物理学报, 2004, 53(2): 337-342. doi: 10.7498/aps.53.337
    [15] 尤云祥, 缪国平. 阻抗障碍物声散射的反问题. 物理学报, 2002, 51(2): 270-278. doi: 10.7498/aps.51.270
    [16] 尤云祥, 缪国平. 三维可穿透目标远场声波反演的一种指示器样本方法. 物理学报, 2002, 51(9): 2038-2051. doi: 10.7498/aps.51.2038
    [17] 赵建林, 杨德兴. 圆柱体的空间圆锥光反射、折射与衍射. 物理学报, 2002, 51(9): 1972-1977. doi: 10.7498/aps.51.1972
    [18] 尤云祥, 缪国平, 刘应中. 用近场声学测量信息可视化多个三维障碍物的一种快速算法. 物理学报, 2001, 50(6): 1103-1109. doi: 10.7498/aps.50.1103
    [19] 周宇峰, 王耀俊, 马力, 高天赋. 流体饱和多孔圆柱体的声波散射. 物理学报, 2000, 49(3): 480-486. doi: 10.7498/aps.49.480
    [20] 胡宁. 圆柱体后及轴对称体后激流平均速度及温度之分布. 物理学报, 1944, 5(1): 1-29. doi: 10.7498/aps.5.1
计量
  • 文章访问数:  697
  • PDF下载量:  25
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-01-17
  • 修回日期:  2021-02-16
  • 上网日期:  2021-08-19
  • 刊出日期:  2021-09-05

水下环形凹槽圆柱体散射声场空间指向性调控

  • 上海交通大学, 海洋工程国家重点实验室, 高新船舶与深海开发装备协同创新中心, 上海 200240
  • 通信作者: 范军, fanjun@sjtu.edu.cn
    基金项目: 国家自然科学基金(批准号: 11774229)资助的课题

摘要: 本文提出一种具有深度梯度的环形凹槽结构, 可用于调控水中有限长刚性圆柱体散射声场空间指向性. 基于声学相位阵列理论分析了环形凹槽圆柱声散射空间指向性改变的机理, 研究表明: 凹槽深度方向相位延迟和凹槽间Bragg散射的相互作用使得平面声波垂直于圆柱方向入射其正横方向散射声波发生偏转. 采用有限元方法讨论了凹槽结构参数如占空比、梯度等对圆柱散射声场空间分布特征的影响规律. 多个不同深度梯度环形凹槽单元组合圆柱体散射声场数值计算和实验结果显示: 具有环形凹槽结构圆柱体正横方向散射声波均匀偏转到预定的空间范围内, 使得圆柱体声散射场空间指向性均衡化, 改变了圆柱整体的散射特征, 这为水下目标声隐身设计和声波定向传播提供了新的方法.

English Abstract

    • 声学超表面打破传统超材料的限制, 在很小的尺度上实现调控, 具有紧凑、成本低、损耗小以及便于集成等优势, 具有广泛的应用前景, 对于声学领域的基础研究和新型声学器件的拓展具有重要的指导意义. 由于声学超表面具有界面相位突变的特点, 可以根据结构的几何尺寸任意调控0—2π的相位分布, 实现声波波阵面的调控[1], 声学超表面的出现为自由调控声波提供了一种新的途径. 近几年, 出现了很多声学超表面理论设计和实验验证, 实现了如声波特殊折射/反射[2-5]、聚焦[6-9]、完美吸收[10-13]、声波全息成像[14-17]等许多新现象. 梯度声学人工结构作为一种特殊的声学超表面, 在调控声波传播方向和能量的空间分布具有重要作用, 研究主要集中于空气声学. Christensen等[18]提出了利用凹槽结构产生的声表面波耦合实现声场调控, 自此凹槽结构在声学超表面中被广泛研究. Zhu等[19]通过在刚性平板上设计具有深度梯度的凹槽阵列调控声波能量的空间分布. Jia等[20]在刚性表面上刻有宽度梯度凹槽阵列以减慢声表面波的传播. 文献[21, 22]通过数值仿真和实验验证了刚性表面刻有圆孔半径梯度结构和圆孔深度梯度结构实现了声波的聚集和成像等. 在水声领域, 关于声学超表面的研究较少. Srivastava等[23]在平面上刻有周期性凹槽结构并利用其Bragg散射特征实现了对水下目标无源声学编码. Wu等[24]通过在黄铜圆板上刻有不同深度同心凹槽结构实现了水下超声波的聚焦. Liu等[25]通过实验验证了水下周期性凹槽表面的Scholte-Stoneley波不是频散的, 为利用凹槽结构进行宽带调控声波提供了理论依据. Lee等[26]通过在刚性平面上雕刻各类凹槽结构利用声波的热损耗实现声波的吸收. 以上这些基于凹槽结构的研究利用的是刚性散射, 要实现微变形表面的相位共振需要声波频率较高(兆赫兹级), 而水下高频声波随距离衰减较快, 因此这类声学结构表面设计在水声领域难以应用.

      常见的水下航行器外形近似为圆柱, 因此本文提出了一种具有深度梯度的环形凹槽结构可用于调控水中圆柱体散射声场空间指向性, 扩展了广义斯涅耳定律[27,28]的应用, 其最初设计用于调控二维平面界面超表面的异常反射、透射、吸声等. 本文针对三维曲面目标散射声场的调控展开了研究. 第2节基于有限元方法(以COMSOL multiphysics商业软件为计算平台)计算了具有环形凹槽结构圆柱的频域空间特性, 并分析了调控环形凹槽圆柱散射声场空间指向性的物理机理. 第3节讨论了环形凹槽结构参数如占空比、深度梯度对圆柱正横方向散射声波偏转的影响规律, 并通过精心设计6个不同梯度环形凹槽单元, 成功地将圆柱体正横方向散射声波偏转到预先设计好的空间范围内, 改变了圆柱体的整体声散射特征, 在水声中具有重要的应用价值. 第4节进行了6个环形凹槽单元组合圆柱收发分置声散射水池实验, 获取了其回波特性, 并得到了与理论预估相一致的结果.

    • 为了实现水下圆柱体散射声场空间指向性调控, 根据广义斯涅耳定律, 在圆柱表面引入一个空间变化的相位突变, 反射角${\theta _{\rm{r}}}$为:

      ${\theta _{\rm{r}}} = \arccos \left( { - \cos {\theta _{\rm{i}}} - \frac{\lambda }{{2{\rm{\pi }}}}\frac{{{\rm{d}}\varphi \left( x \right)}}{{{\rm{d}}x}}} \right),$

      其中${\theta _{\rm{r}}}$${\theta _{\rm{i}}}$分别是反射角和入射角, λ是波长, $\varphi \left( x \right)$为表面相位空间分布. 图1(a)给出了环形凹槽圆柱示意图, 圆柱表面刻有不同深度环形凹槽结构, 当入射波进入深度为$h\left( x \right)$的凹槽后被封闭端反射, 反射波和入射波在凹槽开口端的相移为$\varphi \left( x \right) = {{4{\text{π}}h\left( x \right)} / \lambda }$. 通过在不同空间位置选择适当的凹槽深度, 可以得到反射波相移的空间变化. 由于多基地声纳探测水下目标时, 正横方向散射最强, 因此本文只考虑平面声波垂直于目标正横方向入射的情况, 即${\theta _{\rm{i}}} = {90^ \circ }$, (1)式可以简化为:

      图  1  (a)具有环形凹槽结构圆柱示意图; (b)蓝色虚线框局部放大; (c)圆柱散射声场收发分置示意图

      Figure 1.  (a) Schematic illustration of an annular groove cylinder; (b) details of the blue dotted box; (c) bistatic diagram of simulation

      ${\theta _r} = \arccos \left( { - 2g\left( x \right)} \right),$

      其中$g(x) = {\rm{d}}h(x)/{\rm{d}}x$为环形凹槽阵列的深度梯度. 因此, 当入射声波一定时, 通过设计适当地凹槽深度梯度可以改变声波传播方向.

      由于环形凹槽圆柱体周向均匀, 散射声场与水平方位角无关. 为了减少计算量, 基于有限元方法采用COMSOL multiphysics软件中二维轴对称模型压力声学物理场对平面声波正横入射环形凹槽圆柱散射声场进行数值仿真. 环形凹槽圆柱体半径a = 2.5 cm, 长度L = 6.1 cm, 两个相邻凹槽中心的距离d = 0.6 cm, 凹槽宽度w = 0.5 cm, 凹槽深度为hj, j为凹槽序列号, ${h_j} = $$ \left\{ \begin{aligned}& {h_{\max }} - \left( {N - j} \right)gd, \;\;\;g > 0, \\ &{h_{\max }} + \left( {j - 1} \right)gd, \;\;\;\;g < 0, \end{aligned} \right.$(j = 1, 2, 3 ··· N), 环形凹槽最深的深度hmax = 1.95 cm, 深度梯度g = –0.35, 一个凹槽单元有N = 10个凹槽.

      为了分析环形凹槽圆柱散射声场频率和空间分布特征, 分别计算刚性圆柱和环形凹槽圆柱目标强度的频率-角度谱如图2所示, 计算频率为5—100 kHz, 频率增量为1 kHz. 其中横坐标为接收点角度θ, 即观测点与z轴正向的夹角, 纵坐标为频率, 单位为kHz, 伪彩色图像的亮度表示目标强度, 单位为dB.

      图  2  频率-角度谱 (a)刚性圆柱; (b)环形凹槽结构圆柱

      Figure 2.  Frequency-angle spectra of target strength of the finite: (a) Rigid cylinder; (b) annular groove cylinder

      对于光滑表面圆柱, 当入射声波垂直于圆柱方向入射时, 正横方向散射目标强度最强, 具有很强的方向性, 但由于有限长圆柱两端边界的限制, 图2(a)中以90°为中心两边出现明暗相间的“八”字形条纹, 这是圆柱侧面两个亮点之间的几何干涉效应[29,30]. 而当圆柱表面具有环形凹槽结构时, 散射声波方向发生偏转, 如图2(b)所示, 目标强度最强的位置由正横偏转到45°方向, 且声波偏转方向与频率无关, 但具有一定的宽度, 随着频率升高, 宽度变窄. 为了更好地理解深度梯度环形凹槽圆柱正横散射声波偏转的物理机理, 引入声学相位阵列理论. 如图1(b)中红色箭头所示, 入射声波沿着凹槽深度方向传播在封闭端反射然后到达观测点, 声波的传播可以简单地看作是一系列具有相位延迟的阵元组成的线列阵, 远场散射声压可以近似表示为:

      ${p_{\rm{s}}}\left( {\rho,\theta,\omega } \right) \approx \sum\limits_{j = 1}^N {{A_j}} {{\rm{e}}^{{\rm{i}}\left[ {\omega t + {k_0}\left( {{\rho _j} + 2{h_j}} \right)} \right]}},$

      其中${\rho _j}$为观测点和第j个凹槽中心的距离, ${\rho _j} = {\rho _0} + $$ L\cos \theta /2 + \left( {j - 1} \right)d\cos \theta$, 线列阵总长度$L = Nd$, ${\rho _0}$为观测点到阵中心的距离. 入射波是单位振幅的平面波且每个凹槽宽度相同, 所以每个凹槽幅度因子相同, 即${A_j} \equiv A$, 则(3)式化简为:

      $\begin{split} {p_{\rm{s}}}\left( {\rho,\theta,\omega } \right) \approx \;&A{{\rm{e}}^{{\rm{i}}{k_0}{\rho _0}}}{{\rm{e}}^{{\rm{i}}{k_0}d\left( {\cos \theta + 2g} \right)\left( {N - 1} \right)/2}}\\ &\times\frac{{\sin \left[ {{k_0}d\left( {\cos \theta + 2g} \right)N/2} \right]}}{{\sin \left[ {{k_0}d\left( {\cos \theta + 2g} \right)/2} \right]}}. \end{split}$

      根据(4)式获取散射声压幅值最大时的角度${\theta _0} = $$ \arccos ( - 2 g)$, 散射声压为$p\left( {\rho, {\theta _0}, \omega } \right) \approx NA{{\rm{e}}^{{\rm{i}}{k_0}{\rho _0}}}$, 则(4)式可以表示为$p\left( {\rho, \theta, \omega } \right) = p( \rho, {\theta _0}, $$ \omega )D\left( \theta \right)$, 归一化指向性函数$D\left( \theta \right)$

      $\begin{split} &{}\\D\left( \theta \right) =\;& \frac{1}{N}{{\rm{e}}^{{\rm{i}}{k_0}d\left( {\cos \theta + 2g} \right)\left( {N - 1} \right)/2}}\\ &\times\frac{{\sin \left[ {N{k_0}d\left( {\cos \theta + 2g} \right)/2} \right]}}{{\sin \left[ {{k_0}d\left( {\cos \theta + 2g} \right)/2} \right]}}. \end{split}$

      根据(5)式得到环形凹槽圆柱体归一化指向性函数的频率-角度谱(见图3), 在θ = 45°方向散射声波最强, 且偏转声波具有一定宽度, 主波束宽度与频率相关, 频率越高, 宽度越窄, 与数值计算结果吻合较好.

      图  3  凹槽圆柱归一化指向性函数频率-角度谱

      Figure 3.  Frequency-angle spectra of the normalized directional factors for the annular groove cylinder by Eq. (5)

    • 由于环形凹槽深度方向相位延迟和凹槽间Bragg散射的相互作用使得平面声波垂直入射圆柱正横散射声波发生偏转. 详细讨论环形凹槽结构参数, 如占空比、梯度对环形凹槽圆柱散射声场偏转效果的影响.

    • 定义环形凹槽结构占空比为$\delta ={w / d} \times $100%, d = 0.6 cm, 环形凹槽结构梯度g = –0.35. 基于有限元方法计算占空比δ = 0和83.3%的环形凹槽圆柱在频率f = 80 kHz目标强度的空间指向性如图4所示, 占空比δ = 0(即刚性圆柱, 图4中黑色线)在θ = 90°方向目标强度最强, 当圆柱表面具有凹槽结构δ = 83.3%时, θ = 45°方向目标强度最强, θ = 90°方向相较于圆柱目标强度约减小了15 dB. 图5显示了环形凹槽圆柱在θ = 45°(黑色线)和θ = 90°(红色虚线)方向目标强度随占空比变化规律, 显然随着占空比增大, 环形凹槽圆柱反向散射(θ = 90°)目标强度单调下降; 当δ < 30%时, 随着δ增大, 散射声波偏转方向(即θ = 45°)目标强度单调增加, 当δ > 30%, 其目标强度基本不变. 为了不改变圆柱整体外部结构而使圆柱散射声波偏转且正横方向目标强度较低, 文中数值仿真均选定凹槽宽度w = 0.5 cm, 此时环形凹槽圆柱正横目标强度相较于光滑圆柱降低了15 dB左右.

      图  4  频率f = 80 kHz占空比δ = 0和δ = 83.3%凹槽圆柱目标强度空间指向性

      Figure 4.  Spatial directivity of target strength of the annular groove cylinder with δ = 0 and 83.3% at f = 80 kHz.

      图  5  f = 80 kHz, 不同占空比凹槽圆柱在45°和90°方位目标强度

      Figure 5.  Target strength of the annular groove cylinder with varying δ in the 45° and 90° direction at f = 80 kHz.

    • 环形凹槽深度梯度分别为g = –0.25, g = –0.13和g = 0.25, g = 0.13的凹槽圆柱在频率f = 80 kHz的目标强度空间指向性如图6所示, 凹槽圆柱目标强度最大值分别出现在θ = 60°, 75°, 105°和120°方向, 与公式理论预报的声波偏转方向一致. 因此, 通过设计不同环形凹槽的深度梯度可以使得环形凹槽圆柱正横散射目标强度偏转到任意角度, 从而实现对有限长刚性圆柱散射声场空间分布的调控. 利用此特征, 设计不同深度梯度的环形凹槽单元组合使得圆柱正横方向目标强度在一定空间范围内均衡化, 从而改变圆柱整体的回波特征. 环形凹槽圆柱结构参数为: 半径a = 3.5 cm, L = 38 cm, 每个凹槽单元中凹槽最深的深度相同hmax = 1.95 cm. 6个环形凹槽单元梯度依次为g1 = –0.13, g2 = 0.25, g3 = –0.35, g4 = 0.35, g5 = –0.25, g6 = 0.13. 由公式可知, 6个环形凹槽单元对应声波偏转方向依次为θ1 = 75°, θ2 = 120°, θ3 = 45°, θ4 = 135°, θ5 = 60°和θ6 = 105°.

      图  6  f = 80 kHz, 不同梯度环形凹槽结构圆柱目标强度空间指向性

      Figure 6.  Spatial directivity of target strength of the annular groove cylinder with different g at f = 80 kHz.

      获取了频率f = 80 kHz 6个环形凹槽单元组合圆柱不同接收角度的目标强度(图7中红色虚线), 对比相同尺度光滑圆柱(图7中黑色线)目标强度的指向性可以发现: 环形凹槽圆柱在正横方向(即90°)较光滑圆柱目标强度约降低了15 dB, 正横方向散射声波均衡分布在45°—135°范围内. 利用此结构使得圆柱正横方向较为集中的散射能量被分散到其他区域内, 使得散射声场均衡化, 改变了圆柱整体的散射特征.

      图  7  6个环形凹槽单元组合圆柱目标强度指向性, f = 80 kHz

      Figure 7.  Spatial directivity of target strength of the cylinder with six annular groove units at f = 80 kHz

    • 开展了6个环形凹槽单元组合结构圆柱收发分置声散射特性水池实验, 此次实验是在上海交通大学水声工程所水池完成. 两个实验模型(圆柱和环形凹槽圆柱)均为实心不锈钢材料, 近似为刚性, 几何结构参数与3.2节中6个环形凹槽单元组合圆柱相同. 测试模型水平置于空间尺寸为5 m × 5 m × 5 m的充水水池中, 发射换能器采用刚性连接吊放位置距离目标2.15 m, 水听器为B&K8103, 固定于悬臂梁上距离目标1.05 m, 均匀旋转悬臂梁从而获取收发分置条件下测试模型正横方向入射, 不同接收方向的散射声场. 发射换能器、水听器和目标三者处于同一深度, 距离水面2.4 m, 实验模型及具体布放如图8图9所示. 为排除市电干扰, 测量过程中所有设备全程处于UPS直流供电状态.

      图  8  实验模型

      Figure 8.  Experimental models.

      图  9  实验装置布放图

      Figure 9.  Diagram of experimental system setup.

    • 发射信号为线性调频信号, 中心频率f0 = 80 kHz, 带宽B = 40 kHz, 脉宽τ = 0.5 ms, 发射周期T = 500 ms. 首先对实验获取的时域信号进行脉冲压缩再取对数, 得到时间-角度谱, 横轴表示接收角度, 纵轴表示时间, 单位为ms. 由于圆柱两个端面产生的镜反射回波, 在光滑圆柱时域图(图10(a))中出现两条相交的亮线. 对于环形凹槽圆柱, 不同凹槽单元由于接收角度、距离不同产生时延, 所以在空间上分离出6个强亮区, 在时域图中表现为三个强亮点. 由第2节可知, 声波偏转具有一定宽度, 所以在凹槽圆柱时域图中出现6个具有一定宽度的亮线, 如图10(b)所示, 黑色、红色和绿色虚线框分别对应声波偏转方向θ3 = 45°, θ5 = 60°和θ1 = 75°的环形凹槽单元的散射回波.

      图  10  时间-角度谱 (a)刚性圆柱; (b)环形凹槽单元组合圆柱

      Figure 10.  Time-angle spectra: (a) Rigid cylinder; (b) annular groove cylinder.

      对水听器接收到的凹槽圆柱时域信号进行宽带信号目标强度处理得到其目标强度的频率-角度谱, 如图11(b)所示, 凹槽圆柱散射声场空间指向性发生变化, 6个预设角度在其主瓣偏转方向有所增强, 且存在干涉效应, 这是由于偏转声波具有一定宽度, 且凹槽单元的周期性结构之间的干涉. 由于具有相同深度梯度方向的环形凹槽单元之间存在相位差发生干涉, 在45°—75°和105°—135°角度范围内出现很强的条纹, 从几何关系(图11(c))点A、点B和点C (各凹槽单元几何中心)推导了预测干涉条纹的近似公式:

      图  11  凹槽单元组合圆柱目标强度频率-角度谱 (a)数值计算结果; (b)实验结果; (c) 6个凹槽单元圆柱二维图

      Figure 11.  Frequency-angle spectra of target strength for annular groove cylinder: (a) Numerical result; (b) experimental result; (c) 2-D geometry of the annular groove cylinder.

      $f = \frac{{nc}}{{l\cos \theta }},$

      其中n = 1, 2, 3 ···, $l = {\rm{4}}Nd$, 图11中黑色虚线是根据(6)式预测的干涉条纹. 从图11中提取出频率f = 80 kHz凹槽圆柱目标强度随接收角度变化的曲线, 在45°—75°和105°—135°角度范围内凹槽圆柱散射声场均衡化, 90°方向散射声场较弱, 实验结果(红色虚线)与数值计算结果(黑色实线)一致. 由于发射信号直接到达水听器的信号和水听器接收到的圆柱散射波信号在接近圆柱两端一定范围内发生混叠, 因此实验和数值结果在接近0°和180°范围内存在一定的误差. 由于实验用模型加工精度的限制和实验测量引起的误差使得图12中数值计算结果和实验结果存在一定的误差.

      图  12  频率f = 80 kHz凹槽圆柱目标强度指向性

      Figure 12.  The normalized directivity of target strength for cylinder with six annular groove units at f = 80 kHz

      实验中获取频率f = 80 kHz圆柱和凹槽圆柱不同角度的目标强度, 如图13所示, 黑色线代表圆柱实验结果, 红色虚线为凹槽圆柱实验结果. 当入射波垂直于圆柱方向入射时, 目标强度指向性比较尖锐, 且在90°目标强度值最大, 具有很强的方向性. 环形凹槽单元组合圆柱目标强度在90°方向约减小了12 dB, 在45°—75°和105°—135°的角度范围内增大. 对比圆柱和环形凹槽圆柱目标强度指向性显示: 具有不同梯度环形凹槽单元组合圆柱正横散射声波均匀偏转到预定的空间范围内, 使得散射声场均衡化, 改变了圆柱整体的散射特征.

      图  13  频率f = 80 kHz圆柱和凹槽圆柱目标强度实验结果对比

      Figure 13.  The comparison of target strength in experiment between cylinder and annular groove cylinder at f = 80 kHz

    • 本文基于有限元方法对具有环形凹槽结构圆柱的声散射机理进行了分析. 详细讨论了环形凹槽结构参数如占空比、梯度等对圆柱散射声场空间指向性的影响规律: 圆柱正横散射目标强度随着占空比增大单调下降; 而散射声波偏转方向目标强度在占空比δ < 30%时随着占空比增大而增大, 当占空比$\delta > 30\% $, 其目标强度基本不变. 当入射声波一定时, 通过设计适当的环形凹槽梯度可以使圆柱反向散射声波发生定向偏转, 且声波偏转方向与频率无关. 另外, 对多个深度梯度环形凹槽单元组合圆柱进行数值计算和实验验证, 研究发现通过设计不同梯度环形凹槽单元组合圆柱可以使得圆柱反向散射声波目标强度均匀偏转到一定区域内, 使得散射声场均衡化, 改变了圆柱整体的散射特征. 因此通过改变环形凹槽梯度、多个深度梯度环形凹槽单元间不同组合等结构参数调控目标散射声场, 为水下目标声隐身设计提供新思路.

参考文献 (30)

目录

    /

    返回文章
    返回