有限元概述及一维弹性问题的有限元法.ppt

第二部分一维线弹性体的有限元分析 求解问题： 节点上施加的外载荷，通过对刚度矩阵求逆矩阵的方式，我们能够得到各个节点的位移值。 以上获得的整体刚度矩阵是一个奇异矩阵，即矩阵的行列式为0。 解决办法：加边界条件 第二部分一维线弹性体的有限元分析 第二部分一维线弹性体的有限元分析 第二部分一维线弹性体的有限元分析 高斯消去法 第二部分一维线弹性体的有限元分析 通过第三个方程的求解，可知u4=-0.14186mm。将此数值带入到第二个公式中则可获得u3=-0.15232mm。将这两个值带入第一个公式中得u2=-0.02325mm。 获得了以上的位移结果之后，可以通过原来的未经简化的刚度矩阵计算出相应的另外了两个节点载荷的大小。 如针对F1，应用矩阵方程中第一个方程，有 25 u1-25 u2=F1 其中u1=0，u2=-0.02325，则F1=0.5183N。 同理可以计算出F5=1.4186N。实际上在计算出前4个载荷的条件下，根据系统受力平衡可以直接计算出F5。 第三部分一维简单连续体的分析实例 第三部分一维简单连续体的分析实例 从弹簧跨越到一维弹性体 一维、受力方向与其截面相垂直、载荷状态 第三部分一维简单连续体的分析实例 单元的矩阵方程 第三部分一维简单连续体的分析实例 圆锥台结构，其垂直受力的截面为圆形，在其一端施加载荷，分析结构内部的应力应变情况。 第三部分一维简单连续体的分析实例 每段长度为2m，同时假设每一段的结构的截面积为相应的上下两个截面的和的一半 结构的有限元网格划分 第三部分一维简单连续体的分析实例 相关单元的截面积分别如下 A(1)＝0.6599m2； A(2)＝0.4418m2； A(3)＝0.2672m2。 第三部分一维简单连续体的分析实例 第三部分一维简单连续体的分析实例 计算结果： u2＝3.032?10-5m； u3＝7.561?10-5m； u4＝15.049?10-5m 第三部分一维简单连续体的分析实例 计算精度讨论 有限元分析结果： 理论分析结果： 第三部分一维简单连续体的分析实例 把结构分成10个单元，即每个单元的宽度为0.6m，则相应的有限元伸长量计算结果为15.25?10-5m，此时计算结果误差为0.19％。 最新联系方式 杨建国 浙江工业大学 化工机械设计研究所 杭州，310014 电话：0517 手机:yangjianguo@hit.edu.cn yangjg@zjut.edu.cn 有限元概述及一维弹性问题的有限元法 主要参考书目 李人宪.有限元法基础.国防工业出版社 汪建华.焊接数值模拟技术及其应用.上海交通大学出版社 武传松.焊接热过程数值模拟.哈尔滨工业大学出版社 王勖成. 有限单元法. 清华大学出版社 * 第一部分有限单元方法概述 有限单元方法概述解析方法的局限性 工程技术领域中的许多力学问题和场问题，如固体力学中的位移场、应力场分析、电磁学中的电磁分析、振动特性分析、热力学中的温度场分析，流体力学中的流场分析等，都可以归结为在给定边界条件下求解其控制方程的问题。 虽然人们能够得到它们的基本方程和边界条件，但是能够用解析法求解的只是少数性质比较简单和边界比较规则的问题，实际结构的形状和所受到的载荷往往比较复杂，按解析法求解是非常困难的。 如何解决 有限单元方法概述解决这类问题的两个途径 引入简化假设 将方程、结构几何形状和边界条件简化，使之可能进行解析求解（如：材料力学、结构力学、应用弹性力学的内容）。但其可应用情况有限，并且，简化导致较大偏差。 保留问题的复杂性，数值方法求得近似解 近六十年来，数值方法已成为求解科学技术和工程问题的主要工具。它是在现代数学、力学理论的基础上，借助于计算机技术来获得满足工程要求的数值近似解，是现代工程仿真学发展的重要推动力之一。 有限单元方法概述工程领域常用的数值模拟方法 1、有限单元法FEM（ Finite Element Method） 2、边界元法BEM（Boundary Element Method） 3、有限差分法FDM（ Finite Difference Method） 4、离散单元法DEM（Discrete Element Method） 其中有限单元法是最具实用性和应用最广泛的。 有限单元方法概述有限单元法FEM 虽然近些年才采用了有限元这个名字，有限元的概念在几个世纪以前就已经用过了。例如：古代数学家用多边形逼近圆的办法求出圆周长以及圆的面积；现在人们日常生活中丈量土地的时候也是分成一