一种结合室内空间剖分的声线跟踪法

声线跟踪法是应用于室内声场预测的几何声学模拟经典方法之一,但该方法的计算效率随着室内空间布局复杂度的增加而降低。为了高效地应用于各种室内空间布局,提出了一种结合室内空间剖分的声线跟踪法,该方法将整个室内空间剖分成有限个凸多面体子空间,以减少在声线跟踪过程中与声线求交点的墙面数量,并省略了碰撞点有效性判断,显著提高了计算效率,尤其适用于大型复杂的室内空间布局。     

相干声线跟踪理论中的周期界面迭代散射 仿真方法*

界面声反射模拟是室内复杂声学现象仿真的关键。针对传统声学仿真方法对于周期散射结构存在条件下声场仿真精度较低的问题,本文发展了一种基于迭代散射模型的室内相干声线跟踪法。此方法以经典的相干声线跟踪法为基础,将室内中常见的周期散射结构进行几何形状上的简化处理,然后依据周期散射定理给出声波在界面上的散射方向及能量,并将原始声线迭代分裂为相应的散射子声线,继续对其跟踪处理,此迭代散射模型对周期散射结构上的界面散射现象进行了准确的模拟。数值验证结果表明,本文方法可以有效地在低频段提高室内声场仿真精度,可为具有复杂散射现象的室内仿真提供新思路。     

多层奥氏体钢中超声相控阵辐射声场分析

针对一维线性超声相控阵透过楔块,在多层不同晶轴取向的奥氏体钢中的声场辐射问题,结合高斯声束等效点源模型以及基于时间最小原理的射线追踪法,给出了各层介质中透射声场的计算方法。将奥氏体钢近似为横向各向同性介质,计算了相控阵透过楔块,在晶轴取向为0°的奥氏体钢中的无延时纵波声场,计算结果与COMSOL仿真结果相吻合,从而验证了所用方法的正确性。通过加入不同的延迟法则,仿真计算了三层含有不同晶轴取向(0°,30°,100°)的奥氏体钢中的纵波声场,实现了相控阵声场的偏转与聚焦,分析了偏转声场与聚焦声场的传播特性。仿真结果表明,不同的晶轴取向将导致不同的声束偏转以及聚焦效果。通过延迟激励各阵元,虽然可以控制声束偏转或聚焦到预定位置,但是晶轴取向仍会对声束宽度以及幅值产生一定的影响。      

Variations of mode amplitudes in a range-dependent waveguide

The ray method of calculating the mode amplitudes is used to analyze the sound fields in deep-water acoustic waveguides with two types of inhomogeneities of the refractive index: (i) weak inhomogeneities that cause small perturbations of ray trajectories and (ii) strong inhomogeneities with large spatial scales. Simple analytical relations are derived for describing the variations of the mode structure of the sound field in the presence of the aforementioned inhomogeneities. A new criterion defining the validity of the adiabatic approximation is formulated. To illustrate and test the results obtained, a numerical simulation of the sound fields is performed on the basis of the parabolic equation method.

     

多层介质中超声相控阵声场仿真*

由于良好的声束偏转与聚焦特性,超声相控阵已经广泛应用于多层固体介质的缺陷检测。当超声束经过多层介质时,由于反射、透射以及模式转换的存在,多种声束存在于这种结构中,使得声场分析变得复杂。为了提高多层介质检测的准确性,有必要对超声声场的分布规律进行深入地了解。该文结合高斯声束等效点源模型以及射线追踪法,给出了相控阵声源在多层固体介质中激发声场的仿真方法,并且模拟计算了一维线型相控阵在楔块-铝-黄铜-钢四层固体介质中的辐射声场。通过对不同延时法则的计算,实现了声波在这种复杂介质中的偏转与聚焦,进而研究了不同焦点处聚焦声场的分布。结果表明:相控阵方法能使聚焦点处的声场幅值增大,能量集中,提高了检测分辨率;不同聚焦点处声场聚焦效果不同,实际检测时应根据检测区域结构及位置特点,合理放置相控阵换能器。与瑞利积分法的比较表明,该文的仿真方法适用于多层介质相控阵声场的计算。     

Tomography in observation of inhomogenenities in shallow sea during sounding by a focused high-frequency acoustic field: II. Study of the observation efficiency

We study the efficiency of high-frequency acoustic sounding of small-size spatially localized inhomogeneities in shallow sea by a tomographic method. The method is based on that the focusing of the received and radiated fields to the region of the assumed location of the inhomogeneity is matched to the waveguide. During the observation, the statistical check is performed for hypotheses on the object location that are formed using a priori information on the waveguide model and data on the observation-system structure. To decrease the influence of the interference structure of the acoustic fields in the waveguide, we use accumulation of partial images corresponding to separate ray tomographic projections. Using a simulation model of the considered method, we analyze its accuracy and sensitivity under different observation conditions in shallow sea. Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Radiofizika, Vol. 52, No. 3, pp. 192–207, March 2009.     

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

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

Spatial correlation of the high intensity zone in deep-water acoustic field

The spatial correlations of acoustic field have important implications for underwater target detection and other applications in deep water. In this paper, the spatial correlations of the high intensity zone in the deep-water acoustic field are investigated by using the experimental data obtained in the South China Sea. The experimental results show that the structures of the spatial correlation coefficient at different ranges and depths are similar to the transmission loss structure in deep water. The main reason for this phenomenon is analyzed by combining the normal mode theory with the ray theory. It is shown that the received signals in the high intensity zone mainly include one or two main pulses which are contributed by the interference of a group of waterborne modes with similar phases. The horizontal-longitudinal correlations at the same receiver depth but in different high intensity zones are analyzed. At some positions, more pulses are received in the arrival structure of the signal due to bottom reflection and the horizontal-longitudinal correlation coefficient decreases accordingly.The multi-path arrival structure of receiving signal becomes more complex with increasing receiver depth.     

凹面相控阵聚焦声场在液固界面上的反射和折射

开展了凹面线性相控阵列辐射声场在液固界面上的反射和折射特性研究,利用射线近似方法,得到了凹面相控阵聚焦声场在液固界面上反射和折射后的声场渐近解析表达式,对声波在液固曲面上的反射和折射声场进行了分析和讨论. 利用这个解析表达式,对凹面线性聚焦声场在液固平界面和液固圆柱界面情况下固体中折射纵波和折射横波的相控阵声场进行了分析和讨论,发现凹面线阵声场在液固圆柱界面下比液固平界面具有更好的聚焦效果. 关键词： 凹面阵列 超声相控阵 反射与折射     

浅海空气声入水传播特性分析

结合射线和波数谱积分方法,对空气声入水传播途径进行了分析,利用海上试验数据进行了比较检验。结果表明,在浅海环境中,对水下声场有主要贡献的空气声入水传播途径,主要是透射穿过海面边界的折射直达声以及后续的海底反射声途径,其中折射直达声途径的贡献主要集中在声源正下方附近区域,当距离较远时,由于声线扩展损失效应以及直达声影区两方面的限制,折射直达声传播损失显著增加,对接收声场起主要贡献的是可以到达更远水平距离上的海底反射声,包括海底海面多次反射声。      

Acoustic field structure in the first oceanic convergence zone for different frequencies in the audio range

Experimental data on the angular and energy structures of the acoustic field generated in the first convergence zone by a source of continuous pseudo-noise radiation in the frequency range 0.5–4.0 kHz are presented. The experiments are carried out in a tropical region of the Central Atlantic. The acoustic field characteristics are studied with the omnidirectional and highly directional (～2°) reception in the vertical plane. Particular attention is given to studying the field structure at the entrance to the convergence zone at frequencies of 1.25 and 3.15 kHz. The experimental data are compared with the ray calculations. In the experiments, the origin of the convergence zone, which is determined by a sharp increase in the received signal intensity, is found to be about 1.2 km nearer to the source than in the calculations. At a frequency of 1.25 kHz, the convergence zone begins 300 m nearer than at the higher frequency (3.15 kHz). At distances of several kilometers to the zone origin, weak signals that arrive at the same angles as the signals forming the origin of the convergence zone are detected.     

Restless rays, steady wave fronts

Observations of underwater acoustic fields with vertical line arrays and numerical simulations of long-range sound propagation in an ocean perturbed by internal gravity waves indicate that acoustic wave fronts are much more stable than the rays comprising these wave fronts. This paper provides a theoretical explanation of the phenomenon of wave front stability in a medium with weak sound-speed perturbations. It is shown analytically that at propagation ranges that are large compared to the correlation length of the sound-speed perturbations but smaller than ranges at which ray chaos develops, end points of rays launched from a point source and having a given travel time are scattered primarily along the wave front corresponding to the same travel time in the unperturbed environment. The ratio of root mean square displacements of the ray end points along and across the unperturbed wave front increases with range as the ratio of ray length to correlation length of environmental perturbations. An intuitive physical explanation of the theoretical results is proposed. The relative stability of wave fronts compared to rays is shown to follow from Fermat's principle and dimensional considerations.     

Acoustic band pinning in the phononic crystal plates of anti-symmetric structure

Acoustic bands are studied numerically for a Lamb wave propagating in an anti-symmetric structure of a one-dimensional periodic plate by using the method of supercell plane-wave expansion. The results show that all the bands are pinned in pairs at the Brillouin zone boundary as long as the anti-symmetry remains and acoustic band gaps (ABGs) only appear between certain bands. In order to reveal the relationship between the band pinning and the anti-symmetry, the method of eigenmode analysis is introduced to calculate the displacement fields of different plate structures. Further, the method of harmony response analysis is employed to calculate the reference spectra to verify the accuracy of numerical calculations of acoustic band map, and both the locations and widths of ABGs in the acoustic band map are in good agreement with those of the reference spectra. The investigations show that the pinning effect is very sensitive to the anti-symmetry of periodic plates, and by introducing different types of breakages, more ABGs or narrow pass bands will appear, which is meaningful in band gap engineering.     

三维浅海下弹性结构声辐射预报的有限元-抛物方程法

利用多物理场耦合有限元法对结构和流体适应性强、抛物方程声场计算高效准确的特点,提出了三维浅海波导下弹性结构声振特性研究的有限元-抛物方程法.该方法采用多物理场耦合有限元理论建立浅海下结构近场声辐射模型,计算局域波导下结构声振信息,并提取深度方向上复声压值作为抛物方程初始值;然后采用隐式差分法求解抛物方程以步进计算结构辐射声场.重点介绍了该方法对浅海下结构声辐射计算的准确性、高效性以及快速收敛性后,对Pekeris波导中有限长弹性圆柱壳的声振特性进行了分析.研究得出,当圆柱壳靠近海面(海底)时,其耦合频率比自由场下的要高(低),当潜深达到一定范围时,与自由场耦合频率基本趋于一致;在低频远场,结构辐射场与同强度点源声场具有一定的等效性,且等效距离随着频率增加而增加;由于辐射声场受结构振动模态、几何尺寸和简正波模式影响,结构辐射场传播的衰减规律按近场声影响区、球面波衰减区、介于球面波和柱面波衰减区、柱面波衰减区四个扩展区依次进行.     

Space anisotropy of bottom reverberation in the coastal zone

A method for calculating the bottom reverberation in a coastal zone with allowance for the three-dimensional bottom topography, which affects the spatial anisotropy of the incoming reverberation signal, is proposed. The method is based on the ray approach in combination with the use of the ray invariant. Examples of calculating the monostatic bottom reverberation level for one of the coastal regions with a complex bottom topography including the coastal line and an underwater elevation are presented.     

Dynamical energy analysis for built-up acoustic systems at high frequencies

Standard methods for describing the intensity distribution of mechanical and acoustic wave fields in the high frequency asymptotic limit are often based on flow transport equations. Common techniques are statistical energy analysis, employed mostly in the context of vibro-acoustics, and ray tracing, a popular tool in architectural acoustics. Dynamical energy analysis makes it possible to interpolate between standard statistical energy analysis and full ray tracing, containing both of these methods as limiting cases. In this work a version of dynamical energy analysis based on a Chebyshev basis expansion of the Perron-Frobenius operator governing the ray dynamics is introduced. It is shown that the technique can efficiently deal with multi-component systems overcoming typical geometrical limitations present in statistical energy analysis. Results are compared with state-of-the-art hp-adaptive discontinuous Galerkin finite element simulations.     

一种考虑声线弯曲的温度场重建算法

声波在非均匀温度场中传播时会由于声波的折射而产生声线弯曲现象。为提高非均匀温度场声学CT (Computer Tomography)重建精度,提出一种考虑声线弯曲的重建算法。首先用最小二乘法获得一个不考虑声线弯曲的重建温度场,然后用打靶-插值法确定本征声线出射角,用三角形前向展开法追踪声线,获得此温度场中声波发射/接收器间的本征声线,建立本征声线上声波传播时间与温度分布间的关系,进而实现考虑声线弯曲的最小二乘法温度场重建。单峰、双峰温度场仿真重建结果表明:所采用的本征声线追踪法计算简单运行速度快;考虑声线弯曲后,可明显提高温度场重建精度。因此所提重建算法能更好地适应实际温度场重建需求。      

Regions where transient signals are influenced between a source and receiver

At infinite frequency, the only place that a medium can influence the waveform of a received signal that is emitted from a source is along one or more infinitesimally thin ray paths. For any transient signal at finite frequencies, an exact method is developed to compute the regions in a medium that significantly influence the received signal for any specified window of signal travel time. This window is sometimes chosen to surround a peak. Results at finite frequencies differ from those at infinite frequency because of diffraction. The method has its foundation in the integral theorem of Helmholtz and Kirchhoff. Part of the method involves a filter that yields an imperfect but apparently useful picture of influential regions in the presence of interfering waves. The method is useful for quantifying differences between the region of influence and a ray, and for identifying regions in which medium fluctuations significantly influence signal aberrations at a receiver. Four principal results are found at low frequencies. They are: 1) For propagation in homogeneous media, a significant portion of the received signal is influenced by waves that traverse paths that are approximately integer and half-integer numbers of cycles greater than the straight path between source and receiver. Such paths are called 'constructive and destructive paths of influence', respectively. They correspond to edge-diffracted rays for the geometrical theory of diffraction. 2) For reflection from a flat interface in an otherwise homogeneous medium, the received signal is significantly influenced by constructive and destructive paths of influence whose angles of incidence and reflection differ (non-specular reflection). 3) For acoustic propagation centered at 100 Hz in an oceanic acoustic waveguide, the region of influence markedly departs from a ray path, particularly near the reflective ocean surface. The influential region is flat for O(10) km instead of O(1) km for a ray. 4) The first Fresnel zone is an inappropriate scale to characterize the region of influence for transient signals near a steep ray in inhomogeneous media, as assumed by at least one scattering theory. Modification of that theory may yield a better fit with data.     

Near fields scattered by an underwater finite cylindrical baffle

In this work,acoustic vector characteristics of near fields scattered by an underwater finite cylindrical baffle are investigated theoretically and experimentally.The analytic expressions for the scattered pressure and particle velocity are derived using the elastic thin shell theory.Calculations are presented for the scattered near fields of the pressure,the particle velocity and the intensity.It is found that the pressure and the particle velocity fields near the surface of the cylindrical baffle are characterized by complex interference structure,particle velocity directions and the source bearings are not consistent.The phase difference between the pressure and the particle velocity is not zero and the intensity vector does not reflect the sound bearings.It can be noted that the distortions of the fields will make the original vector signal processing method based on the free space assumption be no longer applicable in the presence of the cylindrical baffle.These results can serve as a basis of the application for the acoustic vector sensor on board.