基于ANSYS的抗性消声器性能仿真分析.pdf
《基于ANSYS的抗性消声器性能仿真分析.pdf》由会员分享,可在线阅读,更多相关《基于ANSYS的抗性消声器性能仿真分析.pdf(6页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、收稿日期:!#$%$&(第!)卷(第*期计(算(机(仿(真!+年*月(文章编号:&+$,)*(!+)*$)+$#基于!#$#的抗性消声器性能仿真分析张乃龙,杨文通,费仁元(北京工业大学机电学院,北京,&!)摘要:运用有限元分析软件-./0/对内插管抗性消声器内部声场及流体进行了有限元分析。在讨论了管道声场流体运动方程的基础上,首先建立抗性消声器的流体有限元分析模型和声场有限元分析模型,以定性分析为目标,通过加载、求解以及后处理等一系列的步骤对消声器内部声场及流体进行计算分析,获得消声器内部不同频率下的声压及流体速度和压力分布情况。通过分析结果可以方便直观地获取消声器内插管中心偏移对消声器插入损
2、失及流体流速和压力的影响,分析结果为消声器的优化设计提供可靠的依据。关键词:抗性消声器;声学仿真;流体动力学仿真;有限元中图分类号:12),&3%(文献标识码:-#%&()*%+,),-!,)(./%/+0 12)3*%42 500(26 7)/2-+,!#$#+0*8)6245-.6.78$9:;$?:;$DE7;(F8G8;JK8?D:L 1MN;:9:8G8;9P;?7;79DK8K K:L?Q7J-./0/8K 7RR98S?:7;79DT K:E;S L89S 7;S L9E8S:L PELL9J3-L?J?79U8;L9:Q M:;?8;E8?D WE7?8:;:L R8R K:E
3、;S L89S,L8JK?9D?N J7M?8I PELL9J L9E8S L8;8?P:S97;S K:E;S L89S L8;8?P:S9 7J VE89?,?N;7M:EK?8M RJKKEJ 7;S L9E8S I9:M8?D 7;S L9E8S RJKKEJ S8K?J8VE?8:;:LJ7M?8I PELL9J E;SJ S8LLJ;?879 LJWE;MD 7J:V?78;S LJ:P 7 KJ8K:L K:E;S L89S 7;S L98ES 7;79DK8KRJ:MSEJ,KEMN 7K 9:7S8;,K:9I8;KK8;M?:;BX 7;S RJKKEJ M:E9S V 8
4、;?E8?8I9D:V?78;S3B?RJ:I8SK 7 L:E;S7?8:;L:J:R?8P79 SK8 PELL9J3;1?#:C7M?8I PELL9J;-M:EK?8M K8PE97?8:;;5DSJ:SD;7P8M K8PE97?8:;;8;8?9P;?A引言抗性消声器主要利用声抗的大小来消声,通过各种不同形状的管道和腔室进行适当地组合,提供管道系统的阻抗失配,使声波产生反射或干涉现象,从而降低由消声器向外辐射的声能。抗性消声器的性能和管道结构形状有关,一般选择性较强,适用于窄带噪声和低、中频噪声的控制。由于其在低中频较好的降噪效果及其较小的压力损失在发动机及内燃机的排气消声处理中广
5、泛使用。但是随着消声量的增加,其引起的压力损失相应增加,造成设备功率损耗增加。有效调整消声与压力损失的关系是在设计消声器时必须考虑的问题。本文利用有限元分析软件-./0/的声学分析模块对某型号的扩张式抗性消声器进行声学及流体动力学分析,以对其进行优化设计提供依据。BA声场流体理论基础在声场的流固耦合问题中,要把结构的动力方程与流体斯托克斯方程中的动量方程和连续性方程综合考虑。声学基本方程是在空气作为特殊流体条件下流体方程的简化。在不改变问题性质的前提下,作如下假设:&)流体是可压缩的,即密度随压力变化而变化。!)媒质为理想流体,即媒质不存在粘滞性,声波在其传播时没有能量的损耗,流体没有扰动和紊
6、流。))声传播是一个绝热过程。)媒质的静态压强 2和静态密度!都是常数。BC A数学模型的建立图&所示为轴对称抗性消声器。在!域的二维平面波动方程E!!$%!(&)其中#介质中声传播速度!#流体平均密度+)图!轴对称抗性消声器示意图!流体体积模量 声压(#($,%,&,)时间E!拉普拉斯算符对于谐波设解为#()*(!)式中#压力幅)#+!#!,简谐振动的圆频率,#振动频率式(!)代入()可得:!-!.E!#($)令!#!-(称为波数),则式($)可表示为:!.E!#(%)式()以矩阵形式表示为-!#!#!+/0(/)#(&)式中:/0#$#%#&通过用迦辽金法对方程(&)离散化即得到单
7、元矩阵,在方程(&)左右同时乘以一个声压变化值,然后在一定区域内对体积积分,)123-!$#!#!4(123).)123(/0$)(/)4(123)#)560$(/)4(5)()其中:123 为一定区域的体积$为一定的声压变化值5 为声压向量所指向的表面6界面 5 的单位法向量在声固耦合界面问题中,面被视为界面,由于简化假设,流体的动量方程中法向声压梯度与结构的法向加速度在界面处遵循以下规律:6 E#+%#6#!7#!(()其中7为结构在界面处的位移向量,用矩阵形式表示,即为:60/#+%#60#!#!7()())把())式代入式(),积分变为:)123-!$#!#!4(123).)123(/
8、0$)(/)4(123)#)5%#$60#!#!7()4(5)(*)方程(*)中包含的变量有:声压,结构位移8,1,*。定义有限单元的近似形函数为:#90((#)7#9:07(()其中9为声压形函数9:为位移形函数(为节点声压向量7(#8(,1(,*(为节点位移向量由方程(#)和(),声压和位移的时间二阶导数及声压变化值可表示为:#!#!#90;(!)#!#!7#;7(($)$#90$((%)将式(#)代入方程(*),声波方程的有限元表达式则为:)123-!$(09 904(123);(.)123($()004(123)(.)=%#$(09 609:04(5);7#(&)其中:#/906为流体
9、边界向量由$(.#,等式两边同时消去$(,并把非变量提出积分,得到:-!)1239 904(123);(.)1230?(;(.?(.%#A(0;7#(()其中:?(#-!)1239 904(123),为声场流体的质量矩阵;?(#)1230 4(123),为流体的刚度矩阵;%#A(0#%#)=9 609:04(5),为流固耦合界面的质量矩阵。#消声器的$%&分析用 B95C5 对两种抗性消声器进行声学及流体动力学分(#$析,观察内部结构变化对消声性能及其压力损失的影响。消声器结构如图!所示。图!抗性消声器结构模型#$%消声器的流体动力学分析运用#$%$软件对消声器内部气体流场进行有限元仿真分析,
10、分析内部压力分布及流速变化情况,寻找对压损失影响较大的部位,为结构优化提供指导。#$%$的流场可以归纳为以下七个主要步骤:&)定义单元类型 对于本文的研究采用二维()*+&,&单元。!)建立流场的范围建立消声器内部的流场模型来约束气体的流动范围,并把出口端延长&-倍出口直径的长度,以模拟无穷远端的压力值。.)生成有限元网格采用映射网格划分模式对流场模型进行网格划分,这样可以得到比较均匀的网格。,)施加载荷和边界条件在消声器的入口处定义来流速度为/-01 2;在流体与消声器内壁接触的边界定义节点速度为零,以固定流场边界;在&-倍出口直径处定义出口压力值为零。/)定义分析类型和求解参数定义流体的物
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 ANSYS 抗性 消声器 性能 仿真 分析
限制150内