纵扭超声电机仿真求助

weibaobao

对纵扭超声电机进行分析，模态分析已经OK，但是在谐波分析时遇到问题，最关键的问题是，扭转PZT的极化问题。这是一个四分极的圆形扭转PZT，两个正极两个负极，我现在已经在两个正极上加了1V的电压，负极加了-1V，但是仿真结果怎么都没有位移，一点儿都没有，我一直找不出问题在哪里，请各位帮我看一看问题出在哪儿吧。代码如下（定义材料的代码是引用的，应该没有问题，应该之前用这个材料代码分析都对的，所以这里没有列出来。另外总共是两片扭振子，但我现在只给一片加了载荷，是想等把程序调试好了再给另一片加的）：


mm = 1e-3
kHz = 1e3
pi = 3.1415926535

TPI = 12.0/2*mm
TPO = 25.0/2*mm
TPT = 1.00*mm
NPI = 10.0/2*mm
NPO = 25.0/2*mm
NPT = 4.0*mm
TV1 = 9.0*mm
TV2 = 2.0*mm
TV3 = 1.0*mm
TV4 = 1.0*mm
TV5 = 9.0*mm
TV6 = 1.0*mm
TV7 = 5.0*mm
TV8 = 3.0*mm

YVI1 = 10.0/2*mm
YVO1 = 25.0/2*mm
YVO2 = 38.0/2*mm
RYVI1 = 13.0/2*mm
RYVO1 = 25.0/2*mm
Z0 = 0*mm
Z1 = Z0 + TV1
P11 = Z1 + TV2
P12 = P11 + NPT
P13 = P12 + NPT
Z2 = P13 + TV3
Z3 = Z2 + TV4
Z4 = Z3 + TV5
P21 = Z4 + TV6
P22 = P21 + TPT
P23 = P22 + TPT
P24 = P23 + TPT
P25 = P24 + TPT
P26 = P25 + TPT
P27 = P26 + TPT
Z5 = P27 + TV7
Z6 = Z5 + TV8

CSYS, 0
WPROTA,,,-90
/PREP7
/INPUT,'matlprop2','inp','C:\Users\Sheng Li\Desktop\code of ANSYS'

RECTNG, Z0, Z6, 0, YVO2
allsel
WPCSYS,-1
WPROTA,,90

WPAVE,,YVI1 $ASBW,ALL
WPAVE,,YVO1 $ASBW,ALL

CSYS,0      $WPCSYS,-1
WPAVE,,,Z1 $ASBW,ALL
WPAVE,,,P11 $ASBW,ALL
WPAVE,,,P12 $ASBW,ALL
WPAVE,,,P13 $ASBW,ALL
WPAVE,,,Z2 $ASBW,ALL
WPAVE,,,Z3 $ASBW,ALL
WPAVE,,,Z4 $ASBW,ALL
WPAVE,,,P21 $ASBW,ALL
WPAVE,,,P22 $ASBW,ALL
WPAVE,,,P23 $ASBW,ALL
WPAVE,,,P24 $ASBW,ALL
WPAVE,,,P25 $ASBW,ALL
WPAVE,,,P26 $ASBW,ALL
WPAVE,,,P27 $ASBW,ALL
WPAVE,,,Z5 $ASBW,ALL

WPCSYS,-1
ASEL,S,LOC,Z, Z0, Z2
ASEL,R,LOC,Y, YVO1, YVO2
ADELE,ALL,,,1
ASEL,S,LOC,Z, Z3, Z6
ASEL,R,LOC,Y, YVO1, YVO2
ADELE,ALL,,,1
ALLSEL

WPCSYS, -1
K, 10001,,, -1*mm
K, 10002,,, Z6 + 1*mm
VROTAT,ALL,,,,,, 10001, 10002,-360,8
ALLSEL
VGLUE,ALL

et, 1, solid45
et, 2, solid5
keyopt, 2, 1, 3

csys, 1
vsel, s, loc, z, 0, Z1
vsel, r, loc, x, -YVO1,YVO1  
vatt, 3, , 1, 0
vsel, s, loc, z, Z1, P11
vsel, r, loc, x, -YVO1, -YVI1
vsel, r, loc, x, YVI1, YVO1
vatt, 9, , 1, 0
vsel, s, loc, z, P11, P13
vsel, r, loc, x, -YVO1, YVI1
vsel, r, loc, x, YVI1, YVO1
vatt, 102, , 2, 0
vsel, s, loc, z, P13, Z2
vsel, r, loc, x, -YVO1,-YVI1
vsel, r, loc, x, YVI1, YVO1
vatt, 9, , 1, 0
vsel, s, loc, z, Z2, Z3
vsel, r, loc, x, -YVO2,-YVI1
vsel, r, loc, x, YVI1,YVO2  
vatt, 9, , 1, 0
vsel, s, loc, z, Z3, Z4
vsel, r, loc, x, -YVO1,-YVI1
vsel, r, loc, x, YVI1, YVO1
vatt, 9, , 1, 0
vsel, s, loc, z, Z4, P21
vsel, r, loc, x, -YVO1,-YVI1
vsel, r, loc, x, YVI1, YVO1  
vatt, 9, , 1, 0
vsel, s, loc, z, P21, P27
vsel, r, loc, x, -YVO1,-TPI
vsel, r, loc, x, TPI, YVO1  
vatt, 100, , 2, 0
vsel, s, loc, z, P27, Z5
vsel, r, loc, x, -YVO1,-YVI1
vsel, r, loc, x, YVI1,YVO1  
vatt, 9, , 1, 0
vsel, s, loc, z, Z5, Z6
vsel, r, loc, x, -YVO1,YVO1  
vatt, 3, , 1, 0
vsel, s, loc, z, Z1, Z6
vsel, r, loc, x, -YVI1,YVI1
vatt, 3, , 1, 0
vsel, s, loc, z,Z6 , Z6 + R1 + R2
vsel, r, loc, x, -RYVO1,RYVO1  
vatt, 3, , 1, 0
ALLSEL

meshsize = 3*mm
esize, meshsize
vsweep,all
VPLOT

allsel

asel,s,loc,z,P11
nsla,s,1
cp,3,volt,all
*get,gnd,node,0,num,min
allsel

asel,s,loc,z,P12
nsla,s,1
nsel,r,loc,x,0.0005,1
nsel,r,loc,y,0.0005,1
cp,4,volt,all
*get,plus1,node,0,num,min
allsel

asel,s,loc,z,P12
nsla,s,1
nsel,r,loc,x,-1,-0.0005
nsel,r,loc,y,-0.0005,-1
cp,5,volt,all
*get,plus2,node,0,num,min
allsel

asel,s,loc,z,P12
nsla,s,1
nsel,r,loc,x,0.0005,1
nsel,r,loc,y,-0.0005,-1
cp,6,volt,all
*get,plus3,node,0,num,min
allsel

asel,s,loc,z,P12
nsla,s,1
nsel,r,loc,x,-0.0005,-1
nsel,r,loc,y,0.0005,1
cp,7,volt,all
*get,plus4,node,0,num,min
allsel

vsel, s, loc, z, P21, P22
vsel, a, loc, z, P23,P24
vsel, a, loc, z, P25,P26
local,13,CART,0,0,0, , ,180
wpcsys,-1,13
eslv,s,1
emodif,all,esys,13
csys,0
wpcsys,-1
allsel

asel, s, loc, z, P22
asel, a, loc, z, P24
asel, a, loc, z, P26
nsla, s, 1
cp, 1, volt, all
*get, plus, node, 0, num, min
allsel

asel, s, loc, z, P21
asel, a, loc, z, P23
asel, a, loc, z, P25
asel, a, loc, z, P27
nsla, s, 1
cp, 2, volt, all
*get, gnd, node, 0, num, min
allsel

d, plus, volt, 1, 0
d, gnd, volt, 0, 0
d, plus1, volt, 1, 0
d, plus2, volt, 1, 0
d, plus3, volt, -1, 0
d, plus4, volt, -1, 0
allsel

FREQSTART = 15*kHz
FREQEND =  35*kHz
DELTAFREQ = 1*KHz
FREQBUSRTNUM = (FREQEND - FREQSTART) / DELTAFREQ
Qfactor = 100
lossfactor = 1/Qfactor/2
/SOLU        
!*  
ANTYPE,3     
HARFRQ,FREQSTART,FREQEND,   
NSUBST,FREQBUSRTNUM,   
KBC,1        
!*  
ALPHAD,0,      
BETAD,0,        
DMPRAT,lossfactor,   
!*
solve

端木方正

不知道。非常感谢楼主的资料