第11章有限元原理及其在CAD中的应用ppt课件

1、第十一章 有限元原理及其在CAD中的运用o 有限元概述o 有限元分析的原理及步骤o 有限元分析的前处置o 有限元分析的后处置11.1 有限元概述o有限元法或称有限单元法、有线元素法是求解复杂工程问题的一种近似数值值分析方法o它的根本概念是将一个外形复杂的延续体的求解区域分解成有限个单元组成的等效组合体，经过将延续体离散化，把求解延续体的场变量应力、位移、压力和温度等问题简化为求解有限个单元节点上的场变量值o此时求解的根本方程将是一个代数方程组，而不是原来的描画真实延续体场变量的微分方程组11.1 有限元概述o大大降低求解的难度，其解的精度取决于所采用的单元类型、数量以及对单元的插值函数o有限元

2、法是近30年来工程计算方法领域中的一项艰苦的成就11.1 有限元概述o它不仅广泛的运用在机械设计的力学分析中，而且也胜利运用在求解复杂的非线性工程中o有限元软件那么是有限元法及其运用的集中和完善的表达，它是CAD/CAM系统中的重要组成部分o利用有限元这一先进的技术，我们在设计阶段就可以预测产品的性能，减少许多原型制造及测试实验任务o这样既可以缩短产品设计周期节省实验费用，又可以优化产品的设计，防止了产品的大贮藏设计及缺乏设计11.2 有限元分析的原理及步骤o经过对弹性力学领域中最简单的平面静力问题的讨论，来阐明有限元方法的根本思绪和主要概念o它推导过程与空间三维等其它问题是一致的11.2.1

3、弹性力学平面问题的类型o实践上任何弹性物体都是处在三维受力形状，因此都是空间问题o但是在一定条件下，相当多的实践问题可以近似地按平面问题处置，从而使问题大大简化o平面应力问题o平面应变问题平面应力问题o如下图的薄板，板的厚度z方向的尺寸相对尺寸很小，板边上受有平行于oxy平面并沿z轴均匀分布的载荷o可以将板内各点上沿z轴方向的应力分量 近似处置为零，而各点上的三个应力分量 都平行于oxy平面，故称之为平面应力问题o在平面应力问题中，虽然应力 ，但应变xyxzx,xyyx,0z0z平面应变问题o例如如下图的圆柱体，其长度z方向比直径大很多，载荷平行于截面oxy，且沿着z轴均匀分布o这时，可以将体

4、内各点上沿着z轴方向的三个应变分量 近似地处置为零，而各点上其他的三个应变分量 都平行与oxy面，故而称为平面应变问题o在平面应变问题中，虽然 ，但应力o有些机械零件如花键轴等，可以简化为平面应变问题来处置xzyzz,xyyx,0z0z11.2.2 平面问题有限元分析的步骤o根本问题：知物体区域的边境上的约束条件所受的作用力，求解区域内各点的位移和应力o处理这个根本问题的有限元分析过程普通遵照以下步骤o离散化o单元分析o整体分析o求解未知量离散化o由无限个质点的延续体简化为有限个单元在节点出衔接而成的集合体，称为离散化o离散化的总目的是：将物体分解成充分小的单元，使得简化的位移模型可以在单元内

5、足够近似地来表示准确解，从而在整体上获得称心的计算结果o另一方面又必需留意到，单元不能分得太细，以免任务量过大，充分小而有限小，这是离散化的根本原那么离散化o有限元分析的结果是在物体的这些离散点出计算出位移和应力的近似值o有限元平面问题中常用的单元方式有三角形三节点单元，矩形四节点单元，三角形六节点单元，等参四边形八节点单元等o三角形三节点单元在平面问题中是最简单最常用的单元方式离散化要留意以下问题o当分析对象的构造具有对称性时，构造的几何外形和支撑条件，外载荷的大小和分布对称于x轴和y轴,这是只需对四分之一的构造进展分析，这样可节省近四分之三的计算任务量，假设只对称于某一轴也可计算一半o恣意

6、三角形单元的顶点，必需是相邻单元的顶点，而不能是相邻单元的内点，三角形单元各边边长不应相差太大，不要出现钝角离散化要留意以下问题o在应力集中或应力变化较大的区域，单元应分得细一点，厚度成资料发生突变的地方，应规定为单元的边境限o节点是一个很重要的概念，单元之间仅在节点处铰接，单元之间的力只经过节点传送，每个单元所受的载荷均应该按精神等效原那么移置到节点上。在位移受约束的节点上，应根据实践情况设置约束条件，当节点沿某一方向上的位移为零时，那么设置相应的连杆支座；当节点为固定点，即在oxy平面上不能挪动时，那么设一个铰链支座离散化要留意以下问题o单元划分之后，要对全部单元和全部节点进展编码，要总体

7、节点编码和部分节点o节点的总体编码是由人们自行选定的，然而编码安排得适当，得到的总体刚度矩阵回紧凑得多，可以接生计算机存储量和计算时间o离散化是一个重要的环节，单元划分的方式、大小及节点陈列对有限元分析计算的结果有很大的影响单元分析o单元分析的根本义务是：对于根本未知量 与其对应量 ，推导出节点位移 和节点力 之间的关系单位刚度矩阵 ，并建立单元的平衡方程式o单元位移方式o几何方程o弹性方程o虚功方程 e eP e eP ek eeekP单元位移方式o对整个弹性体而言，内部各点的位移情况极为复杂，不大能够用某种函数来描画o由于进展了离散化，整个弹性体分割成为细小的单元，在细小的单元内，内部各点

8、的变化情况比较简单，因此可以用单元内点坐标的函数来表示内点的独一，这种函数关系称为单元位移方式单元位移方式o由于弹性体的化整为零，是我们可以对内点位移函数关系进展简化，这是有限元素法的关键技术之一，选择位移方式时，简便的方法是将位移u,v表示为坐标x,y的函数，即采用多项式方式o而最简单的，那么是采用线性函数o选择位移方式时，该当保证有限元法解的收敛性，即当单元划分趋向于无穷小时，有限元素法的结果应趋向于问题的正确解答o可以证明，对于三角形三节点单元，上述线性函数位移方式是满足收敛条件的几何方程o由弹性力学实际可知，平面内的应变为o用矩阵可简写为oB称为应变矩阵，是外形函数经微分算子矩阵作用的

9、结果o由于采用线性位移方式，矩阵B中各个元素均为常数，因此单元内各点的应变亦为常数，所以三角形三节点单元通常称为常应变三角形单元xuyuyuxuxyyx, eB弹性方程o这里讨论单元内点应变 和单元内点应力 的转换关系，即推导出其转换矩阵弹性矩阵D，使oD称为弹性矩阵，它反映了应变与应力的内在转换关系o由于应变是常数，相应应力也当然是常数，这常数可以为是单元形心处的应力值，在相邻的各单元中，应力能够是不同的常数，因此在节点处应力有突变，随着单元减小，突变的程度也相应减小 D虚功方程o单元节点力是作用于单元顶点的外力，必然在单元内部引起相应的应力o设此单元各节点有虚位移，那么在单元内有相应的虚应

10、变o根据虚功原理，节点力在顶点的虚位移上所作的虚功应与单元内部应力在虚应变上所作的虚功相等，以此可以推导出单元内点应力与单元节点应力之间的转换关系整体分析o显然，单元分析得出的仅仅是部分的信息，各个单元靠节点衔接起来组成整体，因此必需从全局进展分析o将各个单元的方程单元刚度矩阵按照坚持节点处位移延续性的方式组合起来，就可以得到整个物体的平衡方程整体刚度矩阵，再按照给顶的位移边境条件修正这些方程，使平衡方程有解求解未知量o求解位移o求解应力、应变和约束力求解位移o求解未知位移向量 的关键在于求解平衡方程组 ，当单元个数越多时，节点未知量个数也越多，由此建立起来的平衡方程组，将是一个高阶的线性代数

11、方程组o求解这样的方程组，即使利用大型计算机来运算，也存在着存储容量大和计算速度慢的问题o针对这两个问题，可以采取以下措施：o利用整体刚度矩阵的特性，寻觅接生存储容量的措施o选择有效的计算方法和优良软件，以抑制计算速度的困难 PK 求解位移o调查整体刚度矩阵可以看出存在如下特点：o对称性：整体刚度矩阵是一个对称矩阵，利用这个特性，可以在计算机中只存储矩阵的上三角部分，从而节省近一半的存储容量o稀疏性：整体刚度矩阵是一个稀疏矩阵，其决大多数元素为零，对于大多数问题来说，非零元素占元素总数的50%，假设设法只存储非零元素，存储信息又大大减少o带状分布性：刚度矩阵的非零元素分布在以主对角线为中心的带

12、型区域内，这个矩阵称为带状矩阵，利用带状矩阵的对称性，可在计算机中只存储半带区域的元素计算应力、应变和约束力o采用位移法，一旦位移计算出以后，那么其他各量均可以计算o由于采用线性位移方式，各单元的应变与应力为常量o在这种情况下，应把算出的应力看成为单元形心处的应力计算应力、应变和约束力o假设欲求出接点上的应力，还必需运用插值计算出每一个单元的节点应力值后，再将相邻单元共用的节点上的应力值进展平均o也可以计算出主应力及按照某种强度实际计算的复合应力，以便为后处置做预备o将计算出的位移代入，即可求出各约束反力11.3 有限元分析的前处置o有限元分析的前置处置本质有限元分析置前的诸多数据处置过程o过

13、去进展有限元分析时，经常采用人工生成有限元模型的数据，既耗费时间又易出错，特别是构造复杂，计算精度要求高时，元素网格的划分时间显著添加11.3 有限元分析的前处置o手工划分有限元网格占有限元分析的任务量的45%，有限元分析结果的判别和输出占50%，而有限元分析本身只占5%o因此为减少数据输入的预备任务和提高分析、处置结果的效率，有限元系统都应配有运用方便、功能齐全的前、后置处置程序有限元分析前处置的根本内容o几何模型生成o有限元网格的自动生成o有限元属性数据的生成o模型检验、错误诊断与修正几何模型生成o这个环节主要是生成有限元计算的构造模型，有了构造模型后，可以按照许多方法生成有限元的计算数据

14、o一种途径是设计公用的有限元模型交互生成软件，实现外型与有限元前处置一体化o另一种方法是利用目前曾经有的商品化的图形软件例如AutoCAD经过的三维实体外型生成模型，在经过模型数据转换与有限元网格划分软件相衔接有限元网格的自动生成o要求生成节点坐标和编号，自动生成单元拓扑信息o可以产生多种类型的单元及其组合而成的网格，网格的疏密分布应能由用户控制o对生成的节点应能优化编号，以减少总刚度矩阵带宽，节省内存有限元属性数据的生成o将载荷、资料和边境约束条件等属性数据放入数据文件模型检验、错误诊断与修正o普通运用计算机图形显示方法来显示构造轮廓线图及网格划分图还可包括节点编号、约束与载荷情况等，这样可

15、直观的检验网格划分正确与否o此过程主要是发现并修正反复单元、反复节点、孤立单元、孤立节点、单元翘曲、单元法线不一致、单元内角与长度比异常、单元方向与自在面定义异常等有限元网格的自动生成o有限元网格的自动生成是建立有限元模型的重要技术，它在有限元前处置中占据极其重要的位置，目前常用的网格生成方法有以下几种：o以几何模型的输入数据为根底的方法o以几何模型数据构造为根底的方法以几何模型的输入数据为根底o可分为自底向上的自顶向下两中方法o自底向上法：是指先定义节点，再定义元素，最后构成有限元模型o自顶向下法：是指在已定义的构造上，根据点、线、面之间的关系，自动生成网格的方法o在实践运用中两者可以结合运

16、用以几何模型数据构造为根底o可分为分块分割法和整体分割法o分块分割法，是指先把组成复杂形体的几何模型分割成假设干适于自动生成网格的根本子形体，然后对子形体进展网格划分，最后组合成完好的有限元模型o整体分割法，是指利用四分法或八分法技术对几何形体实施自动划分，生成二维或三维有限元模型o假设几何模型用四叉树或八叉树构造表示，那么采用整体分割法更为适宜有限元分析的后处置o后处置是有限元分析过程的校核和输出结果的加工阶段o为防止打印输出大量的数据表格，以有效的方式对结果进展处置，将输出结果转化为直观明了的图形描画，经过图形的静、动态显示或绘制相应的图形，有效地检查、校核和输出分析结果后处置实现的主要内容o将节点位移、单元应力或内力等力学数据，转换为设计人员所关怀和熟习的设计参数，如应力集中区、最大应力值、最大挠度、最不利载荷组协作用下的应力等o对浩繁的数据进展编辑、挑选出关键的有用