第四章 平面问题的有限单元法解析.ppt

第四章 平面问题的有限单元法 基本条件 （1）等厚度的薄板； （2）体力作用于体内，平行于板的中面，沿板厚不变； （3）面力作用于板边，平行于板的中面，沿板厚不变； （4）约束作用于板边，平行于板的中面，沿板厚不变。 坐标系： 基本条件 （1）很长的常截面柱； （2）体力作用于体内，平行于横截面，沿柱体长度方向不变； （3）面力作用于柱面，平行于横截面，沿柱体长度方向不变； （4）约束作用于柱面，平行于横截面，沿柱体长度方向不变。 坐标系： 以平面应力问题为例介绍平面问题的有限单元法 以最经典的三角形单元为例 应变的离散过程 根据弹性力学中的几何关系，单元内任一点(x,y)的应变表达式为 应变的离散过程 单元内任一点(x,y)的位移(u,v)可以采用节点位移近似表示： 将其代入应变表达式，则 应变的离散过程 为书写方便，应变分量矩阵可用分块矩阵表示 应变的离散过程 由于形函数 所以刚才的矩阵事实上可表示为 应力的离散过程 根据广义虎克定律，对于平面应力问题： 应力的离散过程 如果令： 物理方程简写为 D又称为“弹性矩阵” 将前面应变的表达式代入，则 我们已经知道了应变与位移的关系 那么很自然的 如果发生了虚位移 则会发生虚应变 外力虚功等于内力虚功。 外力虚功 内力虚功 其具体可计算为： 外力虚功等于内力虚功。 结果： 考虑到节点虚位移的任意性： 上式即为有限元控制方程。 此处K称为“刚度矩阵” 如果将求解域划分为多个单元，则 即 三节点等厚三角形单元中B和D的分量均为常量， 则单元刚度矩阵可以表示为 其具体形式为 对于平面应力问题，其具体可计算如下： 把前面获得的有限元控制方程展开： 那么 事实上就是当节点j产生单位位移时，在节点i上需要施加的 节点力。 更具体一点： 当节点i在垂直方向产生单位位移时，在节点i上需要施加的垂直节点力 当节点j在水平方向产生单位位移时，在节点j上需要施加的水平节点力 当节点j在垂直方向产生单位位移时，在节点i上需要施加的水平节点力 当节点i在水平方向产生单位位移时，在节点j上需要施加的垂直节点力 性质1：对称性 性质2：对角线上元素恒为正 性质3：此矩阵为奇异矩阵 意义：没有对节点施加位移约束，所以单元产生任何的刚性位移 都是可以的，由力得不到位移的唯一解。 性质4：此矩阵的各行元素之和为零，由于对称性，各列元素之和 也为零。 单元刚度矩阵形成后，要将单元组成一个整体结构，即整体分析， 基本方法是刚度集成法，即整体刚度矩阵是单元刚度矩阵的集成。 整体刚度矩阵的集成是按对号入座的方式叠加的。 用下面的三角形薄板作为示例： 各个单元的刚度矩阵为： 设单元节点总数为N，每个节点的自由度数为NDOF。 (对于一维情况，NDOF=1；对于二维情况，NDOF=2；三维，NDOF=3)。 Kij 表示j自由度发生单位位移，其他位移为零时， 第i个自由度上必须施加的节点力。 由于整体刚度矩阵具有对称性、稀疏性和非零元素带状分布的 特点，所以没有必要将全部的整体刚度矩阵进行存储。 各行的半带宽D怎么计算： 可用一维数组A来存储半带宽内的元素，而不必储存所有元素。 本例中：总带宽 则可以采用如下方式存储： 最大半带宽是多少？ 相邻节点的编码最大差值+1）×NDOF 教学程序学习阶段 虚功原理建立控制方程 虚功原理建立控制方程 刚度矩阵 总体刚度矩阵（总刚） 单元刚度矩阵（单刚） 单元刚度矩阵 单元刚度矩阵 单元刚度矩阵的物理意义 单元刚度矩阵的物理意义 单元刚度矩阵的性质 单元刚度矩阵的性质 单元刚度矩阵的性质 整体刚度矩阵的形成 共计4个单元，单元节点编号为： 整体刚度矩阵的形成 整体刚度矩阵的形成 如果是二维问题，则总自由度数为2N个， 相应的整体刚度矩阵大小为2N×2N阶方阵。 整体刚度矩阵的意义与性质 总体刚度矩阵中的元素具有如下性质： （1）主对角元素Kii 0 （2）总体刚度矩阵K是对称的奇异矩阵 （3）总体刚度矩阵K是带状稀疏矩阵。 整体刚度矩阵的存储 （1）利用对称性： 只保存整体刚度矩阵上三角的零元素即可； （2）利用稀疏性： 在用分块表示的整体刚度矩阵中，与相关节点对应的分块才能具有非零的元素，其他位置上的分块矩阵的元素为零 （3）利用带状分布： 整体刚度矩阵的非零元素分布在以对角线为中心的带状区域中，每行具有元素的元素的个数叫做“半带宽”，用D表示。