02平面问题有限元解法.pptVIP

平面问题的有限单元解法;基本思想;有限元单元模型中几个重要概念;1. 研究内容;弹性力学的内容及基本假定;弹性力学的内容及基本假定;弹性力学的内容及基本假定;弹性力学的内容及基本假定;基本概念：;(2) 面力;2. 应力;(2) 一点的应力状态;用矩阵表示：;3. 形变;(2) 一点应变状态; 工程力学问题建立力学模型的过程中，一般从三方面进行简化：;基本变量 u ε σ （位移） （应变） （应力） 基本方程 力的平衡方程 几何方程 物理方程 边界条件 求解方法 经典解析 半解析 传统数值解法 现代数值解法（计算机硬件、规范化、标准化、规模化）;平面问题的基本理论;平面应力问题;x;平面应变问题;水坝; 如图所示三种情形，是否都属平面问题？是平面应力问题还是平面应变问题？;三大基本方程;平衡微分方程;几何方程;物理方程;假定已知任一点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