夹心式换能器的谐响应分析

csys,0
nsel,s,loc,x,-fl
nsel,a,loc,x,-fl-ph-ph
*get,nne1,node,,count
cm,electrode1,node
cp,1,volt,all
*get,nelectrode1,node,0,num,min
nsel,s,loc,x,-fl-ph
*get,nne2,node,,count
cm,electrode2,node
cp,2,volt,all
*get,nelectrode2,node,0,num,min
nsel,all
d,nelectrode1,volt,0
d,nelectrode2,volt,0
finish
/solu
antype,modal
modopt,lanb,15,10000,70000
mxpand,15
solve
finish
/prep7
d,electrode1,volt,0,300
d,electrode2,volt,0
/solu
eqslv,sparse
ANTYPE,HARMIC
hropt,full
outres,all,all
harfrq,30800,31000
nsubst,10
kbc,1
solve
finish
谐响应分析最后得出的位移与频率关系的曲线一直不对，求大神看看那里不对。如果需要全部的代码可以联系我QQ：254702127