无网格精细积分算法在二维结构振动问题中的应用研究.pdfVIP

振动与冲击 OF AND 第26卷第5期 JOURNALVIBRATIONSHOCK + + 楚 十 十 甄 无网格一精细积分算法在二维结构振动问题中的应用研究 任传波，贺光宗 (山东理工大学交通与车辆工程学院，淄博255049) 摘 要 将无网格一精细积分法用于二维结构振动问题的求解，通过瞬时最小势能原理构造与弹性动力学方程等 效的能量泛函。由伽辽金无网格法在空间域内进行离散；在时间域上通过与Romberg积分相结合的精细积分法求解，得 到了二维结构的固有频率和振型以及在不同激励作用下的位移、应力和速度响应。 关键词：无网格法，精细积分，最小势能原理，二维结构，振动问题 中图分类号：TU311．3 文献标识码：A 弹性动力学问题可以表示为 其近似函数““(石) (1) or虮『+Z=pui+c厶i(在n内) orif忍ii=ti (在，。上) (2) J=l Ui=ui (在，。上) (3) (5) 初始条件： 其中：p(咒)是m维完备多项式基，常用线性基、二次基， a(算)是系数。 酩(算，。。)=酩。(引搿∈力 (4) 在局部范围内构造带权重的范数J(x) h(27，to)=凹o(搿) 其中P，c分别为材料的质量密度和阻尼系数； ui。=矿M／∥；uij越‘／0t 其中筇；是石的影响域内的n个节点。对式(6)取极 如上偏微分方程表示的动力学问题的求解在计算 值得 力学中占有重要地位，其解析解不容易得到，因此多采 a(戈)=A一1(戈)B(z)u+(7) 用数值解法。其中有限元法Hoo应用最多。本文即是 其中 将无网格和精细积分结合起来求解此问题的一种数值 算法。 A=∑wi(27)p’(石i)p(戈i) 无网格法是近年来发展比较迅速的一种求解边值 B=[Wi(戈)p(戈i)，…W。(戈)p(搿。)] 问题的有效方法，具有近似函数不依赖于网格，且计算 Wi(戈)=w(x一石i)W：(石)为权函数，具有紧支 精度高，前后处理简单的优点。一些学者也采用无网 性，本文选取指数型函数1为权函数。 格法对动力学问题进行了研究[3’4]。本文采用了精度 于是式(5)可以写为 较高的伽辽金无网格法[51在空间域上进行离散，即采 (8) 用移动最小二乘法构造形函数；采用罚因子引入本质 配“(菇)=∑ni(石)u。=J7v(戈)“+ 边界条件。由钟万勰教授提出的精细积分法【6j具有精 其中N(27)为形函数 度高、无条件稳定、自启动等优点。本文结合精细积分 (9) 的2Ⅳ类算法和Romberg积分在时间域内求解。算例 形函数的导数为 证明将无网格法和精细结合的算法具有较高的精确度 和灵活性。