有限元数值计算方法

A. 学好有限元需要哪些数学基础

高等数学（数学分析）、线性代数（高等代数）偏微分方程、常微分方程、泛函分析、复变函数等。

在数学中，有限元法（FEM，Finite Element Method）是一种为求解偏微分方程边值问题近似解的数值技术。求解时对整个问题区域进行分解，每个子区域都成为简单的部分，这种简单部分就称作有限元。它通过变分方法，使得误差函数达到最小值并产生稳定解。类比于连接多段微小直线逼近圆的思想，有限元法包含了一切可能的方法，这些方法将许多被称为有限元的小区域上的简单方程联系起来，并用其去估计更大区域上的复杂方程。它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的(较简单的）近似解，然后推导求解这个域总的满足条件(如结构的平衡条件），从而得到问题的解。这个解不是准确解，而是近似解，因为实际问题被较简单的问题所代替。由于大多数实际问题难以得到准确解，而有限元不仅计算精度高，而且能适应各种复杂形状，因而成为行之有效的工程分析手段。

B. 什么是有限元法和有限差分法

有限元法（finite element method）是一种高效能、常用的数值计算方法。科学计算领域，常常需要求解各类微分方程，而许多微分方程的解析解一般很难得到，使用有限元法将微分方程离散化后，可以编制程序，使用计算机辅助求解。

有限差分方法（finite difference method）一种求偏微分（或常微分）方程和方程组定解问题的数值解的方法，简称差分方法。

(2)有限元数值计算方法扩展阅读：

有限差分法（FDM）的起源，讨论其在静电场求解中的应用。以铝电解槽物理模型为例，采用FDM对其场域进行离散，使用MATLAB和C求解了各节点的电位。由此，绘制了整个场域的等位线和电场强度矢量分布。同时，讨论了加速收敛因子对超松弛迭代算法迭代速度的影响，以及具有正弦边界条件下的电场分布。

有限差分方法(FDM)是计算机数值模拟最早采用的方法，至今仍被广泛运用。

该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。

该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。

C. 什么是有限元法

有限元法（finite element method）是一种高效能、常用的数值计算方法。科学计算领域，常常需要求解各类微分方程，而许多微分方程的解析解一般很难得到，使用有限元法将微分方程离散化后，可以编制程序，使用计算机辅助求解。有限元法在早期是以变分原理为基础发展起来的，所以它广泛地应用于以拉普拉斯方程和泊松方程所描述的各类物理场中（这类场与泛函的极值问题有着紧密的联系）。自从1969年以来，某些学者在流体力学中应用加权余数法中的迦辽金法(Galerkin)或最小二乘法等同样获得了有限元方程，因而有限元法可应用于以任何微分方程所描述的各类物理场中，而不再要求这类物理场和泛函的极值问题有所联系。基本思想：由解给定的泊松方程化为求解泛函的极值问题。

D. 有限积分法和有限差分法

1.1 概念
有限差分方法(FDM)是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方法将求解域划分为差分网格，用有限个网格节点代替连续的求解域。有限差分法以Taylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。

1.2 差分格式
（1）从格式的精度来划分，有一阶格式、二阶格式和高阶格式。
（2）从差分的空间形式来考虑，可分为中心格式和逆风格式。
（3）考虑时间因子的影响，差分格式还可以分为显格式、隐格式、显隐交替格式等。

目前常见的差分格式，主要是上述几种形式的组合，不同的组合构成不同的差分格式。差分方法主要适用于有结构网格，网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。

1.3 构造差分的方法
构造差分的方法有多种形式，目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式：一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等，其中前两种格式为一阶计算精度，后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合，可以组合成不同的差分计算格式。

2. FEM

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

2.2 原理
有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学、土力学的数值模拟。在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。在河道数值模拟中，常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等。

根据所采用的权函数和插值函数的不同，有限元方法也分为多种计算格式。
（1）从权函数的选择来说，有配置法、矩量法、最小二乘法和伽辽金法；
（2）从计算单元网格的形状来划分，有三角形网格、四边形网格和多边形网格；
（3）从插值函数的精度来划分，又分为线性插值函数和高次插值函数等。
不同的组合同样构成不同的有限元计算格式。

对于权函数，伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数；最小二乘法是令权函数等于余量本身，而内积的极小值则为对代求系数的平方误差最小；在配置法中，先在计算域内选取N个配置点。令近似解在选定的N个配置点上严格满足微分方程，即在配置点上令方程余量为0。插值函数一般由不同次幂的多项式组成，但也有采用三角函数或指数函数组成的乘积表示，但最常用的多项式插值函数。

有限元插值函数分为两大类，一类只要求插值多项式本身在插值点取已知值，称为拉格朗日(Lagrange)多项式插值；另一种不仅要求插值多项式本身，还要求它的导数值在插值点取已知值，称为哈密特(Hermite)多项式插值。单元坐标有笛卡尔直角坐标系和无因次自然坐标，有对称和不对称等。常采用的无因次坐标是一种局部坐标系，它的定义取决于单元的几何形状，一维看作长度比，二维看作面积比，三维看作体积比。在二维有限元中，三角形单元应用的最早，近来四边形等参元的应用也越来越广。对于二维三角形和四边形电源单元，常采用的插值函数为有Lagrange插值直角坐标系中的线性插值函数及二阶或更高阶插值函数、面积坐标系中的线性插值函数、二阶或更高阶插值函数等。

2.3 基本原理与解题步骤
对于有限元方法，其基本思路和解题步骤可归纳为：
(1)建立积分方程，根据变分原理或方程余量与权函数正交化原理，建立与微分方程初边值问题等价的积分表达式，这是有限元法的出发点。
(2)区域单元剖分，根据求解区域的形状及实际问题的物理特点，将区域剖分为若干相互连接、不重叠的单元。区域单元划分是采用有限元方法的前期准备工作，这部分工作量比较大，除了给计算单元和节点进行编号和确定相互之间的关系之外，还要表示节点的位置坐标，同时还需要列出自然边界和本质边界的节点序号和相应的边界值。
(3)确定单元基函数，根据单元中节点数目及对近似解精度的要求，选择满足一定插值条件的插值函数作为单元基函数。有限元方法中的基函数是在单元中选取的，由于各单元具有规则的几何形状，在选取基函数时可遵循一定的法则。
(4)单元分析：将各个单元中的求解函数用单元基函数的线性组合表达式进行逼近；再将近似函数代入积分方程，并对单元区域进行积分，可获得含有待定系数(即单元中各节点的参数值)的代数方程组，称为单元有限元方程。
(5)总体合成：在得出单元有限元方程之后，将区域中所有单元有限元方程按一定法则进行累加，形成总体有限元方程。
(6)边界条件的处理：一般边界条件有三种形式，分为本质边界条件(狄里克雷边界条件)、自然边界条件(黎曼边界条件)、混合边界条件(柯西边界条件)。对于自然边界条件，一般在积分表达式中可自动得到满足。对于本质边界条件和混合边界条件，需按一定法则对总体有限元方程进行修正满足。
(7)解有限元方程：根据边界条件修正的总体有限元方程组，是含所有待定未知量的封闭方程组，采用适当的数值计算方法求解，可求得各节点的函数值。

3. 有限体积法

有限体积法（FiniteVolumeMethod）又称为控制体积法。其基本思路是：将计算区域划分为一系列不重复的控制体积，并使每个网格点周围有一个控制体积；将待解的微分方程对每一个控制体积积分，便得出一组离散方程。其中的未知数是网格点上的因变量的数值。为了求出控制体积的积分，必须假定值在网格点之间的变化规律，即假设值的分段的分布的分布剖面。从积分区域的选取方法看来，有限体积法属于加权剩余法中的子区域法；从未知解的近似方法看来，有限体积法属于采用局部近似的离散方法。简言之，子区域法属于有限体积发的基本方法。有限体积法的基本思路易于理解，并能得出直接的物理解释。离散方程的物理意义，就是因变量在有限大小的控制体积中的守恒原理，如同微分方程表示因变量在无限小的控制体积中的守恒原理一样。限体积法得出的离散方程，要求因变量的积分守恒对任意一组控制体积都得到满足，对整个计算区域，自然也得到满足。这是有限体积法吸引人的优点。有一些离散方法，例如有限差分法，仅当网格极其细密时，离散方程才满足积分守恒；而有限体积法即使在粗网格情况下，也显示出准确的积分守恒。就离散方法而言，有限体积法可视作有限单元法和有限差分法的中间物。有限单元法必须假定值在网格点之间的变化规律（既插值函数），并将其作为近似解。有限差分法只考虑网格点上的数值而不考虑值在网格点之间如何变化。有限体积法只寻求的结点值，这与有限差分法相类似；但有限体积法在寻求控制体积的积分时，必须假定值在网格点之间的分布，这又与有限单元法相类似。在有限体积法中，插值函数只用于计算控制体积的积分，得出离散方程之后，便可忘掉插值函数；如果需要的话，可以对微分方程中不同的项采取不同的插值函数。

4. 比较分析
有限差分法（FDM）:直观，理论成熟，精度可眩但是不规则区域处理繁琐，虽然网格生成可以使FDM应用于不规则区域，但是对区域的连续性等要求较严。使用FDM的好处在于易于编程，易于并行。
有限元方法（FEM）:适合处理复杂区域，精度可眩缺憾在于内存和计算量巨大。并行不如FDM和FVM直观。不过FEM的并行是当前和将来应用的一个不错的方向。
有限容积法:适于流体计算，可以应用于不规则网格，适于并行。但是精度基本上只能是二阶了。FVM的优势正逐渐显现出来，FVM在应力应变，高频电磁场方面的特殊的优点正在被人重视。

比较一下：
有限容积法和有限差分法：一个区别就是有限容积法的截差是不定的（跟取的相邻点有关，积分方法离散方程），而有限差分就可以直接知道截差（微分方法离散方程）。有限容积法和有限差分法最本质的区别是，前者是根据积分方程推导出来的（即对每个控制体积分），后者直接根据微分方程推导出来，所以前者的精度不但取决于积分时的精度，还取决与对导数处理的精度，一般有限容积法总体的精度为二阶，因为积分的精度限制，当然有限容积法对于守恒型方程导出的离散方程可以保持守恒型；而后者直接由微分方程导出，不涉及积分过程，各种导数的微分借助Taylor展开，直接写出离散方程，当然不一定有守恒性，精度也和有限容积法不一样，一般有限差分法可以使精度更高一些。
当然二者有联系，有时导出的形式一样，但是概念上是不一样的。

至于有限容积法和有限元相比，有限元在复杂区域的适应性对有限容积是毫无优势可言的，至于有限容积的守恒性，物理概念明显的这些特点，有限元是没有的。目前有限容积在精度方面与有限元法有些差距。

有限元方法比有限差分优越的方面主要在能适应不规则区域，但是这只是指的是传统意义上的有限差分，现在发展的一些有限差分已经能适应不规则区域。对于椭圆型方程，如果区域规则，传统有限差分和有限元都能解，在求解效率，这里主要指编程负责度和收敛快慢、内存需要，肯定有限差分有优势。

E. 有限元法的介绍

有限元法（finite element method）是一种高效能、常用的数值计算方法。科学计算领域，常常需要求解各类微分方程，而许多微分方程的解析解一般很难得到，使用有限元法将微分方程离散化后，可以编制程序，使用计算机辅助求解。有限元法在早期是以变分原理为基础发展起来的，所以它广泛地应用于以拉普拉斯方程和泊松方程所描述的各类物理场中（这类场与泛函的极值问题有着紧密的联系）。自从1969年以来，某些学者在流体力学中应用加权余数法中的迦辽金法(Galerkin)或最小二乘法等同样获得了有限元方程，因而有限元法可应用于以任何微分方程所描述的各类物理场中，而不再要求这类物理场和泛函的极值问题有所联系。基本思想：由解给定的泊松方程化为求解泛函的极值问题。

F. 数值计算方法、有限元法、无网格法的关系

有限元边界元之类的算法都是用来解带有边界条件的偏微分方程， 数值计算教材一般不会介绍这类特殊问题的算法， 一般只介绍最基本常见的算法
有限元是有网格的算法， 跟无网格的算法明显是不同的， 所谓“交叉”，既然是解同类的问题， 有交叉也有各自特点这是正常的

G. 有限元基于什么理论

有限元方法是求解数学物理问题的一种数值计算方法，起源于固体力学，然后迅速扩展到流体力学、传热学、电磁学等其他物理领域。

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

有限元模型 是真实系统理想化的数学抽象。

H. 有限元法和数值分析法有什么区别

有限元法是数值分析法中的一种，是一套求微分方程的系统化数值计算方法，是解决力学问题比较有效的数值计算方法，是将数值计算转换为矩阵计算，有利于计算机运算。数值分析法就是构造一个比较简单的函数关系，来求解方程的近似值。

I. 谁能解释下什么是有限元。

有限元

有限元法（FEA，Finite Element Analysis）的基本概念是用较简单的问题代替复杂问题后再求解。它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的(较简单的）近似解，然后推导求解这个域总的满足条件(如结构的平衡条件），从而得到问题的解。这个解不是准确解，而是近似解，因为实际问题被较简单的问题所代替。由于大多数实际问题难以得到准确解，而有限元不仅计算精度高，而且能适应各种复杂形状，因而成为行之有效的工程分析手段。

英文：Finite Element 有限单元法是随着电子计算机的发展而迅速发展起来的一种现代计算方法。它是50年代首先在连续体力学领域--飞机结构静、动态特性分析中应用的一种有效的数值分析方法，随后很快广泛的应用于求解热传导、电磁场、流体力学等连续性问题。 有限元法分析计算的思路和做法可归纳如下：
编辑本段1） 物体离散化
将某个工程结构离散为由各种单元组成的计算模型，这一步称作单元剖分。离散后单元与单元之间利用单元的节点相互连接起来；单元节点的设置、性质、数目等应视问题的性质，描述变形形态的需要和计算进度而定（一般情况单元划分越细则描述变形情况越精确，即越接近实际变形，但计算量越大）。所以有限元中分析的结构已不是原有的物体或结构物，而是同新材料的由众多单元以一定方式连接成的离散物体。这样，用有限元分析计算所获得的结果只是近似的。如果划分单元数目非常多而又合理，则所获得的结果就与实际情况相符合。
编辑本段2） 单元特性分析
A、 选择位移模式 在有限单元法中，选择节点位移作为基本未知量时称为位移法；选择节点力作为基本未知量时称为力法；取一部分节点力和一部分节点位移作为基本未知量时称为混合法。位移法易于实现计算自动化，所以，在有限单元法中位移法应用范围最广。 当采用位移法时，物体或结构物离散化之后，就可把单元总的一些物理量如位移，应变和应力等由节点位移来表示。这时可以对单元中位移的分布采用一些能逼近原函数的近似函数予以描述。通常，有限元法我们就将位移表示为坐标变量的简单函数。这种函数称为位移模式或位移函数。 B、 分析单元的力学性质 根据单元的材料性质、形状、尺寸、节点数目、位置及其含义等，找出单元节点力和节点位移的关系式，这是单元分析中的关键一步。此时需要应用弹性力学中的几何方程和物理方程来建立力和位移的方程式，从而导出单元刚度矩阵，这是有限元法的基本步骤之一。 C、 计算等效节点力 物体离散化后，假定力是通过节点从一个单元传递到另一个单元。但是，对于实际的连续体，力是从单元的公共边传递到另一个单元中去的。因而，这种作用在单元边界上的表面力、体积力和集中力都需要等效的移到节点上去，也就是用等效的节点力来代替所有作用在单元上的力。
编辑本段3） 单元组集
利用结构力的平衡条件和边界条件把各个单元按原来的结构重新连接起来，形成整体的有限元方程 （1-1） 式中，K是整体结构的刚度矩阵；q是节点位移列阵；f是载荷列阵。
编辑本段4） 求解未知节点位移
解有限元方程式（1-1）得出位移。这里，可以根据方程组的具体特点来选择合适的计算方法。 通过上述分析，可以看出，有限单元法的基本思想是"一分一合"，分是为了就进行单元分析，合则为了对整体结构进行综合分析。 有限元的发展概况 1943年 courant在论文中取定义在三角形域上分片连续函数，利用最小势能原理研究St.Venant的扭转问题。 1960年 clough的平面弹性论文中用“有限元法”这个名称。 1965年 冯康发表了论文“基于变分原理的差分格式”，这篇论文是国际学术界承认我国独立发展有限元方法的主要依据。 1970年 随着计算机和软件的发展，有限元发展起来。 涉及的内容：有限元所依据的理论，单元的划分原则，形状函数的选取及协调性。 有限元法涉及：数值计算方法及其误差、收敛性和稳定性。 应用范围：固体力学、流体力学、热传导、电磁学、声学、生物力学 求解的情况：杆、梁、板、壳、块体等各类单元构成的弹性（线性和非线性）、弹塑性或塑性问题（包括静力和动力问题）。能求解各类场分布问题（流体场、温度场、电磁场等的稳态和瞬态问题），水流管路、电路、润滑、噪声以及固体、流体、温度相互作用的问题。
编辑本段5）有限元的未来是多物理场耦合
5）有限元的未来是多物理场耦合 随着计算机技术的迅速发展，在工程领域中，有限元分析（FEA）越来越多地用于仿真模拟，来求解真实的工程问题。这些年来，越来越多的工程师、应用数学家和物理学家已经证明这种采用求解偏微分方程（PDE）的方法可以求解许多物理现象，这些偏微分方程可以用来描述流动、电磁场以及结构力学等等。有限元方法用来将这些众所周知的数学方程转化为近似的数字式图象。 早期的有限元主要关注于某个专业领域，比如应力或疲劳，但是，一般来说，物理现象都不是单独存在的。例如，只要运动就会产生热，而热反过来又影响一些材料属性，如电导率、化学反应速率、流体的粘性等等。这种物理系统的耦合就是我们所说的多物理场，分析起来比我们单独去分析一个物理场要复杂得多。很明显，我们现在需要一个多物理场分析工具。 在上个世纪90年代以前，由于计算机资源的缺乏，多物理场模拟仅仅停留在理论阶段，有限元建模也局限于对单个物理场的模拟，最常见的也就是对力学、传热、流体以及电磁场的模拟。看起来有限元仿真的命运好像也就是对单个物理场的模拟。 现在这种情况已经开始改变。经过数十年的努力，计算科学的发展为我们提供了更灵巧简洁而又快速的算法，更强劲的硬件配置，使得对多物理场的有限元模拟成为可能。新兴的有限元方法为多物理场分析提供了一个新的机遇，满足了工程师对真实物理系统的求解需要。有限元的未来在于多物理场求解。 千言万语道不尽，下面只能通过几个例子来展示多物理场的有限元分析在未来的一些潜在应用。 压电扩音器（Piezoacoustic transcer）可以将电流转换为声学压力场，或者反过来，将声场转换为电流场。这种装置一般用在空气或者液体中的声源装置上，比如相控阵麦克风，超声生物成像仪，声纳传感器，声学生物治疗仪等，也可用在一些机械装置比如喷墨机和压电马达等。 压电扩音器涉及到三个不同的物理场：结构场，电场以及流体中的声场。只有具有多物理场分析能力的软件才能求解这个模型。 压电材料选用PZT5-H晶体，这种材料在压电传感器中用得比较广泛。在空气和晶体的交界面处，将声场边界条件设置为压力等于结构场的法向加速度，这样可以将压力传到空气中去。另外，晶体域中又会因为空气压力对其的影响而产生变形。仿真研究了在施加一个幅值200V，震荡频率为300 KHz的电流后，晶体产生的声波传播。这个模型的描述及其完美的结果表明在任何复杂的模型下，我们都可以用一系列的数学模型进行表达，进而求解。 多物理场建模的另外一个优势就是在学校里，学生们直观地获取了以前无法见到的一些现象，而简单易懂的表达方式也获得了学生们的好感。这只是Krishan Kumar Bhatia博士在纽约Glassboro的Rowan 大学给高年级的毕业生讲授传热方程课程时介绍建模及分析工具所感受到的，他的学生的课题是如何冷却一个摩托车的发动机箱。Bhatia博士教他们如何利用“设计－制造－检测”的理念来判断问题、找出问题、解决问题。如果没有计算机仿真的应用，这种方法在课堂上推广是不可想象的，因为所需费用实在是太大了。 COMSOL Multiphysics拥有优秀的用户界面，可以使学生方便地设置传热问题，并很快得到所需要的结果。“我的目标是使每个学生都能了解偏微分方程，当下次再遇到这样的问题时，他们不会再担心，” Bhatia博士说，“这不需要了解太多的分析工具，总的来说，学生都反映‘这个建模工具太棒了’”。 很多优秀的高科技工程公司已经看到多物理场建模可以帮助他们保持竞争力。多物理场建模工具可以让工程师进行更多的虚拟分析而不是每次都需要进行实物测试。这样，他们就可以快速而经济地优化产品。在印度尼西亚的Medrad Innovations Group中，由John Kalafut博士带领着一个研究小组，采用多物理场分析工具来研究细长的注射器中血细胞的注射过程，这是一种非牛顿流体，而且具有很高的剪切速率。 通过这项研究，Medrad的工程师制造了一个新颖的装置称为先锋型血管造影导管（Vanguard Dx Angiographic Catheter）。同采用尖喷嘴的传统导管相比，采用扩散型喷嘴的新导管使得造影剂分布得更加均匀。造影剂就是在进行X光拍照时，将病变的器官显示得更加清楚的特殊材料。 另外一个问题就是传统导管在使用过程中可能会使得造影剂产生很大的速度，进而可能会损伤血管。先锋型血管造影导管降低了造影剂对血管产生的冲击力，将血管损伤的可能性降至最低。 关键的问题就是如何去设计导管的喷嘴形状，使其既能优化流体速度又能减少结构变形。Kalafut的研究小组利用多物理场建模方法将层流产生的力耦合到应力应变分 析中去，进而对各种不同喷嘴的形状、布局进行流固耦合分析。“我们的一个实习生针对不同的流体区域建立不同的喷嘴布局，并进行了分析，” Kalafut博士说，“我们利用这些分析结果来评估这些新想法的可行性，进而降低实体模型制造次数”。 摩擦搅拌焊接（FSW），自从1991年被申请专利以来，已经广泛应用于铝合金的焊接。航空工业最先开始采用这些技术，现在正在研究如何利用它来降低制造成本。在摩擦搅拌焊接的过程中，一个圆柱状具有轴肩和搅拌头的刀具旋转插入两片金属的连接处。旋转的轴肩和搅拌头用来生热，但是这个热还不足以融化金属。反之，软化呈塑性的金属会形成一道坚实的屏障，会阻止氧气氧化金属和气泡的形成。粉碎，搅拌和挤压的动作可以使焊缝处的结构比原先的金属结构还要好，强度甚至可以到原来的两倍。这种焊接装置甚至可以用于不同类型的铝合金焊接。 空中客车（AirBus）资助了很多关于摩擦搅拌焊接的研究。在制造商大规模投资和重组生产线之前，Cranfield大学的Paul Colegrove博士利用多物理场分析工具帮助他们理解了加工过程。 第一个研究成果是一个摩擦搅拌焊接的数学模型，这让空客的工程师“透视”到焊缝中来检查温度分布和微结构的变化。Colegrove博士和他的研究小组还编写了一个带有图形界面的仿真工具，这样空客的工程师可以直接提取材料的热力属性以及焊缝极限强度。 在这个摩擦搅拌焊接的模拟过程中，将三维的传热分析和二维轴对称的涡流模拟耦合起来。传热分析计算在刀具表面施加热流密度后，结构的热分布。可以提取出刀具的位移，热边界条件，以及焊接处材料的热学属性。接下来将刀具表面处的三维热分布映射到二维模型上。耦合起来的模型就可以计算在加工过程中热和流体之间的相互作用。 将基片的电磁、电阻以及传热行为耦合起来需要一个真正的多物理场分析工具。一个典型的应用是在半导体的加工和退火的工艺中，有一种利用感应加热的热壁熔炉，它用来让半导体晶圆生长，这是电子行业中的一项关键技术。 例如，金刚砂在2,000°C的高温环境下可以取代石墨接收器，接收器由功率接近10KW的射频装置加热。在如此高温下要保持炉内温度的均匀，炉腔的设计至关重要。经过多物理场分析工具的分析，发现热量主要是通过辐射的方式进行传播的。在模型内不仅可以看到晶圆表面温度的分布，还可以看到熔炉的石英管上的温度分布。 在电路设计中，影响材料选择的重要方面是材料的耐久性和使用寿命。电器小型化的趋势使得可在电路板上安装的电子元件发展迅猛。众所周知，安装在电路板上的电阻以及其他一些元件会产生大量的热，进而可能使得元件的焊脚处产生裂缝，最后导致整个电路板报废。 多物理场分析工具可以分析出整个电路板上热量的转移，结构的应力变化以及由于温度的上升导致的变形。这样做可以用来提升电路板设计的合理性以及材料选择的合理性。 计算机能力的提升使得有限元分析由单场分析到多场分析变成现实，未来的几年内，多物理场分析工具将会给学术界和工程界带来震惊。单调的“设计－校验”的设计方法将会慢慢被淘汰，虚拟造型技术将让你的思想走得更远，通过模拟仿真将会点燃创新的火花。

J. 有限元分析方法是指什么

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

有限元分析是用较简单的问题代替复杂问题后再求解。它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的（较简单的）近似解，然后推导求解这个域总的满足条件（如结构的平衡条件），从而得到问题的解。

因为实际问题被较简单的问题所代替，所以这个解不是准确解，而是近似解。由于大多数实际问题难以得到准确解，而有限元不仅计算精度高，而且能适应各种复杂形状，因而成为行之有效的工程分析手段。

(10)有限元数值计算方法扩展阅读：

有限元方法与其他求解边值问题近似方法的根本区别在于它的近似性仅限于相对小的子域中。20世纪60年代初首次提出结构力学计算有限元概念的克拉夫（Clough）教授形象地将其描绘为：“有限元法=Rayleigh Ritz法+分片函数”，即有限元法是Rayleigh Ritz法的一种局部化情况。

不同于求解（往往是困难的）满足整个定义域边界条件的允许函数的Rayleigh Ritz法，有限元法将函数定义在简单几何形状（如二维问题中的三角形或任意四边形）的单元域上（分片函数），且不考虑整个定义域的复杂边界条件，这是有限元法优于其他近似方法的原因之一。