小尺度阻尼边界封闭空间声传递函数的有限元素法计算

如何求解阻尼边界封闭空间中声源点到接收点的低频声传递函数已成为目前小尺度封闭空间可听化技术研究的关键技术，能处理任意形状及复杂边界条件的有限元素法可作为求解该问题的适合方法，以室内声声有源Helmholtz方程及其相应边界方程为基础，本文推导出了用于小尺度阻尼边界封闭空间声传递函数的有限元素求解方法，并编制了相应的计算机程序，在算例中，首先通过与模态叠加法计算结果进行比较，验证了该方法的正确性。最后计算了某型车体内腔中任意两点间声传递函数。     

适用于任意形状封闭空间的声场预测相干模型

针对现有几何声学的方法对封闭空间内声场进行预测时在中低频段出现较大误差的问题,该文提出一种近似圆锥声束追踪法和相干反射场理论相结合的声场预测新模型。在近似圆锥声束追踪法基础上,考虑声束轴线在边界多次反射时声压和相位的改变,最后计算不同声波之间的干涉效应,建立一种适用于任意形状封闭空间的声场预测相干模型。利用该模型对某一矩形封闭空间进行声场预测,通过对边界元法、Raynoise软件相干和非相干算法的预测结果和本模型的数值模拟结果对比。结果表明,文中提出的方法和边界元法的计算结果在中低频段非常吻合,两者的计算结果平均绝对误差为1.1 d B。本模型在中低频率下与同样考虑了相位的Raynoise相干算法相比有更好的准确性,在较高频率上,本模型计算结果与Raynoise相干算法计算结果非常吻合。     

飞机机体表面声压及舱内降噪优化设计

本文以某型客机为研究对象,从飞行试验数据分析和声学建模两方面研究机体表面声压分布及其对舱内壁板近场辐射声压的影响。首先根据试飞数据分析了机体表面声压分布,然后利用统计能量法建立飞机客舱中后段的声学模型,以试飞数据作为声源输入,研究机体表面声压分布对客舱内部壁板附近声压分布的影响,并在此基础上提出优化设计方案,通过模型验证优化方案的有效性。试飞数据表明：机体表面声压在后应急门前方、靠近地板处最大;巡航速度升高,声压级较大区域的面积随之增加;巡航高度和发动机N1N2频率变化对机体表面声压级分布无明显影响。仿真数据表明：仅蒙皮结构无法有效降低客舱噪声;对声学降噪包进行优化能增加壁板隔声量,降低舱内声压。     

ARDE算法及其在三维叶栅气动优化设计中的应用

本文提出了一种自适应差分进化算法(ARDE)。通过函数试验,分析比较了该算法、标准遗传算法及差分进化算法的优化性能,证明该算法具有良好的收敛速度和全局寻优能力,满足气动优化设计的要求。耦合该算法与曲面造型方法以及CFD求解技术,本文发展了一种适用于叶轮机械三维气动优化设计的全局自动气动优化算法。利用该算法对小展弦比后加载静叶栅和干式TRT动叶栅进行了气动优化设计,结果表明本文提出的自动气动优化算法具有良好的优化性能和应用前景。     

基于声线跟踪法对闭空间声场预测的接收球半径计算

接收球半径是采用声线跟踪法进行声场预测的一个重要参数。以前所用到的接收球半径一般是通过经验估计得到的,缺乏详细的理论分析。本文从声波传播机理角度分析了根据声线密度来确定接收球半径的原因,并推导了矩形闭空间中声线密度和接收球半径的计算方法。声线密度可以通过原始声线数目,声场空间体积与形状,边界吸声系数来确定。在一给定闭空间里,声线密度可看作是均匀分布的,所以接收球半径与空间位置无关,可看作是一个常数。但对不同的声场空间来说,由于空间体积、形状和吸声系数的变化,声线密度是不同的,因此接收球半径也会不同。声线密度越大,接收球半径越小;声线密度越小,接收球半径越大。实验表明,所提出的接收球模型能用来较准确地预测闭空间里的声压级和混响时间等声学参数。     

某型艇近场噪声和自噪声数值计算

采用结构有限元软件NASTRAN和声学边界元软件SYSNOISE,对某型艇振动与声学特性数值计算问题进行了研究。通过基于全船三维有限元分析模型的频响分析计算,获得该船结构振动响应特性参数。以振动频响分析计算结果作为声学边界条件,考虑球面点噪声源的影响,建立了全船三维边界元声学分析模型;并采用间接边界元法,对某型艇近场和自噪声声学特性进行了计算。所得计算结果与该船实测数据进行了对比,证明本文计算方法是可行的。     

室内声场无网格数值算法性能影响因素分析

基于室内无网格数值算法的理论和实验研究,本文进一步利用控制变量法,结合计算效率对影响计算精度的因素进行了全面分析,包括节点、高斯点的分布方式及其密度,以及构造无网格形函数过程中的基函数类型、权函数类型,并通过对一类基准问题的大量数据实验获得了适用于小尺度封闭空间的最佳支持域半径范围。基于有关结论对一个实际机舱进行了建模分析,将声压级计算值与测量值对比,验证了这些结论的正确性。     

面向多目标动态优化设计的混合遗传算法

多目标动态优化设计问题中，特征值或结构响应等目标函数是一个复杂的非线性隐函数，可行性域态不良、解空间一般是多连通的，缺乏高效率的分析和优化算法。此外，在进行动态设计的同时，往往要求设计结果能够使多个相互矛盾的目标如重量、容积等得到优化。因此，动态设计问题实际上是一个具有高度非线性、解空间复杂的多目标动态优化问题，这类问题既具有一般多目标优化问题的复杂性，又具有动态设计问题固有的困难性。     

大规模声学边界元法的GPU并行计算

提出一种大规模声学边界元法的高效率、高精度GPU并行计算方法.基于Burton-Miller边界积分方程,推导适于GPU的并行计算格式并实现了传统边界元法的GPU加速算法.为提高原型算法的效率,研究GPU数据缓存优化方法.由于GPU的双精度浮点运算能力较低,为了降低数值误差,研究基于单精度浮点运算实现的doublesingle精度算法.数值算例表明,改进的算法实现了最高89.8%的GPU使用效率,且数值精度与直接使用双精度数相当,而计算时间仅为其1/28,显存消耗也仅为其一半.该方法可在普通PC机(8GB内存,NVIDIA Ge Force 660 Ti显卡)上快速完成自由度超过300万的大规模声学边界元分析,计算速度和内存消耗均优于快速边界元法.     

多翼离心风机气动噪声的数值分析

多翼离心风机叶片短、流道窄,叶轮出口流速分布不均,引起叶轮与蜗壳干涉作用加剧。本文探讨流场与声场非定常耦合机理,根据声类比理论分析其偶极子声源产生的气动噪声。利用直接边界元声学求解方法建立以蜗壳为界的内外声学模型,分析蜗壳对声传播的散射作用,内部噪声通过蜗壳的进出口传播到风机外部。结果表明:从监测点声压级频谱及A计权声压级分布观察,声压级分布在低频段呈宽频分布,在基频与其倍频处出现波峰并呈逐渐衰减趋势,说明该多翼离心风机气动噪声受叶片周期性旋转压力脉动影响较大。对比噪声测试结果,相对误差为2%以内,分析计算与试验相符。     

核磁共振测井仪探头设计中的数值方法

核磁共振测井仪探头的优化设计能够增强仪器的探测特性,提高仪器的信噪比, 而探头设计中的数值方法对设计结果至关重要.本文利用电磁场有限元方法对贴井壁型核磁共振测井仪探头静磁场和射频场进行了2D和3D的数值模拟,深入分析了数值模型形状、模型尺寸、 单元形状对数值模拟结果的影响,并将有限元数值模拟结果与实测数据做了对比. 结果显示:数值模拟结果与实测数据符合.在设计核磁共振测井仪探头结构时, 选取与井眼形状一致的圆形模型,模型尺寸范围在10—15倍探头外径, 并采用三角形单元可以有效提高数值模拟方法精度,增强优化设计结果的可靠性.     

Multiobjective muffler shape optimization with hybrid acoustics modeling

This paper considers the combined use of a hybrid numerical method for the modeling of acoustic mufflers and a genetic algorithm for multiobjective optimization. The hybrid numerical method provides accurate modeling of sound propagation in uniform waveguides with non-uniform obstructions. It is based on coupling a wave based modal solution in the uniform sections of the waveguide to a finite element solution in the non-uniform component. Finite element method provides flexible modeling of complicated geometries, varying material parameters, and boundary conditions, while the wave based solution leads to accurate treatment of non-reflecting boundaries and straightforward computation of the transmission loss (TL) of the muffler. The goal of optimization is to maximize TL at multiple frequency ranges simultaneously by adjusting chosen shape parameters of the muffler. This task is formulated as a multiobjective optimization problem with the objectives depending on the solution of the simulation model. NSGA-II genetic algorithm is used for solving the multiobjective optimization problem. Genetic algorithms can be easily combined with different simulation methods, and they are not sensitive to the smoothness properties of the objective functions. Numerical experiments demonstrate the accuracy and feasibility of the model-based optimization method in muffler design.     

Finite element acoustic simulation based shape optimization of a muffler

This paper describes a methodology which combines finite element analysis and Zoutendijk’s feasible directions method for mufflers shape design. The main goal is to obtain the dimensions of the acoustic muffler with the transmission loss (TL), being maximized in the frequency range of interest. The improved four parameters method is used for TL evaluations and the Helmholtz’s equation is solved numerically with the finite element method (FEM). The quadratic triangular finite element meshes are adequately constructed to control the pollution error and the optimization problem is solved using the Zoutendijk’s feasible directions method due to robustness and efficiency for problems with nonlinear constraints. Numeric experiments performed with circular expansion chambers with extended inlet and outlet show results for constrained and unconstrained shape optimization.     

Fast topological imaging

Mathematical optimization methods based on the topological sensitivity analysis have been used to develop innovative ultrasonic imaging methods. With a single illumination of the medium, they have proved experimentally to yield a lateral resolution comparable to classical multiple-illumination techniques. As these methods are based on the numerical simulations of two wave fields, they require extensive computation. A time-domain finite-difference scheme is usually used for that purpose. This paper presents the development of an experimental imaging method based on the topological sensitivity. The numerical cost is reduced by replacing the numerical simulations by simple mathematical operations between the radiation patterns of the array’s transducers and the frequency-domain signals to be emitted. These radiation patterns are preliminary computed once and for all. They were obtained with a finite element model for the anisotropic elastodynamic case and with semi-analytical integrations for the acoustic case. Experimental results are presented for a composite material sample and for a prefractal network immersed in water. A lateral resolution below 2.5 times the wavelength is obtained with a single plane wave illumination. The method is also applied with multiple illuminations, so that objects hidden in complex media can be investigated.     

Sensitivity Analysis and Computations of the Time Relaxation Model

This paper presents a numerical study of the sensitivity of a fluid model known as time relaxation model with respect to variations of the time relaxation coefficient $\chi$. The sensitivity analysis of this model is utilized by the sensitivity equation method and uses the finite element method along with Crank Nicolson method in the fully discretization of the partial differential equations. We present a test case in support of the sensitivity convergence and also provide a numerical comparison between two different strategies of computing the sensitivity, sensitivity equation method and forward finite differences.     

CFD/CSD耦合飞行器多学科优化设计新方法

文章以飞行器巡航外形为设计对象, 构建了一种新的飞行器的气动和结构特性评估方法, 即结构模型反迭代方法.该方法较传统的松耦合静气动弹性方法效率提高了4倍以上.以此为基础建立了一种新的飞行器气动/结构耦合多学科优化设计框架, 将优化效率提高4倍以上.采用数值求解N-S方程和结构有限元方程方法作为气动和结构学科的分析工具, 保证了设计结果的可信性.算例表明以巡航外形作为设计对象能够获得与传统方法一致的飞行器气动与结构特性, 以此为基础开展的无人机气动外形优化设计也获得了良好的设计结果.      

DESIGN SENSITIVITY ANALYSIS AND OPTIMIZATION OF AN ENGINE MOUNT SYSTEM USING AN FRF-BASED SUBSTRUCTURING METHOD

The FRF-based substructuring method is one of the most powerful methods in analyzing the responses of complex built-up structures with high modal density. In this paper, a general procedure for the design sensitivity analysis of a vibro-acoustic system has been presented using the FRF-based substructuring formulation. For an acoustic response function, the proposed method gives a parametric design sensitivity expression in terms of the partial derivatives of the connection element properties and the transfer functions of the substructures. The derived noise sensitivity formula is combined with a non-linear programming module to obtain the optimal design for the engine mount system of a passenger car. The objective function is defined as the area of the interior noise graph integrated over a concerned r.p.m. range. The interior noise variations with respect to the dynamic characteristics of the engine mounts and bushings have been calculated using the proposed sensitivity formulation and transferred to a non-linear optimization software. To obtain the FRFs, a finite element analysis was used for the engine mount structures and experimental techniques were used for the trimmed body including the cabin cavity. The optimization based on the sensitivity analysis gives the ideal stiffness of the engine mount and bushings. The resultant interior noise in the passenger car shows that the proposed method is efficient and accurate.     

Energy flow analysis and design sensitivity of structural problems at high frequencies

The design sensitivity formulation of an energy finite element method is presented using the direct differentiation and adjoint variable methods. The continuum method is used to derive the design sensitivity equation of the energy flow equation, whereas the discrete method is used to calculate the variation of the coupling relation. For design variables, material property, panel thickness, and structural shape are taken into account, in addition to the structural damping factor. The design variable's effect on the power transfer coefficient is discussed in detail. Even if the system matrix equation is not symmetric, the adjoint problem is solved using the same factorized matrix from response analysis. Design sensitivity results calculated from the proposed method are compared to the finite difference sensitivity results with a good agreement.     

二维多流体域耦合声场的光滑有限元解法

针对声学有限元分析中四节点等参单元计算精度低,对网格质量敏感的问题,将光滑有限元法引入到多流体域耦合声场的数值分析中,提出了二维多流体域耦合声场的光滑有限元解法。该方法在Helmholtz控制方程与多流体域耦合界面的声压/质点法向速度连续条件的基础上,得到二维多流体耦合声场的离散控制方程,并采用光滑有限元的分区光滑技术将声学梯度矩阵形函数导数的域内积分转换形函数的域边界积分,避免了雅克比矩阵的计算。以管道二维多流体域耦合内声场为数值分析算例,研究结果表明,与标准有限元相比,对单元尺寸较大或扭曲严重的四边形网格模型,光滑有限元的计算精度更高。因此光滑有限元能很好地应用于大尺寸单元或扭曲严重的网格模型下二维多流体域耦合声场的预测,具有良好的工程应用前景。      

Level set based structural topology optimization for minimizing frequency response

For the purpose of structure vibration reduction, a structural topology optimization for minimizing frequency response is proposed based on the level set method. The objective of the present study is to minimize the frequency response at the specified points or surfaces on the structure with an excitation frequency or a frequency range, subject to the given amount of the material over the admissible design domain. The sensitivity analysis with respect to the structural boundaries is carried out, while the Extended finite element method (X-FEM) is employed for solving the state equation and the adjoint equation. The optimal structure with smooth boundaries is obtained by the level set evolution with advection velocity, derived from the sensitivity analysis and the optimization algorithm. A number of numerical examples, in the frameworks of two-dimension (2D) and three-dimension (3D), are presented to demonstrate the feasibility and effectiveness of the proposed approach.