有限元法基础知识介绍课件

有限元法基有限元法基础础知知识识介介绍绍1A有限元的基本思想2A有限元法中的几个基本概念有限元法是把要分析的连续体假想地分割成有限个单元所组成的组合体，简称离散化离散化。这些单元仅在顶角处相互联接，称这些联接点为结结点点。离散化的组合体与真实弹性体的区别在于：组合体中单元与单元之间的联接除了结点之外再无任何关联。但是这种联接要满足变形协调条件，即不能出现裂缝，也不允许发生重叠。显然，单单元之元之间间只能通只能通过结过结点来点来传递传递内力内力。通过结点来传递的内力称为结结点力点力，作用在结点上的荷载称为结结点荷点荷载载。当连续体受到外力作用发生变形时，组成它的各个单元也将发生变形，因而各个结点要产生不同程度的位移，这种位移称为结结点位移点位移。3A有限元法中的几个基本概念在有限元中，常以结点位移作为基本未知量。并对每个单元根据分块近似的思想，假设一个简单的函数近似地表示单元内位移的分布规律，再利用力学理论中的变分原理或其他方法，建立结点力与位移之间的力学特性关系，得到一组以结点位移为未知量的代数方程，从而求解结点的位移分量。然后利用插值函数确定单元集合体上的场函数。显然，如果单元满足问题的收敛性要求，那么随着缩小单元的尺寸，增加求解区域内单元的数目，解的近似程度将不断改进，近似解最终将收敛于精确解。4A有限元法的基本计算步骤n物体离散化n单元特性分析分析单元的力学性质选择位移模式计算等效节点力n单元组集n求解未知节点位移5A有限元法的基本计算步骤1.物体离散化单元选择：应根据连续体的形状选择最能完满地描述连续体形状的单元。常见的单元有：杆单元，梁单元，三角形单元，矩形单元，四边形单元，曲边四边形单元，四面体单元，六面体单元以及曲面六面体单元等等。单元划分：进行单元划分，单元划分完毕后，要将全部单元和结点按一定顺序编号，每个单元所受的荷载均按静力等效原理移植到结点上，并在位移受约束的结点上根据实际情况设置约束条件。6A有限元法的基本计算步骤1.物体离散化3维实体的4面体单元划分平面的三角形单元划分3维实体的6面体单元划分7A有限元法的基本计算步骤1.物体离散化n将某个工程结构离散为由各种单元组成的计算模型，这一步又称作单元剖分或网格划分。n离散后单元于单元之间利用单元的节点相互连接起来；单元节点的设置、性质、数目等应视问题的性质，描述变形形态的需要和计算进度而定。n用有限元分析计算所获得的结果只是近似的。如果划分单元数目非常多而又合理，则所获得的结果就与实际情况相符合。n网格划分的好坏将直接影响到计算结果的准确性和计算进度，甚至会因为网格划分不合理而导致计算不收敛。n网格的划分主要取决于专业知识和经验积累。一个水平高的FEA工程师，80%的时间是用在网格划分上。n对于一般的问题，各种FEA均能自动的进行合理的网格划分。hypermesh最目前好的划分网格工具。8A有限元法的基本计算步骤2.单元特性分析分析单元的力学性质根据单元的材料性质、形状、尺寸、节点数目、位置及其含义等，找出找出单单元元节节点力和点力和节节点点位移的关系式位移的关系式，这是单元分析中的关键一步。此时需要应用弹性力学中的几何方程和物理方程来建立力和位移的方程式，从而导出单元刚度矩阵，这是有限元法的基本步骤之一。9A有限元法的基本计算步骤2.单元特性分析分析单元的力学性质平面问题的三角形单元的例子：单元有三个结点I、J、M，每个结点有两个位移u、v和两个结点力U、V。节点位移节点力取决于材料性质、形状、尺寸10A有限元法的基本计算步骤2.单元特性分析选择位移模式：在反映力和位移的关系式中，依据那一个量是未知量，可建立不同的模型。位移法：选择节点位移作为基本未知量称为位移法；力法：选择节点力作为基本未知量时称为力法；混合法：取一部分节点力和一部分节点位移作为基本未知量时称为混合法。位移法易于实现计算自动化，所以，在有限单元法中位移法应用范围最广。11A有限元法的基本计算步骤2.单元特性分析计算等效节点力：将外在的负载力等效到各个节点上。物体离散化后，假定力是通过节点从一个单元传递到另一个单元。但是，对于实际的连续体，力是从单元的公共边传递到另一个单元中去的。因而，这种作用在单元边界上的表面力、体积力和集中力都需要等效的移到节点上去，也就是用等效的节点力来代替所有作用在单元上得力。弹性体有限元模型12A有限元法的基本计算步骤3.单元组集利用结构力的平衡条件和边界条件把各个单元按原来的结构重新连接起来，形成整体的有限元方程。对由各个单元组成的整体进行分析，建立节点外载荷与结点位移的关系，以解出节点位移，这个过程为整体分析。i节点的节点力：i节点的平衡方程：集中力单元节点力13A有限元法的基本计算步骤3.单元组集 n最终，将所有单元组合起来得到整体的方程：K=RK整体刚度矩阵；全部结点位移组成的列阵；R全部结点荷载组成的列阵。在位移法中，只有是未知的，求解该线性方程组就可得到各结点的位移。将结点位移代入相应方程中可求出单元的应力分量。n有限元法不仅可以求结构体的位移和应力，还可以对结构体进行稳定性分析和动力分析。例如，结构体的整体动力方程：M+C+K=FM整体质量矩阵；C整体阻尼矩阵；K整体刚度矩阵；整体结点位移向量；F整体结点荷载向量。求出结构的自激振动频率、振型等动力响应，以及动变形和动应力等。14A有限元法的基本计算步骤4.求解未知节点位移 n可以根据方程组的具体特点来选择合适的计算方法。n传统有限元分析的数值计算方法之中，有直接计算法（DirectSolver）与迭代法（Iterative）两种。由于在过去的经验中，迭代法一直无法直接而有效的保证数值计算的收敛性，因此，直接计算法在多数有限元素分析软件中，仍然是一种主流的计算方法。15A有限元法的有限元法的过过程程 16A线线性静力分析的基本矩性静力分析的基本矩阵阵方程方程单单元元刚刚度矩度矩阵阵 K=刚度矩阵F=力向量(已知)u=由F引起的未知位移向量 17A总总体体刚刚度矩度矩阵阵18A线线性静力有限元分析步性静力有限元分析步骤骤 结结构离散构离散形成形成单单元元刚刚度矩度矩阵阵装配装配总总体体刚刚度矩度矩阵阵施加施加约约束束边边界条件界条件施加载荷求解方程求解方程计计算位移、算位移、应应力等力等19AF=10000 磅L=10 英寸YINL=10 英寸面积 A=1.5 英寸2弹性模量 E=10106 磅/英寸2面积 A=2.0 英寸2弹性模量 E=10106 磅/英寸2例子：例子：1 1、建立、建立结结构有限元模型构有限元模型 L=10英寸A=1.5英寸2L=10英寸A=2.0英寸2结点 3#结点 1#结点 2#单元 单元 20A2 2、形成、形成单单元元刚刚度矩度矩阵阵 3 3、总总装装刚刚度矩度矩阵阵 21A4 4、施加、施加边边界条件界条件 5 5、施加作用、施加作用载载荷荷 22A6 6、求解矩、求解矩阵阵方程方程 7 7、计计算算单单元元应应力力23A有限元分析方法的发展n在大力推广CAD技术的今天，从自行车到航天飞机，所有的设计制造都离不开有限元分析计算，FEA在工程设计和分析中将得到越来越广泛的重视。n国际上早在20世纪50年代末、60年代初就投入大量的人力和物力开发具有强大功能的有限元分析程序。其中最为著名的是由美国国家宇航局（NASA）在1965年委托美国计算科学公司和贝尔航空系统公司开发的NASTRAN有限元分析系统。该系统发展至今已有几十个版本，是目前世界上规模最大、功能最强的有限元分析系统。24A有限元分析方法的发展与应用目前应用较多的通用有限元软件如下表所列：另外还有许多针对某类问题的专用有限元软件，例如金属成形分析软件Deform、Autoform，焊接与热处理分析软件SysWeld等。软件名称简介MSC/Nastran著名结构分析程序，最初由NASA研制MSC/Dytran动力学分析程序MSC/Marc非线性分析软件ANSYS通用结构分析软件ADINA非线性分析软件ABAQUS非线性分析软件25A有限元方法的应用转向机构支架的强度分析（刘道勇，东风汽车工程研究院动，用MSC/Nastran完成）26A有限元方法的应用金属成形过程分析（用Deform软件完成）:分析金属成形过程中的各种缺陷型材挤压成形的分析。型材在挤压成形的初期，容易产生形状扭曲螺旋齿轮成形过程的分析T形锻件的成形分析27A有限元方法的应用焊接残余应力分析（用Sysweld完成）结构与焊缝布置焊接过程的温度分布与轴向残余应力28A有限元方法的应用BMW曲轴的感应淬火（用SysWeld软件完成）29A有限元方法的应用复杂形状工件的组织转变预测（石伟，用NSHT3D完成）：预测工件的组织分布和机械性能二分之一工件的有限元模型淬火3.06min时的温度分布淬火3.06min时的马氏体分布30A有限元方法的应用其他的应用还包括：电磁学、流体力学、电磁场等等31A有限元方法的发展趋势l从从单纯结单纯结构力学构力学计计算算发发展到求解展到求解许许多物理多物理场场问题问题有限元分析方法最早是从结构化矩阵分析发展而来，逐步推广到板、壳和实体等连续体固体力学分析，实践证明这是一种非常有效的数值分析方法。有限元方法已发展到流体力学、温度场、电传导、磁场、渗流和声场等问题的求解计算，最近又发展到求解几个交叉学科的问题。例如当气流流过一个很高的铁塔产生变形，而塔的变形又反过来影响到气流的流动这就需要用固体力学和流体动力学的有限元分析结果交叉迭代求解，即所谓流固耦合的问题。32A有限元方法的发展趋势n由求解线性工程问题进展到分析非线性问题线性理论已经远远不能满足设计的要求。例如：航天和动力工程的高温部件存在热变形和热应力，要考虑材料的非线性问题；诸如塑料、橡胶和复合材料等各种新材料的出现，只有采用非线性有限元算法才能解决。非线性的数值计算是很复杂的，很难为一般工程技术人员所掌握。为此近年来国外一些公司花费了大量的人力和投资开发诸如MARC、ABQUS和ADINA等专长于求解非线性问题的有限元分析软件，并广泛应用于工程实践。33A有限元方法的发展趋势n增强可视化的前置建模和后置数据处理功能随着数值分析方法的逐步完善，尤其是计算机运算速度的飞速发展，整个计算系统用于求解运算的时间越来越少，而数据准备和运算结果的表现问题却日益突出。在现在的工程工作站上，求解一个包含10万个方程的有限元模型只需要用几十分钟。工程师在分析计算一个工程问题时有80%以上的精力都花在数据准备和结果分析上。34A有限元方法的发展趋势n增强可视化的前置建模和后置数据处理功能目前几乎所有的商业化有限元程序系统都有功能很强的前置建模和后置数据处理模块。使用户能以可视图形方式直观快速地进行网格自动划分，生成有限元分析所需数据，并按要求将大量的计算结果整理成变形图、等值分布云图，便于极值搜索和所需数据的列表输出。35A有限元方法的发展趋势n与CADCAD软件的无缝集成当今有限元分析系统的另一个特点是与通用CAD软件的集成使用，即：在用CAD软件完成部件和零件的造型设计后，自动生成有限元网格并进行计算，如果分析的结果不符合设计要求则重新进行造型和计算，直到满意为止，从而极大地提高了设计水平和效率。当今所有的商业化有限元系统商都开发了和著名的CAD软件（例如Pro/ENGINEER、Unigraphics、SolidEdge、SolidWorks、IDEAS、Bentley 和AutoCAD 等）的接口。36A有限元方法在设计中的作用n以有限元分析为典型代表的CAE，成为缩短设计流程，减少错误设计从而降低设计和制造的成本的有效途径。在现代化的工业中，计算机辅助设计（CAD）,计算机辅助制造（CAM），已经减少了不少设计者的负担。而原本被视为设计和制造过程中的配角-计算机辅助工程（CAE），已经摆脱了以前可有可无的角色，变成设计过程中不可缺少的一环。计算机辅助工程（CAE）在产业中肩负着五大使命：（1）缩短设计所需的时间和降低设计成本。（2）在精确的分析结果下制造出高质量的产品。（3）能够快速的对设计变更作出反应。（4）能充分的和CAD结合并对不同类型的问题进行分析。（5）能够精确的预测出产品的性能。37ACAE所包含的内容n设计验证：设计是否可行？设计是否按照我们期望的方式运行？n方案择优：那一种设计方案是所有方案中最好的？剔除次要的方案n概念验证在不作实体的情况下验证新的设计思想。设计过程可以跳出各种条框的限制。n失效分析分析方法最传统的用途之一；需要一定水平的专业知识。38A有限元软件使用的基本步骤n三个步骤：前处理（PREPROCESSION），求解（SOLUTION），后处理POSTPROCESSION）n前处理单元选择材料定义几何建模网格划分模型局部调整施加荷载39A有限元软件使用的基本步骤n求解分析问题的类型设定分析参数添加荷载条件建立荷载工况求解40A有限元软件使用的基本步骤n后处理结果的文字输出(ResultList)结果的云图输出(ResultContour)结果的矢量输出(ResultVector)结果的路径输出(ResultMapping)ElementTable的提取LoadCase及其组合41AEnd42A