框架结构屈曲的精确有限元求解.pdfVIP

第 41卷 第 6 期 力 学 学 报 Vo1．41．NO．6 2009 年 11月 ChineseJournalofTheoreticalandAppliedM echanics NOV．，2009 框架结构屈曲的精确有限元求解 陈太聪 马海涛 。) (华南理工大学土木工程系亚热带建筑国家重点实验室，广州 510640) 摘要 基于屈曲微分控制方程的一般解，构造了Euler梁在轴力作用下的精确形函数，建立了用于框架结构屈 曲分析的精确有限单元，得到了单元刚度矩阵和几何刚度矩阵的显式表达，并提出了基于常规特征值计算的迭 代算法 以确定屈曲载荷及相应失稳模态的精确解．研究表明，对于线性稳定性分析而言，常规框架有限单元可 视为精确有限单元的 一种近似．若采用精确单元，无需进行网格细分就可以获得精确的屈曲载荷和失稳模态． 数值算例证明了新单元和算法 的效率和精度． 关键词 有限元法，屈曲分析，框架结构，精确单元，形函数 中图分类号： 034 文献标识码： A 文章编号： 0459—1879(2009)06—0953—08 引 言 法建立杆件屈曲后的柔度方程，再通过矩阵求逆得到 单元刚度方程，然后通过求解超越方程来确定结构的 对 于框架结构 的线性静力分析， Tong[1j和 屈 曲载荷 (例如，见文献 [7])．近年来，Eisenberger[， Kanok—Nukulchai等 _2_已给出证明，若有 限单元的 Kim 等 [。]和 Leung[加』通过微分控制方程 的幂级数 位移形函数能满足Euler控制方程，则无需网格细化 位移解构造了压杆屈 曲分析的近似精确有限单元， 即可解得精确的节点位移．就线性静力分析而言， 结构屈 曲载荷的求解也需要采用异于常规特征值求 常规的三次多项式形函数恰好满足 Euler梁弯 曲问 解的特殊算法进行．由于特殊求解算法的应用，上 题的控制方程，因此可以给出相应的精确单元刚度 述两类方法在通用有 限元软件中不易实现，针对由 矩阵，对于考虑横向剪切变形的Timosheko梁，则 多个杆件构成的框架结构的屈 曲分析也不易进行． 可以采取多种不同方法获得用于线性静力分析的单 Ma[nJ针对弹性杆的轴 向振动分析提出了一种 元精确刚度矩 阵，其中几种方法 的具体过程可见文 新型的有限单元和振型分析的精确算法．单元的位 献 [2～5]_ 移形函数根据轴向振动微分控制方程 的位移通解形 进行线性稳定性分析时，通常的做法是采用与 式确定，由此构造得到了精确的单元刚度矩阵和质 线性静力分析时相同的位移形函数建立几何刚度矩 量矩阵，并提出了一种基于常规特征值解法精确计 阵，然后应用这样得到的几何刚度矩阵和线性静力 算 自振频率和相应振动模态的迭代算法．由于该种 分析所用的单元刚度矩 阵建立特征值方程，最后通 单元具备明确的位移形函数、刚度矩阵及质量矩阵， 过求解广义特征值问题确定屈曲载荷和对应失稳模 且 自振频率的求解基于标准的特征值求解算法，因 态 _6J．在采用这类常规单元进行屈曲分析时，往往 此可以方便地在通用有限元软件中实现对 自振频率 必须对线性静力分析采用 的网格进行细分才能得到 和相应振动模态的精确求解．最近，Ma[J又将这 较准确的结果，而且对高阶模态通常需要更细密的 一 方法推广到框架结构的动力分析． 格划分．