（金）(新)第三讲平面问题有限元分析——精品.pptVIP

* * * * * * * * * * * * * * * * * * * * * * 2-10 支承条件的处理 同理，如果结点n有竖向支杆，则平衡方程的第2n个方程应改为 ，为此，在矩阵[K]中，第2n行的对角线元素改为1，该行全部非对角线元素改为0，同时，第2n列全部非对角线元素也改为0。在{P}中，第2n个元素改为0。 2-10 支承条件的处理 2-8节中的结构，结点1有水平支杆，结点2有两个支杆，结点3有竖向支杆。对支承条件处理后，矩阵修改为： 2-10 支承条件的处理 最后考虑支点n的水平位移 为已知非零值 的情况，这时的支承条件为 对平衡方程的第2n-1个方程作如下修改： 对角线系数 乘以一个大数A ( 例如A=1010 )； 右边自由项换成 ； 其余各项保持不变，即有 实际上，上式中除包含大数A的两项外，其他各项相对地都比较小，可以忽略不计，因此与支承条件等价。 对于支点的竖向位移 的支承条件，也可用同样的方法进行修改。 2-11 整体刚度矩阵的特点 在有限元法中，整体刚度矩阵的阶数通常是很高的，在解算时常遇到矩阵阶数高和存贮容量有限的矛盾。找到整体刚度矩阵的特性达到节省存贮容量的途径。 1、对称性。 只存贮矩阵的上三角部分，节省近一半的存贮容量。 2、稀疏性。 矩阵的绝大多数元素都是零，非零元素只占一小部分。 2-11 整体刚度矩阵的特点 2、稀疏性。 矩阵的绝大多数元素都是零，非零元素只占一小部分。 结点5只与周围的六个结点(2、3、4、6、8、9)用三角形单元相连，它们是5的相关结点。只有当这七个相关结点产生位移时，才使该结点产生结点力，其余结点发生位移时并不在该结点处引起结点力。因此，在矩阵[K]中，第5行的非零子块只有七个(即与相关结点对应的七个子块)。 2-11 整体刚度矩阵的特点 2、稀疏性。 一般，一个结点的相关结点不会超过九个，如果网格中有200个结点，则一行中非零子块的个数与该行的子块总数相比不大于9/200，即在5%以下，如果网格的结点个数越多，则刚度矩阵的稀疏性就越突出。 利用矩阵[K]的稀疏性，可设法只存贮非零元素，从而可大量地节省存贮容量。 2-11 整体刚度矩阵的特点 3、带形分布规律。 上图中，矩阵[K]的非零元素分布在以对角线为中心的带形区域内，称为带形矩阵。在半个带形区域中(包括对角线元素在内)，每行具有的元素个数叫做半带宽，用d表示。半带宽的一般计算公式是： 半带宽 d = ( 相邻结点码的最大差值 + 1 ) * 2 上图中相邻结点码的最大差值为4，故d=(4+1)*2=10 利用带形矩阵的特点并利用对称性，可只存贮上半带的元素，叫半带存贮。 2-11 整体刚度矩阵的特点 图(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列元素 元素 2-11 整体刚度矩阵的特点 同一网格中，如果采用不同的结点编码，则相应的半带宽d也可能不同。如图，是同一网格的三种结点编码，相邻结点码的最大差值分别为4、6、8，半带宽分别为10、14、18。因此，应当采用合理的结点编码方式，以便得到最小的半带宽，从而节省存贮容量。 * * * * * * * * * * * * * * * * * * * * * * * 例题：在均质、等厚的三角形单元ijm的ij边上作用有 沿x方向按三角形分布的载荷，求移置后的结点载荷。 取局部坐标s，在i点s=0，在j点s=l，L为ij边的长度。 在ij边上，以局部坐标表示的插值函数为， ， ， 载荷为 例： 总载荷的2/3移置到结点i，1/3移置到结点j，与原载荷同向 2-5 单元应力矩阵 本节利用几何方程、物理方程，实现用结点位移表示单元的应变和单元的应力。 用结点位移表示单元的应变的表达式为 ，[B]矩阵称为几何矩阵。 2-