消声器设计与声学分析

消声器设计与声学分析第1页，共73页，2023年，2月20日，星期一主要内容消声器设计与声学分析基础GT-POWER中的标准声学分析GT-POWER中的线性声学分析Muffler－专用消声器设计模块GEM3D在消声器设计与声学分析中的应用第2页，共73页，2023年，2月20日，星期一消声器设计与声学分析基础第3页，共73页，2023年，2月20日，星期一消声器的设计要求声学要求减少声能的传递－传递损失插入损失空气流动要求空气流过消声元件时受到阻力流体压力上升机械和材料方面的要求气流、温度对材料性能的影响成本的控制安装空间的限制能量损失增大流速过高，摩擦噪声增大第4页，共73页，2023年，2月20日，星期一消声器的设计步骤对噪声源进行实地了解，主要实测发动机排气噪声的频谱，并且分析各个频带上的噪声声压级，通常用A声级测出噪声的声频谱，采用1/3倍频程带宽或倍频程带宽的频率分析装置作出频率特性曲线图，这是消声器设计的主要依据掌握发动机的基本参数、工作条件，如气流的流量、流速、管道的截面参数等明确要设计的消声器的性能要求，如要求达到的消声量和允许的最大压力损失选择消声器的类型，根据对消声器性能的要求、声源特性、环境的影响空间的限制等条件，选择适当类型的消声器设计计算消声器的各种结构尺寸时，选择合适的消声材料及消声结构，要根据消声原理及性能要求作理论计算，并考虑空间的限制和制造成本验算或作必要的实验，确认消声器能达到各项性能指标。根据加工及安装的工艺性要求来设计消声器，有时试验后还要进行改进才最后完成对选择的消声器，进行实际工况下的测试，对消声器性能的理论计算结果进行测试，分析误差并进行改进第5页，共73页，2023年，2月20日，星期一消声元件的结构形式的选取消声元件的分类主动：电子控制、产生与声源幅值相等、相位相反的次声波。半主动：被动控制系统，消声效果由气流控制，如排气系统中的双模态消声器。被动消声器声能被反射和吸收成本太高，汽车进排气系统中应用较少阻性抗性内部装有消声材料扩张消声器旁支消声器赫姆霍兹消声器¼波长管½波长管只能消除窄频带的噪声第6页，共73页，2023年，2月20日，星期一消元件特性扩张消音器特性：扩张比和扩张腔室的长度扩张比增大：传递损失增加（幅值）长度增加时，最大值的中心频率减少，带宽也减少赫姆霍兹特性：容积、连接管长度、连接管的截面积容积增加，共振频率降低，幅值变化没有规律截面积增加时，共振频率也增加，幅值变化也没有规律管长增加时，共振频率下降，幅值变化没有规律主管截面增加时，传递损失减少，频带变窄，但频率不变¼波长管特性：长度、面积比（旁支截面积比主管截面面积比值）长度越大，频率越低截面积比越大，传递损失幅值增大，带宽也增大，但对频率不影响用来消除较高频率第7页，共73页，2023年，2月20日，星期一排气系统噪音特性排气系统噪音组成：空气噪音（脉动噪音）；发动机的气缸数；几何尺寸；（GT主要关注点）冲击噪音；气流对管路的冲击引起；辐射噪音；机械振动；稳定气流引起；不稳定气流引起；气流摩擦噪音；紊流引起；高频；温度影响：温度高，低频降噪能力降低V型发动机，选择适当的Y型管（决定半阶次大小）第8页，共73页，2023年，2月20日，星期一例：简单扩张室式消声器结构形式的选取总容积、扩张比各腔长的分布等消声器的基本结构发动机具体型号排气噪声频谱，消声器设计的一些基本原则和规定参考原消声器或类似工况下的消声器扩张比理想的m为：7～20尽可能取在10～14之间节数通常以2～4节为宜，除非有低噪声或特殊要求，节数不应该多于4节。总体积Y为消声系数，一般取5－6n0为标定转速；N为气缸数为行程数；Vh为排量第9页，共73页，2023年，2月20日，星期一排气消声器的结构形式套管（框型罐）内装式Helmholtz

消声器三管迷路Helmholtz

消声器同心Helmholtz

消声器旁支Helmholtz

消声器穿孔管消声器Helmholtz消声器¼波长管复合消声器阻性消声器Helmholtz消声器频率作用范围：40～200Hz；三管迷路频率作用范围100～500Hz；套管加入消声材料频率作用范围：≧500HzHelmholtz腔第10页，共73页，2023年，2月20日，星期一常用声学评价指标传递损失(TransmissionLoss)消声器入射声（Lwi）和透射声（Lwt）的声功率之差。传声器声源管路声学元件全消音传声器加入全消音装置，防止声波继续反射回来第11页，共73页，2023年，2月20日，星期一常用声学评价指标插入损失（InsertionLoss）系统中插入消声元件之前和之后，在出口处得到的声功率级（或声压级）的差值。声源声源管路管路声学元件传声器传声器W1：没有安装消声元件的系统在测量点的声功率W2：安装了消声元件后在同一点测量的声功率第12页，共73页，2023年，2月20日，星期一常用声学评价指标声压级差值（leveldifference）系统中任意两点的声压级的差值声源管路声学元件传声器1传声器2第13页，共73页，2023年，2月20日，星期一评价指标的比较

传递损失:只要有截面变化，就存在传递损失，而且永远是正值。只与消声元件的结构、介质的阻抗率以及截面面积有关。而与消声元件在系统中的位置没有关系。用于描述单个元件。插入损失：取决于传递损失和消声元件在系统中的位置有关以及声源和出口处的声学特性。可能是正值也可能是负值。用于描述整个系统。声压级差值：与声源特性无关，但却取决于出口处的声学特性。用于评价出口处的最终的噪声指标。第14页，共73页，2023年，2月20日，星期一GT-POWR中的标准声学分析

——非线性声学第15页，共73页，2023年，2月20日，星期一声学预测工具提供了多种预测整个系统声学特性的工具运用非线性预测工具可以把频域信号转化成时域信号麦克风（传声器）传递损失插入损失传递损失（2Mic）声音文件转换快速傅立叶变换阶次追踪数据提取第16页，共73页，2023年，2月20日，星期一OrderTrackingandHarmonics例:四冲程发动机@6000rpmBasicf:50HzRotationf:100Hz1st，2nd，4th

harmonics:50Hz，100Hz，200Hz1st，2nd，4th

orders:100Hz,200Hz,400Hz.对于二冲程发动机：Harmonics=OrdersEngineCycleEngineRevolutionBasicFrequencyRotationFrequencyordersHarmonics*Mult.*Mult.第17页，共73页，2023年，2月20日，星期一采样率、存储率和频率混叠GT-POWER不能预测频率高于1000Hz的噪音。在进行声学性能分析时要用小的离散长度在使用外部数据进行声学性能分析的时候，采样频率也非常重要，要保证所有的采样的点都被存储起来。根据Nyquist理论，当采样频率低于信号中最高频率的2倍的时候，将发生混叠数字滤波器不能消除混叠现象，它只与采样频率相关第18页，共73页，2023年，2月20日，星期一采样率、存储率和频率混叠由于GT-POWER采用的时间步长很短，所以基本上能消除频率混叠现象。输出数据是离散的如果每个时间步长内的数据没有完全被存储起来将发生aliasing（频率混叠）一般只有在预测高频事件中才会产生频率混叠现象当采样频率≤12691Hz时预测的结果是准确的第19页，共73页，2023年，2月20日，星期一‘AcoustExtMicrophone’相当于一个外置的随机的白噪声麦克风，可根据声源特性将气流的压力波模拟成声波通过在自由场中放置麦克风，计算由于流动脉动引起的噪声假定脉动流动是由于单极子的振动辐射产生计算通过噪声（考虑到多普勒效应）把麦克风和车辆之间的距离设置为定值，即传声器随车辆一起运动。把麦克风的空间位置确定，车辆运动考虑了车辆和外置传声器之间的相对运动，产生多普勒效应当车辆朝着传声器方向运动的时候，由于相对运动，导致传声器测得的波长变短、频率变高当车辆远离传声器方向运动时候，由于相对运动，导致传声器测得波长变长，频率变低Multiplesourcesmayinfluence1microphone(dualexhaust)在V6.2版本以后可以考虑流动噪声（见注释）第20页，共73页，2023年，2月20日，星期一‘AcoustExtMicrophone’该模块用于模拟过程中也可以用于计算结果的后处理中输出结果：dBvs.Frequency（频谱分析）dBOrderTracking（阶次分析）LinearOrderandFrequencyContours（？？）dBOrderandFrequencyContours（Cambell图）Post－processDuringasimulation第21页，共73页，2023年，2月20日，星期一Microphone－计算方法传感器传出每个时间步长的速度信号麦克风把速度信号转换成压力信号GT-POWER对压力信号进行插值用于FFT转换对压力信号进行FFT转换计算每个频率下声波的振幅2134第22页，共73页，2023年，2月20日，星期一‘AcoustExtMicrophone’输入参数推荐频域分析的最大频率≦1000Hz应用24阶低通滤波器瞬态窗函数的选择稳态工况：NO瞬态工况：Hanning窗Coherence（决定声源基频的一致性）：not－Coherence球面声波or半球面声波实际情况应该是介于两者之间预测结果相差3dB把非周期信号转换成周期信号第23页，共73页，2023年，2月20日，星期一采用Chung和Blaser的理论和方法在4个外置麦克风间通过自相关谱和互相关谱得到声功率之差由speaker’提供白噪声声源必须采用无回声的终止端（防止声波的反射）四麦克风方法用于把压力分解到正向和反向部件‘AcoustTransLoss’第24页，共73页，2023年，2月20日，星期一‘AcoustTransLoss’【注意】所有的传感器必须放置在管路的中心位置需要检查离散长度、整个管路的长度以及传感器的位置DistanceinobjectshouldbesameasdistancebetweenSensorConns？？？模型中传感器的位置必须与实验设备的位置一样在做试验对比的时候要选择合适的测点位置传递损失对采样点的多少十分敏感一般设置最大时间步长为0.088deg，对应的FFT中最大采样点为4096。

numberofpoints:assumesthatadriverof360degreesisused.(derivationisleftasexercisefortheuser)（？？？）第25页，共73页，2023年，2月20日，星期一‘AcoustTransLoss’传递损失分析模块用法：

Duringthesimulation（正在进行的模拟）

Post－Processing（后处理）‘data’模块用于从模拟结果中得到数据确保模拟之初所有数据点都能储存到‘data’中第26页，共73页，2023年，2月20日，星期一‘AcoustTf2Mic’在管路内部的两个位置放置两个麦克风用于预测管路内部的动态声压传递函数定义为每个频率下信号2和信号1的振幅之比考虑到声波反射和尾管的影响结果可以用dB表示的降噪量？？signal1signal2第27页，共73页，2023年，2月20日，星期一‘AcoustInsLoss’原理是计算两个外置麦克风之间的声压级之差。由于插入损失与声源特性有关，因此计算时必须要将消声器与发动机相连。提供插入消声元件后的声压级损失分析步骤：不插入消声元件，储存流速信号在‘data’模块中插入消声元件，用存储在‘data’中的信号以及‘AcoustInsLoss’计算插入损失第28页，共73页，2023年，2月20日，星期一Filter（滤波器）采用低通滤波器消除采样数据中的高频信号消除在差值过程中Contamination的可能性？24th

阶次的滤波器应用效果最好（TransLoss、Mic）第29页，共73页，2023年，2月20日，星期一‘AcoustToWAVFile’用于产生‘*.wav’文件稳态计算：重复最后一个循环的模拟结果声音文件的时间由用户控制瞬态计算：声音文件的时间等于整个瞬态过程的时间输入参数推荐：Bits:16SamplingRate:44.1kHzFilter:54thorder第30页，共73页，2023年，2月20日，星期一DataWindows（瞬态窗）当采样的数据为非周期性信号，采用FFT时会产生旁瓣现象。通过对数据加窗函数可以消除旁瓣现象窗函数可以改善信号的频率响应，但加窗的最主要的目的是消除旁瓣（推荐使用Hanning窗）窗函数常用于瞬态分析（转速变化比较剧烈的情况）第31页，共73页，2023年，2月20日，星期一GT－power中的声学评价指标传声损失分析插入损失分析实际试验中测量传递损失是非常困难的。它反应的是消声器本体的声学性能，是消声器的固有特性。通常用与不同结构的消声器进行对比分析。此外特别适合缺乏发动机参数的情况，也可以将由试验取得的发动机出口处的状态参数输入到声源模块且计算速度快。安装消声器前后在发动机排气口某固定测点处测得的计数声级之差，它直接反映了消声器对发动机噪声的消声效果。该指标能综合反应内燃机、消声器和消声器尾管的匹配特性。由于插入损失考虑声源特性，因此必须建立发动机模型。通常情况下计算插入损失必须建立两个模型一个带消声器的模型和一个不带消声器的空管模型，两个模型在排气口处的声压级之差即为插入损失。此声波可以转化成声音文件播放出来，形象直观。第32页，共73页，2023年，2月20日，星期一GT-POWR中的线性声学分析

第33页，共73页，2023年，2月20日，星期一线性声学模型计算方法：传递矩阵（频域分析）输入参数：几何尺寸、边界条件、初始条件适用对象：有声学分析经验的工程师（使用过其他相似的声学分析工具）模型中常用的流动部件：PipeFlowsplitOrificeCatalystBrickEndEnvironmentEndFlowCap第34页，共73页，2023年，2月20日，星期一线性声学模型采用传递矩阵（4－pole）来描述系统的属性对于每个组件都有一个独立的传递矩阵例：从位置1到位置2的传递矩阵为每个单独的传递矩阵合并成一个大的传递矩阵描述整个系统第35页，共73页，2023年，2月20日，星期一线性声学模型计算方法：基于平面声波理论－HelmholtzequationInputRespondsTransferMatrixBlackBox压力、流量的传递矩阵形式正向、反向的传递矩阵形式声波传播、反射的传递矩阵形式第36页，共73页，2023年，2月20日，星期一线性声学模型的应用计算要求：计算速度快、只反应系统的整体性能不在时域上进行计算不能考虑详细细节（部件的流速、压力etc.)不考虑空间位置只反映整个系统的响应适用的频率范围：Airbox——声速为343m/s等效直径为200mmMuffler——声速为500～600m/s等效直径为350mm最大频率≤1000Hz主要应用范围：比较几个消声器的性能差异第37页，共73页，2023年，2月20日，星期一线性声学分析的影响因素平均流动考虑了Mach数的影响质量流量可以通过每个连接件的“InitialMassFlowRate”

获得ViscoThermal考虑流体粘性和热力学性质的影响温度可以通过每个连接件的“WallTemperature”

获得第38页，共73页，2023年，2月20日，星期一线性声学预测工具

Template Description AnalogAcoustSource 线性声源 EndFlowSpeakerEndFlowRadiation 线性声学 EndEnvironment, 边界条件 EndFlowAnechoicAcoustLinEigen 本征频率的计算 None(order/mode (naturalfrequencies) storageflag)AcoustLinTransLoss 传递损失 AcoustTransLoss &transfermatrixAcoustLinExtMic 声压级的计算 AcoustExtMicrophone &moreAcoustLinInsLoss 插入损失 AcoustInsLossTMBlackBox 传递矩阵 Data第39页，共73页，2023年，2月20日，星期一线性声源把声源特性定义以下为频率以及每个阶次下发动机转速的函数一般用‘speaker’模块声源特性不是十分重要（传递损失的计算）一般性的分析特殊的声源模型数据从‘multiload’中获得每个发动机对应一个声源第40页，共73页，2023年，2月20日，星期一RadiationBoundary出口边界条件出口边界可以定义为：Openend（球型)Flange（半球型）Anechoictermination（全消音室）Data（输入边界的声阻抗）可以进行边界声阻抗的建模声阻抗数据可以通过外部文件读入第41页，共73页，2023年，2月20日，星期一Natural(Eigen)Frequencies(固有频率的计算）固有频率可以反应整个系统的响应状况反映的是系统受到激励时自由振动的频率在一个固定激励作用下进行频率扫略可以确定系统的共振频率所有的固有频率都可以在后处理里中显示出来？？？第42页，共73页，2023年，2月20日，星期一例：Helmholtz腔和1/4波长管的固有频率分析固有频率分析系统的固有频率RLT结果显示第43页，共73页，2023年，2月20日，星期一例：1/4波长管传递损失的计算第44页，共73页，2023年，2月20日，星期一例：排气消声器声压级（差）的计算计算插入损失第45页，共73页，2023年，2月20日，星期一例：扩张室插入损失的计算第46页，共73页，2023年，2月20日，星期一‘TMBlackBox’（黑匣子传递矩阵）允许输入端为传递矩阵可以替代一个组件或者系统‘TMBlackBox’获取途径：应用‘TMatrixGenerator’进行标准的GT-POWER计算其他相似软件试验数据理论公式用户子程序第47页，共73页，2023年，2月20日，星期一BridgingtheGapBetweenLinearandNon-LinearAcousticAnalysisNo－LinearAcousticLinearAcoustic‘MultiLoad’'AcoustSource'产生声源特性（p/z）‘TMatrixGenerator’TMBlackBox’产生传递矩阵（TM）

Gap第48页，共73页，2023年，2月20日，星期一TransferMatrixGenerator用带有边界条件的标准（非线性）的声学计算产生用于线性声学的Blackbox第49页，共73页，2023年，2月20日，星期一MultiloadMethodofSourceCharacterization基于时间的流动/压力函数转换成基于频率的压力/阻抗函数间接计算声源特性计算方法：内部计算模型外部计算模型外部测量模型第50页，共73页，2023年，2月20日，星期一管路长度必须可变内部计算模型：需要系统的体积流量信号和压力信号外部计算模型：在系统端点传感出体积流量信号需要加入线性声学分析模块外部计算模型：通过microphone传感出系统的压力信号需要加入线性声学分析模块管路长度必须为可变MultiloadMethodofSourceCharacterization第51页，共73页，2023年，2月20日，星期一线性声学模型利用黑盒子产生传递矩阵系数用于线性声学分析插入损失的计算传递损失的计算通过标准求解器计算传递矩阵本征频率的计算第52页，共73页，2023年，2月20日，星期一Muffler－专用消声器设计模块第53页，共73页，2023年，2月20日，星期一Muffler简介特点：专用的图形化前处理程序图形化的界面建立包括各种零部件的复杂结构消声器的三维几何形状通过其内部数据格式数据文件自动生成GT-POWER模型减少消声器建模的工作量(不必做大量的加工和测试）模型优化而准确为复杂形状的建模提供方便第54页，共73页，2023年，2月20日，星期一Content1.基本概念2.基本操作3.离散4.Interference5.模型的快速检查第55页，共73页，2023年，2月20日，星期一1.基本概念1.1模板第56页，共73页，2023年，2月20日，星期一1.Basicconcept1.2Shell（壳体）定义壳体的厚度用于GT-POWER中的参数用于GT-POWER中的参考模型第57页，共73页，2023年，2月20日，星期一1.Basicconcept1.3Shellbaffle（挡板）Leftendatx=-1Rightendatx=1MiddlefacepositinNote:Bafflethicknessreducestheeffectivechambervolume!第58页，共73页，2023年，2月20日，星期一1.Basicconcept1.4OrificeDefinetheinletandoutletportnumberforassemblylinkageinGT-POWER第59页，共73页，2023年，2月20日，星期一1.Basicconcept1.5Porocity第60页，共73页，2023年，2月20日，星期一1.Basicconcept1.6WOOL（吸声材料）Note: