PZT导纳曲线的问题

lovetoend

按照论坛中之前有关做PZT阻抗曲线的帖子，我试着学习了。但是画出的导纳曲线没有实部，自己找了好久也不能找的错误在哪里，所以希望大家帮忙解答，谢谢大家
/prep7
ET,1,SOLID5,3　　　　!定义三维固体耦合场单元
MP,DENS,1,7500 !压电陶瓷密度

MP,PERX,1,762 ! X和Y轴向的介电常数
MP,PERZ,1,663 ! Z轴向的介电常数

TB,PIEZ,1        ! 定义压电系数                  
TBDATA,16,12.7179                     
TBDATA,14,12.7179                     
TBDATA,3,-5.2028                     
TBDATA,6,-5.2028                     
TBDATA,9,15.0804
                     
TB,ANEL,1     ! 压电陶瓷的弹性系数                    
TBDATA,1,1.39E11,7.7837E10,7.4284E10     
TBDATA,7,1.39E11,7.4284E10
TBDATA,12,1.541E11
TBDATA,16,3.0581E10
TBDATA,19,2.5641E10
TBDATA,21,2.5641E10

!!!!!建立模型划分网格!!!!!!!
block,0,0.01,0,0.01,0,0.0005
vatt,1,1
lsel,s,loc,z,0
lsel,a,loc,z,0.005
lesize,all,,,10
ALLSEL,ALL  
lsel,s,line,,9,12
lesize,all,,,2
ALLSEL,ALL
vmesh,all

!!!!!!电势面耦合!!!!!!!!!
nsel,s,loc,z,0.0005
cp,1,volt,all   ! 将施加相同电位势的面耦合
*GET,vtop,NODE,,NUM,MIN
allsel,all
nsel,s,loc,z,0
cp,2,volt,all
*GET,vbot,NODE,,NUM,MIN
allsel,all

d,vtop,volt,1   ! 在表面加电压1 V
d,vbot,volt,0   ! 在表面加电压0 V
nsel,all
!!!!!!!!!!谐分析!!!!!!!!!!!!!
/solu
eqslv,sparce
antype,harmic
hropt,full
outers,all,all
harfrq,100000,500000
nsubst,80
kbc,1
solve
FINISH

/POST26
RFORCE,2,vtop,AMPS   
CFACT,,,,1
PROD,3,2,1,,MHOS,,,-2*3.1416   ! 计算导纳Y=-jwQ/V
realvar,4,3,,,YR
imagin,5,3,,,YI
PLVAR,4,5
---------------------------------

hy1966

谐响应分析需要加阻尼，例如：
DMPRAT,0.006                     ! 对所有频率段加固定的阻尼值

lovetoend

hy1966 发表于 2012-7-17 17:33
谐响应分析需要加阻尼，例如：
DMPRAT,0.006                     ! 对所有频率段加固定的阻尼值 ...

真的是这样，灰常感谢啊！
我原来以为阻尼可加可不加的，看来影响还真大