第3章〓波导里的微波 3.1为什么要学这一章? 通过第2章的学习,我们知道了常用的几种双导体结构的传输线。面对这种类型的传输线,使用的分析方法都是电路基础中学到的,也就是“路”的方法,虽然学的过程也有点小辛苦,但是相比于电磁场相关课程中学过的那些知识,上一章的打开方式已经是相当友好了。然而,也希望大家不要抱有一种不切实际的幻想,以为微波这门课只搞“路”,不搞“场”。作为“电磁场”课程的后续,这门课就是为了解决微波应用中的实际问题,而很多微波结构的分析,光有“路”的方法还是远远不够的,肯定要一手抓“路”,一手抓“场”,坚持两手都要抓,两手都要硬(图31)。 图31微波工程师必备两大法宝 关于“场”和“路”两种分析方法的区别和联系: ①“场”的方法基于麦克斯韦方程,更接近物理本质,可对三维空间场分布进行分析,能够给出的信息也更丰富,涵盖“路”的方法可以提供的全部信息; 而“路”的方法基于基尔霍夫定律,处理的问题是一维的,能提供的信息较为有限,例如阻抗是由电压电流比值定义的,而电压和电流则来源于电场和磁场的积分; ②秉承“万物负阴而抱阳”的原则,一种东西肯定是两面的,“场”的分析方法缺点就是过程烦琐,计算量大,反观“路”的方法则较为简洁,对于一些简单工程实际问题更适用。 话说到这份上,也就不用再装了,本章就要开始用“场”的方法来分析一些问题了。分析什么问题呢?先从给第2章学过的双导体传输线“找碴儿”开始。比如平行双导线,结构的确相当简单,但是封闭性较差,容易造成辐射损耗,相当于一小部分能量在传输的过程中泄漏到周围空间了,传输线成天线了,其他的平行波导板、带状线以及微带线都有这么个毛病,封闭性不太好,遇到阴天下雨的户外场合,甚至都容易进水。再比如同轴线,密封性很不错,几乎没有辐射损耗,结构也算轻便,但是它最大的一个缺点就是可传输的功率比较小,平时整点功率比较弱的信号还行,但如果信号功率搞到几百上千瓦时,传输线一点点的损耗就会造成热量迅速提升,轻则击穿,重则原地自焚。有人可能觉得夸张了,觉得功率不至于这么大,但如果了解一些大型雷达或者微波武器之后,就知道达到这样的功率简直就是分分钟的事情,甚至连微波炉的功率都可以达到这个量级。此外,同轴线、微带线之类的传输线一般必须在两个导体之间填充介质起到支撑的作用,因此损耗也比较大,特别是随着频率的升高,介质带来的损耗会急剧变大。 那么问题来了,既要传输功率特别大,又要损耗特别小时,该用什么传输设备呢?其实前面说这么多,就是为了硬推本章的主角——波导,它大概长什么样子呢?先上图(图32),形成一个直观印象。 图32两种常见波导 波导的造型说白了就是金属管子,横截面是矩形的就叫矩形波导,横截面是圆形的就叫圆波导,这两种相对来说最为常见,当然还有椭圆波导、脊形波导以及扭波导等一些形状比较怪异的,一般专用于特殊场合或者纯粹是学术时尚,较为非主流,不在我们小菜鸟的考虑范围之内。在矩形波导和圆波导之间,本章要重点对付的是矩形波导,原因有二: 应用更广泛,结构更简单。波导这种憨厚的传输结构,一看就是能承受大功率的,而且封闭性也很好,虽然略显笨重,但是在有些场合还真的没它不行。此外,波导不填充介质,直接做成中空的,因此只有金属损耗。而就算这点儿金属损耗,人类也没打算放过,甚至会在波导的内壁镀上一层铜甚至金来进一步降低金属损耗。这么说来,波导还真有点“憨厚多金”的气质了。 3.2波导的故事 现在我们拿波导来传输电磁波感觉司空见惯,但在一开始这种想法是很前卫的。1885年,英国物理学家赫维赛德考虑过电磁波在封闭的空管 中传播的可能性,就是上面提到过的简化麦克斯韦方程组表达式的那位大佬, 图33双导体传输线TEM模电磁场分布示意图 但他很快放弃了,因为他认为必须用两根导体来传输电磁能量。这也难怪,毕竟当时的思维定式已经形成,电流在两个导体内纵向流动,电压在两个导体间形成,波在传输线中是横电磁(TEM)模式,即电场、磁场以及传播方向三者相互垂直,构成右手螺旋定则,在横截面上,电场从一个导体垂直指向另一个导体,磁场则是分别包围着两个导体与电场相互垂直,如图33所示。 从场的角度来看,这种传输模式的场分布和自由空间中的TEM模没什么本质区别,通过对电场、磁场进行积分也可以得出与第2章中同样的电压、电流以及阻抗等参数。 ※如何理解模式? 一说到模式(mode)的概念,好多同学都又开始慌了,这在数学上倒是很好说,就是麦克斯韦方程在特定的边界条件下的解,但这样说可能会使大部分人陷入更加懵圈的状态。为了更好地理解模式的概念,可以先思考一个简单的问题: 有一个大的正方形,边长是10cm,如果我们知道这个正方形是由一定数量的相同的正方形单元拼接而成的,那么这个大正方形应该是什么样子的呢?但凡上过幼儿园,基本上也可以猜到,可能的情况如下,这里的每种情况都可以称为一种模式。 在这个问题中,我们还要考虑一系列条件: 不能越过10cm×10cm的边界,单元之间又不能重叠,还必须是尺寸相等的正方形单元,等等,最终得出结论: 实际可能存在的构成模式可以是模式1(单元个数: 1个,单元边长: 10cm),或者模式2(单元个数: 4个,单元边长: 5cm),抑或是模式3(单元个数: 9个,单元边长: 10/3 cm)等。每个模式都对应着一种可能性,或者说是这个问题的一个解,至于到底采用哪种模式,就取决于很多实际因素了,比如我们手头有哪种尺寸的正方形单元,或者单元尺寸不能小于多少,或者我们更喜欢什么样的单元数目,等等。 其实吧,上述这样一个问题实在简单到令人发指,但是为什么我们还要煞有其事地扯这么多呢,就是为了让大家更好地理解什么叫模式。这里的模式和波导中的模式在本质上都是一回事,就是一个特定问题的不同可能性答案,或者叫解决方案。不同的是,一个是往大方块中填充小方块,另外一个是往金属管子中填充电磁波,所运用的知识种类及复杂度不同而已。 ※※ 要想用一个金属管子来传输电磁波,首先就需要在理论上证明其可行性,奥利弗·赫维赛德(图34)当时直接拿着TEM模式的场分布就开始往管子里塞,然而,脑补完管子中的场分布之后,他很快就发现了一个比较尴尬的矛盾闭环: ①由于TEM模式是横电磁,因此电场和磁场必须都在横截面上,没有磁单极子的情况下磁场又必须打圈圈,这样一来,横截面上就必须存在一个磁场小漩涡; ②产生磁场小漩涡的方法有二(参见图115),要么有一个传导电流,要么有一个位移电流; ③既然管子是空心或者介质填充的,传导电流就别想了,位移电流则要求有和小漩涡所在横截面垂直的时变电场,也就是说要有纵向电场分量,显然又不是TEM模式了。一圈矛盾论证下来,赫维赛德直接放弃了,认为用管子传点自来水还行,用来传输电磁波实在有点无厘头。这倒真不是赫老哥能力不行,大概是时间和精力实在有限,跑去思考别的人生命题了。 图34奥利弗·赫维赛德(英国,1850—1925) 李宗盛大哥曾经说过,感情说穿了,一人挣脱的,另一人去捡。就像感情一样,波导被赫维赛德放弃了之后,转头就被另一个人捡了起来。1897年,瑞利爵士(本名: John William Strutt,爵位是家族继承的男爵,图35)从数学上证明了电磁波在波导中传播是可能的,不过要用到高次模式: 横电(TE)和横磁(TM)模式 ,他同时指出,有无穷多个TE、TM模式,且存在截止频率。这里的TE模式和TM模式就是往波导中塞电磁波时可能存在的两类场分布情况,本质上就是麦克斯韦方程在金属管子中的解。至于截止频率的概念先往后放一放,人生初见,记着这个名词就好。 图35瑞利爵士(英国,1842—1919) 被瑞利爵士把玩了几天之后,波导在接下来将近40年(1897—1936)中,基本被人们逐渐遗忘了。这个也难怪,虽然其间爆发了第一次世界大战,但那场战争基本上还属于拼火力、拼人头的状态,与信息化、电子化完全不沾边,波导自然没什么用武之地。一般来说,这种新潮科技,如果在军事上都用不上,民用就更没戏了,因此被遗忘也属实正常。 故事的下一个转机发生在1936年,纽约AT&T公司的George C. Southworth在一次会议上宣读了关于波导的论文,同一次会议上,MIT的W.L.Barrow宣读了关于圆波导的论文,并提供了波传播的实验验证。至此,波导这样一个简单的金属管子可以用于传输电磁波的事情算是得到AT&T和MIT两大巨头的双重官方认证了(图36)。大家可以留意一下这个时间点,正是二战打得热火朝天时,而雷达首次应用于战争则是在1935年。 图36来自AT&T和MIT两大巨头的双重背书 至此,在经历了被放弃、被捡起、被遗忘、被证实以及被重用等不同境遇之后,波导终于在整个微波传输领域牢牢占据了一席之地。特别是在大功率微波传输场合,结构简单、憨厚老实的波导简直就是不二之选。 3.3如何拿捏波导? 听完波导的故事,就要开始考虑如何对其进行深入研究了。为什么像赫维赛德这样的大佬级人物都会在一开始就否定波导传输微波的可行性?主要就是因为和第2章中学过的传输线结构相比,波导的结构实在是太不相同了,从双导体结构直接变成单导体,搁谁一时间都有点接受不了。 好在瑞利爵士对其不抛弃、不放弃,硬生生从数学上给波导找了一条活路。既然TEM这种模式在封闭的波导中肯定不能存在,那就看看有没有别的模式。本着“往前一步是黄昏,退后一步是人生”的信念,TEM不行,TE行不行?TM行不行?别误会,后面一句真不是骂人啊,意思是允许电场或磁场含有沿着传播方向的分量。这样一来,格局马上就打开了。其中,TE就是横电(transverseelectric)模式,没有纵向电场,但是允许有纵向磁场; TM是横磁(transversemagnetic)模式,没有纵向磁场,但允许有纵向电场。接下来就可以顺着这个思路开展对波导的分析了。 3.3.1波导分析——纵向场量法 图37矩形波导结构示意图 在本章开头就已经铺垫过,分析波导时要用“场”的方法,毕竟对于这样一个单导体结构,连确定的电压和电流都找不到,基尔霍夫的那两条定律也就只能原地退下了。从几何角度来说,波导本质上是一个封闭的规则管状结构,我们要研究的内容就是搞清楚在这个管子中传输微波时其电场磁场是如何分布的。第一次分析波导,与之前的风格一样,弱水三千,只取一瓢,而且一定是最简单、最常用的那一“瓢”——矩形波导(图37)。 假定波导沿着x轴的边为长边,沿着y轴的边为短边,纵向沿着z轴方向,为传播方向。我们的目标就是找出波导内部传输微波时其电场和磁场的分布: E=Ex(x,y,z)ax+Ey(x,y,z)ay+Ez(x,y,z)az=Et+Ezaz H=Hx(x,y,z)ax+Hy(x,y,z)ay+Hz(x,y,z)az=Ht+Hzaz (3.1) 其中,Et和Ht是横向场分量,Ez和Hz是纵向场分量。注意,这里仍旧通过假定波随时间作简谐变化来分离时间和空间,因此只考虑波在空间上的分布变化。 求波导内场分布的过程其实就是在一个封闭的矩形管子中求解麦克斯韦方程组的过程,对于这样一个管状的空间,我们通常采用的方法称为纵向场量法[7]。为啥叫这个名字呢,因为这种方法是以纵向场分量(Ez和Hz)为突破口,通过一系列数学操作来搞定整个场分布的表达式。具体的过程如下: 为保证推导过程的简洁性和整体思绪的连贯性,这里先直接上一个数学结论: 在正交坐标系中,电场的纵向分量和磁场的纵向分量满足标量亥姆霍兹方程,即 2Ez+k2Ez=0(a)2Hz+k2Hz=0(b)(3.2) 其中,k=ωμε。 至于说为什么会有这个结论,大概是先无脑盲猜矩形波导中的电磁场是波,既然是波,那就归亥姆霍兹管,后来通过推导和实验发现在数学和物理上都说得通,那这事儿就这么定了。 解这样一个方程,需要用到分离变量法,令 Ez(x,y,z)=Ez(x,y)Z(z)=Ez(T)Z(z)(a)Hz(x,y,z)=Hz(x,y)Z(z)=Hz(T)Z(z)(b)(3.3) 将式(3.3a)代入式(3.2a)后,可得 1Z(z)d2Z(z)dz2=-1Ez(T)(2t+k2)Ez(T)(3.4) 观察式(3.4),左右两边函数的变量完全不一样居然还能保持相等,因此二者只能同时等于某一常数,假设此常数为γ2,则有 2tEz(T)+(k2+γ2)Ez(T)=0(3.5) d2Z(z)dz2-γ2Z(z)=0(3.6) 看着γ是不是有点眼熟?别怀疑,这就是第2章中那个传播常数γ,反映波在传播过程中幅度和相位的变化,我们只是从场的角度又把它给推导出来了,再次证明了“路”和“场”的殊途同归。 使用同样的手法,可以得到磁场横向分量满足的表达式, 2tHz(T)+(k2+γ2)Hz(T)=0(3.7) 式(3.6)是一个二阶齐次常微分方程,通解很简单,就是一个波沿着z轴传播的形式, Z(z)=A+e-γz+A-eγz(3.8) 这个通解形式也并不陌生,波导中传输的是沿着+z轴和-z轴两个方向传输的波。鉴于后续整个电磁场表达式的复杂性,我们先主动给自己减负,只考虑沿着+z轴传播的分量,并且假定波导是无耗的,则有 Z(z)=A+e-jβz(3.9) 式(3.5)和式(3.7)则为二维标量亥姆霍兹方程,可以通过再一次分离变量并结合具体的边界条件来求得Ez(T)和Hz(T)。 随时间作简谐变化的场其空间表达式满足微分形式麦克斯韦方程, ×E(x,y,z)=-jωμH(x,y,z)(a)×H(x,y,z)=jωεE(x,y,z)(b)·H(x,y,z)=0(c)·E(x,y,z)=0(d) (3.10) 将式(3.1)代入式(3.10),整理后可得 Ex=-1k2cjωμHzy+jβEzx(a) Ey=1k2cjωμHzx-jβEzy(b)Hx=1k2cjωεEzy-jβHzx(c)Hy=-1k2cjωεEzx+jβHzy(d) (3.11) 其中,k2c=k2+γ2。 这表明,所有的横向场分量(Ex,Ey以及Hx,Hy)都可以由纵向场分量(Ez,Hz)通过式(3.11)求出来,而在前述的推导中,我们已经知道了如何求得Ez,Hz,因此也就可以顺势求出波导中完整的电磁场表达式。 上述纵向场量法的推导过程可以总结为一个流程图,更便于理顺整体逻辑关系,如图38所示。 图38纵向场量法流程图 接下来,就可以根据这个纵向场量法流程图的操作手法,对矩形波导进行进一步剖析了。 3.3.2矩形波导的TM模和TE模 如前所述,虽然TEM模式很简单,电场、磁场和传播方向两两垂直,但在矩形波导这样一个封闭的管子中实在是难以容身,因此我们尝试考虑TM或者TE的模式,也就是允许在传播方向上出现电场分量或者磁场分量[8]。接下来就要对纵向场量法进行实操了。 参照式(3.5)和式(3.7),纵向电磁场分量的横向函数Ez(T)和Hz(T)满足亥姆霍兹方程, 2tEz(T)Hz(T)+k2cEz(T)Hz(T)=0(3.12) 对于直角坐标系,在无耗情况下,可展开为 2x2Ez(x,y)Hz(x,y)+2y2Ez(x,y)Hz(x,y)+k2cEz(x,y)Hz(x,y)=0(3.13) 对Ez(x,y)采取分离变量法, Ez(x,y)=X(x)Y(y)(3.14) 则有 X″(x)X(x)+Y″(x)Y(x)=-k2c(3.15) 欲使式(3.15)恒等,左边两项应各等于一个常数,则有 X″(x)X(x)=-k2x,Y″(x)Y(x)=-k2y(3.16) 由此得出 X(x)=A1e-jkxx+A2ejkxx(3.17) Y(y)=B1e-jkyy+B2ejkyy(3.18) 加入纵向电磁场分量的纵向函数Z(z), Ez(T,z)=Ez(T)Z(z)=Ez(x,y)Z(z)(3.19) 由此可得电磁场的纵向分量完整形式: Ez=Ez(x,y,z)=(A1e-jkxx+A2ejkxx)(B1e-jkyy+B2ejkyy)C1e-jβz(3.20) Hz=Hz(x,y,z)=(A3e-jkxx+A4ejkxx)(B3e-jkyy+B4ejkyy)C2e-jβz(3.21) 其中,A1、B1、C1等以及kx,ky均为待定的系数,由具体的边界条件和信号的大小决定。 有了Ez,Hz,就可以通过式(3.11)来求得完整形式的电磁场表达式了。 以上步骤对于TM和TE模式都是通用的,接下来就开始要分头行动了,先说TM模式。 1. TM模式 TM模式顾名思义就是横磁(Transverse Magnetic)模,磁场一定要在横截面上,电场的话,其分量可以在横截面上,也可以沿着信号传播的方向。即,Ez ≠0,Hz=0。这样就可以通过边界条件来确定系数了。对于电场来说,所谓的边界条件,就是金属内壁上,不能有切向的电场分量,即 (an×E)|s=0(3.22) 具体来说,一共4个内壁,全用上, Ez|x=0,a=0,Ez|y=0,b=0(3.23) 将式(3.20)代入式(3.23),可得 Ez=-4A2B2C1sinkxxsinkyye-jβz=E0sinmπaxsinnπbye-jβz(3.24) 其中,m=1,2,3,…, n=1,2,3,…, E0=-4A2B2C1。E0为某种模式的幅度,由TM模式的激励信号来决定。 因为Hz=0,将式(3.24)代入式(3.11)可以求得电磁场的横向分量, Et=-jβk2cEzxax+Ezyay=Exax+Eyay =-jβk2cE0mπacosmπaxsinnπbyax+nπbsinmπaxcosnπbyaye-jβz (3.25) Ht=-jωεk2c-Ezyax+Ezxay=Hxax+Hyay=jωεk2cE0nπbsinmπaxcosnπbyax-mπacosmπaxsinnπbyaye-jβz (3.26) 至此,TM模式下电磁场的所有场分量都已经到齐了,就是式(3.24)、式(3.25)和式(3.26)所示的那样。不用多说,大家肯定已经开始头皮发麻了,所以也就不再把它们写在一起了,不然看着更难受。 尽管形式上很复杂,但是我们还是硬着头皮看一下,最起码要先有一个整体的感性了解。 首先我们要建立一个认知,TM模并不是一个模式,而是一类模式,说明在这个波导空间中不但可以传输电磁波,而且还用无数种模式花样传输电磁波,每个数字m和n的组合就代表了一种TM模式,所以要具体指出哪种模式时,需要加一个脚标,TMmn。想用哪种模式传输,只要搞清了该种模式的场分布,然后在合适的位置给出一个合适的电场激励或者磁场激励即可。 通过式(3.24)、式(3.25)和式(3.26)可知,m和n只要有一个是0,三个式子都会直接 归0,因此,TMmn最低次模就是TM11。 Ez和Ht二者叉乘后对应的传播方向在横截面上,且并不能产生实的功率密度,因此沿着x轴或者y轴,电磁场呈驻波分布,也就是有的位置振幅可以很大,有的位置振幅就小到0,电场或磁场振幅从0到最大再到0的变化过程为一个半波变化。m代表沿着x轴的半波变化个数,n代表沿着y轴的半波变化个数。 横截面中的电场Et和磁场Ht倒是很般配,二者可以直接叉乘出一个沿着z轴方向传播的行波,且对应的功率密度是实的,这也正是电磁波可以在波导中传播的部分。 2. TE模式 TE模式顾名思义就是横电(Transverse Electric)模,电场一定要在横截面上,磁场的话,可以在横截面上,也可以沿着信号传播的方向,即Hz ≠0,Ez=0。对于TE模,直接通过边界条件确定系数时稍显尴尬,因为对于磁场来说,导体内壁上的边界条件要求的是法向磁场为0,而Hz对于内壁来说是切向磁场,本身就可以不为0。为此,我们采取的策略就是: 脸皮一厚,硬往上凑。既然法向磁场要求是0,对于波导内壁短边来说,Hx就是法向分量,对于内壁长边来说,Hy也是法向分量,用Hz通过式(3.11c,d)把Hx和Hy硬凑出来,再套入边界条件即可。 Hx=-jβk2cHzx=βkxk2c(-A3e-jkxx+A4ejkxx)(B3e-jkyy+B4ejkyy)C3e-jβz(3.27) Hy=-jβk2cHzy=βkyk2c(A3e-jkxx+A4ejkxx)(-B3e-jkyy+B4ejkyy)C3e-jβz (3.28) 这时再利用边界条件就丝滑多了 Hx|x=0,a=0,Hy|y=0,b=0(3.29) 将式(3.27)和式(3.28)代入式(3.29)可得 A3=A4,B3=B4,kx=mπa,ky=nπb(3.30) 其中,m=0,1,2,3,…, n=0,1,2,3,…。由此可得到 Hz=4A3B3C3coskxxcoskyye-jβz=H0cosmπacosnπbe-jβz(3.31) 其中,H0=4A3B3C3,由TE模式的激励信号决定。 有了Hz,就可以通过式(3.11)来求出横向场分量的表达式了, Et=jωμk2c-Hzyax+Hzxay=Exax+Eyay=jωμk2cH0nπbcosmπaxsinnπbyax-mπasinmπaxcosnπbyaye-jβz (3.32) Ht=-jβk2cHzxax+Hzyay=Hxax+Hyay=jβk2cH0mπasinmπaxcosnπbyax+nπbcosmπaxsinnπbyaye-jβz(3.33) 与TM模式相类似,TE模也不是一个单一的模式,而是一类模式,要唯一地指出某一个具体的TE模,也要加上脚标,TEmn。与TM模式不同的是,TE模的脚标可以从0开始,也就是可以有TE10或者TE01模式的存在。当然,m和n也不能同时为0,否则式(3.31)~式(3.33)还是会瞬间归零。 3.3.3矩形波导的高通特性 说完TM和TE模式,用矩形波导这么一个封闭的管子来传微波这件事儿在数学上算是有了定论。不仅能传,还能用无数种花样来传,对应着无数个模式。是不是已经按捺不住激动的心和颤抖的手,准备拿着微波就往波导里面灌了?事情往往没这么简单,当你真正用波导传输微波时,就会发现一个很奇特的现象,对于某一种传输模式来说,只有信号频率够高,才能传,频率低于某个值时,波导是拒绝的,这个值称为截止频率,前面介绍过。也就是说,波导对于电磁波有高通的特性(图39),是不是很离谱? 图39波导的高通特性 至于为什么会产生这种高通的现象呢,这就要从kc、k和γ这几个物理量说起了。 (1) 先说k,就是前面说过的波数(wave number),是由波的频率和传输媒质的属性所决定的,即 k=ωμε(3.34) 这个量从麦克斯韦方程时空分离时就跟着了,在亥姆霍兹方程中更是各种露脸。 (2) 接着说γ,这就是传播常数,其虚部β称为相位常数。从式(3.8)可以看出其实部和虚部分别代表波沿着z轴传播时,幅度和相位随空间位置的变化快慢。无耗的情况下,γ=jβ。 (3) 最后说kc, k2c=k2+γ2=mπa2+nπb2(3.35) 这个值叫截止波数(cutoff wave number),不同的模式对应着不同的kc。之前在TEM模式中是见不到kc的,因为在TEM模式下,k=β,即kc=0。然而,在TE和TM模式下,kc就不能再是0了,不然的话,纵向场量法中的式(3.11)从一开始就要出问题了,因为kc直接出现在了分母的位置。此外,通过式(3.35)还可以看出,kc主要由TE和TM模式的编号m、n以及波导的尺寸a、b来决定。 把式(3.35)改成另外一种形式会更容易解释接下来的问题,即 β2=k2-k2c(3.36) 这个式子给出了一个重要的信息: k一定要大于kc才行。为什么呢?假设k<kc,那么β就变成了虚数,而式(3.9)中就会出现随着z的变化,原本是相位的变化,结果现在成了幅度的变化,而且无源情况下幅度肯定是迅速衰减的(此时β是一个负的虚数),肯定是不能传输的,这就是截止(cutoff)的含义。模式(m,n)和波导尺寸(a,b)定下来了,kc就定下来了,那对应的门槛就定下来了,k要高于这个门槛,即频率f要高于某一个值(ω=2πf),这个值就称为截止频率fc。 fc=ωc2π=kc2πμε=mπa2+nπb22πμε=12μεma2+nb2(3.37) 可见,一旦确定了波导尺寸(a,b),波导内填充媒质(μ,ε)和具体模式(m,n),截止频率就确定下来了,信号只有高于这个频率才能够进入波导进行传输。此处请注意,TEmn和TMmn模式的截止频率计算公式是一样的,只要脚标编号相同就具有同样的截止频率。 有截止频率就有截止波长,可以通过下式求出, λc=2πkc=2ma2+nb2(3.38) 所以下次微波信号再想进入波导,至少得先扪心自问一下,频率够高么?波长够短吗?这里可以通过一个具体的例子将波导“嫌低爱高,嫌长爱短”的特性好好展现一下。 由于在微波传输中的重要地位,波导也是有国标的,不同的型号对应着不同的尺寸和适用频率范围。在工程应用中,考虑到传输损耗、传输功率等实际因素,一般波导的长边和短边的大小关系为a>2b,略大于,这个是常识,需要记住。以型号为BJ32的标准波导为例,其尺寸a×b=(72.14×34.04)mm2,空气填充,由式(3.38)可求出其不同模式下的截止波长λcmn,以波长为横轴,可以画出一个不同模式下截止波长的分布图,由此可以非常直观地看出BJ32波导对于不同传输模式、不同信号波长的态度。 由图310可以看出,对于波长大于2a=14.28cm(频率低于2.1GHz)的微波信号,BJ32波导是彻底拒绝的,用啥模式都不好使,这个叫截止区; 当波长在a和2a之间时,只有TE10模式可以传输,其他的模式不可以,这个叫单模传输区,在实际工程应用中,波导绝大部分情况下都是工作在这个区域,因为这个区域的模式很确定,肯定是TE10,不会有其他高次模式的出现; 当信号的波长小于a时,各种高次模式就可能逐个出现了,显然,m、n的值越大,所对应的截止波长就越短,截止频率就越高。不难看出,各种模式都算在内,TE10模式是最低次的,也是最特殊的,我们把它称为矩形波导的基模或主模。鉴于TE10模在理论和实际中的重要性,我们将在后续用一个小节来细细品它。 图310BJ32波导不同模式截止波长分布图 3.3.4矩形波导的传输参量 微波信号的频率达到准入门槛之后,就可以开心地进入波导进行传输了,这时就不免要关心一下在波导中传输时,速度是多少?波长怎么算?这些问题在传输线时没有深入探究,是因为“路”的分析方法中,默认波就是TEM模式,速度和波长计算起来和空间波没什么差别。然而,涉及波导中,情况就稍显复杂了。 1. 相速度vp和群速度vg 沿着z轴传输的导行波其等相位面移动的速度就称为相速度(phase velocity),根据式(3.36)可以算出微波在波导中传输时的相速度为 vp=ωβ=ωk2-k2c=v1-λ/λc2(3.39) 这个速度由于只是单一频点波的某一相位的传播速度,并不携带能量,因此可以大于媒质中的光速v。 而真正可以表征电磁波能量传输速度的则是群速度(group velocity),其代表了频率和相位都很接近的一群波的传播速度,其计算式为 vg=v1-λ/λc2(3.40) 显然,群速度是小于媒质中的光速的。需要注意的是,在TE和TM模式下,无论相速度还是群速度,都是随着信号的波长(或频率)的变化而变化的,因此是色散波。而对于TEM模式来说,等相位面的传播和能量的传播是同步的,截止波长λc无穷大,其传输速度就是光速,与频率无关,是非色散波。 关于相速度和群速度的概念,可以观察调幅波的动图(别往下找,没有,直接上网搜),盯着某一个波峰或者波谷看,其移动的速度就是相速度; 盯着一个包络看,包络移动的速度就是群速度。如果还是觉得不够直观,那就类比一下搬砖,一堆砖你要从A点全搬到B点,单个砖块移动的速度大概就是相速度,整个砖堆移动的速度就类似于群速度。如果还是不明白,那就安心搬砖,没必要跟这个较劲。 2. 波导波长λg 波导波长定义为波导中的波在一个周期的时间内沿z轴传输的距离,记为λg,其计算式为 λg=vpT=vpf=λ1-λ/λc2(3.41) 可以看出,TE和TM模式的波导波长计算式是相同的,相比于自由媒质中传输的TEM波的波长λ要长一些的,毕竟相速度高于光速,在一个周期内的传播距离也会更远。 3. 波阻抗 阻抗的概念是电压和电流的比值,在第2章中很是常见,但是由于波导是用“场”的方法来分析的,且波导内部的电磁场并不均匀,无法通过积分找到确定的电压或者电流,因此这里采用波阻抗的概念,是横向的电场Et和磁场Ht的比值。 ZTM=EtHt=βkZTEM(3.42) ZTE=EtHt=kβZTEM(3.43) 其中,ZTEM是TEM模式下的波阻抗,在真空中为120π。对于TM模来说,由于磁场都集中在横截面上,导致电场比磁场时,分母更大一些,因此TM模的波阻抗ZTM相比于ZTEM更小; 对于TE模式来说,由于电场都集中在横截面上,因此TE模的波阻抗ZTE相比于ZTEM更大。 3.3.5矩形波导的基模: TE10 前面说到过,所有模式中最靓的那个仔就是TE10模式,也是矩形波导的主模(图311)。这个倒也不是硬捧,TE10模式之于矩形波导就是最重要的,没有之一。原因主要有以下几个: (1) TE10是唯一可以实现单模传输的模式,通俗点讲就是“我能传你们不一定能传,我不能传你们都别想传”。不要小看单模传输,每个模式都有不同的场分布,越高次场分布就越复杂,因此微波工程师是很不待见高次模的,根据图310适当地选择一个波导的工作频段范围,使得其中只有最简单的TE10模式难道不香吗? (2) TE10模式下波导的尺寸可以最小。根据式(3.37)可知,如果工作频率确定了,那么截止频率就要低于这个工作频率以保证可以传输,此时,模次越高,m、n就越大, 图311矩形波导的主模——TE10模 那么就相应地要求a和b也要变大才能使得截止频率足够低,因此,越高次模就越需要越大的波导尺寸。作为最低次模的TE10则可以在相对更小的波导中传播,又省空间又省材料难道不香吗? (3) TE10模的场结构分布最为简单。根据m、n的物理意义,沿着x轴只有一次半波变化,沿着y轴干脆就没有变化,稍微努力思考一下就可以脑补出场分布的情形,难道不香吗? 了解了TE10的江湖地位之后,就可以按照3.3.4节的结果对其进行全方位拿捏了。 1. 截止频率和截止波长及其他传输参量 截止频率: fc10=12με1a2+0b2=v2a(3.44) 截止波长: λc10=21a2+0b2=2a(3.45) 相速度和群速度: vp10=v1-λ/2a2,vg10=v1-λ/2a2(3.46) 波导波长: λg10=λ1-λ/2a2(3.47) 波阻抗: ZTE10=ZTEM1-λ/2a2(3.48) 2. 场结构 将m=1,n=0代入式(3.31)~式(3.33),可得TE10模的各个场分量表达式为 Hx=jβaπH0sinπaxe-jβz Hz=H0cosπaxe-jβz Ey=-jωμaπH0sinπaxe-jβz (3.49) 将时间的信息加入进去的话,就变成 Hx=jβaπH0sinπaxsin(ωt-βz) Hz=H0cosπaxcos(ωt-βz) Ey=-jωμaπH0sinπaxsin(ωt-βz) (3.50) 上面的式子虽然已经尽量简洁,但终归还是有点抽象,需要和更形象的图312和图313配合使用。 图312矩形波导TE10模的电场分布(Ey) 根据式(3.49)可以看出,TE10模电场只有Ey分量,磁场则有Hx和Hz分量,因为n=0,因此电磁场沿着y轴是均匀的,没有变化。如图312所示,电场分量Ey沿着x轴呈正弦变化,且从x=0到x=a呈现从0到最大值再到0的一个半波变化(m=1)。此外,电场分量Ey沿着z轴也呈正弦分布,其分布周期就是我们之前提到的波导波长λg。 如图313所示,磁场分量Hx沿着x轴呈正弦分布,和Ey有点像,也是两头为0,中间最大,经历一个半波变化(m=1)。Hx沿着z轴也呈正弦分布,其分布周期也是波导波长λg。 图313矩形波导TE10模的磁场分布(Hy和Hz) 磁场分量Hz沿着x轴呈余弦分布,两头最大,中间为0,沿着z轴也呈余弦分布,分布周期为波导波长λg。 磁场的两个分量Hx和Hz,一个是正弦分布,另一个是余弦分布,因此表现出“你大我小,你小我大”的对应规律,叠加之后,就会形成打圈圈的磁场,相邻的两个圈圈旋转的方向相反,构成一个波导波长λg,因此在波导中,电场像沿着短边的一根根柱子,顶天立地; 磁场则像一个个轮子,旋转前行。两个磁场“轮子”构成一个波导波长λg,“轮子”之间的交界处电场最强。 将电场和磁场画到一起,就形成了图314所示的透视图,电磁波以TE10模式进入波导之后,电场沿着y轴方向时上时下,磁场则打着圈圈,电场和磁场就这样相互激励着一路前行。 图314矩形波导TE10模电磁场分布透视图 注意,图312~图314相当于给电磁场拍的一幅照片,是定格在某一个时间点的,至于场随着时间是如何变化的,可以根据式(3.50)自行脑补,如果实在脑补不能,也可以自学一下HFSS、CST等电磁仿真软件,软件中可以生成动画,看着还挺带劲的。 为什么要花这么大力气来了解波导内部TE10模式的场分布呢?因为只有了解了场的分布,很多实际的工作才可以开展,比如怎样激励出TE10模式的波?波导和其他类型传输线之间的转接如何设计?等等。这些在实际工程中经常会遇到的问题都要以充分了解波导内部电磁场的分布为基础来解决。以TE10模的激励为例,如果知道了电场在某个地方应该是什么样子的,我们就可以用一个小探针放在相应的位置去激励这样一个电场,然后电场带动磁场,磁场带动电场,相应的模式就形成了。当然也可以用小磁环或者孔缝等结构在合适地方激励出相应的磁场,进而形成我们需要的模式。 3. TE10模式的面电流分布 之前提到过波导内部没有唯一确定的电压和电流,但并不代表完全没有电压和电流,只是因为场结构较为复杂,在一个区域去对电场或者磁场进行积分,处处都能积出不同的电压或者电流。由于趋肤效应,波导内壁上感应出的高频电流只在其内壁表面上流动,接下来将对波导内壁上的面电流分布进行介绍。面电流的分布取决于波导内壁表面附近的磁场分布,因此面电流密度Js可以这样计算: Js=an×Hτ(3.51) 其中,an为内壁表面法向方向的单位矢量,Hτ是内壁表面出的切向磁场分量。 在y=0,y=b,x=0,x=a的各个内壁上,均可以通过式(3.51)来计算相应的表面电流,结果如下, Js|y=0=(ay×H)|y=0=H0cosπaxax-jβaπH0sinπaxaze-jβz Js|y=b=-Js|y=0=-H0cosπaxax+jβaπH0sinπaxaze-jβz Js|x=0=(ax×H)|x=0=-H0e-jβzay Js|x=a=Js|x=0=-H0e-jβzay (3.52) 由式(3.52)可知,波导宽壁上下内表面上的电流由两个分量构成,两个壁上的电流图案相同,方向相反。窄壁上的电流分布都是沿着y轴方向,两个壁上的电流方向是完全相同的。 更直观的面电流分布可以看图315,在宽壁上,有的位置就像“泉眼”,电流咕嘟嘟往外冒, 图315矩形波导TE10模磁场及内壁面电流分布 有的地方就像“黑洞”,电流都往里流。上面宽壁是“泉眼”的位置,正下方的宽壁上就是“黑洞”,电流或者通过宽壁的面从一个“泉眼”流向相邻的“黑洞”,或者通过窄壁的面从一个“泉眼”流向正下方的“黑洞”。 了解了面电流的分布,有一个非常直观的好处,就是可以给波导进行开缝操作了。开缝?这在实际工程应用中还真的是一个常规操作。给波导开缝就好像给人开刀,医生必须深入了解人体的结构,不然的话一不小心割破了血管,后果很严重。面对波导也是如此,我们有时会希望在不影响场结构的情况下,给波导开个缝,然后把探针放进去,探测一下里面的场强,从而可以得到驻波比之类的参数,这时就要沿着电流的方向去开一条细细的缝,这样既不影响波导内的信号传输,又不造成波的泄漏,比如波导驻波测量线就是利用这种办法,在宽壁的中心线上开缝。当然,对于工程师来说,有不想让信号泄漏出来的,就有想让信号泄漏出来的,比如搞天线的工程师,做梦都想让信号更多地泄漏出去,只不过这时不叫泄漏,叫辐射。因此拿到波导之后,天线工程师的开缝原则就是哪里电流强切哪里,横着切,一定要把电流切断的那种,这就是波导缝隙天线,也是目前天线大家族中比较拉风的一类。两种切法,如图316所示,充分体现了同一个波导的两种妙用,颇有些庄子“齐物”的妙处,东西还是同一个东西,成为什么样的“器”完全就是人为决定了。 图316同一个波导的两种妙用 以上基本就把矩形波导的相关内容介绍得差不多了,当然,圆波导在实际中也有一定的应用,鉴于其大体套路和矩形波导没有本质区别,同时鉴于其数学推导的烦琐更甚于矩形波导,因此本书决定避其锋芒, 仅作简单介绍,不再过多赘述。 3.3.6圆波导的浅尝辄止 了解了矩形波导,对于波导的整体认知也就建立起来了,实际中矩形波导的应用也是最为广泛的。当然,既然都讲波导了,不提一下圆波导终归感觉有点儿不太好意思,但是如果讲得太深入也有点儿力不从心,因此本小节对圆波导浅尝辄止,顺着矩形波导的思路,大概了解一下圆波导的特点即可。 图317柱坐标系中的圆波导 既然是圆波导,那么横截面上的场分布再用直角坐标系就不太合适了,不过传播方向上的确还是要用z轴的,因此,分析圆波导一般就用柱坐标系,如图317所示,对于圆波导的数学分析也没什么新的东西,本质上也是在柱坐标系中去解亥姆霍兹方程,遵循上面说过的纵向量法的流程,最终解出的圆波导中传播的模式有两种,分别称为TEmn模(H模)以及TMmn模(E模),因为横截面也是二维的,且沿着R或者φ的方向都有可能变化,因此脚标也是两个数,m和n。这里略去推导过程,直接给出结果。 1. 模式的场表达式与截止波长 对于圆波导中TM模(E模),Hz=0,Ez≠0,其场 的横向分量 表达式为 Et=-jβk2ctEz=-jβk2cEzrar+1rEzφaφ=Erar+Eφaφ=-jβkcE0J′mkcrcosmφsinmφe-jβzar±jβmk2crE0Jmkcrsinmφcosmφe-jβzaφ(3.53) Ht=-jωεk2caz×tEz=-jωεk2c-1rEzφar+Ezraφ=Hrar+Hφaφ=jωεmk2crE0Jmkcrsinmφcosmφe-jβzar-jωεkcE0J′mkcrcosmφsinmφe-jβzaφ(3.54) 其中,Jm(kcR)=0为第一类m阶贝塞尔函数,kc是截止波数,由边界条件决定,这里的边界条件是指圆波导的内壁上,切向电场为零这类已知条件,即r=R时,Ez=0,Eφ=0,根据截止波数可得截止波长为 λcTMmn=2πRνmn(3.55) 其中,νmn= kcR,是m阶贝塞尔函数的第n个根,这种函数虽然工科学生看着复杂,但是理学院的那帮学数学的学生已经贴心地把表都给制好了,可通过查表得到νmn的值,如表 31所示。 表 31部分νmn的值 m n 1 2 3 0 2.405 5.520 8.654 1 3.832 7.016 10.173 2 5.135 8.417 11.620 3 6.379 9.761 12.015 同样地,对于圆波导中TE模(E模),Ez =0,Hz≠0,其场 的横向分量 表达式为 Et=jωμk2caz×tHz=jωμk2c-1rHzφar+Hzraφ=Erar+Eφaφ=±jωμmk2crH0Jmkcrsinmφcosmφe-jβzar+jωμkcE0J′mkcrcosmφsinmφe-jβzaφ(3.56) Ht=-jβk2ctHz=-jβk2cHzrar+1rHzφaφ=Hrar+Hφaφ=-jβkcH0J′mkcrcosmφsinmφe-jβzar±jβmk2crH0Jmkcrsinmφcosmφe-jβzaφ(3.57) 截止波长为 λcTEmn=2πRμmn(3.58) 其中,μmn=kcR,是m阶贝塞尔函数导数的第n个根,其值同样可通过查表得到,如表 32所示。 表 32部分μmn的值 m n 1 2 3 0 3.832 7.016 10.173 1 1.841 5.332 8.536 2 3.054 6.705 9.969 3 4.201 8.015 11.346 2. 最低次模TE11的场结构 既然有截止波长,说明圆波导也有高通特性,根据式(3.55)和式(3.58)可算出不同模式的截止波长,如图318所示。当然,通过截止波长还是以一溜烟儿算出截止频率fc、相移常数β、相速度vp、群速度vg、波导波长λg、波阻抗Z等参数,这里不再赘述。 图318圆波导中不同模式截止波长的分布情况 可见,按照截止波长来看,在所有模式中,TE11模式是最低次模,其截止波长最长,但是和矩形波导不同的是,我们一般不太情愿将圆波导TE11模式称为主模,这是因为就算在实际应用中,TE11模也不能算是最常用的,在TE11之外,大家反而还会选择用TM01和TE01两种,这就和矩形波导的情况非常不一样了,毕竟在矩形波导中,TE10模简直就像神一样的存在, 如果工程师不是为了博出位或者脑子被夹了,一般不会用TE10之外的其他模式。既然这样,为什么到了圆波导,最低次模TE11就不太受待见了呢,可以通过其场分布一探究竟,如图319所示,注意其中图319(c)相当于把波导沿着轴线劈成两半,然后把上半部分铺平了画出其内壁上的面电流分布。 图319圆波导中TE11模场结构图与面电流分布图 从图319(a)可以看出,无论是沿着R的方向,还是φ的方向,场分布都不是均匀的,问题就出在这儿了。可以看出穿过圆心的那条电场线,是水平方向的,其他电场线虽然弯曲,但是也大体和这条电场线保持一致。学过电磁场的同学应该知道,电场线末端划过的方向就是极化方向,也就是说,图319(a)中极化方向是水平的,然而,在圆波导中,这种极化方向是相当不稳定的,激励时手法稍微有点偏差,或者传播过程中波导结构稍微有点不均匀性,都会导致这个极化发生偏转,虽然整体场结构看上去都还是一样的TE11模,但是极化方向已经不一样了,可能变成斜的,甚至是垂直的,这种情况称为极化简并,如图320所示。 图320圆波导中TE11模极化的旋转 也就是说,实际中真要拿圆波导来传输TE11模式,可能传了一段距离之后,横截面上的极化方向到底指哪都不能确定了,这也是为什么在介绍图318时,压根就没有提单模传输区,因为在TE11和TM01的截止波长所包夹的区域中,虽然都是TE11模,但极化方向却是不确定的,并不能看成严格意义上的单模传输,或者说,圆波导就没有所谓的单模工作区。这就是人们不待见圆波导中TE11模的原因,虽然它的确是实实在在的最低次模。 3. TM01及TE01场结构 不喜欢圆波导TE11模的原因,也正是人们会常用圆波导TE01或者TM01模的原因,光看角标就知道, m=0意味着这两种模式下横截面上沿着φ的方向,场分布是均匀的。也就是说,没有极化简并的现象,极化方向爱怎么转就怎么转,怎么转都是一样的,这就是相比于圆波导中TE11模的好处。图321和图322分别绘出了TM01和TE01模式的场结构图,以及内壁上的场分布。注意其中的图(c)相当于把波导整体铺平了画出其内壁上的面电流分布。 图321圆波导中TM01模场结构图与面电流分布图 图322圆波导中TE01模场结构图与面电流分布图 由此可以看到圆波导一个比较尴尬的事实,最低次模存在极化简并,不存在极化简并的模式又不是最低次模,单模传输是不要再想了,因此客户体验相比矩形波导还是差了一点点。但这也并不是说圆波导就一无是处了,周星驰曾经说过,哪怕是一张卫生纸,都有它的用处。圆波导在某些场合还真是矩形波导无法取代的,比如雷达系统中,波导某些部分需要能够360°旋转的活动关节,那么这种关节必须得用圆波导,而且需要使用TM01或者TE01模式,这种情况下矩形波导还真就不灵了。 3.3.7“场”的角度回望微波传输线 前面说到过“场”的方法是最本质的,相比于“路”的方法,不仅可以提供更加丰富的信息,其适用的对象也更广泛。即便对于双导体的传输线,也可以完全用“场”的方法进行分析,不过会更复杂。这一小节我们只是通过一些图片和文字来从“场”的角度回望一下第2章学过的几种典型传输线,看看它们内部的场分布是如何的,不涉及公式推导。 图323同轴线柱坐标结构图 首先建立一个常识性的认知,对于双导体的传输线,其基模一般都是TEM模(微带线除外),但这并不代表它们不可以传输高次模,相反,双导体传输线都是可以传输更高次的TE或者TM模的,只要频率足够高,激励手法到位,没什么高次模是不能传的。 接下来以几种典型的传输线结构为例介绍场分布。 1. 同轴线 同轴线柱坐标结构图如图323所示。 同轴线的主模为TEM模式,柱坐标系下,其电磁场的表达式为 E=-te-jkz=ddre-jkzar=V0rln(b/a)e-jkzar(3.59) H=1ZTEMaz×E=V0rln(b/a)ZTEMe-jkzar(3.60) 式(3.59)以及式(3.60)描绘出来就是图324。 图324同轴线内部电磁场分布图 可以看到,同轴线里面TEM模式的电磁场分布还是比较容易想象的,横截面上,电场像辐条,磁场绕着内导体打圈圈,越靠近内导体场强越强。由于是TEM模,没有截止频率和截止波长,甚至直流都可以传输,此时的波数k和传播常数β是一回事。 当然同轴线也可以传输TE或者TM模式,相较于TEM模式,场分布就没那么简单了。为使读者情绪稳定,下面省略推导过程和相关数学表达式,直接给出几个高次模式的场结构图(图325),浅尝一下就好。 图325同轴线高次模的电磁场分布图 2. 带状线和微带线 同轴线虽然封闭性很好,但是有一个比较尴尬的缺点就是很难做到低剖面,也很难和其他小型化微波器件高度集成,因此不利于实现大规模的微波集成电路。微波工程师给出的替代方案是带状线或微带线。 先说带状线,又称为对称微带,是一种可以用印制电路板(PCB)技术实现的微波传输结构,可以看成同轴线的扁平化的产物,说白了就是把同轴线外导体切开,然后内外导体都压扁成一层薄薄的金属,中间用介质板来支撑。这样就可以实现电路板状的超低剖面,同时也方便微波元器件的集成。因此在微波集成电路(MIC)中,基本上采用的都是带状线或者微带线。注意,图326中带状线的结构只是为了让大家看清楚,真正的带状线的厚度肯定远低于同轴线的直径,也就差不多是普通电路板的厚度。 图326同轴线到带状线的演变过程 带状线的主模是TEM模,其结构示意图如图327(a)所示,其横截面上场分布如图327(b)所示。 图327带状线 微带线的结构就更加简单了,相当于带状线去掉微带上方的介质板和金属地。大家天天爱不释手的手机中就有大量的微带线来传输微波信号。微带线也可以看成平行双导线的扁平化设计结果,其演变过程如图328所示。 图328平行双导线到微带线的演变过程 由于微带线的横截面结构不具备对称性,而且微带的上方是空气,下方是介质和金属地,这导致严格意义上的TEM模并不能存在于微带线之上。微带线的主模是准TEM模,本质上是TE模和TM模组合而成的混合模式,只不过和TEM模的场分布很相似,因此我们称其为准TEM模,就像准新娘、准新郎之类的,很接近了,但还不是。 微带线的主模场分布如图329所示。 图329微带线结构示意图及电磁场分布图 可见在介质板和空气的分界面处,电磁场出现了明显的不连续性,这也决定了微带线的主模不可能是纯净的TEM模式。