有限元分析教案

有限元分析教案1第一章有限元法概述在机械设计中，人们常常运用材料力学、结构力学等理论知识分析机械零构件的强度、刚度和稳定性问题。但对一些复杂的零构件，这种分析常常就必须对其受力状态和边界条件进行简化。否则力学分析将无法进行。但这种简化的处理常常导致计算结果与实际相差甚远，有时甚至失去了分析的意义。所以过去设计经验和类比占有较大比重。因为这个原因，人们也常常在设计中选择较大的安全系数。如此也就造成所设计的机械结构整体尺寸和重量偏大，而局部薄弱环节强度和刚度又不足的设计缺陷。近年来，数值计算机在工程分析上的成功运用，产生了一门全新、高效的工程计算分析学科——有限元分析方法。该方法彻底改变了传统工程分析中的做法。使计算精度和计算领域大大改善。§1.1有限元方法的发展历史、现状和将来一，历史有限元法的起源应追溯到上世纪40年代（20世纪40年代）。1943年R.Courant从数学的角度提出了有限元法的基本观点。50年代中期在对飞机结构的分析中，诞生了结构分析的矩阵方法。1960年R.W.Clough在分析弹性力学平面问题时引入了“FiniteElementMethod”这一术语，从而标志着有限元法的思想在力学分析中的广泛推广。60、70年代计算机技术的发展，极大地促进了有限元法的发展。具体表现在：1）由弹性力学的平面问题扩展到空间、板壳问题。2）由静力平衡问题——稳定性和动力学分析问题。3）由弹性问题——弹塑性、粘弹性等问题。二，现状现在有限元分析法的应用领域已经由开始时的固体力学，扩展到流体力学、传热学和电磁力学等多个传统的领域。已经形成了一种非常成熟的数值分析计算方法。大型的商业化有限元分析软件也是层出不穷，如：SAP系列的代表SAP2000（StructureAnalysisProgram）美国安世软件公司的ANSYS大型综合有限元分析软件美国航天航空局的NASTRAN系列软件除此以外，还有MASTER、ALGO、ABIQUES、ADINA、COSMOS等。三，将来有限元的发展方向最终将和CAD的发展相结合。运用“四个化”可以概括其今后的发展趋势。那就是：可视化、集成化、自动化和网络化。§1.2有限元法的特点机械零构件的受力分析方法总体说来分为解析法和数值法两大类。如大家学过的材料力学、结构力学等就是经典的解析力学分析方法。在这些解析力学方法中，弹性力学的分析方法在数学理论上是最为严谨的一种分析方法。其解题思路是：从静力、几何和物理三个方面综合考虑，建立描述弹性体的平衡、应力、应变和位移三者之间的微分方程，然后考虑边界条件，从而求出微分方程的解析解。其最大的有点就是，严密精确。缺点就是微分方程的求解困难，很多情况下，无法求解。数值方法是一种近似的计算方法。具体又分为“有限差分法”和“有限元法”。“有限差分法”是将得到的微分方程离散成近似的差分方程。通过对一系列离散的差分方程求解，得到最终的力学问题近似解。其优点就是：计算简单收敛性好。缺点是：计算程序无法标准化，在不能获得整个问题的微分方程时，该方法不能运用。由于其是将微分方程转为差分方程，所以它是一种数学近似。有限元分析教案2“有限元法”的基本思想就是“先分后合”或者“化整为零，又积零为整”。与有限差分不同，它是在力学模型上进行近似处理，也就是（分块近似）。具体做法：把连续体模型转为由有限个单元组成的离散体模型，离散体模型之间通过一些节点联系。对于每一个离散体个体选择简单的函数近似表示其中的物理变化规律（如位移等），运用力学方程推导单元的平衡方程组，然后集合所有的方程组形成表征整体结构的方程组，引入边界条件，求取最后问题的解。优点：概念清晰、易于学习理解，适用性强，便于电算化。缺点：计算精度受单元划分的影响较大。§1.3有限元分析的一般过程为了能够了解有限元分析的全貌，我们就一个简单的例子，来分析一下有限元分析的三个过程：结构离散化、单元分析、整体分析。一，结构离散化在该阶段中，要完成把连续结构的力学模型转变为离散的力学模型。处理的好坏，直接影响到最后分析结果的正确与否、计算的精度和计算的效率。根据模型的传力特性和分析的目标，正确选择单元类型。通常单元分为：一维单元、二维单元和三维单元。所谓一维单元就是指所求物理量仅随一个坐标变量而变化的单元。如桁架、平面刚架和空间刚架单元。一维单元：杆单元、梁单元。二维单元：三角形单元、四边形单元（平面类问题）三维单元：四面体单元、六面体单元等（空间问题）计算精度和计算效率：取决于单元划分的形状、大小和分布状况。通常单元愈多、愈密集，计算精度愈高，但计算效率愈低。有限元分析工作就是要在精度和效率两者之间做到有机的统一。二，单元分析进行单元分析的目的是为了到处表征单元力学特性的“单元刚度矩阵”。一般说来该过程有三种方法：1，直接法。2，虚功原理法（变分法）。3，加权余数法。直接法概念浅显，易于理解物理含义。变分法需要泛函的数学知识，其推导过程具有严谨的数学概念。加权余数法适用于泛函不存在的应用范围。本教材将运用虚功原理方程结合弹性力学和材料力学中的知识来推求几种常见单元的单刚计算公式。现在先看一个简单的阶梯轴的轴向拉伸问题例：如图所示的变截面直杆，受拉力P，运用有限元方法分析其变形。取任意单元，长度为l，面积为A，单元内任意一点的轴向位移和其位置坐标成正比，即有限元分析教案3u=a0+a1x其中a0、a1为待定系数。由于杆的两个端点节点1、2是单元上的点，所以它们应该满足上述方程。节点1，x1=0，∴u1=a0+a1×0=a0节点2，x2=l，∴u2=a0+a1×la1=(u2-u1)/l将求出的结果带入方程并整理，就得：eNuuNNulxulxxluluuu2121211211式中：N1、N2是形函数[N]形函数矩阵{δ}e节点位移向量由位移与应变的关系知道：dxdudxuduu将上面推出的位移表达式代入，可得：eeeeBllllxlxdxdNdxddxdudxuduu11上式中的[B]称为应变矩阵或几何矩阵。运用材料力学中的虎克定律，可以将应变和应力联系起来。单向应力状态的虎克定律为eeSuulElEBEE21××其中[S]称为应力矩阵。利用虚功方程可以建立力与位移之间的关系，也就是单刚方程。在后面我们将会推导出它的一般形式如下：eeeKF式中：{F}e为单元节点力向量，对我们这个例子应为[U1U2]T。[K]e为单元刚度矩阵。后面将推导出它的计算公式为dVBDBKvTe[D]矩阵是弹性矩阵。对于一维单元来说，就是E。所以我们这儿讨论的例题：111111110lEAAdxllEllKle求得单刚矩阵，也就完成了单元分析。总结单刚矩阵推导的步骤，应该分为四步：1）假定单元内位移变化的近似规律，即选择位移模式。2）运用几何关系，推求位移与应变的关系。有限元分析教案43）应用物理规律，把应变与应力联系起来。4）运用虚功方程的力与位移关系，求出单刚矩阵。单元分析是整个有限元分析的核心。不同的单元因为其力学特性不同，而具有不同的单元刚度矩阵，我们这本教材就是要学习几种常用单元矩阵的推导和计算。了解各种单元的力学特性，为以后选择单元类型打好基础。三，整体分析1，由各单元刚度矩阵组集成整个结构的总刚度矩阵。1111111lEAK1111222lEAK整个结构有三个节点，首先将单元刚度矩阵扩充为3X3的矩阵，移动各元素使之与单刚矩阵中的元素位置相对应，如下：000011011111lEAK110110000221lEAK然后直接相加。222222221111110110lEAlEAlEAlEAlEAlEAlEAlEAK2，把各单元的节点力向量组集成总的节点载荷向量。00PR3，根据边界条件，修改总刚度矩阵，获得总刚方程组。边界条件修改之前的总刚方程：32122222222111111011000uuulEAlEAlEAlEAlEAlEAlEAlEAP修改以后（采用置“0，1”法）有限元分析教案532122111111100011000uuulEAlEAlEAlEAlEAP4，求解方程组，得出总的节点位移向量。解得的解是：0222211321EAPlEAPlEAPluuu有了节点位移，再回代到前面单元推导过程中的公式×和××，就可以求得每个单元的应变和应力了。从这个简单的例子，我们了解了有限元法求解力学问题三大步骤中的内容，想必很多同学会说，这样复杂，如果运用材料力学的知识，我还来得快些。但是大家不要忘记，有限元的计算很多都是编程完成，而且现在很多的商业软件都已经完成了很多的工作。我们学习有限元主要是了解它的原理，并对常见单元的力学特性有所了解，这样对于以后运用有限元起到帮助作用。所以下面章节的内容，就是围绕这个主题展开。要达到这个目的，我们还必须学习必要的弹性力学知识。对弹性力学知识的学习，也对我们以后把握问题的本质有帮助。有限元分析教案6第二章平面问题平面问题在力学研究的课题中属弹性力学的范畴。该类问题不仅本身具有典型性，而且在机械零构件的分析中，也是应用得非常广泛。所以这类问题也称之为经典的力学问题。我们知道，实际的机械零构件都是具有三维空间尺寸的物体，理应作为三维对象处理，但是当物体的几何形状和受力状态处于某些特定的情况下，近似地简化为平面问题，不仅可以大大简化计算的工作量，而且其精度也完全能够满足所要求。如：直齿圆柱齿轮可在垂直与孔轴线的截平面内作平面应力分析就足以了解整个齿轮的受力状态；大坝的横断面可作平面应变分析来了解整个大坝受力情况等。本章是全书的重点，在这里不仅介绍弹性力学的基本知识。还将系统地讲解有限元的基本概念、原理和方法。是学习以后各章节的基础。§2.1外力、应力、应变和位移外力、应力、应变和位移的概念在材料力学中已经学习过，由于这些概念在弹性力学、有限元法中具有和在材料力学中不同的规定，弹性力学中的规定和有限元法是完全相同的，所以在这里我们将按照弹性力学的习惯表达方法把他们集中的加以阐述。一，外力外界作用在物体上的作用力，可以分为两大类：1）体积力分布在物体体积内的力。如：重力、惯性力和磁性力等。单位体积的体力在坐标轴上的分量X、Y、Z，称为体力分量。符号规定为沿坐标轴正向的为正，沿负向的为负。2）面力作用在物体表面的力。如轮压、水压等。它又可细分为集中力与分布力。面力在坐标轴上的投影，表示为X、Y、Z。符号沿正轴为正，负轴为负。二，应力弹性体受到外力作用后，内部产生的抵抗变形的内力。以弹性体中P点为定点的微单元体来考察。所谓微单元体，就是图中PA、PB、PC的边长分别为dx、dy和dz。以下简称这样的微单元体为微元体。微元体每个面上的应力都可以分解为三个应力分量。以图中红面为例，分别是σx、τxy、τxz。应力命名的规则：以应力所在面垂直的坐标轴为第一个下标，应力指向为第二下标。如果下标相同就用一个下标表示。符号规定：正面上的应力与坐标轴同向为正，反之为负。负面上的应力与坐标轴反向为正。反之为负。所谓正面就是面的外法向与坐标轴同向为正。反之为负面。作用在两个相互垂直的面上，并且垂直与该两面交线的剪应力互等。即：τxy=τyx；τyz=τzy；τxz=τzx如此以来，代表P点应力状态的应力分量应有6个，它们