有限元方法3-平面问题.pptVIP

* * * * 每个子块为 2X2，K为6X6 F KQ ; F T*f, Q T*q Tf KT*q; T*k*q K*T*q * * * * * * * * * * 单元刚度矩阵示例 例题：求下图所示单元的刚度矩阵，设 1、求几何矩阵[B] 2、求弹性矩阵[D] 3、求刚度矩阵 平面应力问题 3-3 整体分析 图示结构的网格共有四个单元和六个节点。在节点1、4、6共有四个支杆支承。结构的载荷已经转移为结点载荷。 整体分析的四个步骤： 1、建立整体刚度矩阵； 2、增加支承条件； 3、解方程组，求节点位移； 4、根据节点位移求出应力。 单元分析得出单元刚度矩阵，下面将各单元组合成结构，进行整体分析。 3-3 整体分析 1、建立整体刚度矩阵 也叫作结构刚度矩阵 上图中的结构有六个节点，共有12个节点位移分量和12个节点力分量。由结构的节点位移向量求结构的节点力向量时，转换关系为： 分块形式为： 其中子向量 和 都是二阶向量，子矩阵 是二行二列矩阵。整体刚度矩阵[K]是12*12阶矩阵。 整体坐标系下单刚矩阵 变换矩阵 矩阵变换 T中的元素 整体坐标系下单元刚度矩阵 3-5 整体刚度矩阵的形式 n为节点数目 例4-7 单元节点编码 单元刚度矩阵 总体刚度矩阵 3-6 整体刚度矩阵的特点 在有限元法中，整体刚度矩阵的阶数通常是很高的，在解算时常遇到矩阵阶数高和存贮容量有限的矛盾。找到整体刚度矩阵的特性达到节省存贮容量的途径。 1、对称性。 只存贮矩阵的上三角部分，节省近一半的存贮容量。 2、稀疏性。 矩阵的绝大多数元素都是零，非零元素只占一小部分。 3-7 整体刚度矩阵的特点 2、稀疏性 续） 矩阵的绝大多数元素都是零，非零元素只占一小部分。 节点5只与周围的六个节点 2、3、4、6、8、9 用三角形单元相连，它们是5的相关节点。只有当这七个相关节点产生位移时，才使该节点产生节点力，其余节点发生位移时并不在该节点处引起节点力。因此，在矩阵[K]中，第5行的非零子块只有七个 即与相关节点对应的七个子块 。 3-7 整体刚度矩阵的特点 2、稀疏性 续） 一般的，一个节点的相关结点不会超过九个，如果网格中有200个节点，则一行中非零子块的个数与该行的子块总数相比不大于9/200，即在5%以下，如果网格的节点个数越多，则刚度矩阵的稀疏性就越突出。 利用矩阵[K]的稀疏性，可设法只存贮非零元素，从而可大量地节省存贮容量。 3-7 整体刚度矩阵的特点 3、带形分布规律 上图中，矩阵[K]的非零元素分布在以对角线为中心的带形区域内，称为带形矩阵。在半个带形区域中 包括对角线元素在内 ，每行具有的元素个数叫做半带宽，用d表示。半带宽的一般计算公式是： 半带宽 d 相邻结点码的最大差值 + 1 * 2 上图中相邻节点码的最大差值为4，故d 4+1 *2 10 利用带形矩阵的特点并利用对称性，可只存贮上半带的元素，叫半带存贮。 3-7 整体刚度矩阵的特点 图 a 中的矩阵[K]为n行n列矩阵，半带宽为d。半带存贮时从[K]中取出上半带元素，按图 b 中的矩阵 的排列方式进行存贮，即将上半部斜带换成竖带。存贮量n*d，存贮量与[K]中元素总数之比为d/n，d值越小，则存贮量约省。 矩阵[K] 矩阵 对角线 第1列 r行 r行 r列 45度斜线 r行s列 r行s-r+1列元素 元素 3-7 整体刚度矩阵的特点 同一网格中，如果采用不同的节点编码，则相应的半带宽d也可能不同。如图，是同一网格的三种节点编码，相邻节点码的最大差值分别为4、6、8，半带宽分别为10、14、18。因此，应当采用合理的节点编码方式，以便得到最小的半带宽，从而节省存贮容量。 3-8 支承条件的处理 无约束结构的整体刚阵是奇异的，即整体平衡方程 的解不唯一，所以，必须引入几何约束，才能求得唯一解。 位移约束常分为：节点固定和给定节点位移两种约束。由于引入位移约束条件通常在整体刚阵及节点载荷形成后进行，即此时[K]、 R 中的元素均已按一定顺序分别储存于相应的数组，故引入位移约束时，要求尽量不要打乱[K]、 R 的储存顺序。 3-9 单元载荷移置 连续弹性体离散为单元组合体时，为简化受力情况，需把弹性体承受的任意分布的载荷都向节点移置 分解 ，而成为节点载荷。如果弹性体承受的载荷全都是集中力，则将所有集中力的作用点取为节点，就不存在移置的问题，集中力就是节点载荷。但实际问题往往受有分布的面力和体力，都不可能只作用在节点上。因此，必须进行载荷移置。如果集中力的作用点未被取为节点，该集中力也要向节点移置。 将载荷移置到节点上，必须遵循静力等效的原则。静力等效是指原载荷与节点载荷在任意虚位移上做的虚功相等。在一定