杆系结构动力问题的直接刚度法董萎： 奎兰丹j函数矩阵；运动频率.PDF

维普资讯 ’ ￡ 3’弓 一 №2 』 西北轻工业学院学报 Jun．1995 Vo1．13 JOURNALOFNORTHWEST INSTITUTEOFLIGHT INDUSTRY ·35 · 7 ／／3 杆系结构动力问题的直接刚度法 土盖 (机械工程磊) ^ 摘 要 ／ 本文从构造含有运动频率 的形函数矩阵入手．建立了杆糸结构动力问题的有 艰元法 经验证，该方法所计算的固有频率有效率高，计算的精度也很高 董萎：奎兰 丹j函数矩阵；运动频率；有限元；壹丝堕堙蓝 t象结≠ 叹 中图法分类号：。。。·。 名 1 引 言 机械工程中，高精度旋转、高速运动的机械要求精确的动力设计。机械中像轴类零件、刚 架的动力计算，连杆机构KED分析等大都采用杆系结构动力学有限元法 一 众所周知 在结构动力学的有限元分析中．单元的形函数矩阵 ]与静力分析中的[Ⅳ] 相同。它的各元素只是点的位置坐标的函数，并依此来建立刚阵[]和质阵[ ]。对于静力 学问题，[Ⅳ]是由静态位移的分布确定的．单元的静态位移场是唯一的。但是，对作振动的结 构单元，其动态位移场与静态位移场不同．与结构频率 相关。因此，用一般的动力学有限 元法求解动力学问题 ．使其计算精度受到影响。 如何用有限元法提高计算精度呢?当考虑频率 影响时，在文献口]中相应的刚阵和质 阵分别为 []一 [o]4- [ ．]+ … [Ⅳ]一 [Ⅳo]4- [ ]+ … 上式 由于含有高次的动力修正项，给计算带来很大的困难。 文献 [2]介绍的动力直接刚度法，在导出动力刚度矩阵之后，均按静力法一样的过程进 行分析。但是，在 自振频率分析中，因动力刚度系数是运动频率的函数，必须在分析过程中反 复调整，影响计算精度的提高，计算量也比较大。 本文从建立含有频率m的形函数矩阵[Ⅳ]人手 ．采用与常规的动力有限元分析大体相 同的方法和步骤，所形成的杆系结构动力 问题的有限元法，试图弥补上述方法中的不足之 力学 处。 2 单元形函数矩阵 2．1 梁元形函数矩 阵 设梁单元为等截面、匀质的直梁，不计剪切变形和转动惯量的响应。梁跨间不承受载荷 其 自由振动微分方程为： ·收稿 日期：l994一l2一l3 维普资讯 · 36 · 西北轻工业学院学报 第 13卷 ，+拿一0 式中E1为抗弯刚度， 为质量分布集度。 利用分离度量法，解得；并且用三角函数和双曲函数等价地替换指数函数后为 )一 A sinax + COSd。4-cshd。4-D ch4 (1) 式中a‘一 C百orm ， 待定常数A，B,C，D决定梁振动的形状和振幅。 由静力有限元法可知，梁元的位移模式的矩阵形式 0)= […N N N ]{}