1. 引言
压电材料作为一种力电转化的功能性材料,已被广泛应用于军事、医疗、工业和生活等方面 [1] [2] [3] [4] [5] ,尤其近些年来超声压电换能器的不断深入研究 [6] [7] [8] [9] [10] ,更使得压电材料应用得到了进一步的推广。目前,大功率压电器件的需求不断提高,然而大量的理论和实验表明,压电器件在大功率条件下会出现明显的非线性效应 [11] [12] [13] [14] ,材料参数变化明显,导致出现共振频率漂移、幅值曲线跳跃,以及器件过热、振幅饱和与性能下降等 [15] [16] ,这对大功率压电器件的应用十分不利。目前,压电材料表征侧重于弱场条件,大功率条件下压电材料的性能表征是一个亟待解决的问题。
机械品质因数Qm作为评估压电材料重要参数,其主要反映压电材料在谐振情况下克服内摩擦而消耗的能量。目前测量压电材料Qm值最常用的方法是定义法 [17] [18] 和3 dB (半功率点)法 [17] [19] ,这两种方法主要通过阻抗分析仪测量压电材料的阻抗谱来计算压电材料Qm值,但这两种方法的不足之处在于仅适用于弱场条件。目前大功率条件下压电材料Qm值测量的主要方法是恒定电流法 [20] [21] 和脉冲驱动法 [22] [23] ,其中恒定电流法可通过搭建相应的电路来测量恒定电流下的压电振子导纳曲线,以而求得相应的Qm值;脉冲驱动法主要是对压电振子进行脉冲串激励,通过分析在脉冲串激励结束后的振动速度衰减曲线来获得相应的Qm值。但上述两种方法的测量过程较为复杂,且在测量精度上也有待提高。
美国宾州州立大学K. Uchino教授课题组首次提出了利用温升和振动测量表征压电振子强场机械品质因数的方法 [24] ,但是对该表征方法的可行性并没有进行理论分析和验证,且以往的实验研究并没有给出Qm值与激励电场强度的关系。本文首先对薄长片压电振子的热传导方程进行有限差分求解分析;其次对推导得到的热传导方程进行数值模拟和有限元仿真,研究了薄长片压电振子的振动模态与应变和损耗的关系,模拟了不同材料参数下压电振子,验证利用温度和振动测量压电长条振子机械品质因数方法的准确性。在此基础上,搭建了实验装置对PZT压电材料的强场机械品质因数进行表征,得到了PZT-4和PZT-8在不同电场强度下的Qm值,验证了利用温升和振动测量表征压电振子机械品质因数的可行性。
2. 薄长片压电振子热传导方程的有限差分求解
本节主要对薄长片压电振子径向共振时的热传导方程进行有限差分求解分析,为以实验数据中拟合自发热参数hg建立理论基础。压电长条振子在共振时的发热主要由于应力应变滞后引起的机械损耗导致,而应变分布只与压电长条振子的长度有关,因此其产热分布也仅沿其长度方向发生变化。压电长条振子的几何形状以及热流分布如图1所示。
![](//html.hanspub.org/file/2-1270785x8_hanspub.png?20230915093836026)
Figure 1. Model of the thin piezoelectric plate and heat distribution in the length direction
图1. 薄长片压电振子模型以及在长度方向的热流分布
由图1所示,为进行热传导方程的有限差分求解,薄长片压电振子沿长度方向被平均分成N等分,图中的T1、Ti和TN分别代表长度被均分之后相应位置处的温度。我们可将压电振子径向振动时的热传导模型分成三部分去研究,即左边界(1位置)、右边界(N位置)以及左边界与右边界之间(1至N之间的位置)。根据图1中i点位置处的热流情况可写出如下的能量守恒方程:
(1)
其中Qc,i是传导进入微元i处的热流,而Qc,i+dx是由于传导流出微元i处的热流,Qg,i为微元i处的产热,Qd,i代表微元i处的散热,Ui代表微元i处的总内能,经推导可得到左边界到右边界之间微元处的热传导方程为 [24] :
(2)
其中λ为压电材料的热导率,cp和ρ分别为压电材料的比热容和密度,Tair为空气温度,P代表压电长条横截面周长,A代表压电长条横截面的面积,ε和σ分别表示为辐射率和斯蒂芬玻尔兹曼常数。图2给出了典型实验条件下对流与辐射耗散功率的比值与样品温度之间的关系。由图可知,随着温度不断下降,对流与辐射耗散功率的比值有所减小,但是两者的比值还是相对较大,因此辐射项可以忽略。
![](//html.hanspub.org/file/2-1270785x11_hanspub.png?20230915093836026)
Figure 2. The relationship between the ratio of convection and radiation power and the temperature of the material
图2. 对流与辐射功率比值与材料温度之间的关系
忽略辐射换热并对公式(2)进行有限差分求解可得到如下形式:
(3)
对于左边界和右边界单元可以写出如下的能量守恒方程:
(4)
(5)
对上述两个方程进行整理以同样的方式整理可得:
(6)
(7)
薄长片压电振子在共振时的产热主要由机械损耗引起,而且产热分布与压电元件应变的平方成正比,因此可写出如下方程:
(8)
其中sin(πx/L)为薄长片压电振子沿长度方向伸缩振动时应变分布函数 [25] ,hg为压电振子的发热参数。
3. Qm与机械能和损耗能之间的关系
机械品质因数Qm的值可定义为压电元件的单个周期内总的机械能与耗散的功率之比,其表达式如下:
(9)
其中Ue代表为压电元件存储的机械能,Pd代表耗散的功率。对于径向振动的薄长片压电振子,机械能可由其最大动能表示,且动能作为位移u的函数为
(10)
其中
(11)
上式中Vmax代表压电长条振子边界处的最大振动速度。将上式代入(10)式中可得最大动能为:
(12)
当薄长片压电振子达到稳态温度后耗散功率等于稳态时的产热功率,因此Pd可以表示为:
(13)
将(12)和(13)式代入到(9)式中可得到Qm与振动速度Vmax和发热参数hg之间的关系式如下:
(14)
4. 数值模拟与有限元仿真
本节建立了与实验样品相同尺寸的薄长片压电振子模型,利用数值模拟和有限元仿真计算其共振频率、应变分布、损耗分布以及温升分布,从而验证上述经有限差分求解后的热传导方程的合理性,同时通过有限元仿真验证了利用温升和振动法表征薄长片压电振子Qm值的可行性。本文模拟使用的材料为PZT-4,其材料性能参数如表1所示:
![](Images/Table_Tmp.jpg)
Table 1. The parameters of PZT-4 thin plate under numerical simulation
表1. 数值模拟时使用的PZT-4薄长片压电材料参数
首先利用有限元仿真软件COMSOL Multiphysics 5.6建立了薄长片压电振子模型,在压电振子上下两端施加电场,且施加自由边界条件,经特征频率计算,可得到其径向振动的特征频率为41.106 kHz,且其振动模态图如图3所示。通过数值软件MATLAB编程计算出的特征频率为40.884 kHz,这与通过有限元仿真得到的特征频率基本一致。通过仿真与数值计算得到的特征频率的数值之间有微小误差是由于COMSOL仿真软件内在计算原理与MATLAB数值计算的方法略有差别而造成的,因此可以忽略不计。
![](//html.hanspub.org/file/2-1270785x24_hanspub.png?20230915093836026)
Figure 3. Radial vibration mode diagram of PZT-4 thin piezoelectricvibrator
图3. PZT-4薄长片压电振子径向振动模态图
图4(a)给出了薄长片压电振子应变分布图,由图可知有限元仿真和数值模拟得到的结果基本重合,其应变分布沿长度方向为正弦曲线,且最大应变位置处在样品中心。图4(b)给出了有限元仿真得到薄长片压电振子的应变功率损耗分布图,将得到的应变分布进行平方处理后与功率损耗进行对比,可以看出两者曲线是相吻合的,从而证明公式(8)中薄长片压电振子的产热分布与应变分布的平方成正比是合理的。
![](//html.hanspub.org/file/2-1270785x25_hanspub.png?20230915093836026)
Figure 4. (a) strain distribution and (b) strain power loss distribution of the thin piezoelectric vibrator
图4. 薄长片压电振子的(a) 应变分布及(b) 应变功率损耗分布图
热导率作为传热过程中一个重要的参数,其对薄长片压电振子的温升分布有重要影响。图5给出了经数值模拟和有限元仿真得到的不同热导率时的压电振子的温升分布图。由图可知,随着热导率的增加,薄长片压电振子中心点的温度不断下降,而边界的温度不断增加,且其中心与边界的温差逐渐变小。这是因为热导率越高,压电振子中心热量更容易传导到产热较低的边界处。其次,薄长片压电振子应变最大处(中心点)的温升最大,而边界处的应变几乎为零,边界温升主要是由热量的传导引起的;最后,当热导率较小时薄长片压电振子的温度分布与应变功率耗散分布基本一致,进一步验证了产热分布与应变分布的平方成正比。
![](//html.hanspub.org/file/2-1270785x26_hanspub.png?20230915093836026)
Figure 5. Influence of different thermal conductivity on temperature rise distribution of the thin piezoelectric vibrator
图5. 不同热导率对压电长条振子温升分布的影响
为验证利用温升和振动数据表征薄长片压电振子Qm值的可行性,通过数值模拟和有限元仿真计算了不同机械损耗因子η时的温度分布,结果如图6所示,从图中可以看出,不同损耗因子时有限元仿真和数值求解得到的温度分布基本重合,从而进一步验证了推导得到的热传导方程及其有限差分求解方法是合理的。由图6和表2可知,η越大,温升越小,或者说发热参数越小,这是因为薄长片压电振子在共振时的发热与应变的平方成正比,且应变引起的损耗密度可表示为:
(15)
其中P1代表应变损耗功率密度,ω代表角频率,η代表机械损耗因子,ε代表应变向量,D代表弹性矩阵。当施加在薄长片压电振子两端的电场强度一定时,η越大,Qm越小,压电振子的振动速度会变小,其对应的应变也会变小。由(15)式可知,应变损耗密度与损耗因子η成正比,与应变的平方成正比,因此随着η的增加,薄长片压电振子的温升或者发热参数也随之变小。
![](Images/Table_Tmp.jpg)
Table 2. Comparison of theoretical value with simulated value of mechanical quality factor with different loss factors
表2. 不同损耗因子下机械品质因数理论值与模拟值对比
同时将有限元仿真得到的不同损耗因子下的振动速度以及温度分布进行计算得出相应条件下的Qm值,其结果如表2所示。由表2可知模拟得到的薄长片压电振子的Qm值与理论值的误差保持在2%以内,进一步说明了此方法的可行性。需要说明的是,上述表格中的发热参数是由压电振子温度分布拟合得到。
![](//html.hanspub.org/file/2-1270785x28_hanspub.png?20230915093836026)
Figure 6. Simulation results of temperature distribution of the thin piezoelectric vibrator with different loss factors
图6. 不同损耗因子薄长片压电振子的温度分布仿真结果
5. 强场下机械品质因数测量的实验研究
5.1. 实验装置及样品参数
![](//html.hanspub.org/file/2-1270785x29_hanspub.png?20230915093836026)
Figure 7. Schematic diagram of experimental apparatus
图7. 实验装置示意图
强场下薄长片压电振子机械品质因数测量的装置如图7所示,该实验装置主要包括函数发生器(Tek AFG1062)、功率放大器(NF HSA4011)、激光测振仪(PolyTec OFV-505)、高压差分探头(PINTECH DP-100)、红外热像仪(FLIR E8)、电流探针(TeKtronix A6032)和示波器(RIGOL MSO5072)。其中,函数发生器输出实验所需的交流电压信号值,通过功率放大器将电压放大并加载在样品上,红外热像仪测量样品侧面的温度分布,电流探针用于获取流经样品的真实电流值,高压差分探头和示波器用来检测样品两侧的实际电压和电流数值,而激光测振仪用来测量样品端面的振动速度。
实验过程中主要使用PZT-4与PZT-8薄长片压电材料两种样品,其中PZT-4和PZT-8薄长片压电振子长宽高分别为40 mm,6 mm和2 mm。
5.2. 实验测量与结果
![](//html.hanspub.org/file/2-1270785x30_hanspub.png?20230915093836026)
Figure 8. Impedance spectra of PZT-4 and PZT-8 thin piezoelectric vibrators
图8. PZT-4与PZT-8薄长片压电振子阻抗谱
在实验中激励电压频率为样品的一阶共振频率,首先通过阻抗分析仪测量出两种样品的阻抗谱,其中图8为样品PZT-4与PZT-8的阻抗曲线,从阻抗谱中可知两种样品的一阶共振频率分别为39.52 kHz和41.96 kHz。
要准确拟合发热参数hg,首先要确定样品的对流换热系数,在实验主要是通过拟合样品达到稳态温度后停止激励时样品自然冷却到室温的温度数据,从而计算出样品的平均对流换热系数 [26] [27] ,图9给出了PZT-4和PZT-8振子的在电压幅值为30 V时对流换热系数测量拟合情况。样品PZT-4的平均对流换热为20.3 W/m2·K,样品PZT-8的平均对流换热系数为21.8 W/m2·K,且测量结果与文献 [28] 中测量得到的空气对流换热系数比较接近。
在完成共振频率和对流换热系数的测量后,可进行自发热温升测量,需要注意的是,在强场激励时压电振子会出现弹性非线性,使共振频率发生漂移,因此在不同电压幅值激励时要以共振频率为基准微调其激励频率,使电流也达到了最大值,如此为真实的共振频率,激励样品达到稳态温度,记录下样品的稳态温度分布和样品边界的振动速度,之后改变电压幅值并重复此过程。通过对上述数值模型进行编程,将所测实验数据与理论数值进行拟合便可达到相应电压幅值下的发热参数hg,将得到的发热参数与振动速度代入到(14)式便可得到相应电场下的机械品质因数。
![](//html.hanspub.org/file/2-1270785x31_hanspub.png?20230915093836026)
Figure 9. Simulation results of convective heat transfer coefficients of (a) PZT-4 and (b) PZT-8 piezoelectric vibrators
图9. (a) PZT-4与(b) PZT-8压电振子对流换热系数的拟合结果
![](//html.hanspub.org/file/2-1270785x32_hanspub.png?20230915093836026)
Figure 10. Schematic diagram of mechanical quality factor variation with electric field
图10. 机械品质因数随电场变化示意图
图10给出了PZT-4与PZT-8薄长片压电振子机械品质因数Qm随加在其两端的电场强度变化的实验结果,从图中可以看出:PZT-8的机械品质因数Qm总体高于PZT-4,且随着两端电场强度的增大,两种材料的机械品质因数Qm值变小。从宏观角度来说,在实验过程中可以明显的观察到随着电场增大,薄长片压电振子C的温升增加较快,而振动速度相对于温升来说增加较慢,所以就会导致Qm值变小。从微观角度来分析,畴壁的运动的难易程度直接影响着机械品质因数Qm的大小,在实验中随着加样品在两端的电场强度不断增大,薄长片压电振子的温度不断升高,两者共同作用促使畴壁的运动,使得Qm值变小。
需要说明的是:由于弱场激励时温升较小,图10中1 V/mm的测量结果是通过测量样品的阻抗谱并用3 dB法计算得到的。此外,由于温度过高会使得压电材料的各项性能参数发生较大的变化,因此在进行上述测试时温升一般不应超过20℃。
6. 结论
本文首先对薄长片压电振子的热传导方程进行有限差分建模和求解,并利用有限元仿真和数值求解验证了此方程的合理性;其次求解计算了不同参数时薄长片压电振子的温升分布图,验证了薄长片压电振子共振时应变与损耗之间的关系以及利用温升和振动法表征强场下薄长片压电振子机械品质因数Qm值的可行性;最后搭建实验装置表征了PZT-4和PZT-8薄长片压电振子的Qm值,结果表明:两种材料的机械品质因数随电场的增大而变小,并且相同电场强度下PZT-8的机械品质因数Qm值均大于PZT-4的Qm值。
基金项目
国家自然科学基金资助项目(11874327)。
NOTES
*通讯作者。