梁的有限元分析原理.ppt

* * 3-D Elastic Beam BEAM4 is a uniaxial element with tension, compression, torsion, and bending capabilities. six degrees of freedom at each node BEAM24 3-D Thin-walled Beam The element has plastic, creep, and swelling capabilities in the axial direction as well as a user-defined cross-section. This element allows a different unsymmetrical geometry at each end and permits the end nodes to be offset from the centroidal axis of the beam BEAM44 3-D Elastic Tapered Unsymmetric Beam * * * BEAM188 3-D Linear Finite Strain Beam BEAM188 is suitable for analyzing slender to moderately stubby/thick beam structures. This element is based on Timoshenko beam theory. Shear deformation effects are included. This element is well-suited for linear, large rotation, and/or large strain nonlinear applications. * * BEAM189 3-D Quadratic Finite Strain Beam BEAM189 is a quadratic (3-node) beam element in 3-D. For a description of the low-order beam, see BEAM188. * * 有限元程序设计方法简介 程序基本框图 1、输入基本数据（结构描述）： （1）控制数据：如结点总数、单元总数、约束条件总数等； （2）结点数据：如结点编号、结点坐标、约束条件等； （3）单元数据：如单元编号、单元结点序号、单元的材料特性、几何特性等； （4）载荷数据：包括集中载荷、分布载荷等。 开始 输入基本数据 计算单元刚度矩阵 形成总体刚度矩阵 形成结点荷载向量 引入约束条件 求解方程组，输出结点位移 计算单元应力，输出结果 结束 * 2、单元分析 （1）各单元的bi,ci(i,j,m) ， 面积A； （2）应变矩阵[B]，应力矩阵[S]； （3）单元刚度矩阵[k]； （4）单元等价载荷列向量[F]。 开始 输入基本数据 计算单元刚度矩阵 形成总体刚度矩阵 形成结点荷载向量 引入约束条件 求解方程组，输出结点位移 计算单元应力，输出结果 结束 3、系统分析 （1）整体刚度矩阵[K]的组装； （2）整体载荷列阵{P}的形成； [K]的存储；约束引入；求解 * 总刚存贮 全矩阵存贮法：不利于节省计算机的存贮空间，很少采用。K[i,j] 对称三角存贮法：存贮上三角或下三角元素。 半带宽存贮法 ：存贮上三角形（或下三角形）半带宽以内的元素 。 一维压缩存贮法 ：半带宽存贮中仍包含了许多零元素。存贮每一行的第一个非零元素到主对角线元素。 * 等带宽形式 UBW UBW 行 号 1 → IR → N→ 1 列 号 JC 行 号 1 → IR → N→ 1 JC-(IR-1) 方阵形式 （1）半带宽存贮法 * 方阵存贮和半带宽存贮地址关系 存贮方式 行号 列号 方阵存贮 IR JC 等带宽存贮 IR JC-IR+1 半带宽计算：设结构单元网格中相邻结点编号的最大差值是d，则最大半带宽为UBW： 结点编号：欲使最大半带宽UBW最小，必须注意结点编号方法，使直接联系的相邻节点的最大点号差最小。 * 举例 B = 2(4-1+1) = 8 B = 2(6-1+1) = 12 Advantages of 2D Storage 1)Space-saving; 2)Easy to be computerized Disadvantages of 2D Storage Enormous storage is required when local bandwidth is large. * 例：计算下图半带宽。 结点数N=91，总刚