4.1 Drátový dipól

Základní teorie

V této kapitole se seznámíme s výpočtem parametrů drátového dipólu pomocí momentové metody. Veškeré informace se snažíme co nejsrozumitelněji vysvětlit ve vrstvě A. Vrstvu B v tomto případě využíváme k uvedení anglické verze této kapitoly. Činíme tak proto, abychom čtenáře seznámili s anglickou terminologií, využívanou v oblasti antén a numerických metod.

Obr. 4.1A.1
Obr. 4.1A.1Drátový dipól

Všechny důležité technické parametry antén, jakými jsou např. zisk, vstupní impedance nebo směrová charakteristika, mohou být relativně snadno vypočteny, pokud známe rozložení proudu na anténním vodiči. Výpočet rozložení proudu je však bohužel dosti komplikovaný, protože při jeho určování je třeba řešit integrální rovnice.

K řešení integrálních rovnic existují dva základní přístupy – iterační a momentový. Iterační metody vycházejí z hrubé aproximace proudového rozložení (uvažujeme např. sinusové rozložení proudu na drátovém dipólu), která je iteračně zpřesňována. Oproti tomu momentové metody převádějí řešení integrální rovnice na řešení soustavy rovnic lineárních, s nimiž si bez problémů poradí např. Matlab.

V této kapitole naší učebnice se budeme zabývat pouze momentovou analýzou drátových antén. Ve všech případech budeme předpokládat, že anténa je tvořena kruhovým válcem o poloměru a a že je dlouhá 2h. Osu anténního vodiče umístíme do osy z (obr. 4.1A.1) válcového souřadného systému (r, ρ, z). Dále předpokládáme, že se anténa nachází ve vakuu (μ = μ0, ε = ε0, σ = 0) a že veškeré možné ztráty jsou nulové.

Obr. 4.1A.2
Obr. 4.1A.2Budicí elektrické pole ve vstupní štěrbině

Uprostřed anténního vodiče (z = 0) budeme uvažovat krátkou štěrbinku. Tuto štěrbinku připojíme k hypotetickému harmonickému zdroji, který vytváří rotačně souměrné budicí pole (obr. 4.1A.2). Napětí ve štěrbince pak můžeme popsat vztahem

V=-gapEzdz. ( 4.1A.1 )

Dále předpokládáme, že toto napětí je rovno jednomu voltu. Ve vztahu (4.1A.1) značí Ez z-ovou složku intenzity budícího elektrického pole ve štěrbině, a tedy i na (válcovém) povrchu této krátké části antény (obr. 4.1A.2). Mimo štěrbinu je Ez nulové, protože předpokládáme dokonalou elektrickou vodivost vodiče antény.



I. Momentová metoda

Uvažujme obecnou integrální rovnici ve tvaru

abf(z,ξ)dξ=g(z), ( 4.1A.2 )

kde f je neznámá funkce (v našem případě rozložení proudu na anténě), <a,b> značí analyzovanou oblast a g je známá funkce, popisující působení zdrojů. Momentové řešení rovnice (4.1A.2) potom můžeme rozepsat do následujících tří kroků:

  1. Neznámou funkci f aproximujeme pomocí lineární kombinace známých bázových funkcí fn a neznámých koeficientů cn
    ff˜=n=1Ncnfn. ( 4.1A.3 )
  2. Formální aproximaci (formální proto, že neznáme koeficienty cn) hledané funkce f~ dosadíme zpět do řešené rovnice a zaměníme pořadí sčítání a integrování
    n=1Ncnabfn(z,ξ)dξ=g(z)+R(z). ( 4.1A.4 )
    V uvedeném vztahu značí R(z) tzv. reziduum, které vyjadřuje skutečnost, že aproximace řešení f~ není identická se zcela přesným řešením rovnice. Vztah (4.1A.4) je jednou rovnicí pro N neznámých koeficientů cn.
  3. Co možná nejpřesnější aproximaci řešení získáme tehdy, když reziduum R bude minimální. Reziduum tudíž minimalizujeme tzv. metodou vážených reziduí: součin vhodné váhové funkce w a rezidua R, integrovaný přes analyzovanou oblast <a, b> musí být nulový [5]. Pokud pro váhování postupně použijeme N různých váhových funkcí, získáme soustavu N lineárních rovnic pro N neznámých koeficientů cn
    abwm(z)R(z)dz=0m=0,1,...N, ( 4.1A.5a )
    n=1Ncnabwm(z)abfn(z,ξ)dξdz=abwm(z)g(z)dz. ( 4.1A.5b )
    Jak váhové funkce tak funkce bázové musejí být lineárně nezávislé na intervalu <a,b>.

II. Bázové funkce

Bázové funkce mohou globální nebo lokální povahu. Globální bázové funkce jsou definovány přes celou analyzovanou oblast <a,b>. Např. soustava funkcí

fn(z)=cos(πnzh) ( 4.1A.6 )
Obr. 4.1A.3
Obr. 4.1A.3Vícebázová aproximace a) po částech konstantní, b) po částech lineární

je lineárně nezávislá na <a,b> a koeficienty cn

f(z)f˜(z)=n=1Ncnfn=n=1Ncncos(πnzh) ( 4.1A.7 )

zde mají význam Fourierových koeficientů proudového rozložení.

Aproximaci založenou na globálních bázových funkcích nazýváme jednobázovou aproximací.

Lokální bázové funkce jsou definovány přes celou analyzovanou oblast také, avšak pouze na určité podoblasti nabývají nenulové funkční hodnoty (obr. 4.1A.3). Pokud jsou bázové funkce normovány (tzn. pokud se jejich funkční hodnota mění od nuly do jedničky), pak koeficienty cn mají význam uzlových hodnot (vzorků) hledané funkce f (obr. 4.1A.3). Aproximaci založenou na lokálních bázových funkcích nazýváme vícebázovou aproximací.


III. Váhové funkce

Mezi nejčastěji používané přístupy k minimalizaci rezidua patří kolokační metoda a metoda Galerkinova.

Kolokační metoda využívá k váhování Diracovy impulsy, umístěné do bodů, v nichž počítáme hodnoty neznámého proudového rozložení

wm(z)=δ(zzm). ( 4.1A.8 )

Kolokační metoda vykazuje velmi nízké výpočetní nároky, jelikož díky filtrační vlastnosti Diracových impulsů je jedna integrace eliminována

n=1Ncnabfn(zm,ξ)dξ=g(zm). ( 4.1A.9 )

Na druhou stranu je minimalizace rezidua vztažena pouze k bodům zm, do nichž byly umístěny váhovací impulsy.

Galerkinova metoda využívá k váhování funkce, které jsou identické s funkcemi bázovými

wm(z)=fm(z). ( 4.1A.10 )

Galerkinova metoda vykazuje ve srovnání s metodou kolokační vyšší výpočetní nároky protože v jejím případě k eliminaci jednoho integrování nedochází. Na druhou stranu jsou však do procesu minimalizace rezidua zahrnuty všechny body analyzované oblasti z ∈ <a,b>.


IV. Drátové antény

Uvažujme drátovou anténu z obr. 4.1A.1. Potom můžeme vyzařované elektromagnetické pole popsat pomocí vektorového potenciálu A (popisuje působení proudů na anténě) a skalárního potenciálu φ [2]

2Az(z)z2+k2Az(z)=μ0Jz(z), ( 4.1A.11a )
2φ(z)z2+k2φ(z)=ρ(z)ε0. ( 4.1A.11b )

Zde Jz značí z-ovou složku proudové hustoty [A.m-2], vnucenou anténě zdrojem, ρ je objemová hustota náboje [C.m-3] na anténním vodiči, Az značí z-ovou složku vektorového potenciálu, φ je skalární potenciál, k=2π/λ značí vlnové číslo a λ vlnovou délku.

Elektrony, které přitékají do antény jako proud, se hromadí na anténním vodiči jako náboj. V druhé polovině periody se směr proudu otočí a náboje z konců anténního vodiče odtékají zpět do zdroje. Jelikož náboje a proudy na anténě spolu souvisejí, musíme vzájemně je svázat. Činíme tak podmínkou kontinuity [2]

Jz(z)z+jωρ(z)=0. ( 4.1A.12a )

Pokud je poloměr anténního vodiče mnohem menší než vlnová délka a << λ, potom můžeme předpokládat, že proudy a náboje jsou soustředěny na ose vodiče. Tento předpoklad je samozřejmě nesprávný (v důsledku povrchového jevu jsou náboje a proudy soustředěny na povrchu vodiče), avšak metoda i přes tento nesprávný předpoklad dává překvapivě dostatečně přesné výsledky [5]. Řešíme-li (s uvážením uvedeného chybného předpokladu) soustavu (4.1A.11), dostáváme [2]

Az(z)=μ4π2hIz(ξ)exp[ jkR(z,ξ) ]R(z,ξ)dξ, ( 4.1A.12b )
φ(z)=14πε2hσ(ξ)exp[ jkR(z,ξ) ]R(z,ξ)dξ. ( 4.1A.12c )

Zde Iz(ξ) značí proud [A] tekoucí v ose anténního vodiče, σ(ξ) je délková hustota náboje [C.m-1] na ose anténního vodiče, R(z,ξ) je vzdálenost mezi pozicí ξ zdroje pole Iz(ξ) a σ(ξ), dále z je místo, v němž počítáme potenciály A(z) a φ(z). Známe-li potenciály A(z) a φ(z), můžeme vypočíst intenzitu vyzařovaného elektrického pole [2]

Ezs(z)=jωAz(z)φ(z)z. ( 4.1A.12d )

Elektrická intenzita musí splňovat okrajovou podmínku na povrchu dokonale elektricky vodivého anténního vodiče S

Ezi+Ezs=0 na S . ( 4.1A.12e )

Ezi značí z-ovou složku (tj. složku tečnou k povrchu anténního vodiče) vektoru elektrické intenzity dopadající vlny. Dopadající vlna je vytvořena zdroji mimo vlastní anténu. Když analyzujeme anténu jako vysílací, je Ezi intenzita vytvořená napájecím zdrojem v budicí štěrbině (Ez v obr. 4.1A.2). Když analyzujeme anténu jako přijímací, je Ezi intenzita přijímaného vlnění (po celé délce vodiče).

Chceme-li určit rozložení proudu na anténě, musíme vyřešit soustavu (4.1A.12).

Abychom se mohli postarat o splnění okrajové podmínky (4.1A.12e), musíme počítat elektrickou intenzitu (a tudíž i potenciály A a φ) na povrchu anténního vodiče. Proto můžeme vzdálenost mezi zdroji pole (na ose) a místy pozorování (na povrchu vodiče) vyjádřit jako

R(z,ξ)=a2+(zξ)2. ( 4.1A.13 )

V dalších odstavcích budeme předpokládat po částech konstantní bázové funkce a Diracovy funkce váhové. S využitím těchto funkcí budeme hledat řešení soustavy (4.1A.12).

Obr. 4.1A.4
Obr. 4.1A.4Po částech konstatní aproximace

V prvém kroku musíme analyzovanou anténu diskretizovat. Tato diskretizace je naznačena na obr. 4.1A.4. Dolní hranice segmentů je označena indexem ”-”, horní hranice indexem ”+”. Dolní hranice prvního segmentu a horní hranice posledního segmentu jsou posunuty o polovinu segmentu za konec anténního vodiče, aby bylo možno modelovat uzel proudu I(-h)=I(h)=0 na koncích antény. Délka všech segmentů je stejná Δ = 2α.

Dosadíme-li po částech konstantní aproximaci do (4.1A.12b,c), dostáváme

Az(z)μ4πn=1NInh+(n0,5)Δh+(n+0,5)Δexp[ jkR(z,ξ) ]R(z,ξ)dξ, ( 4.1A.14b )
φ(z)14πεn=1Nσnh+(n0,5)Δh+(n+0,5)Δexp[ jkR(z,ξ) ]R(z,ξ)dξ. ( 4.1A.14c )

Ve výše uvedených vztazích jsou In a σn uzlové hodnoty proudu a nábojové hustoty.

Jelikož první derivace po částech konstantní aproximace je nulová na konstantní části funkce a neexistuje na hranicích, derivace v (4.1A.12a) a (4.1A.12d) jsou nahrazeny konečnými diferencemi. Uvážíme-li že In = Iz (-h+nΔ), můžeme rovnici kontinuity přepsat do tvaru

Iz(h+(n+1)Δ)Iz(h+nΔ)Δ+jωσ(h+(n+0,5)Δ)0 ( 4.1A.15a )

a vztah pro výpočet intenzity elektrického pole přechází na

Ezs(h+nΔ)jωAz(h+nΔ)φ[ h+(n+0,5)Δ ]φ[ h+(n0,5)Δ ]Δ. ( 4.1A.15d )

Vztahy (4.1A.15a) a (4.1A.15d) odpovídají skutečnosti, že Diracovy impulsy jsou při váhování umístěny do středu segmentu u vektorového potenciálu

Az(h+mΔ)μ4πn=1NInh+(n0,5)Δh+(n+0,5)Δexp[ jkR(h+mΔ,ξ) ]R(h+mΔ,ξ)dξ ( 4.1A.15b )

a do krajů segmentů u potenciálu skalárního

φ[ h+(m+0,5)Δ ]14πεn=1Nσn+h+nΔh+(n+1)Δexp{ jkR[ h+(m+0,5)Δ,ξ ] }R[ h+(m+0,5)Δ,ξ ]dξ. ( 4.1A.15c )

Ve vztahu (4.1A.15c), σ n+ = σ [-h+(n+0.5)Δ].

Vztah, (4.1A.15) může být přepsán do kompaktnějšího tvaru

σn+1jω[ In+1InΔ ], ( 4.1A.16a )
Az(m)μ4πn=1NInΔnexp[ jkR(m,ξ) ]R(m,ξ)dξ, ( 4.1A.16b )
φ(m+)14πεn=1Nσn+Δn+exp[ jkR(m+,ξ) ]R(m+,ξ)dξ, ( 4.1A.16c )
Ezi(m)jωAz(m)φ(m+)φ(m)Δ. ( 4.1A.16d )

Při odvození (4.1A.16d) byla uvážena okrajová podmínka (4.1A.12e).

Nyní se podrobněji podívejme na rovnici kontinuity (4.1A.16a). Tato podmínka vyjadřuje skutečnost, že jednotlivé segmenty antény mohou být nahrazeny elementárními elektrickými dipóly (obr. 4.1A.5). Uvážíme-li tento fakt, můžeme vyjádřit příspěvek n-tého segmentu (elementárního dipólu) k hodnotě skalárního potenciálu na pravé hranici m-tého segmentu díky (4.1A.16c) jako

φ(m+)=1jωε[ InΔn+exp(jkR)4πRdξInΔnexp(jkR)4πRdξ ]1Δ. ( 4.1A.17 )

Dosazením (4.1A.17) a (4.1A.16b) do (4.1A.16d) a vynásobením obou stran rovnice délkou segmentu Δ dostáváme

EziΔ=ZI. ( 4.1A.18 )

Pro prvky Zmn impedanční matice Z platí:

Zmn=jωμΔΔnexp[ jkR(m,ξ) ]4πR(m,ξ)dξ+
+1jωε[ Δn+exp[ jkR(m+,ξ) ]4πR(m+,ξ)dξΔnexp[ jkR(m+,ξ) ]4πR(m+,ξ)dξ ]1Δ
1jωε[ Δn+exp[ jkR(m,ξ) ]4πR(m,ξ)dξΔnexp[ jkR(m,ξ) ]4πR(m,ξ)dξ ]1Δ,
( 4.1A.19 )
Obr. 4.1A.5
Obr. 4.1A.5Anténa jako soubor elementárních elektrických dipólů

Zmn popisuje příspěvek proudu a náboje na segmentu n k napětí indukovanému na segmentu m.

Jelikož složka elektrické intenzity, tečná k anténnímu vodiči, je nulová na všech segmentech vyjma napájecí štěrbiny, prvky sloupcového vektoru napětí (levá strana rovnice 4.1A.18) jsou nulové vyjma případu napájecí štěrbiny (na štěrbince jsme předpokládali napětí 1 V). Z rovnice (4.1A.18) tedy můžeme vyjádřit sloupcový vektor uzlových hodnot rozložení proudu na anténě I. Poměr vstupního napětí a vstupního proudu je potom roven vstupní impedanci antény.

Jako příklad si uveďme výsledek analýzy symetrického dipólu s délkou ramene h = λ a s poloměrem anténního vodiče a = 0.001588 λ . Rozložení proudu na anténě, získané pomocí popsané metody, je nakresleno na obr. 4.1A.6.

Ve vrstvě C uvádíme uživatelský popis programu, s jehož pomocí je možno dosáhnout uvedeného výsledku. Ve vrstvě D pak čtenář nalezne informace o tom, jak je možno program efektivně sestavit v Matlabu.




Obr. 4.1A.6
Obr. 4.1A.6Rozložení proudu na symetrickém dipólu. Po částech konstantní aproximace, Diracovo vážení. Délka dipólu 2λ, průměr 0,001588λ, 64 segmentů.

Copyright © 2010 FEEC VUT Brno All rights reserved.