第2章 杆系结构有限元静力分析程序系统
2.1 概 述
根据第1章的内容,不难掌握杆系有限元结构分析的基本过程,接下来的工作是如何编制计算机程序,完成杆系结构的静力分析。通常,程序的具体结构是多种多样的。它取决于预计的解题规模、程序的通用程度以及计算机的硬件和软件系统购配置状况。因而,并不存在某种**的*优程序结构。另一方面,各种程序系统的基本结构在实质上并无显著的差别。因而,在这里,只着重对程序的核心部分予以简单的分析,希望使读者掌握基本的数值解法和程序设计能力,而不追求程序细节的优化。 线性结构静力分析计算程序的基本内容和步骤可表示为如图 2.1所示的粗略框图。在框图中,矩阵K指的是受约束结构的总刚度矩阵,δ是独立结点位移列矢量,P为与δ相应的结点荷载总矢量。
在整个计算中,对程序的效率影响*大的是结构刚度矩阵的存储方式和组集方法以及求解线性方程组的算法。这部分内容在读者需要编者计算机程序时可参阅结构有限元的相关资料。
在结构刚度矩阵的存储方式和组集方法以及求解线性方程组的算法中,通常采用的技术有以下几点。
(1)结构刚度矩阵是具有高度稀疏性的对称正定矩阵。一般,总刚度矩阵的元素数量是很大的,采用压缩存储形式是**必要的。等带宽存储和变带宽存储这两种*常用的存储方式,它们都适用于线性方程组的直接解法,也适用于迭代解法。
(2)无论是等带宽存储还是变带宽存储,都可利用结构刚度矩阵的对称性,只存储刚度矩阵的上三角部分或下三角部分。
……