工学平面问题有限元解法公式推导ppt课件.ppt

平面问题的有限单元解法;有限元单元法基本思想;基本变量 u ε σ （位移） （应变） （应力） 基本方程 力的平衡方程 几何方程 物理方程 求解方法 经典解析 半解析 传统数值解法 现代数值解法（计算机硬件、规范化、标准化、规模化）;有限元单元模型中几个重要概念;平面问题有限单元法基本概念;有限元单元法分析步骤（一）;有限元单元法分析步骤（二）;有限元单元法分析步骤（三）;弹性力学中的几个基本概念;弹性力学中的几个基本概念;弹性力学中应力的方向规定;弹性力学中的基本假定;平面问题的基本理论;由于板很薄，外力不沿厚度变化，应力沿板的厚度又是连续分布的，所以可以认为在整个薄板的所有各点： 只剩下平行于xy面的三个平面应力分量，即： 这种问题成为平面应力问题。;只剩下平行于xy面的三个形变分量，即： 这种问题成为平面应变问题。;三大基本方程;平衡微分方程;假定已知任一点P处坐标面上的应力分量σx，σ y ，τx y = τ y x 。求经过该点的，平行于z轴而倾斜于x轴和 y轴的任何倾斜面上应力。 从在P点附近取一个平面AB，它平行于上述斜面，并经过P点划出一个微小的三棱柱PAB。当AB无限小而趋于P点时，平面AB上的应力就成为斜面上的应力。;几何方程;物理方程;边界条件;圣维南原理;圣维南原理的应用;圣维南原理的应用(续);圣维南原理的应用(续);有限单元法中基本量的矩阵表示;弹性力学中基本方程的矩阵表示;虚位移原理;虚位移原理(续);（1）取三角形单元的结点位移为基本位置量： (a) 其中， 称为单元的结点位移列阵； （2）应用插值公式，由单元结点位移求出单元的位移函数： (b) 其中，N称为形函数矩阵； （3）应用几何方程，由单元的结点位移求出单元的应变： (c) 其中，B是表示 与 之间关系的矩阵； ; (f) 其中，Fe 是单元的结点力， k称为单元劲度列阵； ;（7）列出各结点的平衡方程，组成整个结构的平衡方程组。由于节??i受有环绕节点的单元移置而来的节点载荷 和节点力 因而i节点的平衡方程为： （i=1,2,…,n） (h);对三角形三个结点i,j,m结点，位移函数应当等于该节点的位移值，即： ;三角形单元的位移模式;三角形单元的位移模式;简写为： 其中矩阵B可写成分块形式： 其子矩阵为：;将D表达式(9)和B表达式(27)代入上式，并写成分块形式，即得到平面应力问中的应力转换矩阵：;由式(26)引起的虚应变为： 由于结点力在虚位移上的虚功应当等于应力在虚应变上的虚功，即： ;从而建立了单元结点力和结点位移之间的关系。对于三角形单元，B中的元素为常量。 并且 ，因此，k可简写为： k称为单元的劲度矩阵。;载荷向节点移置，单元的载荷列阵;载荷向节点移置，单元的载荷列阵（续）;载荷向节点移置，单元的载荷列阵（续）;载荷向节点移置，单元的载荷列阵（续）;注意：式(46)和(48）中的编码i,j,m仅是每个单元的局部编码，对于整个结构，则将结点的平衡方程按整体结点编码1，2，…，n排列起来，就组成整个结构的结点平衡方程组： ;由整体平衡方程组，解出结点位移δ，便可由式(23)和(30)求出每个单元的位移函数和应力。;平面有限元解法（例）;平面有限元解法——划分单元;平