《有限元法-1-泛函与变分.ppt》由会员分享，可在线阅读，更多相关《有限元法-1-泛函与变分.ppt（47页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第五章第五章 有限元法有限元法内容：内容：基于变分原理，介绍有限元法。基于变分原理，介绍有限元法。以线性静态场中一阶有限元的应用为重点，以线性静态场中一阶有限元的应用为重点，引伸到非线性场、时谐场中的分析，以及等参数有引伸到非线性场、时谐场中的分析，以及等参数有限元法的应用。限元法的应用。1特点特点p 能处理复杂区域和复杂边界条件的求解问题。能处理复杂区域和复杂边界条件的求解问题。p是是求解微分方程求解微分方程的系统化数值计算方法。的系统化数值计算方法。p比传统解法具有理论完整可靠比传统解法具有理论完整可靠.p物理意义直观明确物理意义直观明确.p解题效能强。解题效能强。2历史历史p历史历史n1
2、943年年Courant提出有限元思想。提出有限元思想。n20世纪世纪50年代初期，有限元法在复杂的年代初期，有限元法在复杂的航空结构分析航空结构分析中最先得到应用，中最先得到应用，n1960年年Clough(克拉夫克拉夫)在其著作中首先提出在其著作中首先提出有限元法有限元法（finite element method，简称，简称FEM）这个名称。这个名称。3应用应用p以变分原理为基础的有限元法，因其理论依据的以变分原理为基础的有限元法，因其理论依据的普遍性，广泛地被应用于普遍性，广泛地被应用于各种工程领域各种工程领域：n热传导、渗流、热传导、渗流、n流体力学、空气动力学、土壤力学、流体力学、
3、空气动力学、土壤力学、n机械零件强度分析、机械零件强度分析、n电磁场工程问题等。电磁场工程问题等。4p电气工程领域的应用电气工程领域的应用n1965年年Winslow首先将有限元法应用于首先将有限元法应用于电气工程电气工程问题，问题，n1969年年Silvester将有限元法推广应用于将有限元法推广应用于时谐电磁场时谐电磁场问题。问题。p至今至今n有限元法已经成为各类电磁场、电磁波工程问题定量有限元法已经成为各类电磁场、电磁波工程问题定量分析与优化设计的分析与优化设计的主导数值计算方法主导数值计算方法，n是构成各种先进、实用计算是构成各种先进、实用计算软件包的基础软件包的基础。55.1 概述概
4、述p基本思想：基本思想：传统的有限元法传统的有限元法以变分原理为基础以变分原理为基础。n首先把所要求解的微分方程数学模型首先把所要求解的微分方程数学模型边值问题，边值问题，转化为相应的转化为相应的变分问题变分问题，即，即泛函求极值问题泛函求极值问题；n然后利用然后利用剖分插值剖分插值，离散化离散化变分问题为普通变分问题为普通多元函数多元函数的极值问题的极值问题，即最终归结为一组多元的代数方程组；，即最终归结为一组多元的代数方程组；n解之即得待求边值问题的数值解。解之即得待求边值问题的数值解。6p有限元法的核心在于：有限元法的核心在于：剖分插值剖分插值。n将连续场分割为有限个单元，然用比较简单的
5、将连续场分割为有限个单元，然用比较简单的插值函插值函数数来表示每个来表示每个单元的解单元的解，n但是，它并不要求每个单元的试探解都满足边界条件，但是，它并不要求每个单元的试探解都满足边界条件，而是而是在全部单元总体合成后再引入边界条件在全部单元总体合成后再引入边界条件。n这样，就有可能这样，就有可能对于内部和边界上的单元采用同样的对于内部和边界上的单元采用同样的插值函数插值函数，使方法构造极大地得到简化。，使方法构造极大地得到简化。7n此外，由于变分原理的应用，此外，由于变分原理的应用，使使第二第二、三类三类及不同及不同媒媒质分界面质分界面上的边界条件作为上的边界条件作为自然边界条件自然边界条
6、件在总体合成在总体合成时将隐含地得到满足时将隐含地得到满足。p即即自然边界条件将被包含在泛函达到极值的要求之中自然边界条件将被包含在泛函达到极值的要求之中，不必单，不必单独列出，独列出，p惟一需考虑的仅是强制边界条件（第一类边界条件）的处理惟一需考虑的仅是强制边界条件（第一类边界条件）的处理，p进一步简化了方法的构造。进一步简化了方法的构造。8p有限元法的有限元法的主要特点主要特点是：是：n（1）离散化过程保持了）离散化过程保持了明显的物理意义明显的物理意义。p因为变分原理描述了支配物理现象的物理学中的最小作用原理因为变分原理描述了支配物理现象的物理学中的最小作用原理（如力学中的最小势能原理、
7、静电学中的汤姆逊定理等）。（如力学中的最小势能原理、静电学中的汤姆逊定理等）。p故基于问题固有的物理特性而予以离散化处理，列出计算公式，故基于问题固有的物理特性而予以离散化处理，列出计算公式，应当可应当可保证保证方法的正确性方法的正确性、数值解的存在数值解的存在与与稳定性稳定性等前提要素。等前提要素。9n（2）优异的解题能力优异的解题能力。p与其它数值计算方法相比较，有限元法在适应与其它数值计算方法相比较，有限元法在适应场域边界几何形场域边界几何形状状以及以及媒质物理性质变异媒质物理性质变异情况复杂的问题求解上，有突出的优情况复杂的问题求解上，有突出的优点。即方法应用不受上述二个方面复杂程度的
8、限制。点。即方法应用不受上述二个方面复杂程度的限制。p不同媒质分界面上的边界条件是自动满足的；二、三类边界条不同媒质分界面上的边界条件是自动满足的；二、三类边界条件不必作单独的处理件不必作单独的处理。p此外，此外，离散点配置比较随意离散点配置比较随意，并且取决于有限单元，并且取决于有限单元剖分密度剖分密度和和单元插值函数的选取单元插值函数的选取，可以充分保证所需的，可以充分保证所需的数值计算精度数值计算精度。10n(3)可方便地编写通用计算程序可方便地编写通用计算程序，使之构成模块化的子，使之构成模块化的子程序集合，适应计算功能延拓的需要，从而即可构成程序集合，适应计算功能延拓的需要，从而即可
9、构成各种高效能的计算软件包。各种高效能的计算软件包。n(4)从数学理论意义上讲，有限元法作为应用数学的一从数学理论意义上讲，有限元法作为应用数学的一个重要分支，很小有其它方法应用得这样个重要分支，很小有其它方法应用得这样广泛广泛。p它使微分方程的解法与理论面目一新，推动了泛函分析与计算它使微分方程的解法与理论面目一新，推动了泛函分析与计算方法的发展。方法的发展。11p有限元法的内涵也在不断延拓：有限元法的内涵也在不断延拓：n自从自从1969年以来，在流体力学领域中，通过运用年以来，在流体力学领域中，通过运用加权余量法导加权余量法导出的出的伽辽金法或最小二乘法伽辽金法或最小二乘法同样得到了有限元
10、方程。同样得到了有限元方程。n为提高数值解的计算精度，在为提高数值解的计算精度，在高阶有限元法高阶有限元法的应用范畴中，除了的应用范畴中，除了常用的常用的基于拉格朗日多项式构造基函数基于拉格朗日多项式构造基函数的的等参数有限元法等参数有限元法外，还外，还延拓构成了以延拓构成了以B样条函数基为基函数样条函数基为基函数的的B样条有限元法样条有限元法。pB样条有限元法的提出，不仅保证了以样条有限元法的提出，不仅保证了以位函数为待求量位函数为待求量的数值解的高的数值解的高精度，而且保证了与物理场特性相一致的场量数值解的连续性。精度，而且保证了与物理场特性相一致的场量数值解的连续性。12n把有限元法与把
11、有限元法与其它数值方法相结合其它数值方法相结合而构成的组合法，而构成的组合法，经常是解决特定问题的有效途径。经常是解决特定问题的有效途径。p例如，鉴于三维静态磁场分析的需要，由有限元法与数值积分例如，鉴于三维静态磁场分析的需要，由有限元法与数值积分法相组合而成的法相组合而成的单标量磁位法单标量磁位法，校正了三十余年来简化标量法，校正了三十余年来简化标量法有误的构造模式。有误的构造模式。n数学理论的发展也为有限元法注入了新的活力，数学理论的发展也为有限元法注入了新的活力，p1970年，以年，以A.M.Arthurs为代表提出了互补变分原理，为代表提出了互补变分原理，形成了泛函的所谓双边值问题，产
12、生了互补、对偶有限元法。形成了泛函的所谓双边值问题，产生了互补、对偶有限元法。这样，通过泛函极大与极小值问题的近似数值解，简单地求其这样，通过泛函极大与极小值问题的近似数值解，简单地求其算术平均值，即可获得充分逼近真实解的理想计算结果。算术平均值，即可获得充分逼近真实解的理想计算结果。135.2 变分原理变分原理p从介绍有关泛函、变分问题和变分法等数学概念从介绍有关泛函、变分问题和变分法等数学概念着手，阐述着手，阐述有限元法的变分原理有限元法的变分原理，p为有限元法基本原理的讨论提供必要的数学基础。为有限元法基本原理的讨论提供必要的数学基础。p以以加权余量法加权余量法为基础导出的为基础导出的伽
13、辽金有限元法伽辽金有限元法则将则将在矩量法中展述。在矩量法中展述。145.2.1 泛函与变分问题泛函与变分问题p数学上，通常变量与变量间的关系称为数学上，通常变量与变量间的关系称为函数函数，而，而泛函泛函则是函数集合的函数，也就是函数的函数。则是函数集合的函数，也就是函数的函数。p例如，静电场的势函数例如，静电场的势函数f f(r)是定义在坐标空间的是定义在坐标空间的函数集，系统电场总能量函数集，系统电场总能量U(r)则是定义在该则是定义在该函函 数集中的一个泛函，可记数集中的一个泛函，可记I(r)。15p与与多元函数的极值多元函数的极值问题相对应，在几何、力学上的求解问题相对应，在几何、力学
14、上的求解泛泛函极值函极值的问题。的问题。p最速降线问题。最速降线问题。n研究当质点从定点研究当质点从定点A自由下滑到定点自由下滑到定点B时，为使滑行时间最短，试时，为使滑行时间最短，试求质点应沿着怎样形状的光滑轨道求质点应沿着怎样形状的光滑轨道y=y(x)下滑。下滑。n取取A点为坐标原点，点为坐标原点，y轴竖直向下（图轴竖直向下（图5-1）。）。n则沿曲线则沿曲线y=y(x)滑行线段滑行线段ds所需的时间为所需的时间为16n因此滑行的总时间为因此滑行的总时间为n可见，积分值可见，积分值J=Jy(x)不仅取决于定积分的两端点不仅取决于定积分的两端点x1和和x2，而，而且取决于函数且取决于函数y=
15、y(x)的选择。的选择。n对照函数的定义，变量对照函数的定义，变量J值取决于函数关系值取决于函数关系y(x)，因此，因此J是是函数的函数的函数函数，是含义更为广泛的函数，故称之为是含义更为广泛的函数，故称之为函数函数y(x)的泛函的泛函，记作记作Jy(x)。n于是于是最速降线问题最速降线问题，在数学上，就归结为研究，在数学上，就归结为研究泛函泛函Jy(x)的极的极值问题值问题，即，即17p泛函的极值（极大值或极小值）问题就称为变分问题。泛函的极值（极大值或极小值）问题就称为变分问题。p对于一般问题，对应于一个自变量对于一般问题，对应于一个自变量x的的最简形式的泛函最简形式的泛函：n式中，式中，
16、F为单个自变量为单个自变量x、单个函数、单个函数y(x)及其导数及其导数y(x)的已知函的已知函数。数。p对照函数的定义，对照函数的定义，泛函泛函Jy的自变量的自变量不是一般的自变量，不是一般的自变量，而是一个或几个函数所属的函数族而是一个或几个函数所属的函数族y(x)。函数极值，求值；函数极值，求值；泛函极值，求函数泛函极值，求函数18p在最速下降问题，在端点在最速下降问题，在端点x1和和x2给定的无数个函数之中，给定的无数个函数之中，仅有一个函数仅有一个函数 能使式（能使式（5-2a）中的定积分达到极小值）中的定积分达到极小值函数，这一函数函数，这一函数 被称为被称为极值函数极值函数。p所
17、谓变分问题就在于寻求使泛函达到极值的所谓变分问题就在于寻求使泛函达到极值的该极值函数该极值函数，即分析研究泛函的极值问题。即分析研究泛函的极值问题。p物理学各分支都存在有相应的变分问题（变分原理），例物理学各分支都存在有相应的变分问题（变分原理），例如如n分析力学中的哈密顿原理、最小作用原理，分析力学中的哈密顿原理、最小作用原理，n静电学中的汤姆逊定理，静电学中的汤姆逊定理，n光学中的费尔马原理等。光学中的费尔马原理等。n最优控制论等，变分原理也有重要的应用。最优控制论等，变分原理也有重要的应用。195.2.2 泛函的变分与尤拉方程泛函的变分与尤拉方程p变分问题变分问题的的经典解法经典解法可归
18、纳为两大类。可归纳为两大类。n一类称为直接解法一类称为直接解法，其思想在于，其思想在于p把泛函的极值问题近似地转化为把泛函的极值问题近似地转化为一般多元函数的极值问题一般多元函数的极值问题，用，用有穷维子空间中的函数去逼近无穷维空间中的极值函数，从而有穷维子空间中的函数去逼近无穷维空间中的极值函数，从而近似地求得泛函的极值。近似地求得泛函的极值。p例如：瑞利例如：瑞利里兹法、康脱洛维奇法、伽辽金法等；里兹法、康脱洛维奇法、伽辽金法等；n另一类称为间接解法另一类称为间接解法，即，即p把变分问题转化为把变分问题转化为微分方程（所谓尤拉方程）的定解问题（边微分方程（所谓尤拉方程）的定解问题（边值问题
19、）值问题）来求解。来求解。20p有限差分法在理论上有限差分法在理论上没有以变分原理没有以变分原理为基础，因为基础，因而其收敛性和数值稳定性往往得不到保证而其收敛性和数值稳定性往往得不到保证。p有限元法有限元法是是里兹法与有限差分法相结合里兹法与有限差分法相结合的成果，的成果，它在理论上以变分原理为基础，在具体方法构造它在理论上以变分原理为基础，在具体方法构造上又利用了有限差分法网格离散化处理的思想。上又利用了有限差分法网格离散化处理的思想。21p有限元法的变分原理有限元法的变分原理n第一步与传统的变分法一样，首先把待求的第一步与传统的变分法一样，首先把待求的边值问题边值问题转化为等价转化为等价
20、的的变分问题变分问题，n然后通过然后通过有限单元剖分有限单元剖分的离散化处理，构造一个的离散化处理，构造一个分片解析分片解析的有限的有限元子空间，把变分问题近似地转化为有限元子空间中的元子空间，把变分问题近似地转化为有限元子空间中的多元函数多元函数极值极值问题，问题，n由此直接探求变分问题的近似解（极值函数解），以此作为所求由此直接探求变分问题的近似解（极值函数解），以此作为所求边值问题的近似解。边值问题的近似解。p下面从有限元法构造的需要出发，导出变分问题的解答下面从有限元法构造的需要出发，导出变分问题的解答（极值函数）所（极值函数）所必须满足的必要条件必须满足的必要条件尤拉方程尤拉方程。2
21、2最简泛函的变分问题最简泛函的变分问题p分析由式（分析由式（5-3）给定的最简形式泛函的变分问题，即）给定的最简形式泛函的变分问题，即n设想函数设想函数y(x)稍有变动，记作稍有变动，记作y+d dy（d dy称为函数称为函数y(x)的变分的变分，它反映了整个函数的变化量，显然有别于描述同一函数它反映了整个函数的变化量，显然有别于描述同一函数y(x)因因x变化而引起的函数增量变化而引起的函数增量D Dy），），n则泛函则泛函Jy的值也随之变化，其相应于函数的值也随之变化，其相应于函数d dy的的泛函增量泛函增量为为变分问题变分问题(泛函极值泛函极值)变分方程变分方程尤拉方程尤拉方程因讨论的是y
22、的变化23n设函数设函数F充分光滑，则由充分光滑，则由多元函数的泰勒公式多元函数的泰勒公式可将上式展成可将上式展成n式中，式中，p作为泛函增量作为泛函增量D DJ的的线性主部线性主部 称为称为泛函泛函Jy的一次变分（简称变分的一次变分（简称变分）。）。pd d2J、d d3J、分别是函数变分分别是函数变分d dy及其导数及其导数d dy的二次、三次、的二次、三次、齐次齐次式的积分，依次称为式的积分，依次称为二次变分二次变分、三次变分三次变分、24p设待求变分问题（设待求变分问题（5-4）的解答（）的解答（极值函数极值函数）为）为 y=y(x)（5-7）n设想函数设想函数y从极值解（从极值解（5
23、-7）稍稍变动到）稍稍变动到y+d dy，并把变分，并把变分d dy改记为：改记为：eheh(x)，pe e是一个任意给定的是一个任意给定的微量实参数（实变量）微量实参数（实变量）；ph h(x)是定义于区间是定义于区间x1，x2，且满足齐次边界条件的任意，且满足齐次边界条件的任意选定选定的可微函数，即有：的可微函数，即有：h h(x1)=h h(x2)=0。n则泛函则泛函Jy+eheh=Jy(x,e e)=j j(e e)成为参数成为参数e e的函数，的函数，p因为参数因为参数e e的值确定了的值确定了y=y(x,e e)函数族里的曲线，所以同时也就确定了泛函函数族里的曲线，所以同时也就确定
24、了泛函Jy(x,e e)的值，的值，p而且当而且当e e=0时泛函即获得极值函数的解时泛函即获得极值函数的解。e e是对是对y的增量，不是对的增量，不是对x.对于泛函，不妨把对于泛函，不妨把x看成常数看成常数因y是x的函数，但讨论的是y的变化25p根据微积分学可知，函数根据微积分学可知，函数j j(e e)在在e e=0时取得极值的时取得极值的必要条件是必要条件是n由于由于n因此因此n式中式中j即为J26n故可得故可得n简写为简写为n将上式与式（将上式与式（5-6）相比较，只相差一个数值因子）相比较，只相差一个数值因子e e。27n故故（5-8）等价于）等价于变分方程变分方程n也即也即n利用利
25、用分部积分分部积分，根据，根据变分与微分顺序可以互换变分与微分顺序可以互换的原理，即的原理，即d dy=(d dy)，得，得（线性主部）（线性主部）28n在变分问题中，在变分问题中，变分变分d dy在端点保持为零在端点保持为零n于是，必要条件（于是，必要条件（5-12）成为）成为n由于上式对任意的由于上式对任意的d dy都成立，所以极值函数必须满足以下微分方程都成立，所以极值函数必须满足以下微分方程n这个方程就称为这个方程就称为对应于由式（对应于由式（5-4）给定的最简泛函极值问题的尤）给定的最简泛函极值问题的尤拉方程。拉方程。29(1)多变量泛函的变分问题之一多变量泛函的变分问题之一（第一类
26、边值问题）（第一类边值问题）p若设待求极值函数若设待求极值函数u(r)=u(x,y,z)为一多变量的函数，则其为一多变量的函数，则其与式与式（5-4）对应的变分问题可表示为对应的变分问题可表示为n同前理，设同前理，设函数变分函数变分d du(r)=eheh(r)，n并注意到并注意到p如如等等30n则其泛函变分为则其泛函变分为n通过应用如下向量恒等式通过应用如下向量恒等式n以及高斯散度定理，即以及高斯散度定理，即n且考虑到面元且考虑到面元类似于类似于(5-6)由泰勒展开由泰勒展开31n则由（则由（5-15），极值函数），极值函数u(r)须满足如下的须满足如下的变分方程变分方程n因为因为e e为微
27、量实参数，由为微量实参数，由h h(r)的随意性，可知上述变分方程（的随意性，可知上述变分方程（5-16）的解必）的解必等价于等价于以下边值问题的解答，即以下边值问题的解答，即n式式（5-17a）称为对应于式（称为对应于式（5-14）给定的变分问题的）给定的变分问题的尤拉方程尤拉方程，而式（而式（5-17b）即为相应的边界面即为相应的边界面s上所给定的上所给定的边界条件边界条件。32p可见，对于可见，对于“固定端点固定端点”的变分问题的变分问题，即，即n应应等价于等价于如下的如下的第一类边值问题第一类边值问题（此时，（此时，）：）：因变分方程不含第一类边界条件因变分方程不含第一类边界条件,故应
28、加上故应加上.33p例例5-1 静电学中汤姆逊定理的数学描述静电学中汤姆逊定理的数学描述。n变分原理的应用实质上是对物理学定律的一种重新描变分原理的应用实质上是对物理学定律的一种重新描述。述。n对于电、磁场边值问题而言，静电学中的汤姆逊定理对于电、磁场边值问题而言，静电学中的汤姆逊定理即是描述静电现象的即是描述静电现象的“最小作用原理最小作用原理”。n定理指出，定理指出，处于介质中一个固定的带电导体系统，其处于介质中一个固定的带电导体系统，其表面上电荷的分布，应使合成的静电场具有最小的静表面上电荷的分布，应使合成的静电场具有最小的静电能量电能量。34n任一由任一由n个带电导体构成的个带电导体构
29、成的二维二维具平行平面场特征具平行平面场特征静电场静电场的能量积的能量积分式为分式为p其中其中n对照式（对照式（5-14），），（5-20a）所示的能量积分所示的能量积分即是一类取决于二元即是一类取决于二元电位函数电位函数j j(x,y)分布的分布的泛函泛函。n因而，根据因而，根据汤姆逊定理汤姆逊定理，对于二维静电场，其规律性就归结为下，对于二维静电场，其规律性就归结为下述述变分问题变分问题：35n故其极值函数故其极值函数j j(x,y)的解答应满足的解答应满足尤拉方程（尤拉方程（5-19a），），并在并在边界上边界上应满足条件（应满足条件（5-19b）。p将式（将式（5-21a）中函数）中函
30、数F进行尤拉方程（进行尤拉方程（5-19a）相关各项的运算，即得）相关各项的运算，即得与变分问题与变分问题(5-21a)对应的尤拉方程为对应的尤拉方程为n可见，由可见，由条件变分问题条件变分问题（5-21a）、）、(5-21b)给出的极值函数给出的极值函数j j(x,y)应应满足具有给定边值（满足具有给定边值（5-21b）的）的拉普拉斯方程拉普拉斯方程（5-22）。）。n显然，式（显然，式（5-22）和（）和（5-21b）一起即构成众所周知的）一起即构成众所周知的第一类边值问第一类边值问题题，而在物理意义上该边值问题即表征为静电学中的，而在物理意义上该边值问题即表征为静电学中的汤姆逊定理汤姆逊
31、定理。36p与例与例5-1相仿，通过尤拉方程，可知与下述变分问题相仿，通过尤拉方程，可知与下述变分问题n等价的边值问题为等价的边值问题为n上式即为二维具有平行平面场特征的上式即为二维具有平行平面场特征的泊松方程的第一类边值问题泊松方程的第一类边值问题。37(3)多变量泛函的变分问题之二多变量泛函的变分问题之二（第二、三类边值问题）（第二、三类边值问题）p若在式（若在式（5-14）所示变分问题的基础上，再添加一项由）所示变分问题的基础上，再添加一项由面积分面积分表述的表述的泛函泛函，即令，即令p式中，积分面式中，积分面s即为积分域即为积分域v的界面。的界面。n则同前理，参照式（则同前理，参照式（
32、5-16）导出过程，可得其）导出过程，可得其变分方程变分方程为为38n上述变分方程的解必上述变分方程的解必等价于以下边值问题的解答等价于以下边值问题的解答，即，即39p可见，若设变分问题为可见，若设变分问题为n则根据变分方程则根据变分方程d dJ=0，由其对应的尤拉方程的定解问题，由其对应的尤拉方程的定解问题式（式（5-27a）和）和(5-27b)，可知与上述变分问题（，可知与上述变分问题（5-28）等价的边值问题）等价的边值问题为为n上式即为上式即为泊松方程泊松方程的的第三类边值问题第三类边值问题。40p若令第三类边界条件若令第三类边界条件式（式（5-29b）中的中的f1(rb)=f2(rb
33、)=0，则得如下泊松方程的齐次则得如下泊松方程的齐次第二类边值问题第二类边值问题：n不难看出，其等价变分问题可基于式（不难看出，其等价变分问题可基于式（5-28），通过令该式中），通过令该式中f1=f2=0，得知为，得知为41p可见，可见，n第二或第三类边界条件第二或第三类边界条件在变分问题中被包含在泛函达在变分问题中被包含在泛函达到极值的要求之中到极值的要求之中，不必单独列出。，不必单独列出。n还可证明，还可证明，场域内不同媒质分界面上的场域内不同媒质分界面上的边界条件边界条件也包也包含在泛函达到极值的要求之中含在泛函达到极值的要求之中，且系自动满足，不必，且系自动满足，不必另行处理。另行处
34、理。n因此，常称这些边界条件为因此，常称这些边界条件为自然边界条件自然边界条件，而相应的，而相应的变分问题（变分问题（5-28）或（）或（5-31）称为）称为无条件变分问题无条件变分问题。42p但对于但对于第一类边界条件第一类边界条件，则在变分问题中与在边，则在变分问题中与在边值问题中一样，必须作为定解条件列出值问题中一样，必须作为定解条件列出。n换句话说，变分问题的极值函数解必须在满足这一边换句话说，变分问题的极值函数解必须在满足这一边界条件的函数类中去寻求。界条件的函数类中去寻求。n因此，称这类边界条件为因此，称这类边界条件为强制边界条件强制边界条件，而相应的变，而相应的变分问题（分问题（
35、5-18a）、（）、（5-18b）或（）或（5-21a）、（）、（5-21b）或（或（5-23a）、（）、（5-23b）称为）称为条件变分问题条件变分问题。43p可以看出，可以看出，n与与能量积分对应的泛函能量积分对应的泛函Jj j，二次地依赖于函数，二次地依赖于函数j j及其及其偏导数，故又称偏导数，故又称Jj j为函数为函数j j的二次泛函的二次泛函，n而相应的而相应的变分问题变分问题即称为即称为二次泛函的极值问题二次泛函的极值问题。p对于对于磁场问题磁场问题，n同样可构成由标量磁位同样可构成由标量磁位j jm或向量磁位或向量磁位A所描述的二次泛所描述的二次泛函函Jj jm或或JA，得到与
36、给定边值问题等价的变分问题。，得到与给定边值问题等价的变分问题。44函数变量函数变量j j的求取的求取p由上可见，由上可见，n求泛函的极值和解欧拉方程，在数学上都可以代表同求泛函的极值和解欧拉方程，在数学上都可以代表同一个物理问题。一个物理问题。n对两者求得近似解都具有同样的效果。对两者求得近似解都具有同样的效果。n但是在实际计算中，对后者求解往往是困难的，而对但是在实际计算中，对后者求解往往是困难的，而对前者求近似解则常常并不太困难。前者求近似解则常常并不太困难。n因此常利用对泛函离散化的有限元法求解场函数因此常利用对泛函离散化的有限元法求解场函数j j。45p有限元法须先估计有限元法须先估
37、计j j函数的形式。一般来讲函数的形式。一般来讲,精确估计出此精确估计出此泛函在极值情况下泛函在极值情况下j j函数的形式是不可能的。函数的形式是不可能的。p但是原则上可以采用猜测出的函数近似表示但是原则上可以采用猜测出的函数近似表示j j(x,y,x,),其中其中 对应于对应于N个未知的参数个未知的参数q qi,再计算泛函再计算泛函I(j j)，然后，然后用取最小值的条件用取最小值的条件得到得到N个方程个方程,这个方程组可用来求出参数这个方程组可用来求出参数q qi的解。的解。46p有限元法是将网络节点上的函数有限元法是将网络节点上的函数j j的离散值作为参数的离散值作为参数,网络元网络元(
38、剖分元剖分元)内的函数值则采用多项式插值从周围临近节点上的这些参数值求出。内的函数值则采用多项式插值从周围临近节点上的这些参数值求出。p例如例如,选择用三角形元将求解区域划分为子区间的网络，对泛函选择用三角形元将求解区域划分为子区间的网络，对泛函I(j j)求极小值求极小值,就得到节点上未知的势函数的值，然后采用线性插值法，就得到节点上未知的势函数的值，然后采用线性插值法，则可以求出在一个三角形元内的任意一点则可以求出在一个三角形元内的任意一点(x,y)上的函数值。上的函数值。p有限元法的最后解是函数在这些节点上的估计值。有限元法的最后解是函数在这些节点上的估计值。p由于用来求泛函极小值的函数是近似的线性迭代函数由于用来求泛函极小值的函数是近似的线性迭代函数,因而所得到的因而所得到的节点上的函数值并不是精确解。该截断误差可以通过减小元的尺寸或节点上的函数值并不是精确解。该截断误差可以通过减小元的尺寸或提高迭代函数的阶数来降低。提高迭代函数的阶数来降低。47