有限元流体模型和直接边界元分析问题
模型参数如下:矩形房间,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压缩名 有限元计算 加声源的模型好像是有些问题的。
页:
[1]