导电压头作用下的功能梯度压电涂层二维黏附接触问题研究.pdf

应用数学和力学编委会，：导电压头作用下的功能梯度压电涂层二维黏附接触问题研究韩立夫，刘铁军（内蒙古工业大学 理学院，呼和浩特；河套学院，内蒙古 巴彦淖尔）摘要：纳米压痕实验是研究材料的力学性能和表面形貌的重要手段，当接触区尺寸减小时，压头与试件接触表面间的黏附作用将无法忽视，因此，考虑黏附作用对压头作用下的接触问题具有重要的价值 功能梯度压电材料（）兼具梯度材料和压电材料的优点，用作涂层可有效地抑制接触损伤和破坏 该文将针对梯度压电材料在导电压头作用下的黏附接触问题开展研究，假设功能梯度压电涂层的材料参数按照指数形式变化，基于 黏附模型，利用 积分变换获得了功能梯度压电涂层在导电压头作用下的二维无摩擦黏附接触问题的控制奇异积分方程，并采用 的数值方法求解，获得了黏附应力、梯度参数和压头所带电荷对力电耦合响应的影响 研究结果为利用功能梯度压电材料涂层改善材料表面的接触行为提供了理论依据，同时可为压电结构及器件的设计提供帮助 关 键 词：功能梯度压电涂层；黏附；积分变换；奇异积分方程中图分类号：文献标志码：，（，；，）：，（），应用数学和力学 卷 期 年 月 ，收稿日期：；修订日期：基金项目：国家自然科学基金（）；内蒙古自然科学基金（）作者简介：韩立夫（），男，讲师，博士生（：）；刘铁军（），男，教授，博士，博士生导师（通讯作者：）引用格式：韩立夫，刘铁军 导电压头作用下的功能梯度压电涂层二维黏附接触问题研究 应用数学和力学，（）：，：；引 言随着微机电系统和仿生器械的快速发展，对于结构和构件的黏附接触及损伤的研究已经引起了科学家和工程师们的广泛关注 对于均匀弹性材料的黏附接触问题，学者们先后建立了 刚体模型、模型、模型、（）模型、双 模型 近年来，学者们将这些经典模型发展到新型复合材料中，并研究其黏附接触行为 等研究了压电半空间球面刚性冲头的微尺度黏附接触问题，结果表明，压电效应对压电材料的黏附接触行为有显著影响；等研究了球形压头与弹性层状介质之间的无摩擦黏附接触问题；等建立了幂型梯度材料的黏附接触模型，给出了黏附接触界面临界拉脱半径的解析解，研究结果表明，功能梯度材料的黏附接触拉脱力与材料的弹性模量无关，但依赖于材料的梯度变化、球体半径及黏附能；和 等相继建立了幂型梯度材料的轴对称无摩擦 黏附接触模型和受表面粗糙度影响的黏附接触模型以及双 黏附接触模型，并将黏附接触模型扩展到了压电材料，以上研究更详细的内容可参见文献 有学者将功能梯度材料的设计理念引入到压电材料，从而为功能梯度压电材料（）的制备提供了设计思路，已经取得的大量研究成果表明，将 用作涂层能有效改善均匀压电材料接触表面的力学性能和损伤，并实现可设计性 等研究了利用 作为均匀压电材料表面的涂层，有效抑制了器件使用过程中的破坏行为；等研究了参数随指数变化的 与刚性绝缘和导电压头的二维无摩擦接触问题，研究结果表明，材料梯度指数和压头特性对功能梯度压电涂层的接触力学性能产生显著影响；等对材料参数呈指数变化的功能梯度压电涂层在绝缘与导电压头作用下的轴对称无摩擦接触问题进行了研究，发现导电压头作用下的最大接触应力值小于绝缘压头下的值，同时还研究了绝缘压头作用下的功能梯度压电涂层部分滑移接触问题；等深入研究了功能梯度压电涂层在导电压头作用下的部分滑移接触问题；刘兴伟等研究了一维六方压电准晶中正 边形孔边裂纹的反平面问题；马占洲等基于层合板模型研究了梯度压电涂层 问题；代文鑫等研究了导电压头作用下的多层 涂层二维接触问题，更详细的内容可参见文献及相关文献 特别地，等建立了平行弹性长圆柱体（均匀弹性材料）之间的 黏附接触模型，并给出了一个参数 来控制不同接触理论的适用范围；等建立了功能梯度材料相关的 黏附接触模型，并比较系统地研究了功能梯度涂层二维和轴对称黏附接触问题，同时还考虑了尺度效应，详细的研究内容可参见文献 本文在功能梯度压电涂层二维接触问题基本解的基础上，利用 黏附理论，建立了刚性圆柱导电压头与功能梯度压电涂层二维黏附接触模型，给出了刚性圆柱导电压头作用下，功能梯度压电涂层二维无摩擦黏附接触问题的控制方程，并转化为 奇异积分方程 标准化后，采用 的方法进行了数值计算，并定量分析了黏附应力、梯度参数和压头所带电荷对拉脱力、接触应力、电荷分布、压痕及电势等力电参数的影响 本文的研究丰富了 黏附接触理论，对解决功能梯度压电涂层二维黏附接触问题具有一定的理论指导意义 功能梯度压电涂层基底结构二维问题的通解考虑厚度为 的功能梯度压电涂层与均匀压电基底半空间完美黏接，轴位于涂层和基底之间的界面处，轴沿竖直方向向上，涂层上表面（，）处受法向集中线载荷、切向集中线载荷以及正集中线电荷 的共同作用，该问题的力学模型如图（）所示，并假设沿厚度方向极化压电材料，功能梯度压电涂层应 用 数 学 和 力 学 年 第 卷的材料参数沿 轴方向呈指数形式变化，即（），（），（）（，），（）其中，和 分别代表功能梯度压电涂层与均匀压电基底半空间在 处的弹性常数、压电常数和介电常数；代表涂层内材料参数的梯度指数，代表涂层和基底是同种均匀压电材料 关于该问题通解的具体推导过程可参见附录，后文中方程推导和数值计算将会用到（）法向、切向集中线载荷、和正集中线电荷 共同作用（）刚性导电压头作用（），（）图 功能梯度压电涂层压电半空间接触力学分析模型 如图（）所示，假设压头是刚性的，其上作用有法向集中线载荷 和正集中线电荷，并按要求将其安置在 涂层和均匀压电基底半空间上，在 涂层表面将会形成接触区 和电势；在接触区外（），假设接触应力和电荷分布均为零，不计摩擦 此问题可用 涂层和均匀压电基底半空间二维无摩擦接触问题的基本解（附录 中式（）（）进行求解 假设（）和（）（）分别为表面接触区内的法向接触压力和电荷分布，即（，）（）和（，）（），利用叠加原理，对附录 中式（）（）在接触区内积分，可以获得刚性导电压头作用下接触表面的位移分量和电势表达式：（，）（）（）（）（）（）（）（）（），（）（，）（）（）（）（）（）（），（）（，）（）（）（）（）（）（）（）刚性圆柱导电压头作用下无摩擦黏附接触问题的求解 建立黏附接触模型为了便于问题的求解，压头选用圆柱型，并假设 涂层和均匀压电基底半空间在刚性圆柱导电压头作用下的二维无摩擦黏附接触满足 黏附理论，如图 所示的 涂层在导电压头作用下的力学模型 图（）为黏附接触模型，图（）为涂层表面的应力分布 根据 黏附理论，其黏附功可表示为 第 期 韩立夫，等：导电压头作用下的功能梯度压电涂层二维黏附接触问题研究，当两个接触物体表面间的距离小于 时，其接触表面的黏附应力为，于是接触表面的应力分布可重新表示为（）（），（）（）黏附接触模型（）涂层表面的应力分布（）（）图 涂层在导电压头作用下的力学模型 对于圆柱压头，当接触区 远小于压头半径 时，圆柱压头的外形可近似为如下抛物线，：（，），（）式中，表示发生在接触区域中心的最大压痕深度，表示刚性圆柱压头的半径，（，）（）黏附接触在 和 两点处的位移须满足如下条件：（，）（，）（）（）（）涂层接触表面 处的压痕 可表示为 （，）（，）（）涂层接触表面 处的电势 可表示为 （，）（，）（）奇异积分方程的建立考虑黏附时，式（）（）重新表示为（，）（）（）（）（，）（）（，）（）（）（）（，），（）（，）（）（）（，）（，）（）（）（，），（）（，）（）（）（，）应 用 数 学 和 力 学 年 第 卷 （，）（）（）（，），（）其中（，）（），（，）（），（，）（），（，）（），（，）（），（，）（）将式（）（）对 求导，可以得到如下无摩擦黏附接触问题的 奇异积分方程：（，）（）（）（，）（）（）（，）（）（）（），（）（，）（）（）（，）（）（）（，）（）（）（），（）（，）（）（）（，）（）（）（，）（）（）（），（）其中（，）（，），（，）（，），假设压头是一个绝缘体，那么法向电位移在接触表面为零，即（，）（，），从而可得到绝缘压头作用下，涂层和均匀压电基底半空间无摩擦黏附接触问题的奇异积分方程（，）（）（）（，）（）（）（），（）（，）（）（）（，）（）（）（）（）若为均匀压电半空间无摩擦黏附接触问题，则奇异积分方程为（，）（）（），（）（，）（）（），（）（，）（）（）（）第 期 韩立夫，等：导电压头作用下的功能梯度压电涂层二维黏附接触问题研究根据静力学平衡关系，考虑黏附时，涂层表面接触应力（）（），电荷分布（）与法向集中力、总电荷 满足下列关系：（），（）（）（）涂层表面接触应力（）在接触区边缘（）处是光滑的，在接触区内表面电势（，）是一个常数，则有（，）同时，根据文献，刚性圆柱导电压头的电荷（）可以分解为下面两个部分：（）（）（），（）其中，（）是由法向载荷 引起的表面电荷分布，且在接触边缘（）光滑，（）是由电势（，）引起的表面电荷分布，其在接触区边缘具有 奇异性，根据文献，结合式（）可得到（）（）（，）（）（）（，）（）（）（），（）（）（）（，）（）（）（，）（）（）（），（）（），（）（），（）（）（）（，），（）（），（）其中，是与法向荷载 关联的电荷，是与电势（，）关联的电荷 将式（）代入式（），得到求解位移差的表达式：（，），()（）（）（，）（，）（，）（，）（）（）（，）（，）（）积分区间的变换利用如下的变量代换，将式（）（）转化为如下形式：，（）（）（，）（）（）（，）（）（）（），（）（）（）（，）（）（）（，）应 用 数 学 和 力 学 年 第 卷 （）（）（），（）（），（）（），（）（）（）（，），（）（），（）其中（，）（）（），；，式（）变换为（）（）（，）（，）（，）（，）（）（）（，）（，）（）（）式（）变换为（）（）（，）（，）（，）（，）（）（）（，）（，）（）式（）变换为（）（）（，）（，）（，）（，）（）（）（，）（，）（）方程的离散化涂层表面接触应力（）和法向荷载引起的电荷分布（）可表示为，（）（），（）（）（）式（）、（）和（）为第一类 奇异积分方程，考虑 的方法，对其进行离散后可数值求解 首先将式（）（）进行离散，从而得到 （，）（）（，）（）()（）（），（）（，）（）（，）（）第 期 韩立夫，等：导电压头作用下的功能梯度压电涂层二维黏附接触问题研究 （）（）（），（）（），（）（），（）其中 （），（）（），（，）（）（），；，式（）（）所构成的方程组中共有 个方程，待求未知数（），（），（），（），（），（），为 个，其中 需要通过构建位移差方程来确定 根据文献，在求解时，离散点 取偶数，忽略式（）和（）中各自的第 个方程，这样得到的 个方程就可以求解除 之外的 个未知量 由于电荷（）具有 奇异性，可设，（）（）（）同样，使用 的方法对式（）和（）进行离散，从而得到 （，）（），（）（），（）其中 （）（），（），（，）（）（）在求解 值时，需要将方程式（）写成差值形式，然后利用二分法进行迭代，其差值形式如下：（）（，）（，）（）（，）（，）（，）（，）（）（，）（，）（）（）（）对式（）进行离散，从而得到 （）（）（）（）（）（）（）（）应 用 数 学 和 力 学 年 第 卷对式（）进行离散，从而得到 （）（）（）（）（）（）（），（）其中（）（）（），（）（）（），（）（）（）算例分析与讨论首先对本文所建立的二维无摩擦黏附接触模型进行退化，目的在于验证其科学有效 去掉式（）（）中的黏附项，然后采用 的计算方法和文献中相同的材料参数（基底由 压电陶瓷制成，具体参数见表）进行数值求解，计算选取 ，涂层厚度 ，圆柱压头半径 ，法向荷载 ，电荷 从而得到如图（）和（）所示的接触应力、电荷分布曲线 可以看出，本文的数值结果与文献给出的结果完全吻合，说明本文所建立的模型和计算方法是合理可靠的 表 压电陶瓷的材料参数 （）（）（）（）（）（）接触应力分布（）电荷分布（）（）图 本文模型计算结果与文献结果的比较 注 为了解释图中的颜色，读者可以参考本文的电子网页版本，后同 本文计算和分析刚性圆柱导电压头作用下 涂层的二维无摩擦黏附接触问题 计算时黏附功 取 ，圆柱压头半径 取 ，涂层厚度 取 ，数值计算结果如下文所示 图 给出了梯度参数 时，不同的黏附应力 下，比值与接触半径（图（）、接触半径与法向荷载（图（）、应力分布（）（图（）、电荷分布（）（图（）、法向荷载 与压痕深度（图（）、法向荷载 与电势（图（）的关系曲线 从图（）中可以看出：随着接触半径的逐渐增大，黏附区与接触区的比值 逐渐减小并趋近于，这与 等和 等，的结论相同；对于相同的接触半第 期 韩立夫，等：导电压头作用下的功能梯度压电涂层二维黏附接触问题研究径，随着黏附应力的增大，比值 减小，且在接触半径较小时，黏附应力对比值 的影响较大 从图（）中可以看出：当黏附应力较大时，随着接触半径的增大，法向荷载由某一值（由压电效应所产生）趋于一负的极大值（即临界拉脱力，主要由黏附效应引起），之后由拉力（负值）逐渐过渡到压力（正值）并逐渐增大；当黏附应力逐渐减小时，黏附效应引起的临界拉脱力逐渐变小，为 压电接触时，黏附效应完全消失，此时的临界拉脱力为零；当法向荷载为压力时，随着黏附应力的逐渐减小，产生相同的接触半径所需的法向荷载增大（）（）（）（）（）（）（）（）图 当 时，对各参数的影响 图（）表示接触半径 时，涂层接触表面的应力分布曲线，在涂层表面与压头接触区边缘处（）的应力值等于黏附应力值 由于黏附区大于接触区（），且在相同的接触区情况下，不同的黏附应力所对应的黏附区不同，因此本文只给出了接触区的应力 从图（）中可以看出：随着黏附应力的逐渐减小，涂层表面的拉应力在减小，压应力在增大，此结论与 等的研究结论一致 从图（）中应 用 数 学 和 力 学 年 第 卷可以看出：电荷分布几乎不受黏附应力变化的影响，即黏附效应对电荷分布的影响作用甚微 从图（）中可以看出：当黏附应力较大时，压痕深度由某一值（压电效应所引起）趋于一个极值（最大的负值，临界拉脱力所对应的压痕深度），之后由负值变为正值并逐渐增大；随着黏附应力的减小，临界拉脱力所对应的压痕深度（负值）逐渐减小；当法向荷载为压力时，随着黏附应力的逐渐增大，产生相同的压痕深度（正值）所需的法向荷载增大 从图（）中可以看出：当黏附应力较大时，电势由某一值（压电效应所引起）趋于一个极值（最大的负值，临界拉脱力所对应的电势），之后由负值变为正值并逐渐增大；随着黏附应力的减小，临界拉脱力所对应的电势（负值）逐渐减小 当法向荷载为压力时，随着黏附应力的逐渐增大，产生相同的电势（正值）所需的法向荷载增大 图 给出了黏附应力 时，不同的梯度参数 下，比值 与接触半径（图（）、接触半径 与法向荷载（图（）、应力分布（）（图（）、电荷分布（）（图（）、法向荷载 与压痕深度（图（）、法向荷载 与电势（图（）的关系曲线 从图（）中可以看出：随着接触半径的逐渐增大，比值 逐渐减小并趋近于；对于相同的接触半径，随着梯度参数的增大，比值 增大；且在接触半径较小时，梯度参数对比值 的影响较大 从图（）中可以看出：当梯度参数较小时，随着接触半径的增大，法向荷载由某一值（由压电效应所产生）趋于一个负值的极大值（临界拉脱力，主要由黏附效应引起），之后由拉力逐渐过渡到压力并逐渐增大；当梯度参数逐渐增大时，黏附效应引起的临界拉脱力逐渐减小；相同法向压力作用下，接触半径将随着梯度参数的增大而减小 从图（）中可以看出：随着梯度参数的逐渐增大，涂层表面的压应力增大 从图（）中可以看出：受压区（应力为正）的电荷分布随着梯度参数的增大而增大，而受拉区（应力为负）的电荷分布则随着梯度参数的增大而减小 从图（）中可以看出：当梯度参数较小时，压痕深度由某一值趋于一个极值（临界拉脱力所对应的压痕深度），之后由负值变为正值且逐渐增大；随着梯度参数的增大，临界拉脱力所对应的压痕深度（负值）逐渐减小（）（）（）（）（）（）第 期 韩立夫，等：导电压头作用下的功能梯度压电涂层二维黏附接触问题研究（）（）图 当 时，对各参数的影响 从图（）中可以看出：当梯度参数较小时，电势由负值变为正值并逐渐增大；随着梯度参数的增大，临界拉脱力所对应的电势（负值）逐渐减小 图 给出了黏附应力 时，不同的电荷 作用下，比值 与接触半径 （图（）、接触半径 与法向荷载（图（）、应力分布（）（图（）、电荷分布（）（图（）、法向荷载 与压痕深度（图（）、法向荷载 与电势（图（）的关系曲线 从图（）中可以看出：随着接触半径的逐渐增大，比值 逐渐减小并趋近于；对于相同的接触半径，正电荷比负电荷所对应的比值 要大；且在接触半径越小时影响越明显（）（）（）（）（）（）应 用 数 学 和 力 学 年 第 卷（）（）图 当 时，对各参数的影响 从图（）中可以看出：当受足量正电荷（）作用时，随着接触半径的增大，法向荷载由拉力逐渐过渡到压力并逐渐增大；当受足量负电荷（）作用时，法向荷载由某一值（由压电效应所产生）趋于一负的极大值（临界拉脱力，主要由黏附效应引起），之后由拉力逐渐过渡到压力并逐渐增大 从图（）和（）中可以看出：足量正电荷作用下涂层表面的应力和电荷分布比等电量负电荷作用下要大，且对电荷分布的影响更明显 从图（）中可以看出：相同法向压力作用下，带足量正电荷的压头所产生的压痕深度相对大，而带等电量负电荷的压头所产生的压痕深度相对小，说明所带的正电荷助长了压头的压入，而所带的负电荷则抑制了压头的压入；相同法向拉力作用下，结果反之 从图（）中可以看出：相同法向压力作用下，带足量正电荷的压头所产生的电势相对小，而带等电量负电荷的压头所产生的电势相对大；相同法向拉力作用下，结果亦反之 结 论综上所述，可得出如下结论：）对本文所建立的黏附接触模型进行退化求解，其结果与 等给出的结果完全吻合，且本文黏附接触模型的计算结果与 等和 等，的结果一致）当黏附应力逐渐减小时，临界拉脱力及其所对应的压痕深度和电势在减小，涂层表面的拉应力减小，压应力增大；当法向荷载为压力时，随着黏附应力的逐渐增大，产生相同的压痕深度和电势所需的法向荷载增大，而产生相同的接触半径所需的法向荷载减小）当梯度参数逐渐增大时，黏附效应引起的临界拉脱力逐渐减小，临界拉脱力所对应的压痕深度和电势减小，受拉区的电荷分布减小，而涂层表面的压应力和受压区的电荷分布增大；相同法向压力作用下，接触半径将随着梯度参数的增大而减小）足量正电荷作用下涂层表面的接触应力和电荷分布比等电量负电荷作用下要大，且对电荷分布的影响更明显；相同法向压力作用下，带足量正电荷的压头所产生的压痕深度相对大，所产生的电势相对小，而带等电量负电荷的压头所产生的压痕深度相对小，所产生的电势相对大；相同法向拉力作用下，反之）接触半径越小，黏附应力、梯度参数和压头带电量对比值 的影响越显著，说明接触区趋于微尺度时，黏附作用更加明显 因此，可以通过合理设计 涂层的梯度参数和黏附参数来改变压电材料表面的黏附接触行为，进而达到抑制压电器件接触损伤和破坏的目的 附 录 在平面应变状态下，横观各向同性 的本构关系为，第 期 韩立夫，等：导电压头作用下的功能梯度压电涂层二维黏附接触问题研究（，）（）（，）（）（，）（）（，），（）（，）（）（，）（）（，）（）（，），（）（，）（）（，）（，）（）（，），（）（，）（）（，）（，）（）（，），（）（，）（）（，）（）（，）（）（，），（）式中，其中 表示涂层，表示基底；（，），（，），（，），（，）和（，）分别为涂层或基底的应力分量和电位移分量；（，），（，），（，）分别为涂层或基底沿 和 方向的位移分量以及电势 忽略自重与体电荷，平衡方程和 方程可以表示为，（，）（，），（）（，）（，），（）（，），（，）（）将方程（）（）代入到方程（）（）中，可得到如下控制方程：（）()，（）（），（）()，（）其中，依次为（，），（，），（，）的简写形式 式（）（）对 进行 积分变换得到域内表达式如下，：（）()，（）（），（）()，（）其中，依次为（，），（，），（，）的简写形式，上标“”表示 积分变换域内变量，为 积分变换的参变量 将式（）（）联立求解，可得到涂层和基底在变换域内位移分量和电势的矩阵表达如下：（，）（，）（，）（，）（），（）其中，（）是需要求解的未知参数，涂层中的（，）（）（）（，），特征根 包括 个实根和 个复根；在均匀压电基底半空间（）中，当 时，积分变换域内应满足，所以基底中（，）（）（）（，），特征根 只包括实部大于零的 个根，上角标“”表示矩阵转置，求（，）的特征方程如下，：（，），（）其中（，）（）（）（）（）应 用 数 学 和 力 学 年 第 卷矩阵（，）中各元素若用 表示，那么 则代表矩阵（，）中第 行第 列所对应的元素，对于基底还需满足 的条件，矩阵（，）中的（）和（）则可表示为（），（）（）对式（）、（）、（）分别进行 积分变换，然后代入式（），可得到功能梯度压电涂层（）的位移分量、电势、应力分量和电荷分量在 积分变换域内的矩阵表达，：（，）（，）（，）（，）（，）（，）（，）（），（）其中，（）（）（）（）（）（）（），（，）（，）（，）（，）（，）（，）（，）（，）表示矩阵（，）中第 行和第 列的元素，展开形式如下：（，），（，）（），（，）（），（，）（）（）（），（，）（）（）（），（，）（）（）（）同理可得，均匀压电基底半空间的位移分量、电势、应力分量和电荷分量在 积分变换域内的矩阵表达为，（，）（，）（，）（，）（，）（，）（，）（），（）其中 ，；，（）（）（）（），（，）（，）（，）（，）（，）（，）（，），（，），（，）（），（，）（），（，）（）（），（，）（）（），（，）（）（）为了确定式（）和式（）中的未知参变量（）（，）和（）（，），在接触表面（）处，需满足如下边界条件：（，）（），（，）（），（，）（），（）其中，（）为 函数 在界面（）处，位移分量、应力分量、电位移分量和电势需满足下列连续性条件：（，）（，），（，）（，），（，）（，），（，）（，），（，）（，），（，）（，）（）在 积分变换域内，上述边界与连续性条件可表示成如下矩阵形式：（，）（），（）（，）（）（，）（）（）将式（）和（）联立求解，可得到（）和（）的表达式：（），（）（），（）其中 （，）（，），（，）将式（）代入式（）中，然后进行 逆变换，可得到 （，）（，）（，）（，），（）其中，（，）是一个 的矩阵，其表达式为（，）（，）从式（）中提取出涂层表面（）处位移分量和电势的表达式：第 期 韩立夫，等：导电压头作用下的功能梯度压电涂层二维黏附接触问题研究（，），（，），（，）（，），（）其中（，）（，），对矩阵（，）进行渐进性分析，得到（，），（，）（）于是，式（）可表示为（，）（，）（，）（，），（）其中 ，（）（）（）（）（）根据 公式和奇偶性以及下列变换关系：（）（），（，）（）（，），（）（，）（）（，），（，）（）（，），（）（），（）（）（）由式（）可导出表面（）处位移分量和电势的基本解为（，）（）（）（）（）（），（）（，）（）（）（）（），（）（，）（）（）（）（）（）式（）（）为 涂层和均匀压电基底半空间在法向、切向集中力，和正集中线电荷 作用下的二维接触问题的基本解 特别地，当不考虑摩擦时，式（）（）中包含切向集中力 的项将消失，可得到 涂层和均匀压电基底半空间在法向集中力 和正集中线电荷 作用下的二维无摩擦接触问题的基本解为（，）（）（）（）（），（）（，）（）（），（）（，）（）（）（）当 涂层和均匀压电基底为同一材料时，式（）（）中的梯度项将全部消失，可得均匀压电半空间接触问题的基本解为（，）（）（），（）应 用 数 学 和 力 学 年 第 卷（，）（），（）（，）（）（）参考文献（）：，：，：，（）：，（）：，：，（）：，（）：，（）：，（）：，（）：，（）：，（）：，：，（）：，（）：靳凡 先进功能材料的微尺度黏附接触力学研究 大连：大连理工大学，（：，（），（）：，：，（）：，（）：，（）：，：，第 期 韩立夫，等：导电压头作用下的功能梯度压电涂层二维黏附接触问题研究，（）：，：，（）：苏洁 功能梯度压电材料的摩擦接触分析 北京：北京交通大学，（：，（）刘兴伟，李星，汪文帅 一维六方压电准晶中正 边形孔边裂纹的反平面问题 应用数学和力学，（）：（，（）：（）马占洲，刘铁军 基于层合板模型的梯度压电涂层 问题研究 力学季刊，（）：（，（）：（）代文鑫，刘铁军 导电压头作用下的多层功能梯度压电材料涂层二维接触问题研究 应用数学和力学，（）：（，（）：（），（）：，：，（）：，（）：，（）：，（）：，：，：李培兴 考虑尺度效应的梯度纳米涂层的粘附与接触力学行为研究 呼和浩特：内蒙古工业大学，（：，（），：，：应 用 数 学 和 力 学 年 第 卷