第4章平面问题的有限元法-3刚度矩阵.pptVIP

作 业 简述形函数的概念和性质 思考题： 平面问题整体刚度矩阵的推导过程 可见，如此求得的结果与按照静力等效原理将表面力{q}向节点i及j分解所得到的分力完全相同。 图4-7?? ⒊体积力的等效载荷列阵{P} 　　与表面力的情况类似，体积力的等效载荷列阵也是由单元体积力的等效节点力在各节点处合成以后，按节点号码顺序排列而成，即 (i) 式中单元e的体积力的等效节点力为 (j) B：用刚体静力等效原理移置单元载荷。 一、体力的移置： y 0 x i j m 　如果任意三角形单元ijk的重心c上受有自重 则按刚体静力等效原理可把W直接移置到i，j，m三个节点上而组成 ： 　　如果单元的重心c受有惯性力Pu作用，且 ，则Pu移置到i，j，m节点上的等效节点力为： 式中： ——旋转角速度 ——是单元重心处x坐标。 二、面力的移置： y 0 x i j m 已知在ij边受有面力q，则移置到i、j 节点上的等效节点力为： y 0 x i j m 当某一边上有三角形分布的面力时， 可由刚体静力等效直接写出 三、集中力的移置： 如，集中力G做用于其一边界上如图： 先将G分解为 ，后，分别按线段的 比例把 和 分别移置到i，j两点上。 即： 即按静力平衡方法分配。 0 x i j m y * 一. 单元刚度矩阵 第四节 刚度矩阵 为了推导单元的节点力和节点位移之间的关系，可应用虚位移原理对图4-2中的单元e进行分析。单元e是在等效节点力的作用下处于平衡的，而这种节点力可采用列阵表示为 (a) 假设在单元e中发生有虚位移，则相应的三个节点i、j、m 的虚位移为 且假设单元内各点的虚位移为{f *}，并具有与真实位移相同的位移模式。 故有 (c) 参照（4-13）式，单元内的虚应变{? *}为 于是，作用在单元体上的外力在虚位移上所做的功可写为 (d) (f) 而单元内的应力在虚应变上所做的功为 (g) 这里我们假定单元的厚度t为常量。把（d）式及（4-16）式代入上式，并将提到积分号的前面，则有 根据虚位移原理，由（f）和（h）式可得到单元的虚功方程即 注意到虚位移是任意的，所以等式两边与相乘的项应该相等，即得 记 (4-24) 则有 (4-25) 上式就是表征单元的节点力和节点位移之间关系的刚度方程，[k]e就是单元刚度矩阵。如果单元的材料是均质的，那么矩阵[D] 中的元素就是常量，并且对于三角形常应变单元，[B]矩阵中的元素也是常量。当单元的厚度也是常量时,因 ,所以式(4-24)可简写为 [k]e =[B]T [D][B]t? (4-26) 与前面讨论过的情况类似，单元刚度矩阵[k]e中任一列的元素分别等于该单元的某个节点沿坐标方向发生单位位移时，在各节点上所引起的节点力。单元的刚度取决于单元的大小、方向和弹性常数，而与单元的位置无关，即不随单元或坐标轴的平行移动而改变。 将（4-26）式写成分块形式，即可得到平面应力问题中三角形单元的刚度矩阵 (4-27) 其中 ( r = i、j、m；s = i、j、m ) (4-28) 对于平面应变问题，只要将上式中的E、?分别换成E / 1-? 2 和? / 1-? 即可。于是 ( r = i、j、m；s = i、j、m ) (4-29) 二、整体刚度矩阵 讨论了单元的力学特性之后，就可转入结构的整体分析。假设弹性体被划分为N个单元和n个节点，对每个单元按前述方法进行分析计算，便可得到N组形如（4-25）式的方程。将这些方程集合起来，就可得到表征整个弹性体的平衡关系式。 (q) (4-30) 若写成分块矩阵的形式，则 1 2 3 4 2 1 q 图 4-5 组装总刚[k]的一般规则： 1.当[krs]中r=s时，该点被哪几个单元所共有，则总刚子矩阵[krs]就是这几个单元的刚度矩阵子矩阵[krs]e的相加。 3.当[krs]中r和s不同属于任何单元时，则总体刚度矩阵 [krs]=[0]。 下面，我们考查一个组装总刚的实例： 1. 整体刚度矩阵及载荷列阵的组集 根据叠加原理，整体结构的各个刚度矩阵的元素显然是由有关单元的单元刚度矩阵的元素组集而成的，为了便于理解，现结合图4-5说明组集过程。 2.当[krs]中 时，若rs边是组合体的内边，则总体刚度矩