CN109444861B - 一种平面声呐阵阻抗特性近场声全息法校准测量方法 - Google Patents
一种平面声呐阵阻抗特性近场声全息法校准测量方法 Download PDFInfo
- Publication number
- CN109444861B CN109444861B CN201811504271.5A CN201811504271A CN109444861B CN 109444861 B CN109444861 B CN 109444861B CN 201811504271 A CN201811504271 A CN 201811504271A CN 109444861 B CN109444861 B CN 109444861B
- Authority
- CN
- China
- Prior art keywords
- field
- sound
- calibration
- array
- radiation
- 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.)
- Expired - Fee Related
Links
- 238000001093 holography Methods 0.000 title claims abstract description 70
- 238000000691 measurement method Methods 0.000 title claims abstract description 11
- 230000005855 radiation Effects 0.000 claims abstract description 88
- 238000005259 measurement Methods 0.000 claims abstract description 74
- 238000000034 method Methods 0.000 claims abstract description 49
- 238000012360 testing method Methods 0.000 claims abstract description 43
- 238000004088 simulation Methods 0.000 claims abstract description 18
- 238000004364 calculation method Methods 0.000 claims abstract description 14
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 9
- 238000011426 transformation method Methods 0.000 claims abstract description 3
- 238000001228 spectrum Methods 0.000 claims description 21
- 238000004458 analytical method Methods 0.000 claims description 16
- 239000002245 particle Substances 0.000 claims description 14
- 238000005516 engineering process Methods 0.000 claims description 12
- 230000006870 function Effects 0.000 claims description 10
- 230000009466 transformation Effects 0.000 claims description 10
- 230000000694 effects Effects 0.000 claims description 7
- 230000000903 blocking effect Effects 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 3
- 230000007613 environmental effect Effects 0.000 claims description 3
- 235000021384 green leafy vegetables Nutrition 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000011160 research Methods 0.000 abstract description 13
- 238000011161 development Methods 0.000 abstract description 9
- 238000012544 monitoring process Methods 0.000 abstract description 7
- 238000013461 design Methods 0.000 abstract description 4
- 239000000523 sample Substances 0.000 description 15
- 230000008569 process Effects 0.000 description 5
- 238000002360 preparation method Methods 0.000 description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 230000003993 interaction Effects 0.000 description 3
- 238000007781 pre-processing Methods 0.000 description 3
- 230000035945 sensitivity Effects 0.000 description 3
- 241000282414 Homo sapiens Species 0.000 description 2
- 230000003321 amplification Effects 0.000 description 2
- 238000004140 cleaning Methods 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 238000009877 rendering Methods 0.000 description 2
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000004806 packaging method and process Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000006798 recombination Effects 0.000 description 1
- 238000005215 recombination Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000002791 soaking Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52004—Means for monitoring or calibrating
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
一种平面声呐阵阻抗特性近场声全息法校准测量方法,属于大型平面发射阵阻抗特性近场声全息测量方法的研究技术领域,其技术要点是:平面发射阵阻抗特性计算方法;基于近场声全息的辐射阻抗校准方法;基于BAHIM法近场声全息变换方法和算法;基于近场声全息的辐射阻抗校准方法数值仿真;建立低频平面声呐阵阻抗特性近场声全息理论测试模型。本发明可根据水声监听船主动声呐设备声特性校准的需求,有针对性地开展平面声呐阵阻抗特性近场声全息法校准基本原理研究、近场声全息法校准方案设计、低频声呐阵阻抗特性校准软件开发以及近场声全息法校准测量试验的相关研究工作,以提升低频发射换能器校准系统的校准能力,解决低频声呐阵阻抗特性的评价问题。
Description
技术领域:
本发明涉及大型平面发射阵阻抗特性近场声全息测量方法的研究,具体涉及一种平面声呐阵阻抗特性近场声全息法校准测量方法。
背景技术:
海洋占地球总面积的70%以上,它蕴含着大量的能源和资源。随着世界科技水平的不断进步,人类探索海洋步伐的不断加快,海洋的地位日益提高,人们越来越认识到海洋必将成为人类赖以生存的基础环境。人类探索海洋时离不开船舶,而掌握船舶在海洋中的运行情况则主要是根据水声来监听船舶的主动声呐,若船舶处的水声声呐为低频,则不利于对船舶的监听,导致监听结果不准确,影响对船舶运行情况的判断。
发明内容:
本发明提供了一种平面声呐阵阻抗特性近场声全息法校准测量方法,其可根据水声监听船主动声呐设备声特性校准的需求,有针对性地开展平面声呐阵阻抗特性近场声全息法校准基本原理研究、近场声全息法校准方案设计、低频声呐阵阻抗特性校准软件开发以及近场声全息法校准测量试验的相关研究工作,以提升低频发射换能器校准系统的校准能力,解决低频声呐阵阻抗特性的评价问题。
本发明的平面声呐阵阻抗特性近场声全息法校准测量方法,为实现上述目的所采用的技术方案在于,包括以下步骤:
a、平面发射阵阻抗特性计算方法
发射器在辐射声波的同时会受到声场的作用,表现为贴近振动面的介质对其施加的作用力,该作用力与由振动面施加给介质的力在数值上相等而方向相反,取决于贴近该面的介质中的声压,根据发射面上声压分布沿振动面的积分可以得到介质对振动面的作用力,而声场中的声压和该处介质的质点振速成正比,即:
P=u·Z0 (1)
由以上获得,贴近发射面的介质质点的法向振速和发射面的法向振速完全相同,于是发射面上的声压分布表示为:
Pa=ua·Za (2)
其中ua:发射面上法向振速;Za:贴近发射面的介质波阻抗,
假设声场中有n个声源,各声源的振速分别为u1,u2,u3…un,在这种振速分布情况下,1号源辐射面上场的总作用力为F1,显然为1号源自身辐射场声压作用在自身辐射面上产生的力与其他源辐射场声压作用在1号源辐射面上产生的力之和,即
其中:f11——1号辐射器自身辐射场声压作用在自身辐射面产生的力;
fin——n号辐射器辐射场声压作用在i号辐射器辐射面产生的力,
i号辐射器总的辐射阻抗为
由于声源产生的场和声源强度成正比,在均匀振幅分布情况下,辐射场的声压和辐射面的振速成正比,所以可以把fil写成
fil=Zilul (5)
则可以得到
其中Zil表示l单位振速声场作用于i号源面上的力,称之为l号源对i号源的互辐射阻抗,简称互阻抗,
当辐射器以同样的振速和相位振动,那么由n个阵元组成的基阵的辐射阻抗就是
由以上可知,只要测量得到各阵元表面的复振速、复声压分布以及该阵元表面位置处由其他阵元产生的辐射声压,即可求得各阵元的自/互辐射阻抗;
b、基于近场声全息的辐射阻抗校准方法
无穷大平面上的格林函数分别为:
在狄里克利边界条件下,可以表示为
根据二维傅里叶变换,可以得到S面上声压信号的角谱与声场中任意平面上辐射复声压信号的角谱之间的关系:
其中
对式(12)进行傅里叶逆变换,即可得到声场中任意一个平面上的场点声压:
在纽曼边界条件下,可以表示为
根据欧拉方程可得,
将其代入式(15)即可得到如下表达形式:
对上式两边同时进行二维傅里叶变换,可以得到波数域下的角谱关系式:
其中
因此,利用源平面上的质点振速可以得到声场中任意一个平面上的场点声压:
当源平面上声压或质点振速为未知量时,假设全息面与源面平行,位于zH处,且zH>zS,则可以得到源面上声压角谱表达式:
得到源平面上声压角谱后,可以得到空间中其他平面上(z>zS)的声压角谱分布:
而对于法向振速而言,可以利用式(17)得到重构方程:
经过二维傅里叶逆变换后,则可以通过全息面上的复声压分布得到源面及空间声场中复振速分布,如下所示:
利用上式得到源面上的复振速分布后,按照各阵元所在位置对复振速进行分块处理,得到各阵元源面复振速,计算其源面复声压及其他阵元表面位置处由该阵元产生的辐射声压,从而计算各阵元的自/互辐射阻抗;
c、基于BAHIM法近场声全息变换方法和算法
基于声强测量的宽带声全息技术,简称BAHIM,它通过在全息面上的声强测量来重构全息面相位,进而得到全息面复声压,基本重构方式如下:
声波在理想介质中传播时,能量随着介质质点振动状态沿波的传播方向传播,从一个区域流向另一个区域,引入了介质中能流的概念,定义单位时间通过与能量传播方向垂直的单位面积上的声能为声能流密度,以表示,在时域中,成为瞬时声强,定义为:
瞬时声强可以表示为:
对于三维声场,任意点处的声强可以表示为:
把上式可以推广到包含实部和虚部的复数声强,即
对空间声场中非单频的限带,通过傅立叶变换将其分解成一系列简谐信号的叠加,在任意频率ω处,空间声场任意点(x,y,z)处时均声强的频域表达式:
可以得到:
由于有功声强中包含相位梯度信息,如果知道空间声场的有功声强和声压幅值的分布,则可以求出相位梯度:
对(31)式两边作二维FFT变换,由空间域变到波数域后,全息面上声压相位的波数域形式为
其中,k=ω/c为声波波数,c为介质中的波速,βx(kx,ky)、βy(kx,ky)分别为点(x,y,z)处x方向和y方向声强与声压幅值平方的比值的二维FFT变换,其表达式如下
通过以上的数学推导,得到了全息面上声压相位在波数域的分布,再对(32)式做二维FFT反变换就可以得到空间域全息面上相位的分布,进而得到全息面复声压的分布,便可应用正常的NAH技术;
d、基于近场声全息的辐射阻抗校准方法数值仿真
针对基于近场声全息理论进行声呐阵辐射阻抗校准的测量模型进行数值仿真计算,分析测量方法与重构算法的误差,同时针对任意形状下近场测量数据推算远场特性的计算方法开展数值仿真分析;
e、建立低频平面声呐阵阻抗特性近场声全息理论测试模型
通过对基于近场声全息的辐射阻抗校准方法仿真理论模型的数值与解析仿真分析,以及对不同声呐阵布阵结构情况下的校准误差仿真分析,针对大型平面声呐阵,建立相应的低频平面声呐阵阻抗特性近场声全息理论模型;
f、设计平面声呐阵阻抗特性近场声全息法校准方案。
进一步地,所述步骤d中的仿真分析包括以下方面,
d1、声场测量几何参数对校准结果的影响;
d2、不同声呐阵布阵结构情况下的校准误差分析;
d3、环境干扰、电噪声等因素引起的随机测量误差和系统误差引起的校准误差分析;
d4、幅值与相位的不一致性对全息相位重构精度的影响分析。
进一步地,所述步骤f需构建平面声呐阵阻抗特性近场声全息法校准测量系统,具体是根据所建立的低频平面声呐阵阻抗特性近场声全息理论测试模型,建立近场声全息多通道测量系统,测量接收系统由测量水听器阵、信号衰减器、采集系统及计算机构成。
本发明的有益效果是:本发明根据水声监听船主动声呐设备声特性校准的需求,有针对性地开展平面声呐阵阻抗特性近场声全息法校准基本原理研究、近场声全息法校准方案设计、低频声呐阵阻抗特性校准软件开发以及近场声全息法校准测量试验的相关研究工作,以提升低频发射换能器校准系统的校准能力,解决低频声呐阵阻抗特性的评价问题,可利于对水声声呐为低频的船舶进行监听,提高船舶低频声呐监听结果的准确性。
附图说明:
图1为本发明的研究技术路线框图;
图2为近场声全息多通道测量系统框图;
图3为水听器阵扫描示意图;
图4为系统技术架构图。
具体实施方式:
实施例一
参照图1,本实施例中的平面声呐阵阻抗特性近场声全息法校准测量方法,包括以下步骤:
a、平面发射阵阻抗特性计算方法
发射器在辐射声波的同时会受到声场的作用,表现为贴近振动面的介质对其施加的作用力,该作用力与由振动面施加给介质的力在数值上相等而方向相反,取决于贴近该面的介质中的声压,根据发射面上声压分布沿振动面的积分可以得到介质对振动面的作用力,而声场中的声压和该处介质的质点振速成正比,即:
P=u·Z0 (1)
可以认为,贴近发射面的介质质点的法向振速和发射面的法向振速完全相同,于是发射面上的声压分布可以表示为:
Pa=ua·Za (2)
其中ua:发射面上法向振速;Za:贴近发射面的介质波阻抗。
声辐射器的辐射阻抗决定于声压对辐射器的辐射面的作用。因为基阵中的每个辐射器处于整个阵的辐射声场中,它面上的声压是所有阵元也包括自身在内的辐射场声压的叠加。因此基阵中每个阵元的总辐射阻抗包括自身辐射场的作用及其他阵元辐射场的作用。当其他阵元不振动时,单辐射器辐射时的辐射阻抗称之为自辐射阻抗。辐射器在声场中受到其他辐射器声场的作用力,为了保持自身声源强度一定,其机械系统的推动力必然有所变化。既然振幅值保持不变,而辐射面受声场的作用力有变化,这表明辐射器的总辐射阻抗发生变化。由于其他声源作用而引起变化的部分,就是总辐射阻抗中除去自辐射阻抗外的部分就是互辐射阻抗引起的。因此基阵中阵元声复合作用可以用辐射场的相互作用力表示,也可以用互辐射阻抗表示。
假设声场中有n个声源。各声源的振速分别为u1,u2,u3…un。在这种振速分布情况下,1号源辐射面上场的总作用力为F1,显然为1号源自身辐射场声压作用在自身辐射面上产生的力与其他源辐射场声压作用在1号源辐射面上产生的力之和,即
其中:f11——1号辐射器自身辐射场声压作用在自身辐射面产生的力;
fin——n号辐射器辐射场声压作用在i号辐射器辐射面产生的力。
i号辐射器总的辐射阻抗为
由于声源产生的场和声源强度成正比,在均匀振幅分布情况下,辐射场的声压和辐射面的振速成正比,所以可以把fil写成
fil=Zilul(5)
则可以得到
其中Zil表示l单位振速声场作用于i号源面上的力,称之为l号源对i号源的互辐射阻抗,简称互阻抗。
当辐射器以同样的振速和相位振动,那么由n个阵元组成的基阵的辐射阻抗就是
由以上原理可知,只要测量得到各阵元表面的复振速、复声压分布以及该阵元表面位置处由其他阵元产生的辐射声压,即可求得各阵元的自/互辐射阻抗;
b、基于近场声全息的辐射阻抗校准原理研究
近场声全息技术利用在声源近场区域内获得的复声压、复振速或声强等声学量,结合表面及声场中的亥姆霍兹-基尔霍夫(Helmholtz-Kirchhoff)积分公式,实现对声源表面或其近场处其他区域内的声振特性重构,以及远场处的辐射声场、辐射声功率等进行预报。
平面全息变换技术由于算法及操作过程简单,尽管并不完全适用于表面形状复杂的声源识别问题,却仍然是全息变换技术中最具有代表性的,发展也比较成熟。在近场平面声全息变换技术中,如果得到的是全息面上的复声压数据,则采用狄里克利(Dirichlet)边界条件处理亥姆霍兹-基尔霍夫积分公式,如果得到的是法向质点振速,则采用纽曼(Neumann)边界条件处理积分公式,这两种边界条件分别对应于瑞利第一积分方程和瑞利第二积分方程,无穷大平面上的格林函数分别为:
在Dirichlet边界条件下,可以表示为
根据二维傅里叶变换,可以得到S面上声压信号的角谱与声场中任意平面上辐射复声压信号的角谱之间的关系:
其中
对式(12)进行傅里叶逆变换,即可得到声场中任意一个平面上的场点声压:
在Neumann边界条件下,可以表示为
根据欧拉方程可得,
将其代入式(15)即可得到如下表达形式:
对上式两边同时进行二维傅里叶变换,可以得到波数域下的角谱关系式:
其中
因此,利用源平面上的质点振速可以得到声场中任意一个平面上的场点声压:
当源平面上声压或质点振速为未知量时,假设全息面与源面平行,位于zH处,且zH>zS,则可以得到源面上声压角谱表达式:
得到源平面上声压角谱后,可以得到空间中其他平面上(z>zS)的声压角谱分布:
而对于法向振速而言,可以利用式(17)得到重构方程:
经过二维傅里叶逆变换后,则可以通过全息面上的复声压分布得到源面及空间声场中复振速分布,如下所示:
利用上式得到源面上的复振速分布后,按照各阵元所在位置对复振速进行分块处理,得到各阵元源面复振速,计算其源面复声压及其他阵元表面位置处由该阵元产生的辐射声压,从而计算各阵元的自/互辐射阻抗;
c、基于BAHIM法近场声全息变换原理和算法研究
基于声强测量的宽带声全息技术(Broadband acoustic holography fromintensity measurement,简称BAHIM),它通过在全息面上的声强测量来重构全息面相位,进而得到全息面复声压。与传统方法相比,该方法不需要了解激励源及与之有关的参考信号,因而能用于宽带声源,基本重构原理介绍如下:
声波在理想介质中传播时,能量随着介质质点振动状态沿波的传播方向传播,从一个区域流向另一个区域,因此引入了介质中能流的概念。定义单位时间通过与能量传播方向垂直的单位面积上的声能为声能流密度,以表示。在时域中,成为瞬时声强,定义为:
瞬时声强可以表示为:
对于三维声场,任意点处的声强可以表示为:
把上式可以推广到包含实部和虚部的复数声强,即
对空间声场中非单频的限带,通过傅立叶变换将其分解成一系列简谐信号的叠加。在任意频率ω处,空间声场任意点(x,y,z)处时均声强的频域表达式:
可以得到:
由于有功声强中包含相位梯度信息,如果知道空间声场的有功声强和声压幅值的分布,则可以求出相位梯度:
其中,k=ω/c为声波波数,c为介质中的波速。βx(kx,ky)、βy(kx,ky)分别为点(x,y,z)处x方向和y方向声强与声压幅值平方的比值的二维FFT变换,其表达式如下
通过以上的数学推导,得到了全息面上声压相位在波数域的分布,再对(32)式做二维FFT反变换就可以得到空间域全息面上相位的分布,进而得到全息面复声压的分布,便可应用正常的NAH技术;
d、基于近场声全息的辐射阻抗校准方法数值仿真
针对基于近场声全息理论进行声呐阵辐射阻抗校准的测量模型进行数值仿真计算,分析测量方法与重构算法的误差,同时针对任意形状下近场测量数据推算远场特性的计算方法开展数值仿真分析;
e、建立低频平面声呐阵阻抗特性近场声全息理论测试模型
通过对基于近场声全息的辐射阻抗校准方法仿真理论模型的数值与解析仿真分析,以及对不同声呐阵布阵结构情况下的校准误差仿真分析,针对大型平面声呐阵,建立相应的低频平面声呐阵阻抗特性近场声全息理论模型;
f、平面声呐阵阻抗特性近场声全息法校准方案设计。
一、本实施例中还需进行低频声呐阵阻抗特性校准软件开发:
(一)首先进行低频声呐阵阻抗特性校准软件架构设计
基于平面声呐阵阻抗特性近场声全息法校准原理模型,进行软件架构和人机交互界面设计,软件需要完成以下功能:校准过程中的各种参数设置;控制仪器仪表设备的正常工作;完成各种算法运算;校准结果的打印输出以及校准结果的保存等。最终开发的低频声呐阵阻抗特性校准软件包含测量系统一致性校准功能模块、声呐阵阻抗特性校准功能模块、试验仪器设备控制模块,文件输出模块等。
本软件系统,开发语言采用C++语言,软件UI采用VS平台进行开发,方便本地操作以及调用操作系统本地资源。展示校准结果文件渲染采用的是开源图形图像处理库OpenGL,大部分图像生成和操作采用封装OpenGL库的OSG框架实现。在进行仪器仪表控制操作、参数配置过程生成的中间数据、以及测量试验数据均采用轻量化SQLite数据库进行存储,SQLite数据库采用文件数据库,方便安装、操作灵活。软件系统技术架构如图4所示。
(二)、测量系统一致性校准功能模块开发
结合近场声全息法校准测量系统的一致性校准基本原理,采用C++软件进行校准功能模块开发,并编译形成可执行文件,通过界面调用该可执行文件,基于一致性校准测量参数,实现近场声全息法测量系统各通道间的一致性校准,并存储一致性校准结果。
(三)、声呐阵阻抗特性校准功能模块开发
结合平面声呐阵阻抗特性近场声全息法校准基本原理,采用C++软件进行声呐阵阻抗特性校准功能模块开发,并编译形成可执行文件,通过界面调用该可执行文件,读取各通道间一致性校准结果,对各通道采集数据进行预处理,进行预处理结果存储。利用预处理数据进行声呐阵阻抗特性计算分析,并存储结果。
(四)、试验仪器设备控制模块开发
利用VS平台,编写各试验仪器的控制程序,进行参数设置与工作状态控制,实现平面声呐阵阻抗特性近场声全息法一体化校准功能。
(五)、计算结果展示模块开发
计算结果展示模块将利用OpenGL三维图形软件与Visual Studio开发工具相结合,进行VS界面下的OpenGL编程,用以展示结果文件,以及生成结果文件渲染,并利用鼠标或键盘消息实现信息交互。通过调用测量系统一致性校准结果、各测量通道采集数据预处理结果以及声呐阵阻抗特性校准结果等数据文件,实现计算结果的实时显示。
二、本实施例中还需进行平面声呐阵阻抗特性近场声全息法校准测量试验
(一)、系统集成和联调
对所建立的多通道测量系统及测量系统的一致性校准系统在进行系统联调试验,包括发射系统及接收系统;
(二)、多通道测量系统的一致性校准与修正
采取易地式耦合腔互易校准方法,以B&K 8103水听器作为标准水听器,对多通道测量系统通道间的幅度和相位的一致性进行校准,并对所有水听器进行配对,组成16对双水听器探头(或80对双水听器探头);
当所有水听器组装为声全息测量基阵后,针对整体的测量环境,利用非线性最小二乘估计法对全息扫描阵的相位一致性进行校准;
(三)、平面声呐阵阻抗特性近场声全息法校准试验基本方案
通过对平面声呐阵阻抗特性近场声全息法校准的基本原理研究,综合分析解析仿真结果与数值仿真结果,设计实际声全息校准测量试验;
(1)试验场地:开阔水域(自由场);
(2)测试样品:平面发射阵(九元阵);
(3)检验频率:100Hz~630Hz;
(4)样品尺寸:长3m,宽1m;
(5)声场布放:合理放置水听器阵,尽量保证声场测试结果的稳定性;
(6)声场测量参数:
(6.1)测量面尺寸
为确保全息变换场的精确,保证测量数据获得必要多声源辐射能量和声源信息,全息面必须足够大。若已知声源面尺度与波长比和测量距离,即可确定测量面的尺寸。若以测量孔径边缘处声压幅值有30dB衰减为准则,对于L/λ=1的情况,测量距离为1/4λ时,测量面边长Lx应为声源面边长L的3倍;测量距离为1/2波长时,测量面边长Lx应为声源面边长L的3.5倍。对于L/λ≥3的情况,测量距离为1/4波长时,测量面为声源面的2倍即可。即针对低频小尺寸声源(L/λ较小),测量面和声源面为共形面最佳,否则需要更大的测量面和较近的测量距离。当L/λ较大时,一般测量面为声源面的2倍即可;
针对本项目中平面发射阵尺寸及测试频段,测量面选择为声源面的2倍左右,约6m*2m;
(6.2)测量间距
为精确测量声场分布,要求全息面上测点足够密。鉴于表面声场波幅驻波形式分布,类似频率-时间域中采样定理,空间测量点间距应为1/(7~10)λ较为合理;
针对本项目中测试频段,测量点间距选择为20cm~30cm;
(6.3)双水听器探头之间的距离
鉴于实用测量水听器的尺寸过小其灵敏度和可靠性都较差,尺寸较大,相互间散射又会影响到单水听器的自由声场。因此,对于水声宽频带声场的测量以及采用双水听器阵对声场作扫描测量,理想情况下,Δr/λ=1/10;对于工程上可采用Δr/λ=1/7;
针对本项目中测试频段,双水听器探头间距选择为20cm~30cm;
(6.4)发射阵与水听器阵距离
高波数波的波幅沿声源面法线方向减弱很快,因此,要求全息面到源面间距尽可能小。对于高频大声源,其间距取1/3λ较为合适;对于低频小声源,则应当取更小距离与波长比为宜;
针对本项目中测试频段,水听器阵与发射阵之间距离选择为40cm~60cm;
(6.5)声场扫描方式:采用程控方式;
(7)测量信号:足够长(5-10周期),多次采集(5-20)并进行多次平均。
(四)、试验前准备工作
(1)测试平面阵准备
按要求组装测试平面阵,清洗平面阵表面,并充分浸泡,使表面充分浸润。
(2)测试水听器阵准备
按照双水听器配对结果,安装声全息测试水听器阵,并清洗水听器表面,将水听器阵与机械扫描装置连接。
(3)测试仪表设备准备
(3.1)检查各测量用仪表是否工作正常,仪表间线路连接是否正确。
(3.2)测试前,检查发射器、水听器是否能够正常工作,与仪表连接是否正确。
(3.3)检查仪表设置是否正确,如功率放大器是否选择合适的阻抗匹配档位,多通道滤波系统是否设置了合适的高通和低通,放大倍数是否核实,以防仪器过载。
(3.4)测试前,打开仪表设备,预热半个小时以上。
(3.5)检查扫描装置行走是否正常,是否能够达到测试精度要求。
(五)、平面声呐阵阻抗特性近场声全息测量方法
(1)将平面发射阵竖直置于开阔水域中,入水深度约为开阔水域深度的1/2;
(2)按照配好对的声强探头组装测量接收直线阵,每两对水听器探头间距要满足声强测量参数;
(3)将水听器阵安装在扫描装置上;利用机械扫描系统,放置在距离平面发射阵不小于1/6波长的位置处,如图3所示:
(4)将平面发射阵及水听器阵静止放置12小时后,程控机械扫描系统实现进行水平方向扫描,扫描完毕后将水听器阵向下移动,使其测点间距满足全息测试参数及声强测试参数,获得两个全息面上的声压数据;
(5)利用全息数据反演平面发射阵各阵元的质点振速,利用边界元法得到各阵元表面的声压,以及由该阵元产生的其他阵元表面的声压,分析各振元间的互辐射阻抗;
(6)利用双水听器探头测量数据计算双水听器对中心位置处的复振速,将其近似划分为各阵元表面的质点振速分布,进而计算平面发射阵各阵元的自/互辐射阻抗,以及远场辐射声场特性计算。
(六)、测量不确定度分析
在每个工况下,进行多次测量试验。利用多次重复测量的数据,对系统进行不确定度评定。系统的不确定度分为A类不确定度和B类不确定度。
A类不确定度:
选用极差法对A类不确定度进行计算评估:
其中,uA:A类不确定度分量;
s(x):试验标准偏差;
xmax:多次测量试验数据的最大值;
xmin:多次测量试验数据的最小值。
C:极差系数,详见下表:
极差系数表
测量次数 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 |
C | 1.13 | 1.69 | 2.06 | 2.33 | 2.53 | 2.70 | 2.85 | 2.97 |
B类不确定度:uB通过对实际测量环境、测量仪器仪表、水听器灵敏度偏差、测量方法等的综合分析获得。
扩展不确定度:U=k×uC(本试验中取k=2)。
实施例二
本实施例在实施例一的基础上,步骤d中的仿真分析包括以下方面:
d1、声场测量几何参数对校准结果的影响
(1)加窗函数对校准结果的影响;
(2)全息面半径:0.5λ~20λ;
(3)测点间距:0.1λ~0.5λ;
(4)全息面与待校样品的间距:0.1λ~0.25λ;
(5)双水听器之间间距:0.05λ~0.5λ;
d2、不同声呐阵布阵结构情况下的校准误差分析
(1)八元阵(2×4):2m×1.3m;
(2)九元阵(3×3):3m×1m;
(3)九元阵(9×1):1m×3m;
d3、环境干扰、电噪声等因素引起的随机测量误差和系统误差引起的校准误差分析;
d4、幅值与相位的不一致性对全息相位重构精度的影响分析。
实施例三
本实施例在实施例二的基础上,在所述F步骤中需构建平面声呐阵阻抗特性近场声全息法校准测量系统,具体是根据所建立的低频平面声呐阵阻抗特性近场声全息理论测试模型,建立近场声全息多通道测量系统,测量接收系统由测量水听器阵、信号衰减器、采集系统及计算机构成,如图2所示,
(1)、水听器阵
由32个水听器组成测量接收直线阵;每两个水听器构成一组声强探头,共16对声强探头(或者采用160个水听器组成80对声强探头,每对声强探头间距为20cm);拟采用32个B&K8103水听器(或160个715所生产水听器);接收灵敏度高于-210dB;工作压力:0~4MPa;直径尺寸小于15mm;
在100Hz~5kHz频率范围内双水听器幅度一致性偏差:≤±1dB;相位一致性偏差:≤±1°;测量不确定度:U≤3dB。
(2)、声全息测量基阵支架:自研,形式:直线阵(32个水听器组成16对双水听器探头)或方阵(160个组成1680双水听器探头);水听器所在位置散射强度小于-40dB;
(3)、信号衰减器
测试频段:100Hz~5kHz;
32个通道,通道间幅度一致性偏差:≤±1dB;相位一致性偏差:≤±1°;
信号衰减范围:20dB~50dB,可调节。
(4)、测量放大、滤波器:拟采用PF28000,不小于32通道,同时具备滤波、放大功能;前端输入噪声:<5微伏;放大量:>60dB;动态范围:>60dB;
(5)、采集系统:拟采用B&K 3660D采集器,Pulse 14.1以上版本;
分析信号频率范围:0~25.6kHz;
输入信号量程:7mv~25v;
分辨率:24bit;
输入噪声:比所选用量程对应的最大峰值电压低80分贝;
幅度一致性:100Hz~2kHz,典型值0.01dB,最大0.2dB;
相位一致性:100Hz~2kHz,<1.5°。
(6)、机械扫描控制系统
移动自由度:X、Y、Z三个方向;
定位精度:高于0.1mm。
Claims (1)
1.一种平面声呐阵阻抗特性近场声全息法校准测量方法,其特征在于,包括以下步骤:
a、平面发射阵阻抗特性计算方法
发射器在辐射声波的同时会受到声场的作用,表现为贴近振动面的介质对其施加的作用力,该作用力与由振动面施加给介质的力在数值上相等而方向相反,取决于贴近该面的介质中的声压,根据发射面上声压分布沿振动面的积分可以得到介质对振动面的作用力,而声场中的声压和该处介质的质点振速成正比,即:
P=u·Z0 (1)
由以上获得,贴近发射面的介质质点的法向振速和发射面的法向振速完全相同,于是发射面上的声压分布表示为:
Pa=ua·Za (2)
其中ua:发射面上法向振速;Za:贴近发射面的介质波阻抗,
假设声场中有n个声源,各声源的振速分别为u1,u2,u3…un,在这种振速分布情况下,1号源辐射面上场的总作用力为F1,显然为1号源自身辐射场声压作用在自身辐射面上产生的力与其他源辐射场声压作用在1号源辐射面上产生的力之和,即
其中:f11——1号辐射器自身辐射场声压作用在自身辐射面产生的力;
fin——n号辐射器辐射场声压作用在i号辐射器辐射面产生的力,
i号辐射器总的辐射阻抗为
由于声源产生的场和声源强度成正比,在均匀振幅分布情况下,辐射场的声压和辐射面的振速成正比,所以可以把fil写成
fil=Zilul (5)
则可以得到
其中Zil表示l单位振速声场作用于i号源面上的力,称之为l号源对i号源的互辐射阻抗,简称互阻抗,
当辐射器以同样的振速和相位振动,那么由n个阵元组成的基阵的辐射阻抗就是
由以上可知,只要测量得到各阵元表面的复振速、复声压分布以及该阵元表面位置处由其他阵元产生的辐射声压,即可求得各阵元的自/互辐射阻抗;
b、基于近场声全息的辐射阻抗校准方法
无穷大平面上的格林函数分别为:
在狄里克利边界条件下,可以表示为
根据二维傅里叶变换,可以得到S面上声压信号的角谱与声场中任意平面上辐射复声压信号的角谱之间的关系:
其中
对式(12)进行傅里叶逆变换,即可得到声场中任意一个平面上的场点声压:
在纽曼边界条件下,可以表示为
根据欧拉方程可得,
将其代入式(15)即可得到如下表达形式:
对上式两边同时进行二维傅里叶变换,可以得到波数域下的角谱关系式:
其中
因此,利用源平面上的质点振速可以得到声场中任意一个平面上的场点声压:
当源平面上声压或质点振速为未知量时,假设全息面与源面平行,位于zH处,且zH>zS,则可以得到源面上声压角谱表达式:
得到源平面上声压角谱后,可以得到空间中其他平面上(z>zS)的声压角谱分布:
而对于法向振速而言,可以利用式(17)得到重构方程:
经过二维傅里叶逆变换后,则可以通过全息面上的复声压分布得到源面及空间声场中复振速分布,如下所示:
利用上式得到源面上的复振速分布后,按照各阵元所在位置对复振速进行分块处理,得到各阵元源面复振速,计算其源面复声压及其他阵元表面位置处由该阵元产生的辐射声压,从而计算各阵元的自/互辐射阻抗;
c、基于BAHIM法近场声全息变换方法和算法
基于声强测量的宽带声全息技术,简称BAHIM,它通过在全息面上的声强测量来重构全息面相位,进而得到全息面复声压,基本重构方式如下:
声波在理想介质中传播时,能量随着介质质点振动状态沿波的传播方向传播,从一个区域流向另一个区域,引入了介质中能流的概念,定义单位时间通过与能量传播方向垂直的单位面积上的声能为声能流密度,以表示,在时域中,成为瞬时声强,定义为:
瞬时声强可以表示为:
对于三维声场,任意点处的声强可以表示为:
把上式可以推广到包含实部和虚部的复数声强,即
对空间声场中非单频的限带,通过傅立叶变换将其分解成一系列简谐信号的叠加,在任意频率ω处,空间声场任意点(x,y,z)处时均声强的频域表达式:
可以得到:
由于有功声强中包含相位梯度信息,如果知道空间声场的有功声强和声压幅值的分布,则可以求出相位梯度:
对(31)式两边作二维FFT变换,由空间域变到波数域后,全息面上声压相位的波数域形式为
其中,k=ω/c为声波波数,c为介质中的波速,βx(kx,ky)、βy(kx,ky)分别为点(x,y,z)处x方向和y方向声强与声压幅值平方的比值的二维FFT变换,其表达式如下
通过以上的数学推导,得到了全息面上声压相位在波数域的分布,再对(32)式做二维FFT反变换就可以得到空间域全息面上相位的分布,进而得到全息面复声压的分布,便可应用正常的NAH技术;
d、基于近场声全息的辐射阻抗校准方法数值仿真
针对基于近场声全息理论进行声呐阵辐射阻抗校准的测量模型进行数值仿真计算,分析测量方法与重构算法的误差,同时针对任意形状下近场测量数据推算远场特性的计算方法开展数值仿真分析;所述仿真分析包括以下方面;
d1、声场测量几何参数对校准结果的影响;
d2、不同声呐阵布阵结构情况下的校准误差分析;
d3、环境干扰、电噪声等因素引起的随机测量误差和系统误差引起的校准误差分析;
d4、幅值与相位的不一致性对全息相位重构精度的影响分析;
e、建立低频平面声呐阵阻抗特性近场声全息理论测试模型
通过对基于近场声全息的辐射阻抗校准方法仿真理论模型的数值与解析仿真分析,以及对不同声呐阵布阵结构情况下的校准误差仿真分析,针对大型平面声呐阵,建立相应的低频平面声呐阵阻抗特性近场声全息理论模型;
f、设计平面声呐阵阻抗特性近场声全息法校准方案
需构建平面声呐阵阻抗特性近场声全息法校准测量系统,具体是根据所建立的低频平面声呐阵阻抗特性近场声全息理论测试模型,建立近场声全息多通道测量系统,测量接收系统由测量水听器阵、信号衰减器、采集系统及计算机构成。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811504271.5A CN109444861B (zh) | 2018-12-10 | 2018-12-10 | 一种平面声呐阵阻抗特性近场声全息法校准测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811504271.5A CN109444861B (zh) | 2018-12-10 | 2018-12-10 | 一种平面声呐阵阻抗特性近场声全息法校准测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109444861A CN109444861A (zh) | 2019-03-08 |
CN109444861B true CN109444861B (zh) | 2022-05-24 |
Family
ID=65558015
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811504271.5A Expired - Fee Related CN109444861B (zh) | 2018-12-10 | 2018-12-10 | 一种平面声呐阵阻抗特性近场声全息法校准测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109444861B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111505568B (zh) * | 2020-04-14 | 2022-06-21 | 哈尔滨工程大学 | 基于差分优化的四面体阵目标方位估计方法 |
CN111812425B (zh) * | 2020-05-28 | 2023-03-24 | 国网河北省电力有限公司电力科学研究院 | 一种基于近场声全息技术的高压并联电抗器故障检测方法 |
CN112577592B (zh) * | 2020-11-27 | 2022-10-28 | 哈尔滨工程大学 | 基于空间傅里叶变换的有限空间平面近场声全息测量方法 |
CN114136427B (zh) * | 2021-12-02 | 2023-05-26 | 哈尔滨工程大学 | 一种可安装于结构表面的水下法向声能流测量装置 |
CN114543979B (zh) * | 2022-02-17 | 2024-05-03 | 浙江工业大学 | 有界空间中基于近场声全息的声源直接辐射远场声学量的预测方法 |
CN116952355B (zh) * | 2023-07-24 | 2024-05-14 | 中国人民解放军海军工程大学 | 一种浅海环境近场辐射噪声测量系统及终端 |
CN116952356B (zh) * | 2023-07-24 | 2024-08-06 | 中国人民解放军海军工程大学 | 基于浅海环境水下声全息技术的近场辐射噪声测量方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108051800A (zh) * | 2017-12-13 | 2018-05-18 | 贵州航天计量测试技术研究所 | 基于球面近场声全息重构无功声强的室内噪声源定位方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008064639A (ja) * | 2006-09-08 | 2008-03-21 | Ono Sokki Co Ltd | 放射パワー演算装置、音圧演算装置、粒子速度演算装置、放射パワー演算プログラム、音圧演算プログラム、および、粒子速度演算プログラム |
JP5453009B2 (ja) * | 2009-07-30 | 2014-03-26 | 日東紡音響エンジニアリング株式会社 | 音響測定方法及びシステム |
CN103543438A (zh) * | 2013-07-08 | 2014-01-29 | 哈尔滨工程大学 | 配置参考声源的水下运动声源柱面矢量聚焦定位方法 |
CN103743469B (zh) * | 2013-12-11 | 2015-12-02 | 西安交通大学 | 压电水声换能器声辐射模态测量方法及系统 |
CN106813767B (zh) * | 2017-01-03 | 2019-07-16 | 湖北文理学院 | 一种usp探头的灵敏度测量方法 |
CN107478325B (zh) * | 2017-09-07 | 2019-07-26 | 安徽大学 | 一种非稳态平面声源的自由场还原方法 |
CN108318123A (zh) * | 2017-12-14 | 2018-07-24 | 中国船舶重工集团公司第七〇五研究所 | 一种水下声基阵互辐射阻抗的近场声层析测试方法 |
CN108519431B (zh) * | 2018-03-21 | 2020-09-11 | 哈尔滨工程大学 | 一种测量声学覆盖层高频法向反射系数的装置及方法 |
-
2018
- 2018-12-10 CN CN201811504271.5A patent/CN109444861B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108051800A (zh) * | 2017-12-13 | 2018-05-18 | 贵州航天计量测试技术研究所 | 基于球面近场声全息重构无功声强的室内噪声源定位方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109444861A (zh) | 2019-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109444861B (zh) | 一种平面声呐阵阻抗特性近场声全息法校准测量方法 | |
CN109375198B (zh) | 一种低频声呐阵阻抗特性校准方法 | |
Potthast | A study on orthogonality sampling | |
CN106680376B (zh) | 一种基于三维声强阵列的隔声测量系统与测量方法 | |
Holstein et al. | Acoustic tomography on the basis of travel-time measurement | |
Narayanan et al. | High subsonic jet experiments: turbulence and noise generation studies | |
Merino-Martínez et al. | Assessment of the accuracy of microphone array methods for aeroacoustic measurements | |
CN107478325B (zh) | 一种非稳态平面声源的自由场还原方法 | |
KR20070072518A (ko) | 소음원의 원거리-장 분석 | |
CN110186546A (zh) | 基于粉红噪声的水听器灵敏度自由场宽带校准方法 | |
Ziemann et al. | Acoustic tomography in the atmospheric surface layer | |
CN110530510B (zh) | 一种利用线性声阵列波束形成的声源辐射声功率测量方法 | |
Jones et al. | Experimental validation of determining sound power using acoustic radiation modes and a laser vibrometer | |
CN109839627B (zh) | 一种水池条件下声呐探测距离量值计量校准系统及方法 | |
CN103048641B (zh) | 具有联合约束优化形式的矢量阵稳健聚焦处理方法 | |
CN110082431A (zh) | 一种用于材料表面声阻抗测量的方法及装置 | |
CN206563731U (zh) | 一种基于三维声强阵列的隔声测量系统 | |
Gombots | Acoustic source localization at low frequencies using microphone arrays | |
Zhang et al. | Localization of cyclostationary acoustic sources via cyclostationary beamforming and its high spatial resolution implementation | |
CN113639934A (zh) | 气体泄漏位置三维定位方法、定位系统 | |
Wang et al. | A fast irregular microphone array design method based on acoustic beamforming | |
Li et al. | Research of the improved mapping of acoustic correlated sources method | |
Batko et al. | Limitations of decibel algebra in the study of environmental acoustic hazards | |
Knežević et al. | Automated Sound Intensity Measuremant With Robot And Intensity Probe | |
Hald | Beamforming and wavenumber processing |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220524 |