第9章弹性力学平谁面问题的有限元法.pptVIP

第9章 弹性力学平面问题的有限元法 2D连续体问题总体刚度矩阵的半带宽: 9.7 约束条件的处理 位移边界条件在大多数情形下有两类: 第一类:零位移边界 第二类: 给定具体数值的 位移边界 第一类:零位移边界 方法:对角元素置“1”法 不变 不变 不变 不变 第二类: 给定具体数值的位移边界 方法:对角元素乘大数法 9.8 单元等效节点载荷列阵 1、体积分布力 2、分布面力 x q m j i y α s 将水平位移分量和节点坐标公式中的第一式， （下标i，j，m 轮换） 其中: 同理: 令 （下标i，j，m轮换） Ni 称为形态函数 [N]称为形态矩阵 单元内的位移函数用节点位移表示，可以简写成， 单元的结点位移记为 单元内的位移记为 选择单元位移函数应满足以下条件： 1）位移模式必须反映单元的刚体位移：常数项 2）位移模式必须反映单元的常量应变：线性 3）相邻单元在公共边界上的位移连续，单元之间不能重叠，也不能脱离。即位移函数在单元之间连续，称为协调性条件。 满足1）、2）条原则的称为完备性单元。 同时满足三条原则的称为完备协调单元。 形态函数N(x,y)具有以下性质： 1）形态函数在单元节点上形态函数的值，具有“本点为1，它点为0”。 2）形态函数在单元中任一点，三个形态函数之和为1。 3）三角形单元任意一边上的形函数，仅与该边的两端节点坐标有关，而与其他节点坐标无关。 例如：在 ij 边上（即i，j节点），有 例题1：如图所示等腰三角形单元，求其形态矩阵[N]。 三角形积为 形态函数为 形态矩阵为 例2 、证明：对平面三角形单元形函数存在下列关系 节点位移 单元内位移分量 应变 应力 节点力 几何方程 物理方程 平衡方程 [K] 位移模式 [N] 9.3 单元刚度矩阵 用节点位移表示单元内部各点位移： 1 由节点位移求应变 记为 ，[B]矩阵称为几何矩阵。 [B]矩阵可以表示为分块矩阵的形式 由几何方程可以得到单元的应变表达式 关于几何矩阵 1) 几何矩阵中每一元素 bi 的物理意义：单元某一节点有单位位移且其它节点位移皆为0时，所引起单元内部的应变分布。 2) 几何矩阵中各元素为常数（单元位移确定后，节点坐标值为定值），因而应变分量为常量，平面三角形单元为常应变单元。（这是由于采用线性位移函数的结果）。 1、推导杆单元的形态函数、形态矩阵、几何矩阵。 练习： E A l 1 2 x 2、已知三角形单元的形态矩阵为 根据形态矩阵求三角形单元的几何矩阵 2 由应变求应力 由物理方程，可以得到单元的应力表达式 [D]称为弹性矩阵，对于平面应力问题 对于平面应变问题 3 由节点位移求应力 定义 为应力矩阵。 将应力矩阵分块表示为 关于应力矩阵 1) 应力矩阵中每一元素的物理意义：单元某一节点有单位位移且其它节点位移皆为0时，所引起单元内部的应力分布。 2) 应力矩阵中各元素为常数，因而应力分量为常量，平面三角形单元为常应力单元。 4 由应力求节点力——根据虚位移原理 5 由节点位移求节点力 单元刚度矩阵 例题2：如图所示等腰三角形单元，求其刚度矩阵[K]，设μ=0。 解：1）求[B] 2）求[D] 3）求[S] 4）求[K] 9.4 单元刚度矩阵的性质与物理意义 （一）单元刚度矩阵的物理意义 假设单元的节点位移如下： 由 ，得到节点力如下： Kix,ix表示i节点在水平方向产生单位位移时，在节点i的水平方向上需要施加的节点力。 Kiy,ix表示i节点在水平方向产生单位位移时，在节点i的垂直方向上需要施加的节点力。 因此单元刚度矩阵中每个元素都可以理解为刚度系数，即在结点产生单位位移时需要施加的力。 （二）单元刚度矩阵的性质 1) 对称性 2) 奇异性 3) 分块性 9.5 整体刚度矩阵的形成 基本方法是刚度集成法，由单元刚度矩阵中的元素累加得到整体刚度矩阵中的元素，即整体刚度矩阵是单元刚度矩阵的集成。 如何得到整体刚度矩阵？ 刚度集成法即结构中的节点力是相关单元节点力的叠加，整体刚度矩阵的元素是相关单元的单元刚度矩阵元素的集成。节点3在整体刚度矩阵的对应系数，应该是单元（1）、（3）、（4）中对应系数的集成。 i i j i j j m j m m i m 单元(1): i，j，m 单元的局部码与整体的总码的对应关系： 1，2，3 单元(2): i，j，m 2，4，5 单元(3): i，j，m 2，5，3 单元(4): i，j，m 3，5，6 i j m i j m 1 2 3 1 2 3 4 5 6 4 5 6 1 2 3 1 2 3 4 5 6 4 5 6 1 2 3 1 2 3 2 3 2 4 5 2 4 5 5 3