5-有限元法-结构静力学问题幻灯片.ppt

* * * * * * * * * * * * * * * * * 5-7 整体刚度矩阵的特点 3、带形分布规律。 上图中，矩阵[K]的非零元素分布在以对角线为中心的带形区域内，称为带形矩阵。在半个带形区域中(包括对角线元素在内)，每行具有的元素个数叫做半带宽，用d表示。半带宽的一般计算公式是： 半带宽 d = ( 相邻结点码的最大差值 + 1 ) * 2 上图中相邻节点码的最大差值为4，故d=(4+1)*2=10 利用带形矩阵的特点并利用对称性，可只存贮上半带的元素，叫半带存贮。 5-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列元素 元素 5-7 整体刚度矩阵的特点 同一网格中，如果采用不同的节点编码，则相应的半带宽d也可能不同。如图，是同一网格的三种节点编码，相邻节点码的最大差值分别为4、6、8，半带宽分别为10、14、18。因此，应当采用合理的节点编码方式，以便得到最小的半带宽，从而节省存贮容量。 5-8 支承条件的处理 无约束结构的整体刚阵是奇异的，即整体平衡方程 的解不唯一，所以，必须引入几何约束，才能求得唯一解。位移约束常分为：节点固定和给定节点位移两种约束。由于引入位移约束条件通常在整体刚阵及节点载荷形成后进行，即此时[K]、{R}中的元素均已按一定顺序分别储存于相应的数组，故引入位移约束时，要求尽量不要打乱[K]、{R}的储存顺序。 5-9 单元载荷移置 连续弹性体离散为单元组合体时，为简化受力情况，需把弹性体承受的任意分布的载荷都向节点移置(分解)，而成为节点载荷。如果弹性体承受的载荷全都是集中力，则将所有集中力的作用点取为节点，就不存在移置的问题，集中力就是节点载荷。但实际问题往往受有分布的面力和体力，都不可能只作用在节点上。因此，必须进行载荷移置。如果集中力的作用点未被取为节点，该集中力也要向节点移置。 将载荷移置到节点上，必须遵循静力等效的原则。静力等效是指原载荷与节点载荷在任意虚位移上做的虚功相等。在一定的位移模式下，移置结果是唯一的，且总能符合静力等效原则。 5-9 单元载荷移置 载荷移置的原则：能量等效，即单元的实际载荷与移置后的节点载荷在相应的虚位移上所做的虚功相等 载荷移置的条件：圣维南原理 载荷移置的方法： 1）直接法（静力等效法，虚功移置法） 2）普遍公式法 0.5ql 0.5ql 0.5ql 0.5ql M M 静力等效 5-9 单元载荷移置 虚功移置：在线性位移模式下，对于常见的一些载荷，可以通过简单的虚功计算得节点载荷。即移置前后虚功相等。 如均质等厚度的三角形单元所受的重力，把1/3的重力移到每个节点，即 5-9 单元载荷移置 例： 总载荷的2/3移置到节点i，1/3移置到节点j，与原载荷同向 5-9 单元载荷移置 普遍公式法 集中力的移置 体力的移置 分布面力的移置 在线性位移模式下，用直接计算法简单；非线性模式下，要用普遍公式计算。 5-10 应力计算 整体方程在经过边界条件处理后，可选用适当的线性方程组求解方法求得各个节点的位移值 ，而单元内任意一点应力的计算可利用单元的位移函数和弹性矩阵求得： 但是，对于单元公共节点上的应力，由于据不同单元求得的结果是不同的，必须进行处理，常用方法有： 1）绕点平均法 2）单位面积加权平均法 3）精确计算法 本章结束 * * * * * * * * * * * * * * * * * * * * * 5-2 平面问题的常应变(三角形)单元 应力矩阵 由平面问题物理方程得： 应变矩阵[B]反映了单元内任一点的应变与节点位移间的关系 应力矩阵[S]反映了单元内任一点的应力与节点位移间的关系 显然，常应变三角形单元的应变矩阵[B]为常量矩阵，说明在该单元上的应力和应变为常值。由此可见，在相邻单元的边界处，应变及应力不连续，有突变。 5-3 单元刚度矩阵 讨论单元内部的应力与单元的节点力的关系，导出用节点位移表示节点力的表达式。 由应力推算节点力，需要利用平衡方程。用虚功方程表示出平衡方程，即外力在虚位移上所作的