CN109238440B - 基于最优空时聚焦技术的声学覆盖层插入损失测量方法 - Google Patents

基于最优空时聚焦技术的声学覆盖层插入损失测量方法 Download PDF

Info

Publication number
CN109238440B
CN109238440B CN201810788869.5A CN201810788869A CN109238440B CN 109238440 B CN109238440 B CN 109238440B CN 201810788869 A CN201810788869 A CN 201810788869A CN 109238440 B CN109238440 B CN 109238440B
Authority
CN
China
Prior art keywords
signal
focusing
transmitting
insertion loss
optimal
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
CN201810788869.5A
Other languages
English (en)
Other versions
CN109238440A (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201810788869.5A priority Critical patent/CN109238440B/zh
Publication of CN109238440A publication Critical patent/CN109238440A/zh
Application granted granted Critical
Publication of CN109238440B publication Critical patent/CN109238440B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H17/00Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

基于最优空时聚焦技术的声学覆盖层插入损失测量方法,包括以下步骤:1)生成最优空时聚焦发射信号;2)无试样情况下直达信号的采集;3)有试样情况下透射信号的采集;4)插入损失测量值的计算。本发明利用多通道发射阵及接收阵目标函数最优化方法,通过目标函数最优解独立控制发射信号,从而实现最优空时聚焦,获取空时域高分辨聚焦信号,实现直达波与多径反射信号的分离及混响抑制,适用于低频条件下声学覆盖层插入损失的测量;同时,压力罐测量实验验证了本发明在声学覆盖层插入损失测量中的有效性。

Description

基于最优空时聚焦技术的声学覆盖层插入损失测量方法
技术领域
本发明涉及一种水下声学覆盖层大样插入损失的低频测量方法。水下声学覆盖层是水声工程中使用广泛且至关重要的水下部件,为评估声学覆盖层试样的插入损失性能,实验室条件下大件样品的插入损失测量是主要的手段之一。
背景技术
声学覆盖层是指敷设在水下部件上的专用功能性声学材料与结构,是水下部件实现声隐身的一个重要手段。随着复合材料的快速发展以及声隐身技术研究的深入,声学覆盖层产品已较为齐全,性能亦逐步提高。声学覆盖层一方面可以将主动探测声波的能量大幅度的吸收,使主动声纳声波的反射减少,降低水下结构的声目标强度,另一方面可以抑制结构自身辐射噪声。为对不同性能的覆盖层性能进行评估,通常包括仿真计算和实际测量两种手段,由于材料结构和组成较为复杂,计算的精度受到限制,对其在水声环境中进行声学性能测量是必不可少的手段。
插入损失是衡量覆盖层声学性能的重要指标之一。现有的声学覆盖层插入损失测量方法包括空间傅里叶变换方法,宽带脉冲压缩法,时间反转聚焦方法、多通道空时逆滤波法等。有限空间环境下,混响及多径效应严重,无法实现平面波的精确分解,从而限制了空间傅里叶变换方法的应用;宽带脉冲压缩法可实现时域脉冲信号波形聚焦,但低频条件下,其单通道发射的特点导致指向性变差,衍射及混响干扰较严重;时间反转聚焦方法可实现空时聚焦,但由于其没有抵消信道幅度的影响,聚焦的时域信号脉宽受限,无法实现最优聚焦效果,不利于低频情况下插入损失的测量;多通道空时逆滤波技术没有用到接收阵及发射阵各阵元信道传递函数之间的相干信息,对于具有多个接收阵元的测量系统,没有实现最佳空时聚焦。
发明内容
为了克服现有技术的上述缺点,进一步提高声学覆盖层插入损失测量技术在低频段测试中的测量精度,本发明提供一种适用于低频段,可有效减少测量误差、提高测量精度的声学覆盖层大样插入损失测量方法,该方法利用多通道发射阵及接收阵目标函数最优化方法,通过目标函数最优解独立控制发射信号,从而实现最优空时聚焦,获取空时域高分辨聚焦信号,实现直达波与多径反射信号的分离及混响抑制。压力罐实验验证了本发明在声学覆盖层插入损失测量中的有效性。
本发明为解决上述技术问题所采用的技术方案如下。
一种基于最优空时聚焦技术的声学覆盖层插入损失测量方法,包括如下步骤:
1)生成最优空时聚焦发射信号:无试样情况下,各发射换能器依次发射初始导引信号,各水听器依次接收信号,从而获取各个发射-接收信道的响应函数,通过目标函数最优化方法,获得各发射通道的二次发射信号;
2)无试样情况下直达信号的采集:各换能器同步发射1)中计算得到的发射信号,则发射信号将在预设聚焦位置处实现最佳聚焦,水听器记录无测试样品时的直达信号pi
3)有试样情况下透射信号的采集:将敷设声学覆盖层的试样放入测试环境,重复2)中步骤,记录有测试样品时的透射信号pt
4)插入损失测量值的计算:利用插入损失的计算公式可得
Figure BDA0001734346420000021
其中IL表示插入损失值。
本发明的技术构思为:通过初始导引信号发射及水听器逐一接收,得到包含电路信道和水声信道信息的接收信号,进而通过目标函数最优化方法,获得换能器阵每一阵元的二次发射信号,然后各换能器同步发射,获得预设水听器处的最优聚焦信号,从而最优化地实现测量信号的空间聚焦及时域脉冲压缩,达到抑制低频混响的目的。
与现有的插入损失测量方法相比,本发明的优点是:利用各发射阵元及接收阵元之间的信道相干信息,通过目标函数最优化方法,获取各发射阵元最优聚焦发射信号,经同步发射,实现预设水听器处的最优空时聚焦,达到低频高分辨及混响抑制的目的。
附图说明
图1是本发明方法的测量系统整体示意图。
图2a和图2b分别是仿真计算得到的最优空时聚焦方法在
Figure BDA0001734346420000031
和(0m,0rad)位置处聚焦的空间能量分布图。
图3a是初始发射信号时域波形图,图3b和图3c分别是
Figure BDA0001734346420000032
和(0m,0rad)位置处的时域脉冲聚焦波形图。
图4是压力罐环境下5mm厚钢板试样插入损失理论及测量试验结果对比图。
具体实施方式
下面结合附图对本发明作进一步描述。
参照图1~图4,一种基于最优空时聚焦技术的声学覆盖层插入损失测量方法,用于有限空间下声学覆盖层大样插入损失的测量。整套测量方法的技术方案如下:
1)生成最优空时聚焦发射信号;
无试样情况下,各换能器依次发射初始导引信号,各水听器阵元依次接收信号,设发射信号频域可表示为en(ω),n=1,2,…,N,其中ω表示信号角频率,n为发射换能器编号,N为换能器个数;接收信号表示为y(ω,rn;rm),m=1,2,…,M,其中m为水听器编号,M为水听器个数,rn、rm分别表示第n个发射换能器和第m个水听器的位置;
根据公式
Figure BDA0001734346420000033
估计每一发射通道与接收水听器振元之间的信道响应函数G(ω,rn;rm);
得到信道响应后,为求解最优化的二次发射信号,建立如下目标函数
J(s1(ω),s2(ω),…,sN(ω))=∫W|Φ(ω,r)-exp(-iωT)δ(r-rc)|2dr (3)
其中,W为水听器的位置集合,sn(ω),n=1,2,…,N,为第n个发射通道产生的二次发射信号,对应的时域信号为sn(t),其中t表示时间序列,n=1,2,…,N;设sn(t)在t=0时刻发射,经信道传递后在预设聚焦位置rc处和T时刻聚焦,要实现最优聚焦,预设位置处的聚焦信号在频域上的理论值应表示为exp(-iωT)δ(r-rc),即空时均为狄拉克函数;目标函数J(s1(ω),s2(ω),…,sN(ω))中,Φ(ω,rm),m=1,2,…,M,为第m个水听器处的接收信号,可表示为
Figure BDA0001734346420000041
其中,xn(ω,rm),m=1,2,…,M为由第n个发射通道发射,并在第m个水听器位置处接收的信号,表示为
xn(ω,rm)=sn(ω)·G(ω,rn;rm) (5)
通过最小化目标函数,最优二次发射信号可通过下式得到
Γ(ω)s(ω)=exp(-iωT)g*(ω) (6)
其中g*(ω)表示函数g(ω)的共轭,Γ矩阵的组成元素为
Γnm(ω)=∫WG(ω,rn;r)G*(ω,rm;r)dr (7)
向量g(ω)的组成元素为
gn(ω)=G(ω,rn;rc),n=1,2,…,N (8)
即Γ矩阵为各信道响应函数的互相关表示,g(ω)为各发射阵元到预设聚焦位置之间响应函数组成的向量。则最优二次发射信号可表示为
Figure BDA0001734346420000042
其中
Figure BDA0001734346420000043
表示矩阵Γ的伪逆矩阵,对该信号做逆傅里叶变换得到时域最优空时聚焦的发射信号s(t)。
2)无试样时直达信号的采集;
如图1所示,通过计算机辅助处理系统控制多通道独立控制信号发生器,经多通道功率放大器,利用发射换能器阵同步发射由步骤1)生成的最优二次发射信号,同时,计算机辅助处理系统控制聚焦位置处的水听器记录没有测试样品时的直达信号pi
3)有试样时透射信号的采集;
将敷设声学覆盖层的试样放入测试环境,重复2)中步骤,此时记录有测试样品时的透射信号pt
4)插入损失测量值的计算;
通过以上三个步骤可获得计算插入损失所需的测量数据,利用2)和3)采集的入射信号pi及透射信号pt,由公式(1)计算可得插入损失测量值。
实例说明:为验证本发明中最优聚焦技术的空时聚焦效果,开展了仿真分析。设刚性圆柱形水罐直径为4m,轴向长度为半无限长,水体声速为cw=1480m/s,水体密度为ρw=1g/cm3。3个换能器分布在直径为1.5m的圆上,圆阵中心位于圆柱的中轴线上,接收阵为五元十字阵,相邻阵元间距为1m,阵中心同样位于圆柱中轴线上,且与发射阵轴向距离为4.5m。初始发射信号频率为1kHz~1.5kHz。图2a和图2b显示的是利用最优空时聚焦技术进行二次发射时,预设聚焦位置分别在
Figure BDA0001734346420000051
和(0m,0rad)时的空间能量分布图;图3a为初始发射信号波形图,图3b和图3c分别是
Figure BDA0001734346420000052
和(0m,0rad)位置的时域脉冲聚焦波形图。由图可知,最优空时聚焦技术实现了信号的空间能量聚焦和时域脉冲聚焦。
为验证最优空时聚焦技术在声学覆盖层插入损失测量中的有效性,开展了压力水罐条件下的试验验证。实验中三元发射阵与测试样品相距4.5m,3个圆形发射换能器阵元均匀分布在半径为1.5m的圆上。测试试样为钢板,其几何尺寸为1.1m×1.0m×5mm,密度为7.84×103kg/m3,声速5470m/s。该试验开展了0.5kHz~5kHz频率下插入损失测量试验的数据采集和处理工作。图4显示的是钢板试样插入损失理论及试验结果对比图,由图可知,在多径现象非常严重的圆柱腔体波导环境下,本发明的测量结果在低频段基本与理论计算值相一致。在频率大于2kHz时,测量误差增大,产生此种现象的可能原因为换能器阵频率覆盖范围较窄,在高频段性能下降较严重。
本说明书实施例所述的内容仅仅是对发明构思的实现形式的列举,本发明的保护范围不应当被视为仅限于实施例所陈述的具体形式,本发明的保护范围也及于本领域技术人员根据本发明构思所能够想到的等同技术手段。

Claims (1)

1.基于最优空时聚焦技术的声学覆盖层插入损失测量方法,包括如下步骤:
1)生成最优空时聚焦发射信号:无试样情况下,各发射换能器依次发射初始导引信号,各水听器依次接收信号,从而获取各个发射-接收信道的响应函数,通过目标函数最优化方法,获得各发射通道的二次发射信号;生成最优空时聚焦发射信号的过程如下:
无试样情况下,各换能器依次发射初始导引信号,各水听器阵元依次接收信号,设发射信号频域可表示为en(ω),n=1,2,…,N,其中ω表示信号角频率,n为发射换能器编号,N为换能器个数;接收信号表示为y(ω,rn;rm),m=1,2,…,M,其中m为水听器编号,M为水听器个数,rn、rm分别表示第n个发射换能器和第m个水听器的位置;
根据公式
Figure FDA0002401159780000011
估计每一发射通道与接收水听器振元之间的信道响应函数G(ω,rn;rm);
得到信道响应后,为求解最优化的二次发射信号,建立如下目标函数
J(s1(ω),s2(ω),…,sN(ω))=∫W|Φ(ω,r)-exp(-iωT)δ(r-rc)|2dr (2)
其中,W为水听器的位置集合,sn(ω),n=1,2,…,N,为第n个发射通道产生的二次发射信号,对应的时域信号为sn(t),其中t表示时间序列,n=1,2,…,N;设sn(t)在t=0时刻发射,经信道传递后在预设聚焦位置rc处和T时刻聚焦,要实现最优聚焦,预设位置处的聚焦信号在频域上的理论值应表示为exp(-iωTδ)r(-rc,即空时均为狄拉克函数;目标函数J(s1(ω),s2(ω),…,sN(ω))中,Φ(ω,rm),m=1,2,…,M,为第m个水听器处的接收信号,可表示为
Figure FDA0002401159780000012
其中,xn(ω,rm),m=1,2,…,M为由第n个发射通道发射,并在第m个水听器位置处接收的信号,表示为
xn(ω,rm)=sn(ω)·G(ω,rn;rm) (4)
通过最小化目标函数,最优二次发射信号可通过下式得到
Γ(ω)s(ω)=exp(-iωT)g*(ω) (5)
其中g*(ω)表示函数g(ω)的共轭,Γ矩阵的组成元素为
Γnm(ω)=∫WG(ω,rn;r)G*(ω,rm;r)dr (6)
向量g(ω)的组成元素为
gn(ω)=G(ω,rn;rc),n=1,2,…,N (7)
即Γ矩阵为各信道响应函数的互相关表示,g(ω)为各发射阵元到预设聚焦位置之间信道响应函数组成的向量;则最优二次发射信号可表示为
Figure FDA0002401159780000021
其中
Figure FDA0002401159780000022
表示矩阵Γ的伪逆矩阵,对该信号做逆傅里叶变换得到时域最优空时聚焦的发射信号s(t);
2)无试样情况下直达信号的采集:各换能器同步发射1)中计算得到的发射信号,则发射信号将在预设聚焦位置处实现最佳聚焦,水听器记录无测试样品时的直达信号pi
3)有试样情况下透射信号的采集:将敷设声学覆盖层的试样放入测试环境,重复2)中步骤,记录有测试样品时的透射信号pt
4)插入损失测量值的计算:利用插入损失的计算公式可得
Figure FDA0002401159780000023
其中IL表示插入损失值。
CN201810788869.5A 2018-07-18 2018-07-18 基于最优空时聚焦技术的声学覆盖层插入损失测量方法 Active CN109238440B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810788869.5A CN109238440B (zh) 2018-07-18 2018-07-18 基于最优空时聚焦技术的声学覆盖层插入损失测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810788869.5A CN109238440B (zh) 2018-07-18 2018-07-18 基于最优空时聚焦技术的声学覆盖层插入损失测量方法

Publications (2)

Publication Number Publication Date
CN109238440A CN109238440A (zh) 2019-01-18
CN109238440B true CN109238440B (zh) 2020-05-01

Family

ID=65072174

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810788869.5A Active CN109238440B (zh) 2018-07-18 2018-07-18 基于最优空时聚焦技术的声学覆盖层插入损失测量方法

Country Status (1)

Country Link
CN (1) CN109238440B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109991590B (zh) * 2019-02-21 2021-02-02 中国船舶重工集团公司第七一五研究所 一种在有限空间压力罐内测试换能器低频发射特性的系统与方法
CN112834623B (zh) * 2021-01-06 2022-01-14 浙江大学 一种基于压缩感知技术的声学覆盖层回声降低测量方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS63227490A (ja) * 1987-03-16 1988-09-21 Sumitomo Heavy Ind Ltd 海洋調査船等における水中音響インテンシテイ計測治具
WO2004088355A1 (fr) * 2003-04-01 2004-10-14 Thales Procede de traitement de signaux, et sonar actif le mettant en oeuvre
WO2011134890A1 (de) * 2010-04-26 2011-11-03 Atlas Elektronik Gmbh Trägermodul für eine unterwasserantenne, sowie verfahren zum befestigen eines wandlerelements auf derartigem trägermodul
CN102539530A (zh) * 2012-01-12 2012-07-04 浙江大学 基于时间反转聚焦的水声无源材料回声降低/反射系数测量方法
CN104792877A (zh) * 2015-04-03 2015-07-22 浙江大学 水下去耦抑振材料去耦性能测量方法
CN105181800A (zh) * 2015-08-20 2015-12-23 浙江大学 基于多通道空时逆滤波技术的声学覆盖层回声降低测量方法
CN105223271A (zh) * 2015-09-01 2016-01-06 浙江大学 一种有限空间条件下水声无源材料隔声量时反聚焦测量方法
CN107340051A (zh) * 2016-05-03 2017-11-10 中国科学院声学研究所 一种用于水声低频声源测试的单端面互易装置

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102004029714B3 (de) * 2004-06-21 2006-02-02 Carcoustics Tech Center Gmbh Prüfstand und Verfahren zur Messung der Schalldämmung oder Einfügungsdämmung an einem Prüfobjekt
US7408842B2 (en) * 2006-03-30 2008-08-05 Materials Sciences Corporation Sonar dome

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS63227490A (ja) * 1987-03-16 1988-09-21 Sumitomo Heavy Ind Ltd 海洋調査船等における水中音響インテンシテイ計測治具
WO2004088355A1 (fr) * 2003-04-01 2004-10-14 Thales Procede de traitement de signaux, et sonar actif le mettant en oeuvre
WO2011134890A1 (de) * 2010-04-26 2011-11-03 Atlas Elektronik Gmbh Trägermodul für eine unterwasserantenne, sowie verfahren zum befestigen eines wandlerelements auf derartigem trägermodul
CN102539530A (zh) * 2012-01-12 2012-07-04 浙江大学 基于时间反转聚焦的水声无源材料回声降低/反射系数测量方法
CN104792877A (zh) * 2015-04-03 2015-07-22 浙江大学 水下去耦抑振材料去耦性能测量方法
CN105181800A (zh) * 2015-08-20 2015-12-23 浙江大学 基于多通道空时逆滤波技术的声学覆盖层回声降低测量方法
CN105223271A (zh) * 2015-09-01 2016-01-06 浙江大学 一种有限空间条件下水声无源材料隔声量时反聚焦测量方法
CN107340051A (zh) * 2016-05-03 2017-11-10 中国科学院声学研究所 一种用于水声低频声源测试的单端面互易装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《Measurement of the Echo Reduction for Passive Materials Based on Inverse Filtering Technique》;Suxuan Li;《IEEE》;20160609;全文 *
《水声无源材料测试衍射效应影响抑制技术研究》;刘彦森;《电声技术》;20100417;第34卷(第4期);全文 *

Also Published As

Publication number Publication date
CN109238440A (zh) 2019-01-18

Similar Documents

Publication Publication Date Title
Engholm et al. Direction of arrival estimation of Lamb waves using circular arrays
Prada et al. Separation of interfering acoustic scattered signals using the invariants of the time-reversal operator. Application to Lamb waves characterization
CN102590349B (zh) 基于时间反转聚焦的水声无源材料插入损失/透射系数测量方法
CN110186546B (zh) 基于粉红噪声的水听器灵敏度自由场宽带校准方法
CN105301114B (zh) 一种基于多通道空时逆滤波技术的声学覆盖层插入损失测量方法
Engholm et al. Adaptive beamforming for array imaging of plate structures using lamb waves
CN109991590B (zh) 一种在有限空间压力罐内测试换能器低频发射特性的系统与方法
CN102183710A (zh) 变压器局放源个数辨识方法与系统
CN109238440B (zh) 基于最优空时聚焦技术的声学覆盖层插入损失测量方法
CN109238441B (zh) 一种基于最优空时聚焦技术的声学覆盖层回声降低测量方法
CN110412133A (zh) 一种基于合成孔径聚焦成像的超声阵列混凝土检测系统
Xu et al. Lamb wave based damage imaging under nonlinear chirp excitation
CN112834623B (zh) 一种基于压缩感知技术的声学覆盖层回声降低测量方法
CN105181805A (zh) 一种基于时反music的多滤波超声成像方法
Li et al. Research on the imaging of concrete defect based on the pulse compression technique
Hua et al. Simultaneous excitation system for efficient guided wave structural health monitoring
Li et al. Measurement of acoustic properties for passive-material samples using multichannel inverse filter
Sun et al. Secondary channel estimation in spatial active noise control systems using a single moving higher order microphone
Liu et al. A robust Lamb wave imaging approach to plate-like structural health monitoring of materials with transducer array position errors
Ma et al. Measurement of echo reduction for passive-material samples using sparse Bayesian learning and least squares estimation
Anderson High amplitude time reversal focusing of sound and vibration
Engholm et al. Using 2-D arrays for sensing multimodal lamb waves
Wang Far-field performances prediction of high frequency projectors using secondary source array method
Tsysar et al. Experimental verification of phased receiving waveguide array for ultrasonic imaging in aggressive liquids
O'Donoughue et al. Time reversal beamforming of guided waves in pipes with a single defect

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