声学噪声模拟问题
我想做一个发动机加隔声罩后看隔声罩一米远处的噪声情况。噪声源（发动机）采用直径为0.3米的实心圆代替，采用solid45单元，其属性设置为：mp,ex,1,2.10e11
mp,dens,1,7600
mp,nuxy,1,0
mp,sonc,1,4800
在实心圆外加一个厚度为0.03米空心圆（空气），在两圆的接触面上施加声固耦合和声压载荷：
选取接触面上的所有点，我用的是下面的命令施加的：
SF,ALL,FSI 声固耦合
D,ALL,PRES,2.12 声压载荷（声压级在100dB左右）
0.03米的空心圆属性设置为：
et,2,fluid30
mp,dens,2,1.21
mp,sonc,2,344
R,2,20e-6
mp,mu,2,0.02
外面就是一个内空的长方体作为吸声结构，其属性设置为：
et,3,fluid30,,1,0
mp,ex,3,1e9
mp,nuxy,3,0.2
mp,dens,3,30
mp,sonc,3,800
mp,mu,3,0.9
R,3,20e-6
在最外面用一个球体来模拟无限远处的吸声面，球体内部设置为空气，et,3,fluid30,,1,0
在球体最外表面，选取表面上的所有点，设置属性为：
et,4,fluid130
mp,dens,2,1.21
mp,sonc,2,344
mp,mu,2,0.02
R,4,4（球体直径为4米）
运用ESURF命令生成无限远处吸声表面。
我大体上是这么做的，可是结果太离谱了，我查看点的计算结果数值有几百的，我输入的才是2.12啊，不知道问题出在哪，请各位大侠帮忙看看，有什么我做的不对的。
还有就是我看声压值有实部虚部之分，这个都怎么进行查看啊？最终要的结构怎么生成啊？我很疑惑，我是个初学者，请大家多多关照
/PREP7
SPH4,0,0,4 模型结构建立
BLOCK,1.7,-1.7,1.4,-1.4,1.4,-1.4,
/VIEW,1,1,2,3
/ANG,1
/REP,FAST
LPLOT
FLST,2,2,6,ORDE,2
FITEM,2,1
FITEM,2,-2
VOVLAP,P51X
LPLOT
BLOCK,0.7,-0.7,0.4,-0.4,0.4,-0.4,
FLST,2,2,6,ORDE,2
FITEM,2,1
FITEM,2,-2
VOVLAP,P51X
LPLOT
BLOCK,0.65,-0.65,0.35,-0.35,0.35,-0.35,
FLST,2,2,6,ORDE,2
FITEM,2,1
FITEM,2,-2
VOVLAP,P51X
LPLOT
BLOCK,-0.6,-0.47,0.08,0.14,0.35,0.38,
FLST,2,2,6,ORDE,2
FITEM,2,1
FITEM,2,5
VOVLAP,P51X
LPLOT
BLOCK,0.6,0.47,0.08,0.14,0.35,0.38,
FLST,2,2,6,ORDE,2
FITEM,2,5
FITEM,2,-6
VOVLAP,P51X
LPLOT
BLOCK,0.6,0.47,0.08,-0.125,0.365,0.38,
BLOCK,-0.6,-0.47,0.08,-0.125,0.365,0.38,
FLST,2,3,6,ORDE,2
FITEM,2,6
FITEM,2,-8
VOVLAP,P51X
LPLOT
FLST,2,2,6,ORDE,2
FITEM,2,1
FITEM,2,9
VADD,P51X
FLST,2,2,6,ORDE,2
FITEM,2,5
FITEM,2,-6
VADD,P51X
LPLOT
BLOCK,-0.6,-0.47,-0.11,-0.125,0.4,0.38,
BLOCK,-0.6,-0.47,-0.085,-0.1,0.4,0.38,
BLOCK,-0.6,-0.47,-0.075,-0.06,0.4,0.38,
BLOCK,-0.6,-0.47,-0.035,-0.05,0.4,0.38,
BLOCK,-0.6,-0.47,-0.025,-0.01,0.4,0.38,
BLOCK,0.6,0.47,-0.025,-0.01,0.4,0.38,
BLOCK,0.6,0.47,-0.035,-0.05,0.4,0.38,
BLOCK,0.6,0.47,-0.06,-0.075,0.4,0.38,
BLOCK,0.6,0.47,-0.1,-0.085,0.4,0.38,
BLOCK,0.6,0.47,-0.11,-0.125,0.4,0.38,
FLST,2,11,6,ORDE,4
FITEM,2,5
FITEM,2,-6
FITEM,2,8
FITEM,2,-16
VOVLAP,P51X
FLST,2,6,6,ORDE,3
FITEM,2,5
FITEM,2,-9
FITEM,2,11
VADD,P51X
FLST,2,6,6,ORDE,4
FITEM,2,1
FITEM,2,12
FITEM,2,-15
FITEM,2,17
VADD,P51X
FLST,2,3,6,ORDE,3
FITEM,2,2
FITEM,2,5
FITEM,2,10
VADD,P51X
VSEL,S, , , 1
VPLOT
LPLOT
alls
SPH4,0,0,0.32
FLST,2,2,6,ORDE,2
FITEM,2,1
FITEM,2,-2
VOVLAP,P51X
LPLOT
SPH4,0,0,0.3
FLST,2,2,6,ORDE,2
FITEM,2,1
FITEM,2,-2
VOVLAP,P51X
LPLOT
NUMMRG,ALL, , , ,LOW
BLOCK,0.7,0.67,0.1,0.3,0.2,-0.2,
FLST,2,2,6,ORDE,2
FITEM,2,2
FITEM,2,4
VADD,P51X
LPLOT
NUMMRG,ALL, , , ,LOW
BLOCK,0.62,0.65,0.07,0.33,0.23,-0.23,
FLST,2,2,6,ORDE,2
FITEM,2,2
FITEM,2,5
VOVLAP,P51X
LPLOT
FLST,2,2,6,ORDE,2
FITEM,2,2
FITEM,2,7
VADD,P51X
LPLOT
SAVE
NUMMRG,ALL, , , ,LOW
et,1,solid45 ：单元定义
et,2,fluid30
et,3,fluid30,,1,0
et,4,fluid130
mp,ex,1,2.10e11 ：噪声源属性（钢结构）
mp,dens,1,7600
mp,nuxy,1,0
mp,sonc,1,4800
CM,_Y,VOLU
VSEL, , , , 1
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 1, , 1, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
mp,dens,2,1.21 ：空气属性
mp,sonc,2,344
R,2,20e-6
mp,mu,2,0.02
：吸声材料属性
mp,dens,3,30
mp,sonc,3,800
mp,mu,3,0.9
R,3,20e-6
CM,_Y,VOLU
VSEL, , , , 6
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 2, 2, 2, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
CM,_Y,VOLU
VSEL, , , , 4
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 2, 2, 3, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
CM,_Y,VOLU
VSEL, , , , 5
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 3, 3, 3, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
CM,_Y,VOLU
VSEL, , , , 8
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 2, 2, 3, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
CM,_Y,VOLU
VSEL, , , , 3
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 2, 2, 3, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
esize,0.05
MSHAPE,1,3D
MSHKEY,0
!*
CM,_Y,VOLU
VSEL, , , , 1
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
LPLOT
CM,_Y,VOLU
VSEL, , , , 6
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
LPLOT
NUMMRG,ALL, , , ,LOW
CM,_Y,VOLU
VSEL, , , , 4
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
LPLOT
NUMMRG,ALL, , , ,LOW
CM,_Y,VOLU
VSEL, , , , 5
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
LPLOT
esize,0.15
CM,_Y,VOLU
VSEL, , , , 8
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
LPLOT
NUMMRG,ALL, , , ,LOW
esize,0.3
CM,_Y,VOLU
VSEL, , , , 3
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
LPLOT
FLST,5,2,5,ORDE,2
FITEM,5,23
FITEM,5,25
ASEL,S, , ,P51X
APLOT
NSLA,S,1
NPLOT
FINISH
/SOL
!*
ANTYPE,3
FLST,5,516,1,ORDE,2
FITEM,5,1
FITEM,5,-516
CM,_Y,NODE
NSEL,R, , ,P51X
CM,_Y1,NODE
CMSEL,S,_Y
CMDELE,_Y
!*
/GO
!*
D,_Y1,PRES,2.12 在钢球（噪声源表面施加2.12pa的声压）
!*
CMDELE,_Y1
!*
FLST,2,516,1,ORDE,2
FITEM,2,1
FITEM,2,-516
SF,P51X,FSI 在钢球（噪声源）与空气接触表面施加声固耦合（只在这里施加了声固耦合）
alls
LPLOT
FLST,5,2,5,ORDE,2
FITEM,5,1
FITEM,5,-2
ASEL,S, , ,P51X
NSLA,S,1
NPLOT
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,1.08222638492,1
/REP,FAST
FINISH
/PREP7
R,4,4,,,,
!*
TYPE, 4
MAT, 2
REAL, 4
ESYS, 0
SECNUM,
!*
esurf
LPLOT
alls
SAVE
FINISH
/SOL
HARFRQ,2000,2000,
NSUBST, ,
KBC,0
!*
/STATUS,SOLU
SOLVE
[ 本帖最后由 xuruikl 于 2009-2-26 13:58 编辑 ]
页:
[1]