(精选)模态分析与综合技术第6章 FEM一般原理和表达格式课件.ppt

第6章 FEM一般原理和表达格式;6.2 常应变三角形（平面应力问题） 由于三角形单元对复杂边界有较强的适应能力，可以用其逼近任何形状，因此很容易将一个平面结构离散成有限个三角形单元。故在二维问题的分析中，三角形单元是最为流行的单元。; 6.2.2 单元分析 任取一个一般性的单元，如图所示。三个节点的编号为i, j, k。节点位移为( ui, vi )、( uj, vj ) 、( uk, vk )　，单元节点位移为：; a1～ a6称为广义坐标。 在节点处应有: ; 当 i, j, k的位置为逆时针排列时，2Δ恒正，且等于三角形单元面积的两倍。将这些结果代入; 可以合并成; 形函数具有以下性质： (i)在节点上形函数的值为：;也必然要求： Ni=1， Nj= Nk=0。同样利用别的节点的位移协调可以其他两个形函数也具有同样的性质。;第6章 FEM一般原理和表达格式; 2　 单元的应变、应力 利用：; 其中 ; 6.2.3 利用最小位能原理建立有限元方程 ;第6章 FEM一般原理和表达格式;对单元使用最小位能原理，即：;称为单元刚度矩阵。对于3节点三角形单元，应变矩阵B为常量阵，因此有：;称为体力f 和面力T的单元等效节点载荷列阵。;　　单元的单位体积重量为r，如图所示。则有：;自重的等效节点载荷是：;　　沿x方向作用在ij边，分布集度为q,如图所示。则边上的面力为：; 上面例子中单元①和②的边上作用着均布的外载荷，可以把它们的合力平分到两节点，如图所示。 ; 将单元节点位移u用结构节点位移a表示：; 使用最小位能原理，整个离散系统的总位能对结构节点位移的变分为0，可以得到系统的有限元求解方程：;其中：; 6.2.4 解方程获得位移及应力 解上述方程可得到各非约束自由度的位移，再由; 第2节给出了三角形单元的广义坐标有限元法的具体的实施步骤。 但在用有限元法进行分析问题时，所遇到的单元种类繁多，如二维单元有三角形单元、矩形单元，三维单元有四面体、五面体或平行六面体等，即使同样形状的单元还可有不同的节点数目。 如何选择合适的单元进行计算？主要要考虑以下因素：求解问题的类型、计算精度要求及计算量等。 下面给出广义坐标有限元法的一般步骤。;　　1 以广义坐标b为待定参数，给出单元内位移u:;2 用节点位移ae表示广义坐标b:; 3 用单元节点位移ae表示单元位移函数u，得到插值函数N; 5 用最小位能原理建立离散体系的平衡方程; 一管道系统，具体结构和参数见系统结构图，采用单元：Pipe16。计算其前5阶固有频率和振型。