CN108009376B - 基于内嵌信号调节板阵列舷侧阵宽频信号增强及吸声方法 - Google Patents

基于内嵌信号调节板阵列舷侧阵宽频信号增强及吸声方法 Download PDF

Info

Publication number
CN108009376B
CN108009376B CN201711396852.7A CN201711396852A CN108009376B CN 108009376 B CN108009376 B CN 108009376B CN 201711396852 A CN201711396852 A CN 201711396852A CN 108009376 B CN108009376 B CN 108009376B
Authority
CN
China
Prior art keywords
array
sound
signal
normal
covering layer
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201711396852.7A
Other languages
English (en)
Other versions
CN108009376A (zh
Inventor
张燕妮
潘杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Publication of CN108009376A publication Critical patent/CN108009376A/zh
Application granted granted Critical
Publication of CN108009376B publication Critical patent/CN108009376B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Transducers For Ultrasonic Waves (AREA)

Abstract

本发明公开了一种基于内嵌信号调节板阵列的舷侧阵宽频信号增强及吸声方法,包括建立内嵌信号调节板阵列的舷侧阵声呐基阵结构模型,获取相关参数;求解该复合结构在平面声波激励下的各层位移与应力振幅系数向量;计算水听器各单元处的信号增益及结构表面吸声系数;确定水听器单元最佳布放位置。通过在舷侧阵去耦覆盖层及透声覆盖层之间嵌入周期排列的信号调节板阵列,利用各信号调节板单元、阵列与不同性质覆盖层子结构的共振耦合及汇聚作用克服了低频水声信号衰减及低频宽带声吸收差的问题,可实现超宽频率范围内(0–21.3KHz)水听器单元接收信号的显著增强,同时能在多个宽频范围内增加吸声性能从而降低主动声呐的探测距离。

Description

基于内嵌信号调节板阵列舷侧阵宽频信号增强及吸声方法
技术领域
本发明属于结构声与水声交叉技术领域,具体涉及舷侧阵声呐的信号增强和吸声性能的综合设计方法,可从综合考虑宽频信号增强和吸声性能的角度指导舷侧阵结构方案的设计。
背景技术
反潜作战作为当今世界海军公认的最大难题,亟需解决的主要技术问题是探测安静性日益提高的潜艇。随着消声瓦等技术的成熟和运用,潜艇的辐射噪声大大降低,而低频范围内目标舰艇的辐射噪声较大且难以被消声瓦等技术所控制,另一方面,鉴于潜艇对中高频信号的反射回波很小,因此低频声呐成为探测声呐的发展方向。舷侧阵声呐作为近年来大力发展的一种低频声呐,安装于水下航行器两舷侧,在远程范围内具备高探测灵敏性,且具有孔径大、兼顾测向和测距等优点,已成为衡量潜艇先进性的重要指标。
增强舷侧阵声呐水听器拾取信号的方法有基于声压水听器的信号拾取及增强、基于矢量水听器的信号增强以及阵列信号处理的方法。声压水听器法将信号调节板(SCP)大面积或完全敷设在去耦覆盖层表面来增强水听器拾取的声压信号。该方式对SCP厚度要求较高,故使得舷侧阵整阵重量大大增加,且阵元增益局限在6分贝且低频增益衰减较大。矢量水听器法依靠直接拾取柔性障板表面较大的振速(或加速度)矢量信号来增强水听器阵增益,虽省去了刚性信号调节板的引入,但仍存在阵元增益有限且低频信号衰减问题。此外,增加水听器阵信噪比还可借助阵列信号处理的方法,然而舷侧阵信号处理涉及强平台噪声下的远场弱信号检测,由于平台噪声产生机理和传播信道的复杂性,使得其信号检测性能有限。
另一方面,吸声系数是表征水下结构声目标强度的物理量。为保证潜艇的作战性能,在满足其声探测的同时,必须满足最基本的声目标强度或吸声性能的要求。所以,舷侧阵结构设计的核心问题之一就是如何兼顾其信号增强和吸声性能。关于增强水下结构的吸声性能国内外已经发展了大量的方法,包括采用多层覆盖层结构和覆盖层中嵌入无规或周期性填充物,如空腔、金属球及多孔性包含物及声子晶体,其吸声机理可归因为多重散射、波形变换和谐振吸收等。这些方法主要对高频或低频窄带起作用,由于水下声波在低频范围的波长远大于同频率空气声的波长,水下低频宽带范围的声吸收仍是难题。至今还没有一种理想的方法,可在增强超宽频水声信号的基础上同时兼顾低频宽带吸声性能。若能有效降低水听器阵的低频信号衰减、实现宽频水声信号增强且兼顾其宽频吸声性能,同时减轻整阵重量,则可大大增加舷侧阵声呐的探测灵敏性及隐身性能,从而提高其作战性能。
发明内容
针对现有技术所存在的上述技术的不足,本发明目的在于提供一种基于内嵌信号调节板阵列来增强舷侧阵信号增强及吸声方法,从综合考虑增强信号和减低目标强度的角度来指导舷侧阵结构的设计。该方法能够实现超宽频水声信号增益显著增强以大幅改善舷侧振声呐的探测性能,并能增加结构表面声吸收以兼顾隐身性能。
本发明是通过下述技术方案来实现的。
本发明的基于内嵌信号调节板阵列的舷侧阵宽频信号增强及吸声方法,包括如下步骤:
(1)建立内嵌信号调节板阵列的舷侧阵声呐基阵结构模型,获取声激励的振幅、频率及入射角参数,确立入射声压在x-z平面的空间分布;获取复合结构的几何材料参数,并得到对应复纵波波速、复剪切波速、弯曲刚度和吻合频率;
(2)求解该复合结构在平面声波激励下的各层位移与应力振幅系数向量;
运用各层之间应力平衡与位移连续性,并根据周期性结构的简谐波展开法及各简谐波的正交性,得到平面声压入射下的复合结构耦合方程的矩阵形式;利用矩阵求逆求解,可得各层位移与应力的振幅系数向量;
(3)计算水听器各单元处的信号增益及结构表面吸声系数
计算水听器各单元处的法向应力分布,从而可得水听器各单元信号增益;计算结构表面的总声压和法向质点振速,可得结构表面入射声能及吸收声能,从而可得结构表面吸声系数;
(4)确定水听器单元最佳布放位置
绘制信号增益关于空间坐标的高度图,确定最大信号增益的空间坐标,即对应水听器单元的最佳布放位置。
优选的,所述步骤(1)中,结构模型从内到外依次为空气层、艇板、去耦覆盖层、周期信号调节板阵列、水听器阵列及透声覆盖层;
选取坐标x轴为去耦覆盖层与透声覆盖层分界线,由于信号调节板采用薄板结构,其厚度影响在较低频范围可忽略,选取坐标原点在第0个信号调节板的左端点与去耦覆盖层的交点,故第0个信号调节板左右两边分别有N个信号调节板。
优选的,所述步骤(1)中,根据水中入射平面波声激励的振幅Pi、频率f、传播声速co,入射角度θi,可确立入射声压在x–z平面的空间分布
Figure GDA0002229487060000041
其中,j为虚数单位;ko=2πf/co
优选的,所述步骤(1)中,确立复合结构的去耦覆盖层、透声覆盖层的几何材料参数,包括密度ρ、厚度h、纵波波速cd0、剪切波速cs0及对应的纵波损耗因子ηd和剪切波损耗因子ηs,则对应复纵波波速为
Figure GDA0002229487060000042
复剪切波速为
Figure GDA0002229487060000043
其中,j为虚数单位;
确定复合结构的艇板、信号调节板的几何材料参数包括密度ρ、杨氏模量E、泊松比υ及厚度h;则弯曲刚度为D=Eh3/12(1-υ2),吻合频率为
Figure GDA0002229487060000044
co为水中传播声速;设信号调节板长度为L,两两间距为d,则空间周期为l=L+d,信号调节板单元个数为2N+1。
优选的,所述步骤(2)中,求解该复合结构在入射平面声波下的各层位移与应力的系数向量,包括下述步骤:
2a)运用现有的各均匀层分界面即空气层与艇板、艇板与去耦覆盖层、透声覆盖层与水分界面的应力平衡与位移连续性,并运用非均匀层分界面即去耦覆盖层与透声覆盖层分界面z=0法向及剪切位移连续,分别得该界面法向及剪切位移连续方程;
2b)并运用去耦覆盖层与透声覆盖层分界面z=0法向及剪切应力平衡可得该界面法向及剪切应力平衡方程;
2c)根据周期性结构的简谐波展开法,将法向及剪切位移和法向及剪切应力展开;
2d)将法向及剪切位移和法向及剪切应力展开式及ΔHn(x)的傅里叶级数展开式代入各界面法向及剪切位移连续方程及法向及剪切应力平衡方程,并利用简谐波的正交性得到平面声压入射下的复合结构耦合方程的矩阵形式;
2e)利用矩阵求逆法求解方程即可得各层中位移与应力振幅的系数向量。
优选的,所述步骤(3)中,计算水听器各单元处的信号增益及结构表面吸声系数,包括下述步骤:
3a)将求得的应力与位移振幅系数向量值代入应力展开式,可得透声覆盖层中水听器单元所在深度hA的法向应力
Figure GDA0002229487060000051
3b)将
Figure GDA0002229487060000052
代入计算信号增益SR公式;
3c)根据透声覆盖层深度hA的法向应力
Figure GDA0002229487060000053
可得其上表面深度为hc2的总声压;
3d)将求得的应力与位移振幅系数向量值代入法向及剪切位移式,可得透声覆盖层上表面的法向质点振速;
3e)根据所得的结构表面的总声压和法向质点振速,可得结构表面入射和吸收的能量;
据此确定出结构表面吸声系数为αSCPs=Πabin
优选的,所述步骤(4)中,确定水听器单元的最佳布放位置见如下步骤:
4a)用MATLAB绘制信号增益关于空间坐标x的高度图;
4b)根据高度图z坐标的高低,判定最大信号增益位置,z坐标越高,信号增益越大;
4c)最大信号增益对应的x坐标即为水听器单元的最佳布放位置。
本发明与现有技术相比,具有以下优点:
1.针对舷侧阵存在的低频信号衰减及低频宽带声吸收差,本发明利用周期性复合结构的共振耦合及汇聚原理,将现有信号调节板分解优化为周期排列的信号调节板单元,利用各信号调节板单元、阵列与不同性质覆盖层子结构的有效耦合及汇聚作用克服了低频水声信号衰减及低频宽带声吸收差的问题。
2.在进行结构设计时,同时考虑水下信号增强及吸声性能,可同时实现0-21.3kHz超宽频水声信号显著增强及多个宽频范围的声吸收显著改善。本方法可从综合考虑宽频信号增强和吸声性能的角度指导舷侧阵结构方案的设计。
3.本发明在提高水声性能的同时减轻了结构重量。
附图说明
图1为本发明结构参数与坐标图;
图2为水声信号增益曲线(x为距离SCP左端点的距离);
图3为结构表面法向吸声系数频响曲线;
图4为结构表面吸声系数频响曲线随入射角变化的规律;
图5为5Ω=0.05时信号增益在一个整周期内的空间分布,其中x′为距离SCP中心的距离。
图1中:1、空气层,2、艇板,3、去耦覆盖层,4、信号调节板阵列,5、水听器线阵,6、透声覆盖层,7、海水。
具体实施方式
下面结合附图和实施例对发明作进一步的详细说明,但并不作为对发明做任何限制的依据。
参照图1,本发明基于内嵌信号调节板阵列的舷侧阵宽频信号增强及吸声方法,具体步骤如下:
步骤一:建立内嵌信号调节板阵列的舷侧阵声呐基阵结构模型,获取声激励的参数及该复合结构的几何材料参数
1.1)建立内嵌信号调节板阵列的舷侧阵声呐基阵结构模型,从内到外依次为空气层1、艇板2、去耦覆盖层3、周期信号调节板阵列4、水听器阵列5及透声覆盖层6;其中,艇板2上下表面分别与去耦覆盖层3和艇内空气层1接触,去耦覆盖层3胶接在艇板2表面,去耦覆盖层3上表面粘贴有透声覆盖层6,两层覆盖层之间粘贴有周期性分布的信号调节板阵列4,水听器线阵5内嵌于透声覆盖层6内、信号调节板阵列4附近。透声覆盖层6外表面与海水7相连,如图1所示。在入声声波透声至其表面时,利用周期信号调节板阵列4与两层性质迥异覆盖层的共振耦合及汇聚作用放大黏贴其上的水听器阵列5的增益;
1.2)根据水中入射平面波声激励的振幅为Pi、频率f、传播声速co,入射角度θi,可确立入射声压在x-z平面的空间分布
Figure GDA0002229487060000071
其中,j为虚数单位;ko=2πf/co
1.3)确立去耦覆盖层的几何材料参数,采用阻抗较小的软橡胶材料,其动力学采用弹性理论模拟。其密度、厚度、纵波波速、剪切波速及对应的纵波损耗因子和剪切波损耗因子分别为ρ1、h1、cd01、cs01、ηd和ηs,则可得其对应复纵波波速和复剪切波速分别为
Figure GDA0002229487060000072
Figure GDA0002229487060000073
cd01为去耦覆盖层纵波波速;cs01为去耦覆盖层剪切波速;
1.4)确立透声覆盖层的几何材料参数。透声覆盖层采用阻抗与水接近的透声橡胶层,其密度、厚度、纵波波速、剪切波速及对应的纵波损耗因子和剪切波损耗因子分别为ρ3、h3、cd03、cs03、ηd3和ηs3,则可得对应复纵波波速和复剪切波速分别为
Figure GDA0002229487060000081
Figure GDA0002229487060000082
cd03为透声覆盖层纵波波速;cs03为透声覆盖层剪切波速;
1.5)确定艇板的几何材料参数。艇板板采用钢板结构,其密度、杨氏模量、泊松比及厚度分别为ρ、E、υ及h。弯曲刚度D=Eh3/12(1-υ2);其吻合频率为
Figure GDA0002229487060000083
1.6)确定信号调节板SCP的几何物理参数。SCP采用薄钢板,其密度、杨氏模量、泊松比及厚度分别为ρ、E、υ及h。弯曲刚度D=Eh3/12(1-υ2)。设SCP长度为L,两两间距为d,则空间周期为l=L+d。SCP单元个数为2N+1,则第n个SCP对应的位置函数可用ΔHn(x)=H(x-nl)-H(x-L-nl)表征,其中H(x)是Heaviside函数,n=–(2N+1)…–2,–1,0,1,2…2N+1;
1.7)坐标x轴选取在去耦覆盖层与透声覆盖层分界线,由于SCP采用薄板结构,其厚度影响在较低频范围可忽略,选取坐标原点在第0个SCP的左端点与去耦覆盖层的交点,故第0个SCP左右两边分别有N个SCP。
步骤二:求解该复合结构在平面声波入射下的各层位移与应力振幅系数向量
2.1)运用现有的各均匀层分界面应力平衡与位移连续性,包括:
运用空气与艇板分界面z=-(hc1+h)的位移连续可得位移连续方程:
Figure GDA0002229487060000084
其中,wh(x)和pt(x,z)分别为艇板的弯曲位移及空气中的透射声压;ω=2πf;
运用艇板与去耦覆盖层分界面z=-hc1的法向与剪切位移连续可得该界面位移连续方程:
Figure GDA0002229487060000085
Figure GDA0002229487060000091
其中,
Figure GDA0002229487060000092
Figure GDA0002229487060000093
分别为去耦覆盖层的法向及剪切位移;uh(x)为艇板x方向的纵向位移;
运用艇板与去耦覆盖层分界面z方向剪力平衡与弯矩为0可得该界面应力平衡方程:
Figure GDA0002229487060000094
Figure GDA0002229487060000095
其中,
Figure GDA0002229487060000096
分别为去耦覆盖层法向及剪切应力;
运用透声覆盖层与水分界面z=hc2法向及剪切应力平衡分别可得该界面应力平衡方程:
Figure GDA0002229487060000097
Figure GDA0002229487060000098
其中,
Figure GDA0002229487060000099
为透声覆盖层法向及剪切应力;pT(x,z)为水中总声压;
运用透声覆盖层与水分界面z=hc2法向位移连续可得该界面位移连续方程:
Figure GDA00022294870600000910
其中,
Figure GDA00022294870600000911
为透声覆盖层法向位移;
并运用非均匀层分界面即去耦覆盖层与透声覆盖层分界面z=0法向及剪切位移连续分别得该界面法向及剪切位移连续方程:
Figure GDA00022294870600000912
Figure GDA00022294870600000913
其中,
Figure GDA0002229487060000101
Figure GDA0002229487060000102
分别为去耦覆盖层和透声覆盖层法向位移;
Figure GDA0002229487060000103
Figure GDA0002229487060000104
分别为去耦覆盖层和透声覆盖层的剪切位移;ΔHn(x)=H(x-nl)-H(x-L-nl),H(x)是Heaviside函数,n=–(2N+1)…–2,–1,0,1,2…2N+1,l为空间周期;
2.2)并运用非均匀层分界面即去耦覆盖层与透声覆盖层分界面z=0法向及剪切应力平衡可得该界面法向及剪切应力平衡方程:
Figure GDA0002229487060000105
Figure GDA0002229487060000106
其中,
Figure GDA0002229487060000107
分别为去耦覆盖层法向及剪切应力,
Figure GDA0002229487060000108
分别为透声覆盖层法向及剪切应力;ω=2πf,h2、m2、D2分别为信号调节板的厚度、单位长度的质量及弯曲刚度;
2.3)根据周期性结构的简谐波展开法,将法向及剪切位移和法向及剪切应力展开为:
Figure GDA0002229487060000109
Figure GDA00022294870600001010
Figure GDA00022294870600001011
Figure GDA00022294870600001012
其中i=1,2分别代表去耦覆盖层和透声覆盖层;km=ko sinθi+2πm/l;m为简谐波阶数,ko为水中声波数;
Figure GDA0002229487060000111
分别为法向和剪切位移及法向及剪切应力的振幅;
2.4)将法向及剪切位移和法向及剪切应力展开式(13)–(16)及ΔHn(x)的傅里叶级数展开式
Figure GDA0002229487060000112
cn=L/lejnπL/lsin c(nπL/l)代入各界面法向及剪切位移连续方程及法向及剪切应力平衡方程中,并利用简谐波的正交性得到平面声压入射下的复合结构耦合方程的矩阵形式;
Figure GDA0002229487060000113
其中,
Figure GDA0002229487060000114
为耦合方程系数矩阵;
Figure GDA0002229487060000115
为激励向量:
Figure GDA0002229487060000116
为应力和位移振幅系数向量;
Figure GDA0002229487060000117
为信号调节板阵列与复合结构其它子结构的耦合矩阵;
2.5)利用矩阵求逆法求解方程(17)可得
Figure GDA0002229487060000118
即可得各层中位移与应力振幅的系数向量。
步骤三,计算水听器各单元处信号增益及结构表面吸声系数
3.1)将求得的应力与位移振幅系数向量的值代入式(15),可计算得透声覆盖层中水听器所在深度hA的法向应力
Figure GDA0002229487060000119
3.2)水听器单元信号增益SR通过下式计算:
Figure GDA00022294870600001110
3.3)根据透声覆盖层深度hA的法向应力
Figure GDA00022294870600001111
可得其上表面深度为hc2的总声压
Figure GDA00022294870600001112
Figure GDA00022294870600001113
为透声覆盖层上表面的法向应力;
3.4)将求得的应力与位移振幅系数向量值代入式(13),可得透声覆盖层上表面的法向质点振速
Figure GDA00022294870600001114
Figure GDA00022294870600001115
为透声覆盖层上表面的法向位移;
3.5)根据所得的结构表面的总声压和法向质点振速,可得结构表面入射和吸收的能量分别为:
Figure GDA0002229487060000121
Figure GDA0002229487060000122
其中
Figure GDA0002229487060000123
Figure GDA0002229487060000124
的复共轭;ρo为水的密度;
将Πin和Πab代入结构表面吸声系数计算公式:
αSCPs=Πabin (21)
其中θi=90°对应的αSCPs为法向吸声系数,否则为斜向吸声系数。
步骤四,确定水听器单元最佳布放位置
4.1)用MATLAB绘制信号增益关于空间坐标x的高度图;
4.2)根据高度图z坐标的高低,判定最大信号增益位置,z坐标越高,信号增益越大;
4.3)最大信号增益对应的x坐标即为水听器单元的最佳布放位置。
本发明优点可通过以下仿真实验进一步说明:
1.确立声激励的参数及所述复合结构的几何材料参数
(1)本实验以振幅为1个单位强度的平面简谐声波垂直入射为例即Pi=1 Pa,θi=90°。所述艇板及信号调节板采用金属钢板、去耦覆盖层采用多孔软橡胶,透声覆盖层采用透声橡胶。其中多孔橡胶密度ρ1=500kg/m3,纵波波速
Figure GDA0002229487060000125
剪切波速为
Figure GDA0002229487060000126
透声橡胶密度ρ1=990kg/m3,纵波波速为
Figure GDA0002229487060000127
剪切波速为
Figure GDA0002229487060000128
各子结构几何性质分别为:艇板厚度h及去耦覆盖层厚度hc1均为0.01m,透声覆盖层厚度hc2=0.02m,信号调节板SCP单元厚度h2=0.0075m,宽0.1m。
取SCP单元宽度L为0.1m,单元间距d也为0.1m,考虑从213Hz(0.01fc)到基板吻合频率21304Hz(fc)下的宽频响应,频率分辨率为0.01fc。其中入射声信号的性质和内嵌SCP阵列的几何和结构参数如表1所示。
表1本发明的激励参数和结构参数
Figure GDA0002229487060000131
2.求解该复合结构在平面声波入射下的各层位移与应力振幅系数向量
将各结构及工作参数代入平面声压激励下复合结构的耦合方程(17),利用矩阵求逆法即可得各层应力与位移振幅系数向量。
3.计算水听器各单元处的信号增益
(1)将所得应力与位移振幅系数向量代入式(15)可得透声覆盖层中水听器所在深度hA法向应力
Figure GDA0002229487060000132
(2)将得到的法向应力代入式(18),可计算得水听器单元信号增益SR。如图2所示,其中细虚线代表无信号调节板(SCP)阵列时本发明结构的信号增益随归一化频率的变化曲线,其他线条为有SCPs不同空间位置的增益曲线。从图2可以看出,对应所给SCP阵列的参数及各空间位置,该舷侧阵结构的信号增益在超宽频范围内(0.05fc–fc)均有所提高。尤其在每个SCP两个端点附近,如图2中粗点线所示,信号增益在超宽频范围0.05fc–fc平均增加了20分贝以上。其中在Ω=0.05,0.17,0.29等附近增益增幅超过40分贝。
4.计算结构表面吸声系数
(1)根据透声覆盖层深度hA的法向应力值
Figure GDA0002229487060000141
可得其上表面深度为hc2的总声压值
Figure GDA0002229487060000142
(2)将求得的应力与位移振幅系数向量值代入式(13),可得透声覆盖层上表面的法向质点振速
Figure GDA0002229487060000143
(3)将所得结构表面的总声压和法向质点振速值代入式(19)和式(20),可得结构表面入射和吸收的能量Πin和Πab;进而根据式(21)可得结构表面吸声系数。
据此计算所得垂直入射时结构表面法向吸声系数αSCPs如图3所示。其中实线表示无SCP阵列的吸声系数,虚线表示内嵌了周期性SCP阵列的吸声系数。从图3可以看出,覆盖层之间无内嵌SCP阵列时吸声系数存在三个峰值,分别对应频率Ω=0.028,0.35和0.71;引入内嵌SCP阵列后,除了以上三个峰值频率,吸声系数在所取整个频带范围内有所增加,尤其在Ω=0.05,0.16,0.29,0.53及0.87附近增加超过20%,在Ω=0.29附近增加超过50%。
改变入射方向,可得不同入射角θi下的斜吸声系数,如图4所示。从图4可以看出,只要入射角≤60°,斜入射条件下的吸声系数比垂直入射下的吸声系数在中高频(f/fc≥0.35)有显著增加。
5.确定水听器单元最佳布放位置
为确定水听器单元的最佳布放位置,以Ω=0.05为例,用MATLAB绘制信号增益关于空间坐标x的高度图,如图5所示。在图5中,高度图z坐标的高低对应信号增益幅值的大小,z坐标越高,对应信号增益越大。对所选频率,最大信号增益出现在距离SCP的两个端点x=±2x10-4m处,即为水听器单元最佳布放位置。
以上仿真实验可以看出,本发明提出的基于内嵌信号调节板阵列的舷侧阵超宽频信号增强及吸声方法,利用各信号调节板单元、阵列与不同性质覆盖层子结构的有效耦合及汇聚作用,可以同时实现低频及宽频范围内信号增益及吸声性能显著增加。同时,在工程实践中,可运用本发明的方法,确定该复合结构的结构参数尤其是信号调节板阵列的几何材料参数,从综合考虑信号增益及吸声性能的角度指导舷侧阵声呐结构方案的设计。
本发明并不局限于上述实施例,在本发明公开的技术方案的基础上,本领域的技术人员根据所公开的技术内容,不需要创造性的劳动就可以对其中的一些技术特征作出一些替换和变形,这些替换和变形均在本发明的保护范围内。

Claims (7)

1.基于内嵌信号调节板阵列舷侧阵宽频信号增强及吸声方法,其特征在于,包括如下步骤:
(1)建立内嵌信号调节板阵列的舷侧阵声呐基阵结构模型,包括从内到外依次为空气层(1)、艇板(2)、去耦覆盖层(3)、周期信号调节板阵列(4)、水听器阵列(5)及透声覆盖层(6);其中,艇板(2)上下表面分别与去耦覆盖层(3)和艇内空气层(1)接触,去耦覆盖层(3)胶接在艇板(2)表面,去耦覆盖层(3)上表面粘贴有透声覆盖层(6),两层覆盖层之间粘贴有周期性分布的信号调节板阵列(4),水听器线阵(5)内嵌于透声覆盖层(6)内、信号调节板阵列(4)附近;
获取声激励的振幅、频率及入射角参数,确立入射声压在x-z平面的空间分布;获取舷侧阵声呐基阵结构模型的复合结构的几何材料参数,并得到对应复纵波波速、复剪切波速、弯曲刚度和吻合频率;
(2)求解该复合结构在平面声波激励下的各层位移与应力振幅系数向量
运用各层之间应力平衡与位移连续性,并根据周期性结构的简谐波展开法及各简谐波的正交性,得到平面声压入射下的复合结构耦合方程的矩阵形式;利用矩阵求逆求解,可得各层位移与应力的振幅系数向量;
(3)计算水听器各单元处的信号增益及结构表面吸声系数
计算水听器各单元处的法向应力分布,从而可得水听器各单元信号增益;计算结构表面的总声压和法向质点振速,可得结构表面入射声能及吸收声能,从而可得结构表面吸声系数;
(4)确定水听器单元最佳布放位置
绘制信号增益关于空间坐标的高度图,确定最大信号增益的空间坐标,即对应水听器单元的最佳布放位置。
2.根据权利要求1所述的基于内嵌信号调节板阵列舷侧阵宽频信号增强及吸声方法,其特征在于,所述步骤(1)中,结构模型从内到外依次为空气层、艇板、去耦覆盖层、周期信号调节板阵列、水听器阵列及透声覆盖层;
选取坐标x轴为去耦覆盖层与透声覆盖层分界线,由于信号调节板采用薄板结构,其厚度影响在较低频范围可忽略,选取坐标原点在第0个信号调节板的左端点与去耦覆盖层的交点,故第0个信号调节板左右两边分别有N个信号调节板。
3.根据权利要求1所述的基于内嵌信号调节板阵列舷侧阵宽频信号增强及吸声方法,其特征在于,所述步骤(1)中,根据水中入射平面波声激励的振幅Pi、频率f、传播声速co、入射角度θi,可确立入射声压在x-z平面的空间分布
Figure FDA0002229487050000021
其中,j为虚数单位;ko=2πf/co
4.根据权利要求1所述的基于内嵌信号调节板阵列舷侧阵宽频信号增强及吸声方法,其特征在于,所述步骤(1)中,确立复合结构的去耦覆盖层、透声覆盖层的几何材料参数,包括密度ρ、厚度h、纵波波速cd0、剪切波速cs0及对应的纵波损耗因子ηd和剪切波损耗因子ηs,则对应复纵波波速为
Figure FDA0002229487050000022
复剪切波速为
Figure FDA0002229487050000023
其中,j为虚数单位;
确定复合结构的艇板、信号调节板的几何材料参数包括密度ρ、杨氏模量E、泊松比υ及厚度h;则弯曲刚度为D=Eh3/12(1-υ2),吻合频率为
Figure FDA0002229487050000024
co为水中传播声速;设信号调节板长度为L,两两间距为d,则空间周期为l=L+d,信号调节板单元个数为2N+1。
5.根据权利要求4所述的基于内嵌信号调节板阵列舷侧阵宽频信号增强及吸声方法,其特征在于,所述步骤(2)中,求解该复合结构在入射平面声波下的各层位移与应力的系数向量,包括下述步骤:
2a)运用现有的各均匀层分界面即空气层与艇板、艇板与去耦覆盖层、透声覆盖层与水分界面的应力平衡与位移连续性,并运用非均匀层分界面即去耦覆盖层与透声覆盖层分界面z=0法向及剪切位移连续,分别得该界面法向及剪切位移连续方程:
Figure FDA0002229487050000031
Figure FDA0002229487050000032
其中,
Figure FDA0002229487050000033
Figure FDA0002229487050000034
分别为去耦覆盖层和透声覆盖层的法向位移;
Figure FDA0002229487050000035
Figure FDA0002229487050000036
分别为去耦覆盖层和透声覆盖层的剪切位移;ΔHn(x)=H(x-nl)-H(x-L-nl),H(x)是Heaviside函数,n=–(2N+1)…–2,–1,0,1,2…2N+1,l为空间周期;
2b)并运用去耦覆盖层与透声覆盖层分界面z=0法向及剪切应力平衡可得该界面法向及剪切应力平衡方程:
Figure FDA0002229487050000037
Figure FDA0002229487050000038
其中,
Figure FDA0002229487050000039
分别为去耦覆盖层法向及剪切应力,
Figure FDA00022294870500000310
分别为透声覆盖层法向及剪切应力;ω=2πf,h2、m2、D2分别为信号调节板厚度、单位长度的质量及弯曲刚度;
2c)根据周期性结构的简谐波展开法,将法向及剪切位移和法向及剪切应力展开为:
Figure FDA0002229487050000041
Figure FDA0002229487050000042
Figure FDA0002229487050000043
Figure FDA0002229487050000044
其中i=1,2分别代表去耦覆盖层和透声覆盖层;km=ko sinθi+2πm/l,m为简谐波阶数,ko为水中声波数;
Figure FDA0002229487050000045
Figure FDA0002229487050000046
分别为法向和剪切位移、法向和剪切应力的振幅;
2d)将法向及剪切位移和法向及剪切应力展开式(5)–(8)及ΔHn(x)的傅里叶级数展开式
Figure FDA0002229487050000047
cn=L/lejnπL/lsinc(nπL/l)代入各界面法向及剪切位移连续方程(1)、(2)及法向及剪切应力平衡方程(3)、(4)中,并利用简谐波的正交性得到平面声压入射下的复合结构耦合方程的矩阵形式;
Figure FDA0002229487050000048
其中,
Figure FDA0002229487050000049
为耦合方程系数矩阵;
Figure FDA00022294870500000410
为激励向量;
Figure FDA00022294870500000411
为式(5)-(8)中各应力和位移振幅的系数向量;
Figure FDA00022294870500000412
为信号调节板阵列与复合结构子结构的耦合矩阵;
2e)利用矩阵求逆法求解方程(9)即可得各层中位移与应力振幅的系数向量为:
Figure FDA00022294870500000413
6.根据权利要求5所述的基于内嵌信号调节板阵列舷侧阵宽频信号增强及吸声方法,其特征在于,所述步骤(3)中,计算水听器各单元处的信号增益及结构表面吸声系数,包括下述步骤:
3a)将求得的应力与位移振幅系数向量值代入应力展开式(7),可得透声覆盖层中水听器单元所在深度hA的法向应力
Figure FDA0002229487050000051
3b)将
Figure FDA0002229487050000052
代入计算信号增益SR公式:
Figure FDA0002229487050000053
3c)根据透声覆盖层深度hA的法向应力
Figure FDA0002229487050000054
可得其上表面深度为hc2的总声压
Figure FDA0002229487050000055
Figure FDA0002229487050000056
为透声覆盖层上表面的法向应力;
3d)将求得的应力与位移振幅系数向量值代入式(5),可得透声覆盖层上表面的法向质点振速
Figure FDA0002229487050000057
其中ω=2πf,
Figure FDA0002229487050000058
为透声覆盖层上表面的法向位移;j为虚数单位;
3e)根据所得的结构表面的总声压和法向质点振速,可得结构表面入射和吸收的能量:
Figure FDA0002229487050000059
Figure FDA00022294870500000510
其中,
Figure FDA00022294870500000511
Figure FDA00022294870500000512
的复共轭;ρo为水的密度;
据此确定出结构表面吸声系数为
αSCPs=Πabin (12)
其中,θi=90°对应的αSCPs为法向吸声系数,否则为斜向吸声系数。
7.根据权利要求1所述的基于内嵌信号调节板阵列舷侧阵宽频信号增强及吸声方法,其特征在于,所述步骤(4)中,确定水听器单元的最佳布放位置通过如下步骤实现:
4a)用MATLAB绘制信号增益关于空间坐标x的高度图;
4b)根据高度图z坐标的高低,判定最大信号增益位置,z坐标越高,信号增益越大;
4c)最大信号增益对应的x坐标即为水听器单元的最佳布放位置。
CN201711396852.7A 2017-09-12 2017-12-21 基于内嵌信号调节板阵列舷侧阵宽频信号增强及吸声方法 Active CN108009376B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201710819356 2017-09-12
CN2017108193561 2017-09-12

Publications (2)

Publication Number Publication Date
CN108009376A CN108009376A (zh) 2018-05-08
CN108009376B true CN108009376B (zh) 2020-05-12

Family

ID=62060519

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711396852.7A Active CN108009376B (zh) 2017-09-12 2017-12-21 基于内嵌信号调节板阵列舷侧阵宽频信号增强及吸声方法

Country Status (1)

Country Link
CN (1) CN108009376B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109883985B (zh) * 2019-02-18 2021-03-26 西安科技大学 一种基于分子动力学理论的纳米声学效应分析方法
CN112528399B (zh) * 2020-11-25 2022-09-27 西北工业大学 一种基于能量集中汇聚的水下平台壁板振动宽频控制方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101696955A (zh) * 2009-10-16 2010-04-21 西安交通大学 一种变参数温度场下多孔金属材料吸声性能测试装置
WO2015111939A1 (ko) * 2014-01-27 2015-07-30 한국해양과학기술원 신호조절판의 제작 및 이를 통한 효율적 신호계측 방법
CN105181800A (zh) * 2015-08-20 2015-12-23 浙江大学 基于多通道空时逆滤波技术的声学覆盖层回声降低测量方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101696955A (zh) * 2009-10-16 2010-04-21 西安交通大学 一种变参数温度场下多孔金属材料吸声性能测试装置
WO2015111939A1 (ko) * 2014-01-27 2015-07-30 한국해양과학기술원 신호조절판의 제작 및 이를 통한 효율적 신호계측 방법
CN105181800A (zh) * 2015-08-20 2015-12-23 浙江大学 基于多通道空时逆滤波技术的声学覆盖层回声降低测量方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Underwater sound radiation from an elastically coated infinite plate with periodic inhomogeneities of finite width;Zhang, Yanni 等;《JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA》;20170710;第142卷(第1期);第91-102页 *
水声吸声覆盖层理论与实验研究;王曼;《万方学位论文库》;20050713;第1-114页 *

Also Published As

Publication number Publication date
CN108009376A (zh) 2018-05-08

Similar Documents

Publication Publication Date Title
Lin et al. Three-dimensional sound propagation models using the parabolic-equation approximation and the split-step Fourier method
Mu et al. Sound insulation characteristics of multi-layer structures with a microperforated panel
CN108009376B (zh) 基于内嵌信号调节板阵列舷侧阵宽频信号增强及吸声方法
CN110223667A (zh) 一种复合空腔声障板
CN112040382B (zh) 基于声阻抗梯度匹配层的高频宽带水声换能器
CN112162266A (zh) 一种基于凸优化理论的共形阵二维波束优化方法
CN101718869B (zh) 具有宽带、超宽覆盖性能的平面型声基阵
CN103713283B (zh) 参量接收阵组成基阵的实现方法
Yang et al. Flexural waves in a floating ice sheet: Modeling and comparison with data
WO2023000583A1 (zh) 基于对称介质表面的目标声散射隐身覆盖层及其实现方法
CN105301112A (zh) 一种橡胶类阻尼材料黏弹性动态力学参数的测算方法
CN116184409A (zh) 一种敷设消声瓦目标的收发分置目标强度快速算法
AU2013276738B2 (en) Structure for an underwater vehicle such as a submarine
Li et al. A Sandwich Coating Containing Micro-perforated Panel for Underwater Sound Absorption
CN112951188B (zh) 一种有源式微穿孔板吸声器及提升其低频吸声性能的方法
CN112834623B (zh) 一种基于压缩感知技术的声学覆盖层回声降低测量方法
Liu et al. Analysis of Liquid Inclination on the Absorption Performance of the Cavity-containing Coating
CN105136900B (zh) 一种利用有限大薄板获取流固界面波的方法
CN112528399B (zh) 一种基于能量集中汇聚的水下平台壁板振动宽频控制方法
Hahn et al. Effect of the geoacoustic parameters on the range-frequency interference in shallow-water waveguide
CN112433218B (zh) 一种舰用共形阵虚拟障板的实现方法
CN107995557A (zh) 传感与消声一体化水听器及其实现方法
CN116205104A (zh) 一种敷设消声瓦目标的收发合置目标强度的快速算法
Zhao et al. Study on Sound Absorption Characteristics in Seepage State of Anechoic Coating
Zhu et al. Modeling and characterizing the typical under-ice acoustic channel for the arctic

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant