有限元原理加权余量法和变分法ppt课件

1、l第三讲l1.偏微分方程求解有限元法的原理加权余量法和变分法1.解析法2.运用范围有限，适用于实际求解，但有剧烈的物理含义常系数微分方程3.某些复杂问题，很思索根本找不到解析解4.2. 数值法5.工程实践中运用广泛，复杂场域问题，但物理含义不很清楚。任何问题总可以找到数值解数学方法222222tJtAA.l2.数值求解方法2/41. 根本思想： 以偏微分方程的近似解来替代其真解，只需近似解与真解足够以偏微分方程的近似解来替代其真解，只需近似解与真解足够接近，就可以近似解作为问题的解，并满足足够的精度。接近，就可以近似解作为问题的解，并满足足够的精度。2. 根本方法：1.假设一个近似解，该解为一

2、组方式上简单函数假设一个近似解，该解为一组方式上简单函数 的线性组合的线性组合来表示，线性组合的系数就是一组待定系数来表示，线性组合的系数就是一组待定系数2.然后建立一种思索了微分方程和边境条件的关于真解然后建立一种思索了微分方程和边境条件的关于真解 和近似解和近似解间误差的目的函数间误差的目的函数 F3.用适当的算法使得该目的函数最小化用适当的算法使得该目的函数最小化最小化的过程就确定了最小化的过程就确定了待定系数，从而也就得到了问题的近似解。待定系数，从而也就得到了问题的近似解。 i iC尝试函数，基函数，形函数.l2.数值求解方法2/4目的函数最小化的目的：一方面，使得近似解最大程度接近

3、真解；目的函数最小化的目的：一方面，使得近似解最大程度接近真解； 另一方面，求得构成近似解的待定系数。另一方面，求得构成近似解的待定系数。数学上，构成目的函数的方法很多，不同的构成方法就构成了不同的数学上，构成目的函数的方法很多，不同的构成方法就构成了不同的数值解法，电磁场中就常见的是：加权余量法和变分法。数值解法，电磁场中就常见的是：加权余量法和变分法。.l3.电磁场位函数偏微分方程的数值求解方法加权余量法电磁场问题总可以用位函数的偏微分方程和相应的边境条件表述电磁场问题总可以用位函数的偏微分方程和相应的边境条件表述222222tJtAA ) 1(1g)2()2(22ht两个偏微分方程方式一

4、样，故以电位方程的求解过程为例。磁位矢两个偏微分方程方式一样，故以电位方程的求解过程为例。磁位矢量的方程可以分解到各个分量上变为标量方程。量的方程可以分解到各个分量上变为标量方程。.在求解场域内，偏微分方程的真解为在求解场域内，偏微分方程的真解为 ，近似解为，近似解为 它由一组简单函数它由一组简单函数 的线性组合表达，表达中有待定系数的线性组合表达，表达中有待定系数 即：即：l3.电磁场位函数偏微分方程的数值求解方法加权余量法加权余量法加权余量法 i iC 1niiiC简单函数，普通选用简单方式的函数，一旦选定就是知的了待定系数是真正的求解目的问题的自在度近似解.l3.电磁场位函数偏微分方程的

5、数值求解方法加权余量法加权余量法就是一种定义近似解与真解之间误差即余数，并设加权余量法就是一种定义近似解与真解之间误差即余数，并设法使其最小的方法。法使其最小的方法。 : 22）（）（上：边界内场域RR加权余量法误差即余数的定义：加权余量法误差即余数的定义：留意：普通余数并不表示近似解与真解间的代数差场域内，加权余留意：普通余数并不表示近似解与真解间的代数差场域内，加权余量法的采用拉普拉斯算子作用后的差别即余数，来代表近似解量法的采用拉普拉斯算子作用后的差别即余数，来代表近似解整体接近偏微分方程真解的程度。整体接近偏微分方程真解的程度。问题的自在度.l3.电磁场位函数偏微分方程的数值求解方法加

6、权余量法当余数小于要求的精度时，就可以以为近似解就是偏微分方程的解。当余数小于要求的精度时，就可以以为近似解就是偏微分方程的解。要减少余数，我们可以经过寻求适当的待定系数来实现。要减少余数，我们可以经过寻求适当的待定系数来实现。为有效表达减小余数的效果，还选取适当的加权函数，以使余数和该加为有效表达减小余数的效果，还选取适当的加权函数，以使余数和该加权函数的积分为权函数的积分为0。“加权余量法的来由。加权余量法的来由。 ; *jjww设加权函数为：.l3.电磁场位函数偏微分方程的数值求解方法加权余量法 ,.2 , 1 d d *jRwRwjj，目标函数：加权余数的定义：加权余数的定义：加权函数

7、的选取方法很多：如点重合、子域重合、最小二乘法、迦辽金法。加权函数的选取方法很多：如点重合、子域重合、最小二乘法、迦辽金法。效果较好的、运用较多的是迦辽金法：效果较好的、运用较多的是迦辽金法：jjjww*即：迦辽金法选取尝试函数本身为加权函数即：迦辽金法选取尝试函数本身为加权函数.l3.电磁场位函数偏微分方程的数值求解方法加权余量法 0 d d )()(趋于则余数最小，令，RjjjRjFRRF由此构建加权量法的目的函数：由此构建加权量法的目的函数：上述过程中，曾经将偏微分方程转化为上述过程中，曾经将偏微分方程转化为j个代数方程组，便于计算机求解。个代数方程组，便于计算机求解。关于函数的函数，称

8、为：泛函数，或泛函.l3. 加权余量法例1例例1.两极电容板内部电场分布问题：两极电容板内部电场分布问题：根据问题特点将根据问题特点将3维问题简化为维问题简化为2维，维，进一步简化为进一步简化为1维。维。该问题是静态电场问题，该问题是静态电场问题，偏微分方程和边境条件：偏微分方程和边境条件：;10; 0002d.22112211211 xCxCCCxCCiiiniii加权余量法求解：加权余量法求解：1.选取尝试函数、构造近似解：选取尝试函数、构造近似解：1,2)(i iix2.结合问题，写出余数表达式：结合问题，写出余数表达式：22C R :22l3. 加权余量法例1实际上恣意选取，操作中越简

9、单越好 020 )()()(222221122122CxCxCxCiii.2.结合问题，写出余数表达式：结合问题，写出余数表达式： ）（）（：Rl3. 加权余量法例1 )xCxC(dx )xCxC( xdxdx0 xx2211221100处：在处：在）（）（221121xCxCxCiii）（ 10 0 00处：在处：在）（）（dxxdxx )dCdC(Rdx R xdxx100022110处：在处：在.3. 加权余数表达式：加权余数表达式：2 , 1 d d )(jRRFjjRj，l3. 加权余量法例1010110010021221232212222110221102dC)d(dCd)ddCd

10、C(dC d )xCxC( ( x d )xCxC( ( x d)C( x d Rd RF,jdx |dxx |0 xd11)R( 1得到一个代数方程：时.3. 加权余数表达式：加权余数表达式：l3. 加权余量法例1010)32( )10(032 d )10)( ( d )0)( ( d)2( d d ,2223132423132 |x221120 |0 x2211202222)(2dCddCdddCdCdCxCxCxxCxCxCxRRFjdxdxdR又得到一个代数方程：时.4. 求解上述两个代数方程组，得到待定系数，从而确定近似解求解上述两个代数方程组，得到待定系数，从而确定近似解l3. 加

11、权余量法例10 /10 21；解得：CdCxdxCxCxCiii10221121近似解：）（加权余量法求解流程：加权余量法求解流程：1.选取尝试函数、构造近似解选取尝试函数、构造近似解2.结合问题，写出余数表达式结合问题，写出余数表达式3. 写出加权余数表达式写出加权余数表达式4. 令各加权余数表达式为令各加权余数表达式为0，得到代数方程组，解之得到待定，得到代数方程组，解之得到待定系数，从而确定近似解系数，从而确定近似解.该静态电场问题的真解解析解：该静态电场问题的真解解析解：l3. 加权余量法例1真解与近似解一样是由于尝试函数选择的刚好，通常是有差别的，如选用三角函数，但求解过程会复杂，可

12、见尝试函数的选取是有技巧的。.l4. 加权余量法求解普通化偏微分方程的归纳 )( )(sq 普通化偏微分方程：线性微分算子那么其他数为：sRqR)()()()()()( 1niiiC 其中：令加权余数为0，构建代数方程：0d )(d )(0d )(d )(1*1*)(sCwqCwswqwFniiijniiijjjRj .l4. 加权余量法求解普通化偏微分方程的归纳由于是线性微分算子，故微分、求和、积分次序可互换,代数方程变形： d d d)(d)(*1*1swqwCwCwjjniiijniiij 0d )(d )(1*1sCwqCwniiijniiij d swd qwCd)(wd)(w*jj

13、niii*jij1有j个代数方程，通常等于待定系数个数.l4. 加权余量法求解普通化偏微分方程的归纳代数方程写成矩阵方式：d d d)(d)(*1*swqwCwwjjniiijij系数鼓励边境条件bFCK系数矩阵nn待定系数矩阵、源矩阵、边境矩阵n1矩阵元素值： d d ddswbqwFwwKjjjjijijji*)()( 虽然元素值还需求积分、微分的求得，还难以借助计算机求解，但至少化为了代数方程组。经过选择适宜的加权函数和尝试函数可以大大简化矩阵元素的矩阵方程。有限元方法就是如此.l5. 加权余量法的进一步优化边境条件的处置22112 hngq适当的选取加权函数，并对加权余数积分进展处置，

14、可使某些边境条件从加权余数的表达式中消逝，从而简化矩阵方程及其系数的求解。以有源静电场问题为例帕松方程.由近似解表述的加权余数为：l5. 加权余量法求解普通化方法的进一步优化 d d RwRwFjjRj*)(留意余数的本质可使上式第二项消失动满足，则第一类边界条件自，使得常常选取）可简化计算但适当的选取（作限制意选，理论上尝试函数可任其中近似解：gCiniii, 1 d d )()(*）（）（ jjww d d d 2112)()()(*hnwgwqwjjj ）（ d d d 2211122)()()(*）（）（）（ nwwwjjj.经过尝试函数，简化加权余数后：l5. 加权余量法求解普通化方

15、法的进一步优化njhwnwqwwhnwqwFjjjjjjRj,.,)()(*)(321d d d d d d 22222 上式第一项，由格林第一定律得： d d d d 212nwnwwwjjjj 降了微分阶数，等于降了近似解尝试函数的延续性要求，从而扩展了其选择范围.代入后：l5. 加权余量法求解普通化方法的进一步优化，则上式被大大简化如选取加权函数：*)(jjjjjjjjjjjjjjjRjwwhwqwwhwnwqwnwnwwhwnwqwwF d d d d d d d d d d d d d 22221222 由于近似解在1类边境上常数，所以此项为0选取特殊加权函数后，两项和为0第二类边境

16、条件也消逝了，阐明曾经自动满足了.令加权余数为0即可得到求解原微分方程的一组代数方程：l5. 加权余量法求解普通化方法的进一步优化njhwqwwFjjjRj,.3 , 2 , 10d d d 2)(这里加权函数只需一个了，进一步，用迦辽金法，选加权函数为尝试函数本身 1niiijjCw，且有近似解表达式：d d d 21hqCjjniiij）（.l5. 加权余量法求解普通化方法的进一步优化d d d 21hqCjjniiij）（由于是线性微分算子，故微分、求和、积分次序可互换,代数方程变形： d d d 21hqCjjiniijd d d)(d)(*1*swqwCwwjjniiijij对比简化

17、前的代数方程：曾经大大简化，关键是边境条件项全部消逝，微积分计算也降阶、简化.l5. 加权余量法求解普通化方法的进一步优化代数方程写成矩阵方式： d d d 21hqCjjiniij njnjnjnnnjnninijiinjbbbfffccckkkkkkkkkkkk1112121111211 d d d 2hbqfkkjjjjijjiij 对称矩阵，简化计算还有积分求和，梯度差分，有限元将作处置小结：简化后1、2类边境条件自动满足； 尝试函数、加权函数选取 微分降阶，简化计算 对称矩阵，简化计算 根据情况源矩阵、边境矩阵能够为0对拉普拉斯方程和帕松方程问题适宜.l6. 简化后加权余量法 例2例

18、1中的静电场问题，变为两电极板接地，中间充溢电荷。帕松方程.34231201xaxaxaxa加权余量法求解：加权余量法求解：1.初选尝试函数、构造近似解：初选尝试函数、构造近似解：4)31,2(i :3210，、xxxxi利用问题，对近似解进展简化，对尝试函数进展优化利用问题，对近似解进展简化，对尝试函数进展优化l6. 简化后加权余量法 例2经过尝试函数的选取，近似解满足1类边境条件，使得1类边境条件在方程中消逝324110 00|aaaaxx；由此，尝试函数和近似解优化为：由此，尝试函数和近似解优化为：.)()(322312211xxcxxccc 2. 修正尝试函数，以满足修正尝试函数，以满

19、足1类边境条件：类边境条件：)(作好准备、其梯度为：、 3231222132231xxxxxxx l6. 简化后加权余量法 例2.3.代公式计算矩阵元素代公式计算矩阵元素 边境矩阵边境矩阵b为为0l6. 简化后加权余量法 例2.4. 封装矩阵：封装矩阵：l6. 简化后加权余量法 例25. 求解矩阵，得近似解：求解矩阵，得近似解：.该有源静态电场问题的真解解析解：该有源静态电场问题的真解解析解：l6. 简化后加权余量法 例2真解与近似解一样是由于尝试函数选择的刚好，通常有差别。如例3.l7. 简化后加权余量法 求解普通化的微分方程 例3偏微分方程描画的问题如下：212|21 2212 xxdxd

20、xxdxdxdxd| )()( ），（加权余量法求解：加权余量法求解：1.初选尝试函数、构造近似解：初选尝试函数、构造近似解：34231201xaxaxaxa 4)31,2(i 3210，、 xxxxi: .利用问题及其边境条件，对尝试函数进展优化使近似解满足边境条件利用问题及其边境条件，对尝试函数进展优化使近似解满足边境条件经过尝试函数的选取，近似解满足1类边境条件，使得1类边境条件在方程中消逝2243211aaaax | 两个方程，两个独立未知数，消两个方程，两个独立未知数，消a1、a2，重定尝试函数，边境条件自动满足，重定尝试函数，边境条件自动满足，简化求解过程简化求解过程l7. 简化后

21、加权余量法 求解普通化的微分方程 例32124823221432224322 aaaxaxaaxdxdxxx|)(| )( .2. 修正尝试函数，以满足修正尝试函数，以满足1、2类边境条件：类边境条件：)()()()()(111314491111314492212321 xxxCxxCxxxxxxx且近似解为：、 l7. 简化后加权余量法 求解普通化的微分方程 例302102| 021 2212 xxdxdxRxdxdxdxdR| )(:)()( 条件自动满足尝试函数的选取，边界），（余数为：余数为：.l7. 简化后加权余量法 求解普通化的微分方程 例32221222124331441 211

22、1314491xxCxCxdxdxdxdRxxxCxxCx )()()()()()(代入： 结合问题，余数的详细表达式为：结合问题，余数的详细表达式为：问题的加权余数目的泛函为：问题的加权余数目的泛函为： 0 d 24331441 d d 2221 )()()(xxCxCRRFjjjRj .4. j=2,3时得代数方程：时得代数方程：5. 求解矩阵，得近似解：求解矩阵，得近似解：l7. 简化后加权余量法 求解普通化的微分方程 例3 0d 24331441 31 d 24331441 21222121222122 )()()()()()()(xxCxCxxxxCxCFR 0d 24331441

23、111 d 24331441 212221221222133 )()()()()()()(xxCxCxxxxxCxCFR .5. 求解矩阵，得待定系数和近似解：求解矩阵，得待定系数和近似解：l7. 简化后加权余量法 求解普通化的微分方程 例3)(.)(.)(.11134710311378244913471013782221 xxxxxxCC、 .真解解析解：真解解析解：l7. 简化后加权余量法 求解普通化的微分方程 例3.l8. 归纳加权余量求解偏微分方程步骤加权余量法求解流程：加权余量法求解流程：1.初步选取尝试函数、构造近似解初步选取尝试函数、构造近似解2.结合问题的边境条件对尝试函数进展

24、修正，以简化求解结合问题的边境条件对尝试函数进展修正，以简化求解3.写出余数表达式写出余数表达式3. 写出加权余数表达式迦辽金方法选取加权函数写出加权余数表达式迦辽金方法选取加权函数4. 令权余数表达式在各尝试函数下为令权余数表达式在各尝试函数下为0，得到代数方程组，解，得到代数方程组，解之得到待定系数，从而确定近似解之得到待定系数，从而确定近似解.l8. 归纳加权余量求解偏微分方程步骤加权余数法求解普通性偏微分方程的方法：方程的近似解被表示为一系列独立的尝试函数的线性组合，其中包括未知的待定系数。通常用迦辽金原理选取加权函数，即令加权函数等于尝试函数本身，从而完成对加权余数的定义，(尝试函数

25、的选取满足边境条件)经过对加权函数在区域内和在边境上的积分使其平均值为零，也就是说，使近似解与准确解之间的差别在某种目的下到达最小化。如此可以构成一个矩阵方式的代数方程组，求解该矩阵方程可以确定待定系数从而得到偏微分方程的独一近似解。.l9. 变分法简介另外一种求解偏微分方程的普通方法，即变分法。变分法与加权余数法类似，近似解也用一系列线性独立的尝试函数表示包括未知的待定系数。与加权余数法不同的是，变分法用另外的方法来构成求解待定系数的矩阵方程。在变分法中，首先要构成一个近似解的函数，称为泛函。从广义来说，加权余数积分(即平均值)也是一种泛函。然后使该泛函最小化，从而减小近似解的误差。普通说来，要找到一个适宜于偏微分方程及边境条件的泛函是一项难度很大的任务。由于前人已做了许多研讨任务，已找到了适宜于许多常见方式的偏微分方程的泛函。对于电磁场方程来说，偏微分方程常具有拉普拉斯、帕松和赫姆