第3章〓时变电磁场 第2章讨论了静态电磁场的基本规律,归纳如下: 静电场的高斯定理: SD·ds=Vρdv或·D=ρ 静电场的环路定律: ∮LE·dl=0或×E=0 静磁场的磁通连续性原理: SB·ds=0或·B=0 静磁场的安培环路定律: ∮LH·dl=SJ·ds或×H=J 反映媒质特性的结构方程: D=εE,B=μH,J=σE 实验和研究证明,以上方程只能反映静态电磁场的规律和性质,对于随时间变化的电磁场则存在缺陷,应加以修正和推广,而且上述方程没有体现出电场与磁场的相互关系。英国物理学家法拉第1831年发现了法拉第电磁感应定律,揭示了“变磁生电”的规律,开启了对电场、磁场相互作用的研究。之后,英国物理学家麦克斯韦总结前人的研究结论,加以自己创造性的思考,提出了位移电流的概念,证实了“变电生磁”的可能性,麦克斯韦还提出了涡旋电场的概念,使法拉第电磁感应定律更具普适性。麦克斯韦根据前人和自己的研究成果,总结出了全面、准确阐述电磁理论的初始麦克斯韦方程组,这组方程的确立标志着经典电磁理论的建立,它们揭示了电场与磁场的相互关系、变化电场和变化磁场的规律和性质。 本章主要内容是时变电磁场理论,先介绍麦克斯韦方程组,再以此为基础阐述时变电磁场的规律和性质,包括边界条件、辅助位函数、能量与能流密度、波动方程和波动性,最后研究一种具有代表性的时变电磁场——时谐电磁场。 3.1麦克斯韦方程组 3.1.1法拉第电磁感应定律与感应电场 法拉第电磁感应定律表述为: 若穿过闭合导体回路L所张曲面的磁通量Ψ=SB·ds发生变化,回路中会出现感应电动势ε,如图31(a)所示,且感应电动势与磁通量变化率的关系为 ε=-dΨdt=ddtSB·ds(31) 感应电动势ε引起的感应电流I力图阻止回路中磁通的变化。 法拉第电磁感应定律说明变化的磁场与电流之间存在密切关系,将磁现象与电现象联系起来。不过,法拉第电磁感应定律没有说明出现感应电动势的真正原因,以及当时变磁场附近不存在导体回路时会发生什么情况。 图31法拉第电磁感应定律 麦克斯韦在对电磁感应现象进行深入分析后,认识到导体中的电流必然由电场引起。他将这种由变化的磁场激励或者说感应出来的电场称为感应电场,记为Ein。Ein在闭合导体回路上的环量等于回路上的感应电动势ε。不论有、无导体回路,时变磁场都会激励起感应电场,如图31(b)所示。因此,电磁感应现象的实质是时变磁场在周围空间激励起感应电场,如果该电场中有导体存在,就会在导体上引起感应电流。 根据以上分析,考虑仅有磁场随时间变化、闭合路径及其所张曲面均不随时间变化的情况,式(31)可写成 ∮LEin·dl=-dΨdt=-ddtSB·ds=-SBt·ds(32) 式中: L为任意闭合路径,不一定是导体回路; S为L所张的任意曲面。 式(32)是法拉第电磁感应定律的积分形式。 依据旋度定理,有 ∮LEin·dl=S(×Ein)·ds 将上式代入式(32),可得 S(×Ein)·ds=-SBt·ds 要使上式对任意曲面S都成立,必定有 ×Ein=-Bt(33) 这是法拉第电磁感应定律的微分形式。 式(33)说明,时变磁场是感应电场的旋涡源,感应电场是有旋场,也可称为涡旋电场。感应电场的电力线是闭合曲线,与其旋涡源——时变磁场的磁力线相交链。感应电场与静电场的激励源不同,性质也不同,但二者对电荷的作用力相同。 【例31】已知无源空间中感应电场强度E=10-2sin(6.28×109t-20.9z)y^(V/m),求磁感应强度。 解: Bt=-×E=Eyzx^=-20.9×10-2cos(6.28×109t-20.9z)x^ 认为静磁场B0为零,则有 B=∫Btdt=-3.33×10-11sin(6.28×109t-20.9z)x^(T) 3.1.2位移电流与全电流定律 第2章中得出的有关静态电磁场的结论和方程只适用于静态电磁场,在电磁场随时间变化的情况下,这些结论和方程会出现缺陷和矛盾,不再适用。 以静磁场的安培环路定律×H=J为例,对其两边取散度,再依据式(145g)可得 ·J=·(×H)≡0(34) 然而,依据电流连续性方程有·J=-ρ/t,当电荷密度不恒定时有ρ/t≠0,因此·J≠0,与式(34)矛盾。这说明对于电荷密度随时间变化的情况,静磁场的安培环路定律不再适用。也就是说,静磁场的安培环路定律不再适用于时变场,必须加以修正。 考虑到电流连续性方程是依据电荷守恒定律推导出来的,本身不会有错误,因此仍然以电流连续性方程为依据来考虑如何修正静磁场的安培环路定律。电流连续性方程可以写成·J+ρ/t≡0的形式,将电场的高斯定理ρ=·D(实践和研究证明电场的高斯定理适用于静态场和时变场)代入该式可得 ·J+(·D)t=·J+Dt≡0(35) 比较式(34)和式(35)可知,若在安培环路定律×H=J的右边加上一项D/t作为修正的安培环路定律,即 ×H=J+Dt(36) 若对上式等号两边均取散度,则有 0≡·(×H)=·J+Dt≡0 依据场论中的恒等式,左边的恒等号始终成立; 依据电流连续性方程和式(35),右边的恒等号始终成立,不会像静磁场的安培环路定律那样出现矛盾。可见,增加D/t项后,安培环路定律得以修正,适用于时变电磁场的情况,也适用于静态场的情况(此时D/t=0,就是静态场的安培环路定律)。 式(36)中,D/t项与电流密度J等价,麦克斯韦认为它是由变化的电场引起的另一种电流密度,称为位移电流密度Jd,其对应的电流称为位移电流id,它与普通电流一样也是磁场的旋涡源。这说明变化的电场也能激励磁场,这就是“变电生磁”的基本规律。应当注意,位移电流与普通电流不同,它不是实际带电粒子定向运动形成的,而是变化电场形成的一种等效电流,其激励磁场的作用与普通电流等效。位移电流是麦克斯韦引入的一个重要概念,它有助于人们理解“变化电场激励磁场”的原理和规律。 普通电流和位移电流一起称为全电流,J+D/t称为全电流密度矢量。由式(35)可知,全电流密度矢量的散度恒为零,其场是无散场,其矢量线为闭合曲线。这说明在带电粒子形成的电流中断之处必然有位移电流连续,形成闭合的电流线。这是因为在带电粒子形成的电流中断之处电量必然发生变化,会产生时变电场,也就出现了位移电流。 引入位移电流之后的修正的安培环路定律称为全电流定律,式(36)是其微分形式。依据旋度定理可推出全电流定律的积分形式为 ∮LH·dl=SJ+Dt·ds(37) 【例32】在无源空间中,磁场强度H=H0sin(ωt-kz)y^(A/m),k为常数。求位移电流密度和电位移矢量。 解: 由于空间中无源,因此传导电流密度J=0,由全电流定律求出位移电流密度,即 Jd=×H=-Hyzx^=kH0cos(ωt-kz)x^(A/m2) 将位移电流密度Jd对时间t积分,求出电位移矢量D(认为静态不变的D0等于零),即 D=∫Dtdt=∫Jddt=kωH0sin(ωt-kz)x^(C/m2) 3.1.3麦克斯韦方程组的形式及意义 法拉第电磁感应定律和全电流定律适用于静态场、时变场,它们分别给出了电场、磁场的旋度,即场与旋涡源的关系。但对矢量场的研究还必须了解它的散度,即场与散度源的关系。第2章推导出了静电场、静磁场与各自散度源的关系,即高斯定理·D=ρ,磁通连续性原理·B=0,在时变场情况下至今还没有发现与它们相矛盾的事实,因此它们也适用于时变场。 高斯定理、磁通连续性原理与法拉第电磁感应定律、全电流定律一起构成宏观电磁场的基本定律,其数学形式就是电磁场的场方程,称为麦克斯韦方程组。其微分形式为 全电流定律: ×H(r,t)=J(r,t)+D(r,t)t(38a) 法拉第电磁感应定律: ×E(r,t)=-B(r,t)t(38b) 磁通连续性定律: ·B(r,t)=0(38c) 高斯定理: ·D(r,t)=ρ(r,t)(38d) 积分形式的麦克斯韦方程组为 全电流定律: ∮LH(r,t)·dl=SJ(r,t)+D(r,t)t·ds(39a) 法拉第电磁感应定律: ∮LE(r,t)·dl=S-B(r,t)t·ds(39b) 磁通连续性定律: SB(r,t)·ds=0(39c) 高斯定理: SD(r,t)·ds=Vρ(r,t)dv(39d) 这些场方程同样也适用于静态电磁场的情况,此时方程中所有对时间求导的项均等于零,方程组(38)和(39)变成与静态电磁场方程组相同的形式。 麦克斯韦方程组揭示了电磁场的各个场矢量与场源的关系,只要知道了这些关系就可以从场源求出电磁场分布。同时,依据场论知识也可以由麦克斯韦方程组得知电场、磁场的性质。由麦克斯韦方程组可以得到如下结论: (1) 电场的散度等于电荷密度,电荷是电场的散度源; 由电荷产生的电场是有散场,电力线起始于正电荷,终止于负电荷。 (2) 磁场的散度恒为零,磁场没有散度源,科学界至今没有证实磁荷(即磁单极子)的稳定存在; 磁场是无散场,磁力线无头无尾,是闭合曲线。 (3) 时变磁场是感应电场的旋涡源,即时变磁场可以激励出感应电场; 感应电场是有旋场,其电力线是闭合曲线,与磁力线相交链。 (4) 全电流是磁场的旋涡源,电流和时变电场都可以激励出磁场; 磁场是有旋场,磁力线是闭合曲线,与全电流线相交链。 (5) 时变电场、时变磁场在某些情况下可以不断互相激励,且电力线与磁力线相互交链。 (6) 在某些情况下,场源(电荷或电流)一旦激励起了时变电场或时变磁场,即使去掉场源,时变电场、时变磁场也会互相激励,并向周围的空间传播,如图32所示,这就是电磁场的传播。 图32相互交链的闭合电力线、闭合磁力线 (7) 在线性媒质中,麦克斯韦方程组是线性方程组,满足叠加原理,即多个场源各自产生的场可以在空间同时存在,空间任一点的场等于所有场源在该点产生的场的叠加。 以上几点说明了电磁场与源的关系、电磁场的性质以及传播机理,也就是麦克斯韦方程组的物理意义。 麦克斯韦方程组含有5个矢量E、D、H、B、J和一个标量ρ,只由4个麦克斯韦方程无法确定地求出所有未知量,因此还需要辅助方程。在第2章中得出了描述媒质中场矢量E、D、H、B、J之间相互关系的结构方程,即 D=εE(310a) B=μH(310b) J=σE(310c) 推出这些方程时并没有规定场是静态的还是时变的,因此它们也适用于时变电磁场。这三个结构方程是麦克斯韦方程组的辅助方程。 结构方程中系数ε、μ、σ分别是媒质的介电常数、磁导率和电导率,称为媒质的电磁参数。真空中,ε0≈8.854×10-12F/m,μ0=4π×10-7H/m,σ=0。其他媒质的ε、μ、σ参数由媒质的物质成分、微观结构决定。电导率σ=0的媒质称为理想介质,σ=∞的媒质称为理想导体,σ为非零有限值的媒质称为导电媒质。 电磁参数值与空间位置无关的媒质称为均匀媒质; 电磁参数不随电磁场矢量的模值变化的媒质称为线性媒质; 电磁参数不随电磁场矢量的方向变化的媒质称为各向同性媒质,否则称为各向异性媒质,其电磁参数是张量。本书中若非特别指明,讨论的一般是均匀、线性、各向同性媒质。本书中,将σ=0、ε和μ均为实常数的均匀、线性、各向同性媒质称为理想媒质,研究电磁问题时往往先考虑理想媒质中的简单情况,再考虑其他媒质中的情况。 【例33】在无源、均匀、线性、各向同性介质中,E=Acos(ωt-βz)x^(V/m),其中A、β为常数。若要这种场存在,β、ω之间应该满足何种关系?其他场矢量等于多少? 解: 任意一种电磁场存在的条件是它的所有场矢量必须同时满足麦克斯韦方程组。应用麦克斯韦方程组及结构方程,可以由E导出另3个场矢量H、D、B及β、ω之间的关系。 首先由结构方程得电位移矢量: D=εE=Aεcos(ωt-βz)x^(V/m) ·D=[Aεcos(ωt-βz)]x=0 在无源区域中ρ=0,上述D满足·D=ρ=0。 再由法拉第定律求出磁场B、H: ×E=-BtBt=-Exzy^=-Aβsin(ωt-βz)y^ B=∫Btdt=Aβωcos(ωt-βz)y^(T),满足·B=0 H=Bμ=Aβωμcos(ωt-βz)y^(A/m) 由无源情况下(J=0)的全电流定律,D、H应满足 ×H=Dt 将前面求得的D、H代入上式可得 -Aμωβ2sin(ωt-βz)x^=-Aεωsin(ωt-βz)x^ 比较等式两边,可得 β2=ω2με 因此,若β、ω满足这个关系,则以上求出的4个场矢量同时满足麦克斯韦方程组。 3.2边界条件 媒质边界指的是媒质电磁参数ε、μ、σ发生变化处的分界面。媒质电磁参数变化导致边界两侧的电磁场也发生变化,其变化规律也就是边界两侧电磁场之间的关系,称为电磁场的边界条件。 边界两侧的电磁场分布情况都满足麦克斯韦方程组,因此边界条件可由麦克斯韦方程组推导出来。在边界上电磁场发生突变时,只能用麦克斯韦积分方程组来推导边界条件。为推导方便起见,在边界上选取规整的闭曲面或闭曲线,并在这些闭曲面或闭曲线上应用麦克斯韦积分方程组,就可推导出电磁场的边界条件。下面就依照这个思路推导电磁场的一般边界条件,然后再讨论一些特殊情况下的边界条件。 3.2.1边界条件的一般形式 1. 电场的边界条件 图33边界面上的矩形回路 在边界面上取一条规整的闭曲线,在其上应用法拉第电磁感应定律(式(39b)),可推导出电场强度矢量E的切向分量的边界条件。 如图33所示,环绕边界上任意一点O,取高度为Δh (其值趋于0)、宽度足够小的矩形回路L,矩形回路所张的曲面记为S。回路上、下底边ΔL2、ΔL1 分别处于媒质2和媒质1中并无限贴近边界,且与O点处的某个方向的切向单位矢量t^平行。记回路上、下底边的有向线段长度分别为ΔL2、ΔL1,根据几何关系,不妨令ΔL2=ΔL1≈ΔL(ΔL很小)。 在此矩形回路L上应用法拉第电磁感应定律,有 ∮LE·dl=∫ΔL1E·dl+∫ΔL2E·dl+∫侧边E·dl=S-Bt·ds 当矩形回路高度Δh趋于0时,在回路两侧边上的线积分∫侧边Edl趋于0; 又由于此时矩形回路所张曲面S的面积趋于0,且B/t为有限值,故B/t在曲面S上的面积分S(-B/t)·ds也趋于0。因此,上式可以写为 ∮LE·dl=∫ΔL1E·dl+∫ΔL2E·dl+∫侧边E·dl=∫ΔL1E·dl+∫ΔL2E·dl =S-Bt·ds=0 由于ΔL1、ΔL2足够小,可以认为ΔL1、ΔL2上的E分别为常矢量E1、E2,上式可写为 E2·t^ΔL2-E1·t^ΔL1=E2·t^ΔL-E1·t^ΔL=(E2-E1)·t^ΔL=0 即 (E2-E1)t^=0或E2t=E1t(311a) 可见,边界面两侧的E的切向分量相等,即边界面处E的切向分量连续。 设n^为边界面上O点处的法向单位矢量,τ^=n^×t^是与矩形回路平面垂直的另一切向单位矢量,则有t^=τ^×n^。应用矢量恒等式A·(B×C)=B·(C×A),式(311a)又可写为 (E2-E1)·(τ^×n^)=τ^·n^×(E2-E1)=0 因t^是边界面上该点处任意方向的切向单位矢量,故τ^的方向也是任意的。若要上式对任意方向的τ^都成立,则必有 n^×(E2-E1)=0(311b) 上式实际上是(311a)的矢量表示形式。 以上推出的式(311a)、式(311b)给出了边界面两侧E的切向分量的关系,也就是E的切向分量在边界面处的边界条件。 将D1t=ε1E1t,D2t=ε2E2t代入式(311a),可得 D2tε2=D1tε1(312) 可见,在边界面两侧D的切向分量不连续。 图34边界面上的圆柱面 下面再由高斯定律(式(39d))推出电位移矢量D的法向分量的边界条件。如图34所示,n^是边界面上任意点O处的法向单位矢量,由媒质1指向媒质2。环绕该点,跨边界两侧,取高度为Δh(其值趋于零)、横截面足够小的扁圆柱面S,其上、下底面ΔS2、ΔS1分别处于媒质2和媒质1中并无限贴近边界,且都垂直于n^。圆柱面在边界上截取的小曲面为ΔS,可认为ΔS1=ΔS2≈ΔS。 在此圆柱闭曲面S上应用高斯定律可得 SD·ds=ΔS1D·ds+ΔS2D·ds+侧面D·ds=Vρdv(313) 当Δh趋于0时,侧面面积也趋于零,因此侧面D·ds趋于0。且由于闭曲面S无限薄,其包围的体积V中的电荷Vρdv就等于边界上的小曲面ΔS上分布的自由面电荷。由于ΔS1、ΔS2、ΔS足够小,可认为ΔS1、ΔS2上的D分别为常矢量D1、D2,ΔS上的面电荷密度也等于常数ρs,于是式(313)可改写为 D2·n^ΔS2-D1·n^ΔS1=D2·n^ΔS-D1·n^ΔS=ρsΔS 即 n^·(D2-D1)=ρs(314a) D2n-D1n=ρs(314b) 式中: D1n、D2n分别为D1、D2的法向分量。上式就是D的法向分量的边界条件。可见,当边界面上分布有自由面电荷时,D的法向分量不连续,其突变量等于边界面上的自由面电荷密度ρs。 将D1n=ε1E1n,D2n=ε2E2n代入式(314b),可得 ε2E2n-ε1E1n=ρs,ρs≠0 0,ρs=0(315) 可见,无论边界面上有无自由面电荷,E的法向分量都不连续。 2. 磁场的边界条件 图35边界面上的矩形回路 由全电流定律(式(39a))可以推导出磁场强度H的切向分量边界条件。在如图35所示的边界面处的无限小矩形回路L上应用全电流定律,可得 ∮LH·dl=∫ΔL1H·dl+∫ΔL2H·dl+∫侧边H·dl =SJ·ds+SDt·ds(316) 式中: S为L所张的曲面,矩形回路在边界面上截取的曲线段记为ΔL。当Δh趋于0时,S的面积也趋于零,有S(D/t)·ds趋于零。而且,当Δh趋于0时,穿过S曲面的电流等于边界面上流过线段ΔL的面电流,即当媒质分界面上存在自由电流密度J时,有 limΔh→0ΔSJ·ds=∫Δl(limΔh→0ΔhJ)·τ^dl=∫ΔlJS·τ^dl 由于ΔL1、ΔL2、ΔL足够小,可认为ΔL1、ΔL2上的H分别为常矢量H1、H2,ΔL上的面电流密度也等于常矢量Js,再考虑到ΔL1=ΔL2≈ΔL,因此式(316)可改写为 H2·t^ΔL2-H1·t^ΔL1=H2·t^ΔL-H1·t^ΔL=(H2-H1)·t^ΔL=Js·τ^ΔL 即 (H2-H1)·t^=Js·τ^ 将t^=τ^×n^代入上式,可得 (H2-H1)·(τ^×n^)=Js·τ^ 应用矢量恒等式A·(B×C)=B·(C×A),上式可以改写为 [n^×(H2-H1)]·τ^=Js·τ^ 因t^是任意方向的,故τ^也是任意方向的,要使上式对任意方向的τ^均成立,必然有 n^×(H2-H1)=Js(317a) 即 H2t-H1t=Js(317b) 这就是磁场强度的切向分量满足的边界条件,说明当边界上分布有传导面电流时,H的切向分量不连续,突变量等于边界面上的传导面电流密度Js。 将B1t=μ1H1t,B2t=μ2H2t代入式(317b),得 图36边界面上的圆柱面 B2tμ2-B1tμ1=Js,Js≠0 0,Js=0(318) 可见,不管边界面上有没有传导面电流,B的切向分量都不连续。 下面从磁通连续性定律(式(39c))出发,推导磁感应强度B的法向分量边界条件。在图36所示边界处的小圆柱面S上应用磁通连续性定律,可得 SB·ds=ΔS1B·ds+ΔS2B·ds+侧面B·ds=0 与推导电场的法向分量的边界条件的过程类似,当Δh趋于零时,上式可以写为 n^·(B2-B1)=0(319a) 即 B2n=B1n(319b) 可见,边界面处B的法向分量连续。 将B1n=μ1H1n,B2n=μ2H2n代入式(319b),可得 μ2H2n=μ1H1n(320) 显然,边界面处H的法向分量不连续。 由4个麦克斯韦积分方程出发,经过前述推导最终得到电磁场的边界条件为 n^×(H2-H1)=Js n^×(E2-E1)=0 n^·(B2-B1)=0 n^·(D2-D1)=ρs或 H2t-H1t=Js E2t=E1t B2n=B1n D2n-D1n=ρs(321) 式中: Js为边界面上的传导面电流; ρs边界面上的自由面电荷。式(321)是电磁场的一般边界条件。 根据边界条件,可以由边界面一侧的电磁场,结合边界面上的自由电荷、传导电流分布情况,推出边界面另一侧的电磁场; 也可以推出边界面两侧电磁场满足的一些关系式,例如在分析电磁场传播入射到不同媒质交界面上的反射、折射问题时,就可以利用边界条件推出反射系数和折射系数,得到入射场、反射场、折射场的场矢量模值之间的关系式。 下面考虑一些特殊情况下的边界条件。 3.2.2理想导体表面的边界条件 设理想导体为媒质1,理想介质为媒质2。理想导体σ=∞,若理想导体中存在非零电场,必导致J=σE=∞,与“电流强度为有限值”的物理事实矛盾,因此理想导体中必定没有电场,也没有时变磁场(若存在时变磁场,必定会感应出电场),故理想导体中不存在时变电磁场,即E1、D1、H1、B1均为零。而理想介质中的场矢量E2、D2、H2、B2可以去掉下标“2”,记为E、D、H、B。因此,在理想导体与理想介质边界情况下,电磁场的边界条件为 n^×H=Js n^×E=0 n^·B=0 n^·D=ρs或 Ht=Js Et=0 Bn=0 Dn=ρs(322) 式中: n^为理想导体表面上的法向单位矢量,指向理想导体外侧。 以上边界条件说明: 理想导体表面上,电场切向分量、磁场法向分量均等于零,电场矢量、电力线必然垂直于理想导体表面,磁场矢量、磁力线必然平行于理想导体表面; 理想导体表面的自由电荷密度等于电位移矢量的法向分量,传导电流密度的模值等于磁场强度的切向分量。由此,可以根据理想导体表面处的电磁场求出理想导体表面的自由电荷、传导电流密度的分布情况。 【例34】如图37所示,在间距为d的两无限大导电平板之间充满空气,其中电场强度E=E0cos(ωt-βz)x^ (β为常数),求两导电平板表面上的面电荷密度和面电流密度。 图37例34图 解: 空气近似为介电常数、磁导率分别等于ε0、μ0的理想介质。上、下导电板内侧的外法线方向分别为n^d=-x^,n^0=x^。ρs、Js分别由D、H确定,故先由E求出D、H。 D=ε0E=ε0E0cos(ωt-βz)x^ ρs0=n^0·D=ε0E0cos(ωt-βz),x=0 ρsd=n^d·D=-ε0E0cos(ωt-βz),x=d 可见,上、下两导电平板内侧的面电荷密度等值异号,电力线从正电荷指向负电荷。 由×E=-Bt可知B=-∫(×E)dt,再由H=Bμ0可求出 H=E0ε0μ0cos(ωt-βz)y^ Js0=n^0×H=E0ε0μ0cos(ωt-βz)z^,x=0 Jsd=n^d×H=-E0ε0μ0cos(ωt-βz)z^,x=d 可见,上、下两导电平板内侧的面电流大小相等、方向相反。 3.2.3理想介质表面的边界条件 若边界面为两种理想介质的交界面,且交界面上没有人为增添的自由面电荷和传导面电流,此时Js=0,ρs=0,则电磁场边界条件为 n^×(H2-H1)=0 n^×(E2-E1)=0 n^·(B2-B1)=0 n^·(D2-D1)=0或 H2t=H1t E2t=E1t B2n=B1n D2n=D1n(323) 可见,在两种理想介质的交界面处,电场强度矢量、磁场强度矢量的切向分量都是连续的,电位移矢量和磁感应强度矢量的法向分量都是连续的。 图38例35图 【例35】如图38所示,已知y=0的无限大平面为两种电介质的分界面,该分界面上无自由电荷。介质2一侧的电场强度E2=x^10+y^20(V/m),分界面两侧的介电常数分别为ε1=5ε0,ε2=3ε0。求D1、D2和E1。 解: 先由E2求出D2,即 D2=ε2E2=ε0(x^30+y^60) 由题中条件可知,在两种介质分界面上,y^分量是法向分量,x^分量是切向分量。利用电位移矢量的边界条件可得 D1n=D2n=60ε0 由此可得 E1n=D1nε1=12 由电场的边界条件得 E1t=E2t=10 进而得到 D1t=ε1E1t=50ε0 所以 D1=ε0(x^50+y^60)(C/m2) E1=(x^10+y^12)(V/m) 3.3时变电磁场的位函数 麦克斯韦方程组描述了时变电场、时变磁场以及与产生它们的源之间的关系,方程中的电场和磁场是相互关联的,直接求解通常比较困难。为便于求解,可引入电磁场的位函数作为求电磁场场矢量的辅助函数,即先由场源求出位函数,再由位函数求出场矢量。下面先导出位函数满足的方程及位函数与场矢量的关系。 3.3.1标量位与矢量位 在均匀、线性、各向同性的不导电媒质中,麦克斯韦方程组为 ×H=J+Dt(324a) ×E=-Bt(324b) ·B=0(324c) ·D=ρ(324d) 因为·B=0,根据场论中“任何矢量的旋度的散度恒等于0,即·(×A)≡0”的结论,可将B写成一个矢量A的旋度的形式,即 B=×A(325) 将A定义为电磁场的矢量位。显然,这样由式(325)定义的B仍然满足·B=·(×A)=0。 把B=×A代入式(324b),并整理移项,可得 ×E+At=0 根据场论中“任何标量的梯度的旋度恒等于0,即×(u)≡0”,可将E+A/t写成一个标量Φ的负梯度的形式,即 E+At=-Φ(326) 将Φ定义为电磁场的标量位。显然E+A/t仍然满足 ×(E+A/t)=×(-Φ)=0 定义了电磁场的矢量位、标量位之后,场矢量B、E可以用它们来表示,即 B=×A(327a) E=-Φ-At(327b) 显然,只要求出A和Φ,就可以通过式(327)求出场矢量B和E。 3.3.2位函数满足的方程与解 通过麦克斯韦方程组(324)和位函数与场矢量的关系式(327)可以导出位函数满足的方程。在均匀、线性、各向同性的不导电媒质中,应用B=μH,D=εE以及ε、μ均为常数的已知条件,并将式(327)代入式(324a),可得 1μ××A=J-εtΦ+At 应用矢量恒等式××A=(·A)-2A将上式左边展开,并移项整理,可得 2A-εμA2t2=-μJ+·A+εμΦt(328) 为简化上述方程,注意到式(325)中只规定了A的旋度等于B,并没有对它的散度作出任何规定,这就意味着矢量位A 并不唯一。根据矢量场的性质,想要唯一确定某一矢量场,不仅要确定其旋度,还要确定其散度。于是,为使上式更简洁,可以设定A的散度为 ·A=-εμΦt(329) 此时,式(328)等号右边括号中的函数将等于0,方程式得到大大简化。在电磁场问题分析中,通常将这个设定称为洛伦兹条件,最早由丹麦物理学家Ludvig Lorenz于1867年提出。 将洛伦兹条件代入式(328),可得 2A-εμA2t2=-μJ(330a) 这就是均匀、线性、各向同性的不导电媒质中电磁场的矢量位A满足的微分方程。 为导出标量位Φ所满足的方程,将式(327b)代入式(324d),考虑到ε、μ均为常数,并交换算子和(·)/t算子的运算次序,可得 2Φ+t(·A)=-ρε 将洛伦兹条件代入上式,可得 2Φ-εμ2Φt2=-ρε(330b) 这就是均匀、线性、各向同性的不导电媒质中电磁场的标量位Φ满足的微分方程。 1. 标量位的解 首先研究无界、均匀、线性、各向同性、不导电媒质空间中标量位Φ(r,t)的解。 先求解点电荷的位函数,再应用叠加原理求出任意电荷分布产生的位函数。假设无界空间中仅在原点处有点电荷q(t),显然这种电荷分布具有球对称性。因此,空间的位函数分布也具有球对称性,Φ(r,t)仅是球坐标r和时间t的函数。由Φ(r,t)的方程(330b)可知,在除原点以外的无界空间中,在球坐标系下,Φ(r,t)满足的方程为 1r2rr2Φ(r,t)r=1v22Φ(r,t)t2 式中: v=1/εμ。 上式又可改写为 2[rΦ(r,t)]r2=1v22[rΦ(r,t)]t2 若以rΦ(r,t)为待求函数,由偏微分方程知识可知,在无界空间中,该偏微分方程的通解为达朗贝尔解,即 rΦ(r,t)=ft-rv+gt+rv(331) 上式右边第二项不符合实际物理情况,应当舍去,因此,位于原点的时变点电荷q(t)产生的位函数应为 Φ(r,t)=1rft-rv(332) 依据静态场知识,原点处电量为q的恒定点电荷在原点以外的无界空间产生的位函数为 Φ(r)=q4πεr(333) 静态场是时变场的特例,故式(333)是式(332)在q(t)等于常数q情况下的特例,比较这两式就可知函数f的具体形式。 于是,式(332)具体形式应为 Φ(r,t)=qt-rv4πεr(334) 这就是位于原点的时变点电荷q(t)在离自身距离为r(r≠0)的空间点产生的标量位。 对于式(334)可以理解为,位于坐标原点的时变点电荷在时刻t、位置r处产生的标量位由t-r/v时刻的电荷状态决定,即标量位是滞后于点电荷的。这是符合物理学因果定律的,电荷是“因”,标量位是“果”,有“因”才有“果”,“果”必然滞后于“因”。因此,式(334)又称为滞后位。 有了这一概念后,再来看式(331)的第二项为何不符合实际物理情况。根据上述类比过程,式(331)的第二项对应的解应表示为 Φ(r,t)=qt+rv4πεr 上式表明,位于坐标原点的时变点电荷在时刻t、位置r处产生的标量位由t+r/v时刻的电荷状态决定,即标量位是超前于点电荷的。这显然不符合物理学的因果定律,故而需要将这一种解舍去。 图39分布电荷示意图 如图39所示,若时变电荷分布在体积V中,则V中r′处的体积微元dv′中包含的电荷ρ(r′,t)dv′可以看作位于r′处的点电荷,它在空间任意r点(r≠r′)处产生的位函数为 dΦ(r,t)=ρr′,t-|r-r′|vdv′4πε|r-r′| 由于Φ(r,t)满足的方程(330b)是线性方程,故Φ(r,t)满足叠加原理。因此,体积V中所有电荷在r点产生的位函数Φ(r,t)应等于所有体积微元dv′产生的位函数dΦ(r,t)的叠加,即 Φ(r,t)=14πεVρr′,t-|r-r′|v|r-r′|dv′(335) 这就是无界、均匀、线性、各向同性、不导电媒质空间中标量位方程(330b)的解。 2. 矢量位的解 在直角坐标系中,矢量位的波动方程式(330a)可以分解成如下三个标量方程: 2Ai(r,t)-εμAi(r,t)t2=-μJi(r,t)(i=x,y,z)(336) 这三个标量方程与标量位Φ(r,t)的方程(330b)形式相同,故解的形式也相同,即 Ai(r,t)=μ4πVJir′,t-|r-r′|v|r-r′|dv′(i=x,y,z)(337) 写成矢量形式为 A(r,t)=μ4πVJr′,t-|r-r′|v|r-r′|dv′(338) 这就是无界、均匀、线性、各向同性、不导电媒质空间中矢量位的波动方程(330a)的解。 3.4电磁场的能量 在第2章中推出了静电场、静磁场的能量密度,它们分别为 we=12E·D, wm=12H·B 在得到这些能量密度公式时,对场的时变性质没有要求,因此它们也可用于计算时变场的能量。对于时变电磁场,由于时变电场和时变磁场可以相互激励、相互伴生的,时变电磁场的能量密度w是电场能量密度we与磁场能量密度wm之和,即 w=we+wm=12E·D+12H·B(339) 区域V中的总电磁能量为 W=Vwdv=V(we+wm)dv=V12E·D+12H·Bdv(340) 静态电磁场的场矢量不随时间变化,其能量是静态储存的。而对于时变电磁场,空间某一位置处的电磁能量总是在动态变化的,这就意味着不同空间位置之间必然存在着电磁能量的流动。实际上,时变电磁场能够从场源向周围空间传播,电磁场的能量也随之传播,形成电磁能流。任何物理现象及其工程应用都是通过能量传递来实现的,电磁场的应用也依靠电磁能量的传播才得以实现,比如,利用电磁场进行无线通信就是携带信息的电磁能量持续传播,使得远处的电磁场接收装置(包括接收天线和接收机)接收到这些信息。因此,电磁能量的大小和流动方向是研究电磁场时十分关心的特征。 下面推导如何用时变电磁场的场矢量表示电磁能量的大小和流动方向,即时变电磁场的能流。 图310有限无源区域内的 传导电流示意图 先假设在无源区域(媒质电导率为σ)中研究电磁场的能流。如图310所示,该区域内没有外加场源,其中的传导电流由电场E在媒质中引起,电流密度J=σE,它与电场、磁场的关系为J=×H-D/t。该传导电流引起的热损耗功率密度为 J·E=σE2=E·(×H)-E·Dt 将矢量恒等式E·(×H)=H·(×E)-·(E×H)代入上式右边,并根据×E=-B/t,最终有 -·(E×H)=J·E+H·Bt+E·Dt(341) 假设媒质是各向同性的线性媒质,上式的最后两项可改写为 H·Bt=12H·Bt+12B·Ht=t12H·B E·Dt=12E·Dt+12D·Et=t12E·D 式中: 12H·B、12E·D分别是磁场能量密度和电场能量密度。 为更好地理解式(341),依据其等号左边一项,引入一个新矢量——坡印廷矢量S,定义为 S=E×H(342) 因此,式(341)可以写为 -·S=J·E+t12H·B+12E·D 将上式两边在该无源区域中的任意体积V上做体积分,并依据散度定理 V·Sdv=SS·ds 可得 -SS·ds=VJ·Edv+tV12H·B+12E·Ddv(343) 式中: S为体积V的边界面。 式(343)右边第一项表示单位时间内V中电流损耗的焦耳热能量,第二项表示V中电磁能量随时间的增长率,或者说是单位时间内V中增加的电磁能量。由于已经假设研究区域中无外加场源,依据能量守恒定律,上式右边两项之和应等于“单位时间内通过V的边界面S流入V中的电磁能量”,即上式等号左边的项-SS·ds,如图311所示。式(343)描述的电磁能量守恒定律又称为坡印廷定理。 图311电磁能量守恒定律示意图 下面再讨论V内有外加电源的情况。此时V中的电流密度J可看作外加源电流Jg和电场引起的传导电流Jc=σE之和。将式(343)中的J替换成Jg+Jc,移项整理,可得 -SS·ds-VJg·Edv=VσE2dv+tV12H·B+12E·Ddv(344) 这就是V中有外加电源情况下的坡印廷定理。-VJg·Edv表示单位时间内V中的外加电源提供的能量。上式说明: “单位时间内通过边界面S流入V中的电磁能量”(用-SS·ds表示)与“V中外加电源提供的能量”之和,等于“单位时间内V中媒质电流损耗的焦耳热能量”和“单位时间内V中增加的电磁能量”之和,这同样符合能量守恒定律。 由无源、有源两种情况下的坡印廷定理式(343)、式(344)可知,“单位时间内通过V的边界面S流入V中的电磁能量”可用坡印廷矢量S在闭曲面S上的曲面积分-SS·ds来表示,积分微元S·ds表示单位时间内流过面积微元ds的电磁能量。因此,坡印廷矢量S=E×H描述了电磁能量的流动特性: S的方向表示电磁能量流动的方向,该方向垂直于电场强度和磁场强度,并依次构成右手螺旋法则; S的模值等于单位时间内通过与能流方向垂直的单位面积的电磁能量,其单位为瓦特/米2(W/m2)。坡印廷矢量S又称为电磁场的能流密度矢量。 研究电磁问题时,通过计算任意点处的坡印廷矢量S来得知该点处的电磁能流的方向和大小,通过计算坡印廷矢量S在某个曲面上的积分SS·ds来计算单位时间内流过该曲面的总电磁能量,即电磁功率。 【例36】在某无源理想介质区域中,E=Ecos(ωt-kz)x^(V/m),求坡印廷矢量以及坡印廷矢量的时间平均值。 解: 首先应用麦克斯韦方程组,由E求出H。 Bt=-×E=-Eksin(ωt-kz)y^ H=B μH=kωμEcos(ωt-kz)y^(A/m) S=E×H=kωμE2cos2(ωt-kz)z^(W/m2) S只有z分量,可见电磁能量沿z轴方向流动。 由上式可见,S是t的周期函数,其平均值也就是它在一个变化周期T中的平均值。 ωT=2π,T=2π/ω Sav=1T∫T0kωμE2cos2(ωt-kz)dt=k2ωμE2(W/m2) 3.5时变电磁场的波动性 前面在讨论麦克斯韦方程组时已经指出,一旦场源在空间激发起时变电磁场,时变电场、时变磁场就会持续地相互激发,即使去掉场源,时变电磁场也会脱离场源而存在,并以一定速度向远处传播。为了研究时变电磁场传播的具体形式和性质,应当推导出时变电磁场场矢量各自满足的方程。本节将从麦克斯韦方程组出发,推导出时变场场矢量各自满足的方程,这些方程具有波动方程的形式,称为电磁场的波动方程。 3.5.1波动方程 设所讨论的区域中有时变场源电流Jg和时变场源电荷ρ,媒质是均匀、线性、各向同性的,此区域中麦克斯韦方程组为 ×H=(Jg+σE)+εEt(345a) ×E=-μHt(345b) ·H=0(345c) ·E=ρε(345d) 对式(345b)两边取旋度,并交换旋度算子×与偏微分算子(·)/t的次序,可得 ××E=-μt(×H) 利用矢量恒等式××E=(·E)-2E以及式(345a)、式(345d),由上式可得 2E-σμEt2-με2Et2=μJgt+ρε(346a) 用类似的方法,从式(345a)出发又可推导出 2H-σμHt2-με2Ht2=-×Jg(346b) 式(346)具有波动方程的形式,通常将其称为电磁场的波动方程或达朗贝尔方程。 若研究的区域中无源,将Jg=0,ρ=0这两个条件代入式(346),得到无源情况下的电磁场波动方程为 2E-σμEt-με2Et2=0(347a) 2H-σμHt-με2Ht2=0(347b) 若无源区域中的媒质是不导电媒质,其电导率σ=0。将σ=0代入式(347),得到不导电媒质中的无源波动方程为 2E-με2Et2=0(348a) 2H-με2Ht2=0(348b) 依据有源波动方程(式(346)),可以由场源求出辐射的电磁波。依据无源波动方程(式(347)或式(348)),可以研究电磁波脱离场源之后在无源空间中传播的性质。 3.5.2波动性 为了说明电磁场的波动性,先在简单情况下求解波动方程,再分析波动性。 设媒质是均匀、线性、各向同性、不导电的,所研究区域中无源,此时场矢量满足的波动方程为方程(348)。令ψ表示E(r,t)或H(r,t)的任意直角坐标分量。为简化求解,假设ψ仅是坐标z和时间t的函数,即ψ=ψ(z,t),将其代入式(348),得ψ所满足的方程为 2ψ(z,t)2z-με2ψ(z,t)t2=0(349) 由偏微分方程知识可知,在无界空间中,该偏微分方程的通解为达朗贝尔解,表示为 ψ(z,t)=ft-zv+gt+zv(350) 式中: v=1/εμ; f、g为二阶可微函数。 下面分析达朗贝尔解的物理意义。第一项f(t-z/v)有空间坐标z 和时间坐标t两个自变量。假设在z=z1的P1点,函数f随t变化的曲线如图312(a)所示,f在t1、t2、t3时刻的函数值分别为f1=f(t1-z1/v),f2=f(t2-z1/v),f3=f(t3-z1/v)。 图312不同的空间点处函数的时变形式 在z=z2=z1+Δz(z2>z1)的P2点,同样的函数值f1必然出现在t′1=t1+Δt=t1+Δz/v时刻,因为f[t′1-(z1+Δz)/v]=f(t1-z1/v)=f1,所以可以认为P1点在t1时刻的函数值经过一段长为Δt的时间传到了P2点,传播的速度v=Δz/Δt。同理,P1点在其他时刻的函数值也会经过一段长为Δt的时间传到P2点。比如,P1点在t2时刻的函数值为f2,则P2点在t′2=t2+Δt时刻的函数值也等于f2,因为f[t′2-(z1+Δz)/v]=f[t2+Δt-(z1+Δz)/v]=f(t2-z1/v)=f2; P1点在t3时刻的函数值为f3,则P2点在t′3=t3+Δt时刻的函数值也等于f3; …。因此,P2点总是在重复P1点出现过的函数值,只不过时间滞后了Δt。图312(b)中绘出P2点处函数f随时间t的变化曲线,它与P1点的时变曲线完全相同,这说明函数f在P1点的变化(即振动)经过一段时间就传播到了P2点,传播速度为v,传播方向为z值增大的方向。“振动在空间中的传播”被定义为“波动”,由此可知函数f表示一个向正z方向传播的波动。 用同样的方法分析式(350)的第二项g(t+z/v),可知,它表示向z值减小方向,也就是向负z方向传播的波动,传播速度v=1/εμ。因此,波动方程(349)的通解ψ包含两个向相反方向传播的波动。 前面已经假设,ψ表示场矢量E(r,t)或H(r,t)的任意直角坐标分量,同理E(r,t)或H(r,t)的其他直角坐标分量也和ψ一样具有波动性,因此E(r,t)、H(r,t)都具有波动性,这也是将其称为电磁波的原因。电磁波的传播速度v=1/εμ,在空气(其ε、μ值近似为ε0、μ0)中,v0=1/ε0μ0≈2.9979×108(m/s)。 需要说明的是,上述分析过程中并没有给出ψ中函数f(t-z/v)、g(t+z/v)的具体表示式。实际上,函数f(t-z/v)、g(t+z/v)的具体表示式取决于电磁波辐射源的空间分布形式及其时间变化形式,可通过有源波动方程(346)求解。但不管f(t-z/v)、g(t+z/v)的具体表示式如何,都不影响上述波动性分析的结论。 3.6时谐电磁场 通过分析时变电磁场的波动性可知,时变电磁场场矢量的时变形式取决于场源ρ、J的时变形式。理论上来说,时变形式有无穷多种,但在研究电磁学和应用电磁场时往往特别关注时谐电磁场,即场分量随时间做简谐变化(包括正弦变化和余弦变化)的电磁场。时谐电磁场也是一种实际中广泛应用的电磁场,任意时变形式的电磁场理论上都可以用傅里叶级数或傅里叶积分展开成时谐电磁场的叠加。因此,研究时谐电磁场得到的结论可以推广到其他时变形式的电磁场,时谐电磁场具有典型性、代表性,研究时谐电磁场是研究其他时变形式电磁场的基础。 3.6.1时谐电磁场的瞬时表示式和复数表示式 时谐电磁场由时谐场源产生。时谐电磁场的场矢量也随时间t以角频率ω做简谐变化,各个分量也可用时谐函数表示,以电场强度E(r,t)为例: E(r,t)=Ex(r,t)x^+Ey(r,t)y^+Ez(r,t)z^ =Exm(r)cos[ωt+x(r)]x^+Eym(r)cos[ωt+y(r)]y^+ Ezm(r)cos[ωt+z(r)]z^(351) 式中: Exm(r)、Eym(r)、Ezm(r)分别为Ex、Ey、Ez分量在r点处的振幅; x(r)、y(r)、z(r)分别为Ex、Ey、Ez分量在r点处的初始相位。 其他场矢量也可以表示成类似的形式,这种形式的表示式称为时谐电磁场的瞬时表示式。 实际上,时谐场瞬时表示式中与时间有关的角频率ω和时谐函数形式cos(ωt+…)都可以作为已知信息而隐匿,只需要关注与空间坐标r有关的信息,例如振幅Exm(r)、Eym(r)、Ezm(r)和初始相位x(r)、y(r)、z(r)。按照这样的思路,可将时谐场的场矢量、场源都用复数表示,并将场方程写成复数形式,使得时谐场的表示式和计算、分析过程大为简化。 以电场强度E(r,t)为例,其x分量可以写成如下复数形式: Ex(r,t)=Exm(r)cos[ωt+x(r)]=Re[Exm(r)ejx(r)ejωt]=Re[Ex(r)ejωt] (352) 式中: 复数Ex(r)=Exm(r)ejx(r)称为Ex分量在r处的复振幅,它与时间无关,但包含了Ex分量与空间位置r有关的所有信息,即振幅Exm(r)和初始相角x(r)。依据Ex(r)和角频率ω,就可以通过式(352)还原出Ex(r,t),因此用复振幅Ex(r)来代替Ex(r,t)既简洁又完备。 类似地,电场强度E(r,t)的其他分量也可以写成用其复振幅表示的形式: E(r,t)=Exm(r)cos[ωt+x(r)]x^+Eym(r)cos[ωt+y(r)]y^+Ezm(r)cos[ωt+z(r)]z^ =Re[Exm(r)ejx(r)ejωt]x^+Re[Eym(r)ejy(r)ejωt]y^+Re[Ezm(r)ejz(r)ejωt]z^ =Re[Ex(r)ejωt]x^+Re[Ey(r)ejωt]y^+Re[Ez(r)ejωt]z^ =Re{[Ex(r)x^+Ey(r)y^+Ez(r)z^]ejωt} =Re[E(r)ejωt](353a) 式中: Ey(r)=Eym(r)ejy(r)、Ez(r)=Ezm(r)ejz(r)分别是Ey分量、Ez分量的复振幅。而矢量E(r)=Ex(r)x^+Ey(r)y^+Ez(r)z^称为E(r,t)在r处的复振幅矢量,又称为E(r,t)的复数表示式。复数表示式比瞬时表示式更简洁,因此在时谐场条件下往往采用复数表示式来表示场矢量。 场矢量D(r,t)、H(r,t)、B(r,t)、J(r,t)都可按上述方法用其复数形式来表示: D(r,t)=Re[D(r)ejωt], D(r)=Dx(r)x^+Dy(r)y^+Dz(r)z^(353b) H(r,t)=Re[H(r)ejωt], H(r)=Hx(r)x^+Hy(r)y^+Hz(r)z^(353c) B(r,t)=Re[B(r)ejωt], B(r)=Bx(r)x^+By(r)y^+Bz(r)z^(353d) J(r,t)=Re[J(r)ejωt], J(r)=Jx(r)x^+Jy(r)y^+Jz(r)z^(353e) 需要注意以下三个方面: (1) 只有时谐电磁场的场矢量才能写成复数表示式,复数表示式只与空间位置变量有关,与时间变量无关。 (2) 复数表示式仅仅是时谐电磁场的一种简化数学表示形式,真实的电磁场矢量是与之相对应的瞬时表示式。 (3) 只有频率相同的时谐场之间才能直接使用复数表示式进行运算。 分析、研究时为方便起见,往往需要将场矢量的瞬时形式与复数形式相互转换,以下例题说明了转换过程和需注意的事项。 【例37】写出与下列场矢量表示式对应的瞬时表示式或复数表示式: (1) E(r)=E0sin(kxx)sin(kyy)e-jkzzz^ (2) E=2jE0sinθcos(kxcosθ)e-jkzzsinθx^ (3) E=Eymcos(ωt-kx+α)y^+Ezmsin(ωt-kx+α)z^ (4) H=H0sin(kz-ωt)x^+H0cos(kz-ωt)z^ 解: (1) 原表示式是复数表示式,要写出它的瞬时表示式。利用式(353a)可得 E(r,t)=Re[E(r)ejωt] =Re[E0sin(kxx)sin(kyy)e-jkzzejωtz^] =E0sin(kxx)sin(kyy)cos(ωt-kz)z^ (2) 原表示式是复数表示式,要写出它的瞬时表示式。注意复系数j要写成e的指数形式,然后将e的所有指数相加得到时谐函数的相位,于是有 E(r,t)=Re[E(r)ejωt] =Re[2jE0sinθcos(kxcosθ)e-jkzzsinθejωtx^] =Re[2E0sinθcos(kxcosθ)ejπ/2e-jkzzsinθejωtx^] =2E0sinθcos(kxcosθ)cos(ωt-kzzsinθ+π/2)x^ =-2E0sinθcos(kxcosθ)sin(ωt-kzzsinθ)x^ (3) 原表示式是瞬时表示式,要写出它的复数表示式,须将其转化为式(353a)的形式。先将所有时谐函数都转化为cos函数,才能将它们写成复数的实部,于是有 E(r,t)=Eymcos(ωt-kx+α)y^+Ezmsin(ωt-kx+α)z^ =Eymcos(ωt-kx+α)y^+Ezmcos(ωt-kx+α-π/2)z^ =Re(Eymej(-kx+α)ejωty^+Ezmej(-kx+α)e-jπ/2ejωtz^) =Re[(Eymej(-kx+α)y^-jEzmej(-kx+α)z^)ejωt] =Re[E(r)ejωt] E(r)=Eyme-j(kx-α)y^-jEzme-j(kx-α)z^ (4) 原表示式是瞬时表示式,要写出它的复数表示式。将ωt前的负号化成正号,才能转化为式(353c)的形式,于是有 H(r,t)=H0sin(kz-ωt)x^+H0cos(kz-ωt)z^ =-H0cos(ωt-kz-π/2)x^+H0cos(ωt-kz)z^ =Re(-H0e-jkze-jπ/2ejωtx^+H0e-jkzejωtz^) =Re[(jH0e-jkzx^+H0e-jkzz^)ejωt] H(r)=jH0e-jkzx^+H0e-jkzz^ 3.6.2时谐场的复数形式麦克斯韦方程组和结构方程、边界条件 在时谐电磁场前提下,麦克斯韦方程组也可以写成复数形式。以全电流定律为例将其中所有的物理量都用复数表示式来表示,即 ×Re[H(r)ejωt]=Re[J(r)ejωt]+Re[D(r)ejωt]t(354) 交换取实部的Re运算与偏微分算子·、×的运算次序,可得 Re[×H(r)ejωt]=Re[J(r)ejωt+jωD(r)ejωt] 即 Re{[×H(r)-J(r)-jωD(r)]ejωt}=0(355) 上式对任意t都成立的条件为 ×H(r)-J(r)-jωD(r)=0 即 ×H(r)=J(r)+jωD(r) 上式就是时谐电磁场情况下全电流定律的复数形式。 类似地,在时谐电磁场情况下可将麦克斯韦方程组中其他三个方程写成复数形式,得到复数形式的麦克斯韦方程组: ×H(r)=J(r)+jωD(r)(356a) ×E(r)=-jωB(r)(356b) ·B(r)=0(356c) ·D(r)=ρ(r)(356d) 可以看到,与瞬时形式的麦克斯韦方程组相比,复数形式的麦克斯韦方程组唯一的变化就是对时间的偏微分项/t替换为jω,·、×等算子均保持不变。另外,需要特别注意,方程组(356)只适用于时谐电磁场,不适用于其他时变类型的电磁场。 时谐电磁场的结构方程和边界条件的复数形式分别为 D(r)=εE(r) B(r)=μH(r) J(r)=σE(r)(357) n^×[H2(r)-H1(r)]=Js(r) n^×[E2(r)-E1(r)]=0 n^·[B2(r)-B1(r)]=0 n^·[D2(r)-D1(r)]=ρs(r)或 H2t(r)-H1t(r)=Js(r) E2t(r)=E1t(r) B2n(r)=B1n(r) D2n(r)-D1n(r)=ρs(r)(358) 时谐电磁场的结构方程、边界条件的形式与一般结构方程、边界条件完全相同,只是其中的场矢量均写成复数表示式。 【例38】用场矢量的复数表示式重新求解例34中的问题。 解: 例34中给出的是时谐电磁场的电场强度瞬时表示式 E=E0cos(ωt-βz)x^=Re[(E0e-jβzx^)ejωt] 其复数表示式为 E(r)=E0e-jβzx^ 因此,有 D(r)=ε0E(r)=ε0E0e-jβzx^ ρs0(r)=n^0·D(r)=ε0E0e-jβz,x=0 ρsd(r)=n^d·D(r)=-ε0E0e-jβz,x=d 由×E(r)=-jωB(r),H(r)=B(r)/μ0可求出 H(r)=E0=ε0μ0e-jβzy^ Js0(r)=n^0×H(r)=E0ε0μ0e-jβzz^,x=0 Jsd(r)=n^d×H(r)=-E0ε0μ0e-jβzz^,x=d 由上述例题可知,对于时谐电磁场采用复数表示式和复数形式的麦克斯韦方程组可以使计算、分析得到简化。比如,原来利用法拉第电磁感应定律×E(r,t)=-B(r,t)/t由电场强度E(r,t)求磁感应强度B(r,t),要进行旋度运算和积分运算; 但在时谐场情况下,利用复数形式的法拉第电磁感应定律×E(r)=-jωB(r),由电场强度E(r)求磁感应强度B(r),就只需做旋度运算和除法运算了,比积分简单。而且,场矢量的表示式也更加简洁。 3.6.3时谐场的复坡印廷矢量和复坡印廷定理 1. 复坡印廷矢量 坡印廷矢量S(r,t)=E(r,t)×H(r,t)表示时变电磁场在时刻t、空间r点的能流密度,称为时变场的瞬时坡印廷矢量或瞬时能流密度。如果S(r,t)是周期变化的,则它在每个周期T内的时间平均值称为平均坡印廷矢量或平均能流密度 Sav(r)=1T∫T0S(r,t)dt 时谐电磁场情况下,场矢量的瞬时表示式E(r,t)、H(r,t)可以用其复数表示式写为 E(r,t)=Re[E(r)ejωt]=12[E(r)ejωt+E*(r)e-jωt] H(r,t)=Re[H(r)ejωt]=12[H(r)ejωt+H*(r)e-jωt] 式中: E*(r)、H*(r)分别为E(r)、H(r)的共轭。 因此,时谐电磁场的瞬时坡印廷矢量S(r,t)可以表示为 S(r,t)=E(r,t)×H(r,t) =12[E(r)ejωt+E*(r)e-jωt]×12[H(r)ejωt+H*(r)e-jωt] =Re12E(r)×H*(r)+Re12E(r)×H(r)ej2ωt 式中: 第一项与时间t无关,不随时间变化; 第二项以2ω为角频率随t做正弦变化。 图313中画出第一项、第二项及S(r,t)的时变曲线,可见S(r,t)的平均值为 Sav(r)=1T∫T0S(r,t)dt=Re12E(r)×H*(r)(359) 图313瞬时坡印廷矢量的时变曲线 式中Re(·)中的复矢量E(r)×H*(r)/2定义为时谐场的复坡印廷矢量或复能流密度矢量S(r),即 S(r)=12E(r)×H*(r)(360) 由式(359)可知,平均坡印廷矢量与复坡印廷矢量的关系为 Sav(r)=Re[S(r)](361) 应当注意的是,时谐场的复坡印廷矢量S(r)并不是其瞬时坡印廷矢量S(r,t)的复数表示式,即S(r,t)≠Re[S(r)ejωt],只能通过定义式(360)来求复坡印廷矢量S(r)。 2. 复坡印廷定理 时谐场情况下,坡印廷定理也可以写成复数形式。在体积V的边界曲面S上求复坡印廷矢量S(r)的面积分,并应用散度定理,可得 SS·ds=S12E×H*·ds=V·12E×H*dv(362) 将式 ·(E×H*)=H*·×E-E·×H* 及 ×E=-jωB,×H*=J*-jωD* 代入式(362),可得 -S12E×H*·ds=V12E·J*dv+j2ωV14B·H*-14E·D*dv (363a) 这就是时谐场满足的复坡印廷定理。 若媒质参数σ、ε、μ均为实数,则有 V12E·J*dv=V12σ|E|2dv=Pl V14B·H*dv=V14μ|H|2dv=Wm V14E·D*dv=V14ε|E|2dv=We 式中: Pl为体积V中的平均焦耳热损耗功率; Wm、We分别是体积V中的平均磁场能量和平均电场能量。 因此,复坡印廷定理式(362)可写为 -SS·ds=Pl+j2ω(Wm-We)(363b) 它表明,在体积V中无外加电源的情况下,通过边界面S流入体积V中的复功率-SS·ds的实部正好等于体积V中的平均焦耳热损耗功率,而流入的复功率的虚部等于V中平均磁场能量与平均电场能量之差的2ω倍,是流入V中的无功功率。 若V中存在外加电源,则将式(363a)中的J替换成Jg+Jc,此时复坡印廷定理表示为 -S12E×H*·ds-V12E·J*gdv =V12E·J*cdv+j2ωV14B·H*-14E·D*dv(364) 【例39】无源空间中电场强度E=Ecos(ωt-kz)y^(V/m),求复坡印廷矢量和平均坡印廷矢量。 解: 首先要将场矢量的瞬时表示式转化为复数表示式,并求出磁场强度。于是有 E(r,t)=Ecos(ωt-kz)y^(V/m)E(r)=Ee-jkzy^(V/m) ×E(r)=-jωμH(r)H(r)=-kEωμe-jkzx^(A/m) S(r)=12E(r)×H*(r)=kE22ωμz^(W/m2) Sav(r)=Re[S(r)]=kE22ωμz^(W/m2) 3.6.4时谐场的波动方程 时谐场的波动方程可以写成复数形式。将有源波动方程中的时谐场场矢量用其复数表示式表示,算子(·)/t、2(·)/t2分别用jω、(jω)2代替,可得有源区域中的复波动方程,即 2E(r)-jωσμE(r)+k2E(r)=jωμJ(r)+ρ(r)ε(365a) 2H(r)-jωσμH(r)+k2H(r)=-×J(r)(365b) 式中: k=ωεμ。 以上两个方程称为时谐场场矢量的亥姆霍兹方程。 将无源波动方程中的场矢量写成复数形式就可以得到无源、不导电媒质中的复数形式波动方程,即 2E(r)+k2E(r)=0(366a) 2H(r)+k2H(r)=0(366b) 3.6.5时谐场的位函数 1. 位函数方程 在3.3节介绍了时变场的位函数,由式(330)可知,标量位函数和矢量位函数分别由场源ρ(r,t)、J(r,t)决定。时谐电磁场的场源ρ(r,t)、J(r,t)随时间做时谐变化,因此其位函数也是时谐函数,也可以写成复数表示式。 可以由一般时变场的位函数方程直接推出时谐场的位函数满足的方程。将方程(330)中的时谐场场矢量用其复数表示式表示,2(·)/t2用(jω)2代替,即得时谐场的位函数方程 2A(r)+k2A(r)=-μJ(r)(367a) 2Φ(r)+k2Φ(r)=-ρ(r)ε(367b) 式中: k=ωεμ。 2. 位函数 求解方程(367)就可以求出时谐场的复数形式位函数,更简便的方法是直接由一般时变场的位函数推出时谐场的复数形式位函数。 先考虑标量位函数Φ。设ρ(r,t)为时谐函数,表示为 ρ(r,t)=ρm(r)cos(ωt+)=Re[ρ(r)ejωt](368) 式中: ρ(r)=ρm(r)ej是ρ(r,t)的复数表示式。 与此同理,一般时谐场的位函数式(335)中积分号内的ρ(r′,t-|r-r′|/v)可在时谐场情况下表示为 ρr′,t-|r-r′|v=ρm(r′)cosωt-|r-r′|v+ =Reρm(r′)ejejωt-|r-r′|v =Re[ρ(r′)e-jk|r-r′|ejωt](369) 式中 ρ(r′)=ρm(r′)ej,k=ω/v=ωεμ 可见,ρ(r′)e-jk|r-r′|是ρ(r′,t-|r-r′|/v)的复数形式,因此式(335)的复数形式为 Φ(r)=14πεVρ(r′)e-jk|r-r′||r-r′|dv′(370) 这就是无界空间中时谐场的复数形式标量位函数,即方程(367b)的解。 类似地,在无界空间中,时谐场的复数形式矢量位为 A(r,t)=μ4πVJ(r′)e-jk|r-r′||r-r′|dv′(371) 式(370)、式(371)表示的时谐场位函数,仍然是滞后位。相位因子e-jk|r-r′|表示波传播|r-r′|距离后相位滞后了k|r-r′|。在时谐场情况下,相位滞后实际上就意味着时间滞后,或者说,时间滞后可以用相位滞后来体现。 3.7案例与实践 根据3.5节的分析,时变电磁场是以波动的形式在空间传播的。在电磁波传播路径上的不同点处,场矢量的振动规律相同,仅是在时间上有滞后。本节主要对时谐电场的波动性进行仿真和可视化展示。 案例: 假设有时谐电磁场,其电场矢量为 E(r,t)=x^Ex(z,t)=x^Exmcos(ωt-kz+) 绘制该电场矢量在空间的传播特征,并观察z=λ和z=1.5λ这两点处电场矢量的振动状态。电场初始参数: f=1GHz,Exm=1,=0。 图314中,利用quiver3函数模拟了时谐电场沿z轴方向的传播过程,并且单独标示了z=λ和z=1.5λ两点处的电场矢量的振动状态。首先,从左图可以明显看到电场整体上沿+z轴方向推进。其次,可以看到,对于z=λ和z=1.5λ这两点,电场矢量的振动状态完全一致,但是后者较前者在相位上滞后180°,这主要是空间上的滞后引起的相位滞后。从另一个角度分析,可以认为z=1.5λ处的振动在重复z=λ处的振动,只是其时间上滞后了0.5个周期,符合3.5.2节分析的波动性特征。感兴趣的读者请查看本页的动画。 动画 图314时谐电磁场传播示意图 思考题 31什么是位移电流,它与传导电流有何区别与联系。全电流定律所表征的方程是什么? 32什么是边界条件?理想导体表面的电磁场有何特点? 33简述时变电磁场的矢量位与标量位的引入过程。时变场的位函数与静态场的位函数有何区别与联系。 34什么是洛伦兹规范,如果在引入时变场位函数的时候采用的是库仑规范,那么矢量位A和标量位Ф所满足的方程是怎样的? 35如何描述电磁能量流动的大小和方向?简述坡印廷定理。 36瞬时坡印廷矢量、复坡印廷矢量、平均坡印廷矢量各自的区别与联系是什么? 37什么是复坡印廷定理?其中各项的物理内涵是什么? 38时谐场复数形式麦克斯韦方程组与一般瞬时形式麦克斯韦方程有何区别与联系? 39时变场所满足的波动方程是什么?如何理解场的波动性? 310什么是滞后位?滞后位所代表的物理本质是什么? 练习题 31以铜为例,证明在微波频率(300MHz~3000GHz)下,良导体中的位移电流远小于传导电流。假设铜的介电常数ε=ε0,电导率σ=5.8×107S/m,其内部的传导电流为J0cosωtz^(A/m2)。 32设真空中的磁感应强度B=10-3cos(6π×103t)cos(2πz)y^(T),求位移电流密度。 33已知E=E0[cos(ky-ωt)x^+cos(ky-ωt)z^],式中E0、k、ω为常数,由麦克斯韦方程组确定与之相联系的B。 34一段由理想导体构成的同轴腔,内导体半径为a,外导体内半径为b,长为L,两端用理想导体板短路,导体之间为空气,如题34图所示。已知在a≤ρ≤b,0≤z≤L区域内的电磁场为 E=Aρsin(kz)ρ^, B=Bρcos(kz)θ^ (1) 确定A、B之间的关系; (2) 确定k; (3) 求同轴腔所有内壁上的自由电荷密度和传导电流密度。 题34图 35设z=0平面为两种磁介质的分解面,z<0空间的μr1=4,z>0空间的μr2=3,分界面z=0处的面电流为 Js=3sin(2π×106t)y^(A/m) z<0空间中紧挨z=0界面处的磁场为 H1=(2x^+3z^)sin(2π×106t)(A/m) 求z>0空间中紧靠z=0界面处的磁场H2、B2。 36求下列各复场量的瞬时表示式: (1) E=E0e-jβzx^; (2) E=E0sin(βz)x^; (3) H=jH0cos(βz)y^; (4) E=5ej30°x^+6ej220°y^+e-j40°z^。 37已知空气中H=-j2cos(15πx)e-jkzy^,f=3×109Hz,试求E和k。 38分别处于x=0,x=d两个平面上的两块无限大平行理想导体板之间的电场为E=y^Emsinnπdxe-jk2-(nπ/d)2z。 (1) 求H; (2) 写出E、H的瞬时表示式; (3) 求两导体板上的面电流密度。 39无限长理想导体所围区域为0≤x≤a,0≤y≤b,如题39图所示,其中电场为E=Ey0sinnπaxe-jβzy^。 题39图 (1) 求区域中的位移电流密度和磁场; (2) 求瞬时坡印廷矢量和平均坡印廷矢量; (3) 计算穿过该区域任一横截面的平均功率。 310真空中存在的电磁场为 E=jE0sin(kz)x^, H=E0ε0μ0cos(kz)y^ 式中: k=2π/λ=ω/c,λ为波长。 求z=0,λ/4各点的坡印廷矢量的瞬时值和平均值。 311无源(J=0,ρ=0)空气中两个不同频率的电磁波的电场强度分别为 E1=E1e-jω1z/cx^E2=E2e-jω2z/cx^ 式中: c=1/ε0μ0。 证明合成场的平均能流密度等于两个波各自的平均能流密度之和。 312已知E=E0cos(kz-ωt)x^+E0sin(kz-ωt)y^ (1) 求磁场、坡印廷矢量; (2) 对于给定的z(如z=0),确定E矢量随时间变化的运动轨迹; (3) 求电场能量密度、磁场能量密度和平均能流密度。 313证明在无源空间的电磁场可以用矢量位Am表示成E=-1ε×Am,H=-jωAm+(·Am)jωεμ,而Am是方程2Am+ω2εμAm=0的解。 314在麦克斯韦方程中,若忽略Dt和Bt,证明矢量位和标量位满足泊松方程2A=-μJ,2Φ=-ρε。 315证明矢量位A=sin(βy)cos(ωt)x^是无源媒质中方程2A-εμA2t2=0的解,其中β=ωεμ。求与A对应的电场强度、磁场强度。