平面问题有限元解法（公式推导讲解）课件.ppt

南京农业大学工学院机械工程系 平面问题的有限单元解法 有限元单元法基本思想 有限单元法的思想是将物体（连续的求解域）离散成有限个且按一定方式相互联结在一起的单元组合，来模拟或逼近原来的物体，从而将一个连续的无限自由度问题简化为离散的有限自由度问题求解的一种数值分析法。物体被离散后，通过对其中各个单元进行单元分析，最终得到对整个物体的分析。 有限单元法的分析步骤如下： 物体离散化 单元特性分析 单元组集，整体分析 求解未知节点的位移 由节点的位移求解各单元的位移和应力 物体变形及受力情况的描述 基本变量 u ε σ （位移） （应变） （应力） 基本方程 力的平衡方程 几何方程 物理方程 求解方法 经典解析 半解析 传统数值解法 现代数值解法（计算机硬件、规范化、标准化、规模化） 有限元单元模型中几个重要概念 单元 网格划分中每一个小的块体 节点 确定单元形状、单元之间相互联结的点 节点力 单元上节点处的结构内力 载荷 作用在单元节点上的外力 （集中力、分布力） 约束 限制某些节点的某些自由度 弹性模量（杨式模量）E 泊松比（横向变形系数）μ 密度 平面问题有限单元法基本概念 有限单元法(FEM)是20世纪50年代以来随着计算机的广泛应用而发展起来的一种数值解法。简单地说，就是用结构力学方法求解弹性力学问题。 平面问题的有限单元法求解 将连续体变换成为离散化结构。即将连续体划分为有限多个有限大小的单元，这些单元仅在一些结点连接起来，构成一个所谓离散化结构。（对于平面问题，常用的单元是三角形单元） 用结构力学方法进行求解 有限元单元法分析步骤（一） 结构离散化 将结构分成有限个小的单元体，单元与单元、单元与边界之间通过节点连接。结构的离散化是有限元法分析地第一步，关系到计算精度和效率，包括以下三个方面： 单元类型的选择。选定单元类型，确定单元形状、单元节点数、节点自由度数等。 单元划分。网格划分越细，节点越多，计算结果越精确，但计算量越大。网格加密到一定程度后计算精度提高就不明显，对应应力变化平缓区域不必要细分网格。 节点编码。 注意：有限元分析的结构已不是原有的物体或结构物，而是由同样材料、众多单元以一定方式连接成的离散物体。所以，用有限元分析计算所获得的结果是近似的（满足工程要求即可）。 有限元单元法分析步骤（二） 单元特性分析 选择未知量模式 选择节点位移作为基本未知量时，称为位移法； 选节点力作为基本未知量时，称为力法； 取一部分节点位移和一部分节点力作为未知量，称为混合法。 分析单元力学性质 根据单元材料性质、形状、尺寸、节点数目、位置等，找出单元节点力和节点位移关系式，应用几何方程和物理方程建立力和位移的方程式，从而导出单元刚度矩阵。 计算等效节点力 作用在单元边界上的表面力、体积力或集中力都需要等效地移到节点上去，即用等效力来替代所有作用在单元上的力。 有限元单元法分析步骤（三） 整体分析 集成整体节点载荷矢量 F 。结构离散化后，单元之间通过节点传递力，作用在单元边界上的表面力、体积力或集中力都需要等效地移到节点上去，形成等效节点载荷。将所有节点载荷按照整体节点编码顺序组集成整体节点载荷矢量。 组成整体刚度矩阵K ，得到总体平衡方程： 引进边界约束条件，解总体平衡方程求出节点位移。 通过上述分析可以看出有限单元法的基本思想是“一分一合”，分是为了进行单元分析，合是为了对整体的结构进行综合分析。 弹性力学中的几个基本概念 弹性力学中的几个基本概念 弹性力学中应力的方向规定 每一个面上的应力可以分解为一个正应力和两个切应力。 正应力用σ表示，加上一个下标字母，表示作用面和作用方向。 切应力用τ表示，并加上两个下标字母，表示作用面和作用方向。前一个字母表示作用面垂直于哪一个坐标轴，后一个字母表示作用方向沿着哪一个坐标轴。 弹性力学中的基本假定 连续性——假定整个物体的体积都被组成这个物体的介质所填满，不留任何空隙。 完全弹性——假定物体在引起形变的外力被除去之后能恢复原形，而没有任何剩余形变。 均匀性——假定整个物体有同一材料组成的，物体的所有各部分具有相同的弹性。 各向同性——假定物体的弹性在所有各个方向都相同。 小变形——假定位移和形变是微小的，物体受力之后，整个物体所有各点的位移都远远小于物体原来的尺寸，因而应变和转角都远小于1。 平面问题的基本理论 任何一个实际的弹性力学问题都是空间问题，但是如果所考察的弹性体具有某种特殊的形状，并且承受的是某些特殊的外力和约束，就可以把空间问题简化为近似的平面问题。 两种典型的平面问题 平面应力问题 平面应变问题 平面应力问题