一个声学仿真“问题”

captain007

小弟最近学习ＡＮＳＹＳ，想用ＡＮＳＹＳ做一个声学方面的仿真，可是求解完成后发现声压十分的有问题，附上命令流，斑竹和大大留步解答一下。里面的一些warning能不能解决？建模和加载有什么不妥？大家执行的时候可以把载荷步减小一点，网格化粗一些，因为很费时。各位大大和斑竹费心了。
/PREP7  
/TITLE,acoustic TL cal  
ET,1,solid45                             ! structural element   
et,2,fluid30                             !coupled interface
et,3,fluid130                            !infinite element  
r,3,1.9,0,0,0                            !real set for infinite element
et,4,fluid30,,1,0                        !general 3D acoustic element
SAVE
  
MP,EX,1,7e7                              !material property
MP,DENS,1,2700  
MP,NUXY,1,0.33  
MP,SONC,1,3080  
MP,DENS,2,1.21  
MP,SONC,2,344   
MP,MU,1,0   
SAVE

K,1,0.5,0.5,0                            !modal
K,2,0.5,-0.5,0  
K,3,-0.5,-0.5,0
K,4,-0.5,0.5,0

LSTR,       4,       1  
LSTR,       1,       2  
LSTR,       2,       3  
LSTR,       3,       4  
SAVE

FLST,2,4,4  
FITEM,2,1   
FITEM,2,4   
FITEM,2,3   
FITEM,2,2   
AL,P51X

FLST,2,1,5,ORDE,1   
FITEM,2,1   
VEXT,P51X, , ,0,0,0.0008,,,,
FLST,3,1,6,ORDE,1   
FITEM,3,1   
VGEN, ,P51X, , , , ,-0.0004, , ,1   
SAVE

CYL4,0,0,2, , , ,0.0008
FLST,3,1,6,ORDE,1   
FITEM,3,2   
VGEN, ,P51X, , , , ,-0.0004, , ,1   
SAVE

VOVLAP,1,2
SAVE

sphere,0,1.5,0,360
SAVE

SPHERE,1.9,1.5,0,360,   
SAVE  

FLST,2,4,6,ORDE,2   
FITEM,2,1   
FITEM,2,-4  
VOVLAP,P51X
SAVE

VSEL,S, , ,       1                                 !meshing volume 1
VPLOT
TYPE,   1   
MAT,       1
ESYS,       0   
SECNUM,
SAVE

FLST,5,1,4,ORDE,1   
FITEM,5,11  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
LESIZE,_Y1, , ,5, , , , ,1  
MSHAPE,0,3D
MSHKEY,1
CM,_Y,VOLU  
VSEL, , , ,       1
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y
VMESH,_Y1
SAVE

NSLA,S,1                                          !applying boundary to volume 1
NSEL,R,LOC,X,0.5
FLST,2,96,1,ORDE,10
FITEM,2,2   
FITEM,2,17  
FITEM,2,-31
FITEM,2,258
FITEM,2,273
FITEM,2,-287
FITEM,2,513
FITEM,2,-516
FITEM,2,577
FITEM,2,-636
D,P51X, , , , , ,ALL, , , , ,  
SAVE

NSLA,S,1
NSEL,R,LOC,X,-0.5   
FLST,2,96,1,ORDE,14
FITEM,2,1   
FITEM,2,32  
FITEM,2,47  
FITEM,2,-60
FITEM,2,257
FITEM,2,288
FITEM,2,303
FITEM,2,-316
FITEM,2,517
FITEM,2,-520
FITEM,2,637
FITEM,2,-640
FITEM,2,697
FITEM,2,-752
D,P51X, , , , , ,ALL, , , , ,
SAVE
   
NSLA,S,1
NSEL,R,LOC,Y,0.5
FLST,2,96,1,ORDE,6  
FITEM,2,1   
FITEM,2,-16
FITEM,2,257
FITEM,2,-272
FITEM,2,513
FITEM,2,-576
D,P51X, , , , , ,ALL, , , , ,
SAVE
   
NSLA,S,1
NSEL,R,LOC,Y,-0.5   
FLST,2,96,1,ORDE,10
FITEM,2,17  
FITEM,2,32  
FITEM,2,-46
FITEM,2,273
FITEM,2,288
FITEM,2,-302
FITEM,2,577
FITEM,2,-580
FITEM,2,637
FITEM,2,-696
D,P51X, , , , , ,ALL, , , , ,  
SAVE

FLST,5,2,6,ORDE,2                                    !select zhangban                        
FITEM,5,9   
FITEM,5,11  
VSEL,S, , ,P51X
VPLOT
SAVE

TYPE,   1                                            !mesh zhangban
MAT,       1
ESYS,       0   
SECNUM,
ESIZE,0.05,0,   
MSHAPE,1,3D
MSHKEY,0
FLST,5,2,6,ORDE,2   
FITEM,5,9   
FITEM,5,11  
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
VMESH,_Y1   
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
SAVE

FLST,5,2,6,ORDE,2                                  !boundaries of zhang ban  
FITEM,5,9   
FITEM,5,11  
VSEL,S, , ,P51X
VPLOT   
FLST,2,8915,1,ORDE,8
FITEM,2,1   
FITEM,2,-60
FITEM,2,257
FITEM,2,-316
FITEM,2,513
FITEM,2,-752
FITEM,2,1537
FITEM,2,-10091  
D,P51X, , , , , ,ALL, , , , ,
SAVE

FLST,5,4,6,ORDE,2                                  !spacial acoustic et  
FITEM,5,5   
FITEM,5,-8  
VSEL,S, , ,P51X
VPLOT   
SAVE

TYPE,   2   
MAT,       2
ESYS,       0   
SECNUM,
ESIZE,0.05,0,
MSHAPE,1,3D
MSHKEY,0
FLST,5,4,6,ORDE,2   
FITEM,5,5   
FITEM,5,-8  
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
VMESH,_Y1
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
SAVE

esel,s,type,,1                                    !modify acoustic et
nsle,s
esln,s,0
nsle,s
esel,inve
nsle,s
emodify,all,type,4
alls
SAVE

ASEL,S, , ,      23                               !infinite boundary
nsla,s,1
type,3  
real,3  
mat,2   
esurf   
SAVE

ASEL,S, , ,      24
nsla,s,1
type,3  
real,3  
mat,2   
esurf   
SAVE

ASEL,S, , ,      26
nsla,s,1
type,3  
real,3  
mat,2   
esurf   
SAVE

ASEL,S, , ,      27
nsla,s,1
type,3  
real,3  
mat,2   
esurf   
SAVE

VSEL,S, , ,       6                                !apply loads on surface
VPLOT   
FLST,2,2,5,ORDE,2   
FITEM,2,20  
FITEM,2,-21
SFA,P51X,1,PRES,
SAVE

asel,s,,,1                                         !mark coulping surface
asel,a,,,2
asel,a,,,19
asel,a,,,22
asel,a,,,25
asel,a,,,28
NSLA,S,1
esel,s,type,,2
sf,all,fsi,1
alls
SAVE

ANTYPE,3                                           !loader step
HROPT,FULL  
HROUT,ON
LUMPM,0
EQSLV,FRONT,1e-008,
PSTRES,0
FLST,2,1208,1,ORDE,2
FITEM,2,1537
FITEM,2,-2744
SAVE

HARFRQ,0,500,   
NSUBST,5,
KBC,1   
SAVE

GPLOT                                             !choose all and save all  
esel,all
nsel,all
SAVE

finish  

/SOLU                                             !solve
/STATUS,SOLU
SOLVE