基于ODE求解器的高层建筑筒体结构动力特性分析

论文摘要

筒体结构在现代建筑中发挥着重要的作用,它不仅具有好的强度和刚度,且造型美观、使用灵活、受力合理,整体性强,因此在高层建筑中被广泛采用。本文通过有限条（板）元技术理论对高层建筑筒体结构进行竖直方向上的单元划分,取各条元结线的位移函数为基本未知函数,并由最小势能原理建立高层建筑筒体结构动力特性的常微分方程组及相关的边界条件,然后对方程采取一系列的数值变换技巧,将其转换成标准的常微分方程组特征值问题,以便可以直接调用ODE 求解器进行求解。ODE 求解器——COLSYS 具有通用性,可自适应求解,可满足用户预先对解答精度所指定的误差限。其中,变换技巧主要有逆幂迭代法,它具有很好的放大低阶、消减高阶的作用;特征值移位,可以极大地加速逆幂迭代的收敛速度;正交化,即从用户提供的初始解中除去前i-1 个特征成分,余下的再作为逆幂迭代的初始解,其放在逆幂迭代过程中进行。对筒体结构而言,其结构沿高度方向一般为等截面或分段变截面的,有限元线法理论更适用于这种结构。本文在对筒体结构进行动力特性分析的过程中,主要进行了下面的研究: 1) 通过力学和数值分析的方法,实现了实际结构到力学模型的转化; 2) 通过有限条元技术和最小势能原理,完成了力学模型到数学模型的转化; 3) 应用数学变换技巧,将常微分方程组转换成标准的非线性ODE 问题; 4) 应用计算机语言,达到数学模型到计算机模型的转化,完成本课题的研究。在本文的分析中,采用了分布质量的板条单元,与以往的的集中质量模型不同,从算例分析表明,与常规有限元法相比,本文方法具有两大特点:一是采用了有限条单元,单元数少,划分单元直观;二是采用了ODE 求解器解法,该方法是数值解析法,求解精度高。

论文目录

• 中文摘要
• 英文摘要
• 1 绪论
• 1.1 高层筒体结构的发展概况
• 1.2 筒体结构动力特性求解方法比较
• 1.3 0DE 求解器
• 1.4 有限元线法概述
• 1.5 课题的研究目的、主要内容及方案
• 1.5.1 课题的研究目的
• 1.5.2 课题的主要内容
• 1.5.3 课题的研究方案
• 2 筒体结构的有限条模型
• 2.1 筒体结构与有限元线法
• 2.2 筒体结构基本假定
• 2.3 筒体结构的等效弹性参数的处理
• 2.3.1 实体单元
• 2.3.2 框筒单元
• 2.3.3 加肋框条
• 2.3.4 连梁虚条
• 2.4 条单元分析
• 2.4.1 位移函数
• 2.4.2 插值函数的取法
• 2.4.3 正交异性板的受力分析
• 2.4.4 条元应变能
• 2.4.5 条元的惯性力势能
• 2.4.6 坐标转换
• 2.5 楼板的影响
• 2.5.1 楼板刚性条件的引入
• 2.5.2 楼板质量矩阵
• 2.6 角柱势能
• 2.7 整体分析
• 2.7.1 等截面筒体结构动力特性
• 2.7.2 变截面筒体结构动力特性
• 3 高层建筑结构动力特性的求解器解法
• 3.1 常微分方程边值问题通用程序——COLSYS 简介
• 3.1.1 求解问题与范围
• 3.1.2 方法与功能
• 3.1.3 COLSYS 程序的使用
• 3.2 常微分方程特征值问题的求解
• 3.2.1 标准非线性问题的转换
• 3.2.2 逆幂迭代法
• 3.2.3 正交化
• 3.2.4 特征值移位
• 3.2.5 算法流程
• 3.3 高层建筑结构动力特性的求解器解法
• 3.3.1 等截面筒体结构振动控制微分方程的变换
• 3.3.2 变截面筒体结构振动控制微分方程的变换
• 4 程序的设计与使用
• 4.1 程序简介
• 4.1.1 程序设计
• 4.1.2 程序的功能及使用说明
• 4.2 程序框图
• 4.2.1 程序总流程图
• 4.2.2 核心子程序getsolutn 流程图
• 4.2.3 主要子程序简介
• 5 算例分析
• 5.1 算例
• 5.1.1 算例说明
• 5.1.2 算例计算结果及分析
• 5.1.3 ANSYS 计算结果
• 5.1.4 算例小结
• 5.2 算例2
• 5.2.1 算例说明
• 5.2.2 算例计算结果及分析
• 5.2.3 ANSYS 计算结果
• 5.2.4 算例小结
• 5.3 本章小结
• 6 结论与展望
• 6.1 结论
• 6.2 展望
• 致谢
• 参考文献
• 附录
• 独创性声明
• 学位论文版权使用授权书