有限元3-单元劲度矩阵n.pptVIP

形函数及其性质; 形函数是有限单元法中的一个重要函数，它具有以下性质：;（i、j、m）;x;性质3 在三角形单元的边界ij上任一点（x，y），有 ;性质4 形函数在单元上的面积分和在边界上的线积分公式为 ;计算单元位移函数举例;计算单元位移函数举例;计算单元位移函数举例;3、位移模式与解答的收敛性; （1）位移函数的个数 等于单元中任意一点的位移分量个数。本单元中有u和v，与此相应，有2个位移函数；; （4）位移函数中必须包含单元的刚体位移。 ; ◆ 位移函数的形式 一般选为完全多项式。为实现（4）—（6）的要求，根据Pascal三角形由低阶到高阶按顺序、对称地选取；多项式的项数一般应等于单元节点自由度数。;例：平面应力矩形板被划分为若干三角形单元。;①、②、③、④单元的位移函数都是;第三次课 第6章 用有限单元法解平面问题 6－4、5 单元劲度矩阵与相关问题 （单元分析） ; 回顾：单元分析 取结点位移作基本未知量。由结点位移求结点力： 其中，转换矩阵称为单元刚度矩阵。单元分析的主要目的就是要求出单元刚度矩阵。 单元分析的步骤可表示如下： ;6－4、5 单元劲度矩阵与相关问题; 根据几何方程和位移函数可以求得单元应变。;;(6-16a）;上式简写一般式： ; 由于 与x、y无关，都是常量，因此[B]矩阵也是常量。单元中任一点的应变分量是[B]矩阵与单元结点位移的乘积，因而也都是常量。因此，这种单元被称为常应变单元（精度较低!）。;平面应力问题的弹性矩阵; 将应变表达式：;平面应力问题的物理方程;其中：[S]称为单元应力矩阵，并有;(6-20）;(6-21）;3、由单元结点位移表出单元结点力 — 虚功方程;1) 虚功方程（等价于平衡方程和应力边界条件）;2) 结点力与结点位移（实与虚）; ;结点力虚功;单元应力虚功;虚功方程的应用; ; ; 注意到;式(6-24)中子矩阵[krs]为2×2矩阵，且有 ;(6-26）; 平面问题的单元刚度矩阵[k]不随单元（或坐标轴）的平行移动而改变。 由公式(6-25）、(6-26）知，[krs]矩阵和其中的br、cr 、 bs、cs （r、s=i、j、m）有关。 ;i; 平面问题的单元刚度矩阵[k]不随单元的放大或缩小而改变。(板书补充解释) ; （1）单元刚度矩阵中每个元素有明确的物理意义例如，kij表示单元第j个自由度产生单位位移(?j=1)，其他自由度固定（=0）时，在第i个自由度产生的节点力Fi。;（2）[k]的每一行或每一列元素之和为零; ;（3）[k]是对称矩阵 由[k]各元素的表达式，可知[k]具有对称性。 ;单元刚度矩阵性质（ 对称性证明 ）;单元刚度矩阵性质（ 对称性证明 ）;单元刚度矩阵性质（ 对称性证明 ）; 注意到;（4）单元刚度矩阵是奇异矩阵 即[k]的行列式为零（由行列式性质） 。 单元刚度矩阵是在单元处于平衡状态的前提下得出的。单元作为分离体看待，作用在它上面的外力（单元力）必定是平衡力系。然而，研究单元平衡时没有引入约束。承受平衡力系作用的无约束单元，其变形是确定的，但位移不是确定的。所以出现性质（3）中的“平动问题”，即单元可以发生任意的刚体运动。从数学上讲，方程(1-28）的解不是唯一的或不能确定的。由此，单元刚度矩阵一定是奇异的。; ; ; ; 与坐标轴同向为正，反向为负;例：计算平面应力直角三角形单元刚度矩阵; 第一步：计算bi、ci和单元 面积A。 ; 第二步：求子矩阵 ;胆梯档渡号沧酬牌擅党钡针打云睁根织治脐陷汹蚤慷毗彪复缮镑钓藕魂寐有限元3-单元劲度矩阵n有限元3-单元劲度矩阵n