有限元基础教程(绪论)ppt课件

Finite Element Method and ANSYS程程 强强 有限元分析及有限元分析及ANSYS北京工业大学机电学院北京工业大学机电学院Finite Element Method and ANSY课程设置课程设置英文名称：英文名称：Finite element method and ANSYS课程类型：专业选修课程类型：专业选修学时：学时：32 学分：学分：2适用对象：机械制造及自动化专业本科生适用对象：机械制造及自动化专业本科生先修课程：线性代数、理论力学，材料力学先修课程：线性代数、理论力学，材料力学使用教材及参考书：使用教材及参考书：1.有限元分析基础教程有限元分析基础教程，曾攀，清华大学，曾攀，清华大学，20082.Ansys软件软件课程教学基本要求课程教学基本要求本本课课程程教教学学环环节节主主要要包包括括：课课堂堂讲讲授授（20学学时时）、作作业业、上上机机计计算（算（12学时）。学时）。每章布置每章布置1-2道习题，以巩固所学的理论。道习题，以巩固所学的理论。不单设考试，以大作业的报告、平时作业和考勤综合评定成绩。不单设考试，以大作业的报告、平时作业和考勤综合评定成绩。课程设置及要求课程设置及要求课程设置课程设置及要求绪绪 论论1.1概况概况1.2有限元方法的历史有限元方法的历史1.3有限元分析的内容和作用有限元分析的内容和作用1.4有限元分析的一般过程有限元分析的一般过程1.5有限元法的基本概念有限元法的基本概念1.6有限元法的发展趋势有限元法的发展趋势绪 论1.1概况1概况概况 有限元方法(finite element method)或有限元分析(finite element analysis)，是求取复杂微分方程近似解的一种非常有效的工具，是现代数字化科技的一种重要基础性原理。有限元分析必须包含三个方面：(1)有限元方法的基本数学力学原理；有限元方法的基本数学力学原理；(2)基于原理所形成的实用软件；基于原理所形成的实用软件；(3)使用时的计算机硬件。使用时的计算机硬件。随着现代计算机技术的发展，一般的个人计算机就能满足第(3)方面的要求；因此，本课的重点将在以上的第(1)和第(2)方面，将通过一些典型的实例来深入浅出地系统阐述有限元分析的基本原理，并强调原理的工程背景和物理概念。1概况 有限元方法(finite element m1概况概况有限元法：有限元法：也称为有限单元法或有限元素法，基本思想是将求解区域离散离散为一组有限个、且按一定方式相互连接在一起的单元单元的组合体。由于单元可以被分割成各种形状和大小不同的尺寸，所以它能很好地适应复杂的几何形状、复杂的几何形状、复杂的材料特性和复杂的边界条件复杂的材料特性和复杂的边界条件。它是将弹性理论、计弹性理论、计算数学和计算机软件算数学和计算机软件有机地结合在一起的一种数值计算方法。1概况有限元法：也称为有限单元法或有限元素法，基本思想是将求 任何连续体都可以假想地分割成有限个简单形状单单元元体的组合，在有限元法中将这些简单形状的单元体称为单元，把单元与单元之间设置的相互连接点，称为节节点点，如图1.1所示。从理论上说，单元的分割可以是任意的，不过在实际计算中必须根据研究对象的特点，使单元分割既满足力学分析要求，又能使计算简便。图1.1 单元与节点 任何连续体都可以假想地分割成有限个简单左左图图所所示示,为为分分析析齿齿轮轮上上一一个个齿齿内内的的应应力力分分布布,可可分分析析图图中中所所示示的的一一个个 平平面面截截面面内内位位移移分分布布.作作为为近近似似解解,可可以以先先求求出出图图 中中 各各 三三 角角 形形 顶顶 点点 的的 位位 移移.这这 里里 的的三角形就是单元三角形就是单元,其顶点就是节点。其顶点就是节点。从物理角度理解从物理角度理解从物理角度理解从物理角度理解,可把一个连续的齿形截面单元之间在节可把一个连续的齿形截面单元之间在节点处以点处以铰链相链接铰链相链接,由单元组合而成的结构近似代替原连续结由单元组合而成的结构近似代替原连续结构构,在一定的约束条件下在一定的约束条件下,在给定的载荷作用下在给定的载荷作用下,就可以求出各就可以求出各节点的位移节点的位移,进而求出应力进而求出应力.从数学角度理解从数学角度理解从数学角度理解从数学角度理解,把这个求解区域剖分成许多三角形子域把这个求解区域剖分成许多三角形子域,子域内的位移可用相应各节点的待定位移合理插值来表示子域内的位移可用相应各节点的待定位移合理插值来表示.例如例如例如例如:1概况概况左图所示,为分析齿轮上一个齿内的应力分布,可分析图中所示的一真实系统真实系统有限元模型有限元模型 有限元模型有限元模型 是真实系统理想化的数学抽象是真实系统理想化的数学抽象。有限有限元模型由一些简单形状的元模型由一些简单形状的单元单元组成，组成，单元单元之间通过之间通过节点连接，并承受一定载荷。节点连接，并承受一定载荷。定义定义1概况概况真实系统有限元模型 有限元模型 是真实系统理想化的数学抽象。2有限元方法的历史有限元方法的历史数学力学理论研究方面（国外）：数学力学理论研究方面（国外）：1，1870年，英国科学家Rayleigh采用假想的“试函数试函数”来求解复杂的微分方程，1909年Ritz将其发展成为完善的数值近似方法，为现代有限元方法打下坚实基础。2，1943年，Courant发表了第一篇使用三角形区域三角形区域的多项式函数来求解扭转问题的论文；3，1956年波音公司的研究人员在分析飞机结构时系统研究了离散杆、梁、三角形的单元刚度离散杆、梁、三角形的单元刚度表达式；4，1960年Clough在处理平面弹性问题，第一次提出并使用“有限元方法有限元方法”的名称；5，1955年德国的Argyris出版了第一本关于结构分析中的能能量原理量原理和矩阵方法矩阵方法的书，为后续的有限元研究奠定了重要的基础。2有限元方法的历史数学力学理论研究方面（国外）：2有限元方法的历史有限元方法的历史数学力学理论研究方面（国内）：数学力学理论研究方面（国内）：我国的一些学者也在有限元领域做出了重要的贡献：1.胡海昌于1954提出了广义变分原理广义变分原理；2.钱伟长研究了拉格朗日乘子法与广义变分原理拉格朗日乘子法与广义变分原理之间关系；3.钱令希研究了力学分析的余能原理余能原理；4.冯康独立地、并先于西方奠定了有有限限元元分分析析收收敛敛性性的理论基础。2有限元方法的历史数学力学理论研究方面（国内）：2有限元方法的历史有限元方法的历史有限元软件应用及学术论文：有限元软件应用及学术论文：随着计算机技术的飞速发展，基于有限元方法原理的软件大量出现，并在实际工程中发挥了愈来愈重要的作用；目前，专业的著名有限元分析软件公司有几十家，国际上著名的通用有限元分析软件有ANSYS，ABAQUS，MSC/NASTRAN，MSC/MARC，ADINA，ALGOR，PRO/MECHANICA，IDEAS，还有一些专门的有限元分析软件，如LS-DYNA，DEFORM，PAM-STAMP,AUTOFORM，SUPER-FORGE等；国际上著名的主要有限元分析软件状况见表1-1。有关有限元分析的学术论文，每年也不计其数，学术活动非常活跃，表1-2 列出的是刊登有限元分析论文的常见学术期刊。2有限元方法的历史有限元软件应用及学术论文：刊登有限元分析论文的常见学术期刊刊登有限元分析论文的常见学术期刊1-2刊登有限元分析论文的常见学术期刊1-2国际上著名的主要有限元分析软件状况见表国际上著名的主要有限元分析软件状况见表1-1国际上著名的主要有限元分析软件状况见表1-1有限元方法的发展过程有限元方法的发展过程见见图图1-1有限元方法的发展过程见图1-12有限元方法的历史有限元方法的历史有限元方法的主要创始人：有限元方法的主要创始人：1.Richard Courant(美国数学家),1910年在哥廷根大学获得博士学位,1928年创办了应用数学研究所，并在19281933年担任所长，开始应用数学的研究，1936年去纽约大学创立了数学研究所，1964年该研究所被命名为Courant数学科学研究所，Courant的的名名字字还还与与有有限限元元方方法法紧紧密密联联系系在在一一起起，是是他他给给这这种种数数值值求求解解偏偏微微分分方方程程以以坚坚实实的的理理论论基基础础，“CourantFriedrichsLewy条件”以及“Courant最小值原理”都是以他的名字命名的重要定理。2有限元方法的历史有限元方法的主要创始人：2有限元方法的历史有限元方法的历史有限元方法的主要创始人：有限元方法的主要创始人：2.John Argyris(德国人)被誉为有限元分析的发明人和创始者之一，所创立的力力学学分分析析的的矩矩阵阵位位移移方方法法为为有有限限元元方方法法的前身，被称为世界上最先进的数学工具，先后因在有限元方法以及混沌理论方面的贡献而获得菲利普王子奖章、美国最高荣誉的爱因斯坦奖，获得包括5个G8国家在内的由16发达国家授予的杰出科技与工程奖；Argyris于1950年代在斯图加特大学创立计算机应用研究所，开发了几个著名的在欧洲广泛使用的有限元分析商业化软件，在有限元方法的应用方面也作出了杰出的贡献。2有限元方法的历史有限元方法的主要创始人：2有限元方法的历史有限元方法的历史有限元方法的主要创始人：有限元方法的主要创始人：3.Olgierd Cecil Zienkiewiz(英国人)。在1968年，创办了International Journal for Numerical Methods in Engineering杂志并任主编，该杂志至今仍然是该领域的主要刊物。Zienkiewiez教授还是5所科学院的院士，这是对他在有限元方法领域的奠基性发展和贡献的赞誉。1967年，他他出出版版了了有有限限元元领领域域的的第第一一本本专专著著有有限限元元方方法法，该书目前也出版到第5版，可参见有限元方法第1卷：基本原理)(清华大学出版社，2008)。2有限元方法的历史有限元方法的主要创始人：3有限元分析的作用有限元分析的作用 基于功能完善的有限元分析软件和高性能的计算机硬件对设计的结构进行详细的力学分析，以获得尽可能真实的结构受力信息，就可以在设设计计阶阶段段对可能出现的各种问题进行安安全全评评判判和设设计计参参数数修修改改，据有关资料，一个新产品的问题有60以上可以在设计阶段消除，甚至有的结构的施工过程也需要进行精精细细的的设设计计，要做到这一点，就需要类似有限元分析这样的分析手段。下面举出几个涉及机械工程、土木工程、车辆工程、航空工程以及生物工程的实例。3有限元分析的作用 基于功能完善的有限元分析软件和高性ANSYSANSYS中支中支架计算模型架计算模型ANSYS中支架计算模型ANSYSANSYS中计算中计算模型的网格划模型的网格划分图分图ANSYS中计算模型的网格划分图计算得出的计算得出的支架变形图支架变形图计算得出的支架应力支架应力彩图彩图支架应力支架变形支架变形的动画图的动画图支架变形的动画图立体车库变立体车库变形应力彩图形应力彩图有限元基础教程(绪论)ppt课件这是一个对模型飞机翅膀这是一个对模型飞机翅膀的模态分析的模态分析的例题的例题这是一个对模型飞机翅膀3有限元分析的作用有限元分析的作用 奥运场馆的鸟巢是建设中最艰难的，枝蔓重达42000T，其中，顶盖以及周边悬空部位重量为14000T。在施工时，采用了78根支柱进行支撑，也就是产生了78个受力区域，在钢结构焊接完成后，需要将其缓慢而又平稳地卸去，让鸟巢变成完全靠自身结构支撑；如何卸载？图图1-2 北京奥运场馆北京奥运场馆鸟巢的钢铁枝蔓结构鸟巢的钢铁枝蔓结构支撑塔架支撑塔架3有限元分析的作用 奥运场馆的鸟巢是建设中最艰难的奥运鸟巢的有奥运鸟巢的有限元模型限元模型奥运鸟巢的有限元模型3有限元分析的作用有限元分析的作用图1-3给出的是现代列车车厢整体结构的有限元分析模型3有限元分析的作用图1-3给出的是现代列车车厢整体结构的有限3有限元分析的作用有限元分析的作用图1-4为空客A350后机身第19框的设计与有限元分析过程3有限元分析的作用图1-4为空客A350后机身第19框的设计3有限元分析的作用有限元分析的作用图1-5为人体肩部区域的骨胳有限元分析模型及计算结果 3有限元分析的作用图1-5为人体肩部区域的骨胳有限元分析模型3有限元分析的作用有限元分析的作用图1-6为精密卧式加工中心有限元分析模型及计算结果 3有限元分析的作用图1-6为精密卧式加工中心有限元分析模型及3有限元分析的作用有限元分析的作用图1-7为龙门铣床的有限元分析模型及计算结果 3有限元分析的作用图1-7为龙门铣床的有限元分析模型及计算结4 4 4 4、有限元分析问题的一般过程、有限元分析问题的一般过程、有限元分析问题的一般过程、有限元分析问题的一般过程确定位移模式确定位移模式单元特性分析单元特性分析单元刚度矩阵单元刚度矩阵整体分析整体分析整体刚度矩阵整体刚度矩阵解方程组解方程组输出计算结果输出计算结果选择结果表示形式选择结果表示形式结构离散化结构离散化选择单元类型选择单元类型引入边界引入边界条件条件解释结果解释结果直接平衡法；功和能量法；加权残余法消元法；迭代法叠加法4、有限元分析问题的一般过程确定位移模式单元特性分析整体分析5 5 有限元法的基本概念有限元法的基本概念有限元法的基本概念有限元法的基本概念l结构离散（有限元建模）结构离散（有限元建模）l内容：内容：1）网格划分）网格划分-即把结构按一定规则分割成有限单元即把结构按一定规则分割成有限单元 2）边界处理）边界处理-即把作用于结构边界上约束和载荷处理为即把作用于结构边界上约束和载荷处理为节点约束和节点载荷节点约束和节点载荷l要求：要求：1）离散结构必须与原始结构保形）离散结构必须与原始结构保形-单元的几何特性单元的几何特性 2）一个单元内的物理特性必须相同）一个单元内的物理特性必须相同-单元的物理特性单元的物理特性5 有限元法的基本概念结构离散（有限元建模）单元与节点单元与节点单元与节点单元与节点l单元单元：即原始结构离散后，满足一即原始结构离散后，满足一定几何特性和物理特性的最小结构定几何特性和物理特性的最小结构域。域。l节点节点：单元与单元间的连接点。：单元与单元间的连接点。l节点力节点力：单元与单元间通过节点的：单元与单元间通过节点的相互作用力相互作用力l节点载荷节点载荷：作用于节点上的外载。：作用于节点上的外载。注意：注意：1)节点是有限元法的重要概念，节点是有限元法的重要概念，有限元模型中，相邻单元的作用通有限元模型中，相邻单元的作用通过节点传递，而单元边界不传递力，过节点传递，而单元边界不传递力，这是离散结构与实际结构的重大差这是离散结构与实际结构的重大差别；别；2）节点力与节点载荷的差别。）节点力与节点载荷的差别。节点载荷节点载荷节点力节点力单元与节点单元：即原始结构离散后，满足一定几何特性和物理特性非法结构离散不同材料不同材料节点不合法节点不合法非法结构离散不同材料节点不合法典型单元类型单元类型单元类型单元图形单元图形节点数节点数节点自由度节点自由度杆单元杆单元22梁单元梁单元23平面单元平面单元32平面四边形平面四边形42轴对称问题轴对称问题 32板壳单元板壳单元43四面体单元四面体单元43典型单元类型单元类型单元图形节点数节点自由度杆单元22梁单插值函数（或位移函数）插值函数（或位移函数）插值函数（或位移函数）插值函数（或位移函数）l用以表示单元内物理量变化（如位移或位移场）的近似函数。由于该用以表示单元内物理量变化（如位移或位移场）的近似函数。由于该近似函数常由单元节点物理量值插值构成，故称为近似函数常由单元节点物理量值插值构成，故称为插值函数插值函数，如单元，如单元内物理量为位移，则该函数称为内物理量为位移，则该函数称为位移函数位移函数。l选择位移函数的一般原则：选择位移函数的一般原则：1）位移函数在单元节点的值应等于节点位移（即单元内部是连续的）；）位移函数在单元节点的值应等于节点位移（即单元内部是连续的）；2）所选位移函数必须保证有限元的解收敛于真实解。）所选位移函数必须保证有限元的解收敛于真实解。注注：为了便于微积分运算，位移函数一般采用：为了便于微积分运算，位移函数一般采用多项式多项式形式，在单元内形式，在单元内选取适当阶次的多项式可得到与真实解接近的近似解选取适当阶次的多项式可得到与真实解接近的近似解.插值函数（或位移函数）用以表示单元内物理量变化（如位移或位移位移函数的构造方法（广义坐标法）位移函数的构造方法（广义坐标法）位移函数的构造方法（广义坐标法）位移函数的构造方法（广义坐标法）l广义坐标法广义坐标法 一维单元位移函数：一维单元位移函数：为待定系数，也称为广义为待定系数，也称为广义坐标坐标位移函数的构造方法（广义坐标法）广义坐标法位移函数的构造方法（插值函数法）位移函数的构造方法（插值函数法）位移函数的构造方法（插值函数法）位移函数的构造方法（插值函数法）l插值函数法插值函数法 即将位移函数表示为各个节点位移与即将位移函数表示为各个节点位移与已知插值基函数积的和。已知插值基函数积的和。如一维单元如一维单元l二维单元二维单元注：注：Ni可为可为Lagrange、Hamiton多项式或形函数，多项式或形函数，在在+1-1间变化间变化位移函数的构造方法（插值函数法）插值函数法 即将位移函数6有限元方法的发展趋势有限元方法的发展趋势1、与CAD软件的无缝集成在用CAD软件完成部件和零件的造型设计后，能直直接接将将模模型型传传送送到到CAE软软件件中中进进行行有有限限元元网网格格划划分分并并进进行行分分析析计计算算，如果分析的结果不满足设计要求则重新进行设计和分析，直到满意为止，从而极大地提高了设计水平和效率。为了满足工程师快捷地解决复杂工程问题的要求，许多商业化有限元分析软件都开发了和著名的CAD软件（例如Pro/ENGINEER、Unigraphics、SolidEdge、SolidWorks、IDEAS、Bentley和AutoCAD等）的接口。有些CAE软件为了实现和CAD软件的无缝集成而采用了CAD的建模技术，如ADINA软件由于采用了基于Parasolid内核的实体建模技术，能和以Parasolid为核心的CAD软件（如 Unigraphics、SolidEdge、SolidWorks）实现真正无缝的双向数据交换。6有限元方法的发展趋势1、与CAD软件的无缝集成6有限元方法的发展趋势有限元方法的发展趋势2、更为强大的网格处理能力网格生成的质量和效率都有了很大的提高，但在有些方面却一直没有得到改进，如对三维实体模型进行自动六面体网格划分和根据求解结果对模型进行自适应网格划分。自自适适应应性性网网格格划划分分是是指指在在现现有有网网格格基基础础上上，根根据据有有限限元元计计算算结结果果估估计计计计算算误误差差、重重新新划划分分网网格格和和再再计计算算的的一一个个循循环环过过程程。对于许多工程实际问题，在整个求解过程中，模型的某些区域将会产生很大的应变，引起单元畸变，从而导致求解不能进行下去或求解结果不正确，因此必须进行网网格格自自动动重重划划分分。自适应网格往往是许多工程问题如裂纹扩展、薄板成形等大应变分析的必要 条件。6有限元方法的发展趋势2、更为强大的网格处理能力6有限元方法的发展趋势有限元方法的发展趋势3、由求解线性求解线性问题发展到求解非线性求解非线性问题许多工程问题如材料的破坏与失效、裂纹扩展等仅靠线性理论根本不能解决，必须进行非线性分析求解。例如：薄板成形就要求同时考虑结构的大位移、大应变（几何非线性）和塑性（材料非线性）；而对塑料、橡胶、陶瓷、混凝土及岩土等材料进行分析或需考虑材料的塑性、蠕变效应时则必须考虑材料非线性。众所周知，非线性问题的求解是很复杂的，它不仅涉及到很多专门的数学问题，还必须掌握一定的理论知识和求解技巧，学习起来也较为困难。为此国外一些公司花费了大量的人力和物力开发非线性求解分析软件，如ADINA、ABAQUS等。它们的共同特点是具有高效的非线性求解器、丰富而实用的非线性材料库，ADINA还同时具有隐式和显式两种时间积分方法。6有限元方法的发展趋势3、由求解线性问题发展到求解非线性问题6有限元方法的发展趋势有限元方法的发展趋势4、由单一结构场单一结构场求解发展到耦合场问题耦合场问题的求解有限元分析方法最早应用于航空航天领域，主要用来求解线线性性结结构构问问题题。现在用于求解结构线性问题的有限元方法和软件已经比较成熟，发展方向是结结构构非非线线性性、流流体体动动力力学学和耦合场问题和耦合场问题的求解。例如：由于摩擦接触而产生的热问题，金属成形时由于塑性功而产生的热问题,需要结构场和温度场的有有限限元元分分析析结结果果交交叉叉迭迭代代求求解解，即热热力力耦耦合合的问题。当流体在弯管中流动时，流体压力会使弯管产生变形，而管的变形又反过来影响到流体的流动这就需要对结构场和流场的有限元分析结果交叉迭代求解，即所谓流流固固耦耦合合的问题。由于有限元的应用越来越深入，人们关注的问题越来越复杂，耦合场的求解必定成为CAE软件的发展方向。6有限元方法的发展趋势4、由单一结构场求解发展到耦合场问题的6有限元方法的发展趋势有限元方法的发展趋势5、程序面向用户的开放性面向用户的开放性随着商业化的提高，各软件开发商为了扩大自己的市场份额，满足用户的需求，在软件的功功能能、易易用用性性等方面花费了大量的投资，但由于用户的要求千差万别，不管他们怎样努力也不可能满足所有用户的要求，因此必须给用户一个开放的环境，允许用户根据自己的实际情况对软件进行扩充，包括用户自定义单元特性、用户自定义材料本构（结构本构、热本构、流体本构）、用户自定义流场边界条件、用户自定义结构断裂判据和裂纹扩展规律等等。关注有限元的理论发展，采用最先进的算法技术，扩充软件的能力，提高软件性能以满足用户不断增长的需求，是CAE软件开发商的主攻目标，也是其产品持续占有市场，求得生存和发展的根本之道。6有限元方法的发展趋势5、程序面向用户的开放性