按照论坛中之前有关做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 ---------------------------------

按照论坛中之前有关做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 ---------------------------------