suma1989 发表于 2012-5-21 22:28

室内声场仿真求助

背景:长方体房间,中心位置放置一个平板,通过施加一定的频率使其振动,然后在气体中产生噪声,最后要求得到声压分布图
已完成,网格划分,但是流固接触,边界条件,施加频率载荷都还不完全,得不到正确的声压分布图,大家帮忙看下嘛,谢谢
/PREP7/TITLE,Room Acoustic Analysis
ET,1,SOLID45                   !壳单元 类型1
ET,2,FLUID30                   ! acoustic fluid element with ux & uy(0-Structure present at interface (unsymmetric element matrix))
ET,3,130                     ! acoustic infinite line element
!r,3,0.3,0,0
ET,4,FLUID30,,1,0            ! acoustic fluid element without ux & uy(1-No structure at the interface (symmetric element matrix))
! material properties
!材料1为铁
MP,EX,1,2.1E11
MP,DENS,1,7800
MP,NUXY,1,0.3

!材料2为空气介质
MP,DENS,2,1.21
MP,SONC,2,344


MP,DENS,3,2.35E-3   ! (Rho) density of air (lb/ft**3)
MP,SONC,3,1100.0    ! (C) speed of sound in air (ft/sec)
MP,MU,3,0.04      ! (Beta) absorption coefficient of the walls
                  ! ('Beta' should be between 0 - 1)

!building model
wprot,0,-90,0
block,-6,6,-4,4,-2,2
wprot,0,90,0
BLOCK,-0.2,0.2, -0.002, 0.002, -0.1, 0.1    !BLOCK, X1, X2, Y1, Y2, Z1, Z2

!meshing
!each of the length X 100
lesize,14,,,40,1
lesize,16,,,40,1
lesize,17,,,40,1
lesize,19,,,40,1
lesize,21,,,20,1
lesize,22,,,20,1
lesize,23,,,20,1
lesize,24,,,20,1
lesize,13,,,4,1
lesize,15,,,4,1
lesize,18,,,4,1
lesize,20,,,4,1
VSEL,s,,,2!select entity of number 2
VATT,1,,1   !VATT,MAT, REAL, TYPE, ESYS, SECNUM
mshkey,1    !Specifies whether free meshing or mapped meshing should be used to mesh a model.1— Use mapped meshing.
VMESH,2
alls

lesize,2,,,50,1
lesize,4,,,50,1
lesize,5,,,50,1
lesize,7,,,50,1
lesize,1,,,30,1
lesize,3,,,30,1
lesize,6,,,30,1
lesize,8,,,30,1
lesize,9,,,10,1
lesize,10,,,10,1
lesize,11,,,10,1
lesize,12,,,10,1
VSEL,s,,,1
VATT,2,,2   !指定的材料属性为空气介质(与结构耦合)
mshkey,1
VMESH,1
alls

! Boundary conditions
alls
nsel,s,loc,X,6
nsel,a,loc,X,-6
nsel,a,loc,Y,4
nsel,a,loc,Y,-4
nsel,a,loc,Z,2
nsel,a,loc,Z,-2
D,ALL,UX,,,,,UY,UZ! Constrain all displacements to zero at the walls
SF,ALL,IMPD,1

!建立流-固接触面
alls
Asel,s,,,7
Asel,A,,,9
Asel,A,,,10
Asel,A,,,11
Asel,A,,,8
Asel,A,,,12
NSLA,S,1
esel,s,type,,2
sf,all,fsi,1
alls

/solu
antype,harmic
HARF,80,80
SOLVE

耿庆涛 发表于 2012-10-6 13:31

楼主,我也做声场分析 ,加我qq吧:545621944,咱们可以多交流。
页: [1]
查看完整版本: 室内声场仿真求助