请大侠们看看这个声场仿真，有代码和图片

lpc8584

我做了一下压电陶瓷的声场仿真，但是声场分布明显有问题，麻烦各位老师们看一下：
/PREP7
/TITLE

S11=12.3E-12
S12=-4.05E-12
S13=-5.13E-12
S33=15.5E-12
S44=39E-12
S66=32.7E-12

D15=496E-12
D31=-123E-12
D33=289E-12

EP11=370
EP33=635

RHO=7500
!*******************压电陶瓷参数定义
!定义柔顺矩阵
TB,ANEL,1,,,1
TBDATA,1,S11,S13,S12
TBDATA,7,S33,S13
TBDATA,12,S11
TBDATA,16,S44
TBDATA,19,S44
TBDATA,21,S66
!定义压电矩阵
TB,PIEZ,1,,,1
TBDATA,2,D31
TBDATA,5,D33
TBDATA,8,D31
TBDATA,10,D15
TBDATA,15,D15
!定义介电常数
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,PERX,1,,EP11
MPDATA,PERY,1,,EP33
MPDATA,PERZ,1,,EP11
!定义密度和阻尼系数
MP,DENS,1,RHO
MP,DMPR,1,0.0042
!定义材料
ET,1,PLANE223,1001,,1
ET,2,FLUID29,,0,1
ET,3,129,,,1
R,3,0.05,0.05,0
ET,4,FLUID29,,1,1
!****************空气参数定义
MP,DENS,2,1.21
MP,SONC,2,344
!****************建模
R=6.25E-3
T=1E-3
RL=0.05

CYL4,RL,0,0,0,RL,360

RECT,RL-R,RL+R,-T/2,T/2

AOVLAP,1,2

ASEL,S,AREA,,2
AATT,1,1,1,0
LESIZE,5,,,8,1
LESIZE,6,,,2,1
LESIZE,7,,,8,1
LESIZE,8,,,2,1
MSHKEY,0
MSHAPE,0,2D
AMESH,2

ASEL,S,AREA,,3
AATT,2,1,2,0
LESIZE,1,,,80,1
LESIZE,2,,,80,1
LESIZE,3,,,80,1
LESIZE,4,,,80,1
MSHKEY,0
MSHAPE,0,2D
AMESH,3

ESEL,S,TYPE,,1
NSLE,S
ESLN,S,0
NSLE,S
ESEL,INVE
NSLE,S
EMODIF,ALL,TYPE,4
ESEL,ALL
NSEL,ALL

local,11,1,0.05,0,0
CSYS,11
NSEL,S,LOC,X,RL
TYPE,3
REAL,3
MAT,2
ESURF
ESEL,ALL
NSEL,ALL

CSYS,0

LSEL,S,,,5
LSEL,A,,,6
LSEL,A,,,7
LSEL,A,,,8
SFL,ALL,FSI
NSEL,ALL
ESEL,ALL

LSEL,S,,,5
NSLL,S,1
CP,1,VOLT,ALL
*GET,N_TOP,NODE,,NUM,MIN
D,ALL,VOLT,1


ALLSEL,ALL

LSEL,S,,,7
NSLL,S,1
CP,2,VOLT,ALL
*GET,N_BOT,NODE,,NUM,MIN
D,ALL,VOLT,0

allsel
FINISH


/SOLU
EQSLV,SPARSE
ANTYPE,HARMIC
HROPT,FULL
OUTERS,ALL,ALL
HARFRQ,77000,80000
NSUBST,10
KBC,1
SOLVE

/POST1
SET,FIRST
NSEL,ALL
PLNSOL,PRES