有限元方法课件第四章平面三角形单元.pptxVIP

第四章 平面三角形单元第一页，共九十页。第四章 平面三角形单元 §4–1有限元法的基本思想 §4–2三角形常应变单元 §4–3形函数的性质 §4–4刚度矩阵 §4–5等效节点力载荷列阵 §4–6有限元分析的实施步骤 §4–7计算实例第二页，共九十页。§4-1 有限元法的基本思想一、有限元法的基本思想 假想的把一连续体分割成数目有限的小体（单元），彼此间只在数目有限的指定点（结点）出相互连结，组成一个单元的集合体以代替原来的连续体，再在结点上引进等效力以代替实际作用于单元上的外力。选择一个简单的函数来近似地表示位移分量的分布规律，建立位移和节点力之间的关系。 有限元法的实质是：把有无限个自由度的连续体，理想化为只有有限个自由度的单元集合体，使问题简化为适合于数值解法的结构型问题。第三页，共九十页。 微分 数目增到∞ 建立一个描述连续体经 典 解 法 ——（解析法） 大小趋于 0 性质的偏微分方程 有限单元 离散化 集合 总体分析解有限元法——连续体——单元——代替原连续体（近似法） （单元分析） 线性方程组二、经典解与有限元解的区别：第四页，共九十页。yx三、有限元法算题的基本步骤1. 力学模型的选取（平面问题，平面应变问题，平面应力问题，轴对称问题，空间问题，板，梁，杆或组合体等，对称或反对称等）例如：为平面应力问题，由于结构的对称性可取结构的1/4来研究，故所取的力学模型第五页，共九十页。 根据题目的要求，可选择适当的单元把结构离散化。对于平面问题可用三角元，四边元等。2. 单元的选取、结构的离散化例如：第六页，共九十页。——单元内任一点的位移列阵；——单元的结点位移列阵；——单元的形函数矩阵；（它的元素是任一点位置坐　　　　　标的函数）3. 选择单元的位移模式　结构离散化后，要用单元内结点的位移通过插值来获得单元内各点的位移。在有限元法中，通常都是假定单元的位移模式是多项式，一般来说，单元位移多项式的项数应与单元的自由度数相等。它的阶数至少包含常数项和一次项。至于高次项要选取多少项，则应视单元的类型而定。（4-1） 第七页，共九十页。——单元内任一点应变列阵；——单元的应变矩阵；（它的元素仍为位置坐标的　 函数）4. 单元的力学特性分析 把（4-1）式代入几何方程可推导出用单元结点位移表示的单元应变表达式：（4-2）式中： 再把（4-２）式代入物理方程，可导出用单元结点位移列阵表示的单元应力表达式：（4-3）第八页，共九十页。——单元内任一点的应力列阵；——单元的弹性矩阵，（它与材料的特性有关）式中：　最后利用弹性体的虚功方程建立单元结点力阵与结点位移列阵之间的关系，即形成单元的刚度方程式：（4-4）——单元刚度矩阵式中：（4-5）第九页，共九十页。 用直接刚度法将单刚　　组集成总纲　　，并将　　组集成总载荷列阵　　，形成总体结构的刚度方程：　解出整体结构的结点位移列阵　　后，再根据单元结点的编号找出对应于单元的位移列阵　　，将　　代入（4-3）式就可求出各单元的应力分量值。5. 建立整体结构的刚度方程（4-6）　考虑整体结构的约束情况，修改整体刚度方程之后，（4-6）式就变成以结点位移为未知数的代数方程组。解此方程组可求出结点位移。6. 求解修改后的整体结构刚度方程7. 由单元的结点位移列阵计算单元应力第十页，共九十页。8. 计算结果输出 求解出整体结构的位移和应力后，可有选择地整理输出某些关键点的位移值和应力值，特别要输出结构的 变形图、应力图、应变图、结构仿真变形过程动画图及整体结构的弯矩、剪力图等等。第十一页，共九十页。§4-2 三角形常应变单元一、离散化 在运用有限单元法分析弹性力学平面问题时，第一步就是要对弹性体进行离散化，把一个连续的弹性体变换为一个离散的结构物。对于平面问题，三角形单元是最简单、也是最常用的单元，在平面应力问题中，单元为三角形板，而在平面应变问题中，则是三棱柱。 假设采用三角形单元，把弹性体划分为有限个互不重叠的三角形。这些三角形在其顶点（即节点）处互相连接，组成一个单元集合体，以替代原来的弹性体。同时，将所有作用在单元上的载荷（包括集中载荷、表面载荷和体积载荷），都按虚功等效的原则移置到节点上，成为等效节点载荷。由此便得到了平面问题的有限元计算模型，如图4-1所示。 第十二页，共九十页。??图4-1 弹性体和有限元计算模型第十三页，共九十页。??图4-2 平面三角形单元第十四页，共九十页。二、位 移 首先，我们来分析一下三角形单元的力学特性，即建立以单元节点位移表示单元内各点位移的关系式。设单元e的节点编号为i、j、m，如图4-2所示。由弹性力学平面问题可知，每个节点在其单元平面内的位移可以