功能梯度压电材料半空间中的Rayleigh波.docVIP

功能梯度压电材料半空间中的Rayleigh波 曹小杉，金峰，王子昆 西安交通大学航天航空学院，西安 710049 电子邮箱: jinzhao@mail.xjtu.edu.cn，传真: 86-29电话: 86-29摘要 对于材料参数沿深度方向连续变化的横观各向同性压电材料半无限体中Rayleigh波的传播问题，应用Wentzel-Kramers-Brillouin (WKB) 方法得到了渐进解析解答。对于自由表面电学开路和电学短路边界条件，得到了问题的相速度方程。讨论了压电材料参数梯度变化对Rayleigh波频散特性、自由表面质点位移幅值及相应衰减特性的影响。研究结果表明：在压电梯度半无限体的表面可传播耦合的Rayleigh面波，其频散形式和自由表面质点位移椭圆度取决于材料参数的变化趋势，衰减特性与均匀材料相似，梯度变化对传播性能的影响在电学开路和电学短路两种情况下基本一致。通过Rayleigh面波相速度变化量，梯度变化参数与波数的关系，可为借助实验测量确定其材料参数梯度变化参数提供理论依据。 关键词: 功能梯度压电材料；瑞利表面波；频散特性； WKB 方法; 材料性能表征 1 引言 功能梯度压电材料（FGPM）是借鉴功能梯度概念制备的压电材料，由于其固有的机电耦合性质[1]，已经被广泛应用于各个领域。目前已有的研究多限于静力学问题[2]或振动问题[3]。随着动力测量、超声波技术、无损检测和人工智能等高科技领域的迅速发展，研究这种介质中的弹性波的传播问题极为必要。 研究连续变化的非均匀介质及其结构中的弹性波传播问题，由于控制方程的复杂性，使得问题求解相当困难，只有对于某些特殊的情形，可以求得解析解。通常采用的数值解法，是将材料划分成许多层，每层中的材料参数设为均匀，再采用有限元法[4-8]，传递矩阵法[9]等进行计算。在已有的研究报道中，也有采用WKB法研究复杂介质中导波传播特性的，如功能梯度层状结构的love波传播[10-11]，沿曲面传播的Rayleigh波等[12]，而文献[13]则采用摄动法得出功能梯度材料半空间Rayleigh面波的波函数解答和频散曲线。 对于材料参数沿厚度变化的横观各向同性梯度压电半无限体Rayleigh面波，其控制方程是无法解耦的三元二阶微分方程组。本文利用待定系数法将其变换为含一个未知函数的六阶变系数微分方程，再采用WKB法求解其波函数解答与自由表面分别为电学短路和电学开路情况的相速度方程。以具体压电材料为例，在材料参数沿深度方向按指数规律连续变化的情况下，分析了压电材料各项性能参数梯度变化对Rayleigh面波频散性能，自由表面质点位移幅值和衰减性能的影响。为借助Rayleigh面波的实验测量确定其材料参数梯度变化指数提供了理论依据。 2 基本方程 设横观各向同性梯度压电介质半无限体如图1所示，表面外力自由，极化方向沿轴。假定弹性系数、压电系数及介电系数是坐标z的确定函数。 考察上述半无限体耦合Rayleigh面波的传播特性，可设波沿x轴正向传播，位移函数和电势函数可表示如下： （1） 问题的控制方程可以表示为： （2） 这里“′”表示对坐标z的导数。 一般来说，半无限体上表面空间通常是气体介质，电势函数应当满足： （3） 自由表面电学边界条件考虑电学开路和电学短路两种情况。 Figure 1 功能梯度压电材料半空间无限体 3 问题求解 控制方程（2） （4） 其中为波数，为波的相速度。，和为深度的待定常数。问题的控制方程是无法解耦的三元二阶微分方程组本文利用待定系数法将其变换为含一个未知函数的六阶变系数微分方程，然后再借助于WKB法求解其波函数解答与自由表面分别为电学短路和电学开路的相速度方程 其中为梯度指数，并定义材料参数相对变化量为 和分别代表和时的性能参数，算例中取值为PZT-1陶瓷的材料参数。 4. 1 材料的非均匀性对弥散特性的影响 为了揭示材料的非均匀性对弥散效应的影响，分析中采用一个材料参数变化而其他材料参数不变的简化方案。详细研究了不同材料参数变化对Rayleigh波频散曲线的影响。数值算例中分别选取材料参数相对变化量，梯度参数分别为1， 3，5。 结果表明，无论是材料性能参数递增变化，还是递减变化，在梯度压电半无限体表面均可传播耦合Rayleigh面波，