la propagation_OEM-guide d'ondes
la propagation_OEM-guide d'ondes
la propagation_OEM-guide d'ondes
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
D2 Propagation d’une <strong>OEM</strong> dans un <strong>guide</strong> d’onde<br />
Dans le domaine des hyperfréquences (de 10 8 à 10 12 Hz), les courants deviennent<br />
superficiels et les pertes par effet Joule dans une ligne de transmission deviennent<br />
trop importantes. On les supprime donc ainsi que le diélectrique central et on ne<br />
conserve que le conducteur extérieur pour forme un « tube », à section rectangu<strong>la</strong>ire,<br />
circu<strong>la</strong>ire ou elliptique, qui va servir à transporter l’énergie de ces ondes<br />
hyperfréquences entre les émetteurs et les antennes.<br />
A l'intérieur de ce tube à parois conductrices le champ électrique E r et l’excitation<br />
magnétique H r<br />
de l'onde électromagnétique qui s’y propage doivent satisfaire à<br />
certaines conditions aux limites: <strong>la</strong> composante tangentielle de E r et <strong>la</strong> composante<br />
normale de H r sont nulles à <strong>la</strong> surface du conducteur. Pour étudier <strong>la</strong> <strong>propagation</strong> de<br />
l’onde dans le <strong>guide</strong>, il faut donc trouver <strong>la</strong> combinaison de E r et de H r satisfaisant<br />
d'une part aux conditions aux limites et d'autre part aux équations de Maxwell.<br />
Dans <strong>la</strong> suite de cet exposé, on va étudier <strong>la</strong> <strong>propagation</strong> d’une onde<br />
électromagnétique dans un <strong>guide</strong> rectangu<strong>la</strong>ire (Figure 36) et on assimilera l’air au vide. De plus, on supposera que les parois<br />
sont parfaitement conductrices (conductivité infinie) de sorte que l’onde n’y pénètre pas et n’y dissipe aucune énergie. Le <strong>guide</strong><br />
est alors considéré sans pertes.<br />
D2.1. Les équations de <strong>propagation</strong> dans le <strong>guide</strong><br />
D.2.1.1. Les équations de Maxwell dans le <strong>guide</strong><br />
Puisque l’air contenu dans le <strong>guide</strong> est assimilé au vide, on peut utiliser les équations de Maxwell dans le vide.<br />
div E = 0<br />
Eq.100<br />
∂H<br />
rot E = −μ 0 Eq.101<br />
∂t<br />
div H = 0<br />
Eq.102<br />
∂E<br />
rot H = ε 0 Eq.103<br />
∂t<br />
D.2.1.2. Les équations de <strong>propagation</strong><br />
Ce sont les mêmes que celle obtenues pour <strong>la</strong> <strong>propagation</strong> dans le vide.<br />
2<br />
∂ E<br />
ΔE − μ0 ε0<br />
= 0<br />
Eq.104<br />
2<br />
∂t<br />
2<br />
∂ H<br />
ΔH − μ0 ε0<br />
= 0<br />
Eq.105<br />
2<br />
∂t<br />
Les grandeurs E r et de H r sont des fonctions qui dépendent du temps et des trois variables d’espace x, y et z (en coordonnées<br />
cartésiennes). Si on désigne par Ox l’axe du <strong>guide</strong>, qui est aussi <strong>la</strong> direction de <strong>propagation</strong> de l’onde dans le <strong>guide</strong>, on fait jouer<br />
un rôle différent à cette variable. Aussi, on va écrire l’opérateur Lap<strong>la</strong>cien comme <strong>la</strong> somme du Lap<strong>la</strong>cien longitudinal<br />
2<br />
2 2<br />
∂<br />
∂ ∂<br />
Δ x = et du Lap<strong>la</strong>cien transversal. Δ<br />
2<br />
⊥ = + .<br />
2 2<br />
∂x<br />
∂y ∂z<br />
Les deux équations de <strong>propagation</strong> dans le <strong>guide</strong> s’écrivent alors :<br />
Δ ⊥<br />
2<br />
x =<br />
∂ E<br />
E + Δ E − μ0<br />
ε0<br />
2<br />
∂t<br />
2<br />
0<br />
Eq.106<br />
∂ H<br />
Δ ⊥ H + Δx H − μ0<br />
ε0<br />
= 0<br />
Eq.107<br />
2<br />
∂t<br />
D.2.1.3. Les équations de <strong>propagation</strong> pour une onde sinusoïdale<br />
Le <strong>guide</strong> est supposé illimité suivant Ox (ou, ce qui est physiquement équivalent, fermé sur son impédance itérative), <strong>la</strong><br />
dépendance en x est contenue dans l’argument de l’exponentielle quand on utilise <strong>la</strong> notation complexe. Ainsi, les expressions<br />
générales de E r et de H r sont :<br />
E(x, y,z,t)<br />
Figure 36: Schéma d’un <strong>guide</strong><br />
d’ondes rectangu<strong>la</strong>ire<br />
jωt−K.x<br />
= E0(y,z)<br />
e<br />
Eq.108<br />
H(x, y,z, t)<br />
jωt<br />
−K.x<br />
0(y,z)<br />
e<br />
= H<br />
Eq.109<br />
41
avec E 0 (y,z) = E0x<br />
ex<br />
+ E0y<br />
ey<br />
+ E0z<br />
e z , H 0 (y,z) = H0x<br />
ex<br />
+ H0y<br />
ey<br />
+ H0z<br />
e z et K=α+jk g (α est une constante<br />
d’atténuation et k g le vecteur d’onde dans le <strong>guide</strong>)<br />
2<br />
2<br />
∂<br />
On en déduit<br />
E K<br />
2 ∂ E 2<br />
= E et = −ω E .<br />
2<br />
2<br />
∂x<br />
∂t<br />
Quand <strong>la</strong> <strong>propagation</strong> se faisait dans l’air assimilé au vide, le vecteur d’onde vaut alors μ 0 ε 0 ω. Pour différencier le vecteur<br />
d’onde dans le <strong>guide</strong> (plus précisément dans le milieu limité par le <strong>guide</strong>) du vecteur d’onde dans le même milieu illimité, nous<br />
poserons :<br />
2<br />
km<br />
0<br />
0<br />
2<br />
= μ ε ω<br />
Eq.110<br />
De sorte que les équations de <strong>propagation</strong> pour une onde sinusoïdale prennent <strong>la</strong> forme suivante :<br />
Δ ⊥<br />
E + (K<br />
2<br />
+ k<br />
2<br />
m<br />
)E = 0<br />
Eq.111<br />
2<br />
2<br />
Δ ⊥ H + (K + km)H<br />
= 0<br />
Eq.112<br />
Pour trouver <strong>la</strong> structure de l’onde, il serait nécessaire de résoudre les équations de <strong>propagation</strong> pour chacune des composantes<br />
de E r et de H r en tenant compte des conditions aux limites imposées par les parois du <strong>guide</strong>.<br />
Mais les composantes ne sont pas totalement indépendantes entre elles puisqu’elles sont liées par les équations de Maxwell. Pour<br />
réduire le nombre de variables, on choisit généralement de faire jouer un rôle particulier à E0x<br />
et H 0x<br />
, les composantes de E r<br />
et de H r dans <strong>la</strong> direction de <strong>propagation</strong>. On utilise alors les équations de Maxwell pour trouver les autres composantes en<br />
fonction de E0x<br />
et H 0x<br />
.<br />
D.2.1.4. Re<strong>la</strong>tions entre les composantes<br />
Il faut exprimer E 0y<br />
, E 0z<br />
, H 0y<br />
et H0z<br />
en fonction de E0x<br />
et H 0x<br />
.<br />
∂H<br />
Les deux équations de Maxwell qui couplent les champs sont rot E = −μ 0 et<br />
∂t<br />
⎧∂E<br />
∂E<br />
z y ∂Hx<br />
⎪ − = −μ0<br />
⎪<br />
∂y<br />
∂z<br />
∂t<br />
∂H<br />
⎪∂E<br />
∂ ∂H<br />
x Ez<br />
y<br />
rot E = −μ0<br />
⎨ − = −μ0<br />
∂t<br />
⎪ ∂z<br />
∂x<br />
∂t<br />
⎪∂Ey<br />
∂Ex<br />
∂Hz<br />
⎪ − − = −μ0<br />
⎩<br />
∂x<br />
∂y<br />
∂t<br />
Ce qui donne trois équations:<br />
∂E<br />
∂E<br />
0z<br />
0y<br />
− = −jμ0ω<br />
H0x<br />
∂y<br />
∂z<br />
∂E<br />
0 x<br />
+ KE0z<br />
= −jμ0ω<br />
H 0y<br />
∂z<br />
∂E0x<br />
− KE0y<br />
− = −jμ0ω<br />
H0z<br />
∂y<br />
De <strong>la</strong> même façon :<br />
⎧∂H<br />
∂H<br />
z y ∂E<br />
x<br />
⎪ − = ε0<br />
⎪ ∂y<br />
∂z<br />
∂t<br />
∂E<br />
⎪∂H<br />
∂ ∂E<br />
x Hz<br />
y<br />
rot H = ε0 ⎨ − = ε0<br />
qui permet d’obtenir :<br />
∂t<br />
⎪ ∂z<br />
∂x<br />
∂t<br />
⎪∂H<br />
y ∂H<br />
x ∂E<br />
z<br />
⎪ − − = ε0<br />
⎪⎩<br />
∂x<br />
∂y<br />
∂t<br />
rot H<br />
∂E<br />
= ε 0 .<br />
∂t<br />
Eq.113<br />
Eq.114<br />
Eq.115<br />
42
∂H<br />
∂y<br />
∂H<br />
∂z<br />
− KH<br />
∂H<br />
0z<br />
0y<br />
− = jε0<br />
ω E0x<br />
∂z<br />
0 x<br />
+ KH0z<br />
= jε0<br />
ω E 0y<br />
∂H<br />
∂y<br />
0x<br />
0y<br />
− = jε0ω<br />
E0z<br />
Eq.116<br />
Eq.117<br />
Eq.118<br />
La combinaison des équations précédentes permet d’écrire l’ensemble des composantes des champs en fonction E0x<br />
et H 0x<br />
:<br />
Par exemple, <strong>la</strong> combinaison de l’Eq.117et de l’Eq.115 donne :<br />
∂H0x<br />
∂H0x<br />
K ∂E<br />
x<br />
j E KH<br />
(KE0y<br />
0<br />
2<br />
2<br />
∂H0x<br />
∂E0x<br />
ε 0 ω 0y = + 0z = +<br />
+ ) et − ( K + ε0<br />
μ0<br />
ω ) E0y<br />
= jμ0ω<br />
+ K<br />
∂z<br />
∂z<br />
jμ0ω<br />
∂y<br />
∂z<br />
∂y<br />
Et finalement :<br />
−1<br />
⎡ ∂H<br />
⎤<br />
0x ∂E0x<br />
E0y<br />
= ⎢jμ0ω<br />
+ K ⎥<br />
Eq.119<br />
2 2<br />
(K + k m)<br />
⎣ ∂z<br />
∂y<br />
⎦<br />
Pour les autres composantes, les expressions sont :<br />
1 ⎡ ∂H<br />
⎤<br />
0x ∂E0x<br />
E0z<br />
= ⎢jμ<br />
ω − K<br />
2 2 0<br />
⎥<br />
Eq.120<br />
(K + km)<br />
⎣ ∂y<br />
∂z<br />
⎦<br />
1 ⎡ ∂E<br />
⎤<br />
0x ∂H0x<br />
H0y<br />
= ⎢jε<br />
ω − K<br />
2 2 0<br />
⎥<br />
Eq.121<br />
(K + km)<br />
⎢⎣<br />
∂z<br />
∂y<br />
⎥⎦<br />
−1<br />
⎡ ∂E<br />
⎤<br />
0x ∂H0x<br />
H0z<br />
= ⎢jε<br />
ω + K<br />
2 2 0<br />
⎥<br />
(K + km)<br />
⎣ ∂y<br />
∂z<br />
⎦<br />
Eq.122<br />
Les équations précédentes, ajoutées aux conditions aux limites et aux équations de <strong>propagation</strong> pour les composantes et<br />
H<br />
0x<br />
permettent en principe de trouver <strong>la</strong> structure de l’onde qui se propage dans le <strong>guide</strong>.<br />
Le paragraphe suivant rappelle les conditions aux limites pour le passage d’une onde électromagnétique du vide à un conducteur.<br />
D.2.1.5.<br />
Conditions aux limites<br />
Les conditions générales sont:<br />
∧ (E − E ) 0<br />
n12 2 1 =<br />
n<br />
12 .(D2<br />
− D1)<br />
= σ<br />
n =<br />
12 ∧ (H2<br />
− H1)<br />
jS<br />
n 1<br />
12 .(B2<br />
− B ) = 0<br />
s<br />
E0x<br />
Appliquées à un <strong>guide</strong> rectangu<strong>la</strong>ire de dimensions a et b<br />
(Figure 37), on obtient des conditions pour certaines des<br />
composantes. Comme on ne connaît pas les charges et les<br />
courants surfaciques, on utilise les composantes tangentielles<br />
de E r ainsi que les composantes normales de H r .<br />
Pour les composantes tangentielles de E r , on obtient :<br />
E x (0,z) 0 E x (a,z) E x (y,0) E x (y,b)<br />
0 = , 0 = 0 , 0 = 0 , 0 = 0 , E 0 y (y,0) = 0 , E 0 y (y,b) = 0 , E 0 z (0,z) = 0 , E 0 z (a,z) = 0<br />
Pour les composantes normales de H r , on obtient :<br />
H 0 y (0,z) = 0 , H 0 y (a,z) = 0 , H 0 z (y,0) = 0 , H 0 z (y,b) = 0<br />
Figure 37: Composantes des champs E r et H r à considérer<br />
pour les conditions aux limites.<br />
43
D2.2. La re<strong>la</strong>tion de dispersion dans le <strong>guide</strong><br />
L’équation de <strong>propagation</strong> pour <strong>la</strong> composante selon Ox du champ électrique est :<br />
2 2<br />
Δ ⊥ E + (K + k )E 0 .<br />
0 x<br />
m 0x<br />
=<br />
Ce qui s’écrit encore :<br />
2<br />
∂ E<br />
∂y<br />
0x<br />
2<br />
2<br />
∂ E<br />
+<br />
∂z<br />
0x<br />
2<br />
= −(K<br />
2<br />
+ k<br />
2<br />
m<br />
) E<br />
0x<br />
Eq.123<br />
Remarque : Cette équation peut aussi être obtenue en utilisant l’ensemble des équations liant les composantes des champs entre<br />
elles et en écrivant que div E = 0 .<br />
L’Eq.123 ne fait intervenir que les variables d’espace y et z. Sa résolution, en utilisant les conditions aux limites, permet de<br />
trouver (K<br />
2 2<br />
+ km)<br />
, une grandeur positive qui ne dépend que de <strong>la</strong> géométrie du <strong>guide</strong> et est indépendante de <strong>la</strong> pulsation de<br />
l’onde qui s’y propage.<br />
On pose :<br />
2<br />
c<br />
2<br />
2<br />
m<br />
k = K + k<br />
Eq.124<br />
L’équation de dispersion peut alors s’écrire :<br />
Où<br />
2<br />
k<br />
m<br />
et<br />
K<br />
2<br />
2<br />
c<br />
2<br />
k<br />
c<br />
sont deux constantes positives.<br />
2<br />
m<br />
= k − k<br />
Eq.125<br />
K 2 dépend à <strong>la</strong> fois de <strong>la</strong> pulsation et de <strong>la</strong> géométrie du <strong>guide</strong>, par l’intermédiaire de<br />
pour <strong>la</strong> géométrie.<br />
D2.3. Les conditions de <strong>la</strong> <strong>propagation</strong> guidée<br />
D.2.3.1. <strong>la</strong> fréquence de coupure<br />
2 2<br />
k<br />
m<br />
et k<br />
c<br />
étant deux constantes positives, K 2 est un réel qui est soit positif, soit négatif.<br />
Puisque K=α+jk g :<br />
2<br />
k<br />
m<br />
2<br />
= μ0<br />
ε0ω<br />
pour <strong>la</strong> pulsation et<br />
• K 2 >0 implique K=α. L’onde ne se propage pas. Dans ces conditions, une onde introduite dans le <strong>guide</strong> s’amortit,<br />
même dans un <strong>guide</strong> sans pertes. Quand le milieu est l’air assimilé au vide, le coefficient d’atténuation est :<br />
2<br />
2 2 ω<br />
α = kc<br />
−<br />
Eq.126<br />
2<br />
c<br />
• K 2 ω c avec :<br />
c<br />
2<br />
2 ωc<br />
kc<br />
= Eq.128<br />
2<br />
c<br />
2<br />
2<br />
k<br />
c<br />
44
Un <strong>guide</strong> d’onde se comporte comme un filtre passe-haut.<br />
Pour qu’une onde puisse s’y propager, sa fréquence doit être<br />
ωc<br />
supérieure à <strong>la</strong> fréquence de coupure fc = . 2 π<br />
La courbe de dispersion kg(ω) est représentée sur <strong>la</strong> Figure<br />
38. Le vecteur d’onde dans le <strong>guide</strong> k g est une fonction<br />
croissante de <strong>la</strong> pulsation de l’onde. Il tend vers c<br />
ω quand <strong>la</strong><br />
pulsation est très grande devant <strong>la</strong> pulsation de coupure.<br />
k g<br />
ωω c<br />
<strong>propagation</strong><br />
D.2.3.2. <strong>la</strong> re<strong>la</strong>tion fondamentale de <strong>la</strong><br />
<strong>propagation</strong> guidée<br />
Elle est obtenue à partir de l’équation de dispersion en<br />
introduisant les trois longueurs d’onde suivantes :<br />
2π<br />
• λ m , <strong>la</strong> longueur d’onde dans le milieu : km<br />
=<br />
λ<br />
• λ g , <strong>la</strong> longueur d’onde dans le <strong>guide</strong> : k<br />
g<br />
2π<br />
=<br />
λ<br />
g<br />
m<br />
Figure 38: Re<strong>la</strong>tion de dispersion dans un <strong>guide</strong>. La<br />
courbe en pointillé représente k g =ω/c.<br />
λ g<br />
ω c<br />
ω>ω c<br />
<strong>propagation</strong><br />
ω
2 ω − ωc<br />
2<br />
La vitesse de groupe est obtenue en différenciant kg<br />
= : 2c<br />
k dk = 2ω<br />
dω<br />
2<br />
g g .<br />
c<br />
2<br />
dω<br />
2<br />
kg<br />
c<br />
Ce qui donne pour expression de <strong>la</strong> vitesse de groupe = c = .<br />
dk ω v<br />
v<br />
g c 1<br />
2<br />
c<br />
2<br />
2<br />
g<br />
2<br />
ϕ<br />
ω<br />
= −<br />
Eq.131<br />
ω<br />
La vitesse de groupe est plus faible que <strong>la</strong> vitesse de l’onde dans le milieu illimité et le produit v g .v ϕ est égal à c 2 .<br />
Les évolutions de <strong>la</strong> vitesse de phase et de <strong>la</strong> vitesse de groupe<br />
en fonction de <strong>la</strong> pulsation de l’onde sont représentées sur <strong>la</strong><br />
Figure 40.<br />
Quand <strong>la</strong> fréquence de l’onde est très grande devant <strong>la</strong><br />
fréquence de coupure, <strong>la</strong> vitesse de phase et <strong>la</strong> vitesse de<br />
groupe tendent vers une même valeur, <strong>la</strong> vitesse de l’onde<br />
dans le milieu illimité.<br />
Pour une fréquence élevée devant <strong>la</strong> fréquence de coupure, le<br />
<strong>guide</strong> n’est plus dispersif.<br />
v ϕ<br />
v g<br />
c<br />
ωω c<br />
<strong>propagation</strong><br />
D.2.3.4.<br />
Les différents modes de <strong>propagation</strong><br />
Si dans un câble coaxial, on peut faire se propager des ondes Figure 40: Vitesse de phase et vitesse de groupe dans le<br />
transverses électromagnétiques (TEM : le champ électrique E <strong>guide</strong> en fonction de <strong>la</strong> pulsation de l’onde<br />
et l’excitation magnétique H sont perpendicu<strong>la</strong>ires à <strong>la</strong><br />
direction de <strong>propagation</strong>), il n’en est pas de même dans un<br />
<strong>guide</strong>.<br />
En effet, en utilisant E 0x<br />
=0 et H 0x<br />
=0 dans le système d’équations qui relie les composantes entre elles, on trouve que toutes<br />
les composantes s’annulent sauf si on pouvait avoir k g =k m , ce qui n’est possible que dans un milieu illimité et pas dans un <strong>guide</strong><br />
où k g
2 2<br />
⎧∂<br />
E0x<br />
∂ E0x<br />
2<br />
⎪ + = −k<br />
E<br />
2 2 c 0x<br />
⎪ ∂y<br />
∂z<br />
⎪ 2<br />
2<br />
∂ H0x<br />
∂ H0x<br />
2<br />
⎪ + = −k<br />
H<br />
2<br />
2 c 0x<br />
⎪ ∂y<br />
∂z<br />
⎪<br />
⎪ − j ⎡ ∂H0x<br />
∂E<br />
E0y<br />
=<br />
⎪ ⎢μ<br />
ω + k<br />
2 0<br />
g<br />
k<br />
c ⎣ ∂z<br />
∂y<br />
⎨<br />
⎪ j ⎡ ∂H0x<br />
∂E<br />
E0z<br />
=<br />
⎪<br />
⎢μ<br />
ω − k<br />
2 0<br />
g<br />
kc<br />
⎣ ∂y<br />
∂z<br />
⎪<br />
⎪ j ⎡ ∂E0x<br />
∂H<br />
⎪H0y<br />
= ⎢ε<br />
ω − k<br />
2 0<br />
g<br />
⎪ kc<br />
⎣ ∂z<br />
∂y<br />
⎪<br />
⎪ − j ⎡ ∂E0x<br />
∂H<br />
H0z<br />
=<br />
⎪<br />
⎢ε0ω<br />
+ k<br />
2<br />
g<br />
⎩ kc<br />
⎣ ∂y<br />
∂z<br />
0x<br />
0x<br />
0x<br />
0x<br />
⎤<br />
⎥<br />
⎦<br />
⎤<br />
⎥<br />
⎦<br />
⎤<br />
⎥<br />
⎦<br />
⎤<br />
⎥<br />
⎦<br />
Eq.132<br />
Dans le cas d’une onde TE ( E =0), le système d’équation précédent devient :<br />
2<br />
2<br />
0x<br />
∂ H0x<br />
∂ H0x<br />
+ = −k<br />
2<br />
2<br />
∂y<br />
∂z<br />
− jkg<br />
∂H0x<br />
H0y<br />
=<br />
2<br />
k ∂y<br />
c<br />
2<br />
c<br />
H<br />
0x<br />
E<br />
H<br />
0y<br />
0z<br />
− jμ0ω ∂H0<br />
=<br />
2<br />
k ∂z<br />
c<br />
x<br />
E<br />
0z<br />
jμ0ω ∂H0<br />
=<br />
2<br />
k ∂y<br />
c<br />
− jkg<br />
∂H<br />
0 x<br />
= Eq.133<br />
2<br />
k ∂z<br />
c<br />
x<br />
D.2.4.1.<br />
Séparation des variables<br />
H 0x est une fonction ψ des variables y et z. Compte tenu de <strong>la</strong> géométrie du <strong>guide</strong>, on peut considérer que H 0x<br />
dépend<br />
indépendamment de y et z, c'est-à-dire que <strong>la</strong> dépendance en y (en z) est <strong>la</strong> même quel que soit z (ou y). Ceci est une bonne<br />
approximation sauf lorsque on approche les coins du <strong>guide</strong>.<br />
Avec ψ(y,z)=f(y).g(z), on peut trouver une nouvelle expression du Lap<strong>la</strong>cien transversal :<br />
2<br />
2<br />
2<br />
2<br />
∂ Ψ( y,z) ∂ Ψ(y,z)<br />
∂ f (y) ∂ g(z) 2<br />
+ = g(z) + f (y) = −kc<br />
f (y).g(z)<br />
2<br />
2<br />
2<br />
2<br />
∂y<br />
∂z<br />
∂y<br />
∂z<br />
En divisant par ψ(y,z), on obtient :<br />
2<br />
2<br />
1 ∂ f (y) 1 ∂ g(z) 2<br />
+<br />
= −k<br />
2<br />
2 c<br />
f (y) ∂y<br />
g(z) ∂z<br />
Cette somme d’une fonction de y et d’une fonction de z est une constante quels que soient y et z. C’est donc que chacun des<br />
termes est une constante.<br />
2 2 2<br />
En posant k = k + k , on peut écrire deux équations indépendantes :<br />
1<br />
f (y)<br />
2<br />
c<br />
cy<br />
∂ f (y)<br />
= −k<br />
2<br />
∂y<br />
2<br />
cy<br />
et<br />
cz<br />
1<br />
g(z)<br />
2<br />
∂ g(z)<br />
= −k<br />
2<br />
∂z<br />
2<br />
cz<br />
dont les solutions sont :<br />
0<br />
jk cy y<br />
0<br />
− jk cy y<br />
f (y) = A e + B e<br />
Eq.134<br />
0<br />
jk czz<br />
0<br />
− jk czz<br />
g(z) = C e + D e<br />
Eq.135<br />
D.2.4.2. Les vecteurs d’onde de coupure<br />
2 2<br />
Les constantes A0, B0, C0, D0 ainsi k<br />
cy<br />
et k<br />
cz<br />
peuvent être déterminées à partir des conditions aux limites. Celles-ci ont déjà<br />
été précisées lors de l’étude générale. Pour une onde TE :<br />
E 0 y (y,0) = 0<br />
E 0 y (y,b) = 0<br />
E 0 z (0,z) = 0<br />
E 0 z (a,z) = 0<br />
H 0 y (0,z) = 0<br />
H 0 y (a,z) = 0<br />
H 0 z (y,0) = 0<br />
H 0 z (y,b) = 0<br />
47
On va utiliser les conditions aux limites pour H 0 y .<br />
jk<br />
= −<br />
2<br />
k<br />
∂H<br />
∂y<br />
jk<br />
g 0x g<br />
H0y<br />
= −<br />
2<br />
c kc<br />
df(y)<br />
g(z)<br />
dy<br />
df (y) ⎞ df (y)<br />
−k<br />
cyy<br />
k cyy<br />
df (y) ⎞<br />
H 0 y (0,z) = 0 impose alors ⎟ = 0 , = jkcy<br />
A0<br />
e − jkcy<br />
B0<br />
e et ⎟ = jkcy(A0<br />
− B0)<br />
dy ⎠ dy<br />
dy<br />
y=<br />
0<br />
⎠y=<br />
0<br />
Les deux constantes A 0 et B 0 sont égales.<br />
La même condition en y=a donnera :<br />
df (y) ⎞<br />
k cya<br />
−k<br />
cya<br />
⎟ = jkcy<br />
A0<br />
(e + e ) = −2kcy<br />
A0<br />
sin(kcya)<br />
= 0<br />
dy ⎠y=<br />
a<br />
Ce qui permet d’obtenir <strong>la</strong> constante k cy .<br />
mπ<br />
k cy<br />
= Eq.136<br />
a<br />
De <strong>la</strong> même façon, on trouverait :<br />
k cz<br />
nπ<br />
= Eq.137<br />
b<br />
Pour chacun des couples de valeurs (m, n), il existe un vecteur d’onde de coupure qui s’écrit :<br />
mπ<br />
2 nπ<br />
2<br />
kc (m,n) = ( ) + ( )<br />
Eq.138<br />
a b<br />
Il lui correspond une fréquence de coupure ω c(m,n) au dessus de <strong>la</strong>quelle une onde TE peut se propager suivant le mode TE mn .<br />
D.2.4.3. Les longueurs d’onde et les fréquences de coupure<br />
La fréquence de coupure est liée au vecteur d’onde de coupure par <strong>la</strong> re<strong>la</strong>tion<br />
On en déduit l’expression des longueurs d’onde de coupure λ c(m,n) .<br />
2<br />
λc(m,n)<br />
=<br />
m 2 n 2<br />
( ) + ( )<br />
a b<br />
k<br />
c<br />
2π<br />
=<br />
λ<br />
c<br />
Eq.139<br />
On en déduit les fréquences de coupure pour les modes TE mn :<br />
c m 2 n 2<br />
f c (m,n) = ( ) + ( )<br />
Eq.140<br />
2 a b<br />
Il peut exister une infinité de modes de <strong>propagation</strong> caractérisés par les nombres entiers m et n. La fréquence de coupure est liée<br />
à ces nombres et, plus l'onde TE mn est d'ordre élevé, plus <strong>la</strong> fréquence de coupure est élevée.<br />
La fréquence de coupure <strong>la</strong> plus faible est celle qui est associée au mode TE 10 si a>b (ou au mode TE 01 si a
Les autres composantes se déduisent aisément de cette expression.<br />
− jμ0ω ∂H0x<br />
jμ0ω<br />
nπ<br />
y z<br />
E0y<br />
= = H0<br />
cos(mπ<br />
).sin(nπ<br />
)<br />
Eq.142<br />
2<br />
2<br />
k ∂z<br />
k b<br />
a b<br />
c<br />
c<br />
jμ0ω ∂H0x<br />
− jμ0ω<br />
mπ<br />
y z<br />
= = H sin(mπ<br />
).cos(n )<br />
Eq.143<br />
2<br />
2<br />
k ∂y<br />
k a<br />
a b<br />
E0z<br />
0<br />
π<br />
c<br />
c<br />
− jkg<br />
∂H<br />
jk<br />
0x g mπ<br />
y z<br />
= = H sin(mπ<br />
).cos(n )<br />
Eq.144<br />
2<br />
2<br />
k ∂y<br />
k a<br />
a b<br />
H0y<br />
0<br />
π<br />
c<br />
c<br />
− jk g ∂H<br />
jk<br />
0x g nπ<br />
y z<br />
= = H cos(mπ<br />
).sin(n )<br />
Eq.145<br />
2<br />
2<br />
k ∂z<br />
k b<br />
a b<br />
H0z<br />
0<br />
π<br />
c<br />
c<br />
Les composantes du champ électriques sont toutes nulles si m et n sont nuls.<br />
Par contre, m ou n peut être nul, ce qui donne les ondes TE 0n ou TE m0 .<br />
D.2.4.5.<br />
La structure d’une onde TEmn.<br />
A partir des expressions précédentes, on observe que<br />
dépendance en y et z. Le rapport<br />
E<br />
H<br />
0y<br />
0z<br />
E 0y<br />
et 0z<br />
est égal à l’opposé du rapport<br />
H d’une part ainsi que E 0z<br />
et H 0y<br />
d’autre part ont <strong>la</strong> même<br />
E<br />
H<br />
0z<br />
0y<br />
et vaut<br />
μ0ω<br />
.<br />
k<br />
g<br />
On peut donc écrire ainsi les composantes E et H :<br />
⎧<br />
⎪<br />
⎪0<br />
⎪ jμ0ω<br />
nπ<br />
y<br />
E : H0⎨<br />
cos(mπ<br />
).sin(nπ<br />
2<br />
⎪ kc<br />
b a<br />
⎪−<br />
jμ0ω<br />
mπ<br />
y<br />
⎪ sin(mπ<br />
).cos(nπ<br />
2<br />
⎩ kc<br />
a a<br />
H : H<br />
0<br />
⎧ y<br />
⎪cos(mπ<br />
).cos(nπ<br />
⎪ a<br />
⎪−<br />
kg<br />
⎨ E0z<br />
⎪ μ0ω<br />
⎪ kg<br />
⎪ E0y<br />
⎪⎩<br />
μ0ω<br />
z<br />
)<br />
b<br />
×<br />
e<br />
z<br />
)<br />
b<br />
z<br />
)<br />
b<br />
j( ωt<br />
−k<br />
x)<br />
g<br />
×<br />
e<br />
j( ωt<br />
−k<br />
gx)<br />
Eq.146<br />
Eq.147<br />
Le produit sca<strong>la</strong>ire E . H est donc nul. Les vecteurs E et H sont perpendicu<strong>la</strong>ires. La structure d’une onde TE est représentée<br />
à un instant donné et pour une abscisse donnée sur <strong>la</strong> Figure 42. E est transversal, <strong>la</strong> composante transversale de B est<br />
perpendicu<strong>la</strong>ire à E de même que le champ total.<br />
D.2.4.6. La structure d’une onde TE10.<br />
Pour une onde TE10 (m=1 et n=0), les composantes de l’onde s’écrivent :<br />
⎧ y<br />
⎧<br />
⎪cos(<br />
π )<br />
⎪<br />
a<br />
0<br />
⎪<br />
⎪<br />
⎪kg<br />
y<br />
E0 : H0⎨0<br />
et H0 : H0⎨<br />
sin( π )<br />
Eq.148<br />
⎪−<br />
jμ0ω<br />
y ⎪kc<br />
a<br />
⎪ sin( π ) ⎪0<br />
⎪⎩<br />
kc<br />
a ⎪<br />
⎩<br />
49
Figure 42: Représentation schématique de <strong>la</strong><br />
structure d’une onde TE.<br />
Figure 42: Représentation schématique des lignes de champ<br />
électrique et magnétique pour le mode TE10.<br />
Le champ électrique a une seule composante, suivant OZ, et l’excitation magnétique a une composante longitudinale et une<br />
composante transversale suivant OY et toutes ces composantes sont indépendantes de <strong>la</strong> variable z.<br />
Les lignes de champ sont représentées sur <strong>la</strong> Figure 42.<br />
• Dans un p<strong>la</strong>n YOZ, les lignes de champ électrique sont des segments de droites parallèles à OZ et les lignes de champ<br />
magnétique sont des segments de droites parallèles à OY.<br />
D.2.4.7.<br />
• Dans le p<strong>la</strong>n XOZ, les lignes de champ électrique sont des segments de droite parallèles à OZ qui se reproduisent avec<br />
<strong>la</strong> périodicité λ g .<br />
• Dans le p<strong>la</strong>n XOY, les lignes de champ magnétique sont des ellipses centrées sur l’axe du <strong>guide</strong> réparties avec <strong>la</strong><br />
périodicité λ g . Ces ellipses se déforment quand on approche de <strong>la</strong> paroi.<br />
Impédance itérative du <strong>guide</strong><br />
Dans un <strong>guide</strong> d’onde, l’impédance est le rapport des composantes transversales<br />
TE 10 , ce rapport est Z<br />
μ ω<br />
μ<br />
μ ε<br />
ω<br />
μ<br />
0 0 0 0<br />
0 m<br />
i = =<br />
= , que l’on peut encore écrire :<br />
kg<br />
ε0<br />
kg<br />
ε0<br />
kg<br />
k<br />
E<br />
H<br />
T<br />
T<br />
pour une onde progressive. Pour le mode<br />
λg<br />
Z i = Z 0<br />
Eq.149<br />
λm<br />
Où Z 0 est l’impédance itérative du vide.<br />
Lorsque le <strong>guide</strong> est fermé sur une impédance Z r , il se produit des ondes stationnaires. Si on introduit le coefficient de réflexion<br />
j R<br />
γ = Γ e θ pour le champ électrique, on montre que l’expression des champs est :<br />
R<br />
R<br />
j( ωt<br />
−k<br />
x)<br />
j( ωt<br />
k x)<br />
g<br />
g<br />
ET(x,<br />
t) = E0(e<br />
+ γRe<br />
) et HT<br />
(x, t) = H0(e<br />
− γRe<br />
) .<br />
Et l'impédance en un point z s'écrit:<br />
− jkz jkz<br />
e + γ R e<br />
Z(x) = Zi<br />
− jkz jkz<br />
e − γ e<br />
g +<br />
R<br />
j( ωt<br />
−k<br />
x)<br />
2<br />
g +<br />
j( ωt<br />
k x)<br />
Eq.150<br />
Le module du champ résultant peut être mis sous <strong>la</strong> forme E = E.E = A [1 + Γ + 2 Γ cos( θ + 2kz)]<br />
qui varie entre<br />
A0 [1 + ΓR<br />
] et A0 [1 − ΓR<br />
] .<br />
Le rapport d’ondes stationnaires (ROS) est alors :<br />
*<br />
2<br />
0<br />
2<br />
R<br />
R<br />
R<br />
50
ρ =<br />
E<br />
E<br />
max<br />
min<br />
1 + Γ<br />
=<br />
1 − Γ<br />
R<br />
R<br />
Eq.151<br />
La détermination de ρ à partir de <strong>la</strong> mesure de E max et E min permet de connaître Γ<br />
R<br />
. Les positions a M des maximums du champ<br />
vérifient θR = 2kg<br />
a M − 2Kπ<br />
tandis que les positions a m des minimums vérifient θ R = 2kg<br />
a m − (2K + 1)<br />
π . On peut en<br />
2π<br />
déduire, en utilisant kg<br />
= , <strong>la</strong> valeur du déphasage θ R introduit par <strong>la</strong> réflexion.<br />
λg<br />
A partir de Γ R et θ R , il est alors possible de déterminer l’impédance du <strong>guide</strong> en tout point.<br />
D.2.4.8. Interprétation de <strong>la</strong> <strong>propagation</strong> guidée<br />
Pour le mode TE 10 , l’expression du champ électrique est<br />
− jμ0ω<br />
y j( ωt<br />
−k<br />
g .x)<br />
E = H0<br />
sin( π )e ez<br />
.<br />
kc<br />
a<br />
On peut encore l’écrire :<br />
y y<br />
2Z 0 k<br />
jπ<br />
− jπ<br />
m<br />
a a j( ωt<br />
−k<br />
g .x)<br />
H0<br />
(e − e ) e e z<br />
kc<br />
π<br />
π<br />
j( ωt<br />
−(k<br />
g .x + y)) j( ωt<br />
−(k<br />
g .x − y))<br />
a<br />
a<br />
0 ( e<br />
− e<br />
) ez<br />
.<br />
E = =<br />
E<br />
π<br />
π<br />
On introduit alors les deux vecteurs k' = kgex<br />
+ ey<br />
et k" = kg<br />
ex<br />
− ey<br />
et le<br />
a<br />
a<br />
vecteur r = xe + ye<br />
(Figure 43).<br />
x<br />
y<br />
Le champ électrique se met alors sous <strong>la</strong> forme :<br />
E = E<br />
j( ωt<br />
−k'.r)<br />
j( ωt<br />
−k".r<br />
+π)<br />
0 ( e − e ) ez<br />
Il s’écrit comme <strong>la</strong> somme de deux ondes p<strong>la</strong>nes, déphasées de π et qui se propagent dans deux directions symétriques par<br />
rapport à l’axe du <strong>guide</strong>.<br />
Si on se restreint à une seule dimension pour le guidage, c'est-à-dire si on prend un <strong>guide</strong> infiniment <strong>la</strong>rge suivant OZ, alors le<br />
π<br />
vecteur d’onde de coupure est kc<br />
= ey<br />
. Le vecteur k ' a une<br />
a<br />
2 2 2<br />
norme telle que k ' = kg<br />
+ kc<br />
, qui est égale à <strong>la</strong> norme de k m ,<br />
le vecteur d’onde dans le milieu.<br />
Dans ce cas, on peut interpréter l’onde qui se propage dans le<br />
<strong>guide</strong> comme <strong>la</strong> somme de deux ondes qui se propagent dans des<br />
directions symétriques par rapport à l’axe avec une vitesse égale à<br />
<strong>la</strong> vitesse de l’onde dans le milieu illimité.<br />
On peut alors donner une interprétation graphique des résultats<br />
obtenus pour <strong>la</strong> vitesse de phase et <strong>la</strong> vitesse de groupe (Figure<br />
44). Pendant un intervalle de temps τ, une onde parcourt <strong>la</strong><br />
distance L 1 dans le milieu illimité. La distance entre les deux<br />
p<strong>la</strong>ns de phase sur l’axe du <strong>guide</strong> est L 3 .<br />
La vitesse de phase dans le <strong>guide</strong> est donc égale à <strong>la</strong> vitesse dans le milieu (c dans l’air assimilé au vide) divisée par le rapport<br />
kc<br />
ωc<br />
c c<br />
L 1 /L 3 qui vaut cosθ. Comme sinθ<br />
= = sinθ, on trouve comme précédemment: vϕ = = .<br />
km<br />
ω<br />
2<br />
2<br />
1−<br />
sin θ ωc<br />
1−<br />
2<br />
ω<br />
Pendant l’intervalle de temps τ, un paquet d’onde parcourt L 1 dans le milieu illimitée et seulement L 2 dans le <strong>guide</strong>. La vitesse<br />
de groupe est donc égale à <strong>la</strong> vitesse dans le milieu multipliée par le rapport L 2 /L 1 qui vaut cosθ. On trouve donc :<br />
2<br />
c<br />
2<br />
2 ω<br />
vg = c 1−<br />
sin θ = c 1−<br />
.<br />
ω<br />
Figure 43: Directions des vecteurs<br />
k ' et k " .<br />
Figure 44: Interprétation des vitesses de phase et de<br />
groupe<br />
51
D2.5. Conclusion<br />
La <strong>propagation</strong> d’une onde dans un <strong>guide</strong> n’est possible que si <strong>la</strong> fréquence de l’onde est supérieure à une fréquence de coupure<br />
qui dépend de <strong>la</strong> géométrie du <strong>guide</strong> et du mode de <strong>propagation</strong>. Quelle que soit sa fréquence, une onde TEM ne peut se<br />
propager dans un <strong>guide</strong>. D’une manière générale, il se propage un mé<strong>la</strong>nge d’ondes TE et d’ondes TM. Toutefois, quand on veut<br />
<strong>guide</strong>r une onde d’une fréquence donnée, on choisit les dimensions pour faire se propager le mode TE 10 .<br />
Même lorsque les pertes sont négligeables, un <strong>guide</strong> est dispersif. Cependant, si on se p<strong>la</strong>ce assez loin de <strong>la</strong> fréquence de<br />
coupure, <strong>la</strong> vitesse de phase et <strong>la</strong> vitesse de groupe tendent vers <strong>la</strong> vitesse dans le milieu illimité.<br />
Remarques :<br />
• L’étude des ondes TM dans un <strong>guide</strong> rectangu<strong>la</strong>ire est semb<strong>la</strong>ble à celles des ondes TE. On trouve par exemple les mêmes<br />
expressions pour les fréquences de coupure. Mais les couples (m, n) ne sont pas les mêmes. Il n’existe pas d’ondes TM(0,n)<br />
ou TM(m,0) car <strong>la</strong> composante longitudinale du champ électrique est un produit de sinus.<br />
• Dans les <strong>guide</strong>s à section circu<strong>la</strong>ires, <strong>la</strong> résolution des équations de Maxwell conduit à des fonctions de Bessel plutôt qu’à<br />
des fonctions trigonométriques. Il existe toujours des fréquences de coupure qui dépendent du mode de <strong>propagation</strong> et du<br />
rayon du <strong>guide</strong>. Par contre, le mode fondamental qui correspond à <strong>la</strong> longueur d'onde de coupure <strong>la</strong> plus grande n’est pas le<br />
mode TE 11 mais le mode TE 10 .<br />
• Pour les fréquences plus élevées dans le domaine optique, <strong>la</strong> <strong>propagation</strong> des ondes se fait à l’aide de fibres optiques. La<br />
fibre optique <strong>la</strong> plus simple consiste en deux cylindres concentriques de matériaux diélectriques d'indices de réfraction<br />
différents. Une infinité d’ondes lumineuse peut se propager dans une fibre optique, ce<strong>la</strong> correspond à des modes de<br />
<strong>propagation</strong>. Lorsque le diamètre de <strong>la</strong> fibre est très grand devant <strong>la</strong> longueur d’onde ( de 50 à 100 μm), on peut traiter <strong>la</strong><br />
<strong>propagation</strong> de l’onde en utilisant les rayons lumineux. Dans une fibre à saut d’indice, les trajets sont des segments de<br />
droite, les rayons subissent une réflexion totale à l’interface entre les deux diélectriques. Dans une fibre à gradient d’indice,<br />
les trajets des rayons lumineux ne sont plus rectilignes. Quelquefois le type de fibre, chacun des trajets correspond à un<br />
mode, ces fibres de fort diamètre sont des fibres multimodes. Si on utilise des fibres optiques à saut d’indice dont le<br />
diamètre est de l’ordre de grandeur de <strong>la</strong> longueur d’onde (5 à 10 μm), il faut traiter <strong>la</strong> <strong>propagation</strong> à l’aide des équations de<br />
Maxwell, comme pour le <strong>guide</strong> d’onde. Mais les conditions aux limites à l’interface des deux diélectriques sont plus<br />
complexes car le champ n’est pas nul dans le diélectrique de <strong>la</strong> gaine. Comme dans le <strong>guide</strong> d’onde, une fibre de diamètre<br />
fixée possède des fréquences de coupure et on peut y faire propager un seul mode en choisissant judicieusement <strong>la</strong><br />
fréquence de l’onde.<br />
52