压电陶瓷测量原理.docx

1、压电陶瓷测量原理压电陶瓷及其测量原理近年来，压电陶瓷的研究发展迅速， 取得一系列重大成果，应用范围不断扩大，已深入到国民经济和尖端技术的各个方面中， 成为不可或缺的现代化工业材料之一。 由于压电材料的各向异性，每一项性能参数在不同的方向所表现出的数值不同， 这就使得压电陶瓷材料的性能参数比一般各向同性的介质材料多得多。 同时，压电陶瓷的众多的性能参数也是它广泛应用的重要基础。（一）压电陶瓷的主要性能及参数（1）压电效应与压电陶瓷在没有对称中心的晶体上施加压力、 张力或切向力时，则发生与应力成比例的介质极化，同时在晶体两端将出现正负电荷，这一现象称为正压电效应；反之，在晶体上施加电场时， 则将产

2、生与电场强度成比例的变形或机械应力， 这一现象称为逆压电效应。 这两种正、逆压电效应统称为压电效应。 晶体是否出现压电效应由构成晶体的原子和离子的排列方式， 即晶体的对称性所决定。 在声波测井仪器中， 发射探头利用的是正压电效应， 接收探头利用的是逆压电效应。（2）压电陶瓷的主要参数1、介质损耗介质损耗是包括压电陶瓷在内的任何电介质的重要品质指标之一。在交变电场下，电介质所积蓄的电荷有两种分量：一种是有功部分（同相） ，由电导过程所引起；另一种为无功部分（异相），由介质弛豫过程所引起。介质损耗是异相分量与同相分量的比值，如图 1所示，IC为同相分量，h为异相分量，IC与总电流I的夹角为& ,其

3、正切值为图1交流电路中电压-电流矢量图（有损耗时）2、机械品质因数机械品质因数是描述压电陶瓷在机械振动时， 材料内部能量消耗程度的一个参数， 它也是衡量压电陶瓷材料性能的一个重要参数。 机械品质因数越大， 能量的损耗越小。产生能量损耗的原因在于材料的内部摩擦。机械品质因数 Qm的定义为:C _ 谐振时振子储存的机械能Qm =谐振时振子每周所损失的机械能父2兀机械品质因数可根据等效电路计算而得式中R1为等效电阻（Q）, 0s为串联谐振角频率（Hz）, Ci为振子谐振时的等效电容（F）, L1为振子谐振时的等效电感。 Qm与其它参数之间的关系将在后续详细推导。不同的压电器件对压电陶瓷材料的 Qm值

4、的要求不同，在大多数的场合下（包括声波测井的压电陶瓷探头），压电陶瓷器件要求压电陶瓷的 Qm值要高。3、压电常数压电陶瓷具有压电性， 即在其外部施加应力时能产生额外的电荷。 其产生的电荷与施加 的应力成比例，对于压力和张力来说，其符号是相反的，电位移 D （单位面积的电荷）和应力仃 的关系表达式为： d =Q =drA式中Q为产生的电荷（C）, A为电极的面积（m2）, d为压电应变常数（C/N）。在 逆压电效应中，施加电场 E时将成比例地产生应变 S,所产生的应变 S是膨胀还是收缩，取决于样品的极化方向。S=dE i, r D S两式中的压电应变常数 d在数值上是相同的，即 d= 二 E另一

5、个常用的压电常数是压电电压常数 g,它表示应力与所产生的电场的关系， 或应变与所引起的电位移的关系。常数 g与d之间有如下关系： g = d式中名为介电系数。在声波测井仪器中， 压电换能器希望具有较高的压电应变常数和压电电压常数，以便能发射较大能量的声波并且具有较高的接受灵敏度。4、机电耦合系数当用机械能加压或者充电的方法把能量加到压电材料上时，由于压电效应和逆压电效应，机械能（或电能）中的一部分要转换成电能（或机械能） 。这种转换的强弱用机电耦合系数k来表示，它是一个量纲为一的量。机电耦合系数是综合反映压电材料性能的参数，它表示压电材料的机械能和电能的耦合效应。机电耦合系数的定义为：k2=电

6、能转变为机械能 或者k2=机械能转变为电能输入电能 输入机械能机电耦合系数不但与材料参数有关， 还与具体压电材料的工作方式有关。 对于压电陶瓷来说，它的大小还与极化程度相关。它只是反映机、电两类能量通过压电效应耦合的强弱，并不代表两类能量之间的转换效率。 压电材料的耦合系数在不同的场合有不同的要求， 当制作换能器时，希望机电耦合系数越大越好。（二）压电换能器的等效电路压电换能器的等效电路表示法， 是利用电学网络术语表示压电陶瓷的机械振动特性， 即把某些力学量模拟为电学量的方法。把压电换能器用等效电路来表示，有很多优点：其一，可以把力学上复杂的振动问题有效地进行简化； 其二，为了得到换能器的各个

7、参数，从而定量地分析或筛选换能器； 其三，实际应用的需要，因为在实际的应用当中， 压电换能器也是接入到具体的电子线路中的，得到压电换能器的等效电路能够更好地对其外围电路进行匹配 设计。由此可见，得到压电换能器的等效电路是十分必要的。2.3压电换能器的谐振特性将压电换能器按照图 2-2所示线路连接。当改变信号频率时，可以发现，通过压电陶瓷换能器的电流也随着发生变化，其变化规律如图 2-3 （a）所示。从图2-3 （a）可以看出,当信号为某一频率 fm时，通过压电陶瓷换能器的电流出现最大值 Imax ；而当信号变到另 频率fn时，传输电流出现最小值 Imin。由流经它的电流随频率的变化可以看出，压

8、电陶瓷2-3 (b)所示。图2-2 压电陶瓷换能器谐振特性接线示意图换能器的阻抗是随频率的变化而变化的，其变化规律同电流相反，如图图2-3 压电陶瓷换能器电流、阻抗同频率的关系曲线从图中可以看出，当信号频率为 fm时，通过压电陶瓷换能器的电流最大，即其等效阻抗最小，导纳最大；当信号频率为 fn时，通过压电陶瓷换能器的电流最小，即其等效阻抗最大，导纳最小。因此把 fm称为最大导纳频率或最小阻抗频率；而把 fn称为最小导纳频率或最大阻抗频率。而当信号频率继续增大时，还会出现一系列的电流的极大值和极小值， 如 图2-4所示。fin fh fml fnl % % f图2-4 压电陶瓷换能器电流随频率变

9、化示意图（多谐振模式）2.2.4压电换能器的等效电路根据交流电路相关知识，对于图 2-5所示好的LC电路来说，其阻抗 Z也随着频率的变化而变化。在图 2-2所示的线路中，用 LC电路代替压电陶瓷换能器，可以发现，在压电陶瓷换能器的谐振频率处，只要选择合适的 Li、Ci、Ri和C0，通过LC电路的电流和LC电路的阻抗的绝对值随频率的变化曲线，分别同图 2-1中的（b）和（c）的关系曲线非常相似。也就是说，在串联谐振频率附近， 压电陶瓷换能器的阻抗特性和谐振特性同LC电路的阻抗特性和频率特性非常相似。 因此，利用机电类比的方法， 可以用一个 LC电路来表示压电陶瓷换能器白参数和特性，这个 LC电路

10、即为压电陶瓷换能器的等效电路。图2-5 LC电路对压电陶瓷换能器来说， 在任何串联谐振频率附近， 其电行为可以用图 2-3所示的LC电路来表示。在压电陶瓷换能器的串联谐振频率附近， 如果值存在一种振动模式， 即没有其它寄生响应，则在串联谐振频率附近很窄的频率范围内， 可以认为压电陶瓷换能器的等效参数R、Ci、Ri和Co与频率无关。在实际中通过选择合适的尺寸进行加工处理，是可以将所需要的振动模式同其他模式充分隔离开来的。另外，考虑到在实际中，在通电之后，压电陶瓷换能器必然会存在能量的损耗，这一能量损耗可用一个并联电阻 R0来等效。所以其最终等效电路图如图 2-6所示。Ci4C11R-iCD图2-

11、6压电陶瓷换能器等效电路图图中串联支路中的 L称为压电陶瓷换能器的动态电感， g称为动态电容，R称为动态电阻。这三个参数用来表征压电陶瓷换能器在工作 (加电源激励产生振动) 的情况下，振动联电阻 R又称静态电阻，表征换能器的电损耗的大小。2.2.5压电换能器的导纳特性根据已得到的压电换能器的等效电路图， 来进一步分析其导纳特性。 为了简化推导,2-5假定压电陶瓷换能器没有电损耗， 即R0=0,此时其等效电路即为一个 LC电路，如图式中：Y为换能器的总的导纳值,(2-1)Y) = jBo =孜0为并联支路的导纳值，Y =G = jB1为串联支路的导纳值。先对串联支路进行分析。Y0 = jB0 =

12、 j 02 G 2 2 G 2 2 Gx = , B1 =（ : x -R1 ） = （x - RI ） - -G1 Gi Ri Ri R2 2G c /12 1 2 2 1 2所以，B1 *g m=两边同时加上（大，可得（G1一与丁） +B1 =H（2-3） R1 2 R1 2 R1 2 R11 1若以电导为横坐标， 电纳为纵坐标，则式（2-3）表示一个以（不,0）为圆心，xg为2 R1 2 R1半径的圆，也即是我们所说的导纳圆。如图 2-7中虚线所示图2-7 导纳圆图对于串联支路进行分析，根据串联谐振频率的定义，令 B1 =0,则由式（2-3）可得到G1 =0或 G =1/R。由于实际的压

13、电陶瓷换能器的动态电阻 R0不可能为零，根据式（2-2）1中G1的表达式可以知道，只有 G1 =1/R满足串联谐振的条件。即： 切匕工=0，所以 , C11可以得到串联支路的谐振频率（又称机械共振频率） ：/（2-4）L1C1接着考虑加入静态电容后的情况。由式（ 2-1）可知，考虑静态电容后换能器的导纳相当于在串联支路的电纳（虚部）加上 丫。鉴于一般情况下，压电陶瓷换能器的机械品质因 数都较大，也即在串联谐振频率 6s附近，丫0 = j0C0的值随频率的变化很小，可以近似认为是一个常数。因此，只需将串联支路所得到的导纳圆的纵坐标向上平移一个常数， 而横坐标保持不变即可得到加入静态电容后换能器的

14、导纳关系图，如图 2-7中点划线所示。若再考虑到换能器的静态电阻并不为零， 则实际中的导纳圆不可能与纵轴相切， 而是向横轴的正,如图2-7中实线圆所示向平移一定的量（平移距离的大小取决于静态电阻的阻值）对导纳圆图进行简要的分析可知：当 f fs即0 0S时，电纳值大于零，当 f A fs即3A8s时，电纳值小于零。所以，随着频率的增加，导纳圆是沿顺时针方向变化的。另外，在串联谐振频率的附近， 还存在着两个频率点使得换能器总的电纳为零， 此时电源信号经过换能器之后只有幅值的改变，而没有相位的变化，也即电压和电流信号同相位。这两个频率 中，值较小的那个频率 fr称为谐振频率，较大 fa的称为反谐振

15、频率。另外还存在使得换能器的导纳值取得最大的频率 fm ,导纳值最小的频率fn。连接原点和串联谐振频率点， 与导纳圆的交点处的频率 fp称为并联谐振频率。另外，需要特别指出的是，上述讨论是在一个振动模态谐振频率 8s附近较小的频率变化范围内进行的，并且只有在导纳圆的直径远大于这个频率范围内 6c0的变化时才是正确的，否则换能器的导纳曲线将变得十分复杂，具有蔓叶曲线的特征。根据以上导纳圆图的推导过程， 下面介绍一下压电陶瓷换能器等效电路中各个参数和导纳圆图的关系，并给出各自的计算公式。在换能器的导纳圆图中作平行于纵轴的直径， 交导纳圆于两点，分别记作f1、f2。在f1点处，串联支路的动态电导和电

16、纳值相等，即 Gi = Bi。由式（2-2）可得:/I 1 、-( iLi -) 11c1Ri2 (/Li -7)2,1C1在f2点处，串联支路的动态电导和电纳值相等， 但符号相反，即Gi = -Bi。由式（2-2）结合式（2-5）和式（2-6）,可得:再由式（2-4）可得:(2-8)C 一二5 一 2 ,s动态电阻的值可以通过导纳圆的直径求得:(2-11)Ro =1 D式中V。为圆心的纵坐标。静态电阻Ro的值可由导纳圆偏离纵轴的距离（或圆心的横坐标）来确定:(2-13)Ro=TX。一 一2R式中X。为圆心的横坐标。至此，我们已得到压电陶瓷换能器等效电路中所有参数的计算公式。2.3测量原理在上

17、一节中，得到的压电陶瓷换能器等效电路参数的计算公式都是基于导纳圆的， 也即是基于各个频率下的电导和电纳值的， 因此我们需要得到每个频率点的导纳值。 为此采用图 2-8所示的测量原理图进行测量。图2-8压电陶瓷换能器测量原理示意图图2-8中，AC为频率可控的交流信号源， R表示源内阻，Rm称为精密电阻，Ua为U B相对于U a会有一个加在压电陶瓷换能器山的电压信号 Ub为经过换能器之后的电压信号。根据前面章节所介绍的压电陶瓷的导纳特性可以知道，在经过换能器之后的电压信号 幅度和相位的变化。不失一般性，在这里设定:UA = UAmej/ UB=UBmej或城 (2-14)其中：UAm, UBm分别

18、表示两路信号的幅值， 0为信号的角频率，e为信号的初始相位，日为两路信号的相位差。按照习惯表达，先求压电陶瓷换能器的阻抗，再取倒数得到导纳。7 U Ua -Ub一 I 一 UbRm将式（2-14）代入得= Rm(Umcos -1) - j Umsin二 R jX UBm UBm对应得到：R=Rm(UmcOS0 -1) , X =-RmUAm-sin9 (2-15)U Bm U Bm再由导纳和阻抗的关系可得11 -Y G jBZ R jX由以上推导可以看出，换能器的导纳和阻抗值仅与加在其两端的电压信号的幅值比和相位差有关，因此只需要得到两路信号的幅值和相位信息即可得到换能器等效电路的各个参数。而

19、实际中，只需要对两路信号进行采样， 再通过对采样所得数据进行处理便可得到幅值和相位信息。2.4正弦信号的测量方法根据上一节介绍的测量原理可知， 要得到压电换能器在测试频率下的电导和电纳值， 就需要测得其两端正弦信号的幅值比和相位差。 但是实际中，硬件电路实现的仅是对两路信号的A/D转换采集，也即是得到的是两路正弦信号的一系列的离散的点。在这一节中，将介 绍从这些采集到的离散的点计算其幅值和相位的方法。2.4.1数字相关法随着微处理器和大规模集成电路的迅速发展， 在测试系统中，越来越多的传统的测量方法被数字化测量方法所取代。近年来，由于相关函数法具有提高测试精度， 减少或简化硬件设计，能够充分利

20、用测试系统中的数据采集系统和微型计算机， 提高测试系统的可靠性和可维护性的诸多优点，使得相关技术原理在相位差的测量及数字信号处理中得到了广泛应用，并展现出良好的应用前景1、相关函数法原理相关函数法利用两同频正弦信号的延时为零时的互相关函数值与其相位差的余弦值成正比的原理获得相位差。设两路被测信号为：Nx(t)、 Ny(t)分别表x(t)和y(t)是相关的,x(t) = Asin(2兀ft) + Nx(t), y(t) = Bsing +甲)+ Ny(t) (2-17)其中：A、B分别表示两路信号的幅值， f表示信号的频率,示两路信号的干扰噪声信号， 十表示两路信号的相位差。显然，信号 则两路信

21、号的互相关函数为:1 丁Rxy( ) 0 x(t)y(t )dtT (2-18)1 TRxy( ) = 0Asin(2二ft) Nx(t) Bsin(2二f(t ) ) Ny(t )dt， ，、， 1式中T为信号的周期，即T =- f当七二0时，有1 TRxy(Q)=- 10 Asin(2nft) +Nx(t) Bsin(2nft +中)+ Ny(t)dt由于噪声信号之间不相关，噪声和信号之间也不相关，将上式进一步展开得:1 TRxy = - o ABsin(2二ft)sin(2二ft )dt1 T=丁 。 ABsin(2二：ft)sin(2：ft)cos+cos(2：ft)sin dt1 T

22、 =0 ABsin(2:ft)sin(2:ft)cos dt =所以，可以得到相位差的计算公式：1 _AB cos :22Rxy(Q) =arccos( )AB(2-19)而信号幅值的大小可由信号的自相关函数求得1 ， 1 -Rx( ) 2rx(t)x(t )dt 2T Asin(2 ：ft)Asin(2 ：f (t )dt (2-2Q)T W T -2当。二Q时，有1 c c cRx(Q) 2r A2sin2(2 ft)dt =T wA22所以可得信号幅值的计算公式：A i J2Rx(Q) (2-21)B = .2Ry(Q)将上式代入式(2-19),可得相位差计算公式的另一种表达式:邛 =a

23、rccos( . Rxy(。) =) Rx(Q)Ry(Q)(2-22)而在实际中，是没有完整精确的信号的表达式的， 有的是对信号的模数转换所得到的离散的数据，离散序列的自相关和互相关的计算公式如下:1 n Rxy(0)- x(i)y(i) n t 1 n1 2Rx(0) = v x(i) (2-23)n i包 1 n 2Ry(0) J y(i)2 n -式中：n表示采样个数，i表示第i个采样点，x(i)、y(i)分别表示两路信号的第 i个点的转换得到数值。由式(2-23)分别求出两路信号的自相关和互相关函数值之后， 再由式(2-21)和式(2-22)即可得到两路信号各自的幅值和它们之间的相位差

24、。 但是，需要指出的是，由数字相关法求得的相位差，并不能区分是超前还是滞后， 这就需要采用其他方法来确定相位差符号的正负号。根据前面测量原理中的介绍，由式( 2-15)和式(2-16)可知，压电陶瓷换能器的电导值仅取决于两路信号相位差的余弦值， 而电纳的值是在电导值取得最大的时候发生变号。 由此，可以先求得电导的值， 再通过循环找其最大值， 并从使电导取得最大值时的相位差开始，把相位差变号，得到新的相位差序列，再由新的相位差序列求电纳的值即可。图 2-9表示的为采用数字相关法对一号压电换能器测量数据的处理结果， 其中(a)表示的是导纳圆图， (b)表示的是电导和电纳值随测试频率的变化曲线。(a

25、) (b)图2-9数字相关法处理结果2、相关函数法的特点及误差分析通过上面对相关函数法测量原理的理论推导过程可以看出， 相关函数法测量信号的幅值和相位差与信号的频率无关。也即是说相关函数法不受频率的影响， 可以用来测量未知频率的信号的相位差。同时，相关函数法测量原理的推导都是基于正弦函数的， 因此，它只能用于测量正弦或余弦信号，并不能测量一般的周期信号。由于噪声干扰信号和原信号并不相关， 所以相关函数法能够有效的抑制噪声干扰。 但是,如果在系统中存在相关性较强的干扰信号， 并且信噪比又比较低的情况下，相关函数法测量误差就会比较大。由相关函数法离散序列的最终计算公式可以看出， 其计算结果与采样的

26、点数有关，也即是说测量误差的大小与采样点数是相关的， 采样点数越大，计算结果越接近真实值，测量误差也就越小。综合以上对相关函数法的特点的分析， 可知相关函数法对于采样转换信号中的直流偏移和噪声等干扰具有很强的抑制能力， 它的误差主要是因为采用有限长度的样本代替了高斯白噪声和均匀分布的 A/D量化误差，使得被检正弦信号与噪声信号并非完全不相关。所以， 相关函数法的测量误差与 A/D转换的位数、信号的信噪比和采集点数有关。2.4.2快速离散傅里叶变换法现代信号分析采用数字化方式实现， 其核心是离散傅立叶变换，它完成了从时域到频域的转换，不仅可以实现线性谱分析，而且还是均方谱分析的关键。离散傅立叶变

27、换( DFT)实现了信号首次在频域表示的离散化， 使得频域也能够用计算机进行处理， 但由于用于实际时计算量太大而使应用受到限制。 直到1965年由Cooly和Tukey建立了一种快速傅立叶变换一一FFT时，DFT的应用才成为现实1、FFT获取正弦波幅值和相位的原理设采集正弦信号得到的离散序列为 x(n), n=1,2,KN。则该序列的离散傅里叶变换为：N 口洛X(k) =DFTx(n) x(n)e Nn -0N 2-：k 2-：k=x(n)cos( n) - j sin( n) = ReX(k) Im X(k) n=0 N Nk =0,1,2,N -1 (2-24)则其初始相位为：r + /R

28、eX(k)、| f da - arctan( ), k 二 - 1ImX(k) fs其中：fs是信号的采样频率，N是采样长度。在对时域离散序列进行傅立叶变换之后， 可以得到其离散的幅度谱和相位谱， 在幅度谱和相位谱中找到对应时域波形的频率的谱线就可以得到时域的正弦波形的幅值和相位信息。图2-10所示的是采用快速离散傅里叶变换法对采集到的数据处理的结果。(a) (b)图2-10快速离散傅里叶变换法处理结果2、FFT的特点及误差分析通过傅里叶变换可以只提取基波参数， 因此谐波的存在并不影响基波成分， 所以谐波的存在对应用这种方法测量相位差几乎没有影响； 对于噪声干扰，只有当高斯白噪声接近基波的频率

29、分量时才会影响到基波的相位， 所以应用FFT法测量相位差也能有效地抑制高斯白噪声干扰。但是，实际上信号是连续的无限长的序列，用 FFT对其进行谱分析时，必须截 短形成有限长序列，再进行周期延拓，这样就不可避免的造成信号频谱的泄漏， 由此便产生了相位差测量误差。误差现象主要是：混叠现象、栅栏效应和截断效应。要想减小相位差测 量误差，就必须提高谱分辨率。实际中可通过提高采样频率或者增加采样数据长度来提高谱 分辨率，进而达到减小相位差测量误差的目的。2.4.3正弦曲线参数拟合法设被测的正弦信号为：f (t) =A0Sin(2 ：ft ) D (2-25)其中：f表示信号频率，Ao表示被测信号幅值，中

30、表示被测信号的初始相位角， D表示被测信号的直流分量。 由于被测信号的频率为已知的， 故只需对测得的数据进行三参数的正弦曲线拟合，即可得到被测信号的幅值和相位信息。为此，进一步将上式展开可得：f = Asin(2ft) Bcos(2 二 ft) D (2-26)其中：A = A0 cos ,B = Aq sin ：从而将被测信号的幅值和初始相角转化为对参数 A、B的求取。其基本思想就是寻找合适的 A、B和D的值，使得其测量残差的平方和取得最小。 ,、 _ ， 1 , 、设每个频率下测量的时间序列为 ti =,（i =1,2,3n）,fsn为测量数据的个数， fs为采样频率，fi为每个点的测量值。则测量残差的表达式为:n一 一 一 _ 2=， f (t) 一 Asin(2 二fti) 一 B cos(2 二fti) 一 D (2-27)i 4要使得上式取得最小值，可对其参数求偏导，并令其为零。即:cz dA cAd- = 田 田n n/ f (i) - Asin(2Tifti) - Bcos(27if