水声共形阵辐射指向性计算方法及其实验验证

用边界元方法,给出了水声发射换能器共形阵各阵元的单元指向性以及基阵的自然波束指向性和束控波束指向性。首先建立了换能器共形阵的边界元模型,并根据换能器结构和声传播介质建立了模型的边界条件,计算了其辐射指向性。然后在消声水池中对共形阵各阵元以及基阵的辐射指向性进行了实验测量。将理论计算结果与实验测量结果进行了对比,二者吻合得比较好。研究结果表明:基于边界元方法建立的模型及计算方法,对基阵的辐射指向性进行理论计算和预报是可行的。     

宽带密排平面阵声辐射分析

应用声有限元-边界元法计算了宽带基阵声场分布.建立了平面9元阵有限元-边界元模型,给出了阵元之间互阻抗数值计算方法,分析了阵元间互辐射对声性能影响的变化规律,计算了9元阵3个谐振频率点处指向性和波束宽度.利用3谐振换能器研制出一种宽带密排平面9元阵,并在消声水池中对其输入阻抗和指向性进行了实验测量.研究结果表明,按阵元...     

非近轴近似多高斯声束模型的相控阵换能器声场计算

针对近轴近似的多高斯声束模犁无法准确计算离轴区域声场的问题,我们提出了一个矩形源的非近轴近似多高斯声束模型,用于计算偏转聚焦时线性相控阵换能器辐射的声场.利用该方法计算了由矩形阵元组成的相控线阵的辐射声场,并与瑞利积分、近轴近似方法及夫琅和费近似方法进行了分析和对比,计算结果表明该方法克服了常规近轴近似条件的限制,可以准确、快速地计算较大偏转范围内的换能器辐射声场,在计算相控阵换能器偏转聚焦声场中具有独特的优势.     

cMUT阵元的声学分析与环形扫描成像

电容式微机械超声换能器(cMUT)具有宽带宽、灵敏度高、机械阻抗低和与电子电路集成制造等优点.由于超声换能器声场设计的好坏决定了成像质量的优劣,为了明确阵元参数与辐射声场的关系,该文对不同的阵元参数进行了计算仿真.与压电陶瓷超声换能器的结构不同,cMUT阵元是由多个电容单元(cell)并联构成.因此分析了cell的半径...     

带楔块二维面阵列超声相控阵声场特性分析

开展了带楔块二维面阵脉冲超声相控阵辐射声场特性的研究。将安装在一定角度斜楔块上二维面阵相控阵换能器声场问题简化为液固界面情况进行讨论。以单阵元在液固平界面条件下的辐射声场为基础,推导了聚焦偏转法则,给出了带楔块二维面阵超声相控阵声场计算方法。以检测材料为钢板,安装在倾斜角为36°有机玻璃楔块上的频率为5 MHz、8×8二维面阵相控阵的辐射声场计算为例,分析了在不同偏转角和不同聚焦深度下检测材料中的辐射声场特性。计算结果表明该方法可有效的分析带楔块二维面阵超声相控阵声场特性并用于指导二维超声面阵角束探头的设计。      

二维超声阵列换能器声场的两种时域仿真模型比较

针对二维超声阵列换能器的声场仿真,以xMATRIX超声换能器结构为原型,设计仿真实验对分割阵元空间脉冲响应模型和离散点阵代表模型的计算精度、计算效率、影响参数和适用场合进行了分析比较.结果表明,分割阵元空间脉冲响应模型计算速度快,精度高,对阵元尺寸、阵元间距、聚焦深度及阵元高宽比例敏感,广泛应用于均匀介质中的声场仿真....     

热疗用环形平面换能器声场特性的研究

引入负声源，导出了环形平面换能器声场的简洁表示式，计算了单个圆环、多个同心圆环同时激励及环阵相控聚焦声源产生的声场，分析了环的宽度、环数及聚焦点位置等因素对声场特征的影响。提出根据各阵元在声场中某点产生的声压的相位来确定各阵元电激励信号的延时，此方法较根据各阵元到该点的平均声程确定各阵元电激励信号的延时，能更难确地在该点实现聚焦。     

求解随机振动结构声辐射的统计边界元方法

机器的机械振动在音频范围内大多是具有明显随机性的随机过程，普通边界元方法不能计算这类具有随机边界条件的结构声辐射问题。本文提出将边界元方法和统计方法相结合的统计边界元方法，导出了求解随机振动物体所辐射噪声声场的计算公式。通过求解随机声场算例的计算，表明本文所提出的方法是正确的。     

Pekeris波导中垂直发射阵的阵增益快速估计

为了利用单换能器的声强级估计垂直发射阵的声强级,研究了Pekeris波导中快速估计发射阵增益的方法。首先,基于等效深度简正波理论,计算了Pekeris波导中垂直发射阵的平均声强,进而得到用模态阶数求和表示的发射阵增益的表达式。然后对模态阶数用积分代替求和,并对引入的误差进行修,得到3/2次方衰减场区发射阵增益的快速估计表达式;有效声影区声场由于只剩1阶简正模态,可直接简化得到发射阵增益的快速估计表达式,在海底参数未知时可用海深近似代替等效深度来估计有效声影区的发射阵增益。计算机仿真和海上实验得到的阵增益快速估计误差均小于1 dB,验证了快速估计方法的有效性,这一发射阵增益快速估计方法可以为工程应用提供参考。      

微纳相控线阵超声换能器参数的理论分析*

微纳相控线阵超声换能器参数(阵元数目、阵元宽度及阵元间距)直接影响其横向声场分布,而其横向声场分布是能否实现高成像分辨率、大探测深度的决定性因素,也是制备换能器的主要依据。该文利用数值模拟研究微纳相控线阵超声换能器阵元参数对其横向声场中主瓣强度、-3 dB主瓣宽度、第一级旁瓣及栅瓣的影响。结果表明,主瓣强度随着阵元数目增加而增大,随阵元间距减小而增大,随着阵元宽度的增大呈现先增大后减小的趋势;-3 dB主瓣宽度随着阵元数目和阵元间距的增大而减小,随着阵元宽度的减小而减小;此外,减小阵元数目、减小阵元间距或增大阵元宽度均可以抑制旁瓣;栅瓣在阵元间距满足一定条件时可以完全消除。通过这些研究为微纳相控线阵超声换能器的优化设计与制备提供理论参考。     

Directivity calculation with experimental verification for a conformal array of underwater acoustic transducers

The radiation directivity of a complicated conformal array of underwater acoustic transducers is presented based on the boundary ＂element method. It includes the element directivity of each transducer, the natural beam pattern and the controlled beam pattern of the transducer array. At first, the boundary element model of the conformal array is built up, and then the boundary condition is exerted on the model according to the design and environment in which the transducer array is used, and the radiation directivity of the conformal array is calculated. An experiment has been done to measure the directivity in an anechoic water tank. The calculated and the experimental results are compared and analyzed. They are consistent to each other. It is shown that the boundary element method together with the detailed calculations is successful to simulate and predict the radiation directivity of an underwater acoustic transducer array.     

Optimization of transmitting beam patterns of a conformal transducer array

A method is presented to calculate the driving-voltage weighting vector of a conformal array of underwater acoustic transmitting transducers to obtain a low-sidelobe beam pattern based on the measured receiving array manifold. The relationship among three quantities is given, which are, respectively, the radiated acoustic field, the measured receiving array manifold matrix and the driving-voltage weighting vector of the transducer array. Then, the driving-voltage weighting vector of the array is calculated using the optimization method to obtain a low-sidelobe transmitting beam pattern. At the frequency of 12.5 kHz, the receiving array manifold matrix of a 27-element conformal array is measured in an anechoic water tank. The driving-voltage weighting vector of the array is calculated using the proposed method. In addition, the computer simulation and experiments are carried out. The results agree well and show that the proposed method can obtain a low-sidelobe transmitting beam pattern and at the same time provide the largest amplitude of pressure in the axial direction when the maximum amplitude of the driving voltages of the array elements keeps unchanged.     

A numerical study on the flow and sound fields of centrifugal impeller located near a wedge

Centrifugal fans are widely used and the noise generated by these machines causes one of the serious problems. In general, the centrifugal fan noise is often dominated by tones at blade passage frequency and its higher harmonics. This is a consequence of the strong interaction between the flow discharged from the impeller and the cut-off in the casing. However, only a few researches have been carried out on predicting the noise because of the difficulty in obtaining detailed information about the flow field and considering the scattering effect of the casing. The objective of this study is to understand the generation mechanism of sound and to develop a prediction method for the unsteady flow field and the acoustic pressure field of the centrifugal impeller. A discrete vortex method is used to model the centrifugal impeller and a wedge and to calculate the flow field. The force of each element on the blade is calculated by the unsteady Bernoulli equation. Lowson's method is used to predict the acoustic source. In order to consider the scattering and diffraction effects of the casing, Kirchhoff-Helmholtz boundary element method (BEM) is developed. The source of Kirchhoff-Helmholtz BEM is newly developed, so the sound field of the centrifugal fan can be obtained. A centrifugal impeller and wedge are used in the numerical calculation and the results are compared with the experimental data. Reasonable results are obtained not only for the peak frequencies but also for the amplitudes of the tonal sound. The radiated acoustic field shows the diffraction and scattering effect of the wedge.     

一种矢量水听器基阵模拟器的设计

矢量水听器的每个阵元同时测量声场中声矢量和质点振速的3个分量,相对于声压水听器阵来说,矢量阵获取声场中更多的信息。利用矢量阵所获得速度场的信息可以去除目标方位估计中的 180°模糊。模拟器可以模拟实际海洋环境中目标的辐射特性和噪声特性,应用模拟器可以有效地缩短声纳的研制周期。本文提出一种矢量水听器基阵模拟器的设计方案,该方案解决了矢量阵中宽带信号的90°移相问题、时延精确控制问题和宽带噪声的谱状控制问题。     

流体加载下加肋板结构的声辐射特性研究

本文采用有限元和边界元方法对加肋板结构的声辐射进行了计算分析,研究分析了加肋板结构的板厚、板面积、板边长比、肋骨惯性矩和间距、边界条件以及板材和流体介质等对结构辐射声功率的影响,得到了一些有意义的结论,从而为揭示加板结构声辐射规律,降低结构的声辐射提供了一定的依据。     

An approximate method to acoustic radiation problems: element radiation superposition method

An approximate method is brought forward to predict the acoustic pressure based on the surface velocity. It is named Element Radiation Superposition Method （ERSM）. The study finds that each element in Acoustic Transfer Vector （ATV） equals the acoustic pressure radiated by the corresponding surface element vibrating in unit velocity and other surface elements keep still, that is the acoustic pressure radiated by the corresponding baffled piston vibrating in unit velocity. So, it utilizes the acoustic pressure radiated by a baffled piston to establish the transfer relationship between the surface velocity and the acoustic pressure. The total acoustic pressure is obtained through summing up the products of the surface velocity and the transfer quantity. It adopts the regular baffle to fit the actual baffle in order to calculate the acoustic pressure radiated by the baffled piston. This approximate method has larger advantage in calculating speed and memory space than Boundary Element Method. Numerical simulations show that this approximate method is reasonable and feasible.     

Numerical analysis of immersed finite cylindrical shells using a coupled BEM/FEM and spatial spectrum approach

This study numerically analyzes submerged cylindrical shells using a coupled boundary element method (BEM) with finite element method (FEM) in conjunction with the wave number theory, in which the spatial Fourier transform of surface velocity for cylinders is directly related to pressure in a far field. The acoustic loading is formulated using a symmetric complex matrix derived from a boundary integral equation where the symmetry is based on an acoustic reciprocal principle for surface acoustics. In this formulation the acoustic loading matrix is a large acoustic element whose degree of freedom is connected to the normal displacement of the vibrating structures. The coupled BEM/FEM equation is a banded, symmetric matrix, and thus its bandwidth can be minimized using a proper algorithm. This formulation significantly increases numerical efficiency. The computed normal velocity is thus transformed to wave number representation to examine acoustic radiation. A finite plane cylindrical shell, without attached stiffeners, and a shell with internal ring stiffeners are chosen to demonstrate the present analysis procedure. The far field pressure computed directly from the integral equation and predicted by wave number theory correlates closely with increasing vibrating frequency. Meanwhile, the influences of the internal ring structures on acoustic radiation are examined using the wave number theory, which helps in understanding how internal structures influence radiated noise.     

Uncertainty of array-based measurement of radiated and absorbed sound intensity

Patch near-field acoustic holography (NAH) coupled with an array of sound intensity probes allows separating the sound field incident on a surface from the one radiated by the surface itself. Although the measurement principle has been successfully used to separate the noise source contribution from disturbing sources and/or noise reflections, the method accuracy has not been investigated in the literature. We describe the results of experiments meant to evaluate the uncertainty in the identification of noise radiated by vibrating panels with different absorption characteristics in presence of an incident acoustic radiation using the statistically optimized near-field acoustic holography. Measurement errors were evaluated through tests performed in controlled acoustic conditions. Results evidenced that the measurement uncertainty depends on the accuracy of the microphone array positioning and on the incident sound field. These conclusions were in agreement with the results obtained by simulations in the phase of instrument optimization.     

Analytical prediction of break-out noise from a reactive rectangular plenum with four flexible walls

This paper describes an analytical calculation of break-out noise from a rectangular plenum with four flexible walls by incorporating three-dimensional effects along with the acoustical and structural wave coupling phenomena. The breakout noise from rectangular plenums is important and the coupling between acoustic waves within the plenum and structural waves in the flexible plenum walls plays a critical role in prediction of the transverse transmission loss. The first step in breakout noise prediction is to calculate the inside plenum pressure field and the normal flexible plenum wall vibration by using an impedance-mobility approach, which results in a compact matrix formulation. In the impedance-mobility compact matrix (IMCM) approach, it is presumed that the coupled response can be described in terms of finite sets of the uncoupled acoustic subsystem and the structural subsystem. The flexible walls of the plenum are modeled as an unfolded plate to calculate natural frequencies and mode shapes of the uncoupled structural subsystem. The second step is to calculate the radiated sound power from the flexible walls using Kirchhoff-Helmholtz (KH) integral formulation. Analytical results are validated with finite element and boundary element (FEM-BEM) numerical models.     

A three-dimensional hybrid finite element/spectral analysis of noise radiation from turbofan inlets

This paper describes a new three-dimensional (3D) analysis of tonal noise radiated from non-axisymmetric turbofan inlets. The novelty of the method is in combining a standard finite element discretisation of the acoustic field in the axial and radial coordinates with a Fourier spectral representation in the circumferential direction. The boundary conditions at the farfield, fan face and acoustic liners are treated using the same spectral representation. The resulting set of discrete acoustic equations are solved employing the well-established BICGSTAB or QMR iterative algorithms and a very effective specialised preconditioner based on the axisymmetric mean geometry and flow field. Numerical examples demonstrate the suitability of the new method to engine configurations with realistic 3D features, such as relatively large degrees of asymmetry and spliced acoustic liners. The examples also illustrate the two advantages of the new method over a traditional 3D finite element approach. The new method requires a significantly smaller number of unknowns as relatively few circumferential Fourier modes in the spectral solution ensure an accurate field representation. Also, due to the effective preconditioner, the spectral linear solver benefits from stable iterations at a high rate of convergence.