问一个关于施加耦合场的问题
上图为一个双片的悬臂梁压电振子，我应该如何施加耦合场做分析，才不会使我的计算出错。我的命令流如下，
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
请问我错在哪了？谢谢！
页:
[1]