djh713 发表于 2006-11-9 12:54

求助：消声器问题

用SYSNOISE做消声器传递损失分析，DBEM、IBEM、FEM都试过了，可结果就是和理论求解不一样，趋势差不多，最小传损处也一样，可峰值差了近两倍，而且在300Hz以下也不一致，到底是什么原因啊，我做的就是一个简单的一个扩张腔的消声器。

FEM轴对称：

OptionFEM Frequency Fluid Return
Import Mesh Format NASTRAN File femaxi.DAT Return
AxiSymmetry Return
Set 1 Name "inlet"
    Face 54
    Return
Set 2 Name "outlet"
    Face 59
    Return
Material Fluid
    Name 'air'
    Sound 3.4000e+002Rho 1.2250e+000
    Elements All
    Return
Boundary VelocityReal -1 Imag 0
    Faces Set 1
    Return
Boundary AdmittanceReal 2.401e-3 Imag 0
    Faces Set 2
    Return
SolveFrequency 50 To 3000 LinStep 50 Return
Combine
Read Node 11 Pressure Return
Amplitude
Read Node 1 Pressure Return
Read Node 1 VelocityVector 0 1 0 Return
Constant Real 416.5 Imag 0
Multiply
Add
Constant Real 2 Imag 0
Divide
Amplitude
Divide
Log10
Constant Real 20 Imag 0
Multiply
Write Table 1 name tl Return
Return


BEM:

Import Mesh Format NASTRAN File BEM.DAT Return
Check Mesh Merge Return
Material Fluid
    Name 'air'
    Sound Real 3.4000e+02Rho Real 1.2250e+00
    Return
Boundary VelocityReal 1 Imag 0
    Elements X=0.06
    Return
Boundary ImpedanceReal 416.5 Imag 0
    Elements X=-0.06
    Return
SOLVE
   Frequency 50 to 3000 linstep 50
   return
import point format nastran file point.dat return
postrocess
   points all
   frequency 50 to 3000 linstep 50
   Near 2
   Far 5
   Quadrature 2 2 1
   Save Results Step 1
   Return
Return
combine
read point 176 pressure return
amplitude
read point 92 pressure return
amplitude
divide
amplitude
log10
constant real 20 imag 0
multiply
return

最后得出的传递损失曲线比理论值大了两倍左右。
急切需要得到帮助　，谢谢各位了。
附图中的结果是我除了2以后得出的。

小刚 发表于 2007-1-10 21:03

我的问题和您一样，而且我是用发动机进行实测的，结果更是大相径庭。如果你或你的朋友是这方面的高手不妨指导一下。谢了

ZHANGJIAWEI 发表于 2007-6-19 16:48

我粗略看了一下你的程序，好像没有加BOUNDARY CONDITIONS,要先在import and define tables for relation BCs中导入table

wrg 发表于 2007-6-20 09:35

cmd flie中几点建议：
FEM：1、最好Check mesh
         2、我认为你的进出口端的set应该不止一个face,（Face 54 Face 59），应该是一序列   
            的，网格不止于这么粗。
         3、 求解的步长也太大了 50Hz
         4、 TL 计算公式没有什么大问题，有个地方建议修改下Node 1 VelocityVector 0 1 0    ，
               直接Constant Real 1 Imag 0 就行,没必要那么复杂。
BEM：1、需要复制节点
         2、定义自由边的声压跳动
         3、步长太大
         4、TL计算公式最好是用你上面FEM中使用的，简单的声压比肯定是不准的。
如果有log file可能能给你更针对的方案，你先按上面的再多做几遍，看结果怎么样。不行的话，你把问题再描述下。

jianyuan 发表于 2007-7-21 18:04

定义自由边的声压跳动
这个定义为多少是不是0啊

jianyuan 发表于 2007-7-22 10:26

有限元搞定.用边界元怎么计算.我算的结果刚好是有限元的一半!!怎么回事
页: [1]
查看完整版本: 求助：消声器问题