高应变率下多孔未极化PZT95/5铁电陶瓷的非线性力学行为*

李成华，蒋招绣，王贝壳，张 振，王永刚

(宁波大学冲击与安全工程教育部重点实验室，浙江 宁波 315211)

Pb(Zr0.95Ti0.05)O3(PZT95/5)铁电陶瓷处于铁电相与反铁电相相界附近，具有丰富的相结构、外场诱导畴变和相变特性[1]。在外加电场作用下，PZT95/5铁电陶瓷容易被极化，而且撤除电场后仍然保留较强的剩余极化强度。在冲击波加载下，极化PZT95/5铁电陶瓷发生相变而去极化，从而迅速释放表面束缚的电荷，在外电路中产生脉冲高电流或高电压[2]。因此，PZT95/5铁电陶瓷常被用来制备小巧但高功率的脉冲电源，在工业和国防科技领域具有广泛的应用前景[3]。

针对冲击波加载下PZT95/5铁电陶瓷的电学响应、力学响应以及力电耦合响应，人们开展了大量研究工作，集中探讨了极化状态、外载荷条件(包括力学、电学及力电耦合加载)以及负载电阻等因素的影响[4-6]。除了冲击波加载，人们还对准静态单轴压缩和静水压缩下PZT95/5铁电陶瓷的相变行为和去极化行为开展了许多工作，重点讨论了应力状态的影响[7-9]。近年来，人们通过添加造孔剂制备出具有不同孔隙微结构特征的PZT95/5铁电陶瓷，研究发现：与密实PZT95/5铁电陶瓷相比，在冲击波压缩下多孔PZT95/5铁电陶瓷具有更好的界面阻抗匹配性、更高的抗电击穿强度和更低的相变压力，但是随着孔隙率增大，PZT95/5铁电陶瓷的放电性能稍有弱化[10-11]。准静态单轴压缩下多孔PZT95/5铁电陶瓷的畴变应力、相变压力、压缩强度都随着孔隙率的增大而减小，但断裂应变基本不依赖于孔隙率[12-13]。一般来说，脆性材料的力学特性对加载应变率非常敏感。目前，有关冲击压缩下多孔PZT95/5铁电陶瓷力学性能及其应变率敏感性的研究工作比较匮乏。

本文中，采用添加造孔剂的方法制备多孔未极化PZT95/5铁电陶瓷；采用改进的分离式霍普金森压杆(split Hopkinson pressure bar, SHPB)测试系统，对多孔PZT95/5铁电陶瓷开展一系列高应变率单轴压缩；采用超高速相机实时监测PZT95/5铁电陶瓷试样在动态单轴压缩下的变形与破坏，再结合数字图像相关性(digital image correlation, DIC)方法得到试样动态全场应变信息，分析SHPB实验中试样应变均匀性问题，重点讨论应变率和孔隙率对PZT95/5铁电陶瓷非线性力学行为的影响。

1 实 验

1.1 试样制备

以聚甲基丙烯酸甲酯(PMMA)作为造孔剂制备多孔PZT95/5铁电陶瓷。用球磨法将质量分数(wPMMA)为0、1%、2%和3%的PMMA分别与PZT95/5粉末球磨6 h(以锆球作为球磨介质，酒精作为溶剂，wPMMA=0表示未添加PMMA)。球磨后的浆料于100 ℃下干燥24 h，在干燥后的粉末中添加质量分数为6%的PVA液体黏合剂，混合均匀后在200 MPa压力作用下压缩成型。成型的材料在850 ℃下对PMMA进行排塑，再将排塑后的素胚放入富含PbO环境的密封氧化铝坩埚中，并在1 300 ℃下烧结2 h。对烧结好的样品表层进行打磨和清洗，最终得到∅6 mm×6 mm和∅30 mm×6 mm两种尺寸的具有不同孔隙率的未极化PZT95/5铁电陶瓷试样。

制备出4种不同孔隙率的未极化PZT95/5铁电陶瓷试样，分别标记为1#、2#、3#和4#。试样的低倍率(500倍)扫描电镜(SEM)形貌如图1所示。可见，造孔剂在试样中产生了较均匀的球形微孔洞，微孔洞数量随着造孔剂添加量的增大而增多。图2显示了典型的球形微孔洞局部放大SEM图像，球形微孔洞的直径约为30 μm，平均晶粒尺寸约12 μm。

1.2 实验方法

式中：εr(t)和εt(t)分别为入射杆上实测的反射波应变和透射杆上实测的透射波应变，Eb、Ab和c0分别为杆的弹性模量、截面面积和弹性波波速，l和As分别为试样的原始长度和截面面积。

对于脆性材料PZT95/5铁电陶瓷而言，SHPB实验的难点在于其破坏应变小，在有效测量时间内，试件难以达到应力/应变均匀性要求。有关试件应力均匀性问题的讨论已开展非常多[15-16]。理论分析[17]指出，弹性波在试样中传播4个来回后，试件两端的相对应力差将低于5%，此时可近似认为试件达到了应力均匀化要求。事实上，除了关注试件应力均匀性问题，应变均匀性更重要，尤其对脆性材料而言。为了监测试件应变均匀性，本研究中在传统SHPB实验装置上搭建了基于超高速相机与DIC分析技术的试件全场应变测量系统。其中，超高速相机采用英国Specialised-imaging公司的Kirana相机(最高拍摄速度为5×106帧/秒，图像分辨率为924×768像素，可连续拍摄180张)，在SHPB实验中的拍摄速度为1×106帧/秒，曝光时间为1 μs。对实时拍摄的照片进行DIC分析即可获得试样全场应变信息。DIC技术[18]是一种材料表面变形场测量技术，其基本原理为：首先获取试样表面散斑图像，散斑图像上的所有特征点以像素点为坐标，以像素灰度作为信息载体；然后，在试样表面划分若干个图像子区，类似于有限元计算中划分网格，在图像移动或变形过程中，基于图像子区灰度值不变的假定，追踪每个图像子区在变形后图像的位置，即可获得所有图像子区中心点处的位移矢量，从而获得物体表面的变形信息。实验前，在试件外表面喷涂哑光黑漆和哑光白漆以制作高质量散斑。

在SHPB实验中，采用波形整形技术实现恒应变率加载[19]。这里采用无氧铜垫圈(内、外径分别为4和8 mm，厚度为1.5 mm)作为波形整形器，将矩形入射波调整为近似三角形入射波。图4给出了实测的入射波、反射波和透射波的原始电压信号。可见，在试样发生破坏前，反射波出现明显的平台段，表明实现了试样的恒应变率加载。另外，为了对比分析DIC全场应变测量的有效性，选取一发实验，在试样中部粘贴一片应变片，通过应变片直接测量试样应变，实测的电压信号也在图4中显示。

2 实验结果与讨论

2.1 应变均匀性讨论

采用DIC方法得到某时刻试样全场应变云图和不同时刻试样沿轴向的应变分布，如图5和图6所示，可以看出：(1) 试样沿轴向的应变分布并不均匀，两端应变明显增大，甚至超过中部应变的两倍，这归因于与杆接触的试样端部存在应力集中；(2) 试样中部约3 mm长度范围内的应变分布相对均匀。

图7给出了应变片(粘贴在试样中部)实测应变时程与DIC方法得到的试样中部平均应变时程的对比。两者的一致性很好，验证了DIC全场应变测量方法的有效性。另外，图6也给出了采用SHPB原理(式(2))计算得到的试样应变。可见，计算结果明显偏大，原因在于：SHPB原理成立的前提是试样满足应变均匀性条件，而对于脆性材料PZT95/5，应力集中所引起的试样应变不均匀性比较严重。因此，为了准确获取PZT95/5铁电陶瓷的动态应力-应变关系曲线，选择试样中部DIC平均应变或应变片实测值作为试样应变比较合理，必须摒弃或修正SHPB实验所采用的传统数据处理方法。

2.2 应变率的影响

图8给出了多孔未极化PZT95/5铁电陶瓷(孔隙率φ约为11%)在不同应变率下的应力-应变曲线，其中准静态实验的详细介绍参见文献[12]。与准静态结果相比，高应变率下PZT95/5铁电陶瓷的压缩强度和断裂应变都显著提高，暗示多孔未极化PZT95/5铁电陶瓷的力学行为具有强烈的应变率效应。然而，SHPB实验中脆性材料压缩强度的显著提高是否可简单地归因于材料的应变率效应一直存在争议。众所周知，在SHPB实验中，试样在轴向快速变形的同时，径向也快速向外膨胀；径向膨胀会在试样中产生惯性应力，惯性应力相当于给试样额外施加了惯性侧向约束。许多实验结果已揭示出脆性材料的动态压缩强度对侧向约束应力非常敏感。由此看来，在SHPB实验中径向惯性效应与材料内禀的应变率效应是相互耦合的。Li等[20]认为混凝土或岩石类脆性材料的动态压缩强度显著增大主要是由径向惯性引起的侧向约束所导致，不能错误地归因于材料本构关系中的应变率效应。基于弹性变形假设，Forrestal等[21]给出了径向惯性引起圆柱形试样惯性应力的理论计算公式：

除了关注动态压缩强度，从图8中还观察到无论准静态压缩还是冲击压缩，多孔未极化PZT95/5铁电陶瓷都具有明显的非线性变形特征。通常认为，畴变是PZT铁电陶瓷材料宏观非线性行为的内在机理[22]。下面基于试样全场应变测量结果(由DIC方法获得)，讨论多孔未极化PZT95/5铁电陶瓷非线性变形的内在机理，重点关注应变率的影响。图11给出了典型的准静态加载(10-4s-1)和冲击加载(220 s-1)下试样轴向应变和径向应变随着加载应力的变化曲线(以压缩应变为正)。在加载初期，多孔未极化PZT95/5铁电陶瓷发生弹性变形，轴向应变和径向应变随着轴向应力的增大呈线性增长；在B、B′点和E、E′点之后，轴向应变和径向应变开始偏离初始直线，进入非线性变形阶段，此时PZT95/5铁电陶瓷发生畴变，畴变改变晶格形状，同时促进轴向应变和径向应变快速增长，B点和E点对应的应力可定义为畴变临界应力；随着压缩应力进一步增大，轴向应变继续增长，但径向应变在C′点和F′点出现反转，即轴向压缩下试样从径向膨胀转变为径向收缩，这种异常现象显然无法用畴变解释。在压缩应力的作用下，处于FE-AFE相界附近的PZT95/5铁电陶瓷不仅会发生畴变，还会发生相变。在菱方晶系FE相转变为正交斜方晶系AFE相的相变过程中，体积将收缩，从而导致试样发生径向收缩，因此C′点和F′点所对应的应力即为相变临界应力。由此看来，多孔未极化PZT95/5铁电陶瓷非线性变形的内在机理是畴变和相变的共同作用。从图11容易确定：在10-4s-1应变率下，多孔未极化PZT95/5铁电陶瓷的畴变临界应力σqd和相变临界应力σqf分别为104和136 MPa；而当应变率提高到220 s-1时，畴变临界应力σdd和相变临界应力σdf分别为207和289 MPa。由此看来，σdd和σdf都随着应变率的增加而增大。按照上述方法，另外获取了其他应变率下的实验结果，并以10-4s-1应变率下的实验结果作为参考进行归一化处理，处理后σdd和σdf随应变率的变化曲线如图12所示。可见，可以采用统一公式对曲线进行拟合，即：

2.3 孔隙率的影响

采用高精度电子密度计对每个试样的密度进行测量。由密实PZT95/5铁电陶瓷密度(8.0 g/cm3)以及试样实测密度，换算出每个试样的孔隙率。根据孔隙率大小进行筛选，选取孔隙率约为5%、11%、15%、20% 4种试样。在近似相同的应变率(约200 s-1)下对4种试样进行3组重复性SHPB实验，所得动态应力-应变曲线如图13所示。可以看出：在试样变形初期，实验结果的重复性较好；而在变形后期，实验结果的重复性变差，特别是断裂强度。虽然实验结果具有一定分散性，但是对于陶瓷类脆性材料而言，图13所示实验结果的重复性是可以接受的。下面基于这些实验数据，讨论孔隙率对未极化PZT95/5铁电陶瓷动态力学性能的影响。

孔隙率对未极化PZT95/5铁电陶瓷的动态压缩强度有显著的影响，如图14所示，其中还给出了准静态压缩强度随孔隙率的变化关系。图14显示：准静态压缩强度随着孔隙率的增大线性衰减，而动态压缩强度则随着孔隙率的增大呈非线性衰减。还有一点值得注意，孔隙率为20%的未极化PZT95/5在高应变率下的压缩强度与准静态实验结果基本接近，没有表现出应变率效应。事实上，脆性材料的宏观压缩强度在细观上取决于试样中与拉伸应力相关的微裂纹成核和扩展。对于内部含有较多微孔洞的高孔隙率试样，微裂纹一般成核于微孔洞周围，微孔洞之间容易借助微裂纹扩展而相互贯通，形成贯穿性宏观裂纹，最终引起试样在较低应力下发生破坏。

采用2.2节中关于畴变临界应力和相变临界应力的确定方法，根据图13中的实验结果，得到未极化PZT95/5铁电陶瓷的畴变临界应力σdd和相变临界应力σdf随着孔隙率φ的变化曲线，如图15所示。可以看到，σdd和σdf均随着φ的增大基本上呈线性衰减，线性拟合后得到直线在纵轴上的截距，分别定义为基体材料的畴变临界应力σdd0(305 MPa)和相变临界应力σdf0(514 MPa)。以σdd0和σdf0为基础进行归一化处理后，发现可以采用统一的线性表达式拟合数据，即：

3 结 论

利用超高速相机结合DIC方法，在SHPB实验中实时获取了试样的全场应变信息，发现脆性材料PZT95/5铁电陶瓷的轴向应变在试样端部的分布极不均匀，只在中部分布较均匀，中部的平均应变与应变片的实测值一致，而由SHPB原理计算的试样应变则明显偏大。在高应变率单轴压缩下，多孔未极化PZT95/5铁电陶瓷表现出明显的应变率效应和非线性变形特征。由轴向应变及径向应变与轴向应力之间的关系，揭示出未极化PZT95/5铁电陶瓷非线性变形的物理机理是畴变和相变的共同作用。未极化PZT95/5铁电陶瓷压缩强度的动态增强因子随着应变率的增大呈幂指数增长，归一化的畴变临界应力和相变临界应力也随着应变率的增大呈幂指数增长。未极化PZT95/5铁电陶瓷的动态压缩强度随着孔隙率的增大呈非线性衰减，而畴变临界应力和相变临界应力则随着孔隙率的增大呈线性衰减。

参考文献:

[1] 钟纬烈.铁电体物理学[M].北京:科学出版社,1996.

[2] 杜金梅,张毅,张福平,等.冲击加载下PZT95/5铁电陶瓷的脉冲大电流输出特性[J].物理学报,2006,55(5):2584-2589.

DU Jinmei, ZHANG Yi, ZHANG Fuping, et al. Large current out-put of PZT 95/5 ferroelectric ceramics under shock loading[J]. Acta Physica Sinica, 2006,55(5):2584-2589.

[3] LYSNE P C, PERCIVAL C M. Electric energy generation by shock compression of ferroelectric ceramics: Normal-mode response of PZT 95/5[J]. Journal of Applied Physics, 1975,46(4):1519-1525.

[4] SETCHELL R E. Shock wave compression of the ferroelectric ceramic Pb0.99(Zr0.95Ti0.05)0.98Nb0.02O3: Depoling currents[J]. Journal of Applied Physics, 2005,97(1):013507.

[5] SETCHELL R E. Shock wave compression of the ferroelectric ceramic Pb0.99(Zr0.95Ti0.05)0.98Nb0.02O3: Hugoniot states and constitutive mechanical properties[J]. Journal of Applied Physics, 2003,94(1):573-588.

[6] ZHANG F P, HE H L, LIU G M, et al. Failure behavior of Pb(Zr0.95Ti0.05)O3ferroelectric ceramics under shock compression[J]. Journal of Applied Physics, 2013,113(18):183501.

[7] ZEUCH D H, MONTGOMERY S T, HOLCOMB D J. Uniaxial compression experiments on lead zirconate titanate 95/5-2Nb ceramic: Evidence for an orientation-dependent, “maximum compressive stress” criterion for onset of the ferroelectric to antiferroelectric polymorphic transformation[J]. Journal of Materials Research, 2000,15(3):689-703.

[8] ZEUCH D H, MONTGOMERY S T, HOLCOMB D J. The effects of nonhydrostatic compression and applied electric field on the electromechanical behavior of poled lead zirconate titanate 95/5-2Nb ceramic during the ferroelectric to antiferroelectric polymorphic transformation[J]. Journal of Materials Research, 1999,14(5):1814-1827.

[9] AVDEEV M, JORGENSEN J D, SHORT S, et al. Pressure-induced ferroelectric to antiferroelectric phase transition in Pb0.99(Zr0.95Ti0.05)0.98Nb0.02O3[J]. Physical Review B, 2006,73(6):064105.

[10] FENG N, NIE H, CHEN X, et al. Depoling of porous Pb0.99(Zr0.95Ti0.05)0.98Nb0.02O3ferroelectric ceramics under shock wave load[J]. Current Applied Physics, 2010,10(6):1387-1390.

[11] SETCHELL R E. Shock wave compression of the ferroelectric ceramic Pb0.99(Zr0.95Ti0.05)0.98Nb0.02O3: Microstructural effects[J]. Journal of Applied Physics, 2007,101(5):053525.

[12] 蒋招绣,辛铭之,申海艇,等.多孔未极化Pb(Zr0.95Ti0.05)O3铁电陶瓷单轴压缩力学响应与相变[J].物理学报,2015,64(13):134601.

JIANG Zhaoxiu, XIN Mingzhi, SHEN Haiting, et al. Mechanical properties and phase transformation of porous unpoled Pb(Zr0.95Ti0.05)O3ferroelectric ceramics under uniaxial compression[J]. Acta Physica Sinica, 2015,64(13):134601.

[13] 蒋招绣,申海艇,辛铭之,等.多孔极化PZT95/5铁电陶瓷单轴压缩力学响应与放电特性[J].固体力学学报,2016,37(1):50-57.

JIANG Zhaoxiu, SHEN Haiting, XIN Mingzhi, et al. Mechanical properties and depoling of porous poled PZT95/5 ferroelectric ceramics under uniaxial compression[J]. Chinese Journal of Solid Mechanics, 2016,37(1):50-57.

[14] 王礼立.应力波基础[M].2版.北京:国防工业出版社,2010.

[15] YANG L M, SHIM V P W. An analysis of stress uniformity in split Hopkinson bar test specimens[J]. International Journal of Impact Engineering, 2005,31(2):129-150.

[16] 朱珏,胡时胜,王礼立.SHPB试验中粘弹性材料的应力均匀性分析[J].爆炸与冲击,2006,26(4):315-322.

ZHU Jue, HU Shisheng, WANG Lili. Analysis on stress uniformity of viscoelastic materials in split Hopkinson bar tests[J]. Explosion and Shock Waves, 2006,26(4):315-322.

[17] 毛勇建,李玉龙.SHPB试验中试件的轴向应力均匀性[J].爆炸与冲击,2008,28(5):448-454.

MAO Yongjian, LI Yulong. Axial stress uniformity in specimens of SHPB tests[J]. Explosion and Shock Wave, 2008,28(5):448-454.

[18] SCHREIER H, ORTEU J J, SUTTON M A. Image correlation for shape, motion and deformation measurements[M]. New York: Springer, 2009.

[19] FREW D J, FORRESTAL M J, CHEN W. Pulse shaping techniques for testing brittle materials with a split Hopkinson pressure bar[J]. Experimental Mechanics, 2002,42(1):93-106.

[20] LI Q M, MENG H. About the dynamic strength enhancement of concrete-like materials in a split Hopkinson pressure bar test[J]. International Journal of Solids and Structures, 2003,40(2):343-360.

[21] FORRESTAL M J, WRIGHT T W, CHEN W. The effect of radial inertia on brittle samples during the split Hopkinson pressure bar test[J]. International Journal of Impact Engineering, 2007,34(3):405-411.

[22] 杨卫.力电失效学[M].北京:清华大学出版社,2000.