认知不确定二维声场的响应特性数值分析

针对声学参数存在认知不确定性的问题,为实现认知不确定声场声压响应的预测。提出了解决二维认知不确定声场的有限元法(Evidence Theory-based Finite Element Method,ETFEM),引入证据理论,采用焦元和基本可信度的概念来描述认知不确定参数,基于摄动法的区间分析技术,推导了认知不确定声场声压响应的标准差、期望求解公式。为验证本文方法的可行性。以认知不确定参数下的二维管道声场模型和某轿车二维声腔模型为例进行了数值计算,对比离散随机变量得到认知不确定参数的声场分析结果和相应的随机声场所得分析结果,研究表明:该方法能够有效的处理认知不确定参数下的二维声场,为工程问题中噪声预测提供可靠的分析模型。      

概率盒框架下混合认知不确定声场的响应预测

针对含概率盒-证据混合认知不确定参数声场的响应预测问题,提出了一种概率盒框架下的改进区间蒙特卡洛方法。该方法首先将混合认知不确定参数转换为纯概率盒形式,然后结合有限元方法推导出混合认知不确定声场的盖根鲍尔多项式代理模型,再采用蒙特卡洛方法求解代理模型得到声压响应。以含概率盒-证据混合认知不确定参数的二维管道声场模型和卡车乘客舱声腔模型为例,计算结果表明混合认知不确定参数影响下的声压响应为概率盒形式,其包括声压响应极值和相应的概率信息,并且所提方法较常规混合离散方法效率更优,较基于一阶摄动法的区间蒙特卡洛方法准确性更高。研究结果表明:所提方法可以有效预测混合认知不确定声场的声压响应,并可进行声学性能的风险和保守估计。      

不确定声场分析的二阶区间摄动有限元法

针对一阶区间摄动有限元法在声场参数不确定程度增大时误差过大的缺陷,在二阶Taylor展开的基础上推导了声学二阶区间摄动有限元法,并将其应用于区间不确定声场的声压响应分析。该方法先对声学区间有限元方程的声压响应向量进行二阶Taylor展开,获取声压响应的二阶近似响应向量;再根据二次函数极值定理获得声压响应向量的上下界。二维管道声场与轿车声腔模型的数值分析算例表明,与一阶区间摄动有限元法相比,二阶区间摄动有限元法有效提高了计算精度。因此二阶区间摄动有限元适合不确定度更大的区间不确定声场声压响应分析,具有良好的工程应用前景。      

Numerical analysis of the 2D acoustic field with fuzzy-random parameters

Regard for the fuzziness and the randomness in some acoustic fields,a method for the numerical analysis of the 2D acoustic field with Fuzzy-Random parameters was proposed based on the equivalent conversion of information entropy.In the proposed method,a fuzzyrandom acoustic field was treated as a pure fuzzy acoustic field or a pure random acoustic field by transforming all the variables into fuzzy variables or random variables.Perturbation finite element methods for analyzing the two-dimensional acoustic fuzzy and random field are deduced.The sound pressure response of a 2D acoustic tube and the 2D acoustic cavity of a car with fuzzy-random parameters were analyzed by the proposed method and the Monte Carlo method,the results show that the proposed method can be well applied to the numerical analysis of the 2D acoustic field with fuzzy-random parameters,and has good prospect of engineering application.     

模糊随机参数二维声场数值计算

解决声场参数同时具有模糊性和随机性的问题,实现模糊随机声场声压响应的预测,引入了信息熵理论,利用信息熵的等效转换,将模糊随机声场转化为纯随机声场或者纯模糊声场进行求解,推导了基于摄动法的二维随机声场和模糊声场的有限元计算公式。以模糊随机参数下的二维管道声场模型和某轿车二维声腔模型为例进行了数值计算,所得结果与蒙特卡洛法(Monte Carlo Method)所预测声压变化范围基本一致,同时,转化为纯随机声场和纯模糊声场所求得声压响应变化范围也基本一致,说明了本文方法计算结果的准确性。因此本文方法能很好地应用于模糊随机参数下二维声场的预测,具有重要的工程应用价值。      

Modified sub-interval perturbation finite element method for 2D acoustic field prediction with large uncertain-but-bounded parameters

Based on the sub-interval perturbation analysis, a modified sub-interval perturbation finite element method is proposed to determine the bounds of sound pressure in the 2D acoustic field with large uncertain-but-bounded parameters. In the proposed method, the inversion of the invertible sub-interval dynamic stiffness matrix is approximated by a modified approximate interval-value Sherman–Morrison–Woodbury formula to overcome the drawbacks arising from the dependency phenomenon of parameters and the unpredictable effect of neglecting the higher order terms in Neumann series. The modified sub-interval perturbation procedures are implemented in a numerical finite element framework. Numerical results on a 2D acoustic tube and a 2D acoustic cavity of a car with large uncertain-but-bounded parameters evidence the remarkable accuracy and effectiveness of the proposed method.     

The rapid uncertainty prediction of the ocean-acoustic coupled model

Focusing on the rapid prediction of acoustic field uncertainty in environment with temporal and spatial sound speed perturbation, evolvement of sound speed structure over time is predicted based on the ocean-acoustic coupled model to obtain the uncertainty distribution of the vertical structure of sound speed. Further, a method combining the arbitrary polynomial chaos expansion with the empirical orthogonal function is proposed to reduce the dimensionality of uncertain parameters and to obtain the uncertainty distribution of the acoustic field. Simulations have shown that the computational complexity can be reduced by 2 orders of magnitude compared to the conventional polynomial chaos expansion while ensures the same precision.Moreover, the computational complexity is not influenced by the complexity of the sound speed profile. The acoustic field and uncertainty predicted in uncertain environment by proposed method also have been tested with the experimental data.     

海洋—声学耦合过程不确定性快速预报

针对海洋声速空间非均匀和时变条件下的声场不确定性快速预报问题,根据海洋-声学耦合模式预报声速场时空变化过程,获取声速垂直结构不确定性分布规律,提出经验正交函数-随机多项式展开方法,以降低不确定参数维度,得到声场不确定性分布.数值计算表明该方法在保证同等精度的同时可大幅减少计算量,相比常规随机多项式展开方法计算效率可提高2个数量级,且计算量不随声速剖面复杂程度变化而改变.经海上实验验证,结果表明该方法可实现环境不确定性条件下的声场及其不确定性的快速预报。      

浅海不确定声场的随机多项式展开法研究

为获得求解浅海不确定声场的普适模型,建立了随机多项式展开法与Helmholtz方程的非嵌入式耦合模型,其间运用概率配点法求解多项式系数。针对仅当海水深度不确定时的Pekeris波导、声速剖面和海深均不确定时的Pekeris波导以及下限深度不确定温跃层等几种情形,计算了传播损失概率密度分布。结果表明所建模型对声场计算模型普适性强,计算精度和计算效率高,可用于研究含多个不确定环境参数、声速剖面复杂的浅海环境中声传播的不确定性。      

Numerical approach for quantification of epistemic uncertainty

In the field of uncertainty quantification, uncertainty in the governing equations may assume two forms: aleatory uncertainty and epistemic uncertainty. Aleatory uncertainty can be characterised by known probability distributions whilst epistemic uncertainty arises from a lack of knowledge of probabilistic information. While extensive research efforts have been devoted to the numerical treatment of aleatory uncertainty, little attention has been given to the quantification of epistemic uncertainty. In this paper, we propose a numerical framework for quantification of epistemic uncertainty. The proposed methodology does not require any probabilistic information on uncertain input parameters. The method only necessitates an estimate of the range of the uncertain variables that encapsulates the true range of the input variables with overwhelming probability. To quantify the epistemic uncertainty, we solve an encapsulation problem, which is a solution to the original governing equations defined on the estimated range of the input variables. We discuss solution strategies for solving the encapsulation problem and the sufficient conditions under which the numerical solution can serve as a good estimator for capturing the effects of the epistemic uncertainty. In the case where probability distributions of the epistemic variables become known a posteriori, we can use the information to post-process the solution and evaluate solution statistics. Convergence results are also established for such cases, along with strategies for dealing with mixed aleatory and epistemic uncertainty. Several numerical examples are presented to demonstrate the procedure and properties of the proposed methodology.     

Multiple-array passive acoustic source localization in urban environments

In many situations of interest, obstacles to acoustic wave propagation such as terrain or buildings exist that provide unique challenges to localization. These obstacles introduce multiple propagation paths, reflections, and diffraction into the propagation. In this paper, matched field processing is proposed as an effective method of acoustic localization in a two dimensional scattering environment. Numerical techniques can be used to model complex propagation in a space where analytical solutions are not feasible. Realistically, there is always some uncertainty in model parameters that in turn can adversely affect localization ability. In particular, uncertainty in array location, sound speed, and various parameters affecting inter-array coherence only are investigated. A spatially distributed, multiarray network is shown to mitigate the effects of uncertainty. Multiarray inverse filter processing techniques are evaluated through perturbation of uncertain model parameters. These techniques are more accurate and flexible to implement than other matched field processing methods such as time reversal.     

基于声压场描述的扩散声场多自由度互易原理研究

应用于复杂结构中频声振分析的扩散场多自由度互易原理采用位移变量描述系统,实质为弹性波场互易原理,应用于声波场时会造成模型自由度数不必要的增加。建立基于声压描述的扩散声场受挡模型,利用声辐射模态描述扩散声场中结构的表面受挡声压;据此提出基于声压描述的扩散声场多自由度互易原理,发现扩散声场中结构表面受挡声压的互谱矩阵与该结构在自由空间中振动辐射声波的声阻矩阵成正比。该互易原理与传统的单自由度互易原理表达形式相似,但适用于任意自由度结构。该互易原理可用于扩散声场中复杂结构的表面受挡声压的自谱及相关分析,仿真研究表明当边界元网格尺寸小于声波波长的1/6(线性单元)或1/3(二次单元)时,数值解与理论解完全吻合。      

Hybrid uncertainty propagation of coupled structural–acoustic system with large fuzzy and interval parameters

Based on the finite element framework and uncertainty analysis theory, this paper proposes a first-order subinterval perturbation finite element method (FSPFEM) and a modified subinterval perturbation finite element method (MSPFEM) to solve the uncertain structural–acoustic problem with large fuzzy and interval parameters. Fuzzy variables are used to represent the subjective uncertainties associated with the expert opinions; whereas, interval variables are adopted to quantify the objective uncertainties with limited information. By using the level-cut strategy and subinterval methodology, the original large fuzzy and interval parameters are decomposed into several subintervals with small uncertainty level. In both FSPFEM and MSPFEM, the subinterval matrix and vector are expanded by the Taylor series. The inversion of subinterval matrix in FSPFEM is approximated by the first-order Neumann series, while the modified Neumann series with higher order terms is adopted in MSPFEM. The eventual fuzzy interval frequency responses are reconstructed by the interval union operation and fuzzy decomposition theorem. A numerical example evidences the remarkable accuracy and effectiveness of the proposed methods to solve engineering structural–acoustic problems with hybrid uncertain parameters.     

多项式混沌展开的浅海水声环境参数敏感度分析

水下声场是海洋环境参数的非线性函数,因此环境参数的不确定性必然会引起水下声场预测的不确定性,从而导致与声场预测相关的声呐探测和水声通信等设备性能的下降.该文将Sobol敏感度指数的计算与多项式混沌展开方法相结合,将敏感度指数表达为环境条件的函数,利用Q范数约束的双曲截断方案来减少多项式项数,有效地分析了设定的"浅海负梯...     

Pekeris waveguide comparisons of methods for predicting acoustic field amplitude uncertainty caused by a spatially uniform environmental uncertainty (L)

Acoustic field calculations in underwater environments are often uncertain because the environmental parameters required for such calculations are uncertain. This letter compares the accuracy of direct simulations, the field shifting approximation, and polynomial chaos expansions for predicting acoustic amplitude uncertainty in 100-m-deep Pekeris waveguides having spatially uniform uncertain water-column sound speed. When this sound speed is Gaussian-distributed with a standard deviation of 1 m/s, direct simulations and polynomial chaos expansions, based on 21 field calculations, are more accurate than the field shifting approximation, based on two field calculations. This ranking reverses as the sound-speed standard deviation increases to 20 m/s.     

Proper orthogonal decomposition and cluster weighted modeling for sensitivity analysis of sound propagation in the atmospheric surface layer

Outdoor sound propagation predictions are compromised by uncertainty and error in the atmosphere and terrain representations, and sometimes also by simplified or incorrect physics. A model's predictive power, i.e., its accurate representation of the sound propagation, cannot be assessed without first quantifying the ensemble sound pressure variability and sensitivity to uncertainties in the model's governing parameters. This paper describes fundamental steps toward this goal for a single-frequency point source. The atmospheric surface layer is represented through Monin-Obukhov similarity theory and the acoustic ground properties with a relaxation model. Sound propagation is predicted with the parabolic equation method. Governing parameters are modeled as independent random variables across physically reasonable ranges. Latin hypercube sampling and proper orthogonal decomposition (POD) are employed in conjunction with cluster-weighted models to develop compact representations of the sound pressure random field. Full-field sensitivity of the sound pressure field is computed via the sensitivities of the POD mode coefficients to the system parameters. Ensemble statistics of the full-field sensitivities are computed to illustrate their relative importance at every down range location. The central role of sensitivity analysis in uncertainty quantification of outdoor sound propagation is discussed and pitfalls of sampling-based sensitivity analysis for outdoor sound propagation are described.     

Distinguishing between the interior pressures induced from two independent sources within a room using the probabilistic approach

In this paper, a probabilistic approach is introduced and used to distinguish between the interior pressures that are induced from two independent sound sources within a rectangular room. One source is a vibrating wall of the room and the other is an interior point source. The model is set up using the modal analysis method. In the probabilistic method, the probability densities of the uncertain acoustic model parameters are computed. The values of highest probability density are identified and assigned to the uncertain parameters in the acoustic model. The contributions of the sound sources to the total pressure are then obtained from the acoustic model input with the identified parameters. The effects of the modal truncations, the modeling errors, and the measurement locations on the accuracies of the identification process are studied.     

基于时间反转法的相控换能器声场的仿真研究

相控换能器具有焦距可调的优势。本文以82阵元相控换能器建立的3D数值仿真模型为例,基于时间反转法提取阵元的激励信号,利用时域有限差分(FDTD)法对Westervelt声波非线性传播方程进行声场数值仿真,研究不同阵元分布、偏离声轴的距离、设定焦距大小对形成声场的影响,可调控范围及其消除旁瓣方法。研究结果表明,随机分布相控阵可明显降低声场中的旁瓣;随着偏离声轴距离的增加,主瓣声压幅值逐渐减小,旁瓣与主瓣的最大声强比值r逐渐增大,且沿声轴的可调控范围逐渐减小;随声轴方向上设定焦距的增加,主瓣声压幅值先增大后减小,r值先减小后增大;基于时间反转法的高声压旁瓣消除法可在一定程度上扩大相控阵声场的可调控范围。     

Vibro-acoustic behaviors of an elastically restrained double-panel structure with an acoustic cavity of arbitrary boundary impedance

An analytical study on the vibro-acoustic behaviors of a double-panel structure with an acoustic cavity is presented. Unlike the existing studies, a structural–acoustic coupling model of an elastically restrained double-panel structure with an acoustic cavity having arbitrary impedance on sidewalls around the cavity is developed in which the two dimensional (2D) and three dimensional (3D) modified Fourier series are used to represent the displacement of the panels and the sound pressure inside the cavity, respectively. The unknown expansions coefficients are treated as the generalized coordinates and the Rayleigh–Ritz method is employed to determine displacement and sound pressure solutions based on the energy expressions for the coupled structural–acoustic system. The effectiveness and accuracy of the present model is validated by numerical example and comparison with finite element method (FEM) and existing analytical method, with good agreement achieved. The influence of key parameters on the vibro-acoustic behaviors and sound transmission of the double-panel structure is investigated, including: cavity thickness, boundary conditions, sidewall impedance, and the acoustic medium in the cavity.     

Acoustic characteristics of pulse detonation engine with ellipsoidal reflector

Acoustic characteristics of a pulse detonation engine(PDE) with and without an ellipsoidal reflector are numerically and experimentally investigated. A two-dimensional(2 D) non-splitting unstructured triangular mesh Euler solver based on the space-time conservation element and solution element(CE/SE) method is employed to simulate the flow field of a PDE.The numerical results clearly demonstrate the external flow field of the PDE. The effect of an ellipsoidal reflector on the flow field characteristic near the PDE exit is investigated. The formation process of reflected shock wave and reflected jet shock are reported in detail. An acoustic measurement system is established for the PDE acoustic testing. The experimental results show that the ellipsoidal reflector changes the sound waveform and directivity of PDE sound. The reflected shock wave and reflected jet shock result in two more positive pressure peaks in the sound waveform. The ellipsoidal reflector changes the directivity of PDE sound from 20 to 0. It is found that the peak sound pressure level(PSPL) and overall sound pressure level(OASPL) each obtain an increment when the PDE is installed with a reflector. The maximum relative increase ratio of PSPL and OASPL are obtained at the focus point F2, whose values are 6.1% and 6.84% respectively. The results of the duration of the PDE sound indicate that the reflecting and focusing wave generated by the reflector result in the increment of A duration and B duration before and near focus point F2. Results show that the ellipsoidal reflector has a great influence on the acoustic characteristic of PDE sound. The research is helpful for understanding the influence of an ellipsoidal reflector on the formation and propagation process of PDE sound.