5有限元解平面问题.ppt

在整理应力成果时，读者应注意，应用三角形单元时， （1）采用两单元平均法和绕结点平均法的 应力成果比较接近，但前者的精度略 好于后者。 （2）边界面的应力，宜采用向外插值的方 法求出。 例题4 图（a）所示的深梁，在跨中受集中力F的作用，若取 试用有限单元法求解跨中的位移。 解 1.将图 划分网格，化为离散化结构，如图（b）所示。由于结构具有对称性，可取 部分进行分析，如 所示。 （a） 图（c） 2. 中，只有两个未知结点位移 其余的结点位移均为零。 未知的结点位移列阵是 对应的结点荷载列阵是 3.下面我们直接来建立对应于未知结点 位移的平衡方程式， 图（c） 应用几何方程，求出单元的应变列阵： 3.单元应变和应力 S——称为应力转换矩阵，写成分块形式为 再应用物理方程，求出单元的应力列阵： B ——称为应变矩阵，用分块矩阵表示， 对于线性位移模式，求导后得到的应变和应力，均成为常量，因此，称为常应变（应力）单元。相邻单元的应力是跳跃式的。 4.单元的结点力列阵与刚度矩阵 考虑其中一个单元： （2）单元与周围的单元在边界上已没有联系，只在结点处互相联系。 （1）将作用于单元上的各种外荷载，按静力等效原则移置到结点上去，化为等效结点荷载。故单元内已没有外荷载。 结点作用于单元上的力，称为结点力， 假设发生一组结点虚位移 　　 则单元内 任一点（x,y）的虚位移为 单元内任一点（x,y）的虚应变为 代入虚功方程：在单元中，外力（结点力 ）在虚位移（结点虚位移　 　）上的虚功，等于应力 在虚应变 上的虚功，即： 其中 与 无关，故式(a) 成为 式(b)是由应力求结点力的一般公式。 因为 是独立的任意的虚位移，虚功方程对任意的 均应满足，可得出 代入 (b) 式（c）是由结点位移求结点力的一般公式， －－称为单元的劲度矩阵 K 其中： 再将应力公式代入上式，得 （c） （d） 对于三角形单元，B矩阵内均为常数， 有 代入B，D，即得平面应力问题中三结点三角形单元的刚(劲)度矩阵，可写成如下分块矩阵的形式： 其中， 结点力与结点位移的关系（c)式可写为（子块形式） (1)k是6×6的方阵, k中元素 表示仅在单元结点s沿n方向产生单位位移时引起结点r沿l方向的结点力。 (2)由反力互等定理， 所以k是对称矩阵，以对角线为对称轴。 单元劲度矩阵k的性质： (3)当单元作刚体平移时，如 三角形内不产生应力和应变，结点力也为0。 (4)由(3)可导出行列式 (即k为奇异矩阵 )。 (5)k的元素与　　 单元的形状和方位等有关，但与单元的大小和刚体的平动以及作 度转动无关。 即,k中每一行（或列）元素之和为0（其中第1、3、5元素之和(对应x向)或2、4、6元素之和(对应y向)也为0）。 某等腰直角三角形单元ijm如图所示,已知在所选取的坐标系中,单元结点坐标分别为: 例题 应用 可得 应用公式 可得该单元的应力转换矩阵为 可得该单元的应力转换矩阵为 应用教材式(6-37)及式(6-38)可得该单元的单元刚度矩阵为 现考察结点力与单元中的应力之间的关系。为了简单起见，假定只有结点i发生位移ui，如右图(a)所示。由上面的单元刚度矩阵得相应的结点力为： 其中， 。相应的结点位移及结点力如图所示 另一方面，由于发生了位移ui,则根据上面得到的应力转换矩阵，可得应力分量为 如上图(b)中单元的两直角面所示。根据单元的平衡条件,还可得出斜面上的应力,如图所示。若将这三个面上的应力分别按静力等效原则移置到结点上去，可以得到图(a)中相同的结点力。 ■ §5-4 荷载向结点移置?单元的结点荷载列阵 在FEM中，须将作用于单元中的外荷载向结点移置，化为等效结点荷载， （2）变形体静力等效原则－－在任意的虚位移上，使原荷载与移置荷载的虚功相等。 1、等效原则 （1）刚体静力等效原则－－使原荷载与移置荷载的