压电材料的变分原理课件

压电材料的变分原理

和有限元分析方法赵寿根压电材料的变分原理

和有限元分析方法赵寿根1前言压电材料由于其机电耦合特性，受到使用者的欢迎。当将压电材料应用到结构中时，由于结构形式的多样性、边界条件的多样性和外界激励环境的复杂性，能得到解析解的情形很少，因而大多数情况需要采用数值的方法来分析结构。变分原理是进行数值计算的基础，因而研究压电材料的变分原理为建立压电材料的有限元模型和方程提供了依据。1前言压电材料由于其机电耦合特性，受到使用者的欢迎。当将压2基本方程压电材料具有力电耦合特性，根据连续弹性介质理论和电介质理论，基于线弹性、小变形假设，基本方程及条件如下。

(1)运动方程(2)电学方程

若不存在自由电荷则式子等于零2基本方程压电材料具有力电耦合特性，根据连续弹性介质理论和2.1本构方程(3)力学耦合方程或(4)电学耦合方程或2.1本构方程或(4)电学耦合方程或2.3几何方程(5)变形方程(6)电场方程2.3几何方程(6)电场方程2.3边界条件(7)力学边界条件(8)电学边界条件在Sσ上在Su上在Sq上在Sv上2.3边界条件(8)电学边界条件在Sσ上在Su上在Sq上3力电耦合系统的能量泛函(1)动能(2)应变能(3)电势能3力电耦合系统的能量泛函(1)动能(2)应变能(3)(4)外力功(5)外电荷功符合说明：Ve、Vp和V=Ve+

Vp分别为弹性材料体积、压电材料体积和总体积。(4)外力功(5)外电荷功符合说明：Ve、Vp和V=5系统广义泛函由Hamilton原理，系统广义泛函为：

将上面的本构方程、几何方程带入得到系统的能量泛函为：5系统广义泛函由Hamilton原理，系统广义泛函为：能量泛函写成矩阵形式有：上式即为分析压电耦合结构、建立各种位移形式的运动微分方程的变分形式方程。能量泛函写成矩阵形式有：上式即为分析压电耦合结构、建立各种位6有限元方法有限元分析，即有限元方法（冯康首次发现时称为基于变分原理的差分方法），是一种用于求解微分方程组或积分方程组数值解的数值技术.这一解法基于完全消除微分方程,即将微分方程转化为代数方程组(稳定情形);或将偏微分方程(组)改写为常微分方程(组)的逼近,这样可以用标准的数值技术(例如欧拉法,龙格-库塔方法等)求解.有限元法最初起源于土木工程和航空工程中的弹性和结构分析问题的研究.它的发展可以追溯到AlexanderHrennikoff(1941)和RichardCourant(1942)的工作.这些先驱者使用的方法具有很大的差异,但是他们具有共同的本质特征:利用网格离散化将一个连续区域转化为一族离散的子区域,通常叫做元.Hrennikoff的工作离散用类似于格子的网格离散区域;Courant的方法将区域分解为有限个三角形的子区域,用于求解来源于圆柱体转矩问题的二阶椭圆偏微分方程.Courant的贡献推动了有限元的发展,绘制了早期偏微分方程的研究结果.从有限元的基本方法派生出来的方法很多，则称为三维单元。如有限条法、边界元法、杂交元法、非协调元法和拟协调元法等，用以解决特殊的问题。6有限元方法有限元分析，即有限元方法（冯康首次发现时称为基6.1本构方程压电材料的线性本构方程为：用矩阵形式表示为：6.1本构方程压电材料的线性本构方程为：用矩阵形式表示为：6.2有限元列式(四面体单元为例)对于每一个单元，机械应变可以表示为：6.2有限元列式(四面体单元为例)对于每一个单元，机械应变又位移u、v、w可以用单元节点位移和形函数表示：又位移u、v、w可以用单元节点位移和形函数表示：从而可得：从而可得：上式用矩阵形式表示为：Bu为包含形函数微分的矩阵：上式用矩阵形式表示为：Bu为包含形函数微分的矩阵：对于每个单元的x、y、z方向的位移向量表示为：

同样对于每个单元，电场向量可以表示为为：对于每个单元的x、y、z方向的位移向量表示为：同样对于又，电势可以用单元节点电势和形函数表示：从而：又，电势可以用单元节点电势和形函数表示：从而：上式用矩阵形式表示为：Bφ为包含形函数微分的矩阵：上式用矩阵形式表示为：Bφ为包含形函数微分的矩阵：由压电材料的虚功原理：由前面有：由压电材料的虚功原理：由前面有：又：将上两式代入虚功原理有：又：将上两式代入虚功原理有：对上面的式子进行化简，就可得到单元的有限元方程：式中：对上面的式子进行化简，就可得到单元的有限元方程：式中：对单元有限元方程进行组装可得压电材料的整体有限元方程：对单元有限元方程进行组装可得压电材料的整体有限元方程：7

ANSYS中的压电分析压电分析只能用下列单元类型之一：(1)PLANE13,KEYOPT(1)=7，耦合场四边形板单元(2)SOLID5,KEYOPT(1)=0或3，耦合场六面体单元(3)SOLID98,KEYOPT(1)=0或3，耦合场四面体单元KEYOPT选项激活压电自由度：位移和电压。对于压电分析，必须激活位移和电压自由度，即对于SOLID98要选择Degrees

of

Freedom选项的值为UX,

UY,

UZ,

VOLT。对于PLANE13要选择Degrees

of

Freedom选项的值为UX,

UY,

VOLT。对于SOLID5和SOLID98，KEYOPT(1)=3仅激活压电选项。7ANSYS中的压电分析压电分析只能用下列单元类型之一：压电材料材料参数的输入介电常数是反映材料的介电性质，或极化性质的，通常用ε来表示。不同用途的压电陶瓷元器件对压电陶瓷的介电常数要求不同。例如，压电陶瓷扬声器等音频元件要求陶瓷的介电常数要大，而高频压电陶瓷元器件则要求材料的介电常数要小。压电陶瓷极化处理之前是各向同性的多晶体，这是沿1(x)、2(y)、3(z)方向的介电常数是相同的，即只有一个介电常数。经过极化处理以后，由于沿极化方向产生了剩余极化而成为各向异性的多晶体。此时，沿极化方向的介电性质就与其他两个方向的介电性质不同。设陶瓷的极化方向沿3方向，则有关系

ε11=ε22≠ε33

即经过极化后的压电陶瓷具有两个介电常数ε11和ε33压电材料材料参数的输入介电常数是反映材料的介电性质，或极化性介电系数矩阵(介电常数)用MP命令(MainMenu>Preprocessor>MaterialProps>MaterialModels>Electromagnetics>RelativePermittivity>Orthotropic)定义PERX、PERY和PERZ。这些常数分别表示的是介电系数矩阵

[ε]s（上标“s”表示常数值是在常应变条件下计得到的）的对角分量ε11,ε22,ε33。介电系数矩阵(介电常数)压电陶瓷具有压电性，即施加应力时能产生额外的电荷。其所产生的电荷与施加的应力成比例，对于压力和张力来说，其符号是相反的，用介质电位移D（单位面积的电荷）和应力T（单位面积所受的力）表示如下D=Q/A=dT

式中，d的单位为库仑/牛顿（C/N），这正是正压电效应。还有一个逆压电效应，既施加电场E时成比例地产生应变S，其所产生的应变为膨胀或为收缩取决于样品的极化方向。

S=dE

式中，d的单位为米/伏（m/v）。上面两式中的比例常数d称为压电应变常数。对于正和逆压电效应来讲，d在数值上是相同的。对于企图用来产生运动或振动(例如，声纳和超声换能器)的材料来说，希望具有大的压电应变常数d。另一个常用的压电常数是压电电压常数g，它表示内应力所产生的电场，或应变所产生的电位移的关系。常数g与常数d之间的关系如下：g=d/e

此外，还有不常用的压电应力常数e和压电劲度常数h；e把应力T和电场E联系起来，而h把应变S和电场E联系起来，既T=-eE

;

E=-hS压电陶瓷具有压电性，即施加应力时能产生额外的电荷。其所产生的压电矩阵可以定义[e]型（压电应力矩阵）或[d]型（压电应变矩阵）的压电矩阵。[e]型矩阵典型地与刚度矩阵[c]的各向异性弹性输入有关，而[d]矩阵与柔度矩阵[s]的输入相关。

[e]矩阵和[d]矩阵使用下列数据表输入：2D情况3D情况压电矩阵2D情况3D情况大多数已公布的压电材料的[e]矩阵