凯威特型单层球面网壳的特征值屈曲特性分析.pdfVIP

凯威特型单层球面网壳的特征值屈曲特性分析 胡瑞，张宁宁 辽宁工程技术大学研究生院，辽宁阜新（123000 ） E-mail：kaiyuanrain@163.com 摘 要：本文主要以施威德勒型单层球面网壳为研究对象,对其进行特征值屈曲分析，特征 值屈曲分析用于预测一个理想弹性结构的理论屈曲强度，也即弹性屈曲分析方法。由于初始 缺陷和非线性使得很多实际结构的屈曲行为不是在弹性屈曲强度处发生，所以特征值屈曲分 析的结构过于保守。经典的屈曲分析是采用特征值屈曲分析法，它适用于对一个理想弹性结 构的理想屈曲强度（歧点）进行预测，主要是使用特征值公式计算造成结构负刚度的应力刚 度阵的比例因子。特征值屈曲主要特点是在结构未变形位置建立结构总体弹性刚度阵和几何 刚度阵，最后把稳定分析转化为求解矩阵广义特征值问题。 关键词：特征值屈曲分析，稳定，屈曲，变形 中图分类号：TU688.6 1. 引言 近年来, 网壳结构发展迅速、形式多样，网壳结构在大跨度建筑中已越来越多地被利用, 而且具有广阔的发展前景。网壳结构在达到屈曲载荷之前其位移—变形曲线表现出线性关 系，达到屈曲以后曲线将跟随另外的路径[1] 。网壳结构失稳有两种性质根本不同的失稳形式， 即分枝点失稳和极值点失稳，对应着两种分析方法:第 1 类为线弹性最小特征值屈曲问题， 用于确定一个理想弹性结构的理论屈曲强度;第 2 类为极点问题，即考虑了结构几何非线性 和材料非线性情况下的极限承载力问题[2] 。第1 类稳定问题无论在理论分析中还是在工程应 用上都占有重要的地位，这不仅是由于第1 类稳定问题的力学概念比较明确，在数学上可归 结为特征值问题，因而求解相对较容易，更重要的是因为第1 类和第2 类稳定问题有着良好 的相关性，往往代表着第2 类稳定问题的上限，所以工程中通常以第 1 类稳定问题的计算结 果作为设计的依据。 2. 结构特征值屈曲分析 特征值屈曲分析方法是包括材料非线性、大变形等非线性因素的静力分析法，计算过程 可以一直进行到结构的限制载荷或最大载荷[3] 。对网壳这种大型空间结构, 稳定性能研究有 着重要的意义,稳定问题一般分为三类，第一类失稳是理想化情况，即达到某个荷载时，除 结构原来的平衡状态可能存在外，出现第二个平衡状态，所以又称平衡分岔失稳或分枝点失 稳，而数学处理上是求解特征值问题，故又称特征值屈曲。第二类失稳是结构失稳时，变形 将大大发展，而不会出现新的变形形式，即平衡状态不发生质变，也称极值点失稳。还有一 种跳跃失稳是当荷载达到某值时，结构平衡状态发生一明显的跳跃(运动) ，突然过渡到非邻 近的另一具有较大位移的平衡状态，可称为第三类失稳[4] 。 第一类稳定问题在有限元分析中为特征值屈曲(Eigenvalue buckling)，其基本步骤为：建 立几何模型或有限元模型；打开应力效应选项，以获得屈曲分析用的应力刚度，进行静力求 解；设置屈曲模态数目和扩展选项，进行特征值屈曲求解；观察结果，得到各个屈曲模态的 解[5] 。需要注意的是特征值屈曲为线弹性，对所定义的材料非线性将被忽略。 - 1 - 3. 数值模拟分析 3.1 工程概况 为了精确模拟实际结构，该模型采用空间三维模型。本文采用一个施威德勒型网壳结构， 其俯视平面形状为圆形，低平面的直径为80m，球面直径为80m，矢高为5.3m，球面中心角 为60度，具体模型前视图、底视图，见图1，图2 。 图1 模型前视图 图2 模型后视图 Fig.1 Model front view Fig.2 Model back view 3.2 数值模拟结果及分析 网壳结构为— 多自由度系统, 将网壳简化为空间梁体系。 (1) 运动方程 网壳结构体系的振动方程为[6]