有限元流体模型和直接边界元分析问题

榫眼

模型参数如下：
矩形房间,x向为6m,y向为5m,z向为4m,顶面铺设吸声材料,界面声阻抗为37ρc,其他墙面为刚性壁，声速340，空气密度为1.21，
分别用有限元和边界元(BEM Direct Collocation Node Interior Uncoupled Frequency)求解声压级响应，在房间内a（2.0，1.5，1.0）点布置单位点声源，激励频率为30hz，b（4.0，1.5，2.5）和c（3.0，1.5，2.0）分别为接收点
1.边界元法
New Name 'bem_2' Model 1 File bem_2.sdb Return
Option
BEM Direct Collocation Node Interior Uncoupled Frequency
Return
Import Mesh Format Patran File damp_sen\myroom.out Return
Check Mesh Return
Material Fluid Name 'air' Sound Real 3.4000e+002Rho Real 1.21
Return
Boundary ImpedanceReal 15221.8 Imag 0
Nodes BetweenZ= 4,4
Return
Source Name 'p' Spherical
Amplitude Real 1 Imag 0
Position 2.0 1.5 1.0
Return
Point 4.0 1.5 2.5 Return
Point 3.0 1.5 2.0 Return
Solve
Frequency 1 To 100 LinStep 1
Return
PostProcess
Points All
Frequency 1 To 100 LinStep 1
Near 2 Far 5 Quadrature 2 2 1
Save Results Step 1
Return
Store
Case Active
Results
Frequency 1 To 100 LinStep 1
Points All
Return
Return


[ 本帖最后由 榫眼 于 2008-6-18 17:45 编辑 ]

榫眼

2.有限元法（FEM Frequency Fluid）（模态法求解） 节点330（2.0，1.5，1.0）为声源；
  节点763 (4.0，1.5，2.5）为b点，节点618（3.0，1.5，2.0）为c点
New Name 'fem_001' Model 1 File fem_001.sdb Return
Option
FEM Frequency Fluid Return
Import Mesh Format Patran File damp_sen\room.out Return
Material Fluid Name 'air' Sound Real 3.4000e+002 Rho Real 1.210e+000
Elements All
Return
Modes
Vector 10 Tolerance 1.0000e-6 Iteration 100
BLanczos
Return
Save Return
Set1 Name"a" FacesZ= 4
Return
Boundary ImpedanceReal 15221.8 Imag 0
Faces Set 1
Return
Set Name
"Envelope" Envelope Elements All
Return
Set 3 Name"g" Difference2 1
Return
Boundary VelocityReal 0 Imag 0 Faces Set 3
Return
Source Name 'p' Spherical Amplitude Real 10 Imag 0 Node 330
Return
Solve
Frequency 1 To 100 LinStep 1
Return
Store Case Active Potentials
Frequency 1 To 100 LinStep 1
Nodes 763
618
Return
Return

这两种方法得到的结果却完全不一样？（其中边界元的结果为正确的结果）
我不知道有限元方法中那处错了，请大家帮我看看

[ 本帖最后由 榫眼 于 2008-6-18 17:47 编辑 ]

榫眼

下面是模型文件，去掉txt后缀，变为rar压缩名

pengweicai

有限元计算 加声源的模型好像是有些问题的。