声呐导流罩内声场特性研究

为研究机械激励下声呐导流罩壳内噪声预报问题,采用数值计算(有限元法-边界元法)结合实验的方法对一导流罩模型的壳内声场进行研究。首先依据实际导流罩模型建模,并加载机械激励进行数值计算,得到壳内声场的均方声压,然后进行湖内实验测量,最后对比了数值计算和实验的结果。研究表明:有限元法加间接边界元法是分析导流罩内声场这种封闭结构流-固-声耦合问题的有效方法;导流罩壳内声场的数值计算结果与实验测量结果吻合较好;声场均方声压与结构均方振速之比反映了罩壳结构内平均声场特性,可用来对比分析数值计算和实验的结果。     

边界面法数值结果的云图表征

边界面法继承了传统边界元法的优点,并将几何实体的边界曲面离散为参数空间里的曲面单元,在处理一些特殊问题如移动边界、高梯度、大变形等方面显示出特殊的优越性。但是也使得计算结果的后处理遇到困难。提出了一种基于黎曼度量推进波前法生成三角背景网格的实用边界面法计算结果后处理方法。该法对求解域剖分成三角背景网格然后将计算结果映射到网格节点上,通过区域填充获得计算结果的云图。温度场的数值算例表明该方法可靠实用。     

谱元法与透射边界的配合使用及其稳定性研究

在无限域波动模拟中引入透射边界条件时,目前多将边界上的透射公式与内域的有限元法结合使用,其计算精度由有限元方法决定,而谱元法因结合有限元和频谱法的优势则比有限元空间域积分具有更高的计算精度。该文基于谱元法非等距网格划分特性,研究了内域的谱元法与边界上的透射公式结合的理论方法,给出了相应的透射公式使用方法,并基于建立的谱元法波动数值模型探讨了透射公式的稳定性问题。研究表明:空间域插值系数需控制在一个合理范围内,空间域插值方法相对于时间域插值方法更为稳定,高频失稳出现可能性相对较小;Gamma算子的使用可提高模拟的精度,采用Gamma算子后对于高阶透射公式仍可出现低频漂移现象,可结合降阶消漂的方式实现稳定精度高的透射边界应用。     

结构声辐射计算的边界无网格法

研究函数有限维逼近插值形函数的一般要求,介绍采用移动最小二乘构建无网格插值形函数的方法与步骤;通过配点法将Kirchhoff-Helmholtz边界积分方程离散为受边界条件约束的线性方程组;最后通过分块矩阵法求解约束方程组,得到离散后的声辐射传输模型数值表达式。在计算实例中,分别用边界无网格法和边界元法建立声辐射传输模型进行声场计算,计算声场值与解析值相对比的结果表明,由于边界无网格法插值形函数根据求解情况自行构建,因此更灵活,具有更高的插值和计算精度。     

半无限域中结构体辐射声场计算

本文从无限域Ｈｅｌｍｈｏｌｔｚ积分方程入手，用边界元法（ＢＥＭ）计算任意形状结构体辐射声场，对用ＢＥＭ计算声场出现的奇异积分，特征频率，边角点处向不连续，高频段计算误差较大等问题进行了有效，简便的处理，并用理论算例验证了这些方法有较精度。     

消声水池声场仿真研究

采用网格计算方法建立了消声水池声场的低频仿真分析模型,利用该模型对消声水池中三维声场进行了仿真分析,验证了其可靠性,并将分析结果与镜像声学方法计算进行了比较.通过分析得到不同频率不同吸声效果下消声水池声场的起伏程度等参数.     

求解远场声辐射问题的一种快速边界元法

针对传统边界元法奇异性的固有缺陷,研究一种求解结构辐射声场的间接边界元法。通过在内域虚拟曲面上配置一系列单偶极子源,等效模拟结构外辐射声场,并根据单元体积速度匹配法的边界条件方案,给出数值计算辐射声场的间接边界元模型。结合球冠形辐射器的仿真算例,表明该间接边界元法在求解远场声特性时具有较高的计算效率和精度,为分析设计具有大空间指向性的辐射器提供参考。     

城市轨道车辆车内噪声特性研究

采用有限元法结合边界元法研究按额定速度行驶的城市轨道客车车内噪声特性。车辆行驶工况下车身所受到的结构激励力以车底架上的载荷谱描述,采用有限元计算该激励下车体表面的振动速度,并采用边界元法预报相应的车内声场。最后,分析车身表面不同位置对声场最大声压级点的贡献量及其随频率和空间分布的规律。     

管道声传播问题的一种高精度数值模拟

采用高精度的切比雪夫谱元法对管道声学模型的声场进行数值模拟。该解法对基于线性化欧拉方程的声场传播控制方程在空间上进行谱元离散,并推导了隐式时间积分的公式。通过Gaussian扰动波问题进行了算法测试,然后对带有管壁反射边界和进出口吸收边界的管道声学模型进行数值模拟,并给出具体算例进行求解分析,其结果符合线性声学理论,表明谱元方法对求解和分析计算气动声学问题具有良好的适应性和有效性。     

柴油机油底壳辐射声场的仿真研究

采用有限元与边界元联合求解的方法进行了柴油机油底壳辐射声场的仿真研究,利用ANSYS软件进行了油底壳的谐响应分析为SYSNO ISE软件的边界元法声场模拟提供边界条件,计算结果与实验结果的高度一致表明该方法是完全可行的,并且可以为低噪声设计提供指导。     

Iterative coupling algorithms for large multidomain problems with the boundary element method

A new parallel Robin-Robin adaptive iterative coupling algorithm with dynamic relaxation parameters is proposed for the boundary element method (BEM), and relaxation parameters are derived for other existing iterative coupling algorithms. The performances of the new algorithm and of the modified existing algorithms are investigated in terms of convergence properties with respect to the number of subdomains, mesh density, interface mesh conformity, and BEM element types. Results show that the number of subdomains and the refinement level of the mesh are the two dominant factors affecting the performances of the considered algorithms. The proposed parallel Robin-Robin algorithm shows the best overall convergence behavior for the tested large problems, thanks to its effectiveness in handling complex boundary conditions and large number of subdomains, thus resulting to be very promising for efficient parallel BEM computing and large coupling problems. Source code is available at https://github.com/BinWang0213/PyBEM2D .     

A meshless singular boundary method for three‐dimensional elasticity problems

This study documents the first attempt to extend the singular boundary method, a novel meshless boundary collocation method, for the solution of 3D elasticity problems. The singular boundary method involves a coupling between the regularized BEM and the method of fundamental solutions. The main idea here is to fully inherit the dimensionality and stability advantages of the former and the meshless and integration‐free attributes of the later. This makes it particularly attractive for problems in complex geometries and three dimensions. Four benchmark 3D problems in linear elasticity are well studied to demonstrate the feasibility and accuracy of the proposed method. The advantages, disadvantages, and potential applications of the proposed method, as compared with the FEM, BEM, and method of fundamental solutions, are also examined and discussed. Copyright © 2015 John Wiley & Sons, Ltd.     

Dynamic analysis of large-scale SSI systems for layered unbounded media via a parallelized coupled finite-element/boundary-element/scaled boundary finite-element model

An algorithm for a parallelized coupled model based on finite element method (FEM), boundary element method (BEM), and scaled boundary FEM (SBFEM) for harmonic and transient dynamic response of large-scale 2D structures embedded in or on layered soil media is presented. The BEM and SBFEM are used for modelling the dynamic response of the unbounded media. The standard FEM is used for modelling the finite region and the embedded structure. The objective of the development of this parallelized coupled model is to use the power of high performance computing, and to take into account the advantages and evade the disadvantages of the above mentioned numerical methods for modelling of the unbounded media in soil-structure interaction (SSI) systems. The development of the parallel algorithm for this model is essential for solving arbitrarily shaped large-scale SSI problems, which cannot be solved within reasonable elapsed times by a serial algorithm. The efficiency of the proposed parallel algorithm and the validity of the coupled model are shown by means of three numerical examples, indicating the excellent accuracy and applicability of the parallel algorithm with considerable time-savings in large-scale problems.     

Computational algorithms and applications of element-free Galerkin methods for nonlocal damage models

In using complex material models, especially the strain-gradient-dependent damage models, the convergence of the finite element computation increasingly becomes a problem. Due to large strains in damaging elements the computation may often result in non-convergence. For the higher-order gradient plasticity the special element formulation would often be necessary, which causes additional difficulties in implementation and computations. In recent years, meshless methods have been developed as an alternative to the finite element method (FEM) and can overcome some known shortcomings of FEM. In the present paper an algorithm of element-free Galerkin (EFG) methods for strain-gradient based nonlocal damage models has been developed and used to simulate ductile material damage. The method provides a reliable and robust results for material failure with large damage zones. The strain gradient-dependent terms can be evaluated from the direct differentiation. The investigation confirms that the nonlocal damage model with element-free Galerkin method is suitable for computing the damage problems and predicting the size effects. With the help of the meshless method, material failure in specimens as well as the size effects are predicted accurately.     

声学-结构设计灵敏度分析

声学 -结构设计灵敏度分析揭示了结构振动引起的辐射声压与结构设计变量之间的关系。分别用有限元法和边界元法计算结构设计灵敏度和声学灵敏度。将两个灵敏度结合得到最终的声学 -结构设计灵敏度。在边界元计算中 ,采用退化元处理奇异积分问题 ,对特征频率不唯一问题采用CHIEF方法处理。以脉动球和箱体为例 ,验证了算法的可行性和精确性。     

Topology optimization of exterior acoustic-structure interaction systems using the coupled FEM-BEM method

This paper presents an efficient topology optimization procedure for exterior acoustic-structure interaction problems, in which the coupled systems are formulated by the boundary element method (BEM) and the finite element method (FEM). So far, the topology optimization based on the coupled FEM-BEM still faces several issues needed to be addressed, especially the efficient design sensitivity analysis for the coupled systems. In this work, we contribute to these issues in two main aspects. Firstly, the adjoint variable method (AVM) formulations are derived for sensitivity analysis of arbitrary objective function, and the feedback coupling between the structural and acoustic domains are taken into consideration in the sensitivity analysis. Secondly, in addition to the application of fast multipole method (FMM) in the acoustic BEM response analysis, the FMM is now updated to adapt to the arising different multiplications in the AVM equations. These accelerations save considerable computing time and memory. Numerical tests show that the developed approach permits its application to large-scale problems. Finally, some basic observations for the optimized designs are drawn from the numerical investigations.     

基于拓扑优化的结构加强筋布局降噪方法研究

基于变密度拓扑优化方法,提出了通过优化加强筋布局来降低谐振结构辐射声功率的策略。优化中将加强筋单元的伪密度作为设计变量,约束为加强筋质量上限,所用寻优算法为移动渐近线(MMA)系列算法中的MMA-GC-MMA混合算法。对结构的振动响应使用有限元方法求解,对声辐射使用边界元方法求解,并在此基础上分析了目标函数对设计变量的灵敏度。以加筋箱体为例进行优化,验证了所提方法在降噪设计中的可行性和有效性,并对结果进行了讨论。     

低压母线系统振动与声辐射特性分析

采用有限元与边界元法,以某型号母线桥为例,对低压大电流母线设备正常运行中噪声偏大的问题进行研究。建立母线桥有限元结构分析及边界元声学分析的三维模型。在求解系统固有频率及振型的基础上,应用有限元法完成了母线桥结构振动响应的数值计算。然后,将母线桥外表面的节点振动位移作为声学边界条件,采用边界元法计算了母线桥的辐射噪声特性,得到母线桥辐射声场的分布规律,确定系统中主要振源和声源的部位。根据分析结果提出了结构改进方法,计算和实验结果表明该方法控制噪声的可行性和有效性。     

自然元与无限元耦合方法在岩土工程粘弹性分析中的应用

自然单元法是一种新的偏微分方程数值解法,由于其位移插值函数采用无网格的方式构造且形函数满足插值性质,从而克服传统有限元方法对单元网格信息的依赖,大大简化数值计算的前处理过程,同时又能像有限元那样准确施加边界条件,在岩土工程中具有广阔的应用前景;介绍了自然元与无限元的基本原理,针对在处理岩土工程无限域或半无限域问题时需要人为确定边界条件而带来计算误差的问题,引入无限元模拟无穷远处边界条件,与自然元相结合形成耦合分析方法;并根据粘弹性理论,采用Laplace插值,编制了基于自然元与无限元耦合方法的二维粘弹性分析程序,通过算例验证了算法的正确性,结果也表明相对于纯自然单元法,耦合方法能够显著提高分析结果的精度,在此基础上拓展了自然单元法在岩土工程中的应用范围.     

基于有限元素法的小尺度封闭空间低频声传递函数求解方法研究

如何求解封闭空间中声源点到接收点的低频声传递函数已成为目前小尺度封闭空间可听化技术研究的关键问题。具有较高计算精度并能够处理任意结构形状的有限元素法可作为求解该问题的合适的方法。本文推导出了小用于小尺度封闭空间低频声传递函数的有限元素求解方法,并编制了相应的计算机程度。最后的算例中通过与经典模态迭加法给出的结果进行比较,验证了该方法的正确性。