压电陶瓷测量基本知识.docx
《压电陶瓷测量基本知识.docx》由会员分享，可在线阅读，更多相关《压电陶瓷测量基本知识.docx（28页珍藏版）》请在冰豆网上搜索。
压电陶瓷测量基本知识
压电陶瓷及其测量原理
近年来，压电陶瓷的研究发展迅速，取得一系列重大成果，应用范围不断扩大，已深入到国民经济和尖端技术的各个方面中，成为不可或缺的现代化工业材料之一。
由于压电材料的各向异性，每一项性能参数在不同的方向所表现出的数值不同，这就使得压电陶瓷材料的性能参数比一般各向同性的介质材料多得多。
同时，压电陶瓷的众多的性能参数也是它广泛应用的重要基础。
（一）压电陶瓷的主要性能及参数
（1）压电效应与压电陶瓷
在没有对称中心的晶体上施加压力、张力或切向力时，则发生与应力成比例的介质极化，同时在晶体两端将出现正负电荷，这一现象称为正压电效应；反之，在晶体上施加电场时，则将产生与电场强度成比例的变形或机械应力，这一现象称为逆压电效应。
这两种正、逆压电效应统称为压电效应。
晶体是否出现压电效应由构成晶体的原子和离子的排列方式，即晶体的对称性所决定。
在声波测井仪器中，发射探头利用的是正压电效应，接收探头利用的是逆压电效应。
（2）压电陶瓷的主要参数
1、介质损耗
介质损耗是包括压电陶瓷在内的任何电介质的重要品质指标之一。
在交变电场下，电介质所积蓄的电荷有两种分量：
一种是有功部分（同相），由电导过程所引起；另一种为无功部分（异相），由介质弛豫过程所引起。
介质损耗是异相分量与同相分量的比值，如图1所示，IC为同相分量，Ir为异相分量，IC与总电流I的夹角为，其正切值为
图1交流电路中电压-电流矢量图（有损耗时）
2、机械品质因数
机械品质因数是描述压电陶瓷在机械振动时，材料内部能量消耗程度的一个参数，它也是衡
量压电陶瓷材料性能的一个重要参数。
机械品质因数越大，能量的损耗越小。
产生能量损耗
的原因在于材料的内部摩擦。
机械品质因数Qm的定义为:
谐振时振子储存的机械能c
Qm谐振时振子每周所损失的机械能2兀
机械品质因数可根据等效电路计算而得
sR1C1
sL1
式中R1为等效电阻（Q）,s为串联谐振角频率（Hz）,C1为振子谐振时的等效电容
（F），L1为振子谐振时的等效电感。
Qm与其它参数之间的关系将在后续详细推导。
不同的压电器件对压电陶瓷材料的Qm值的要求不同，在大多数的场合下（包括声波
测井的压电陶瓷探头），压电陶瓷器件要求压电陶瓷的Qm值要高。
3、压电常数
压电陶瓷具有压电性，即在其外部施加应力时能产生额外的电荷。
其产生的电荷与施加
D（单位面积的电荷）和
的应力成比例，对于压力和张力来说，其符号是相反的，电位移
Q
应力的关系表达式为：
Ddr
A
式中Q为产生的电荷（C）,A为电极的面积（m2）,d为压电应变常数（C/N）。
在逆压电效应中，施加电场E时将成比例地产生应变S,所产生的应变S是膨胀还是收缩，
取决于样品的极化方向。
S=dE
DS
两式中的压电应变常数d在数值上是相同的，即d
E
另一个常用的压电常数是压电电压常数g，它表示应力与所产生的电场的关系，或应
变与所引起的电位移的关系。
常数g与d之间有如下关系：
g-
式中为介电系数。
在声波测井仪器中，压电换能器希望具有较高的压电应变常数和压电电压常数，以便能发射较大能量的声波并且具有较高的接受灵敏度。
4、机电耦合系数
当用机械能加压或者充电的方法把能量加到压电材料上时，由于压电效应和逆压电效应，机械能（或电能）中的一部分要转换成电能（或机械能）。
这种转换的强弱用机电耦合
系数k来表示，它是一个量纲为一的量。
机电耦合系数是综合反映压电材料性能的参数，它表示压电材料的机械能和电能的耦合效应。
机电耦合系数的定义为：
k2=电能转变为机械能或者k2=机械能转变为电能
输入电能或输入机械能—
机电耦合系数不但与材料参数有关，还与具体压电材料的工作方式有关。
对于压电陶瓷
来说，它的大小还与极化程度相关。
它只是反映机、电两类能量通过压电效应耦合的强弱，并不代表两类能量之间的转换效率。
压电材料的耦合系数在不同的场合有不同的要求，当制
作换能器时，希望机电耦合系数越大越好。
（二）压电换能器的等效电路
压电换能器的等效电路表示法，是利用电学网络术语表示压电陶瓷的机械振动特性，即
把某些力学量模拟为电学量的方法。
把压电换能器用等效电路来表示，有很多优点：
其一，可以把力学上复杂的振动问题有效地进行简化；其二，为了得到换能器的各个参数，从而定
量地分析或筛选换能器；其三，实际应用的需要，因为在实际的应用当中，压电换能器也是
接入到具体的电子线路中的，得到压电换能器的等效电路能够更好地对其外围电路进行匹配设计。
由此可见，得到压电换能器的等效电路是十分必要的。
2.3压电换能器的谐振特性
将压电换能器按照图2-2所示线路连接。
当改变信号频率时，可以发现，通过压电陶
瓷换能器的电流也随着发生变化，其变化规律如图2-3（a）所示。
从图2-3（a）可以看
出，当信号为某一频率fm时，通过压电陶瓷换能器的电流出现最大值Imax；而当信号变到
另一频率fn时，传输电流出现最小值Imin。
由流经它的电流随频率的变化可以看出，压电
陶瓷换能器的阻抗是随频率的变化而变化的，其变化规律同电流相反，如图2-3（b）所示。
IE电陶瓷
图2-2压电陶瓷换能器谐振特性接线示意图
图2-3压电陶瓷换能器电流、阻抗同频率的关系曲线
（a）电流-频率关系曲线（b）阻抗-频率关系曲线
从图中可以看出，当信号频率为fm时，通过压电陶瓷换能器的电流最大，即其等效阻
抗最小，导纳最大；当信号频率为fn时，通过压电陶瓷换能器的电流最小，即其等效阻抗
最大，导纳最小。
因此把fm称为最大导纳频率或最小阻抗频率；而把fn称为最小导纳频率
或最大阻抗频率。
而当信号频率继续增大时，还会出现一系列的电流的极大值和极小值，如
图2-4所示。
图2-4压电陶瓷换能器电流随频率变化示意图（多谐振模式）
2.2.4压电换能器的等效电路
根据交流电路相关知识，对于图2-5所示好的LC电路来说，其阻抗Z也随着频率
的变化而变化。
在图2-2所示的线路中，用LC电路代替压电陶瓷换能器，可以发现，在
压电陶瓷换能器的谐振频率处，只要选择合适的Li、Ci、R和Co，通过LC电路的电流
和LC电路的阻抗的绝对值随频率的变化曲线，分别同图2-1中的（b）和（c）的关系曲
线非常相似。
也就是说，在串联谐振频率附近，压电陶瓷换能器的阻抗特性和谐振特性同LC
电路的阻抗特性和频率特性非常相似。
因此，禾U用机电类比的方法，可以用一个LC电路
来表示压电陶瓷换能器的参数和特性，这个LC电路即为压电陶瓷换能器的等效电路。
图2-5LC电路
对压电陶瓷换能器来说，在任何串联谐振频率附近，其电行为可以用图2-3所示的LC
电路来表示。
在压电陶瓷换能器的串联谐振频率附近，如果值存在一种振动模式，即没有其
它寄生响应，则在串联谐振频率附近很窄的频率范围内，可以认为压电陶瓷换能器的等效参
数R、Ci、Ri和C0与频率无关。
在实际中通过选择合适的尺寸进行加工处理，是可以将所需要的振动模式同其他模式充分隔离开来的。
另外，考虑到在实际中，在通电之后，压电陶瓷换能器必然会存在能量的损耗，这一
能量损耗可用一个并联电阻Ro来等效。
所以其最终等效电路图如图2-6所示。
0-
Ro
CD
图2-6压电陶瓷换能器等效电路图
图中串联支路中的Li称为压电陶瓷换能器的动态电感，Ci称为动态电容，Ri称为动态
电阻。
这三个参数用来表征压电陶瓷换能器在工作
（加电源激励产生振动）的情况下，振动
部分所受到的力阻抗和介质对振动的反作用的强弱。
并联电容C0又称静态电容，表征压电
陶瓷换能器在未加激励的情况下等效为一个纯电容，它的值的大小与换能器的形状有关。
并
联电阻R0又称静态电阻，表征换能器的电损耗的大小。
225压电换能器的导纳特性
根据已得到的压电换能器的等效电路图，来进一步分析其导纳特性。
为了简化推导,
为串联支路的导纳值。
先对串联支路进行分析。
YojBOj0
为半径的圆，也即是我们所说的导纳圆。
如图2-7中虚线所示
图2-7导纳圆图
对于串联支路进行分析，根据串联谐振频率的定义，令Bi=0，则由式（2-3）可得到
接着考虑加入静态电容后的情况。
由式（2-i）可知，考虑静态电容后换能器的导纳相
当于在串联支路的电纳（虚部）加上丫。
。
鉴于一般情况下，压电陶瓷换能器的机械品质因
数都较大，也即在串联谐振频率s附近，丫0jCo的值随频率的变化很小，可以近似认
为是一个常数。
因此，只需将串联支路所得到的导纳圆的纵坐标向上平移一个常数，而横坐
标保持不变即可得到加入静态电容后换能器的导纳关系图，如图2-7中点划线所示。
若再
考虑到换能器的静态电阻并不为零，则实际中的导纳圆不可能与纵轴相切，而是向横轴的正
向平移一定的量（平移距离的大小取决于静态电阻的阻值），如图2-7中实线圆所示
对导纳圆图进行简要的分析可知：
当ffs即s时，电纳值大于零，当ffs即
s时，电纳值小于零。
所以，随着频率的增加，导纳圆是沿顺时针方向变化的。
另外,
换能器之后只有幅值的改变，而没有相位的变化，也即电压和电流信号同相位。
这两个频率中，值较小的那个频率fr称为谐振频率，较大fa的称为反谐振频率。
另外还存在使得换能器的导纳值取得最大的频率fm，导纳值最小的频率fn。
连接原点和串联谐振频率点，
与导纳圆的交点处的频率fp称为并联谐振频率。
另外，需要特别指出的是，上述讨论是在一个振动模态谐振频率s附近较小的频率变化范围内进行的，并且只有在导纳圆的直径远大于这个频率范围内Co的变化时才是正确的，否则换能器的导纳曲线将变得十分复杂，具有蔓叶曲线的特征。
根据以上导纳圆图的推导过程，下面介绍一下压电陶瓷换能器等效电路中各个参数和导
纳圆图的关系，并给出各自的计算公式。
在f2点处，串联支路的动态电导和电纳值相等，但符号相反，即GBi。
由式（2-2）
Li
结合式（2-5）和式（2-6），可得:
（2-7）
再由式（2-4）可得:
1
（2-8）
Ci2
sLi
c1sL11L1
机械品质因数：
QmSC1R1R1RG（2-9）
动态电阻的值可以通过导纳圆的直径求得:
（2-11）
Ro1D
静态电容Co的值也可由导纳圆偏离横轴的距离来确定：
Coyo（2-12）
s
式中yo为圆心的纵坐标。
Ro
静态电阻Ro的值可由导纳圆偏离纵轴的距离（或圆心的横坐标）来确定:
（2-13）
式中xo为圆心的横坐标。
至此，我们已得到压电陶瓷换能器等效电路中所有参数的计算公式。
2.3测量原理
在上一节中，得到的压电陶瓷换能器等效电路参数的计算公式都是基于导纳圆的，也即
是基于各个频率下的电导和电纳值的，因此我们需要得到每个频率点的导纳值。
为此采用图2-8所示的测量原理图进行测量。
图2-8压电陶瓷换能器测量原理示意图
图2-8中，AC为频率可控的交流信号源，R表示源内阻，Rm称为精密电阻，Ua为
加在压电陶瓷换能器山的电压信号Ub为经过换能器之后的电压信号。
根据前面章节所介
绍的压电陶瓷的导纳特性可以知道，在经过换能器之后的电压信号Ub相对于Ua会有一个
幅度和相位的变化。
不失一般性，在这里设定：
UaUAmej，UbUBmej（）（2-14）
其中：
UAm，UBm分别表示两路信号的幅值，为信号的角频率，为信号的初始
相位，为两路信号的相位差。
按照习惯表达，先求压电陶瓷换能器的阻抗，再取倒数得到导纳。
UUaUb
ZTUb
Rm
将式（2-14）代入得
再由导纳和阻抗的关系可得
由以上推导可以看出,
换能器的导纳和阻抗值仅与加在其两端的电压信号的幅值比和相位差有关，因此只需要得到两路信号的幅值和相位信息即可得到换能器等效电路的各个参
数。
而实际中，只需要对两路信号进行采样，再通过对采样所得数据进行处理便可得到幅值和相位信息。
2.4正弦信号的测量方法
根据上一节介绍的测量原理可知，要得到压电换能器在测试频率下的电导和电纳值，就
需要测得其两端正弦信号的幅值比和相位差。
但是实际中，硬件电路实现的仅是对两路信号
的A/D转换采集，也即是得到的是两路正弦信号的一系列的离散的点。
在这一节中，将介
绍从这些采集到的离散的点计算其幅值和相位的方法。
2.4.1数字相关法
随着微处理器和大规模集成电路的迅速发展，在测试系统中，越来越多的传统的测量方
法被数字化测量方法所取代。
近年来，由于相关函数法具有提高测试精度，减少或简化硬件
设计，能够充分利用测试系统中的数据采集系统和微型计算机，提高测试系统的可靠性和可
维护性的诸多优点，使得相关技术原理在相位差的测量及数字信号处理中得到了广泛应用，并展现出良好的应用前景
1、相关函数法原理
相关函数法利用两同频正弦信号的延时为零时的互相关函数值与其相位差的余弦值成
正比的原理获得相位差。
设两路被测信号为:
x（t）Asin（2ft）Nx（t）,y（t）Bsin（2ft）Ny（t）（2-17）
则两路信号的互相关函数为:
1T
Rxy（）-0x（t）y（t）dt
T（2-18）1T
Rxy（）T0[Asin（2ft）Nx（t）][Bsin（2f（t））Ny（t）]dt
式中T为信号的周期，即T1
当=0时，有
1T
Rxy（O）0[Asin（2ft）Nx（t）][Bsin（2ft）Ny（t）]dt
由于噪声信号之间不相关，噪声和信号之间也不相关，将上式进一步展开得:
所以，可以得到相位差的计算公式:
而信号幅值的大小可由信号的自相关函数求得
当T=0时，有
所以可得信号幅值的计算公式:
A.2&（0）
B,'2Ry（0）（2-21）
将上式代入式（2-19）,可得相位差计算公式的另一种表达式:
arccos（Ry⑼）（2-22）
£Rx（0）Ry（0）
而在实际中，是没有完整精确的信号的表达式的，有的是对信号的模数转换所得到的离散的数据，离散序列的自相关和互相关的计算公式如下:
1n1
Rxy（0）—x（i）y（i）
ni0
1n12
Rx（0）x（i）（2-23）
ni0
1n12
Ry（0）—y（i）
ni0
式中：
n表示米样个数，i表示第i个米样点，x（i）、y（i）分别表示两路信号的第i个点的转换得到数值。
由式（2-23）分别求出两路信号的自相关和互相关函数值之后，再由式（2-21）和式
（2-22）即可得到两路信号各自的幅值和它们之间的相位差。
但是，需要指出的是，由数
字相关法求得的相位差，并不能区分是超前还是滞后，这就需要采用其他方法来确定相位差
符号的正负号。
根据前面测量原理中的介绍，由式（2-15）和式（2-16）可知，压电陶瓷
换能器的电导值仅取决于两路信号相位差的余弦值，而电纳的值是在电导值取得最大的时候发生变号。
由此，可以先求得电导的值，再通过循环找其最大值，并从使电导取得最大值时的相位差开始，把相位差变号，得到新的相位差序列，再由新的相位差序列求电纳的值即可。
图2-9表示的为采用数字相关法对一号压电换能器测量数据的处理结果，其中（a）表示的是导纳圆图，（b）表示的是电导和电纳值随测试频率的变化曲线。
（b）
图2-9数字相关法处理结果
2、相关函数法的特点及误差分析
通过上面对相关函数法测量原理的理论推导过程可以看出，相关函数法测量信号的幅值
和相位差与信号的频率无关。
也即是说相关函数法不受频率的影响，可以用来测量未知频率
的信号的相位差。
同时，相关函数法测量原理的推导都是基于正弦函数的，因此，它只能用
于测量正弦或余弦信号，并不能测量一般的周期信号。
由于噪声干扰信号和原信号并不相关，所以相关函数法能够有效的抑制噪声干扰。
但是,
如果在系统中存在相关性较强的干扰信号，并且信噪比又比较低的情况下，相关函数法测量
误差就会比较大。
由相关函数法离散序列的最终计算公式可以看出，其计算结果与采样的点
数有关，也即是说测量误差的大小与采样点数是相关的，采样点数越大，计算结果越接近真
实值，测量误差也就越小。
综合以上对相关函数法的特点的分析，可知相关函数法对于采样转换信号中的直流偏移
和噪声等干扰具有很强的抑制能力，它的误差主要是因为采用有限长度的样本代替了高斯白
噪声和均匀分布的A/D量化误差，使得被检正弦信号与噪声信号并非完全不相关。
所以，相关函数法的测量误差与A/D转换的位数、信号的信噪比和采集点数有关。
2.4.2快速离散傅里叶变换法
现代信号分析采用数字化方式实现，其核心是离散傅立叶变换，它完成了从时域到频域
的转换，不仅可以实现线性谱分析，而且还是均方谱分析的关键。
离散傅立叶变换（DFT）
实现了信号首次在频域表示的离散化，使得频域也能够用计算机进行处理，但由于用于实际
时计算量太大而使应用受到限制。
直到1965年由Cooly和Tukey建立了一种快速傅立叶变换一一FFT时，DFT的应用才成为现实
1、FFT获取正弦波幅值和相位的原理
设采集正弦信号得到的离散序列为x（n）,n=1,2，…KN。
则该序列的离散傅里叶变换为:
.2kjnx（n）eN
n0
k0,1,2,N-1（2-24）
则其初始相位为：
丄，Re[X（k）]、_f。
‘
arctan（）,k1
Im[X（k）]fs
其中：
fs是信号的采样频率，N是采样长度。
在对时域离散序列进行傅立叶变换之后，可以得到其离散的幅度谱和相位谱，在幅度谱
和相位谱中找到对应时域波形的频率的谱线就可以得到时域的正弦波形的幅值和相位信息。
图2-10所示的是采用快速离散傅里叶变换法对采集到的数据处理的结果。
图2-10快速离散傅里叶变换法处理结果
2、FFT的特点及误差分析
通过傅里叶变换可以只提取基波参数，因此谐波的存在并不影响基波成分，所以谐波的
存在对应用这种方法测量相位差几乎没有影响；对于噪声干扰，只有当高斯白噪声接近基波
的频率分量时才会影响到基波的相位，所以应用FFT法测量相位差也能有效地抑制高斯白
噪声干扰。
但是，实际上信号是连续的无限长的序列，用FFT对其进行谱分析时，必须截
短形成有限长序列，再进行周期延拓，这样就不可避免的造成信号频谱的泄漏，由此便产生
了相位差测量误差。
误差现象主要是：
混叠现象、栅栏效应和截断效应。
要想减小相位差测量误差，就必须提高谱分辨率。
实际中可通过提高采样频率或者增加采样数据长度来提高谱分辨率，进而达到减小相位差测量误差的目的。
243正弦曲线参数拟合法
设被测的正弦信号为：
f（t）Aosin（2ft）D（2-25）
其中：
f表示信号频率，A0表示被测信号幅值，表示被测信号的初始相位角，D表
示被测信号的直流分量。
由于被测信号的频率为已知的，故只需对测得的数据进行三参数的正弦曲线拟合，即可得到被测信号的幅值和相位信息。
为此，进一步将上式展开可得:
（2-26）
A、B的求取。
fAsin（2ft）Bcos（2ft）D
其中：
AAocos,BAosin
从而将被测信号的幅值和初始相角转化为对参数
其基本思想就是寻找合适的A、B和D的值，使得其测量残差的平方和取得最小。
设每个频率下测量的时间序列为ti-1,（i1,2,3n）,n为测量数据的个数，fs为采样
fs
频率，
fi为每个点的测量值。
则测量残差的表达式为：
n
2
[f（t）Asin（2fti）Bcos（2ftjD]（2-27）
i1
要使得上式取得最小值，可对其参数求偏导，并令其为零。
即:
A
A
n
[f（i）
i1
Asin（2fti）
Bcos（2
fti）
D]2
0
B
B
n
[f（i）
i1
Asin（2ftj
Bcos（2
fti）
D]2
0
（2-28）
D
D
n
[f（i）
i1
Asin（2fti）
Bcos（2fti）
D]2
0
进一步化简得到:
nn
sin（2ftJ[Asin（2ft」Bcos（2ftjD]sin（2ftjf（i）
i0i0
n
cos（2fti）[Asin（2ftj
i0
Bcos（2fti）
D]
n
cos（2fti）f（i）
i0
n
[Asin（2fti）
i0
Bcos（2fti）
n
D]
i0
f（i）
对于式（2-29），
构造如下三个矩阵
:
f1（tj
sin2ft1
sin2ft2
sin2ftn
A
T.
f1（t2）
cos2ft1
cos2ft2
・・・
cos2ftn
X
B,F
1
1
1
1
D
■■■
f1（tn）
（2-29）
式（2-29）可写成如下矩阵形式:
TXtf
上式中X的解为:
则被测信号幅值的计算公式为
初始相角的计算公式为:
2-11为参数拟合法的处理结果。
步便可得到每个频率下换能器的电导和电纳值。
图2-11正弦曲线参数拟合法处理结果
2.5导纳圆的带约束最小二乘曲线拟合
通过以上章节的介绍，我们已经得到了各个测试频率下压电换能器的电导和电纳值，绘
制出了导纳圆图，但这还是不够的。
由压电换能器等效电路的各个参数的计算公式可以看出，我们还需要得到导纳圆的圆心和半径的值。
为此，就需要对所得到的离散点进行圆曲线拟合。
拟合圆的方法有很多种，常用的有平均值法、加权平均法和最小二乘法。
平均值法的思想是分别计算各个离散点的横、纵坐标的平均值，作为圆心的横、纵坐标，
将圆心到各个离散点的距离的平均值作为半径。
这种方法计算简单，适用于离散点分布较均
匀的情况，但对于分布不均的情况，所计算的圆心位置会偏向离散点分布较密集的一侧，半
径的计算值也会偏小，误差较大。
加权平均法是对平均值法的改进，它在计算圆心坐标时加
入一个与两相邻点间弧长相关的系数，降低了离散点分布不均的影响，减小了误差。
但是，由于两相邻点间的弧长是无法精确得到的（实际中采用两点间的距离来代替的），所以其误
差仍较
升级会员 