看看这两种做法为什么会有这样的结果差别？

fxdl · 发表于 2009-5-18 21:47

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

这是消音器的那个例子，用两种方法求传递损失，第一种按照软件自带的步骤做的，得到的结果，
第二种是按照网上有个兄弟的处理方法做的，两种方法结果不一样，为什么呢？

登录/注册后可看大图

上面是按照消音器的例子一步不差的做出来的结果。
具体步骤:
Environment Section SETUP BTNFILE 'muffler.BTN'
Return
Environment BTNFILE=mufflerNOFLOW.BTN Return
New Name 'noflow' Model 1 File noflow.sdb Return
Option  FEM Frequency Fluid Return
Import Mesh Format Free File muffler.fre Return
Material Fluid
Name 'air'
Sound 3.4000e+02  Rho 1.2250e+00
Elements all
Return
Set 10 Name  "inlet"
Faces Z=-.08
Return
Boundary Velocity  Real -1 Imag 0
Faces Set 10
Return
Set 20 Name  "outlet"
Faces Z=0.6
Return
Table 1 Name 'admittance' File admittance.tbl Return
Boundary Admittance Table 1
Faces Set 20
Return
Solve  Frequency 20 To 300 LinStep 2 Return
Combine
Read Node 925 Pressure Return
Read Node 925 Pressure Return
Conjugate
Multiply
Root
Read Node 3244 Pressure Return
Read Node 3244 Pressure Return
Conjugate
Multiply
Root
Divide
Real
Log10
Constant Real 20 Imag 0
Multiply
Write Table 2 Name 'attenuation' Return
Return

下面的图是按照这个步骤做的。
Option  FEM Frequency Fluid Return
Import Mesh Format Free File muffler.fre Return
AxiSymmetry Return
Set 1 Name "inlet"
Face Z=-0.08
Return
Set 2 Name "outlet"
Face Z=0.6
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 20 To 300 LinStep 2 Return
Combine
Read Node 925Pressure Return
Amplitude
Read Node 3244Pressure Return
Read Node 3244 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

为什么结果不太一样呢？差在哪了呢？
另外在下面的步骤中设定导纳边界条件Boundary Admittance  Real 2.401e-3 Imag 0
其中为什么在Real 中设置2.401e-3 ？怎么来的这个数据？

账号		自动登录	找回密码
密码			我要加入

[综合讨论] 看看这两种做法为什么会有这样的结果差别？

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。