基于Matlab的大型梁结构固有振动边界元法仿真.pdfVIP

第15卷第2期 湖南工程学院学报 Vbl．15．Nb．2 2005年6月 JoumalofHunanInstituteof June．2005 Engineering 吴莉华 (湖南工程学院机械电子工程系，湖南湘潭411101) 摘要：简要分析了边界元法在梁类结构振动仿真研究中的优点和缺点，应用边界元法对矩形梁弯曲 振动进行了固有振动仿真分析，当梁的弹性常数改变时，结构的固有振动频率会相应地改变，为通过改 变结构材料参数来改善梁结构的动力学性能进行理论探索． 关键词：梁；固有振动；边界元法；仿真 中图分类号：THl22文献标识码：A 系，即可求出领域Q内的内点值． 0前言 1边界元法的特点 Element 边界元法(BEM，BoundaryMeihod)作 为边界积分方程的离散数值法有其特有的优点，如 与有限元法相比，边界元法有如下特点： 适合于求解无穷域数值问题，具有降维特点等，比如 (1)据准备工作量小、解析时间短、计算费用低； 可以将三维问题降为二维问题，将二维问题降为一 (2)析计算精度高； Element 维问题．边界元法与有限元法(FEM，Finite (3)可以实现降维，如一维问题工程问题可变成 MethOd)不同，它只在对象体的边界上进行离散，分零维问题，其边界就是点；二维、三维工程问题可以 割单元，属于边界型数值解析方法．对于三维振动问 分别降为一维和二维问题，可见边界元法使解析问 题尤其是长宽比很大的梁类零件的振动问题，采用 题降了一维； 边界元法较为简便．对长宽比很大的仿真模型， (4)处理一些特殊问题非常方便，如对长宽比很 FEM仿真精度很低，甚至因发散而无法实现仿真分 大的三维仿真模型，边界元法能够克服有限元法的 析，BEM则能够克服FEM的这类缺点．从解析原理缺点，适应这类仿真模型的特点和要求，顺利完成仿 和建模过程来看，边界元法实际上是一种数理解析 真分析； 。 和数值解析相结合的解析方法，其求解过程包括数 (5)边界元法建模较复杂，需进行数理解析；与 理数值解析和方程求解． 有限元法相比，多了数理解析和求基本解的工作． (1)数理解析：首先应用三维拉普拉斯基本解和 (6)联立方程求解麻烦； 积分原理，建立领域与边界的关系式即积分方程式； (7)现有应用软件功能尚不够强大，没有现成的 然后将领域问题转化为边界问题，建立边界上的关 商品化软件可供使用． 系式即边界积分方程式． 可见，在长宽比很大的梁结构的振动问题求解 (2)数值解析：边界积分方程式是在边界上进行 中，应用边界元法较好．本文在分析了大型梁结构的 积分，所以边界元法建模只需要在边界上进行离散， 将边界分割成单元，写出边界节点方程式，再用求和 序，并对不同弹性常数的大型梁结构进行了边界元 的方法写出边界联立方程式． 法固有振动仿真和对比分析，揭示出这类结构的固 (3)求解：对边界联立方程式求解，即可求出对 有振动特征与弹性常数的内存关系，为通过改变结 象体边界上的节点位移；然后根据领域与边界的关 构的材料