ansys 压电单晶阻抗分析

例子是从一篇论文中看到的，基本情况如下：模拟一个长度振动模式的柱体的阻抗谱
长度振动模式的谐振器：4.8mm×4.8mm×30mm
基本参数：(单位MKS制)
弹性柔顺系数s11=83.5e-12，s12=-19.9e-12，s13=-31.1e-12，s33=70.2e-12，s44=172.9e-12，s66=237.2e-12
压电常数：d31=-469e-12，d33=1094e-12，d15=4.9e-12
介电常数：ept11=16.4，ept33=2581
密度：dens=5336
单元类型：SOLID5（自由度X，Y，Z，Volt）
边界条件：长度方向（30mm）上表面1V耦合电压，下表面0V耦合电压，无其他边界条件
谐响应分析Harmonic
damping=0.007
频率范围：1kHz~80kHz
问题：自己按上述条件做出来之后，谐振频率倒是对的（约为29kHz）,但是反谐振频率与论文中的有很大出入（我模拟出来的是36.6kHz,论文中是54kHz）。
不知道到底哪里出错了，跪求做过类似模拟的大神指导！谢谢！