有限元 整体分析.ppt

3-3 单元刚度矩阵 3-4 单元刚度矩阵的物理意义及其性质 3-5 整体分析 3-6 整体刚度矩阵的形式 * 根据虚功原理，有 则 建立了单元的结点力与结点位移之间的关系， 称为单元刚度矩阵。它是6*6矩阵，其元素表示该单元的各结点沿坐标方向发生单位位移时引起的结点力，它决定于该单元的形状、大小、方位和弹性常数，而与单元的位置无关，即不随单元或坐标轴的平行移动而改变。 由于[D]中元素是常量，而在线性位移模式下，[B]中的元素也是常量，且 因此 可以进一步得出平面应力问题和平面应变问题中的单元刚度矩阵。 已经求出了下列关系 结点力和结点位移的关系：(以简单平面桁架为例) 平面问题中，离散化的单元组合体极为相似，单元组合体在结点载荷的作用下，结点对单元、单元对结点都有作用力与反作用力存在，大小相等方向相反，统称为结点力。 结点力和结点位移的关系前面已经求出： 单元刚度矩阵的物理意义： 将 写成分块矩阵 写成普通方程 其中 表示结点S(S=i,j,m)产生单位位移时，在结点r(r=i,j,m)上所需要施加的结点力的大小。 将结点力列矩阵 与结点位移列矩阵 均展开成(6*1)阶列矩阵，单元刚度矩阵相应地展开成(6*6)阶方阵： 元素K的脚码，标有“-”的表示水平方向，没有标“-”的表示垂直方向。 单元刚度矩阵的物理意义： 单元刚度矩阵的物理意义： 单元刚度矩阵的每一个元素都有明显的物理意义。 表示结点S(S=i,j,m)在水平方向、垂直方向产生单位位移时，在结点r(r=i,j,m)上分别所要施加的水平结点力和垂直结点力的大小。例如 表示结点j在垂直方向产生单位位移时，在结点i所需要施加的水平结点力的大小。 单元刚度矩阵的性质： 1)对称性： 是对称矩阵 2)奇异性： 是奇异矩阵， 单元刚度矩阵所有奇数行的对应元素之和为零，所有偶数行的对应元素之和也为零。由此可见，单元刚度矩阵各列元素的总和为零。由对称性可知，各行元素的总和也为零。 单元刚度矩阵的性质： 例题：求下图所示单元的刚度矩阵，设 1、求[B] 2、求[D] 3、求[S] 4、求 图示结构的网格共有四个单元和六个结点。在结点1、4、6共有四个支杆支承。结构的载荷已经转移为结点载荷。 对单元的分析得出单元刚度矩阵，下面，将各单元组合成结构，进行整体分析。 将各单元组合成结构，进行整体分析。 整体分析分4个步骤 1、建立整体刚度矩阵； 2、根据支承条件修改整体刚度矩阵； 3、解方程组，求出结点位移；(消去法与叠加法) 4、根据结点位移求出应力。 不需作坐标转换 ； 结构刚度矩阵的形成与杆系相同； 1、建立整体刚度矩阵(也叫作结构刚度矩阵) 上图中的结构有六个结点，共有12个结点位移分量和12个结点力分量。由结构的结点位移向量求结构的结点力向量时，转换关系为： 分块形式为： 其中子向量 和 都是二阶向量，子矩阵 是二行二列矩阵。整体刚度矩阵[K]是12*12阶矩阵。 2、根据支承条件修改整体刚度矩阵。 建立整体刚度矩阵时，每个结点的位移当作未知量看待，没有考虑具体的支承情况，因此进行整体分析时还要针对支承条件加以处理。 在上图的结构中，支承条件共有四个，即在结点1、4、6的四个支杆处相应位移已知为零： 建立结点平衡方程时，应根据上述边界条件进行处理。 3、解方程组，求出结点位移。 通常采用消元法和迭代法两种方法。 4、根据结点位移求出应力。 整体刚度矩阵 是单元刚度矩阵 的集成。 1、刚度集成法的物理概念： 刚度矩阵中的元素是刚度系数，即由单位结点位移引起的结点力。 由3-5节的例题可见，与结点2和3相关的单元有单元①和③，当结点3发生单位位移时，相关单元①和③同时在结点2引起结点力，将相关单元在结点2的结点力相加，就得出结构在结点2的结点力。由此看出，结构的刚度系数是相关单元的刚度系数的集成，结构刚度矩阵中的子块是相关单元的对应子块的集成。 2、刚度矩阵的集成规则： 先对每个单元求出单元刚