平面压电结构的径向基函数无网格法.pdfVIP

第24卷第4期 苏 州 大 学 学 报 (自然科学版) V0L24 No．4 2008年 10月 JOURNALOFSUZHOUUNIVERSITY(NATURALSCIENCEEDITION) 0ct．2008 平面压电结构的径向基函数无网格法 刘 芳，姚林泉 (苏州大学数学科学学院，江苏 苏州215006) 摘 要：利用压电材料的正、逆压电效应制成的传感器和驱动器在很多领域有广泛的应用．压电结构的机电耦 合效应，给问题的求解带来了一定的困难．本文利用无网格法对压电平面问题的控制方程进行求解．利用径向 基函数进行插值近似，直接配点法对控制方程进行离散，得到无网格离散的线性控制方程组，最后通过数值计 算与分析，验证了方法的可行性以及不同径向基函数对结果的影响． 关键词：径向基函数；配点型无网格法；压 电结构 中图分类号：0343 文献标识码：A 文章编号：1000—2073(2008)04—0018—06 0 引 言 无网格法是近年来发展起来的一种新的数值计算方法．它采用基于点的近似，可以彻底或部分地消除网 格，不需要网格的初始划分和重构，不仅可以保证计算的精度，而且可以减少计算的难度．它克服了有限元法 中网格划分带来的缺点和不便，实际上是对有限元法的补充．无网格法的基本过程与有限元法非常相似，主要 有区域离散、形函数构造、方程离散、求解代数方程以及后处理等步骤．形成形函数的方法有：移动最小二乘法 (MLS)、重构核近似法(RKM)、径向基函数法(RBF) 等． 径向基函数法，是一类以计算点 到节点 ，的距离为 自变量的函数，具有形式简单、各向同性等优点，因 此能成功应用在数值计算中．而具有紧支特性的径向基函数(CSRBF)，具有带状稀疏系数矩阵的特点，适用 于求解大型科学与工程问题．本文将径向基函数引入配点法中用于分析压电结构． 1 压电控制方程 考虑平面压电材料，建立 坐标系．设压电材料的极化方向沿z轴，．z平面下二维压电材料的连续性方 程为[23 (i)力学平衡方程 ： + …+ =0，丁… + + =0． (1．1) (ii)电学平衡方程：D +D =0． (1．2) CI1 C13 0 0 一e3l C13 C33 0 0 一e33 0 0 。55 一e15 0 (1．3) 0 0 e15 l 0 e31 e33 0 0 螽 ㈣ 几何方程 ： =u ；= =1 ， + u， )． (1．4) 电场与电势的关系：E =一咖 E：=一咖 (1．5) 力边界条件： n+7-n：一 ：0，Ixzn+ n：一 =0． (I．6) 收稿 日期：2007—08—15 基金项 目：国家自然科学基金，江苏省 自然科学基金 (BK2006725)，暨南大学重点实验室开放课题基金资助项 目 作者简介：刘 芳 (1983一)，女，江苏泰兴人，硕士研究生，主要从事计算力学的研究． 第4期 刘 芳，姚林泉：平面压电结构的径向基函数无网格法 19 位移、电势和电场边界条件：=。，= 。，=咖。，咖 =