请看看电压分布对不对

lchylinchy

有一个压电晶片，放在空气中，在晶片上施加一声压，需要求得压电晶片上的电压，麻烦帮我看看对不对，附录是电压分布图
fini
/clear
/prep7
/title,quartz vibration anslysis
!et,1,PLANE13,7  !结构单元类型
et,1,PLANE223,1001
et,2,29         !流体单元类型与结构接触
et,3,129   !吸收边界
et,4,29,,1,0      !流体单元类型与结构没有接触 在xyz方向没有位移
r,3,6e-3,0,0
!press=1e-8     !力的大小
E1=2.0e9                        ! Young's modulus, N/m^2
NU12=0.29                       ! Poisson's ratio
G12=0.775e9                     ! Shear modulus, N/m^2
d31=2.2e-11                     ! Piezoelectric strain coefficients, C/N
d32=0.3e-11
d33=-3.0e-11
ept33=12                        ! Relative permittivity at constant stress
tb,ANEL,1,,,1                 ! Elastic compliance matrix（tb定义非线性材料单元表，同tbdata一起连用）
!最后一个1表示输入的弹性矩阵为弹性形式，如果是0表示输入的弹性矩阵为刚度形式
tbda,1,1/E1,-NU12/E1,-NU12/E1
tbda,7,1/E1,-NU12/E1
tbda,12,1/E1
tbda,16,1/G12
tbdata,19,2.09e10  
tbdata,21,2.09e10
tb,PIEZ,1,,,1                  ! Piezoelectric strain matrix
tbda,2,d31
tbda,5,d33
tbda,8,d32
tb,DPER,1,,,1                   ! Permittivity（介电常数） at constant stress
tbdata,1,ept33,ept33
!tblist,all
mp,dens,1,2650        !石英密度
mp,dens,2,2.35E-3    !空气密度
MP,SONC,2,1100       !空气声速
mp,mu,2,0.75           !石英英叉的吸收系数
mp,dens,3,2.35E-3    !空气密度
MP,SONC,3,1100       !空气声速
!结构模型
blc4,-0.75e-3,-2.7e-3,1.5e-3,1.8e-3
blc4,-0.75e-3,-0.9e-3,0.6e-3,3.6e-3
blc4,0.15e-3,-0.9e-3,0.6e-3,3.6e-3
aadd,all
cyl4,0,0,6e-3,0,,360
aovlap,1,4
!esize,0.2e-3,0,  !指定单元划分大小
lesize,1,,,15
lesize,15,,,3
lesize,7,,,6
lesize,11,,,6
lesize,4,,,18
lesize,2,,,18
lesize,14,,,36
lesize,13,,,36
lesize,12,,,36
lesize,6,,,36
mshape,0,2d
mshkey,0
!smrt,1
asel,s,,,4
aatt,1,,1
amesh,all
asel,s,,,2
aatt,2,,2
smrt,1
mshkey,0
mshape,0,2d
amesh,all
save,pingm122_mesh,db
!将不与结构相连的单元类型改为4
esel,s,type,,1
nsle,s
esln,s,0
nsle,s
esel,inve
nsle,s
emodif,all,type,4
esel,all
nsel,all
!定义吸收边界
csys,1
nsel,s,loc,x,6e-3
type,3
real,3
mat,2
esurf
esel,all
nsel,all
csys,0
!压电耦合
!定义上层电极
lsel,s,,,6,12,6
NSLL,S,1
cp,1,volt,all
*get,n_top,node,0,num,min
!d,n_top,volt,0
esel,all
nsel,all
!定义下层电极
nsel,s,loc,x,-0.75e-3
nsel,r,loc,y,-2.7e-3,2.7e-3
cp,2,volt,all
*get,n_bot1,node,0,num,min
d,n_bot1,volt,0
esel,all
nsel,all
nsel,s,loc,x,0.75e-3
nsel,r,loc,y,-2.7e-3,2.7e-3
cp,3,volt,all
*get,n_bot2,node,0,num,min
d,n_bot2,volt,0
esel,all
nsel,all
!定义流固界面
csys,0
nsel,s,loc,x,-0.75e-3,0.75e-3
nsel,r,loc,y,-2.7e-3
esel,s,type,,2
sf,all,fsi,1i
nsel,s,loc,x,-0.75e-3
nsel,r,loc,y,-2.7e-3,2.7e-3
esel,s,type,,2
sf,all,fsi,1
nsel,s,loc,x,0.75e-3
nsel,r,loc,y,-2.7e-3,2.7e-3
esel,s,type,,2
sf,all,fsi,1
Lsel,s,,,6,7,1
Lsel,a,,,11,12,1
Lsel,a,,,15
NSLL,S,1
esel,s,type,,2
sf,all,fsi,1
esel,all
nsel,all

!边界条件
nsel,s,loc,x,-0.75e-3,0.75e-3
nsel,r,loc,y,-2.7e-3
d,all,ux,0,,,uy
alls
fini
/solu
antype,harmic
hropt,full
sfl,6,pres,1
sfl,12,pres,1
harfrq,1000,1000
NSUBST,5
KBC,1
solv
fini
/post1
set,first
plnsol,press
plesol,volt
prnsol,volt
prnsol,press