第五章静态场的边值问题课件.ppt(70页)
- 【下载声明】
1. 本站全部试题类文档,若标题没写含答案,则无答案;标题注明含答案的文档,主观题也可能无答案。请谨慎下单,一旦售出,不予退换。
2. 本站全部PPT文档均不含视频和音频,PPT中出现的音频或视频标识(或文字)仅表示流程,实际无音频或视频文件。请谨慎下单,一旦售出,不予退换。
3. 本页资料《第五章静态场的边值问题课件.ppt(70页)》由用户(ziliao2023)主动上传,其收益全归该用户。163文库仅提供信息存储空间,仅对该用户上传内容的表现方式做保护处理,对上传内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(点击联系客服),我们立即给予删除!
4. 请根据预览情况,自愿下载本文。本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
5. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007及以上版本和PDF阅读器,压缩文件请下载最新的WinRAR软件解压。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 第五 静态 边值问题 课件
- 资源描述:
-
1、1第五章 静态场的边值问题2标题添加点击此处输入相关文本内容点击此处输入相关文本内容前言点击此处输入相关文本内容标题添加点击此处输入相关文本内容3边值问题研究方法解析法数值法分离变量法镜像法复变函数法有限差分法有限元法边界元法矩量法模拟电荷法45.1 唯一性定理和解的叠加原理 一.唯一性定理1、表述 在给定的区域内,泊松方程(或拉普拉斯方程)满足所给定的全部边界条件的解是唯一的。2、边界条件的形式给定全部边界上的函数值给出全部边界上函数的法向导数值给定部分边界上的函数值,而其余边界上给出函数的法向导数值“狄利赫利”边界条件“聂曼”边界条件混合边界条件1Cs2Cns2121,CnCss53、唯一
2、性定理的证明证明:考虑泊松方程,用反证法 在区域内存在两个不同的函数 和 都满足相同的泊松方程 ,并且在区域边界S上满足同样的边界条件。12f2令210)(22122122ffuFFuFu)(22()()则有利用Fu,取则对上式两边在区域内作体积分,然后运用散度定理,得dsds2)()(6ddsns2)(dsnsd)(将 代入上式得 0)(112121CCssss0)(222121CCnnnnssss0,021ssn1Cs对于第一类边界条件2Cns对于第二类边界条件对于第三类边界条件不论对哪类边界条件,面积分 都等于零 sdsn 因此有 0)(2d7由于 恒为正值,故上式成立的条件是2)(0C
3、解之可得讨论:对于第一类边界条件,0S因为所以C=0第二类和第三类边界条件的情况 C21对于求解场函数来说,解是唯一的。4、应当明确只有在区域的所有边界上给出唯一的边界条件时,边值问题的解才是唯一确定的。唯一性定理 给出了求解电磁场问题的理论依据不论采用什么方法,只要得到的解能够在区域内满足方程而在边界上满足边界条件,这个解就是该边值问题的唯一正确解。8二.解的叠加原理解的可叠加性是方程线性的必然结果。1、表述对拉普拉斯方程,若 、都是满足方程 的解 n1202nnaaa2211则其中ai为任意常数。也是方程 的解 02并且 、都是满足方程 的解 n12201122pnnaaa则其中ai为任意
4、常数。也是方程 的解 2f对泊松方程,若 是满足方程 的一个任意解 2fp92、证明221122()pnnaaa讨论泊松方程的情况对叠加得到的结果两边作 运算,得222221122pnnaaa12000nfaaaf 因此 是方程 的解 2f令 f=0,即可得到拉普拉斯方程情况的证明3、应用求解边界问题时,可以先将复杂边界条件分解成便于求解的几个边界条件,则总的边界问题解就是这些解的叠加。1012321230sssCCC 111222333222123112312231233000000000sssssssssCCC 例:分解为三个边界问题分解后每个边值问题都只有一个非齐次边界值,求解变得容易。
5、123原问题的解应该是三个问题解的叠加115.2 拉普拉斯方程的分离变量法一.直角坐标系中的分离变量法 直角坐标系中拉普拉斯方程的表达式为 0222222zyx)()()(),(zhygxfzyx令 ,并代入上式0)()(1)()(1)()(1222222zzhzhyygygxxfxf)()()(zhygxf并两边同除以 得 1、方法介绍2xk2yk2zk变量的分离12则上式分解成三个独立的全微分方程,即 0)()(222xfkdxxfdx0)()(222ygkdyygdy0)()(222zhkdzzhdz222,zyxkkk称为分离常数,分离常数之间满足约束关系 0222zyxkkk12()
6、f xAxA12()sincosxxf xAk xAk xxkAxkAxfxxchsh)(2112xxk xk xAeA e全微分方程解的选取(以 为例)xk20 xk 20 xk 220 xxkk13对 和 也都有与上述相同形式的解。在 、三组可取解中各取其一并相乘,即可得到一个解的表达式。)()()(),(zhygxfzyx)chsh)(cossin)(cossin(21121zkCzkCykBykBxkAxkAzzyyyxx)(yg)(zh)(xf)(yg)(zh解的选取并不是任意的,因为存在约束条件2220 xyzkkk20 xk 20yk 20zk 如果取 ,则必须取 a.对于有两个
7、零值边界的方向,其对应的函数取三角函数;b.对于单零值边界方向,对应的函数一般取双曲函数形式;c.而有无限远边界的方向,一般取指数函数形式。d.若位函数与某一坐标变量无关,则该变量对应的函数取成常数,并取作1。根据边界条件来选择函数:14拉普拉斯方程的通解),3,2,1(,nikkkziyixi本征值:满足齐次边界条件的分离常数可以取一系列特殊值 本征值对应的函数称为本征函数或本征解。)()()(),(),(11niiiiniizhygxfzyxzyx所有本征解的线性叠加构成满足拉普拉斯方程的通解在许多问题中,单一本征函数不能满足所给的边界条件,而级数形式的通解则可以满足单个解函数所无法满足的
8、边界条件。152、例题例5.1 边界条件如图,求电位分布 o y 0U b a 0UU 图51 长方形截面的导体长槽 0U 0U x 解:内部电位与z无关,是二维问题0),(),(2222yyxUxyxU)()0(00abyUx)()0(0bbyUax)()0(00caxUy)()0(0daxUUby 根据边界条件写出通解x方向电位有两个零值边界,应取三角函数形式;y方向电位为单零值边界,应取双曲函数形式。)(xf)(yg16222kkkyx令分离常数 则电位通解为iiiiiiiiiykBykBxkAxkAyxUchshcossin),(2121求待定系数将边界条件(a)代入得)0(chsh0
9、212byykBykBAiiiiii若对任意y成立,则 02iA112(,)sinshchiiiiiiiU x yAk x Bk yBk yiiiiiiiykBykBakAchshsin0211再利用边界条件(b)得 通解变为若对任意y成立,则),3,2,1(iaiki17再利用边界条件(c)得 iiiiBxkA21sin0ykxkCykBxkAyxUiiiiiiiiishsinshsin),(11若对任意x成立,则 02iBiiiBAC11其中 ,若用整数n代替i,则上式表示为 yanxanCyxUnnshsin),(1把边界条件(d)代入上式,得)0(shsin10axbanxanCUnn
10、将U0 在(0,a)区间内展开为 的傅立叶级数sinnxa于是0sinnnUfxan=100041,3,5.2sin02,4,6.anUnn xfUdxnaan其中18),6,4,2(0),5,3,1(sh140nnbannUCn对比两边各项系数可得xanbanyannUyxUnsinshsh4),(05,3,1因此,在槽内的电位分布是 0UU 0U 0U 0U 19 y 图53 缝隙导体板间的电位 0UU 0U x b 2b o 例5.2 有两块一端弯成直角的导体板相对放置,中间留有一小缝,如图所示。设导体板在x轴和z轴方向的长度远远大于两导体板间的距离b,上导体的电位为U0,下导体接地。求
11、两板间的电位分布。解:电位分布与z无关,这是一个二维拉普拉斯问题 00(0)()yUxa 0(0)()y bUUxb 0),(),(2222yyxUxyxU边界条件20002()002()()xbUybUcby0(0)()xUUyybdb边界条件的分解 0UUbya ybUUxa00 ybUUUxb000 ybUUxb00 0 xbU(A)(B)图54 图53的场分解 y y x x 00yaU 0bybU 00ybU 21对于(A)问题,情况与平行板电容器相同,两极间为匀强电场ybUyxUa0),()sin()(bynBygnnexp)(bxnAxfnn对于(B)问题沿y轴方向存在两个零电位
12、边界,取 沿x轴正方向边界无界且电位趋于零,取 因此(B)问题的电位通解应为 ybneCyxUxbnnnbsin),(1将边界条件(c)代入,得)()(220sin0001bybybUUbyybUybnCnn22ybnenUyxUxbnnnb2sin)1(),(210利用傅立叶展开并比较系数,可求得(B)问题的解利用叠加原理所以原问题的电位表达式为 ybnenUybUyxUxbnnn2sin)1(),(2010原问题的解应该是A 问题和B 问题解的线性叠加23二.圆柱坐标系中的分离变量法1、方法介绍变量的分离圆柱坐标中拉普拉斯方程为22222110rrrrrz01122222)(zhhrggr
13、frrfr)(2fghr)()()(),(zhgrfzr令 ,代入上式,并在两边同乘以得得上式第二项只与 有关,可以先分离出来,令其等于常数 2n0222gndgd将上式代回到式(A),各项同除以 r2,得(A)240112222)(zhhrnrfrrfr令上式左边最后一项等于常数 ,则上式分离成为两个方程2zk0222hkzdhdz01)()(222frnkrdfdrrddrz解的选取02nnAnAgcossin)(21 时,对于()g x对于()h xzkzkzzeBeB2102zkzBBzh21)(时,022zzkkzkBzkBzhcossin)(21 时,02n21)(AAg 时,02
14、zkzkBzkBzhzzchsh)(21 时,250)(22222fnrkrdfdrrdfdrz对于 ,是一个n阶贝塞尔方程()f r0,022zknnnrCrCrf21)(时 022zknrCCrfln)(21 时 0,022zkn)(N)(J)(21rkCrkCrfznzn 时 0,0222zzkkn)(K)(I)(21rkCrkCrfznzn 时 )(Jrkzn称为n阶第一类贝塞尔函数)(Nrkzn称为n阶第二类贝塞尔函数,也叫聂曼函数 第一类和第二类变形贝塞尔函数()()nznzIk rKk r其中本征解的叠加构成通解)()()(),(),(zhgrfzrzriiiiii26例5.3
15、电场强度为 的均匀静电场中放入一半径为a的电介质长圆棒,棒的轴线与电场相垂直,棒的电容率为 ,外部电容率为 ,求任意点的电位。20E1 P a r y 2 图55 均匀电场中的电介质棒 x 1 2 1解:在柱坐标系中求解 根据边界条件,写出通解(,)U rz与z无关,取 ,则()1h z 20zk 故(,)(,)U rzU r对于 ,根据问题的形式,可知)(g 2gggg所以取 n2 0,并且 n为正整数)(gncos只取 项27对于 ,因为20,02zn k()f r所以 取()f r12nnC rC r于是通解可以表示为)(cos),(211nnnnnrCrCnrUnrCrUnnncos)
16、,(1112021(,)coscosnnnUrE rC rn 求分区通解r a 时,总电位包括外电场和介质棒两部分的贡献 cos00rExE外电场电位介质棒电位)(cos),(211nnnnnrCrCnrU规定了200 xU考虑到当 时,介质棒产生的电位是有限值,故r 20nC28利用静电场的边界条件求系数ararrUrUaUaU221121,),(),(将U1和U2的通解表达式代入上面两个边界条件,得 naCaEnaCnnnnnncoscoscos12011nanCEnanCnnnnnncoscoscos1122021111比较上面两式两边 的系数,得 cos111021211120221C
17、 aE aC aCEC a n=1时02212121021211;2EaCEC以上两式联立求解,得 29n 1时12111122nnnnnnnnC aC anC anC a 120,0nnCC联立解上面两式求解,得)(2cos2),(021202121arxErErUcoscos),(20212102raErErU求出圆柱内外电位求出圆柱内外电场强度)(2sin2cos202120212021211arExEErUE)(sin1cos10222121022212122)()(arEraErarUE212012cos()aErr()ra30 1mU2mU中的无限长圆柱体是和上面例子完全类似的磁场
18、在均匀外磁场 边值问题。用0H和 分别表示圆柱体内外的磁标位,则它们的边界 条件也和电位的形式相同,即 r 时,cos02rHUmar 时,rUrUUUmmmm221121;此时,1mU和2mU的解也与静电场1U和2U的解有相同的形式,只需把1U和2U解中的1和 2分别用1和2代替,0E用0H代替,便得到1mU和2mU的解。同样可以得知圆柱体内的磁场强度1H是均匀的。31 a x b o z y 图56 导体圆筒内的电位 0000arbzzUUUU 例 5.4 一导体圆筒的高度为b,半径为a,所给边界条件为 求圆筒内的电位分布函数。解:由问题对称性可知,电位与变量无关,因此应取0n和1)(g;
19、)(zh应该取成在零点处为零值)(shzkz02zk;并考虑到电位在r=0处为有限值,故)(rf的形式为)(J0rkz。因此电位函数的本征解为:)(J)(sh),(0rkzkAzrUzz的双曲正弦函数,有将边界条件0arU代入上式,得0)(J)(sh),(0akzkAzaUzz所以akz应是零阶贝塞尔函数的根ip032),3,2,1(iapkzizi 把所有本征值)(J)(sh),(0001rapzapAzrUiiii所对应的本征解相加,得到电位函数的通解将边界条件0)(UbzU代入上式,得)(J)(sh00010rapbapAUiiii上式两边同乘以)(J00arprj,从0a对r积分,得d
20、rraprrapbapAdrraprUjiaiiija)()()()(00000010000JJshJ 对上式两边分别应用贝塞尔函数积分公式 aaxxdxxx0100)(J)(J 和正交公式(540),可以得到33 )(J2sh)(J0212001020)(iiiiipabapAppaU)(Jsh201000)(iiiipbappUA 所以 101000000)(JshJsh2),()()()(iiiiiipbapprapzapUzrU 例5.5 将例题5.4的边界条件改成0000UUUUarbzz求圆筒内的电位分布函数。解:由问题对称性可知,电位与变量无关,因此应取0n和1)(g)(zh应该
展开阅读全文