热机电耦合智能板结构的一类有限元模型.pdfVIP

维普资讯 第25卷 第2期 应 用 力 学 学 报 Vo1．25 No．2 2008年 6月 CHINESEJoURNALOFAPPLIEDMECHANICS Jun．2008 文章编号：1000—4939(2008)02-0263—05 热机 电耦合智能板结构的一类有限元模型 王小兵 陈建军 谢永强 梁震涛 陈永琴 (西安电子科技大学机电工程学院 710071 西安) 摘要：针对热机 电压电薄板结构，提出了一种包含4个位移节点、2个电势节点和 8个温度节点的 有限元模型，其位移场按平板壳单元描述，而电势场和温度场则按线性插值方法描述，并基于虚功 原理导出了热机 电有限元方程，通过比较，表明所得有限元方程的正确性。最后以悬臂板为例，对 位移响应、温度响应和输出电压进行了数值仿真，凸现了在小挠度变形下，变化幅值较大的温度对 输 出电压 的重大影响。 关键词：热一机一电耦合；有限单元法；计算机仿真 中图分类号：0242．21 文献标识码：A 的三阶剪切变形理论，而其电势场和温度场都采用 1 引 言 双线性拉格朗日插值函数加以描述，该建模方法也 被HeungSooKim[4]应用到了曲面壳结构中，相比 压电智能板壳结构在热一机耦合研究上的蓬勃 第一类模型而言，该模型较为简洁，但其不易同时满 发展带动了热机电耦合的相关研究，而其有限元方程 足复杂的三维热边界条件。 推导的一般原理和方法在早期文献[1—2]中就已比较 上述各个模型的共同点是认为整个压电层表面 成熟，并被广泛地运用到了不同的有限元模型中。 上各点电势是连续分布的，然而实际使用中，整个压 利用有限元法分别描述位移场 ，电势场和温度 电层是由多个压电片铺设而成，而且每一个压电片 场，就需要有相应的位移节点，电势节点和温度节 的上下表面都镀有一层金属电极[5]，使得整个压电 点，学者们一般偏向于将所设立的节点同时作为位 层表面上各点电势是台阶式离散分布的，在文E6]中 移节点、电势节点和温度节点，实践表明，这种建模 的图6以及文E7]中的图3都体现了这种电势连续 方式有利于复杂结构的有限元编程实现。对于压电 分布与离散分布这二者间的明显差异。这时，文[2— 薄板结构，目前有两类比较典型的节点共用的模型， 3]等使用的节点共用单元模型就不能直接加以使 第一类为位移节点、电势节点和温度节点共用体单 用，从而有必要研究在单元中分别设立位移节点，电 元的8个节点，第二类为共用中性面上四边形单元 势节点和温度节点的有限元模型。由于有诸多成熟 的4个节点。例如 H．S．TzouE]采用的就是第一类 的位移单元模型、电势单元模型、单纯的热分析单元 模型，其好处为便于描述复杂热边界条件下的三维 模型可供选择，并能利用能量原理将它们有机地结 温度场，但该类有限元模型对应的状态方程阶数较 合起来，可使这种建模方式具有更多的自由性、多样 高，而且用于薄板结构日寸需引人位移非协调项；而 性和灵活性，然而，由于该类模型从节点编号到有限 HaozhongGu[。]J针对仅在板上下表面施加热载荷 元编程实现要麻烦得多，使得至今该类模型仍很鲜 下使用的模型就是属于第二类 ，其位移场采用薄板 见。本文就针对压电薄板结构，扬长避短地选用了 *基金项目：国防预研基金(51421060505DZ0155)；陕西省自然科学基金 (2005A009) 来稿 日期：2006-09