用ANSYS软件分析压电陶瓷的振动状态
用ANSYS软件分析压电陶瓷的振动状态
第30卷第1期2008年3月湖北大学学报(自然科学版)Journal of Hubei University (Natural Science )Vol.30 No.1 Mar.,2008 收稿日期:2006206202作者简介:雷辉(19812 ),男,硕士生文章编号:100022375(2008)0120029205用ANSYS 软件分析压电陶瓷的振动状态雷辉,周双娥(湖北大学数学与计算机学院,湖北武汉430062)摘 要:近年来,压电陶瓷的应用日趋广泛.但是由于压电陶瓷片的边界条件和应力状况比较复杂,利用传统实验手段对其研究不仅耗时费力,而且其结果具有很强的局部性,因此利用大型通用仿真软件ANSYS 8.0来进行计算机仿真.通过对压电陶瓷片中的耦合效应进行计算机模拟分析,得出压电陶瓷的振动状态图.实验结果表明ANSYS 8.0在处理压电耦合场这方面有很强的处理能力,大大简化了建模和计算,强大的后处理功能更是让研究者能够很直观地获得数据结果和模拟图像. 关键词:仿真;压电陶瓷;振动状态 中图分类号:TP302 文献标志码:A1 引言计算机仿真技术是以多种学科和理论为基础,以计算机及其相应构件为工具,通过虚拟试验的方法来分析和解决问题的一门综合性技术[1].近年来,压电陶瓷的应用日趋广泛,而在实际应用中,特别是将压电陶瓷技术应用于混凝土结构的监测中,由于压电陶瓷片的边界条件和应力状况比较复杂,利用传统实验手段对其研究不仅耗时费力,而且其结果具有很强的局部性[2].因此利用计算机仿真技术对压电陶瓷进行研究具有较好的理论与实际意义.本文中利用大型通用有限元分析软件ANS YS 8.0,对压电陶瓷片中的耦合效应进行模拟分析,并得出其模态和谐振态,实验表明ANS YS 8.0能很好地解决压电陶瓷片的压电耦合问题.图1 处理器模型2 ANSYS 仿真原理ANS YS 软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件,由世界上最大的有限元分析软件公司之一的美国ANS YS 开发,它能与多数CAD 软件接口,实现数据的共享和交换.20世纪90年代该软件开始在我国的机械制造、航空航天、汽车交通、铁道、石油化工、能源等领域得到应用,为各领域中产品设计、科学研究作出了很大的贡献[3].ANS YS 软件使用统一的集中式数据库来存储所有模型数据和求解结果(见图1)[4].模型数据(包括实体模型和有限元模型、材料等)通过前处理器写入数据库;载荷和求解结果通过求解器写入数据库;后处理结果通过后处理器写入数据库.3 处理过程3.1 定义材料参数 材料参数包括定义单元类型,这里选取了solid226,并在它的option 选项里选择压 湖北大学学报(自然科学版)第30卷30电这个选项.然后定义压电陶瓷的密度、介电常数、刚度系数和压电常数,后两个参数是用矩阵的形式来表示的.为了方便后期的网格划分处理,添加了一个单元类型mesh200,它主要用来进行面划分,下一节将详细介绍.定义材料参数的部分代码及注释如下所示:ET,1,SOL ID226,101 !定义solid226单元类型KEYOP T,1,1,1001!在solid226选项中选择压电选项ET,2,M ESH200!3!定义mesh200单元类型KEYOP T,2,1,7KEYOP T,2,2,0MPTEMP,,,,,,,,MPTEMP,1,0MPDA TA,DENS,1,7600!定义压电陶瓷密度TB,AN EL,1,1,21,0!定义压电陶瓷的刚度系数TB TEMP,0TBDA TA,,1.32e11,7.1e10,7.3e10,,,TBDA TA,,1.32e11,7.3e10,,,,1.15e11TBDA TA,,,,,3.0e10,,TBDA TA,,2.6e10,,2.6e10,,,MPTEMP,,,,,,,,!定义压电陶瓷的介电常数MPTEMP,1,0MPDA TA,PERX,1,,7.124e-9MPDA TA,PER Y,1,,7.124e-9MPDA TA,PERZ,1,,5.841e-9TB,PIEZ,1,,,0!定义压电陶瓷的压电系数TBMODIF,1,1,TBMODIF,1,2,TBMODIF,1,3,-4.1TBMODIF,2,1,TBMODIF,2,2,TBMODIF,2,3,-4.1TBMODIF,3,1,TBMODIF,3,2,TBMODIF,3,3,14.1TBMODIF,4,1,TBMODIF,4,2,TBMODIF,4,3,TBMODIF,5,1,TBMODIF,5,2,10.5TBMODIF,5,3,TBMODIF,6,1,10.5TBMODIF,6,2,TBMODIF,6,3,3.2 建立模型及网格划分 首先新建一个长为0.005m,宽为0.001m的矩形.然后选取单元类型为mesh200,对该矩形进行面划分,其中长和宽分别划分8份和10份,结果如图2所示.然后用操作选项中的Ext rude命令将其扩展成一个已划分好的圆盘体,其中要在Extrude选项中将单元类型选择为solid226,并在要划分的数目中填入10.接着用Extrude命令将其扩展成一个圆盘,其扩展结果如图3所示.第1期雷辉等:用ANSYS 软件分析压电陶瓷的振动状态31 图2 面划分图 图3 体划分图3.3 添加约束条件和负载 添加的约束条件是在上、下表面的圆心处分别添加位移约束,使其只能沿纵向方向移动.具体代码及注释如下:nsel ,s ,loc ,y ,0nsel ,r ,loc ,z ,0nsel ,r ,loc ,x ,0!选取下表面圆心上的点d ,all ,ux ,0!使其不能沿x 方向运动d ,all ,uz ,0!使其不能沿y 方向运动nsel ,s ,loc ,y ,0.001!选取上表面圆心上的点nsel ,r ,loc ,z ,0!同理nsel ,r ,loc ,x ,0d ,all ,ux ,0d ,all ,uz ,0给圆盘添加的电压负载是在圆盘的上下表面的两个电极上加上耦合电压,其电压值分别为5伏和0伏.这样使得圆盘上下电势差为5伏.其代码及注释如下:nsel ,s ,loc ,y ,0!选取下表面所有节点cp ,1,volt ,all !为下表面添加耦合电压3get ,n_bot ,node ,0,num ,min !定义下表面的电极节点nsel ,s ,loc ,y ,0.001!选取上表面所有节点cp ,2,volt ,all !为上表面添加耦合电压3get ,n_top ,node ,0,num ,min !定义上表面的电极节点nsel ,all !选取所有节点d ,n_bot ,volt ,0!下表面加0伏电压d ,n_top ,volt ,5!上表面加5伏电压3.4 静态和模态下的处理 首先选择处理状态为静态,并在此状态下得出其静态电容,代码及注释如下:/SOL !进入处理环节AN T YPE ,0!选择静态处理/STA TUS ,SOL U SOL V E !求解3get ,cs ,node ,n_top ,rf ,chrg !得到上电极电量值fini !处理结束3SET ,cs ,abs (cs )/5!得到静态电容接着进行模态分析,设定它有20个子步,频率范围为02200000Hz ,并且将上表面短路,代码及注释如下:/SOL U !进入处理环节AN T YPE ,2!选择模态处理MODOPT ,L ANB ,20!定义处理方法和子步数EQSL V ,SPAR !选择处理器MXPAND ,20,,,1!要显示的子步数L UMPM ,0PSTRES ,0MODOPT ,L ANB ,20,0,2000000,,OFF !定义频率范围32 湖北大学学报(自然科学版)第30卷d ,n_top ,volt ,0!上表面短路nsel ,all!选择所有节点/STA TUS ,SOL USOL V E !求解3.5 后处理 后处理是指检查并分析求解的结果的相关操作.这是分析中最重要的环节之一,因为任何分析的最终目的都是为了研究作用在模型上的的载荷是如何影响设计的[5].检查分析结果可使用两个后处理器:POST1(通用后处理器)和POST26(时间历程后处理器).POST1允许检查整个模型在某一载荷步和子步(或对某一特定时间点或频率)的结果.POST26可以检查模型的指定节点的某一结果项相对于时间、频率或其它结果项的变化.在求解时,ANS YS 将结果写入结果文件,进行后处理时,结果文件必须存在且可用.结果文件名的后缀取决于分析类型,对于结构分析的结果文件的后缀为RST [6].本次实验只用到了POST1.在后处理中利用一个循环语句可以得出不同频率下的动态电容、动态电感、电量以及阻抗值.其代码和注释如下:/POST13SET ,nmodes ,20!定义nmodes =203dim ,C ,array ,nmodes!定义电容数组3dim ,L ,array ,nmodes!定义电感数组3SET ,PI2,233.14159!定义PI2=233.14159set ,first!设定第一个子步/com ,3do ,i ,1,nmodes !定义一个循环,从1到20步3get ,Fi ,mode ,i ,freq !得到该步的频率值3get ,Qi ,node ,n_top ,rf ,chrg !得到该步的电量值3SET ,Omi ,Pi23Fi !将线速度转化为角速度3SET ,C (i ),(Qi/Omi )332!计算相关的动态电容3SET ,L (i ),1/(Omi 3323C (i ))!计算相关的动态电感/com ,Mode %i %!在输出窗口中显示第几步/com ,Resonant f requency F =%Fi %Hz !在输出窗口中显示频率值/com ,Dynamic capacitance C =%C (i )%F !在输出窗口中显示动态电容值/com ,Dynamic inductance L =%L (i )%H !在输出窗口中显示动态电感值/com ,charge Q =%Qi %C!在输出窗口中显示电量值/com ,Impedance R =%5/(233.141593Fi 3Qi )%!在输出窗口中显示阻抗值/com ,set ,next !进入下一个子步3enddo!循环结束FINISH同时可以在主菜单的Animation 选项下看一下压电圆盘在各个不同频率下的振动状态,这里截取了几个振动图,如图4所示:(a )频率为29728Hz (b )频率为65741Hz (c )频率为53476Hz图4 振动状态图第1期雷辉等:用ANSYS软件分析压电陶瓷的振动状态33 4 结束语本文对压电陶瓷圆盘添加了约束和对称负载,解决了压电陶瓷压电场与结构场的耦合问题,并且最终得到了压电圆盘在静态下的电容值,以及它在模态下的动态电容、动态电感、电量和阻抗,而且还可以观察圆盘在各个不同频率下的振动状态,以便于今后对材料性能的研究以及对材料的改进.同时,不难发现Ansys8.0在处理压电耦合场这方面有很强的处理能力,像Solid226就是专门针对压电分析而定义的,而且Ansys8.0在其材料库中建立了相关的压电材料,因此大大简化了有限元的建模和计算,强大的后处理功能更是让研究者能够很直观的获得数据结果和模拟图像.参考文献:[1]林书玉,张福成.压电陶瓷圆片振子的多模耦合振动[J].电子学报,1994,12:43249.[2]姜德义,郑拯宇.压电陶瓷片耦合振动模态的ANSYS模拟分析[J].传感技术学报,2003,12:9216.[3]陈大任.压电陶瓷微位移驱动器概述[J].电子元件与材料,1994,2:33240.[4]邵蕴秋.ANSYS8.0有限元分析实例导航[M].北京:中国铁道出版社,2004.[5]刘涛.精通ANSYS[M].北京:清华大学出版社,2002.[6]任重.ANSYS实用分析教程[M].北京:北京大学出版社,2003.Using ANSYS to analyze the vibration state of piezoelectric ceramicL EI Hui,ZHOU Shuang2e(School of Mathematics and Computer Science,Hubei University,Wuhan430062,China) Abstract:The application of piezoelect ric ceramic becomes more and more extensive,However,it needs to take much more time to st udy it by t raditional experiment s and t he result s are often unilateral for complex prezoelect ric ceramics boundary co ndition and st ress state.In t his view,we use big2scale general piezoelect ric ceramic’s boundary condition and st ress state simulating software ANS YS8.0to carry t hrough comp uter simulation.We gain t he pict ure of piezoelect ric ceramic’s vibration state by using comp uter to analyze t he coupling effect of piezoelect ric ceramic.The experiment result indicates t hat ANS YS8.0can competently deal wit h p roblems about piezoelectric coupling field.It optimizes t he time of creating model and comp uting largly,and it s st rong f unction of post dealing makes researchers can directly obtain data result s and simulating images.K ey w ords:simulation;piezoelect ric ceramic;vibratio n state(责任编辑肖铿,胡小洋)。
基于ANSYS的全陶瓷电主轴动态分析及振动性能测试
ｍ ｅ ｈ ｆ Ａｎｓ ｓ ｉ ｄｉａ ｅ ｈ ｔ ｔ ｅ ｉ ｎｅ ｒ ｎ ｐｅ ｄ ｏ ｐ ｎｄ ｅ ｗａ ａ ｗａ ｏ ｉｓｌ ｉｅ ｏ ａｅ ｓ ｅ ｄ ｉ ｃ ｌ ｖ ｉ ｈｅ ｒ ｇ ｏ ｆ ｔ ｏｄ ｏ ｙ ｎ ｃ ｔ ｄ ｔ ａ ｈｅ ｄ ｓ ｇ ｄ ｗｏ ｋｉ ｇ ｓ ｅ ｆｓ ｉ ｌ ｓ ｆ ｒａ ｙ ｆ ｍ ｔ ｉ ｔ ｄ ｒ ｔ ｔ ｐ ｅ ． ｏｕ ｄ ａ ｏ ｄ ｔ ｅ ｉ ｎ ｏ ｒ ｍ ｔ ｓ ｙｍｐａｈｅ ｉ ｉ ｒ ｔｏ ｓ ｅ ｅ ｔ ｅ ｙ ｔ ｃ ｎ ｂ ｅ ｉ ｅ ｈ ｔ ｔ ｉ ｒ ｔｏ ｏ ｏｒｚ ｄ ｓ ｎｄ ｅ ｗ ａ ｔ ｂ ｅ ｕｎ ｒ ｔ ｅ ｄ ｓｇ ａ ｇ ｈ ｏ ｈ ｔ ｔ ｔｃ ｖ ｂ ａｉ ｎ ｆ ｃ ｉ ｌ ；ｉ ａ ｅ ｖ ｒ ｆ ｄ ｔ ａ ｈｅ ｖ ｂ ａｉ ｎ ｏｆｔ ｍ ｔ ｉ ｅ ｐｉ ｌ ｓ ｓａ ｌ ｄｅ ｈ ｅ ｉ ｎ ｒ ｎ ｅ ｔ ｒ ｕｇ ｈｅ ｖ ｉ ｈｅ ｒ ｓ ｔｏ ｈ ｂ ａ ｉ ｎ ｐｅ ｆ ｒ ｎ ｅ ｔ ｓ ， ｈ ｃ ｌｏ ｖ ｒｆｅ ｈ ｅ ｉ ａ ａ ｉ ａ ． ｅ ｕｌ ｆｔ ｅ ｖｉ ｒ ｔｏ ｒ ｏ ｍａ ｃ ｅ ｔ ｗ ｉ ｈ ａｓ ｅ ｉ ｄ ｔ ｅ ｄ ｓ ｉ ｎ ｇ ｗ ｓｒ ｔ ｏｎ １
摘 要 ：对设 计 的 全 陶瓷 电主 轴 进行 结 构 分析 ，在 Ａ ＹＳ中建立 了其轴 承 一主轴 转 子 系统 三 维 有 限元 模 ＮＳ 型 。采 用 Ｓ ｂｐ ｃ ｕ ｓａｅ法计 算 了前 ６阶 固有频 率 和振 型 。计 算 出临界 转 速 ．然后 通过 Ｈａｍｏ ｉ 析 ｒ ｎｃ分
基于ANSYS的压电圆片振子振动模式分析
Ｃ ＨＥ Ｓ ｅ ｇ．ＹＡＮＧ ｕ Ｎ ｈ ｎ Ｊｎ
（ ｏ ｅ ｅ ｏ ａｉｅ ｏｔｗ ｓ ｒ ｏ ｔｃ ｎｃ ｌＵ ｉ ｒ ｔ ，Ｘ ｎ ７ ０ ７ ，Ｃ ｉ ） Ｃ ｌ ｇ ｆ Ｍ ｒ ，Ｎ ｒ ｅｔ ｎ Ｐ ｌ ｅｈ ｉ ｎｖ ｓ ｙ ｉ １ ０ ２ ｈｎ ｌ ｎ ｈ ｅ ｙ ａ ｅｉ ａ ａ
电 声 基 础 匡０ 闶＠响 囿０ 匡０ ⑥囿＠ ＠ ＠ ， ＠＠ ⑨凹
文 章 编 号 ：０ ２ ８ ８ （ ０ ８ １ — ０ ３ ０ １ ０ — ６ ４ ２ ０ ）２ ０ ４ — ３
基 于 Ａ ＹＳ的压 电圆片振子 振动模 式分析 论 文 ・ ＮＳ
・
陈 晟 ．杨 俊
器 和压 电厚度 计等。所有这些应用 的发展也极大地 促 进 了对压 电振 子特性 的深入研究 。 中， 其 对于纵 向极 化 压 电陶瓷 圆片振子的 自由耦合振动 已引起 很多学者 的 关注【］ ｌ。研究 思路很 多 ， － ３ 但总 的来说 ， 统的解决方 法 传
是 建 立 振 子 的 自 由振 动 方 程 用 数 学 物 理 的 方 法 解 答 。 然 而 。 些 方 法 数 据 处理 起 来 比较 繁 琐 。 且 处 理 的结 这 而 果 不 够 直 观 。 文 运用 有 限 元 理 论 ， 助 Ａ Ｓ Ｓ 件 ， 本 借 ＮＹ 软 对 压 电 圆片 振 子 的 自由振 动 进 行研 究 。
ｌ ｔｄ Ｍａ ｙ ｋ ｎ ｓ ｏ ｉ ｒ ｔ ｎ ｍｏ ｅ ｒ ｒ ｓ ｎ ｅ ｙ ｔ ｅ ｆ ｉ ｌ ｍｅ ｔ ｍｅｈ ｄ Ｔ ｅ ｆ ｎ ａ ｎ ａ ｒ ｑ ｅ ｃ ｓ ｏ ａｅ ． ｎ ｉ ｄ ｆ ｖ ｂ ａｉ ｄ ｓ ａ ｅ ｐ ｅ ｅ ｔ ｄ ｂ ｈ ｉ ｔ ｅｅ ｎ ｔ ｏ ． ｈ ｕ ｄ ｍｅ ｔｌ ｆｅ ｕ ｎ ｙ ｆ ｏ ｎｅ ｌｎ ｉ ｄ ｎ ｉ ｒ ｔ ｎ ａ ｄ ｒ ｄ ａ ｉ ｒ ｔ ｎ ａ ｅ ｖ ｒ ｅ ｒ ｔ ｈ ｈ ｏ ｎ ｅ ｔ ｖ ｌ ｅ ． Ｃｏ ａ ｅ ｏ ｔ ｅ ｔａ ｉ ｏ ａ ｏ ｇｔ ｉ ａ ｖｂ ａｉ ｎ ａ ｉ ｌ ｖｂ ａｉ ｒ ｅ ｙ ｎ ａ ｏ ｔ ｅ ｔ ｅ ｒ ａ ｄ ｔｓ ａｕ ｓ ｕ ｌ ｏ ｏ ｙ ｍｐ ｒ ｄ ｔ ｈ ｒ ｄ ｔ ｎ ｌ ｉ ｍｏ ａ ｎ ｌｓｓ ｍｅｈ ｄ ｆ ｖｂ ａｉ ｎ ｍｏ ｅ ｈ ｒ ｐ ｓ ｄ ｍｅ ｈ ｄ ｉ ｒ ｉ ｌ ｎ ｉ ｃ ． ｄ ｌ ａ ａｙ ｉ ｔ ｏ ｓ ｏ ｉｒ ｔ ｄ ，ｔ ｅ ｐ ｏ ｏ ｅ ｔ ｏ ｓ ｍｏ ｅ ｓｍｐ ｅ ａ ｄ ｄｒ ｔ ｏ ｅ
基于ANSYS的纯压电振子瞬态特性分析
（ ｈ ｎ ｏ ｇＵ ｉｅｓｙｏ ｃ ｎ ｅａ ｄＴ ｃ ｎｌｇ ， ｉｇ ａ ６ ５ ０ Ｃ ｉａ Ｓ ａ ｄ ｎ ｎｖｒ ｔ ｆ ｉ ｃ ｎ ｅｈ ｏ ｙ Ｑ ｎ ｄｏ２ ６ １ ， ｈｎ ） ｉ Ｓ ｅ ｏ
在 理 想 频 率 厂＝ ０ ５ｋ ｚ下 对 振 子 表 面施 加 ４ ． Ｈ
Байду номын сангаас２ ０Ｖ的正 弦变 压载荷 （ ０ 载荷分 布 如 图 ３所示 ） 行 进 瞬态分 析 。
ｌ ｌ ２ ｓ Ｉ ＝ ０ｃ ｌ ＝０ ｉ２ ２ ｏ２ Ｖ ０ ｎ ０ ｓ １ ３２０ｏ２ １ ＝０ｓ １ ＝０ｃ ２０ｎ Ｖ ｓ ｉ２
ｊ 薪 ｉ
计 斩
． … …
ｃ，
／ ， 《 ｚ ≯
＇ ；
… 一 … … … … … … … … … 一 … 一 … … … … … … … … … － … … … … … … － ＿ － … … … … … … － －＿
触持电棚 ２２ 第 ｏ 第 期 ０ 年 ４卷 ４ ｌ
（ ）一 阶 纵 向振 动 Ｌ 振 型 ａ ｌ （ ）二 阶弯 曲振 动 Ｂ ｂ ２振 型
图 ２ 纯压 电振 子 的模 态 分 析 结 果
根 据模 态简 并原则 ， 选定 ３ ４ Ｈｚ ８～ ３ｋ 作为 频率 范围， 并对 压 电振子 两 表 面 施加 幅值 为 ２０Ｖ的 正 ０ 弦 电压 进 行 谐 响 应 分 析 ， 到 振 子 的 理 想 频 率 为 得
并能在各场间实现有限耦合的 Ｓ ＬＤ 单元。 Ｏ Ｉ５ 压 电振 子 的性能 参 数 Ｊ压 电振 子材 料 为 锆钛 ：
酸铅 （ ｉ 缩 写 ：Ｚ ４ ， 料 密 度 为 ７５０ ＰＺＴＯ ）（ ＰＴ ） 材 ０
压电陶瓷振动的有限元分析ansys
1.1压电材料概述 1.1.1压电效应1880年法国物理学家皮埃尔和雅各居里兄弟在研究石英晶体的物理性质时 发现：当沿着晶片的某些方向施加作用力使晶片发生变形后， 晶片上相对的两个 表面会出现等量的正负电荷，电荷的密度与施加的力的大小有关， 这种现象称为 压电现象，具有压电现象的介质称为压电体。
压电效应反应了晶体的弹性性能与介电性能之间的耦合。
当对压电陶瓷施加一个与极化方向平行的压力F ,如图1.1（ a ）所示，陶瓷片将产生压缩变形，片 内的正、负束缚电荷之间的距离变小，极化强度也变小。
因此，原来吸附在电极 上的自由电荷，有一部分被释放，片内的正、负电荷之间的距离变大，极化强度 也变大，因此电极上又吸附一部分自由电荷而出现充电现象。
这种由机械效应转变为电效应的现象就是压电效应。
压电效应包括正压电效应和逆压电效应。
如图 所示： 图1.1压电效应示意图：（a ）正压电效应（b ）负压电效应正压电效应：当压电晶体在外力作用下发生形变时，在它的某些相对应的面 上产生异号电荷，这种没有电场作用，只是由于形变产生的极化现象称为正压电 效应。
逆压电效应：当压电晶体施加一电场时，不仅产生了极化，同时还产生了形 变，这种由电场产生形变的现象称为逆压电效应。
1.1.2压电陶瓷的诞生与发展具有压电效应性能的陶瓷称为压电陶瓷， 1942年美国麻省理工学院绝缘研究室发现，在钛酸钡铁电陶瓷上施加直流高压电场，使其自发极化沿电 场方向择优取第一章绪论向，除去电场后仍能保持一定的剩余极化，使它具有压电效应，从此诞生了压电陶瓷。
钛酸钡（BaTiO s ）陶瓷的发现促进了压电材料的发展，它不但使压电材料从一些单晶体材料发展到压电陶瓷等多晶体材料，而且在压电性能上也有了大幅度提高。
当今广泛应用的压电陶瓷是PZT，即Pb Zr,Ti O3压电陶瓷，其压电效应强，稳定性好。
它是由美国学者B.贾菲等人于1954年发现的PbZrO3 - PbTiO3二元系固溶体压电陶瓷，其机械品质因数约为钛酸钡（BaTiO 3）陶瓷的两倍。
基于Ansys的压电传感器设计与振动分析-毕业答辩
三.研究内容
压电传感器的工作原理是基于某些材料的压电效应。如 图1（a）所示。压电力，速度或加速度传感器可以通过 使用正向压电效应来制造。当对压电晶体施加电场时， 晶体不仅产生极化而且产生应变和应力的现象称为反压 电效应，如图1（b）所示。当电场不强时，应变与外部 电场线性相关。压电位移或力跟随器可以通过反压电效 应制成。
三.研究内容
通过实验可以验证哪种方法在模拟中是准确的。
图5 前六阶振型
三.研究内容
当弹性膜片的弹性模量恒定时，研究了弹性膜片宽度对 固有频率的影响。内圈弹性膜片的宽度为L1，外圈弹性 膜片的宽度为L2作为设计变量，石英晶片的极限强度为 SMAX作为状态变量，传感器的一阶固有频率是F作为目 标函数。这里采用的实际目标函数是10000-f。原因是目 标函数将最小值作为优化目标，从而可以最大化传感器 的一阶固有频率。通过优化分析得到弹性膜片宽度与固 有频率之间的关系曲线。见图所示。
基于Ansys的压电传感器设 计与振动分析
----硕士毕业答辩
目录
1. 研究背景及意义 2.国内外研究进展 3.研究内容 4.总结与展望
一.研究背景及意义
随着科学技术的进步，压电材料将应用于人类生产和 生活的各个角落。在最初发现石英晶体具有一定的压 电特性后，经过进一步的研究和探索，人们开发出了 各种压电晶体，然而，如何使石英压电传感器在动态 模式下稳定工作一直是一个重要的课题。传感器分析 的基本问题是建立传感器的数学模型，揭示传感器弹 性体的应力和应变关系，并解释传感器的机械特性和 耦合尺寸之间的本质关系。理论上，对于某石英压电 传感器的力学模型，可以通过三维运动方程加边界条 件精确求解振动规律，但极其复杂的运动方程和边界 条件使得难以获得解析解。
ANSYS在压电智能板振动主动控制中的应用
收稿日期:2006-02-27作者简介:郑骥(1979-),男,安徽滁州人,硕士研究生,从事结构主动控制研究。
文章编号:1007-6743(2006)02-0014-03ANSY S 在压电智能板振动主动控制中的应用郑 骥,赵东东,张京军(河北工程大学土木学院,河北邯郸 056038)摘要:利用ANSY S 软件对粘贴有压电材料的智能板结构建立有限元模型,通过瞬态动力学分析(时间历程分析)并编写APD L (ansys parameter design language )程序实现了对智能板的振动控制仿真。
首先定义程序所需要的荷载步,并利用每一荷载步下的节点应变作为反馈信号传给压电片,考虑到板是二维结构,其节点应变可用节点挠度的二阶差分来代替。
仿真结果表明了该方法的合理性,为进一步研究智能结构振动主动控制奠定了良好的基础。
关键词:ANSY S ;智能结构;主动控制中图分类号:TH113.1 文献标识码:A 振动主动控制是当前振动工程领域内的高新技术,是动力学、控制、计算机与材料科学等诸多学科的综合。
近年来,随着航空航天、机器人技术的发展,柔性结构的振动控制已经成为人们研究的热点和迫切需要解决的问题之一。
用于振动控制的智能结构中,多用压电材料作为传感器和驱动器,将其贴在梁、板等结构元件的表面或嵌入结构内部。
与传统的振动主动控制方法相比,智能结构可以在不明显改变受控结构的质量和体积的条件下,达到自适应调节减振的目的。
目前对采用压电材料进行主动控制研究的结论通常采用数值仿真[1]分析和实验[2]模型进行验证其有效性。
由于数值仿真需要建立系统的状态空间方程,对于复杂结构很难实现且过程繁琐,而实验模型则需要高昂仪器设备和实验费用。
K arag ülle 等[3]对两自由度弹簧系统振动分别用数值方法计算和ANSY S 软件仿真,得到的结果是一致的,并利用ANSY S 软件对智能结构建模完成了智能梁振动主动控制的仿真过程。
基于Ansys的压电陶瓷材料振动特性仿真与研究
基于Ansys的压电陶瓷材料振动特性仿真与研究【摘要】本文以智能材料中压电陶瓷为研究对象，对其进行的主动振动控制。
本论文总结分析了压电陶瓷的特性和相关的物理学方程，依据理论分析的结果，确定PZT材料的具体参数的选择，分析压电陶瓷特性和力学模型。
利用Ansys 仿真软件对梁进行模态分析得到各阶模态下的振动频率，观察振动情况。
并分析PZT在静态、动态电压下的电致伸缩特性。
然后给梁上的压电陶瓷一个激振，观察和分析梁的振动情况。
最后在梁的某个位置粘贴第二片压电陶瓷，对其施加与第一片相反的电压，抑制梁的振动，观察分析抑制情况。
【关键词】压电陶瓷;Ansys仿真;PZT1.引言压电材料是指具有压电效应，能够实现电能与机械能相互转换的晶体材料，受到压力作用时会在两端面间出现电压，进而表现出压电效应[1]。
压电效应的机理是：具有压电性的晶体对称性较低，当受到外力作用发生形变时，晶胞中正负离子的相对位移使正负电荷中心不再重合，导致晶体发生宏观极化，而晶体表面电荷面密度等于极化强度在表面法向上的投影，所以压电材料受压力作用形变时两端面会出现异号电荷[2]。
反之，压电材料在电场中发生极化时，会因电荷中心的位移导致材料变形[3]。
压电陶瓷具有价格低廉、易于批量生产等优点，已被广泛应用于社会生产的各个领域，尤其是在超声领域及电子科学技术领域中，压电陶瓷材料已逐渐处于绝对的优势支配地位，如医学及工业超声检测、水声探测、压电换能器、超声马达、显示器件、电控多色滤波器等[4]。
随着现代高科技的迅猛发展，智能结构和器件广泛应用于信息技术、新材料技术和航天等高技术领域，并日益显示出其巨大的优越性[5]。
近年来，各国都在积极研究功能陶瓷，研究的重点大都是从老材料中发掘新效应，开拓新应组织和结构入手，寻找新的压电材料。
特别值得重视的是随着材料展，目前国际上对压电材料的应用研究十分活跃[6]。
2.压电陶瓷材料的仿真与分析选取了solid226 ，Solid226单元需要介电常数ε，压电常数数（压电应力矩阵[e]或者压电应变矩阵[d]），弹性常数（柔度矩阵[S]或刚度矩阵[C]），以及密度。
基于ANSYS的全陶瓷电主轴动态分析及振动性能测试
从而提高国产加工中心等数控机床在国际市场上的竞争 力， 使我国在该领域的研究和应用达到世界领先水平[1，2]。
电主轴部件是整个加工中心的核心部件， 其动态特 性将对机床的加工精度、 被加工表面的质量和生产率产 生很大影响。 因此， 对高速电主轴进行动态分析， 具有 十分重要的意义。
笔 者 应 用 有 限 元 软 件 ANSYS 对 实 验 室 最 新 设 计 的 全陶瓷电主轴进行了动态分析， 确定其固有频率和振 型， 得到在不同激励下的动态响应， 其结果为评价电主 轴的动态特性提供了依据。 同时对组装好的电主轴进行 了振动性能测试， 测试结果表明在设计转速范围内电主 轴不会发生共振， 符合要求， 同时验证了有限元分析结 果是合理可行的。
5 在 ANSYS 中 的 谐 响应分析
of motorized spindle
固有频率 临界转速
阶
/Hz
/r·min-1
2 2145.8 128748
3 2147.1 128826
4 3141.4 188484
5 3180.8 190848
6 4284.3 257058
谐响应分析用于分析连续的周期性载荷在结构系统 中产生的持续的周期响应 （谐响应）， 以及确定线性结 构承受 （随时间） 按正弦 （简谐） 规律变化的载荷时的 稳态响应。 这种只计算结构的稳态受迫振动， 发生在激 励开始时的瞬态振动则不在谐响应分析中考虑。 谐响应 分析是一种线性分析， 但是也可以分析存在预应力的结 构 [8]。 谐响应主要采用三种方法求解： Full 法、 Reduced
为： n=60×f。 将主轴的固有频率转化为临界转速见表 2。
主轴设计的最大工作转 表 2 电主轴的临界转速
速 为 30000 ~40000r/min， 从 Tab.2 Limited rotate speeds
压电陶瓷ansys分析
你好，我正在做一个ANSYS的压电仿真的瞬态分析，我的问题是怎么在耦合的电极上加载正弦电压载荷？2010-09-09 16:012405401032|分类：电脑/网络|浏览2094次我的命令流：/filename,fengming/prep7et,1,plane42,,,1et,2,plane13,7,,1MP,ex,1,9.8e10 !黄铜MP,nuxy,1,.32MP,dens,1,7500mp,dens,2,7750 !PZT材料密度mp,perx,2,919 !介电常数mp,pery,2,826mp,perz,2,919tb,piez,2 !压电常数tbdata,2,-5.3512tbdata,5,15.784tbdata,8,-5.3512tbdata,10,12.295tbdata,15,12.295tb,anel,2,,,0 !弹性常数tbdata,1,1.2035e11,7.509e10,7.5179e10tbdata,7,1.1087e11,7.509e10tbdata,12,1.2035e11tbdata,16,2.1053e10tbdata,19,2.1053e10tbdata,21,2.2584e10k,1,0,0 !建模k,2,1e-2,0k,3,1.5e-2,0k,4,1.5e-2,1e-4k,5,1e-2,1e-4k,6,1e-2,2e-4k,7,0,2e-4k,8,0,1e-4l,1,2l,2,3l,3,4l,4,5l,5,6l,6,7l,7,8l,8,5l,8,1l,5,2al,1,10,8,9al,2,3,4,10al,5,6,7,8type,1mat,1mshape,0,2Dmshkey,1amesh,1,2type,2mat,2mshape,0,2Dmshkey,1amesh,3nsel,s,loc,x,1.5e-2 !位移约束d,all,ux,0d,all,uy,0/solu !瞬态分析antype,transtrnopt,fulltintp,,0.25,0.5,0.5outpr,all,alloutres,all,allkbc,1nsel,s,loc,y,1e-4 !耦合电压cp,1,volt,alld,all,volt,0*do,i,1,50,1t=i/400000time,tnsubst,1v=40*sin(2*3.14159*42881*t)nsel,s,loc,y,2e-4cp,2,volt,alld,all,volt,vsolve*enddofinish这样加载问题在哪里，对不对？谢谢分享到：2010-09-09 18:30提问者采纳我先看看然后晚上回复你！方法是对这呢，不过在写程序的过程出现了一些小失误你看下我这样改下可以不？与你表达的意思是一致的，只是后面的东西有些改正，前面的材料以及模型尺寸楼主你就自己拿捏一些精度与单位就可以了allsel,allfinish/solu !瞬态分析antype,transtrnopt,fulltintp,,0.25,0.5,0.5outpr,all,alloutres,all,allkbc,1nsel,s,loc,y,1e-4 !耦合电压cp,1,volt,allnsel,s,loc,y,2e-4cp,2,volt,allnsel,s,loc,y,1e-4d,all,volt,0allsel,all*do,i,1,50,1t=i/400000time,tnsubst,1v=40*sin(2*3.14159*42881*t)nsel,s,loc,y,2e-4d,all,volt,vallsel,allsolve*enddofinish提问者评价太谢谢你了，对了的，没有警告了！来自：求助得到的回答asel,s,,,56 ！选择面736cp,1,volt,all ！为该表面施加耦合电压*get,n1,node,,num,min ！定义该表面电极节点d,n1,volt,0 ！加0v电压ALLSEL,ALL ！选择全部asel,s,,,896 ！选择面896nsla,s,1 ！选择面上所有节点cp,2,volt,all ！施加耦合电压*get,n2,node,,num,mind,n2,volt,200 ！加幅值200电压ALLSEL,ALL !选择全部。
有限元计算-压电振子(实验报告)
tbdata,12,11.5E10 ! c33
tbdata,16,3.0E10 ! c44
tbdata,19,2.6E10 ! c55
tbdata,21,2.6E10 ! c66
tb,PIEZ,1 ! 压电（应力）常数, C/m^2
tbdata,3,-5.2 ! e13
实验项目名称
有限元计算-压电振子
实验成绩
实验者
蔡铭辉
专业班级
材物1401
组别
同组者
实验日期
年 月 日
第一部分：实验分析与设计
一、实验内容描述（问题域描述）
利用Ansys通过命令行方式进行长方形压电振子建模，并对其进行加载电压的多场耦合模拟，分析研究振动频率、最大位移（形变量）与模型尺寸之间的关系。
二、实验基本原理与设计（包括实验方案设计，实验手段的确定，试验步骤等，用硬件逻辑或者算法描述）
NSEL,ALL
FINISH
！求解
！短路状态下（谐振状态）
/SOLU
ANTYP,MODAL
MODOPT,LANB,10,50000,150000 ！定义分块LANCZOS法
MXPAND,10
/COM,
NSEL,S,LOC,Z,0
D,ALL,VOLT,0
NSEL,S,LOC,Z,H
D,ALL,VOLT,0
MXPAND,10
/COM
NSEL,S,LOC,Z,H
DDELE,ALL,VOLT
!删除VOLT约束
CP,1,VOLT,ALL
!耦合VOLT约束
NSEL,ALL
SOLVE
*GET,F3,MODE,4,FREQ
压电陶瓷的自振频率分析及详细过程
压电变换器的自振频率分析及详细过程所有图片大小可以自己调节1.在开始菜单中启动Mechanical APDL Product launcher，如图1所示其中Working Directory代表你的工作路径，也就是所有ansys的生成文件都存在这个文件加下。
而job Name 就代表你所做实例的名字，实例中我们将目录设为F盘，工作名设为model2.最后点击run运行进入ansys界面。
如图2所示。
图22.定义单元类型拾取菜单Main M--Preprocessor-Element Type-Add/Edit/Delete,如图3所示。
然后单击，图3然后单击add按钮。
如图4所示，单击ok ,选择solid5单元，这样就完成了对单元的选择。
然后再单击图3中的close按钮。
图43.定义材料属性拾取菜单Main Menu-Preprocessor--Material Props-Material Models.弹出对话框，如图5所示。
定义相对介电常数。
在右侧列表中依次双击"electromagnetic", "relative permittivity", "orthotropic".输入相关参数。
图5在右侧再点击Density,如图6所示。
定义材料的密度。
密度为7500，单击ok。
图6在右侧在依次点击“Piezoelectrics”“Piezoelectric matrix”定义压电应力矩阵，如图7所示。
输入相关参数。
单击ok。
图7在右侧列表中依次双击"Structural", "Linear", "Elastic", "Anisotropic"如图8所示，定义各向异性弹性矩阵设置，输入参数，单击ok。
图8点击Material-Exit完成对材料的定义。
4.建模选取菜单Parameters-Scalar Parameters 在selection中输入L=10E-3然后点击accept，再输入h=20e-3,点击accept，按照相同的步骤输入，w=10e-3,a3=1000.如图9所示。
ANSYS压电分析
A N S Y S压电分析-CAL-FENGHAI-(2020YEAR-YICAI)_JINGBIAN压电分析是一种结构—电场耦合分析，压电效应是石英和陶瓷等压电材料的自然属性。
给压电材料加电压会产生位移，反之使压电材料振动则产生电压。
一个典型的压电分析的应用是压力换能器。
ANSYS压电分析类型（仅在ANSYS Multiphysics或ANSYS Mechanical产品中）有静态、模态、预应力模态、谐响应、预应力谐响应和瞬态分析。
压电分析需要用下列单元类型之一：PLANE13，KEYOPT（1）=7，耦合场四边形单元；SOLID5，KEYOPT（1）=0或3，耦合场六面体单元；SOLID98，KEYOPT（1）=0或3，耦合场四面体单元；PLANE223，KEYOPT（1）=1001，耦合场8节点四边形单元；SOLID226，KEYOPT（1）=1001，耦合场20节点六面体单元；SOLID227，KEYOPT（1）=1001，耦合场10节点四面体单元；KEYOPT选项激活压电自由度、位移和电压。
对于SOLID5和SOLID98，设置KEYOPT（1）=3激活压电仅有选项。
通过使用NLGEOM、SSTIF和PSTRE S命令，KEYOPT选项可用于大变形和大应变的影响。
对于PLANE13单元，大变形和大应变在KEYOPT（1）=7时有效。
对于SOLID5和SOLID98单元，大变形和大应变在KEYOPT（1）=3时有效，小变形和小应变在KEYOPT（1）=0时有效。
自动求解控制对于压电分析是无效的，SOLCONTROL 缺省设置仅对纯结构或纯热分析是有效的。
对于大变形压电分析，必须使用非线性求解命令指定设置。
1.? ? 注意要点分析可以是静态、模态、预应力模态、谐响应、预应力谐响应和瞬态分析。
应注意以下要点：（1）对于模态分析，推荐使用Block Lanczos法（缺省）求解。
（2）对静态、全谐响应或全瞬态分析，选用稀疏矩阵（SPARSE）求解器或Jacobi共轭梯度（JCG）求解器。
Ansys压电分析
ansys压电分析压电效应分析是一种结构－电场耦合分析。
当给石英和陶瓷等压电材料加电压时，它们会产生位移，反之若使之振动，则会产生电压。
压力传感器就是压电效应的一种典型的应用。
压电分析（ANSYS/Multiphysics或ANSYS /Mechanical软件包提供这种分析）类型可以是静力、模态、预应力模态、谐波、预应力谐波和瞬态分析。
压电分析只能用下列单元类型之一：PLANE13, KEYOPT(1)=7，耦合场四边形实体单元SOLID5, KEYOPT(1)=0或3，耦合场六面体单元SOLID98, KEYOPT(1)=0或3，耦合场四面体单元KEYOPT选项激活压电自由度：位移和电压。
对于SOLID5和SOLID98，KEYOPT(1)=3仅激活压电选项。
注意：如果模型中激活了至少一个带有压电自由度（位移和VOLT）的单元，则需要用到VOLT自由度的所有单元必须是上面三种压电单元其中之一。
而且，所有的这些单元均需激活压电自由度。
如果不希望在这些单元中存在压电效应，则需给材料定义非常小的压电特性。
压电KEYOPT用NLGEOM，SSTIF，PSTRES命令可用大挠度和应力刚化作用（参见《ANSYS Commands Reference》对这些命令的更多信息，参见《ANSYS Structural Analysis Guide》及《ANSYS，Inc. Theory Reference》的第三章关于大挠度及应力刚化功能的更多信息）。
对PLANE13，通过设置KEYOPT（1）＝7可用大挠度及应力刚化功能。
对SOLID5及SOLID98通过设置KEYOPT（1）＝3可用大挠度及应力刚化功能。
而且小挠度及应力刚化选项可以通过KEYOPT（1）＝0使用。
注意－对压电分析不能使用自动求解控制。
SOLCONTROL缺省设置只能对纯结构或纯热分析使用。
对大挠度压电分析，必须用非线性求解命令定义有关设置。
用ANSYS分析PZT压电方片振动的模态特性
用ANSYS分析PZT压电方片振动的模态特性张君善【摘要】基于有限元软件ANSYS 14.0的多物理场耦合仿真技术,通过建模对PZT 压电陶瓷方形薄片的振动特性进行了模态分析.通过介绍压电性和多物理场耦合的原理,设置了材料的压电常数和弹性常数矩阵,选用solid 226多物理场单元进行有限元建模和仿真.对几何尺寸为1.5 mm×1.5 mm×0.1 mm的PZT-5H压电陶瓷方片做了振动模态分析,提取了前十阶的模态频率和振型,此分析对高频压电超声换能器的制备有一定的帮助.【期刊名称】《宁夏师范学院学报》【年(卷),期】2015(036)003【总页数】5页(P63-66,85)【关键词】有限元;ANSYS;模态分析;压电振动【作者】张君善【作者单位】武警工程大学装备工程学院,陕西西安710086【正文语种】中文【中图分类】TH164压电效应是有关机械能和电能的耦合转换效应，最初由皮埃尔·居里和雅克·居里兄弟于1880年发现.具有压电效应的材料是压电材料，大致可分为压电单晶、压电陶瓷、高聚物压电材料和压电复合材料等种类.压电材料广泛应用于压电传感器、谐振器、换能器和致动器等压电器件中，其中PZT(Pb(Zr1-xTix)O3,锆钛酸铅)系列压电陶瓷在压电超声换能器领域中的应用尤为广泛.压电超声换能器是利用压电效应进行电能和机械能转换的器件，其中的核心元件为压电有效层，视换能器用途的不同有效层一般为不同形状、尺寸和厚度的压电材料；压电材料在电场和应力作用下，将电势能和机械能进行耦合变换，并以电信号或者超声信号输出.研究压电层的自由振动特性及瞬态电压激励下的振动特性对于超声换能器的制备有重要意义[1,2].有限单元法是在数学理论和工程应用结合推动下发展起来的一种数值计算方法，其基于变分原理和加权余量法来求解偏微分方程.随着计算机处理性能的不断发展，有限元法的功能也愈加强大，工程应用非常广泛.商业化的大型有限元软件种类繁多，其中ANSYS是一款广泛使用的大型有限元工程软件，其包含模块丰富，功能强大，版本升级迅速，市场占有率高.在压电换能器的制作过程中，器件功能设计和模拟仿真是前期的重要环节，通过功能仿真可以优化换能器的参数，并预判器件的最终性能.据报道，姜德义、雷辉等使用ANSYS做过压电陶瓷应用于混凝土无损检测的振动状态研究[3,4]，林书玉等对压电陶瓷圆片的振动做了深入理论研究[5]，向阳等研究了压电陶瓷圆环的径向振动[6].在高频压电超声换能器的应用中，压电陶瓷有效层的厚度往往要减薄到几十微米，且出于加工方便，其几何形状通常设计为方片形.本文使用ANSYS 14.0对PZT-5H压电陶瓷方片的振动模态进行了仿真，设置了其压电常数、介电常数和弹性常数，得到了前十阶的模态解，分析了满足换能器设计参数的压电片的自由振动特性.压电性是关于力学量和电学量之间的相互作用，压电方程描述的是压电材料的力学量和电学量之间相互关系.根据压电材料所受的机械边界条件和电学边界条件的不同，选择不同的自变量和因变量，可以得到四类不同类型的压电方程.分别为如下四种：a.机械自由和电学短路条件；b.机械夹持和电学短路条件；c.机械自由和电学开路条件；d.机械夹持和电学开路条件.在进行压电换能器的有限元仿真中通常选用第二类压电方程，取应变量S和电场强度E为自变量，应力量T和电位移D为因变量，其方程的简写形式为其中，cE为电场强度恒定时压电材料PZT的弹性常数矩阵，e为压电应力系数矩阵，et为e的转置矩阵，εS为PZT材料的夹持介电常数(恒应变时).ANSYS中的多物理场耦合分析功能可实现压电分析.项目开始时，在过滤菜单preferences的选项卡中设置，同时选取Structural 和Electric两个优选项.多物理场耦合分析中的耦合是直接耦合，而非顺序耦合.多物理场耦合提供的分析单元有solid 226，solid 227等，考虑到建立的模型是规则的方形陶瓷片，选用三维六面体20节点的solid 226单元作为分析单元.压电常数是描述压电材料把电能转换为机械能，或把机械能转换为电能的转换系数.压电常数反映了压电材料所具有的压电效应的强弱，直接关系到压电换能器输出的灵敏度.以厚度模式工作的压电超声换能器，经常采用极化方向沿厚度方向的PZT-5H压电陶瓷，因为其压电常数相对较大，机电耦合系数较高，有利于增强换能器性能.模态分析一般用于确定结构的振动特性，如固有频率和振型(模态).对于无阻尼自由振动，动力学基本方程中的阻尼项和激励项均为零，方程的形式为式中，M是质量矩阵，K为刚度矩阵.任何弹性体的自由振动都可以视为一系列简谐振动的叠加.方程(1)的简谐振动解形式为式中，ω为固有圆频率，δ0为各节点振幅值，也称为振型.将(2)式代入(1)式得，由于自由振动的各节点振幅不全为零，所以(3)式的系数矩阵的行列式必为零.即由于刚度阵和质量阵均为n阶方阵，故(4)式是关于ω2的n次方程，可以求解出n个固有频率.ω2称为广义特征值.对应于每个固有频率，由(3)式可得一组节点振动的振幅值δ0，称为结构的主振型.最小特征值对应的主振型称为基本振型.ANSYS做模态分析时，所用的模态提取方法有分块法(Block Lanczos)、子空间法(Subspace)、缩减法(Reduced)等，常用的是分块法和子空间法.压电常数是描述压电材料把电能转换为机械能，或把机械能转换为电能的转换系数.压电常数反映了压电材料所具有的压电效应的强弱，直接关系到压电换能器输出的灵敏度.本文中，PZT-5H压电陶瓷方片的极化方向沿厚度z方向，在软件中对其压电常数和弹性常数进行设置输入.首先，选取单元类型为solid 226单元，在单元选项option中设置分析类型为piezoelectric.在材料属性中依次设置，密度DENS=7500，磁导率MURX=1；相对介电常数PERX=919，PERY=919，PERZ=826.需要注意的是，在弹性系数矩阵和压电常数矩阵设置时，ANSYS中的阵元排列顺序和IEEE标准的顺序不同，IEEE标准中矩阵元的下标顺序是X，Y，Z，YZ，XZ，XY；在ANSYS中，程序要求的是按照X，Y，Z，XY，YZ，XZ的输入顺序[7].在设定材料参数时应按顺序进行转换后再输入到程序中.程序中，压电常数矩阵采用压电应力矩阵[e]的形式，其单元值如表1所示.在ANSYS中，压电系数矩阵和弹性系数矩阵各有两种形式，各自对应，也可以相互转换.压电应变系数矩阵与弹性顺度矩阵相对应，与压电应力系数矩阵相应的，压电材料的弹性系数矩阵采用劲度矩阵[c]形式，其形式如表2所示.在ANSYS前处理中进行实体建模，首先，在工作平面上创建一个几何尺寸为1.5 mm×1.5 mm的矩形面，然后，用extrude命令沿法向拉伸0.1 mm生成三维实体.在网格划分时，对六面体的边线设定单元划分数目，分别为z方向的短边划分为4份，xy平面内的长边划分为10份，而后，选用六面体形状以掠扫(mapped)的方式对整个几何实体进行网格划分，程序最终的划分单元数为400.边界条件的施加为，选取压电方片上下表面的中心点，节点编号分别为221和552，对其X和Y方向自由度进行约束.电压加载到压电方片的上下表面，加载电压大小为10 V.分析类型设为Modal，求解方法设定为分块法(Block Lanczos)，模态扩展及提取数目均设为10，频率范围设为0 MHz～1 MHz.模态分析结果中，列举了几个频率和振型图，如图2所示.同样，对压电材料的振动状态变形情况，截取了几个图，如图3所示.使用有限元软件ANSYS 14.0，对几何尺寸为1.5 mm×1.5 mm×0.1 mm的PZT-5H压电陶瓷方形薄片进行了原理阐述和振动模态分析，得到了陶瓷方片前十阶的模态频率和振型，以及陶瓷方片的变形情况.在已有的报道中，适合高频领域的正方形压电片的分析较少，而高频压电超声换能器的阵元尺寸小，且通常加工为方形，所以本文的仿真结果对于高频压电超声换能器的制备具有一定的帮助. [1] Qi W K，Cao W W.Finite Element Analysis and Experimental Studieson the Thickness Resonance of Piezocomposite Transducers[J].Ultrasonic Imaging，1996，18:1-9.[2] Ramesh R，Kara H，Bowen C.R.Finite element modelling of dense and porous piezoceramic disc hydrophones[J].Ultrasonics，2005，43(3):173-181.[3] 姜德义，郑拯宇，李林，等.压电陶瓷片耦合振动模态的ANSYS模拟分析[J].传感技术学报，2003，16(4):452-456.[4] 雷辉，周双娥.用ANSYS软件分析压电陶瓷的振动状态[J].湖北大学学报，2008，30，(1)：29-33.[5] 向阳，王丽坤.压电陶瓷圆环径向厚度振动分析[J].压电与声光，2012，34(4)：561-564.[6] 林书玉，张福成.压电陶瓷圆片振子的多模耦合振动[J].电子学报，1994，12(12)：64-69.[7] GB/T 176-1987 IEEE Standard on Piezoelectricity[S].New York:The Institute of Electrical and Electronics Engineers,1987.ANSI/IEEE Std 176-1987.(SchoolofEquipmentandEngineering，EngineeringUniversityofChineseArmedPoliceForce，Xi'an，Shaanxi 710086)Key words FEM； ANSYS； Modal analysis； Piezoelectric vibration。
ansys压电分析实例Ansys压电分析
ansys压电分析实例Ansys压电分析导读：就爱阅读网友为您分享以下“Ansys压电分析”的资讯，希望对您有所帮助，感谢您对的支持! ansys压电分析压电效应分析是一种结构－电场耦合分析。
当给石英和陶瓷等压电材料加电压时，它们会产生位移，反之若使之振动，则会产生电压。
压力传感器就是压电效应的一种典型的应用。
压电分析（ANSYS/Multiphysics或ANSYS /Mechanical软件包提供这种分析）类型可以是静力、模态、预应力模态、谐波、预应力谐波和瞬态分析。
压电分析只能用下列单元类型之一：PLANE13, KEYOPT(1)=7，耦合场四边形实体单元SOLID5, KEYOPT(1)=0或3，耦合场六面体单元SOLID98, KEYOPT(1)=0或3，耦合场四面体单元KEYOPT选项激活压电自由度：位移和电压。
对于SOLID5和SOLID98，KEYOPT(1)=3仅激活压电选项。
注意：如果模型中激活了至少一个带有压电自由度（位移和VOLT）的单元，则需要用到VOLT自由度的所有单元必须是上面三种压电单元其中之一。
而且，所有的这些单元均需激活压电自由度。
如果不希望在这些单元中存在压电效应，则需给材料定义非常小的压电特性。
压电KEYOPT用NLGEOM，SSTIF，PSTRES命令可用大挠度和应力刚化作用（参见《ANSYS Commands Reference》对这些命令的更多信息，参见《ANSYS Structural Analysis Guide》及《ANSYS，Inc. Theory Reference》的第三章关于大挠度及应力刚化功能的更多信息）。
对PLANE13，通过设置KEYOPT（1）＝7可用大挠度及应力刚化功能。
对SOLID5及SOLID98通过设置KEYOPT（1）＝3可用大挠度及应力刚化功能。
而且小挠度及应力刚化选项可以通过KEYOPT （1）＝0使用。
注意－对压电分析不能使用自动求解控制。
基于ANSYS模态分析的压电陶瓷参数数据的转换及其输入
面向 分析时压电材料 � 可以进行下面的处理 参数输入与转换 � � � +
�
� � � 首先列出了压电陶瓷的基本压电方程 , 为下面的
�
-1 -1 � � 讨论提供了一个框架 . 然后详细讨论了把制造商提供 . (7 ) � � � � 的数据转换为 A 需要的数据的方法 , 这些数据 由于方程 (2 ) 反映的是应力与电位移而不是应变 包括了刚度矩阵 ,介电常数矩阵和压电常数矩阵 . � 与电位移的关系,所以把方程(7 )代入方程 (2 )可得 压电方程组 � � � � + � , (8 ) 压电 体是结构与电学的耦合体 , 具有压电效应 � � 和 逆压电 效应 , 是一 些材料 (� 例 如 , 石英 ,压 电陶 瓷 ) � � � � � ( ) + � , (9 ) � � 的 天然属 性. 压电体具有机电转换功能 , 可以把力效 -1 -1 � � � � � . (10 ) � � � � + (� ) � � 应转换为电效应 . 对于压电体来说 ,电场作用能使它 把方程 ( 7) 和 (10) 组成方 程组与方程 (3),( 4) 对 产生电 位移 ,一些 力也能在某些方向产生电 位移 同 -1 -1
收 稿日 期� 2 009-0 7-10 基 金项 目� 国 家自 � 然科 学基金 项目 " � 基于 压电 转换 的振 动能 量回 收及 整流 变换 技术 的研 究" (N . 50 57702 3) 作� 者简 介� 罗 平 � (1983-), 男, 安徽 亳 州市人 , 广东 技术 师范 学院 教务 处科 员 研 究方 向 计 算机 测试 与控 制, 数字 资源 建设 及自 动化 技术
引
言
O3 ,19 55 年 ,
基于ANSYS的矩形压电振子瞬态特性分析
孽ｉ积分方法求解结构和其他二阶系统方程。 差分方法 替换 。 ；用于求解显示的瞬态分析； 牛顿法用于隐式瞬态问 根据 辛凯维 奇理 论 ， 用式 （ ） 式 （ ） 利 ２和 ３ 得到 牛
ｉ题的求解 。由于振子压电陶瓷机电耦合瞬态分析过 顿法求解的无条件稳定必须满足 ：
皇；程中， 涉及到 材料非线性、 压电 振子与转子接触非线
｛ ｝ “｝ 。 ， 。 ＝｛。 ＝｛ 。 ≠ 通过 Ｉ ！ ｕ ｛。 ｛ ｝ ｛ ｝ ｝ ｕ ｝ ０ 。 ｃ
命令输入或在瞬态本身开始之前执行一个静态分析 ｉ
载荷步 。进行瞬态分析 中的静载荷步时要关闭瞬态 ！
｛ １ ：（ 一Ｏ ） ｕ＋｝＋（ ｛ ｕ＋ ｝ １ ｔ ｛ １ ｉ “ ｝ 一 ｍ Ｘ 时间积分的影响 ， 结束静载荷分析后会 自动开启 时 间积分 的影 响 。瞬态 时间积 分 的影 响 默认 形 式是开
ｉ０ 引 言
ｉ
１ １牛顿法 理论 介绍 ．
对 于 非线 性 问题 ， Ｎ Ｙ Ａ Ｓ Ｓ用 中心差 分 方法 整合 运动方 程 。牛顿 法 使 用 有 限 差 分 法 ， ｔ ｔ 时 在 到
间 间隔 内 ， ： 有
ｉ
起动 和关 断瞬 态过程所 表 现的特 性是评 价超声
饭 持电棚 ２２ ． 卷 ３ 一 ｌ！第 ０ 第 ． 期 ＿ … ０ 年 ４妻曼 ＾ ． …
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第30卷第1期2008年3月湖北大学学报(自然科学版)Journal of Hubei University (Natural Science )Vol.30 No.1 Mar.,2008 收稿日期:2006206202作者简介:雷辉(19812 ),男,硕士生文章编号:100022375(2008)0120029205用ANSYS 软件分析压电陶瓷的振动状态雷辉,周双娥(湖北大学数学与计算机学院,湖北武汉430062)摘 要:近年来,压电陶瓷的应用日趋广泛.但是由于压电陶瓷片的边界条件和应力状况比较复杂,利用传统实验手段对其研究不仅耗时费力,而且其结果具有很强的局部性,因此利用大型通用仿真软件ANSYS 8.0来进行计算机仿真.通过对压电陶瓷片中的耦合效应进行计算机模拟分析,得出压电陶瓷的振动状态图.实验结果表明ANSYS 8.0在处理压电耦合场这方面有很强的处理能力,大大简化了建模和计算,强大的后处理功能更是让研究者能够很直观地获得数据结果和模拟图像. 关键词:仿真;压电陶瓷;振动状态 中图分类号:TP302 文献标志码:A1 引言计算机仿真技术是以多种学科和理论为基础,以计算机及其相应构件为工具,通过虚拟试验的方法来分析和解决问题的一门综合性技术[1].近年来,压电陶瓷的应用日趋广泛,而在实际应用中,特别是将压电陶瓷技术应用于混凝土结构的监测中,由于压电陶瓷片的边界条件和应力状况比较复杂,利用传统实验手段对其研究不仅耗时费力,而且其结果具有很强的局部性[2].因此利用计算机仿真技术对压电陶瓷进行研究具有较好的理论与实际意义.本文中利用大型通用有限元分析软件ANS YS 8.0,对压电陶瓷片中的耦合效应进行模拟分析,并得出其模态和谐振态,实验表明ANS YS 8.0能很好地解决压电陶瓷片的压电耦合问题.图1 处理器模型2 ANSYS 仿真原理ANS YS 软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件,由世界上最大的有限元分析软件公司之一的美国ANS YS 开发,它能与多数CAD 软件接口,实现数据的共享和交换.20世纪90年代该软件开始在我国的机械制造、航空航天、汽车交通、铁道、石油化工、能源等领域得到应用,为各领域中产品设计、科学研究作出了很大的贡献[3].ANS YS 软件使用统一的集中式数据库来存储所有模型数据和求解结果(见图1)[4].模型数据(包括实体模型和有限元模型、材料等)通过前处理器写入数据库;载荷和求解结果通过求解器写入数据库;后处理结果通过后处理器写入数据库.3 处理过程3.1 定义材料参数 材料参数包括定义单元类型,这里选取了solid226,并在它的option 选项里选择压 湖北大学学报(自然科学版)第30卷30电这个选项.然后定义压电陶瓷的密度、介电常数、刚度系数和压电常数,后两个参数是用矩阵的形式来表示的.为了方便后期的网格划分处理,添加了一个单元类型mesh200,它主要用来进行面划分,下一节将详细介绍.定义材料参数的部分代码及注释如下所示:ET,1,SOL ID226,101 !定义solid226单元类型KEYOP T,1,1,1001!在solid226选项中选择压电选项ET,2,M ESH200!3!定义mesh200单元类型KEYOP T,2,1,7KEYOP T,2,2,0MPTEMP,,,,,,,,MPTEMP,1,0MPDA TA,DENS,1,7600!定义压电陶瓷密度TB,AN EL,1,1,21,0!定义压电陶瓷的刚度系数TB TEMP,0TBDA TA,,1.32e11,7.1e10,7.3e10,,,TBDA TA,,1.32e11,7.3e10,,,,1.15e11TBDA TA,,,,,3.0e10,,TBDA TA,,2.6e10,,2.6e10,,,MPTEMP,,,,,,,,!定义压电陶瓷的介电常数MPTEMP,1,0MPDA TA,PERX,1,,7.124e-9MPDA TA,PER Y,1,,7.124e-9MPDA TA,PERZ,1,,5.841e-9TB,PIEZ,1,,,0!定义压电陶瓷的压电系数TBMODIF,1,1,TBMODIF,1,2,TBMODIF,1,3,-4.1TBMODIF,2,1,TBMODIF,2,2,TBMODIF,2,3,-4.1TBMODIF,3,1,TBMODIF,3,2,TBMODIF,3,3,14.1TBMODIF,4,1,TBMODIF,4,2,TBMODIF,4,3,TBMODIF,5,1,TBMODIF,5,2,10.5TBMODIF,5,3,TBMODIF,6,1,10.5TBMODIF,6,2,TBMODIF,6,3,3.2 建立模型及网格划分 首先新建一个长为0.005m,宽为0.001m的矩形.然后选取单元类型为mesh200,对该矩形进行面划分,其中长和宽分别划分8份和10份,结果如图2所示.然后用操作选项中的Ext rude命令将其扩展成一个已划分好的圆盘体,其中要在Extrude选项中将单元类型选择为solid226,并在要划分的数目中填入10.接着用Extrude命令将其扩展成一个圆盘,其扩展结果如图3所示.第1期雷辉等:用ANSYS 软件分析压电陶瓷的振动状态31 图2 面划分图 图3 体划分图3.3 添加约束条件和负载 添加的约束条件是在上、下表面的圆心处分别添加位移约束,使其只能沿纵向方向移动.具体代码及注释如下:nsel ,s ,loc ,y ,0nsel ,r ,loc ,z ,0nsel ,r ,loc ,x ,0!选取下表面圆心上的点d ,all ,ux ,0!使其不能沿x 方向运动d ,all ,uz ,0!使其不能沿y 方向运动nsel ,s ,loc ,y ,0.001!选取上表面圆心上的点nsel ,r ,loc ,z ,0!同理nsel ,r ,loc ,x ,0d ,all ,ux ,0d ,all ,uz ,0给圆盘添加的电压负载是在圆盘的上下表面的两个电极上加上耦合电压,其电压值分别为5伏和0伏.这样使得圆盘上下电势差为5伏.其代码及注释如下:nsel ,s ,loc ,y ,0!选取下表面所有节点cp ,1,volt ,all !为下表面添加耦合电压3get ,n_bot ,node ,0,num ,min !定义下表面的电极节点nsel ,s ,loc ,y ,0.001!选取上表面所有节点cp ,2,volt ,all !为上表面添加耦合电压3get ,n_top ,node ,0,num ,min !定义上表面的电极节点nsel ,all !选取所有节点d ,n_bot ,volt ,0!下表面加0伏电压d ,n_top ,volt ,5!上表面加5伏电压3.4 静态和模态下的处理 首先选择处理状态为静态,并在此状态下得出其静态电容,代码及注释如下:/SOL !进入处理环节AN T YPE ,0!选择静态处理/STA TUS ,SOL U SOL V E !求解3get ,cs ,node ,n_top ,rf ,chrg !得到上电极电量值fini !处理结束3SET ,cs ,abs (cs )/5!得到静态电容接着进行模态分析,设定它有20个子步,频率范围为02200000Hz ,并且将上表面短路,代码及注释如下:/SOL U !进入处理环节AN T YPE ,2!选择模态处理MODOPT ,L ANB ,20!定义处理方法和子步数EQSL V ,SPAR !选择处理器MXPAND ,20,,,1!要显示的子步数L UMPM ,0PSTRES ,0MODOPT ,L ANB ,20,0,2000000,,OFF !定义频率范围32 湖北大学学报(自然科学版)第30卷d ,n_top ,volt ,0!上表面短路nsel ,all!选择所有节点/STA TUS ,SOL USOL V E !求解3.5 后处理 后处理是指检查并分析求解的结果的相关操作.这是分析中最重要的环节之一,因为任何分析的最终目的都是为了研究作用在模型上的的载荷是如何影响设计的[5].检查分析结果可使用两个后处理器:POST1(通用后处理器)和POST26(时间历程后处理器).POST1允许检查整个模型在某一载荷步和子步(或对某一特定时间点或频率)的结果.POST26可以检查模型的指定节点的某一结果项相对于时间、频率或其它结果项的变化.在求解时,ANS YS 将结果写入结果文件,进行后处理时,结果文件必须存在且可用.结果文件名的后缀取决于分析类型,对于结构分析的结果文件的后缀为RST [6].本次实验只用到了POST1.在后处理中利用一个循环语句可以得出不同频率下的动态电容、动态电感、电量以及阻抗值.其代码和注释如下:/POST13SET ,nmodes ,20!定义nmodes =203dim ,C ,array ,nmodes!定义电容数组3dim ,L ,array ,nmodes!定义电感数组3SET ,PI2,233.14159!定义PI2=233.14159set ,first!设定第一个子步/com ,3do ,i ,1,nmodes !定义一个循环,从1到20步3get ,Fi ,mode ,i ,freq !得到该步的频率值3get ,Qi ,node ,n_top ,rf ,chrg !得到该步的电量值3SET ,Omi ,Pi23Fi !将线速度转化为角速度3SET ,C (i ),(Qi/Omi )332!计算相关的动态电容3SET ,L (i ),1/(Omi 3323C (i ))!计算相关的动态电感/com ,Mode %i %!在输出窗口中显示第几步/com ,Resonant f requency F =%Fi %Hz !在输出窗口中显示频率值/com ,Dynamic capacitance C =%C (i )%F !在输出窗口中显示动态电容值/com ,Dynamic inductance L =%L (i )%H !在输出窗口中显示动态电感值/com ,charge Q =%Qi %C!在输出窗口中显示电量值/com ,Impedance R =%5/(233.141593Fi 3Qi )%!在输出窗口中显示阻抗值/com ,set ,next !进入下一个子步3enddo!循环结束FINISH同时可以在主菜单的Animation 选项下看一下压电圆盘在各个不同频率下的振动状态,这里截取了几个振动图,如图4所示:(a )频率为29728Hz (b )频率为65741Hz (c )频率为53476Hz图4 振动状态图第1期雷辉等:用ANSYS软件分析压电陶瓷的振动状态33 4 结束语本文对压电陶瓷圆盘添加了约束和对称负载,解决了压电陶瓷压电场与结构场的耦合问题,并且最终得到了压电圆盘在静态下的电容值,以及它在模态下的动态电容、动态电感、电量和阻抗,而且还可以观察圆盘在各个不同频率下的振动状态,以便于今后对材料性能的研究以及对材料的改进.同时,不难发现Ansys8.0在处理压电耦合场这方面有很强的处理能力,像Solid226就是专门针对压电分析而定义的,而且Ansys8.0在其材料库中建立了相关的压电材料,因此大大简化了有限元的建模和计算,强大的后处理功能更是让研究者能够很直观的获得数据结果和模拟图像.参考文献:[1]林书玉,张福成.压电陶瓷圆片振子的多模耦合振动[J].电子学报,1994,12:43249.[2]姜德义,郑拯宇.压电陶瓷片耦合振动模态的ANSYS模拟分析[J].传感技术学报,2003,12:9216.[3]陈大任.压电陶瓷微位移驱动器概述[J].电子元件与材料,1994,2:33240.[4]邵蕴秋.ANSYS8.0有限元分析实例导航[M].北京:中国铁道出版社,2004.[5]刘涛.精通ANSYS[M].北京:清华大学出版社,2002.[6]任重.ANSYS实用分析教程[M].北京:北京大学出版社,2003.Using ANSYS to analyze the vibration state of piezoelectric ceramicL EI Hui,ZHOU Shuang2e(School of Mathematics and Computer Science,Hubei University,Wuhan430062,China) Abstract:The application of piezoelect ric ceramic becomes more and more extensive,However,it needs to take much more time to st udy it by t raditional experiment s and t he result s are often unilateral for complex prezoelect ric ceramics boundary co ndition and st ress state.In t his view,we use big2scale general piezoelect ric ceramic’s boundary condition and st ress state simulating software ANS YS8.0to carry t hrough comp uter simulation.We gain t he pict ure of piezoelect ric ceramic’s vibration state by using comp uter to analyze t he coupling effect of piezoelect ric ceramic.The experiment result indicates t hat ANS YS8.0can competently deal wit h p roblems about piezoelectric coupling field.It optimizes t he time of creating model and comp uting largly,and it s st rong f unction of post dealing makes researchers can directly obtain data result s and simulating images.K ey w ords:simulation;piezoelect ric ceramic;vibratio n state(责任编辑肖铿,胡小洋)。