大家帮忙看看,比较mu＝1单元和无限吸收单元的例子

melanie

根据网上一个比较mu＝1单元和无限吸收单元的例子做的。应该是蛮简单的例子，但怎么也得不到和例子一样的图。不知道问题在哪里呢 大家帮忙看看吧。
这是我做的图，衰减也大了，有40db。这是网上的例子，正确的应该做成这样吧，可我怎么都做不出来。
我用GUI做的，命令流有点乱，大家将就看看吧。抱歉抱歉～～～
/PREP7  
!*  
ET,1,FLUID29         单元设置
!*  
KEYOPT,1,2,0
KEYOPT,1,3,0
!*  
ET,2,FLUID29
!*  
KEYOPT,1,2,1
KEYOPT,1,3,0
!*  
ETDEL,2
!*  
R,1,0.00002,       参考声压
!*  
!*  
MPTEMP,,,,,,,,               材料属性
MPTEMP,1,0  
MPDATA,DENS,1,,1.22
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,SONC,1,,344  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,MU,1,,0  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,2,,1.22
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,SONC,2,,344  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,MU,2,,1  
K,1,0.5,0.5,,                  建立模型
K,2,0.5,-0.5,,  
K,3,-0.5,-0.5,,
K,4,-0.5,0.5,,  
LSTR,       4,       1  
LSTR,       1,       2  
LSTR,       2,       3  
LSTR,       3,       4  
FLST,2,4,4  
FITEM,2,1   
FITEM,2,2   
FITEM,2,3   
FITEM,2,4   
AL,P51X
TYPE,   1   
MAT,       1
REAL,       1   
ESYS,       0   
SECNUM,
!*  
FLST,5,4,4,ORDE,2   
FITEM,5,1   
FITEM,5,-4  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,20, , , , ,1   网格划分
!*  
MSHAPE,0,2D
MSHKEY,0
!*  
CM,_Y,AREA  
ASEL, , , ,       1
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
FLST,5,4,4,ORDE,2   
FITEM,5,1   
FITEM,5,-4  
LSEL,S, , ,P51X
NSEL,ALL
NSLL,S,1
ESLN,S  
FLST,2,76,2,ORDE,38
FITEM,2,1   
FITEM,2,-21
FITEM,2,40  
FITEM,2,-41
FITEM,2,60  
FITEM,2,-61
FITEM,2,80  
FITEM,2,-81
FITEM,2,100
FITEM,2,-101
FITEM,2,120
FITEM,2,-121
FITEM,2,140
FITEM,2,-141
FITEM,2,160
FITEM,2,-161
FITEM,2,180
FITEM,2,-181
FITEM,2,200
FITEM,2,-201
FITEM,2,220
FITEM,2,-221
FITEM,2,240
FITEM,2,-241
FITEM,2,260
FITEM,2,-261
FITEM,2,280
FITEM,2,-281
FITEM,2,300
FITEM,2,-301
FITEM,2,320
FITEM,2,-321
FITEM,2,340
FITEM,2,-341
FITEM,2,360
FITEM,2,-361
FITEM,2,380
FITEM,2,-400
EMODIF,P51X,MAT,2,  
/PNUM,KP,0  
/PNUM,LINE,0
/PNUM,AREA,0
/PNUM,VOLU,0
/PNUM,NODE,0
/PNUM,TABN,0
/PNUM,SVAL,0
/NUMBER,0   
!*  
/PNUM,MAT,1
/REPLOT
!*  
NSEL,ALL
ESEL,ALL
/PNUM,KP,0  
/PNUM,LINE,0
/PNUM,AREA,0
/PNUM,VOLU,0
/PNUM,NODE,0
/PNUM,TABN,0
/PNUM,SVAL,0
/NUMBER,0   
!*  
/PNUM,MAT,1
/REPLOT
!*  
FLST,5,1,1,ORDE,1   
FITEM,5,261
CM,_Y,NODE  
NSEL,R, , ,P51X
CM,_Y1,NODE
CMSEL,S,_Y  
CMDELE,_Y   
!*  
/GO
!*  
D,_Y1,PRES,1
!*  
CMDELE,_Y1  
!*  
FINISH  
/SOL
!*  
!*  
ANTYPE,3
HARFRQ,300,300,           频率
NSUBST,1,   
KBC,0   
!*  
/STATUS,SOLU
SOLVE   
FINISH  
/POST1  
SET,LAST
!*  
PLESOL, NMISC,4, 0        查看声压级

[ 本帖最后由 perch88 于 2008-9-3 00:20 编辑 ]