压电马达 频率好低，怎么回事

916754385

命令流见下：
/title,3wei
/units,mks
/prep7

!PZT材料特性参数
mp,dens,1,7500    !密度 kg/m^3
mp,perx,1,804.6
mp,pery,1,804.6
mp,perz,1,659.7

!压电常数  C/m^2
tb,piez,1
tbdata,3,-4.1
tbdata,6,-4.1
tbdata,9,14.1
tbdata,14,10.5
tbdata,16,10.5

!弹性柔顺系数 N/m^2
tb,anel,1
tbdata,1,13.2e10,7.1e10,7.3e10
tbdata,7,13.2e10,7.3e10
tbdata,12,11.5e10
tbdata,16,3.0e10
tbdata,19,2.6e10
tbdata,21,2.6e10

!定义主结构的材料参数
mp,dens,2,7800
mp,ex,2,209e9
mp,nuxy,2,0.3  

!定义压电复合模型的几何参数
a=9e-3
b=2.5e-3
c=6.5e-3
d=7e-3
e=14e-3
f=2.5e-3
h=10e-3
thick=6e-3
wide=8e-3

!定义主结构模型
k,1,a,b,0
k,2,f,c,0
k,3,-f,c,0
k,4,-a,b,0
k,5,-a,-b,0
k,6,-f,-c,0
k,7,f,-c,0
k,8,a,-b,0
k,9,e,0,0
k,10,0,d,0
k,11,-e,0,0
k,12,0,-d,0

a,9,10,11,12
lfillt,1,2,2e-3
lfillt,2,3,2e-3
lfillt,3,4,2e-3
lfillt,1,4,2e-3

ASBL,       1,       5  
ASBL,       3,       8  
ASBL,       4,      11  
ASBL,       5,      14  

adele,1,4

numcmp,area

a,1,2,3,4,5,6,7,8

asba,2,1

numcmp,area

k,100,0,0,0
k,101,0,0,h
l,100,101

VDRAG,       1, , , , , ,      25  !NV1
!建立压电陶瓷模型

block,-a-thick,-a,-b,b,(h-wide)/2,(h+wide)/2        !NV2
block,a,a+thick,-b,b,(h-wide)/2,(h+wide)/2           !NV3

vglue,all

numcmp,volume

!赋予材料属性
et,1,solid5,3
et,2,solid186

!采用映射划分网格连接相邻面
!asel,s,loc,x,-a
!cm,cm_1,area
!cmplot,cm_1
!accat,cm_1

!asel,s,loc,x,a
!cm,cm_2,area
!accat,cm_2

!进行网格划分
!LOCAL, KCN, KCS, XC, YC, ZC, THXY, THYZ, THZX, PAR1, PAR2
local,11,0,-a-thick,0,h/2,0,0,90
local,12,0,a+thick,0,h/2,0,0,270
!lesize,all,l_esize,,,,1,,,1,
mat,1
type,1
esys,11
vmesh,2

mat,1
type,1
esys,12
vmesh,3

mat,2
type,2
esys,0
MSHKEY,0
mshape,1
vmesh,1

csys,0
!定义约束
!D, NODE, Lab, VALUE, VALUE2, NEND, NINC, Lab2, Lab3, Lab4, Lab5, Lab6
!*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM
!加载边界条件
nsel,s,loc,x,-a-thick
cp,1,volt,all
*get,n1,node, ,num,min

nsel,s,loc,x,a+thick
cp,2,volt,all
*get,n2,node, ,num,min

nsel,s,loc,x,-a
cp,3,volt,all
*get,n3,node, ,num,min

nsel,s,loc,x,a
cp,4,volt,all
*get,n4,node, ,num,min

nsel,all

!定义电压边界
d,n3,volt,0
d,n4,volt,0

finish


/solu

antype,trans
trnopt,full
tintp,,0.25,0.5,0.5
outpr,all,all
outres,all,all
kbc,1

!加载边界条件
nsel,s,loc,x,-a-thick
d,all,ux,0, , , , uy,uz

nsel,s,loc,x,a+thick
d,all,ux,0, , , , uy,uz

nsel,all

!定义电压边界
d,n3,volt,0
d,n4,volt,0


pi=3.1415926
freq=50
omega=freq*2*pi
cishu=20
*do,i,1,cishu+1,1
t=i/(cishu*omega)
time,t
nsubst,1
v1=10*sin(2*3.14159*omega*t)
v2=10*sin(2*3.14159*omega*t+pi/2)

d,n1,volt,v1             ! apply v Volt
d,n2,volt,v2             ! apply v Volt

allsel,all
solve
*enddo

finish

/post26
num=462
NSOL,2,num,U,x, Ux_1,
NSOL,3,num,U,y, Uy_1,
xvar,2
plvar,3

916754385

这个压电马达定子是按照一个硕士论文做的，
《南京航空航天大学》 2008年
压电型直线电机的设计
李艳林  

916754385

不知道为什么做不出她的效果，我这个频率很低啊，如果频率放高 了，接触点的运动轨迹就变得杂乱无章，只有放的很低才能看到较好的椭圆运动轨迹