有限单元法-杆.pptVIP

杆件系统的有限元法（九） 例（续） * 杆件系统问题的求解步骤： 1. 对杆件系统划分单元，对单元和节点进行编号； 2. 写出各个单元在整体坐标系下的刚度矩阵（单刚）； 3. 根据各单元刚度矩阵根据节点关系，装配成整体刚度矩阵； 5. 写出系统的边界条件，代入系统方程，求出各节点位移； 杆件系统的有限元法（十） 4. 根据力-刚度-位移的关系，写出系统方程； 6. 求解个节点支反力、各单元应力和应变等未知量。 * 整体刚度矩阵的集成 整体刚度矩阵是由在整体坐标系下，矩阵按照结点编号的顺序组成的行和列的原则，将全部单元刚度矩阵扩展成n×n方阵后对号入座叠加得到。 对于单元1 对于单元2 对于单元3 单元刚度矩阵集成得出整体刚度矩阵 杆件系统的有限元法（十一） * 杆件系统的有限元法（十二） 整体刚度矩阵的性质 整体刚度矩阵 中位于主对角线上的子块 ，称为主子块，其余 为副子块。 a. 中主子块 由结点i的各相关单元的主子块扩展之后叠加求得，即 b. 当结点i、 j为单元e的相关结点时， 中副子块 为该单元e相应的副子块，即 。 c. 当结点i、 j为非相关结点时， 中副子块 为零子块，即 。 d. 仅与各单元的几何特性、材料特性，即A、I、l、E等因素有关。 e. 为对称方阵， f. 为奇异矩阵，其逆矩阵不存在，因为建立整体刚度矩阵时没有考虑结构的边界约束条件。 * g. 为稀疏矩阵，整体刚度矩阵中的非零元素分布区域的宽度与结点编号有关，非零元素分布在以对角线为中心的带状区域内，称为带状分布规律，见图3-10(a)。在包括对角线元素在内的区域中，每行所具有的元素个数叫做把半带宽，以d表示。 最大半带宽等于相邻结点号的最大差值加 1 与结点自由度数的乘积，结点号差越大半带宽也就越大。计算机以半带宽方式存储，见图3-10(b)。半带宽越窄，计算机的存储量就越少，而且可以大幅度减少求解方程所需的运算次数。其效果对大型结构显得尤为突出。 图3-10 整体刚度矩阵存储方法 h. 整体刚度矩阵稀疏阵。 故整体刚度矩阵不能求逆，必须作约束处理方能正确地将结点位移求出，进而求出结构的应力场。 (a) 带状分布规律 (b) 带状存储 * 杆件系统的有限元法（十三） 总体刚阵存储与节点排列 1个节点有2个自由度，n个节点的系统有2n个自由度。 总体刚阵是2n×2n阶。n＝100,000时，总刚有4×1010个元素。全部存储，数据量非常大。 总刚的特点 对称——存储上（下）三角，有n×(2n+1)个 元素，还是太多 带状、稀疏——数据中大量是零，在对角线两 侧带状分布 带状存储——仅存储上（下）三角的带内元素， 数据量减少很多 最大带宽——由A－B的元素数 最大半带宽——由O－B的元素数 最大半带宽愈小，存储量愈小 * 杆件系统的有限元法（十四） 总体刚阵存储与节点排列（续） 单元自由度的总体拓扑关系 非零元素只出现在由有限单元 连接的节点的自由度的位置 最大半带宽 节点号 自由度号 单元 总体 单元 总体 1 i 1 2 i－1 2 2 i 2 j 3 2 j －1 4 2 j 总刚 * 杆件系统的有限元法（十五） 总体刚阵存储与节点排列（续） 每行存储8个元素，少于8×12＝96个元素 每行存储6个元素，少于6×12＝72个元素 大型结构，采用等带宽存储，较上（下）三角存储，少80％～90％，更节省存储空间的方法是变带宽存储，每行的半带宽不同，存储的元素数也不同 * 杆件系统的有限元法（十六） 总体刚阵存储与节点排列（续） 每行存储8个元素，少于8×12＝96个元素 每行存储6个元素，少于6×12＝72个元素 一般来讲，当一个结构较长时，应先顺其较窄的方向编号然后向较长的方向移动编号。 * 小结 弹簧单元和杆单元 基本参数： 单元刚度矩阵 总体刚度矩阵 总刚的组装：单元节点和总体节点的拓扑关系 边界条件的处理：消元法和大数法 反力的计算：回代 总刚的存储，对称、带状、稀疏，带宽、半带宽 * 矩阵代数复习 1、矩阵定义 一组元素按行、列次序排列成的矩形阵列称为矩阵。若矩阵的元素排列为m 行和n列，称为m?n 阶矩阵。 A = a a a a a a a a a n n m m mn 11 12 1 21 22 2 1 2 L L M O M L é ? ê