空间杆系屈曲分析的非线性有限元法研究初探.pdfVIP

钢结构工程研究⑧ ‘钢结构》2010增刊 空间杆系屈曲分析的非线性 有限元法研究初探 钱若军1林智斌1袁行飞2 (1．同济大学，上海200092；2．浙江大学，杭州310058) 提要：对屈曲分析的研究，一部分学者致力于经典的屈曲理论研究，另一部分学者采用非线性有限单元 法对结构的屈曲进行数值分析。经典的屈曲理论推理严密，概念清晰，但是它无法解决复杂的结构系统， 特别是离散系统的屈曲分析；而传统的有限元法虽能求解结构系统，但由于其数学力学模型不足、理论支 持不够严密，并不能获得准确的求解。作者按照经典的屈曲理论采用有限单元法对简单的空间杆系结构的 屈曲进行了初步探讨，按照屈曲理论和有限单元法给出了空间杆系结构屈曲分析的有限元控制方程，对于 非线性有限元屈曲分析进行了初步探讨。 关键词：屈曲分析；非线性分析；势能函数；有限元； 自从欧拉(EuteO研究了著名的直杆屈曲问题以来，屈曲研究已有二百多年的历史。很 多学者在柱子和壳体的屈曲研究分析中奠定了现代屈曲研究的基础。冯。卡fq(VonKarman) 理论分析结果有较大差异现象给出了解释，他们认为差异是平衡点处的强非线性和初始缺 陷引起的。对一些简单结构如单杆和梁杆体系的屈曲问题进行了研究，得到了具有里程碑 式的结论，他们认为系统势能函数中应该包括未知量的高阶量即三阶以上的高阶量，继他 价值直到1960才被人们发现，Koite理论仍被称为屈曲理论研究中的经典理论。对屈曲问 题的研究有两部分学者分别进行研究，一部分学者致力于届曲理论研究，而上世纪七八十 年代以来，另一部分学者采用非线性有限单元法对结构的屈曲进行数值分析。虽然他们的 论文中也引述了经典理论的概念，但是所采用的方法都是非线性有限元的增量和自动增量 过程，这部分的研究结果比较大的影响了结构设计中的屈曲分析。然而，前一部分研究经 典理论的学者认为有限元法无法用于屈曲分析。 本文按照经典的屈曲理论采用有限单元法对简单的空间杆系结构的屈曲进行了研究和 探讨，由于屈曲理论涉及到结构、数学等很多领域，所以本文的研究是初步的，而且是概 念性的，希望以此纠正作者本人在屈曲分析研究中错误和模糊的概念。 系统势能函数的构造 按照经典的屈曲理论，单参数系统的势能函数应包括位移三阶以上的高阶量，即 525 钢结构工程研究⑧ 《钢结构》2010增刊 V=V(u，“2，“3，…，名) 其中，“为系统的位移；五为系统的荷载参数，这里仅考虑单参数。 对于结构体系，单元的势能函数为： ^■ y=毒，eodv-uTp 一V 式中， 占为单元的应变； 盯为单元的应力；P为作用于单元上的外荷载。 单元的应变可分为线性和非线性两部分，即 ￡2sL+￡NL 其中， 气，‰分别为单元应变的线性和非线性部分。 对于空间杆系 ，， 、z 吒2瓦Ou，‰2j1I瓦OuJ’… 传统的有限元平衡方程 在一般的结构分析中，总是忽略非线性部分中位移二阶以上的高阶量，而对于非线性 应变引入迭代格式从而使得非线性部分简化为线性，这样当对式(2)引入物理关系，按照能 量原理，可得平衡方程 8V=0 (5) 经过整理以后可得系统的有限元基本方程 (E+E)衄=2p-f。 (6) 式中， 置，为单元的线弹性刚度矩阵，它表示了材料对单元