下载文档原格式
第3章有限单元法在工程技术领域内，工程师常常运用数学和力学的知识将实际问题抽象成它们应遵循的基本方程（常微分方程或偏微分方程）和相应的边界条件。
对于大多数的工程技术问题，由于物体的几何形状和载荷作用方式是很复杂的，除了方程性质比较简单且几何边界相当规则的少数问题之外，试图按经典的弹性力学和塑性力学方法获得解析解是十分困难的，甚至是不可能的。
为了克服这种困难，有两条解决途径：一是引入简化假设，将方程和边界条件简化为能够处理的问题，从而得到它在简化状态下的解答。
这种方法只在有限的情况下可行，因为过多的简化将可能导致不正确的甚至错误的答案。
另一条解决途径就是数值解法，如有限差分法、边界元法、有限单元法和离散元法等。
对于非线性问题，有限单元法更为有效，且已经出现了许多通用程序。
有限单元法的主要优点是：①建立于严格理论基础上的可靠性。
因为用于建立有限元方程的变分原理或加权余量法在数学上已被证明是微分方程和边界条件的等效积分形式。
只要原问题的数学模型是正确的，同时用来求解有限元方程的算法是稳定、可靠的，如果单元满足收敛准则，则近似解最后收敛于原数学模型的精确解；②适应性强，应用范围广，不仅能成功地分析具有复杂边界条件、非线性、非均质材料、动力学等难题，而且还可以推广到解答数学方程中的其它边值问题，如热传导、电磁场、流体力学等问题；③适合计算机实现的高效性。
由于有限元分析的各个步骤可以表达成规范化的矩阵形式，最后导致求解方程可以统一为标准的矩阵代数问题，特别适合计算机的编程和执行。
已经出现了许多大型结构分析通用程序，如：NASTRAN、ASKA、ADINA、ANSYS、ABAQUS等，可以直接应用。
这些优点使有限单元法得到了广泛的应用和发展。
3.1有限单元法分析的基本步骤在工程或物理问题的数学模型（基本变量、基本方程、求解域和边界条件等）确定以后，有限单元法作为对其进行分析的数值计算方法的基本步骤如下：(1) 离散化一个复杂的弹性体可以看成是由无限个质点组成的连续体，它具有无限个自由度。
有限单元法有限单元法，是一种有效解决数学问题的解题方法。
其基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。
采用不同的权函数和插值函数形式，便构成不同的有限元方法。
有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学的数值模拟。
在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。
在河道数值模拟中，常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等。
根据所采用的权函数和插值函数的不同，有限元方法也分为多种计算格式。
从权函数的选择来说，有配置法、矩量法、最小二乘法和伽辽金法，从计算单元网格的形状来划分，有三角形网格、四边形网格和多边形网格，从插值函数的精度来划分，又分为线性插值函数和高次插值函数等。
不同的组合同样构成不同的有限元计算格式。
对于权函数，伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数 ;最小二乘法是令权函数等于余量本身，而内积的极小值则为对代求系数的平方误差最小;在配置法中，先在计算域内选取N个配置点。
令近似解在选定的N个配置点上严格满足微分方程，即在配置点上令方程余量为0。
插值函数一般由不同次幂的多项式组成，但也有采用三角函数或指数函数组成的乘积表示，但最常用的多项式插值函数。
有限元插值函数分为两大类，一类只要求插值多项式本身在插值点取已知值，称为拉格朗日(Lagrange)多项式插值;另一种不仅要求插值多项式本身，还要求它的导数值在插值点取已知值，称为哈密特(Hermite)多项式插值。
