《压电材料的变分原理ppt课件.ppt》由会员分享，可在线阅读，更多相关《压电材料的变分原理ppt课件.ppt（30页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、School of Aeronautical Science and Engineering, BUAASchool of Aeronautical Science and Engineering, BUAA赵寿根赵寿根School of Aeronautical Science and Engineering, BUAASchool of Aeronautical Science and Engineering, BUAAQ 压电材料由于其机电耦合特性，受到使用者的欢压电材料由于其机电耦合特性，受到使用者的欢迎。当将压电材料应用到结构中时，由于结构形迎。当将压电材料应用到结构中时，由于结构形
2、式的多样性、边界条件的多样性和外界激励环境式的多样性、边界条件的多样性和外界激励环境的复杂性，能得到解析解的情形很少，因而大多的复杂性，能得到解析解的情形很少，因而大多数情况需要采用数值的方法来分析结构。数情况需要采用数值的方法来分析结构。Q 变分原理是进行数值计算的基础，因而研究压电变分原理是进行数值计算的基础，因而研究压电材料的变分原理为建立压电材料的有限元模型和材料的变分原理为建立压电材料的有限元模型和方程提供了依据。方程提供了依据。School of Aeronautical Science and Engineering, BUAASchool of Aeronautical Sci
3、ence and Engineering, BUAAQ 压电材料具有力电耦合特性，根据连续弹性介质压电材料具有力电耦合特性，根据连续弹性介质理论和电介质理论，基于线弹性、小变形假设，理论和电介质理论，基于线弹性、小变形假设，基本方程及条件如下。基本方程及条件如下。 (1) 运动方程运动方程iij ,ijuf (2) 电学方程电学方程fi , iD 若不存在自由电荷则式子等于零School of Aeronautical Science and Engineering, BUAASchool of Aeronautical Science and Engineering, BUAAQ2.1 本构
4、方程本构方程(3) 力学耦合方程力学耦合方程kijkklijklijEdsEkijkklijklijEeCEE或(4) 电学耦合方程jEdDijklikli或jEeDijklikliSchool of Aeronautical Science and Engineering, BUAASchool of Aeronautical Science and Engineering, BUAAQ2.3 几何方程几何方程(5) 变形方程变形方程(6) 电场方程电场方程j ,iji ,ijijuu21 i ,igradESchool of Aeronautical Science and Enginee
5、ring, BUAASchool of Aeronautical Science and Engineering, BUAAQ2.3 边界条件边界条件(7) 力学边界条件力学边界条件(8) 电学边界条件电学边界条件ijijTn 在S上iiuu 在Su上qnDii在Sq上V在Sv上School of Aeronautical Science and Engineering, BUAASchool of Aeronautical Science and Engineering, BUAAQ (1) 动能动能dVuu21dVuu21TiipviievpeQ (2) 应变能应变能dV21dV21Uii
6、viivpeQ (3) 电势能电势能dVDE21UmmvppSchool of Aeronautical Science and Engineering, BUAASchool of Aeronautical Science and Engineering, BUAAQ (4) 外力功外力功viiviTdSTufdVuWeQ (5) 外电荷功外电荷功qsqqdSW符合说明：符合说明： Ve、Vp和和V= Ve+ Vp分别为弹性材料体积、分别为弹性材料体积、压电材料体积和总体积。压电材料体积和总体积。School of Aeronautical Science and Engineering,
7、BUAASchool of Aeronautical Science and Engineering, BUAAQ 由由Hamilton原理，系统广义泛函为：原理，系统广义泛函为：t1t2qTp0dtWWUUT 将上面的本构方程、几何方程带入得到系统的能将上面的本构方程、几何方程带入得到系统的能量泛函为：量泛函为： ssiiVViiijijViiqp0dsqdSuTdvdEDdvdvuu School of Aeronautical Science and Engineering, BUAASchool of Aeronautical Science and Engineering, BUAA
8、能量泛函写成矩阵形式有：能量泛函写成矩阵形式有： qpppssTvvTTVVTTTVTdSdqdsTudVEEdVeEdvEedvCdvuu 上式即为分析压电耦合结构、建立各种位移形式的运动上式即为分析压电耦合结构、建立各种位移形式的运动微分方程的变分形式方程。微分方程的变分形式方程。School of Aeronautical Science and Engineering, BUAASchool of Aeronautical Science and Engineering, BUAAQ有限元分析，即有限元方法（冯康首次发现时称为基于变分原理的差分方法），有限元分析，即有限元方法（冯康首次
9、发现时称为基于变分原理的差分方法），是一种用于求解微分方程组或积分方程组数值解的数值技术是一种用于求解微分方程组或积分方程组数值解的数值技术. 这一解法基于完全这一解法基于完全消除微分方程消除微分方程, 即将微分方程转化为代数方程组即将微分方程转化为代数方程组(稳定情形稳定情形); 或将偏微分方程或将偏微分方程(组组)改写为常微分方程改写为常微分方程(组组)的逼近的逼近, 这样可以用标准的数值技术这样可以用标准的数值技术(例如欧拉法例如欧拉法,龙格龙格-库库塔方法等塔方法等)求解求解.Q有限元法最初起源于土木工程和航空工程中的弹性和结构分析问题的研究有限元法最初起源于土木工程和航空工程中的弹性
10、和结构分析问题的研究. 它的它的发展可以追溯到发展可以追溯到Alexander Hrennikoff(1941)和和Richard Courant(1942)的工作的工作. 这这些先驱者使用的方法具有很大的差异些先驱者使用的方法具有很大的差异, 但是他们具有共同的本质特征但是他们具有共同的本质特征: 利用利用网格网格离散化将一个连续区域转化为一族离散的子区域离散化将一个连续区域转化为一族离散的子区域, 通常叫做元通常叫做元.Hrennikoff 的工作的工作离散用类似于格子的网格离散区域离散用类似于格子的网格离散区域; Courant 的方法将区域分解为有限个三角形的方法将区域分解为有限个三角
11、形的子区域的子区域, 用于求解来源于圆柱体转矩问题的二阶椭圆偏微分方程用于求解来源于圆柱体转矩问题的二阶椭圆偏微分方程. Courant 的贡的贡献推动了有限元的发展献推动了有限元的发展, 绘制了早期偏微分方程的研究结果绘制了早期偏微分方程的研究结果.Q从有限元的基本方法派生出来的方法很多，则称为三维单元。如从有限元的基本方法派生出来的方法很多，则称为三维单元。如有限条法、边有限条法、边界元法、杂交元法、非协调元法和拟协调元法界元法、杂交元法、非协调元法和拟协调元法等，用以解决特殊的问题。等，用以解决特殊的问题。School of Aeronautical Science and Engine
12、ering, BUAASchool of Aeronautical Science and Engineering, BUAA 压电材料的线性本构方程为：压电材料的线性本构方程为：) 3 , 2 , 1()6 , 2 , 1(lESeDiEeSCjljklkljEijkEiki用矩阵形式表示为：用矩阵形式表示为： ESeDEeSCSchool of Aeronautical Science and Engineering, BUAASchool of Aeronautical Science and Engineering, BUAAQ 对于每一个单元，机械应变可以表示为：对于每一个单元，机械
13、应变可以表示为：xvyuxwzuywzvzwyvxuSxyxzyzzyxeSchool of Aeronautical Science and Engineering, BUAASchool of Aeronautical Science and Engineering, BUAAQ 又位移又位移u、v、w可以用单元节点位移和形函数表示：可以用单元节点位移和形函数表示：niiiuNu1niiivNv1niiiwNw1School of Aeronautical Science and Engineering, BUAASchool of Aeronautical Science and Eng
14、ineering, BUAAQ 从而可得：从而可得： xvNyuNxwNzuNywNzvNzwNyvNxuNSniiiniiiniiiniiiniiiniiiniiiniiiniiixyxzyzzyxe111111111School of Aeronautical Science and Engineering, BUAASchool of Aeronautical Science and Engineering, BUAA上式用矩阵形式表示为：上式用矩阵形式表示为： eueuBSBu为包含形函数微分的矩阵：为包含形函数微分的矩阵：000000000 xyxzyzzyxuNNNNNNNNNBS
15、chool of Aeronautical Science and Engineering, BUAASchool of Aeronautical Science and Engineering, BUAA对于每个单元的对于每个单元的x、y、z方向的位移向量表示为：方向的位移向量表示为： ieiiuuv 同样对于每个单元，电场向量可以表示为为：同样对于每个单元，电场向量可以表示为为： zyxEEEEzyxeSchool of Aeronautical Science and Engineering, BUAASchool of Aeronautical Science and Engineer
16、ing, BUAA又，电势可以用单元节点电势和形函数表示：又，电势可以用单元节点电势和形函数表示：niiiN1从而：从而： zNyNxNEEEEniiiniiiniiizyxe111School of Aeronautical Science and Engineering, BUAASchool of Aeronautical Science and Engineering, BUAA上式用矩阵形式表示为：上式用矩阵形式表示为： eeBEB为包含形函数微分的矩阵：为包含形函数微分的矩阵：zyxNNNBSchool of Aeronautical Science and Engineering
17、, BUAASchool of Aeronautical Science and Engineering, BUAAQ 由压电材料的虚功原理：由压电材料的虚功原理：Q 由前面有：由前面有： eeueeeueBuBeDBeuBC eTeeVTeeTeeVTeQudvDEFudvSSchool of Aeronautical Science and Engineering, BUAASchool of Aeronautical Science and Engineering, BUAA又：又：Q 将上两式代入虚功原理有：将上两式代入虚功原理有： eTeVeTTeuTTeeTeVeTTueuTuTe
18、QdvBBuBeBFudvBeBuBCBuTTeTeTuTeTeBEBuSSchool of Aeronautical Science and Engineering, BUAASchool of Aeronautical Science and Engineering, BUAA对上面的式子进行化简，就可得到单元的有限元方程：对上面的式子进行化简，就可得到单元的有限元方程：式中： eeeeeeueueuuFuKKKKQ dvBCBKuTvueuu dvBeBKKTTvuTeueu dvBBKTveSchool of Aeronautical Science and Engineering,
19、BUAASchool of Aeronautical Science and Engineering, BUAA对单元有限元方程进行组装可得压电材料的整体有限元对单元有限元方程进行组装可得压电材料的整体有限元方程：方程： QFuKKKKuuuuSchool of Aeronautical Science and Engineering, BUAASchool of Aeronautical Science and Engineering, BUAAQ 压电分析只能用下列单元类型之一：压电分析只能用下列单元类型之一：(1) PLANE13, KEYOPT(1)=7，耦合场四边形板单元，耦合场四边
20、形板单元(2) SOLID5, KEYOPT(1)=0或或3，耦合场六面体单元，耦合场六面体单元(3) SOLID98, KEYOPT(1)=0或或3，耦合场四面体单元，耦合场四面体单元Q KEYOPT选项激活压电自由度：位移和电压。对于压电分选项激活压电自由度：位移和电压。对于压电分析，必须激活位移和电压自由度，即对于析，必须激活位移和电压自由度，即对于SOLID98要选择要选择Degrees of Freedom选项的值为选项的值为UX, UY, UZ, VOLT。对于。对于PLANE13要选择要选择Degrees of Freedom选项的值为选项的值为UX, UY, VOLT。对于。对
21、于SOLID5和和SOLID98，KEYOPT(1)=3仅激活压电选项。仅激活压电选项。School of Aeronautical Science and Engineering, BUAASchool of Aeronautical Science and Engineering, BUAAQ 介电常数是反映材料的介电性质，或极化性质的，通常用介电常数是反映材料的介电性质，或极化性质的，通常用来表来表示。不同用途的压电陶瓷元器件对压电陶瓷的介电常数要求不示。不同用途的压电陶瓷元器件对压电陶瓷的介电常数要求不同。例如，压电陶瓷扬声器等音频元件要求陶瓷的介电常数要同。例如，压电陶瓷扬声器等音频
22、元件要求陶瓷的介电常数要大，而高频压电陶瓷元器件则要求材料的介电常数要小。大，而高频压电陶瓷元器件则要求材料的介电常数要小。Q 压电陶瓷极化处理之前是各向同性的多晶体，这是沿压电陶瓷极化处理之前是各向同性的多晶体，这是沿1(x)、2(y)、3(z)方向的介电常方向的介电常 数是相同的，即只有一个介电常数。经过数是相同的，即只有一个介电常数。经过极化处理以后，由于沿极化方向产生了剩余极化而成为各向异极化处理以后，由于沿极化方向产生了剩余极化而成为各向异性的多晶体。此时，沿极化方向的介电性质就与其他两个方向性的多晶体。此时，沿极化方向的介电性质就与其他两个方向 的介电性质不同。设陶瓷的极化方向沿的
23、介电性质不同。设陶瓷的极化方向沿3方向，则有关系方向，则有关系 11=2233 即经过极化后的压电陶瓷具有两个介电常数即经过极化后的压电陶瓷具有两个介电常数11和和33School of Aeronautical Science and Engineering, BUAASchool of Aeronautical Science and Engineering, BUAAQ 介电系数矩阵介电系数矩阵(介电常数介电常数) 用用MP命令命令(Main MenuPreprocessorMaterial PropsMaterial ModelsElectromagnetics Relative Pe
24、rmittivityOrthotropic)定义定义PERX、PERY和和PERZ。这些常数分别表示的是介电系数。这些常数分别表示的是介电系数矩阵矩阵 s（上标（上标“s”表示常数值是在常应变条件下表示常数值是在常应变条件下计得到的）的对角分量计得到的）的对角分量11,22,33。Q 压电陶瓷具有压电性，即施加应力时能产生额外的电荷。其所产生压电陶瓷具有压电性，即施加应力时能产生额外的电荷。其所产生的电荷与施加的应力成比例，对于压力和张力来说，其符号是相反的电荷与施加的应力成比例，对于压力和张力来说，其符号是相反的，用介质电位移的，用介质电位移D（单位面积的电荷）和应力（单位面积的电荷）和应力
25、T（单位面积所受（单位面积所受的力）表示如下的力）表示如下D=Q/A=dT 式中，式中，d的单位为库仑的单位为库仑/牛顿（牛顿（C/N），这正是正压电效应。还有一个逆压电效应，既施加电场，这正是正压电效应。还有一个逆压电效应，既施加电场E时成比时成比例地产生应变例地产生应变S，其所产生的应变为膨胀或为收缩取决于样品的极其所产生的应变为膨胀或为收缩取决于样品的极化方向化方向。 S=dE 式中，式中，d的单位为米的单位为米/伏（伏（m/v）。上面两式中的比）。上面两式中的比例常数例常数d称为压电应变常数。对于正和逆压电效应来讲，称为压电应变常数。对于正和逆压电效应来讲，d在数值在数值上是相同的。上
26、是相同的。Q 对于企图用来产生运动或振动对于企图用来产生运动或振动(例如，声纳和超声换能器例如，声纳和超声换能器)的材料来的材料来说，希望具有大的压电应变常数说，希望具有大的压电应变常数d。另一个常用的压电常数是压电。另一个常用的压电常数是压电电压常数电压常数g，它表示内应力所产生的电场，或应变所产生的电位移，它表示内应力所产生的电场，或应变所产生的电位移的关系。常数的关系。常数g与常数与常数d之间的关系如下：之间的关系如下：g=d/e 此外，还有不常用此外，还有不常用的压电应力常数的压电应力常数e和压电劲度常数和压电劲度常数h；e把应力把应力T和电场和电场E联系起来，联系起来，而而h把应变把
27、应变S和电场和电场E联系起来，既联系起来，既T=-eE ; E=-hSSchool of Aeronautical Science and Engineering, BUAASchool of Aeronautical Science and Engineering, BUAAQ 压电矩阵压电矩阵 可以定义可以定义e型（压电应力矩阵）或型（压电应力矩阵）或d型（压电应变矩阵）的型（压电应变矩阵）的压电矩阵。压电矩阵。e型矩阵典型地与刚度矩阵型矩阵典型地与刚度矩阵c的各向异性弹性输的各向异性弹性输入有关，而入有关，而d矩阵与柔度矩阵矩阵与柔度矩阵s的输入相关。的输入相关。 e矩阵和矩阵和d矩矩阵
28、使用下列数据表输入：阵使用下列数据表输入： 111213212223313233414243515253616263eeexeeeyeeezeeeexyeeeyzeeezx 1112212231324142eexeeyeeezeexy2D情况情况3D情况情况School of Aeronautical Science and Engineering, BUAASchool of Aeronautical Science and Engineering, BUAAQ 大多数已公布的压电材料的大多数已公布的压电材料的e矩阵数据都是基于矩阵数据都是基于IEEE标准标准按照按照x,y,z,yz,xz,
29、xy的顺序，而的顺序，而ANSYS的输入数据是按照的输入数据是按照x,y,z,xy,yz,xz的顺序。也就是说，输入该参数时必须通过改的顺序。也就是说，输入该参数时必须通过改变剪切项的行数据以转换到变剪切项的行数据以转换到ANSYS数据格式。数据格式。 111213212223313233616263414243515253eeexeeeyeeezANSYSeeeexyeeeyzeeezx中Ansys中中e矩阵的输入矩阵的输入School of Aeronautical Science and Engineering, BUAASchool of Aeronautical Science and Engineering, BUAAQ 刚度矩阵和柔度矩阵的输入与其它材料相同。刚度矩阵和柔度矩阵的输入与其它材料相同。