功能梯度压电板的样条有限点分析闻泽徐松波徐永伦合肥水泥研究.docVIP

功能梯度压电板的样条有限点分析 闻泽 徐松波 徐永伦 （合肥水泥研究设计院设计分院 邮编230051） 摘要：利用样条有限点法，对功能梯度压电梁板的力学性能进行了三维分析。用3 次B样条函数来拟合中性面的位移（电势）函数，通过引入位移模式，得到结构中任意一点的位移，将力学和电学所有基本量用这个假定的函数表示。对功能梯度压电板的静力问题和屈曲问题进行了公式推导。通过算例和已有的相关计算结果作以比较，证明本文方法的合理性和准确性；并且通过算例分析了不同梯度模式变化下，功能梯度压电板的力学量和电学量沿梯度方向的变化情况。 关键词：功能梯度压电材料 物性参数 样条有限点法 屈曲分析 引言 压电材料在变形时产生电场，在受电场作用时产生变形，因此根据这种内在的机电耦合性质，使得压电材料在诸多领域有着广泛应用。正压电效应被广泛运用于智能结构的感知元件，以实现对结构的形状控制、振动和噪声控制及结构的损伤检测等；逆压电效应被广泛的应用于智能结构的驱动器。然而传统的压电材料在两种不同的界面交接处由于材料性质的突然变化，容易产生力学性能的突变，这对压电材料的受力是很不利的。现在已经开始了对功能梯度压电材料的力学性能进行研究。功能梯度压电材料的组成和性能在空间上的连续变化，减少了传统压电材料在分界面上了应力集中，大大提高了压电材料的可靠性和寿命。功能梯度压电结构中，物性参数随坐标变化而变化，弹性场与电场并存，基本方程包含力场和电场，求解十分困难。近年来，为了实际应用的需要，关于功能梯度压电材料的理论在国内外都取得了广泛的成果；在精确解，简化理论解，有限元解，边界元解等多种分析方法上都取得了不少的进展。Chen和Ding[1]基于状态方程分析了四边简支压电功能梯度板的弯曲问题。 样条有限点法是一种半解析半数值方法, 充分发挥了样条函数的紧凑、精度高的优点，克服了有限元方法中由于场函数的局部化近似所引起的误差，场函数及其梯度在整个求解域内是连续的。功能梯度材料的显著特点是物性参数沿着梯度方向变化，而样条有限点法可以将变化的物性参数直接在梯度方向积分求得整体的物性参数矩阵。在求解的过程中力学概念非常清晰，也易于编制程序。本文将利用样条有限点法对功能梯度压电材料的力学性能进行一些研究。 1 力学模型 不同的梯度变化模型决定了组成和性能的具体变化方式，不同的变化方式又决定了材料内部应力的不同。因此，对材料参数按不同梯度模型的变化的研究是十分有意义，也为计算功能梯度压电材料的力学性能提供了可能。本文所选取的功能梯度计算模型如图1所示，材料的物性参数功能梯度变化是沿z轴方向变化的。假定材料的物性参数从底面到顶面按照种函数变化，材料的功能梯度变化只是z的函数。如图1所示取坐标原点位于板的底面，z为垂直中面的法向坐标。当顶面和底面物性参数已知情况下，板的材料遵循如下梯度分布规律[14]]： 其中可代表、、，分别为弹性刚度系数、压电系数和介电系数， 为功能梯度压电材料板壳底部(或内部)的物参，为材料的体积份数分布函数。 一般情况下，三种函数形式： 1．指数函数形式： 2．幂函数形式： 3．正弦函数形式： 当处材料的物参，已知，通过计算，可以确定形状因子、、的大小.。 2.基本方程: 2.1本构方程: 其中 E,B分别为应力向量,应变向量,电场强度向量及电位移向量;D,e,分别为弹性矩阵,压电应力常数矩阵,介电常数矩阵. 2.2位移模式 在文献[15]中，对水平位移偏离直法线部分的位移作了沿截面几何相似的假定，使得这个偏离位移可以进行分离变量，再引入Bnacob假定，将厚板的位移模式描述如下： 其中：为中性面的位移，为中面剪切角，表征了截面的非线性变形，h为板的厚度。本文将采用（4）的形式作为位移模式进行求解。 电势是标量，也是一个相对量，为了研究的方便，可以选取板的上表面的电势为0，中面的电势为，在梯度方向上线性分布，可以得到任意一点的电势为： = 2.3建立能量泛函： 利用最小势能