发动机辐射噪声分析

（研究生课程论文）振动与噪声控制论文题目：基于LMS Virtual. Lab边界元法发动机辐射声分析指导老师: 学院班级：学生姓名： 学 号:2015年5月基于LMS Virtual.Lab边界元法发动机辐射噪声分析摘要：在国家经济保持快速增长的背景下，国内汽车工业发展迅速。随着汽车保有量增加, 汽车噪声污染问题越来越受到人们的重视。发动机的运行噪声是车辆产生环境噪声的主要因 素，对其辐射噪声的数值分析能够为控制噪声提供良好的理论参考。本文主要介绍了外声场 分析的边界元法的基本理论，利用LMS Vutual.Lab声学模块计算了发动机辐射外声场及其 频率响应，为之后的研究学习提供参考依据。关键词：边界元法，辐射噪声，声固耦合1引言在现代汽车设计过程中，CAE分析起到越来越重要的作用，在汽车设计初期即可快速 的取得结果，从而取代后期人量的试验，使得汽车设计周期人大缩短，降低研发成本。而作 为汽车性能重要指标的NVH (Noise Vibration and Harshness)在现代汽车市场中越来越受到 人们的重视，也成为许多厂家核心竞争力的一部分，涉及车辆的振动噪声问题已经成为汽车 技术领域的一个研究热点。随着国内整机厂汽车CAE技术的成熟，利用CAE技术模拟汽车NVH问题已经不仅 仅局限于零部件及子系统的模态，基于整车模型的整车振动和噪声响应的模拟预测技术也已 经逐渐被掌握。在设计的虚拟样机阶段即可预测振动噪声水平，以便及时的更改设计，达到 可接受的振动噪声水平。发动机是汽车主要的振动和噪声源。发动机怠速时产生的振动与噪 声水平是汽车用户对汽车NVH性能的第一感觉。本文用直接边界元法计算了发动机的辐射 噪声。2数值方法的基础理论2.1边界元法的基本理论有限单元法的基本思想是将连续的求解区域离散为一组有限个、按一定方式相互联结在 一起的单元的组合体。出于单元能按不同的联结方式进行组合，且单元本身又可以有不同形 状，因此可以模型化几何形状复杂的求解域。有限单元作为数值计算方法的另一个重要特点 是利用在每一个单元内假设的插值函数来分片地表示全求解域上待求的未知场函数。由于插 值函数是已知的一个简单函数，那么有限元分析的基本未知量就是未知场函数的节点值。一 经求解出这些未知屋，就可以通过插值函数计算出各个单元内场函数的近似值，从而得到整 个求解域上的近似解。显然随着单元数目的增加，也即单元尺寸的缩小，或者随着单元自由 度的增加及插值函数精度的提高，解的近似程度将不断改进。如果单元是满足收敛要求的, 近似解最后将收敛于精确解。尽管有限元法所取得的成就与口俱増，但有限元法还不是十全十美的。改进有限元法的 努力一直在进行着，但是有限元法的某些不足是无法克服的。例如有限元法需全域离散，导 致问题的自由度和原始信息量人；对无限域只能人为地取成有限域；有限元法的离散技术本 身也存在缺陷，它把本来是连续的介质用仅在节点处连接的有限单元的集合来模拟，这样不 仅带进了离散的误差，而且在单元之间连续的要求较高时，有限单元的构造也很困难；对有 限元法的精度和町靠性也常常会提出疑问，因为对同问题采用不同的程序计算时可能会得出 不同的结果。有限元法的不足用边界元法可以弥补。边界元法仅在边界上离散，使数值计算的维数降 低一维，从而减少了问题的自由度和原始信息量。边界元法采用无限域的基本解，用边界元法求解无限域问题可称是天衣无缝。边界元法的离散误差只产生在边界上，边界元法中部分 采用数值法，部分采用解析式，它的准确度和可靠性已公认是高于有限元法。边界元法是继 有限元法之后的一种别具特色的新的数值方法，它是将描述弹性力学问题的偏微分方程边值 问题化为边界积分方程并吸收冇限元法的离散化技术而发展起来的。将弹性力学问题归结为 求解一组边界积分方程，若在边界上已知三个位移分量和三个面力分量中的三个分量，则由 边界积分方程可以确定其余三个未知分量，而任意内点的位移和应力可由六个边界分量通过 边界积分来确定，这就是边界积分方程方法。边界枳分方程有奇异性，解析求解极为困难。 有限元法所取得的成就吸引人们对边界枳分方程在边界上划分单元进行离散，然后由全部边 界节点的三个已知边界量求出全部边界节点的另外三个边界量，这就是边界元法的由来。边 界元法中包含有有限元法的思想，它把有限元法的按求解域划分单元离散的概念移植到边界 积分方程方法中，但边界元法不是有限元法的改进或发展，边界元法与有限元法存在着本质 的差异。边界枳分方程是边界元法的基本出发点，根据建立边界积分方程时，对基本解的利用方 式的不同，边界单元法可分为直接法和间接法两人类。直接法利用数学上各种积分等式，通 过基本解直接把边界上的待解边界函数与已知边界条件联系起来建立积分方程，从这个方程 解出来的就是未知边界值。间接法则不用边界的待解边界值作为未知函数，而是无限人区域 内沿着该问题的计算边界配置某种点源分布函数作为间接的待解未知量，它对计算区域的影 响，是一系列点源的影响（基本解）的叠加。由于基本解是自动满足控制微分方程的，因而 只要这些间接点源函数在边界各处产生的影响，刚好与给定的边界条件一致时，则根据定解 问题解的唯一性原理，在计算区域范I制内和边界上的影响也就是该边值问题的特解。2.2边界条件车辆噪声考虑小振幅的线性声学问题，外声场满足Helmholtz方程 V2p + k3p = 0式中，沪为拉普拉斯算子；p为声场内任一点的声压；k为波数，k = O/m, 0为角频率, m为声速。声场的边界条件包括（1）声压边界条件P= P（2）法向速度边界条件j dpPqO)dn（3）法向阻抗边界条件P= Zvn= rPqCo onPo为空气密度，j2 =式中，尸为声压函数，“为法向速度函数，n表示边界表面的法向, -b 2为法向阻抗函数。方程还应满足Somerfield远场辐射条件，以解决无限人的外声场问题。|r| （搭 + 旳）二式中，r为结构表面到空间中任意点的矢量，G为结构表面上任意点的矢量, 几为空间中任意点的矢量。Somerfield远场辐射条件保证Helmholtz方程的解在远场收敛，且没有反射波。结构的外声场问题即求解以上方程。边界元法基于Helmholtz积分方程，将区域内的偏微分方程转换成边界上的积分方程, 分割边界为有限人小的边界元素，将边界枳分方程离散成代数方程。边界元方法自动满足 Somerfield辐射条件，可以有效地求解外声场问题。边界元方法利用第二格林公式，由 Helmholtz方程和Green函数，得到边界积分方程。3模型的建立与分析求解（1）运用软件建立了发动机的三维模型，并对其进行网格划分图1发动机网格模型（2）将模型转换成由面网格构成的边界元网格图2发动机直接边界元网格模型（3）导入发动机的振动加速度由于振动加速度位于结构网格上，不再声学网格上，不能计算发动机振动引起的外声场,故需要将数据进行映射转移计算。图3为数据转移后的振动加速度云图。图3转移后的振动加速度（4）设置场点网格场点网格是以发动机的几何中心为原点的八分之一球。图4场点网格（5）定义边界条件并开始计算。利用转移后的振动加速度定义一个振动边界条件，并 计算发动机表面声压、场点声压以及场点声压频率响应函数。Pressure Amplitude dB. 1Occurre nee 986.483.881.178.575.862.559.9On Boundary图5发动机表面声圧dB云图（100Hz）Pressure Amplitude dB(RMS.lOccurrence 2974j671.1图6 300Hz场点上的声压dB云图（300Hz）参考文献1李丽君驾驶室外声场分析的边界元法与无限元法几机械强度.2012, 34 (3): 351-3542邢玉涛,李燕，洪燕等基于LMS VirtualLab的整车振动及声学响应分析J. 20123林肖辉高速汽车外声场气动噪声分析与控制D.长沙：湖南人学.2011.4刘月辉车用发动机表面辐射噪声的研究J.汽车工程200： 24 (3):213-216.