求助：消声器问题

djh713

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

FEM轴对称：

Option  FEM 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+002  Rho 1.2250e+000  
    Elements All  
    Return  
Boundary Velocity  Real -1 Imag 0  
    Faces Set 1  
    Return  
Boundary Admittance  Real 2.401e-3 Imag 0  
    Faces Set 2  
    Return  
Solve  Frequency 50 To 3000 LinStep 50 Return  
Combine  
Read Node 11 Pressure Return  
Amplitude  
Read Node 1 Pressure Return  
Read Node 1 Velocity  Vector 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+02  Rho Real 1.2250e+00  
    Return  
Boundary Velocity  Real 1 Imag 0  
    Elements X=0.06
    Return  
Boundary Impedance  Real 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以后得出的。

小刚

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

ZHANGJIAWEI

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

wrg

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

jianyuan

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

jianyuan

有限元搞定.用边界元怎么计算.我算的结果刚好是有限元的一半!!怎么回事