三角形平面问题摘要.pptVIP

* * * * * * * * * * * * * 10 支承条件的处理 2）对角元素置1法 设边界条件为： 整体方程记为： 为了不改变[K]列数，将第i列的元素乘以已知位移移到等号右，而第i列位置补0；同时，为了不改变[K]的行数，使其体现约束条件，将第i行对角元素置1，其他元素置0。如右： 10 支承条件的处理 图示结构中，节点1有水平约束，节点4有两个方向约束，节点6有竖向约束。对支承条件处理后，整体刚度矩阵修改为： 10 支承条件的处理 3）对角元素乘大数法 将第i行主对角元素乘一个大数，相应载荷分量改为 ，其余元素不变。如右： 原因为：展开第i行有： 两边同除大数得： 近似有： 10 应力计算 整体方程在经过边界条件处理后，可选用适当的线性方程组求解方法求得各个节点的位移值 ，而单元内任意一点应力的计算可利用单元的位移函数和弹性矩阵求得： 但是，对于单元公共节点上的应力，由于据不同单元求得的结果是不同的，必须进行处理，常用方法有： 1）绕点平均法 2）单位面积加权平均法 3）精确计算法 * * * * * * * * * * * * * * * * * * * * * * * * * * * * 4 单元刚度矩阵的物理意义及其性质 上述推导单元刚度矩阵的过程可归纳为 单元刚阵[K]的物理意义是单元受节点力作用后抗变形的能力。其元素 的意义为：当第j个自由度发生单位位移，而其他自由度的位移为0时，在第i个自由度上所施加的力。若按节点来说明，则刚阵中每个子块 表示：当节点j处发生单位位移，而其他节点固定时，在节点i上所施加的力。 { } s [ ] tA B T [ ] tA ] B ][ D [ ] B [ K T e = { } e d { } e { } e F [ ] D [ ] B [ ] ] B ][ D [ S = (6) (3) (3) (6 ╳ 3) (3 ╳ 3) (3 ╳ 6) (3 ╳ 6) (6 ╳ 6) 4 单元刚度矩阵的物理意义及其性质 节点力和节点位移的关系：(以简单平面桁架为例) 平面问题中，离散化的单元组合体极为相似，单元组合体在节点载荷的作用下，节点对单元、单元对节点都有作用力与反作用力存在，大小相等方向相反，统称为节点力。 节点力和节点位移的关系前面已经求出： 4 单元刚度矩阵的物理意义及其性质 单元刚度矩阵的物理意义： 将 写成分块矩阵 写成普通方程 其中 表示节点S(S=i,j,m)产生单位位移时，在节点r(r=i,j,m)上所需要施加的节点力的大小。 4 单元刚度矩阵的物理意义及其性质 单元刚度矩阵的物理意义： 将节点力列矩阵 与节点位移列矩阵 均展开成(6*1)阶列矩阵，单元刚度矩阵相应地展开成(6*6)阶方阵： 元素K的脚码，标有“-”的表示水平方向，没有标“-”的表示垂直方向。 4 单元刚度矩阵的物理意义及其性质 单元刚度矩阵的物理意义： 单元刚度矩阵的每一个元素都有明显的物理意义。 表示节点S(S=i,j,m)在水平方向、垂直方向产生单位位移时，在节点r(r=i,j,m)上分别所要施加的水平节点力和垂直节点力的大小。例如 表示节点j在垂直方向产生单位位移时，在节点i所需要施加的水平节点力的大小。 4 单元刚度矩阵的物理意义及其性质 1）单元刚度矩阵是对称阵，(只要证明 ) 2）单元刚阵主对角线元素恒为正值；因为主对角元素 表示力的方向和位移方向一致，故功总为正值。 3）单元刚阵是奇异阵，即|K|=0，这是因为计算单元刚阵时没有对单元的节点加以约束，虽然，单元处于平衡状态，但容许单元产生刚体位移，故从单元刚度平衡方程不可能得到唯一位移解 ，只能得到唯一的节点力解。 4）单元刚阵所有奇数行的对应元素之和为零，所有偶数行的对应元素之和也为零。由此可见，单元刚阵各列元素的总和为零。由对称性可知，各行元素的总和也为零。 4 单元刚度矩阵的物理意义及其性质 例题：求下图所示单元的刚度矩阵，设 1、求[B] 2、求[D] 3、求[S] 4、求 4 单元刚度矩阵的物理意义及其性质 几点说明： 1）单元刚度方程是满足节点力平衡条件而建立的，即有限元方程是一组节