6第六章 用有限单元法解平面问题2011s.pptVIP

第六章 用有限单元法解平面问题;第六章 用有限单元法解平面问题;前言;;;;;;;;二、有限元法处理弹性问题的基本思路;;;;三、有限元法的主要优点;;四、有限单元法发展趋势和研究方向 ;;五、有限元法在土木工程中的应用;;;第一节 基本量及基本方程的矩阵表示;第一节 基本量及基本方程的矩阵表示;第一节 基本量及基本方程的矩阵表示;第一节 基本量及基本方程的矩阵表示;第一节 基本量及基本方程的矩阵表示;第二节 有限单元法的概念 ;一、基本思路 本节将用结构力学方法的思路介绍有限元法。 有限元法，简单的说，就是用结构力学的方法求解弹性力学问题。首先将连续体变换成为离散化结构，然后再用结构力学方法进行求解。 有限元法的基本思想： 假想的把一连续体分割成数目有限的小体（单元），彼此间只在数目有限的指定点（结点）出相互连结，组成一个单元的集合体以代替原来的连续体，再在结点上引进等效力以代替实际作用于单元上的外力。选择一个简单的函数来近似地表示位移分量的分布规律，建立位移和节点力之间的关系。;经典解与有限元解的区别：;第二节 有限单元法的概念 ;有限元法解题的基本步骤;; 将连续体划分为有限多个有限大小的单元，这些单元仅在一些结点连接起来，构成一个所谓的离散化结构。 根据题目的要求，可选择适当的单元把结构离散化。对于平面问题，最简单而常见的单元式三角形单元，此外还有矩形单元、任意四边形单元。; （1）取各结点位移 为基本未知数 （2）对每个三角形单元，由结点位移 求出单元中的位移分布函数 。 N---形函数矩阵 （3）由单元的位移分布 ，通过几何方程，求出单元的应变列阵，得 ; （4） 由单元的应变 ，通过物理方程，求出单元的应力列阵，得 S—应力转换矩阵 （5）由单元的应力 ，通过虚功方程，求出单元的节点力列阵。 以三角形单元ijk为例：一是结点力 。二是单元ijk内部的应力 。由虚功方程得， ，由此导出 k----单元劲度矩阵 ;;;;;第三节 单元的位移模式与解答的收敛性 ; 在有限元单元法的分析中，首先我们将整个连续弹性体变换为离散化结构，即把一个连续的弹性体变换为一个离散的结构物。但每个单元仍然是一个连续的、均匀的、各向同性的完全弹性体。 对于每个单元，只要求得单元的位移函数，就可以应用几何方程、再应用物理方程求得应力。在有限元法中，结点位移是基本未知量。为此，假定一个位移模式，用来表示单元中的位移函数，此函数在单元结点上等于结点的位移值。 位移模式是根据结点位移值在单元中作出的位移插值函数。 ; 在有限元分析中，我们可以同时使用各种不同形状的单元来离散复杂结构。对于平面问题，三角形单元是最简单、也是最常用的单元；在平面应力问题中，常用单元为三角形板单元；而在平面应变问题中，则常用三棱柱单元。;例如：如图所示的悬臂梁受理问题 分析采用三角形单元，把弹性体划分为有限个互不重叠的三角形。这些三角形在其顶点（即节点）处互相连接，组成一个单元集合体，以替代原来的弹性体。同时，将所有作用在单元上的载荷（包括集中载荷、表面载荷和体积载荷），都按虚功等效的原则移置 ;; 首先，我们来分析一下三角形单元的力学特性，即建立以单元节点位移表示单元内各点位移的关系式。设单元e的节点编号为i、j、m，如图所示。由弹性力学平面问题可知，每个节点在其单元平面内的位移可以有两个分量，所以整个三角形单元将有六个节点位移分量，即六个自由度。用列阵可表示为：;? ? 图 平面三角形单元; 从弹性力学平面问题的解析解法中可知，如果弹性体内的位移分量函数已知，则应变分量和应力分量也就确定了。但是，如果只知道弹性体中某几个点的位移分量的值，那么就不能直接求得应变分量和应力分量。因此，在进行有限元分析时，必须先假定一个位移模式。由于在弹性体内，各点的位移变化情况非常复杂，很难在整个弹性体内选取一个恰当的位移函数来表示位移的复杂变化，但是如果将整个区域分割成许多小单元，那么在每个单元的局部范围内就可以采用比较简单的函数来近似地表示单元的真实位移，将各单元的位移式连接;起来，便可近似地表示整个区域的真实位移函数。这种化繁为简、联合局部逼近整体的思想，正是有限单元法的绝妙之处。;(c);? 平面三角形单元;其中;