自由液面和刚性壁面对结构振动声辐射的影响

在边界元理论的基础上,建立了半无限流体域中、特别是与无限大平面接触的结构流固耦合振动方程,计算了方箱的固有频率、振动响应及声辐射,探讨了自由液面和刚性壁面对结构的固有频率、振动响应和有关声学物理量的影响,并讨论了板厚、结构阻尼、结构与无限大平面的距离等有关因素。研究结果表明:自由液面和刚性壁面对结构的固有频率、振动响应及声辐射的影响是不容忽视的。     

进一步改进边界元方法以克服振动声辐射计算中解的非唯一性

章利用基于三次B样条插值的边界元方法，对振动体外部声辐射问题进行了研究，对CHIEF法及其改进方法作了进一步的改进，提出在加权余量意义下，通过把内部Helmholtz积分方程与其对内点坐标取导后的方程式作线性叠加，在域外构作的一个小体积块上进行积分以形成补充方程，经与表面Helmholtz积分方程相结合，来求解任意频率下的声辐射问题，并以脉动球和摆动球作为算例，说明本提出的方法能够有效地克服在特殊频率处解的非唯一性问题。     

半球封闭圆柱体表面绕射声场研究

基于一致性几何绕射理论,对半球封闭圆柱体上半球表面的高频绕射声场进行了数值计算与分析。获得了单个点声源作用下,半球表面声场的声压分布图,其中包括一些干涉条纹,干涉条纹的区域及形状与点声源的位置和频率有关;在多个点声源作用下,发现在接近半球顶点的附近,多个声波经多途传播后在该区域的叠加声压由衰减转为增强的现象。     

改进的非负声强和有用声强识别表面声辐射*

识别结构表面对远场声辐射有实质贡献的区域,对控制远场声辐射有重要意义。该文以非负声强法为基础得到改进的声强计算方法,这种方法比非负声强法计算量小,能更好地识别出表面对远场有贡献的区域。该文对有用声强法从特征值和截断准则的角度进行了一些改进。最后,以四边简支薄板为例,分别采用这三种方法计算几种不同频率下的噪声辐射区域。数值计算结果表明,采用该文提出的方法可以较好地识别出结构表面对远场声辐射有实质贡献的区域,而且计算量相对较小。     

压电晶体表面激发的二维声表面波Green函数

最近,我们研究了在压电晶体半空间表面激发的二维广义Green函数。在本文将仔细讨论广义Green函数的表面波部分。用稳相法给出声表面波远场渐近解。但是由于这样的渐近解在接近慢度曲线的拐点处失效,我们发展了用变型Airy函数表示的普适渐近解,从而解决了在拐点附近的解的衔接问题。对于Y切LINbO_3,进行了声表面波Green函数分量和功率流密度的数值计算,结果表明在Y切LiNbO_3的Z方向上出现聚焦现象.     

强弱耦合板-腔系统的声辐射模态计算与分析EI北大核心CSCD

针对板-腔耦合系统的声辐射模态(ARM)计算问题,提出了一种基于能量原理的声辐射模态计算方法,该方法从能量原理的动力学方程构建起声压模态幅值和结构模态幅值的关系,通过将声势能表示为结构模态幅值向量的二次型形式,得到板-腔耦合系统的声辐射模态,弥补了前人理论在解决声腔为阻抗壁面和结构-声为强耦合条件时的不足。通过数值算例验证了本文计算方法的正确性和有效性,在此基础上分析了壁面和结构-声耦合条件变化对声辐射模态特性的影响。结果表明:声辐射模态辐射效率曲线会在声腔模态频率处产生峰值,阻抗壁面的引入会降低声辐射模态辐射效率在峰值处的幅值,并且阻抗值越小,幅值衰减效应越明显,具体表现为声势能曲线在辐射效率峰值频率处幅值会下降;强耦合条件下低频段声势能响应主要由弹性板结构模态激发,响应峰值密度更高,幅值更低。低频同频宽的声辐射模态辐射效率峰值数更少,峰值频率更高。     

Gauss声束对离轴椭圆柱的声辐射力矩

利用部分波展开法求解得到了Gauss声束入射下刚性和非刚性椭圆柱的声散射系数,推导了一般情况下的声辐射力矩表达式.在此基础上,通过一系列数值仿真详细分析了离轴距离、入射角度和束腰半径对声辐射力矩的影响.结果表明:正向与负向声辐射力矩均可以在一定条件下存在;低频情况下刚性椭圆柱比非刚性椭圆柱更容易产生较强的声辐射力矩;特定频率的入射声场可以激发出非刚性椭圆柱不同阶的共振散射模式,因而非刚性椭圆柱的声辐射力矩峰值与频率的关系更密切;增加束腰半径有利于扩大散射截面,进而增加椭圆柱的声辐射力矩.该研究结果预期可以为利用声辐射力矩实现粒子的可控旋转和流体黏度的反演提供一定的理论指导.     

激光空泡刚性半球面内运动

采用甚高速照相技术与建立激光空泡在刚性半球壁面内的运动模型相结合的方法,确定了激光空泡在刚性半球面壁内的运动特性与无量纲距离的关系,提出了最佳无量纲距离概念。结果表明:半球反射面的半径与激光空泡最大半径之比小于1.1时,激光空泡在第1次膨胀时就会产生严重的变形并弹出半球面,并产生空化泡和空蚀,它们均会严重影响激光声的传播。该比值在1.1~3.3时,激光空泡将在第3次收缩之前接触半球面,容易对壁面造成空蚀。在该比值大于3.3的情况下,激光空泡在第3次收缩之前不会接触半球面,对激光声的传播和反射特性影响较小。如果考虑把空泡第1次溃灭时产生的激光声的声学中心控制在击穿点时,需要把该比值控制在5以上。     

具有弹性障板的有限长充液管道的管口声辐射

求解了下端自由、上端带弹性障板的有限长充液管道,在下端受到流体脉动激励时,上端管口向水下半无限空间的声辐射。求解过程中,管道部分采用行波法求解,其中管壁与管内流体的耦合归结到求解充液管道的色散曲线;弹性障板的振动部分则采用模态法来求解,其中弹性障板及管口与外界流体间的声振耦合最终归结到求解模态的辐射阻抗。数值研究表明:管口辐射的流体声和弹性障板辐射的结构声对系统的声辐射均有较大的贡献。     

积分方程法与求解谐振频率的声散射

表面Ｈｅｌｍｈｏｌｔｚ积分公式可以有效地求解物体的辐射声场和散射声场，但并不是在任何频率下都能得到满意的结果。当波数等于或接近目标内部问题的特征值时，将产生非唯一解，严重影响求解的准确性。本文根据混合Ｈｅｌｍｈｏｌｔｚ积分公式，用最小平方正交法有效地求解谐振频率的声散射。文中以圆柱和椭圆柱的声散射为例进行了计算，并讨论了目标内附加计算点的选取问题。     

Solving the hypersingular boundary integral equation in three-dimensional acoustics using a regularization relationship

Regularization of the hypersingular integral in the normal derivative of the conventional Helmholtz integral equation through a double surface integral method or regularization relationship has been studied. By introducing the new concept of discretized operator matrix, evaluation of the double surface integrals is reduced to calculate the product of two discretized operator matrices. Such a treatment greatly improves the computational efficiency. As the number of frequencies to be computed increases, the computational cost of solving the composite Helmholtz integral equation is comparable to that of solving the conventional Helmholtz integral equation. In this paper, the detailed formulation of the proposed regularization method is presented. The computational efficiency and accuracy of the regularization method are demonstrated for a general class of acoustic radiation and scattering problems. The radiation of a pulsating sphere, an oscillating sphere, and a rigid sphere insonified by a plane acoustic wave are solved using the new method with curvilinear quadrilateral isoparametric elements. It is found that the numerical results rapidly converge to the corresponding analytical solutions as finer meshes are applied.     

基于Burton-Miller方程的轴对称结构声学边界元方法

提出了一种求取轴对称结构任意边界条件下声辐射特性的边界元方法。采用Burton和Miller改进型公式将高阶奇异项转化为弱奇异项之和,保证声辐射参数的唯一性,且计算简单精确。将结构表面声压与振速按照旋转轴角度进行Fourier级数展开,利用级数的正交性建立各项待定系数的求解公式;然后转化格林函数的法向偏导为切向偏导,方便直接计算各项积分,并将面积分公式表示为沿结构边界的线积分和沿旋转角度的积分;进一步采用二次等参单元离散结构边界线,建立声压与振速的关系矩阵,从而确定结构声辐射参数。以脉动球源和横向振动球源为例计算,与解析解和传统边界元法结果作对比,说明该方法的有效精确性。     

弯曲圆盘换能器与镜像虚源的互辐射作用

建立了弯曲圆盘换能器镜像虚源等效模型,利用脉动球源互作用原理对弯曲圆盘与其镜像虚源的互辐射作用进行理论分析,给出了互辐射阻抗及系统谐振频率的数学表达。提出了虚源互作用的低频换能器设计思想,将刚性反射板引入弯曲圆盘换能器临近辐射面的声场中,通过理论分析、有限元模拟和样机实验研究了低频换能器谐振频率与主要结构参数之间的关系。结果表明,弯曲圆盘与镜像虚源间的互辐射作用可以有效降低换能器的谐振频率,当反射板直径与弯曲圆盘辐射面直径相当时,谐振频率可降低至其自身谐振频率的50%以下;当反射板直径为弯曲圆盘辐射面直径2倍时,谐振频率可降低至37%。      

Sound wave scattering by a spherical scatterer located near an ice surface

An echo signal is simulated, which is reflected from a spherical scatterer located near an ice surface. The homogeneous water medium in which the scatterer is located is assumed semi-infinite. For the scattering coefficients of the sphere, asymptotic formulas are obtained by the saddle point method, which can be used for sufficiently large distances between the source emitting a spherical wave and the scatterer. For the occurring branch cut integrals using the steepest descent method, asymptotic expressions are also obtained. Numerical results are obtained for an acoustically rigid sphere and an ice sphere. The density of the ice medium and speed of longitudinal waves in it coincide with the analogous parameters of the ice cover. In a wide frequency range of 8–12 kHz, echo signals are compared that have been calculated for two models of media: a water half-space bordering an ice half-space and an ice-covered homogeneous waveguide with a fluid bottom under the assumption that the source placed in the water layer is directional. It is shown that in a large distance interval between the source and the spherical scatterer, the half-space model sufficiently accurately describes the echo signal while substantially reducing calculation time (by approximately a factor of 10 for the waveguide with a depth of 200 m and a sandy bottom considered in the paper).     

An assessment of the high frequency boundary element and Rayleigh integral approximations

This paper revisits the popular Rayleigh integral approximation and also considers a second approximation, the high frequency boundary element method, which is similar to the Rayleigh integral. The Rayleigh integral approximation under consideration is enhanced so that only the elements visible to a particular point in the field are used to calculate the sound pressure at that point. It is demonstrated how both the Rayleigh integral and high frequency boundary element method are approximations to the boundary integral equation so that similarities between the two methods are recognized. Several test cases were conducted in order to assess and compare both methods. The first set of test cases was the pulsating and oscillating sphere. Both methods were then compared on more applied examples including a running engine, construction cab, and transmission housing. It was concluded that though both methods can reliably predict the sound power for some problems, the high frequency boundary element method is the more robust.     

Dictums for problem solving and approximation in mathematical acoustics: examples involving low-frequency vibration and radiation

A sequence of dictums for mathematical acoustics is given representing opinions intended to be regarded as authoritative, but not necessarily universally agreed upon. The dictums are presented in the context of the detailed solution for a class of problems involving the forced vibration of a long cylinder protruding half-way into a half-space bounded by a compliant surface (impedance boundary) characterized by a spring constant. One limiting case corresponds to a cylinder vibrating within an infinite rigid baffle, and another limiting case corresponds to a vibrating cylinder on the compliant surface of an incompressible fluid. The second limiting case is identified as analogous to that of a floating half-submerged cylinder whose vibrations cause water waves to propagate over the surface. Attention is focused on vibrations at very low frequencies. Difficulties with insuring a causal solution are pointed out and dictums are given as to how one overcomes such difficulties. Various approximation techniques are described. The derivations involve application of the theory of complex variables and the method of matched asymptotic expansions, and the results include the apparent entrained mass in the near field of the cylinder and the radiation resistance per unit length experienced by the vibrating cylinder.     

Indentation on one-dimensional hexagonal quasicrystals: general theory and complete exact solutions

This paper presents a general account of the indentation responses of a one-dimensional hexagonal quasicrystal half-space pressed by an axisymmetric rigid punch. Based on Green's functions of the half-space subjected to point sources on the surface, the mixed boundary value problem is transformed to integral equations and solved exactly using the results of the potential theory method. Explicit expressions for the generalised pressures and indentation forces are derived for three common indenters (cylinder, cone and approximate sphere) in a systematic manner. For conical and spherical indenters, relations between the contact radius and indentation loads are determined. The coupling phonon–phason fields in the half-space under indentation are accurately expressed in terms of elementary functions. Numerical calculations are performed and discussions on related physical phenomena are given. The present exact solutions can serve as benchmarks for approximate or numerical analyses and can guide the experimental characterisation of material properties of quasicrystals.     

利用声辐射模态对声场进行重构

声辐射模态是定义在振动表面上的一组互相独立的基函数,描绘了振动表面的多极子辐射模式。本文利用声辐射模态对平板辐射的声场进行重构，并同时就频率和模态数对误差的影响进行了分析。针对复杂振动表面幅射的声场，重构时需正则化处理。数值计算结果和对应的理论值达到很好的吻合。     

Prediction of sound level at high-frequency bands by means of a simplified boundary element method

A simplified boundary element method (BEM) for dealing with high-frequency sound is proposed. The boundary integral equation is modified into a quadratic form to enable the prediction of sound levels in the one-third octave band analysis. Monopole and dipole source terms in the conventional BEM are transformed into the auto- and cross-spectra of two vibrating sources, in which the cross-spectra are eventually neglected by assuming that the correlation coefficients involved are negligible. The present method is compared with the Rayleigh integral for calculating the sound pressure radiated from a baffled panel, in terms of the application limit. The characteristic length of the boundary element and the applicable frequency range can be determined by the lower limit value of the correlation coefficient. As a test example, the field pressure radiated from a partially vibrating sphere is predicted and the resultant trend is in good agreement with the analytic solution as far as the related correlation coefficient satisfies the assumption. The overdetermination process for overcoming nonuniqueness in exterior radiation problems is unnecessary in the present method because phase information can be ignored. The results of the calculation show that the proposed method is acceptable for solving the exterior radiation problem at a high-frequency range in a timely manner.     

VIRTUAL REALITY OF SOUND GENERATED FROM VIBRATING STRUCTURES

The advancement of virtual reality (VR) technology in cyberspace is amazing, but its development is mainly concentrated on the visual part. In this paper, the development of VR technology to produce sound based on the exact physics is studied. Our main concern is on the sound generated from vibrating structures. This may be useful, for example, in apprehending sound field characteristics of an aircraft cabin in design stage.To calculate sound pressure from curved surface of a structure, a new integration scheme is developed in boundary element method. Several example problems are solved to confirm our integration scheme. The pressure distributions on a uniformly driven sphere and cylinders are computed and compared with analytic solutions, and radiation efficiency of a vibrating plate under one-dimensional flow is also calculated. Also, to realize sound through computer simulation, two concepts, “structure-oriented analysis” and “human-oriented analysis”, are proposed. Using these concepts, virtual sound field of an aircraft cabin is created.