【精选】有限元分析理论基础

有限元分析概念

有限元法：把求解区域看作由许多小的在节点处相互连接的单元（子域）所构成，其模型给出基本方程的分片（子域）近似解，由于单元（子域）可以被分割成各种形状和大小不同的尺寸，所以它能很好地适应复杂的几何形状、复杂的材料特性和复杂的边界条件

有限元模型：它是真实系统理想化的数学抽象。由一些简单形状的单元组成，单元之间通过节点连接，并承受一定载荷。

有限元分析：是利用数学近似的方法对真实物理系统（几何和载荷工况）进行模拟。并利用简单而又相互作用的元素，即单元，就可以用有限数量的未知量去逼近无限未知量的真实系统。

线弹性有限元是以理想弹性体为研究对象的，所考虑的变形建立在小变形假设的基础上。在这类问题中，材料的应力与应变呈线性关系，满足广义胡克定律；应力与应变也是线性关系，线弹性问题可归结为求解线性方程问题，所以只需要较少的计算时间。如果采用高效的代数方程组求解方法，也有助于降低有限元分析的时间。

线弹性有限元一般包括线弹性静力学分析与线弹性动力学分析两方面。

非线性问题与线弹性问题的区别：

1）非线性问题的方程是非线性的，一般需要迭代求解；

2）非线性问题不能采用叠加原理；

3）非线性问题不总有一致解，有时甚至没有解。

有限元求解非线性问题可分为以下三类：

1）材料非线性问题

材料的应力和应变是非线性的，但应力与应变却很微小，此时应变与位移呈线性关系，这类问题属于材料的非线性问题。由于从理论上还不能提供能普遍接受的本构关系，所以，一般材料的应力与应变之间的非线性关系要基于试验数据，有时非线性材料特性可用数学模型进行模拟，尽管这些模型总有他们的局限性。在工程实际中较为重要的材料非线性问题有：非线性弹性（包括分段线弹性）、弹塑性、粘塑性及蠕变等。

2）几何非线性问题

几何非线性问题是由于位移之间存在非线性关系引起的。

当物体的位移较大时，应变与位移的关系是非线性关系。研究这类问题一般都是假定材料的应力和应变呈线性关系。它包括大位移大应变及大位移小应变问题。如结构的弹性屈曲问题属于大位移小应变问题，橡胶部件形成过程为大应变问题。

3）非线性边界问题

在加工、密封、撞击等问题中，接触和摩擦的作用不可忽视，接触边界属于高度非线性边界。

平时遇到的一些接触问题，如齿轮传动、冲压成型、轧制成型、橡胶减振器、紧配合装配等，当一个结构与另一个结构或外部边界相接触时通常要考虑非线性边界条件。

实际的非线性可能同时出现上述两种或三种非线性问题。

有限元理论基础

有限元方法的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。采用不同的权函数和插值函数形式，便构成不同的有限元方法。

1.加权余量法：

是指采用使余量的加权函数为零求得微分方程近似解的方法称为加权余量法。（Weighted residual method WRM ）是一种直接从所需求解的微分方程及边界条件出发，寻求边值问题近似解的数学方法。加权余量法是求解微分方程近似解的一种有效的方法。 设问题的控制微分方程为:

在V 域内 在S 边界上 式中 :

L 、B ——分别为微分方程和边界条件中的微分算子； f 、g ——为与未知函数u 无关的已知函数域值； u ——为问题待求的未知函数 ()0

L u f -=(5.1.1)()0

B u g -=(5.1.2)

混合法对于试函数的选取最方便，但在相同精度条件下，工作量最大。对内部法和边界法必须使基函数事先满足一定条件，这对复杂结构分析往往有一定困难，但试函数一经建立，其工作量较小。

无论采用何种方法，在建立试函数时均应注意以下几点：

(1)试函数应由完备函数集的子集构成。已被采用过的试函数有幂级数、三角级数、样条函数、贝赛尔函数、切比雪夫和勒让德多项式等等。

(2)试函数应具有直到比消除余量的加权积分表达式中最高阶导数低一阶的导数连续性。

(3)试函数应与问题的解析解或问题的特解相关联。若计算问题具有对称性，应充分利用它。

显然，任何独立的完全函数集都可以作为权函数。按照对权函数的不同选择得到不同的加权余量计算方法，主要有：配点法、子域法、最小二乘法、力矩法和伽辽金法。其中伽辽金法的精度最高。

2、虚功原理

——平衡方程和几何方程的等效积分“弱”形式

虚功原理包含虚位移原理和虚应力原理，是虚位移原理和虚应力原理的总称。他们都可以认为是与某些控制方程相等效的积分“弱”形式。虚功原理：变形体中任意满足平衡的力系在任意满足协调条件的变形状态上作的虚功等于零，即体系外力的虚功与内力的虚功之和等于零。

虚位移原理是平衡方程和力的边界条件的等效积分的“弱”形式；

虚应力原理是几何方程和位移边界条件的等效积分“弱”形式。

虚位移原理的力学意义：如果力系是平衡的，则它们在虚位移和虚应变上所作的功的总和为零。反之，如果力系在虚位移（及虚应变）上所作的功的和等于零，则它们一定满足平衡方程。所以，虚位移原理表述了力系平衡的必要而充分条件。一般而言，虚位移原理不仅可以适用于线弹性问题，而且可以用于非线性弹性及弹塑性等非线性问题。

虚应力原理的力学意义：如果位移是协调的，则虚应力和虚边界约束反力在他们上面所作的功的总和为零。反之，如果上述虚力系在他们上面所作的功的和为零，则它们一定是满足协调的。所以，虚应力原理表述了位移协调的必要而充分条件。

虚应力原理可以应用于线弹性以及非线性弹性等不同的力学问题。但是必须指出，无论是虚位移原理还是虚应力原理，他们所依赖的几何方程和平衡方程都是基于小变形理论的，他们不能直接应用于基于大变形理论的力学问题。

3、最小总势能法

应变能：作用在物体上的外载荷会引起物体变形，变形期间外力所做的功以弹性能的形式储存在物体中，即为应变能。

由n 个单元和m 个节点组成的物体的总势能为总应变能和外力所做功的差：

()11

=n m e i i e i Fu ==∏Λ

-∑∑ 最小势能原理：对于一个稳定的系统，相对于平衡位置发生的位移总会使系统的总势能最小，即：

()110n m e i i e i i i i Fu u u u ==?∏??=Λ-=???∑∑，i=1,2,3,……,n

有限元法的收敛性

有限元法是一种数值分析方法，因此应考虑收敛性问题。

有限元法的收敛性是指：当网格逐渐加密时，有限元解答的序列收敛到精确解；或者当单元尺寸固定时，每个单元的自由度数越多，有限元的解答就越趋近于精确解。

有限元的收敛条件包括如下四个方面：

1）单元内，位移函数必须连续。多项式是单值连续函数，因此选择多项式作为位移函数，在单元内的连续性能够保证。

2）在单元内，位移函数必须包括常应变项。每个单元的应变状态总可以分解为不依赖于单元内各点位置的常应变和由各点位置决定的

变量应变。当单元的尺寸足够小时，单元中各点的应变趋于相等，单元的变 …… 此处隐藏：2470字，全部文档内容请下载后查看。喜欢就下载吧 ……