变截面压杆稳定计算的有限元法.pdfVIP

第 3 期 (总第 108 期) 山　西　机　械 N o 3 2000 年 9 月 SHAN X I　M A CH IN ER Y Sep 文章编号:(2000) 0300 1602 变截面压杆稳定计算的有限元法 1 1 2 戴保东 , 郑　荣 , 吉子轩 ( 1 太原重型机械学院, 山西　太原　030024; 2 太原冶金工业学校, 山西　太原　030002) 摘要: 利用有限元法对变截面压杆进行弹性稳定性分析, 可以计算变截面压杆的临界荷载及相应的计算长度系 数。计算结果表明, 这种方法具有计算简便、便于计算机处理等特点, 是对经典理论成果的补充与完善。 关　键　词: 临界荷载; 有限单元法; 屈曲 中图分类号: TB 115　　　文献标识码: A 引　言 但是变截面压弯杆单元的刚度矩阵没有统一的形式, 变截面压杆由于其节省材料且具有良好的受力特 不便于计算机编程。本文在此采用近似方法, 即将变 性, 故在起重机械、建筑结构中被大量使用。但是与 截面杆单元近似用一个以其中点截面为准的等截面杆 之相关的稳定性分析理论却极不完善, 这主要是由于 单元代替。显然, 这样处理会带来一定的误差, 但是 在理论上推导变截面压杆的稳定方程是很繁难的, 而 只要增加单元数 目, 就可以大大降低这种近似所引起 且对应的微分方程又是超越方程, 要想进一步求解这 的误差。 些方程往往要遇到数学上的困难。从现有的文献资料 在 由单元刚度矩阵集结形成原始总体刚度矩阵 和设计规范来看, 传统设计方法在引用经典理论成果 后, 还应引入边界条件对原始总体刚度矩阵作进一步 对这类问题进行计算时, 尚有疑难之处, 因此对这一 的修改方可求解。因为本文只考虑求解临界荷载, 不 类问题进行深入研究是很有必要的。本文将有限元法 计算结点位移, 所以在此可以采用 “划零置一”的方 引入变截面压杆的稳定计算中, 应用计算机求解, 获 法对原始总刚度矩阵进行修改, 即将原始总刚度矩阵 ( ) 得了令人满意的结果。 [ ] 和 [ ] 中与边界位移为零 如 = 0 相对应 KE KG j 1　基本原理 的元素 j j 换为 1, 而将对应的第 行、第 列的其它