有限元分析理论基础.docx

1、有限元分析理论基础2 有限元法的基本原理2.1有限元简介有限元法：把求解区域看作由许多小的在节点处相互连接的单元（子域）所构成，其模型给出基本方程的分片（子域）近似解，由于单元（子域）可以被分割成各种形状和大小不同的尺寸，所以它能很好地适应复杂的几何形状、复杂的材料特性和复杂的边界条件有限元模型：它是真实系统理想化的数学抽象。由一些简单形状的单元组成，单元之间通过节点连接，并承受一定载荷。有限元分析：是利用数学近似的方法对真实物理系统（几何和载荷工况）进行模拟。并利用简单而又相互作用的元素，即单元，就可以用有限数量的未知量去逼近无限未知量的真实系统。线弹性有限元是以理想弹性体为研究对象的，所考

2、虑的变形建立在小变形假设的基础上。在这类问题中，材料的应力与应变呈线性关系，满足广义胡克定律；应力与应变也是线性关系，线弹性问题可归结为求解线性方程问题，所以只需要较少的计算时间。如果采用高效的代数方程组求解方法，也有助于降低有限元分析的时间。线弹性有限元一般包括线弹性静力学分析与线弹性动力学分析两方面。非线性问题与线弹性问题的区别：1）非线性问题的方程是非线性的，一般需要迭代求解；2）非线性问题不能采用叠加原理；3）非线性问题不总有一致解，有时甚至没有解。有限元求解非线性问题可分为以下三类：1）材料非线性问题材料的应力和应变是非线性的，但应力与应变却很微小，此时应变与位移呈线性关系，这类问题

3、属于材料的非线性问题。由于从理论上还不能提供能普遍接受的本构关系，所以，一般材料的应力与应变之间的非线性关系要基于试验数据，有时非线性材料特性可用数学模型进行模拟，尽管这些模型总有他们的局限性。在工程实际中较为重要的材料非线性问题有：非线性弹性（包括分段线弹性）、弹塑性、粘塑性及蠕变等。2）几何非线性问题几何非线性问题是由于位移之间存在非线性关系引起的。当物体的位移较大时，应变与位移的关系是非线性关系。研究这类问题一般都是假定材料的应力和应变呈线性关系。它包括大位移大应变及大位移小应变问题。如结构的弹性屈曲问题属于大位移小应变问题，橡胶部件形成过程为大应变问题。3）非线性边界问题在加工、密封、

4、撞击等问题中，接触和摩擦的作用不可忽视，接触边界属于高度非线性边界。平时遇到的一些接触问题，如齿轮传动、冲压成型、轧制成型、橡胶减振器、紧配合装配等，当一个结构与另一个结构或外部边界相接触时通常要考虑非线性边界条件。实际的非线性可能同时出现上述两种或三种非线性问题。2.2刚塑性有限元法及其变分原理介绍2.2.1刚塑性有限元法简介刚塑性有限元法是在1973 年提出来的，这种方法虽然也基于小应变的位移关系，但忽略了材料塑性变形时的弹性变形部分，而考虑了材料在塑性变形时的体积不变条件。它可用来计算较大变形的问题，所以近年来发展迅速，现已广泛应用于分析各种金属塑性成形过程。 刚塑性有限元法的理论基础是

5、变分原理，它认为在所有动可容的速度场中，使泛函取得驻值的速度场是真实的速度场。根据这个速度场可以计算出各点的应变和应力。 对于实际的金属成形加工过程，弹性变形部分远小于塑性变形部分( 弹性应变与塑性应变之比通常在1/1001/1000 )，因而可忽略弹性变形，将材料模型简化为刚塑性模型。采用刚塑性模型可大大简化有限元列式和求解过程。与弹塑性有限元法相比较，可采用较大的时间增量步长。在保证足够的工程精度的前提下，可提高计算效率。 由于刚塑性有限元法采用率方程表示，材料变形后的构形可通过在离散空间对速度的积分而获得，从而避开了应变与位移之间的几何非线性问题。由于忽略了弹性变形，刚塑性有限元法仅适合

6、于塑性变形区的分析，不能直接分析弹性区的变形和应力状态，也无法处理卸载和计算残余应力与变形。由于刚塑性模型假设，对一般的体积不可压缩材料，因为其静水压力与体积应变率无关，如要计算应力张量，还必须进行应力计算的处理。 从数学的角度来讲，有限元法是解微分方程的一种数值方法。它的基本思想是：在整个求解区域内要解某一微分方程很困难(即求出原函数)时，先用适当的单元将求解区域进行离散化，在单元内假定一个满足微分方程的简单函数作为解，求出单元内各点的解；然后，再考虑各单元间的相互影响，最后求出整个区域的场量。1.刚塑性有限元法的求解过程（1）离散化处理；（2）单元分析的基础上集合成总体方程组； （3）刚塑

7、性有限元法集合成的总体方程组为一非线性方程组，还须线性化处理并采用迭代方法求解。 2.刚塑性有限元法按照处理方法的不同分成如下5种： （1）流函数法；（2）拉格朗日乘子法；（3）罚函数法；（4）泊松系数v 接近0.5法；（5）材料可压缩性法。3.刚塑性材料基本假设对于大变形金属塑性成形问题，将变形体视为刚塑性体，即把变形中的某些过程理想化，便于数学上处理。此时，材料应满足下列假设：（1）不计材料的弹性变形； （2）材料的变形流动服从 Levy- Mises 流动法则；（3）材料是均质各向同性体；（4）材料满足体积不可压缩性；（5）不计体积力与惯性力；（6）加载条件(加载面)给出刚性区与塑性区的

8、界限。2.2.2变分原理1.第一变分原理刚塑性材料的第一变分原理又称为马尔柯夫(Markov)变分原理, 其为：在满足：（1）速度应变关系 (2.1) (2）体积不可压缩条件 (2.2)（3）速度边界条件(在Su上) (2.3)的一切动可容场u*ij,*ij中使泛函： (2.4)可以证明，在一切满足应变速率与速度关系和速度边界条件的ui中，使泛函式取驻值的ui为真实解。 2.完全广义变分原理在第一变分原理中，所选择的速度场必须满足(1)，(2)和(3)式，实际问题中，有些条件比较容易满足，而有些条件则不易满足。为了容易选择速度场，应用条件变分的概念，引用拉格朗日(Lagrangian) 乘子

9、，i和 ，将运动许可解所必须满足的条件引入泛函中，则得到新的泛函： (2.5) 在任意选取的ui、中，真实解使泛函取驻值，这就是刚塑性完全广义变分原理。3.第一变分原理和完全广义变分原理对比第一变分原理所选择的ui和只要求满足运动许可条件，而静力许可条件是通过变分近似满足的。据广义变分原理预选的ui和不受任何约束，所有的方程均由变分近似满足。所以，由第一变分原理计算的近似解较广义变分原理得到的解更精确。但前者在预选满足运动许可条件的速度场时比后者困难。4.不完全的广义变分原理 在选取运动许可解ui和时，可将其应满足的三个条件中的任意两个或一个事先得到满足，而将其余的一个或两个，通过拉格朗日乘子

10、引入泛函中，组成新的泛函，真实解使泛函取驻值，这就是不完全广义变分原理。在选择速度场时应变速率与速度的关系(1)式和速度边界条(3)式容易满足，而体积不可压缩条件(2)式难于满足。因此，可以把体积不可压缩条件用拉格朗日乘子入引入到泛函中，得到新泛函： (2.6)可以证明，在一切满足应变速率与速度关系和速度边界条件的ui中，使泛函(2.6)式取驻值的ui为真实解。 按照Markov变分原理求解时，面临速度场选取的困难。因而在实际求解时常采用不完全广义变分原理求解塑性变形过程。对刚塑性体按Markov变分原理确定的泛函为： (2.7)解决的问题是寻找某种方式将体积不可压缩条件式引入泛函(2.7)中

11、，构成新的泛函，使问题转变成对新泛函的无条件的驻值问题。通常采用拉格朗日乘子法、罚函数法及修正罚函数法来构造新的泛函。通过这样的方法将体积不可压缩条件引入后，便能求静水压力m， 从而解决了因忽略材料的弹性变形而带来的应力计算的困难。2.2.3刚塑性材料流动的基本方程欲求解变形体在塑性变形时的场变量，首先要建立基本方程组。设变形体的体积为v，在v内给定体力pi；表面积为S，在S的一部分力面S上给定面力qi，在St的另一部分速度（位移）面Sv上给定速度v0i，则材料在流动过程中满足下列力学基本方程1力平衡方程 (2.8)2力边界条件 即在St上 (2.9)3几何方程 (2.10)4速度边界条件即在

12、Sv上 (2.11) 5体积不可压缩方程 (2.12)6屈服准则假设材料符合密赛斯屈服准则，即: (2.13)式中k是变形过程的函数，如材料是理想刚塑性体时，k=const。得到，符合密赛斯屈服准则的应力应变关系式： （2.14）刚塑性有限元法借助于变分原理可求出近似解，对变形场的位能泛函进行变分，当变分取得驻值时，变形场满足平衡微分方程和力学边界条件。处理体积不变条件的方法有两种： 一是在假设初始速度场时，除了满足速度边界条件以外，还应严格满足体积不变条件，这种方法给假设初始速度场带来困难。另一种方法是假设初始速度场只满足速度边界条件，而对体积不变引入约束条件，即拉格朗日乘子进行有条件变分。

13、这种方法在运算中较易实现，目前已得到广泛应用。2.2.4虚功原理变形体的虚功原理可以叙述如下：变形体中满足平衡的力系在任意满足协调条件的变形状态上作的虚功等于零，即体系外力的虚功与内力的虚功之和等于零。虚功原理是虚位移（功率）原理和虚应力（率）原理的总称，它们都是与某些控制方程相等效的积分“弱”形式，虚位移（功率）原理是平衡方程和力的边界条件的等效积分“弱”形式；虚应力原理则是几何方程和位移（速度）边界条件的等效积分“弱”形式。下面来推导虚功率原理。首先考虑平衡方程 (2.15)以及力的边界条件 (2.16)我们可以采用相应的方法建立与他们等效的积分形式，在这里权函数不失一般地取速度的变分vi

14、及其边界值（取负值）。这样就可得到上面两式的等效积分形式 (2.17)对上式体积分中的第一项进行积分，并注意到应力张量是对称张量，以及由于vi是速度的变分，因而有在速度边界上vi=0，再考虑体积内部满足几何方程，则可以得到 （2.18）将上式代入(2.11)式，就得到经分部积分后的“弱”形式虚功率方程 (2.19)上式第一项是变形体内应力在虚应变率上所作之功，即内力虚功率；第二、第三项分别为体积力、面力所作的虚功率。外力和内力的虚功率和为零。这就是虚功率原理。应当指出虚功率原理是力系平衡的必要和充分条件。还应指出的是，在推导虚功效率方程时，并未涉及物理方程（应力应变率）关系，所以虚功率方程不仅

15、可以用于线弹性问题，而且还可用于非线性问题。所以虚功方程是建立塑性加工过程中有限元法的一个重要工具。2.3 UG简介UG（Unigraphics NX）是Siemens PLM Software公司出品的一个产品工程解决方案，它为用户的产品设计及加工过程提供了数字化造型和验证手段。Unigraphics NX针对用户的虚拟产品设计和工艺设计的需求，提供了经过实践验证的解决方案。UG同时也是用户指南（user guide）和普遍语法（Universal Grammer）的缩写。这是一个交互式CAD/CAM(计算机辅助设计与计算机辅助制造)系统，它功能强大，可以轻松实现各种复杂实体及造型的建构。它在诞生之初主要基于工作站，但随着PC硬件的发展和个人用户的迅速增长，在PC上的应用取得了迅猛的增长，已经成为模具行业三维设计的一个主流应用。UG的开发始于1969年，它是基于C语言开发实现的。UG NX是一个在二和三维空间无结构网格上使用自适应多重网格方法开发的一个灵活的数值求解偏微分方程的软件工具。其一个给定过程的有效模拟需要来自于应用领域(自然科学或工程)、数学(分析和数值数学)及计算机科学的知识。然而，所有这些技术在复杂应用中的使用并不是太容易。这是因为组合所有这