联合波叠加法的全息理论与实验研究

当空间声场中同时存在多个相干声源时，运用常规近场声全息方法无法重建每个相干声源表面的声学信息，当然也无法预测每个声源单独产生的空间声场，相干声场的全息重建与预测已成为全息技术推广应用过程中亟待解决的问题.在提出联合波叠加法并将其应用于空间声场变换的基础上，对其进行了实验研究.通过对实际相干声场的全息重建与预测，验证了常规波叠加法在相干声场重建中的局限性、联合波叠加法在相干声场全息重建与预测过程的可行性和准确性，还研究了Tikhonov正则化方法在抑制声学逆问题的非适定性中的有效性和滤波系数的选择原则的可行性，以提高全息重建与预测的精度. 关键词： 近场声全息 联合波叠加 相干声场 Tikhonov正则化     

单层传声器阵列重建相干声场的近场声全息方法研究

针对传声器阵列两侧存在相干声源的非自由声场重建问题,提出基于球面谐波函数扩展近场声全息理论的相干声场重建方法。该方法在已知测量面两侧声源几何位置时,使用单层传声器阵列获取测量面处的声压分布,通过最小二乘法获得与目标声源和干扰噪声源响应对应的最优球波函数扩展项数和最优系数向量,结合测点位置的空间坐标进行声波分解,并分别重建出各声源在测量面上的声压分布。为了验证方法的有效性,分别给出了相干噪声源为球形声源和非球形声源的仿真验证,并在全消声室内对双扬声器产生的相干声场的重建进行了实验验证。结果表明:该方法对球形声源和非球形声源干扰下的声场重建都具有较好的效果,球形声源干扰下的重建精度更高。      

基于波叠加法的近场声全息空间分辨率增强方法

提出一种基于波叠加法的近场声全息空间分辨率增强方法. 该方法在波叠加法的基础上,利用全息面声压信号求得布置在全息面附近的虚源面上的简单源源强,再根据求得的简单源源强实现对全息面声压的插值,进而利用插值后的全息面声压数据进行重建. 该方法可以提高近场声全息重建图像的空间分辨率,减少测量工作量,简化测量过程. 通过仿真对影响插值结果的参数进行了分析,给出了合理的选取范围;通过仿真和实验研究验证了该方法的有效性. 关键词： 波叠加法 近场声全息 空间分辨率     

Computation of acoustic radiation from vibrating structures in motion

The noises from a vibrating structure in motion are often encountered in engineering practice, for example, tire noise and pass-by noise of moving vehicles. Consequently, research on the radiation characteristics of moving acoustic source is of significance. In this paper, a new computational method based on the wave superposition approach is developed for the acoustic field from a vibrating structure in motion. It inherits the advantages of the wave superposition approach in the acoustic computation, and in which a method of moving simple sources is used to eliminate the influence of the Doppler effect. By the proposed method, the acoustic radiation from the moving vibrating structure can be calculated easily with the same implementation process as the conventional wave superposition approach performed in the stationary acoustic field. Finally, the validity and feasibility of the proposed method are verified by the numerical results.     

空间声场全息重建的波叠加方法研究

提出了基于波叠加法的近场声场全息技术，并将其用于任意形状物体的声辐射分析.在声辐射计算问题中，边界元法是通过离散边界面上的声学和位置变量来实现，而波叠加方法则通过叠加辐射体内部若干个简单源产生的声场来完成.因而，基于波叠加法的声全息就不存在边界面上的参数插值和奇异积分等问题，而这些问题是基于边界元法的声全息所固有的.与基于边界元法的声全息相比较，基于波叠加法的声全息在原理上更易于理解，在计算机上更容易实现.实验结果表明：该种全息技术在重建声场时，具有令人满意的重建精度. 关键词： 声全息 逆问题 波叠加方法 正则化方法     

波叠加声场重构精度的影响因素分析

波叠加法是利用等效源的思想来求解声辐射问题的一种数值方法.其原理并没有确切地规定如何配置等效源.正由于这些不确切性的存在,当将它应用到声场重构时,有必要检验重构的准确性,分析影响重构精度的各种因素,目的是通过分析而获得一些提高重构精度的规律.通过理论分析和数值仿真,分析了各因素对重建准确性的影响.证明等效源必须配置在辐射体内部,而且最好能与辐射体表面共形;但测量面并不需要与辐射体表面共形,而且也不需要在近场测量;还发现,Tikhonov和TSVD正则化方法很适合与波叠加法结合用于声场重构.依据这些规律能更好地应用波叠加法进行声场重构,也将拓宽其应用范围. 关键词： 声场重构 波叠加法 误差影响因素 正则化方法     

基于波叠加法的patch近场声全息及其实验研究

提出了基于波叠加法的patch近场声全息技术(PNAH).该方法首先通过线性叠加在振动体内部“虚置”的若干不同强度的简单源产生的波场来拟合实际声场，实现近场外推，然后采用计算简便的FFT方法进行全息重建.该方法解决了现有PNAH技术近场外推过程中存在的外推精度和计算效率问题，具有外推精度高、计算速度快、实现方法简单等特点，更适用于工程计算.为了验证本方法的正确性与有效性，将其用于单点激励固支板的全息实验中，取得了令人满意的结果.     

采用球面波叠加法还原自由声场的方法

为解决非自由声场中近场声全息重建时,干扰声在目标声源表面产生的散射影响,提出一种基于球面波叠加法的自由场还原技术。该技术首先采用基于球面波叠加法的声场分离技术获得向内和向外传播的声场,然后以目标声源的表面导纳作为边界条件,实现目标声源辐射声和散射声的分离,从而获得等效于自由声场的测量条件。该技术为准确实现非自由声场中的噪声源识别创造了条件。文中首先详细描述了该技术的基本原理,并提出一种最优球面波展开项数选取方法,最后通过数值仿真说明了该技术的有效性。结果表明:在频率较低时,散射声影响较小,采用声场分离技术和自由场还原技术效果相当;但随着频率升高,散射声影响逐步增强,必须采用自由场还原技术才能准确获得目标声源辐射声。      

Recovery of the free field in noisy environment by using the spherical wave superposition method

To remove the scattering effect of the disturbing sound on the target source when implementing nearfield acoustic holography in a non-free field, a free field recovery technique based on the spherical wave superposition method is proposed. In the method, the sound field separation technique based on the spherical wave superposition method is first used to separate the incoming and outgoing fields, and a further step for separating the radiated and scattered fields is performed by utilizing the surface admittance of the target source as the boundary condition. The technique makes it possible to correctly identify noise sources in a non-free sound field. The basic principle of the technique is described firstly, a method for choosing the optimal number of spherical wave expansion terms is given, and two numerical simulations are used to demonstrate the validity of this technique. It is shown that, for the lower frequency, the scattering effect can be neglected, and the radiated field of the target source can be obtained by the sound field separation technique, however, as the increasing of the frequency, the scattering effect cannot be neglected, and the free field recovery technique has to be used to obtain the radiated field of the target source.     

基于联合波叠加法的浅海信道下圆柱壳声辐射研究

针对浅海信道下弹性结构声辐射预报尚无高效可靠的研究方法,提出了一种浅海信道下弹性结构声辐射快速预报的联合波叠加法.该方法结合了浅海信道传输函数、多物理场耦合数值计算法和波叠加法理论,运用该方法可对浅海信道下弹性结构辐射声场进行快速预报.经数值法和解析解法验证后,从信道下辐射源、环境影响和辐射声场测量的角度研究分析了浅海信道下弹性圆柱壳的声辐射特性,阐释了进行浅海信道下结构声辐射研究的必要性.研究结果表明,仅在低频浅海信道下弹性结构可近似等效为点源,信道上下边界对声场产生显著的耦合影响,高频段的空间声场指向性分布尤为明显,垂直线列阵进行信道下结构辐射声功率测量时,测量结果受到信道环境边界和潜深的影响较大.     

采用振速测量的改进统计最优近场声全息方法

针对常规统计最优近场声全息在空间多源声场重建过程中所需波函数项数多、重建精度不理想的问题,本文提出了一种基于振速测量的改进统计最优近场声全息算法。与常规算法不同,改进算法主要根据声源特点选取合适的波函数组合来计算声场传递矩阵。通过数值仿真初步验证了该方法的准确性和有效性,并与常规算法进行了详细的对比分析。仿真结果表明,改进算法重建精度高,随频率的变化相对误差波动较小,且随着频率的增大相对误差有逐渐减小的趋势;此外,不同的波函数组合,重建效果差异很大,当选取的波函数与声源共形且数量一致时重建效果最好。     

三维复杂温度场高精度声学测量方法

声学温度场检测技术通过多路径声波传播时间数据,反演被测区域的温度分布。提供了一种高精度的三维复杂温度场的声学测量方法。首先从射线声学角度给出了三维非均匀温度场中声波传播路径的数学模型。在此基础上,将三维温度场的重建问题转化为声波传播路径的求解和温度场的反演问题,建立了基于多项式修正径向基函数(RBF-PR)和改进的Tikhonov正则化三维温度场重建算法。采用两种典型的炉膛三维温度场模型,在信噪比SNR=35 dB下进行了数值模拟,分析了声波传播路径在非均匀温度场中的弯曲特性、算法的重建质量和抗噪性,同时进行了实际炉膛内二维温度场的重建。结果表明了提出的考虑声线弯曲的温度场重建算法具有精度高,抗噪性强、适用性好的特点,为实现高精度的复杂温度场的声学测量提供了有效方法。      

利用宽带引导声源重构拷贝声场的目标声源定位

为了避开传统的匹配场目标定位技术对环境先验知识的依赖性,提出了在均匀浅海环境中只知道少量环境参数的情况下,利用垂直接收阵和不同距离上的两枚宽带引导声源重构声场对目标声源进行定位的一种方法。这种方法主要基于简正波估计和声场重构两种关键技术,同时省去了匹配场定位技术中大量的拷贝声场计算。数值仿真主要采用线性Bartlett匹配处理器分析了目标定位效果,在信噪比高于10 d B的情况下,定位效果良好。     

基于统计最优和波叠加的联合局部近场声全息

提出了一种基于统计最优近场声全息和波叠加法的联合局部近场声全息技术。首先利用两次统计最优近场声全息的声源定位结果来指导配置等效源,其后利用波叠加法进行局部声场重构。该技术适合于中低频声场的局部重建,计算快速,重建精度高;可以在测量数据有缺失的情况下重建声场。进行了脉动球声源模型的数值仿真,并在半消声室内对电机噪声源进行了实验,仿真实验都准确地重构了声源所辐射的外部声场。该技术可以重建任意类球形声源辐射的声场。      

A wave field synthesis approach to reproduction of spatially correlated sound fields

This article discusses an open-loop wave field synthesis (WFS) approach for the reproduction of spatially correlated sound fields. The main application concerns laboratory reproduction of turbulent boundary layer wall pressure on aircraft fuselages and measurement of their sound transmission loss. The problem configuration involves reconstruction of random sound pressure distributions on a planar reproduction surface using a planar array of reproduction monopoles parallel to the reproduction plane. In this paper, the WFS formulation is extended to sound fields with imposed time and spatial correlation properties (or equivalently imposed cross-spectral density in the frequency and wave number domains). Numerical examples are presented for the reproduction of a propagating plane wave, diffuse acoustic field and wall pressure in subsonic or supersonic turbulent boundary layers. The reproduction accuracy is examined in terms of the size of the source plane and reproduction plane, their separation, and the number of reproduction sources required per acoustic wavelength. While the reproduction approach cannot reconstruct sub-wavelength correlation scales of subsonic turbulent boundary layers, it effectively reconstructs correlation scales larger than the acoustic wavelength, making it appropriate for diffuse acoustic field and supersonic turbulent layers.     

适用于循环平稳声场的基于波叠加法的近场声全息技术

在原有的平面循环平稳近场声全息基础上，提出一种基于波叠加法的循环平稳近场声全息技术，可以对具有复杂表面的声源进行全息重建，重建的声源表面声压谱相关密度函数能反映出调制信号的信息.声源表面声压谱相关密度函数全息图形象地反映了调制信号在表面的强弱分布情况，可由此确定调制信号源的产生位置.仿真分析和实验验证表明，基于波叠加法的循环平稳近场声全息技术可以更准确地反映循环平稳声场的调制特性.该方法继承了波叠加法的优点，无需计算边界奇异积分，计算效率高、精度好.     

适用于循环平稳声场的基于波叠加法的近场声全息技术

在原有的平面循环平稳近场声全息基础上,提出一种基于波叠加法的循环平稳近场声全息技术,可以对具有复杂表面的声源进行全息重建,重建的声源表面声压谱相关密度函数能反映出调制信号的信息.声源表面声压谱相关密度函数全息图形象地反映了调制信号在表面的强弱分布情况,可由此确定调制信号源的产生位置.仿真分析和实验验证表明,基于波叠加法的循环平稳近场声全息技术可以更准确地反映循环平稳声场的调制特性.该方法继承了波叠加法的优点,无需计算边界奇异积分,计算效率高、精度好. 关键词： 近场声全息 循环平稳信号 波叠加     

Direct analysis of dispersive wave fields from near-field pressure measurements

Flexural waves play a significant role for the radiation of sound from plates. The analysis of flexural wave fields enables the detection of sources and transmission paths in plate-like structures. The measurement of these wave fields can be carried out indirectly by means of near-field acoustic holography, which determines the vibrational wave field from pressure information measured in a plane close to the plate under investigation. The reconstruction of the plate vibration is usually obtained by inverting the forward radiation problem, i.e., by inversion of an integral operator. In this article, it is shown that a pressure measurement taken in the extreme near-field of a vibrating plate can directly be used for the approximate analysis of the dispersive flexural wave field. The inversion step of near-field acoustic holography is not necessarily required if such an approximate solution is sufficient. The proposed method enables fast and simple analysis of dispersion characteristics. Application of dispersion compensation to the measured field allows for visualizations of propagating wavefronts, such that sources and scatterers in the plate can be detected. The capabilities of the described approach are demonstrated on several measurements.     

Wave envelopes method for description of nonlinear acoustic wave propagation

A novel, free from paraxial approximation and computationally efficient numerical algorithm capable of predicting 4D acoustic fields in lossy and nonlinear media from arbitrary shaped sources (relevant to probes used in medical ultrasonic imaging and therapeutic systems) is described. The new WE (wave envelopes) approach to nonlinear propagation modeling is based on the solution of the second order nonlinear differential wave equation reported in [J. Wójcik, J. Acoust. Soc. Am. 104 (1998) 2654-2663; V.P. Kuznetsov, Akust. Zh. 16 (1970) 548-553]. An incremental stepping scheme allows for forward wave propagation. The operator-splitting method accounts independently for the effects of full diffraction, absorption and nonlinear interactions of harmonics. The WE method represents the propagating pulsed acoustic wave as a superposition of wavelet-like sinusoidal pulses with carrier frequencies being the harmonics of the boundary tone burst disturbance. The model is valid for lossy media, arbitrarily shaped plane and focused sources, accounts for the effects of diffraction and can be applied to continuous as well as to pulsed waves. Depending on the source geometry, level of nonlinearity and frequency bandwidth, in comparison with the conventional approach the Time-Averaged Wave Envelopes (TAWE) method shortens computational time of the full 4D nonlinear field calculation by at least an order of magnitude; thus, predictions of nonlinear beam propagation from complex sources (such as phased arrays) can be available within 30-60 min using only a standard PC. The approximate ratio between the computational time costs obtained by using the TAWE method and the conventional approach in calculations of the nonlinear interactions is proportional to 1/N2, and in memory consumption to 1/N where N is the average bandwidth of the individual wavelets. Numerical computations comparing the spatial field distributions obtained by using both the TAWE method and the conventional approach (based on a Fourier series representation of the propagating wave) are given for circular source geometry, which represents the most challenging case from the computational time point of view. For two cases, short (2 cycle) and long (8 cycle) 2 MHz bursts, the computational times were 10 min and 15 min versus 2 h and 8 h for the TAWE method versus the conventional method, respectively.