【物理】单元分析.pptVIP

4 弹性力学平面问题的有限元分析 三角形单元中的节点位移如下： 根据克莱姆法则，可求出?1， ?2 ， ?3 同理可得 vi = ?4+?5xi+?6yi vj = ?4+?5xj+?6yj vm = ?4+?5xm+?6ym 解出?4 , ?5 , ?6 ? 写成矩阵形式 单元的应变 由于 单元的应力 根据弹性方程 对于平面应力情况 形函数的性质 Ni 在其余二节点上的值等于零 把面积的行列式以第一行展开乘第二行的代数余子式 把面积的行列式以第一行展开乘第三行的代数余子式 在单元上任一点的三个形函数之和等于 1 面积坐标 面积坐标与形函数的关系 * 杆系问题以结点作为分割单元的“结点”是很自然的，但对于平面问题，待分析物体是连续的，并不存在实际结点。要将物体“拆”成单元，必须用一些假想的线或面作人为地分割。将物体进行分割时，必须保证相邻单元具有公共边界。假定相邻单元仅在一些点（顶点或顶点加边中点）相连接。这些点即为“结点”。实际计算时，可将连续体分成多种形状单元，为讨论简单，现暂时规定只用一种单元来分割。 以位移为未知量的有限元法，最关键的工作是建立单元位移场，因此本节主要介绍各种单元位移场的建立。 引 言 平面问题有限元法可用的单元很多，先介绍最简单的单元：三角形。 第四章 平面问题的有限元分析 §4-1 有限元法基本思想和解题步骤 一、有限元法的基本思想 假想的把一连续体分割成数目有限的小体（单元），彼此间只在数目有限的指定点（结点）出相互连结，组成一个单元的集合体以代替原来的连续体，再在结点上引进等效力以代替实际作用于单元上的外力。选择一个简单的函数来近似地表示位移分量的分布规律，建立位移和节点力之间的关系。 有限元法的实质是：把有无限个自由度的连续体，理想化为只有有限个自由度的单元集合体，使问题简化为适合于数值解法的结构型问题。 二、经典解与有限元解的区别： 微分 数目增到∞ 建立一个描述连续体 经 典 解 法 —— （解析法） 大小趋于 0 性质的偏微分方程 有限单元 离散化 集合 总体分析解 有限元法——连续体——单元——代替原连续体 （近似法） （单元分析） 线性方程组 结构物的离散 假定三角形单元的位移模式 （i，j，m 轮换） (a) 4-1 三角形单元 建立单元内任意点的位移与节点位移的关系，单元节点位移坐标为( xi , yi ), ( xj , yj ), ( xm , ym ) 每一点的位移由下列方程给出，在 i点上 的水平位移方程为： ui=?1+?2 xi +?3 yi uj=?1+?2 xj +?3 yj um=?1+?2 xm +?3 ym 其中 其中 展开后 令 (i, j, m) Ni — 单元的形函数 可得 可以写成 所以，单元的位移模式： [N] — 形态矩阵 根据几何方程 得出 写成矩阵形式 矩阵[B]称为几何矩阵 因此单元内任一点的应变是节点位移的函数, [B] 是常数，所以三角形单元是常应变单元。 其中 （i= i, j, m） 令[S]=[D][B] [S] — 应力矩阵 把[S]矩阵分块，得 其中[Si]如下 (i=i, j, m) （i=i, j, m) 对于平面应变情况 （i=i, j, m) 可知，三角形单元中的应力各处相等 §4-3 形函数的性质及面积坐标 形函数在节点i上的值=1 形函数在节点上的值 在三角形面积表示的行列式中以第一行展开 同理可得 当 当 所以 在三角形单元任一边如 i j 边上的形函数 利用形函数的这一性质可以证明，相邻单元的位移分别进行线性插值之后，在其公共边上将是连续的。 面积坐标的定义 在三角形内任意一点p定义 ?i, ?j, ?m分别表示节点i, j, m所对应的三角形面积。 ?i+ ?j+