问一个关于施加耦合场的问题

cyf7895123

上图为一个双片的悬臂梁压电振子，我应该如何施加耦合场做分析，才不会使我的计算出错。我的命令流如下，
finish
/clear
/filnam,yadianliang
/title,static ansys
/prep7
mp,DENS,1,7600                 
mp,perx,1,7.124e-9               
mp,pery,1,7.124e-9
mp,perz,1,5.841e-9

tb,ANEL,1                  
tbdata,1,13.2E10,7.1E10,7.3E10  
tbdata,7,13.2E10,7.3E10      
tbdata,12,11.5E10            
tbdata,16,3E10            
tbdata,19,2.6E10              
tbdata,21,2.6E10            

tb,PIEZ,1                    
tbdata,3,-4.1            
tbdata,6,-4.1              
tbdata,9,14.1              
tbdata,14,10.5              
tbdata,16,10.5              

mp,dens,2,8920
EX,2,10.6e10
nuxy,2,0.35

mp,dens,3,390
EX,3,9.8e8
nuxy,3,0.49

L=50e-3
W=20e-3
H =0.2e-3   

L1=70e-3
H1=0.25e-3

H2=5e-3

local,11               
local,12,,,,,,,180                        
csys,11                        
block,0,L1,0,W,O,H1
block,0.01,0.01+L,0,W,0,-H
block,0.01,0.01+L,0,W,H1,H1+H

block,0,L1,0,W,0.01,0.01+H1
block,0.01,0.01+L,0,W,0.01,0.01-H
block,0.01,0.01+L,0,W,0.01+H1,0.01+H1+H

block,L1,L1+0.003,0,W,-0.01,0.02

vglue,all                  
et,1,solid5,3            
et,2,solid45
et,3,solid45

asel,s,loc,z,0
cm,CM_1,area
cmplot,CM_1
accat,CM_1

asel,s,loc,z,H1
cm,CM_2,area
accat,CM_2

LESIZE,ALL,1e-3, , , ,1, , ,1,
mat,1 $ type,1 $ esys,11
Vsweep,2
mat,2 $ type,2 $ esys,11
Vsweep,8
mat,1 $ type,1 $ esys,12
Vsweep,3

mat,1 $ type,1 $ esys,11
Vsweep,5
mat,2 $ type,2 $ esys,11
Vsweep,10
mat,1 $ type,1 $ esys,12
Vsweep,6

mat,3 $ type,3 $ esys,11
Vsweep,9

nsel,s,loc,z,0
nsel,r,loc,x,0.01,0.01+L           !定义下层电极
cp,1,volt,all
*get,n_bot,node,0,num,min
nsel,s,loc,z,H+H1         !定义上层电极
cp,2,volt,all
*get,n_top,node,0,num,min
nsel,s,loc,z,0            !压电片中间面电压耦合
nsel,r,loc,x,0.01,0.01+L
cm,CM_3,node
nsel,s,loc,z,H1
nsel,r,loc,x,0.01,0.01+L
cm,Cm_4,node
cmsel,s,cm_3,node
cmsel,a,cm_4,node
cp,3,volt,all

nsel,s,loc,z,0.01
nsel,r,loc,x,0.01,0.01+L           !定义下层电极
cp,4,volt,all
*get,n_bot,node,0,num,min
nsel,s,loc,z,0.01+H+H1         !定义上层电极
cp,5,volt,all
*get,n_top,node,0,num,min
nsel,s,loc,z,0.01            !压电片中间面电压耦合
nsel,r,loc,x,0.01,0.01+L
cm,CM_3,node
nsel,s,loc,z,0.01+H1
nsel,r,loc,x,0.01,0.01+L
cm,Cm_4,node
cmsel,s,cm_3,node
cmsel,a,cm_4,node
cp,6,volt,all

fini

/solu
antype,static
nsel,s,loc,x,0
d,all,all,0
nsel,s,loc,z,H1        
nsel,r,loc,x,0,0.01  
d,all,uz,0
nsel,s,loc,z,0.01+H1        
nsel,r,loc,x,0,0.01  
d,all,uz,0
nsel,s,loc,x,L1
nsel,r,loc,y,W/2
nsel,r,loc,z,0.02
f,all,fz,-1
nsel,s,loc,z,H1
nsel,r,loc,x,0.01,0.01+L       !设置0电势面
d,all,volt,0
nsel,s,loc,z,0.01+H1
nsel,r,loc,x,0.01,0.01+L       !设置0电势面
d,all,volt,0
nsel,all
solve
fini

请问我错在哪了？谢谢！