2 程序设计与软件应用.ppt

第2章 有限元程序设计方法 2.1 程序基本框图 1、输入基本数据（结构描述）： （1）控制数据：如结点总数、单元总数、约束条件总数等； （2）结点数据：如结点编号、结点坐标、约束条件等； （3）单元数据：如单元编号、单元结点序号、单元的材料特性、几何特性等； （4）载荷数据：包括集中载荷、分布载荷等。 2、单元分析 （1）各单元的bi,ci(i,j,m) ， 面积A； （2）应变矩阵[B]，应力矩阵[S]； （3）单元刚度矩阵[k]； （4）单元等价载荷列向量[F]。 总刚存贮 全矩阵存贮法：不利于节省计算机的存贮空间，很少采用。K[i,j] 对称三角存贮法：存贮上三角或下三角元素。 半带宽存贮法 ：存贮上三角形（或下三角形）半带宽以内的元素 。 一维压缩存贮法 ：半带宽存贮中仍包含了许多零元素。存贮每一行的第一个非零元素到主对角线元素。 例：计算下图半带宽。 4、引入约束条件 手算时采用去行列法 ，而计算机编程时采用乘大数法。即：指定结点位移对应的主对角元素乘上一个大数 ，同时将{P}中对应元素换为结点位移指定值与扩大了的主对角线元素的乘积 。 5、线性方程组求解 求解方法常用：GAUSS消元法，QR分解法等。其程序在一些专著中列出（例如见：徐士良编。FORTRAN常用算法程序集。清华大学出版社，1992）。在此不作详细介绍，其方法参阅[数值分析]有关书籍。 6、单元应力