21
4 Transição entre Guias de Onda Superquadráticos Este capítulo é dedicado ao estudo da matriz de espalhamento de descontinuidades em guias de onda superquadráticos através do método do casamento de modos [21]. Os campos modais em guias de onda superquadráticos foram determinados no Capítulo 3 utilizando o método variacional de Rayleigh-Ritz com funções de base trigonométricas. A aplicação deste tipo de função de base apresenta duas vantagens ao cômputo da matriz de espalhamento: 1. Permite computar um grande número de modos no guia; 2. Reduz o tempo de execução do programa, visto que a integração dupla inerente ao método do casamento de modos é resolvida analiticamente em uma das dimensões. O estudo dos guias de onda superquadráticos resulta em um eficiente algoritmo que pode ser aplicado à análise e projeto de diversos dispositivos de microondas como, por exemplo, transições entre guias de diferentes seções. Uma breve discussão do método do Casamento de Modos será apresentada a seguir. 4.1. O Método do Casamento de Modos Seja uma descontinuidade entre dois guias de onda mostrada na figura 4.1. A região à esquerda da descontinuidade, será denominada região I, seção S I , e a região à direita, região II, seção S II . Considera-se S I contida em S II . O procedimento para determinação da matriz de espalhamento da descontinuidade é apresentado em [21], e será detalhado abaixo.

4 Transição entre Guias de Onda Superquadráticos · (4.19b), e efetuando-se as diferenciações indicadas, obtém-se o seguinte conjunto de equações: j. z M m N n TM m n TM TM

Embed Size (px)

Citation preview

Page 1: 4 Transição entre Guias de Onda Superquadráticos · (4.19b), e efetuando-se as diferenciações indicadas, obtém-se o seguinte conjunto de equações: j. z M m N n TM m n TM TM

4 Transição entre Guias de Onda Superquadráticos

Este capítulo é dedicado ao estudo da matriz de espalhamento de

descontinuidades em guias de onda superquadráticos através do método do

casamento de modos [21].

Os campos modais em guias de onda superquadráticos foram determinados

no Capítulo 3 utilizando o método variacional de Rayleigh-Ritz com funções de

base trigonométricas. A aplicação deste tipo de função de base apresenta duas

vantagens ao cômputo da matriz de espalhamento:

1. Permite computar um grande número de modos no guia;

2. Reduz o tempo de execução do programa, visto que a integração dupla

inerente ao método do casamento de modos é resolvida analiticamente

em uma das dimensões.

O estudo dos guias de onda superquadráticos resulta em um eficiente

algoritmo que pode ser aplicado à análise e projeto de diversos dispositivos de

microondas como, por exemplo, transições entre guias de diferentes seções.

Uma breve discussão do método do Casamento de Modos será apresentada a

seguir.

4.1. O Método do Casamento de Modos

Seja uma descontinuidade entre dois guias de onda mostrada na figura 4.1.

A região à esquerda da descontinuidade, será denominada região I, seção SI, e a

região à direita, região II, seção SII. Considera-se SI contida em SII.

O procedimento para determinação da matriz de espalhamento da

descontinuidade é apresentado em [21], e será detalhado abaixo.

DBD
PUC-Rio - Certificação Digital Nº 0124855/CA
Page 2: 4 Transição entre Guias de Onda Superquadráticos · (4.19b), e efetuando-se as diferenciações indicadas, obtém-se o seguinte conjunto de equações: j. z M m N n TM m n TM TM

Transição entre Guias de Onda Superquadráticos 43

Figura 4.1 – Seção longitudinal de uma descontinuidade entre guias.

Uma onda eletromagnética, composta por um somatório de modos TE e

TM, com amplitudes conhecidas, incide na estrutura através da região I ou da II,

propagando-se ao longo do eixo dos z, e incidindo na descontinuidade. Deseja-se

determinar os campos espalhados (refletidos e transmitidos).

Os campos elétricos e magnéticos transversais, nas regiões I e II, em z = 0-

e z = 0+, serão expressos através de suas expansões modais:

( )∑=

+=J

jjIjIjII eBAE

1

rr (4.1a)

( )∑=

−=J

jjIjIjII hBAH

1

rr (4.1b)

( )∑=

+=I

iiIIiIIiIIII eBAE

1

rr (4.2a)

( )∑=

−=I

iiIIiIIiIIII hBAH

1

rr (4.2b)

onde os índices j e i estão associados aos pares ordenados (m1,n1) e (m2,n2), que

caracterizam os modos ( )11nmTM/TE e ( )

22nmTM/TE , nas regiões I e II,

respectivamente. jIer , jIhr

, iIIer , e iIIhr

são as componentes transversais dos campos

modais nas regiões I e II; AjI e AiII, são as amplitudes dos campos incidentes e BjI e

BiII são as amplitudes dos campos espalhados.

As condições de contorno em z = 0 (continuidade dos campos elétrico e

magnético transversais) impõem:

O

SII

SI z

Região I

Região II

DBD
PUC-Rio - Certificação Digital Nº 0124855/CA
Page 3: 4 Transição entre Guias de Onda Superquadráticos · (4.19b), e efetuando-se as diferenciações indicadas, obtém-se o seguinte conjunto de equações: j. z M m N n TM m n TM TM

Transição entre Guias de Onda Superquadráticos 44

Para pontos no interior de SI

=

=

(4.3b)

(4.3a)

III

III

HH

EErr

rr

Para pontos no interior de SII - SI

=

==

(4.4b)0

(4.4a)0

I

III

H

EEr

rr

Combinando-se as equações (4.1), (4.2), (4.3) e (4.4) resulta:

- Para pontos no interior de SI:

( ) ( )∑∑==

+=+I

iiIIiIIiII

J

jjIjIjI eBAeBA

11

rr (4.5a)

( ) ( )∑∑==

−=−I

iiIIiIIiII

J

jjIjIjI hABhBA

11

rr (4.5b)

- Para pontos no interior de SII - SI:

( ) 01

=+= ∑=

I

iiIIiIIiIIjI eBAe rr (4.6a)

( ) 01

=−= ∑=

I

iiIIiIIiIIjI hABhrr

(4.6b)

Multiplicando-se vetorialmente ambos os membros de (4.5a) e (4.6a) por

iIIhr

, i = 1, 2, ..., I, integrando-se sobre a superfície SII, e lembrando-se que:

∫ =⋅×II

IIII

S

ii sdhe 021rrr , se i1 ≠ i2 (propriedade de ortogonalidade dos modos),

obtém-se:

( ) ( )

I...,,,i

BAqBApJ

jiIIiIIiijIjIij

211

=

+=+∑= (4.7)

onde,

DBD
PUC-Rio - Certificação Digital Nº 0124855/CA
Page 4: 4 Transição entre Guias de Onda Superquadráticos · (4.19b), e efetuando-se as diferenciações indicadas, obtém-se o seguinte conjunto de equações: j. z M m N n TM m n TM TM

Transição entre Guias de Onda Superquadráticos 45

∫ ⋅×=IS

iIIjIij sdhep rrr (4.8)

∫ ⋅×=IIS

iIIiIIii sdheq rrr (4.9)

De forma análoga, multiplicando-se vetorialmente ambos os membros de

(4.5b) e (4.6b) por jIer , j = 1, 2, ..., J, e integrando-se sobre a superfície SI, tem-se:

( ) ( )

J...,,,j

ABrpBAI

ijIjIjjijiIIiII

211

=

−=−∑=

(4.10)

onde,

∫ ⋅×=IS

jIjIjj sdher rrr (4.11)

As equações (4.7) e (4.10) formam um sistema de J+I equações com J+I

incógnitas, que sob a forma matricial é expresso por:

[ ] [ ] [ ]{ } [ ] [ ] [ ]{ }IIIIII BAQBAP +=+ (4.12a)

[ ] [ ] [ ]{ } [ ] [ ] [ ]{ }IIIIIIT BARABP −=− (4.12b)

onde,

- [AI] e [BI] são matrizes colunas J x 1 contendo os coeficientes AjI e BjI;

- [AII] e [BII] são matrizes colunas I x 1 contendo os coeficientes AiII e BiII;

- [P] é uma matriz I x J, com elementos pij definidos em (4.8);

- [Q] é uma matriz diagonal I x I, com elementos qii definidos em (4.9);

- [R] é uma matriz diagonal J x J, com elementos rjj definidos em (4.11).

O sistema definido em (4.12) pode ser reescrito sob a forma:

[ ][ ] [ ][ ] [ ][ ]IIII APBQBP −=− (4.13a)

DBD
PUC-Rio - Certificação Digital Nº 0124855/CA
Page 5: 4 Transição entre Guias de Onda Superquadráticos · (4.19b), e efetuando-se as diferenciações indicadas, obtém-se o seguinte conjunto de equações: j. z M m N n TM m n TM TM

Transição entre Guias de Onda Superquadráticos 46

[ ][ ] [ ] [ ] [ ][ ]IIIT

I ARBPBR =+ (4.13b)

que resolvido resulta em:

[ ] [ ][ ]ASB = (4.14)

onde,

[ ] [ ][ ]

=

II

I

BB

B , [ ] [ ][ ]

=

II

I

AA

A e [ ] [ ] [ ][ ] [ ]

=

22211211

SSSS

S

sendo [S] a matriz de espalhamento desejada.

Após manipulações algébricas nas equações (4.13a) e (4.13b), obtém-se:

[ ] [ ] [ ] [ ] [ ]{ } [ ] [ ] [ ] [ ]{ }PQPRPQPRS TT 11111 −−− −+= (4.15a)

[ ] [ ] [ ] [ ] [ ]{ } [ ]TT PPQPRS11212

−−+= (4.15b)

[ ] [ ] [ ][ ] [ ]{ } [ ]PPRPQS T 11221−−+= (4.15c)

[ ] [ ] [ ][ ] [ ]{ } [ ] [ ][ ] [ ]{ }TT PRPQPRPQS 11122 −−− −+−= (4.15d)

Dessa forma, para o cálculo da matriz [S] basta que se determinem os

elementos das matrizes [R], [Q] e [P], e se aplique às expressões (4.15a) a (4.15d).

4.2. Obtenção das Matrizes [P], [Q] e [R]

Nesta seção serão calculados os elementos das matrizes [P], [Q] e [R], para

o caso em que os guias das regiões I e II são superquadráticos, como representado

na figura 4.2.

DBD
PUC-Rio - Certificação Digital Nº 0124855/CA
Page 6: 4 Transição entre Guias de Onda Superquadráticos · (4.19b), e efetuando-se as diferenciações indicadas, obtém-se o seguinte conjunto de equações: j. z M m N n TM m n TM TM

Transição entre Guias de Onda Superquadráticos 47

Figura 4.2 – Seção transversal de uma descontinuidade entre guias de onda

superquadráticos e sistema de coordenadas utilizado.

4.2.1. Determinação dos Campos Modais

As expressões para os campos modais nas regiões I e II, a serem utilizadas

no cálculo das matrizes [R], [Q] e [P], são obtidas através dos potencias vetores

elétrico e magnético.

Para a região I, os potenciais vetores elétrico e magnético, TEjIψ e TM

jIψ , que

são obtidos a partir das equações (3.3) e (3.7), respectivamente, são dados por:

∑∑=

=

=

maxTEjI

max

jI

M

m

zjN

n

TEnm

TEjI e

byncos

axmsenC

1 0 1

1

1

1

1 1

11 22βππ

ψ (4.16a)

∑∑=

=

+

=

maxTMjI

max

jI

M

m

zjN

n

TEnm

TMjI e

bynsen

axmcos

by

axC

0 1 1

1

1

1

01011 1

11

11 221 β

γγππ

ψ (4.16b)

Mmax e Nmax, são os valores máximos dos índices m1 e n1, respectivamente; γ1 é o

parâmetro γ relativo ao primeiro guia (seção SI). 200

2jIcTE

TEjI k−= εµωβ e

200

2jIcTM

TMjI k−= εµωβ , sendo ω a freqüência angular em rad/s; 0µ e 0ε são a

permeabilidade e permissividade do vácuo, respectivamente.

DBD
PUC-Rio - Certificação Digital Nº 0124855/CA
Page 7: 4 Transição entre Guias de Onda Superquadráticos · (4.19b), e efetuando-se as diferenciações indicadas, obtém-se o seguinte conjunto de equações: j. z M m N n TM m n TM TM

Transição entre Guias de Onda Superquadráticos 48

Os números de onda de corte dos modos TE (TM), jIcTEk (

jIcTMk ), e os

coeficientes jITEnmC

11 e jITM

nmC11

, são obtidos através do método de Rayleigh-Ritz

(Capítulo 3). Observa-se que o índice j corresponde ao par ordenado (m1,n1).

A partir das equações (4.16a) e (4.16b), de acordo com [23], obtém-se as

componentes dos campos modais transversais na região I.

- Modos TE

ye

TEjITE

x jI ∂∂

−=ψ

(4.17a)

xe

TEjITE

y jI ∂∂

(4.17b)

TEy

TEjITE

x jIjIeh

−=

0ωµβ

(4.17c)

TEx

TEjITE

y jIjIeh

=

0ωµβ

(4.17d)

Substituindo-se a expressão de TEjIψ dada em (4.16a) nas equações (4.17a) e

(4.17b), e efetuando-se as diferenciações indicadas, obtém-se o seguinte conjunto

de equações:

∑∑=

=

=

maxTEjI

max

jI

jI

M

m

zjN

n

TEnm

TEx e

bynsen

axmsen

bnCe

1 0 1

1

1

1

1

1

1 1

11 222βπππ (4.18a)

∑∑=

=

=

maxTEjI

max

jI

jI

M

m

zjN

n

TEnm

TEy e

byncos

axmcos

amCe

1 0 1

1

1

1

1

1

1 1

11 222βπππ (4.18b)

∑∑=

=

−=

maxTEjI

max

jI

jI

M

m

zjN

n

TEnm

TEjITE

x eb

yncosa

xmcosa

mCh1 0 1

1

1

1

1

1

01 1

11 222βπππ

ωµβ

(4.18c)

∑∑=

=

=

maxTEjI

max

jI

jI

M

m

zjN

n

TEnm

TEjITE

y eb

ynsena

xmsenb

nCh1 0 1

1

1

1

1

1

01 1

11 222βπππ

ωµβ

(4.18d)

DBD
PUC-Rio - Certificação Digital Nº 0124855/CA
Page 8: 4 Transição entre Guias de Onda Superquadráticos · (4.19b), e efetuando-se as diferenciações indicadas, obtém-se o seguinte conjunto de equações: j. z M m N n TM m n TM TM

Transição entre Guias de Onda Superquadráticos 49

- Modos TM

zxje

TMjITM

x jI ∂∂∂

ωε

2

0

1 (4.19a)

zyje

TMjITM

y jI ∂∂∂

ωε

2

0

1 (4.19b)

TMyTM

jI

TMx jIjI

eh

−=

βωε0 (4.19c)

TMxTM

jI

TMy jIjI

eh

=

βωε0 (4.19d)

Substituindo-se a expressão de TMjIψ dada em (4.16b) nas equações (4.19a) e

(4.19b), e efetuando-se as diferenciações indicadas, obtém-se o seguinte conjunto

de equações:

zj

M

m

N

n

TMnm

TMjITM

x

TMjI

max max

jI

jI

eb

ynsen

axm

senby

ax

am

byn

sena

xmcos

ax

aCe

βγγ

γ

πππ

ππγωεβ

= =

+

= ∑∑

1

1

1

1

01011

1

0 1 1

1

1

1

1

0101

1

0

221

2

22

11

1 1

1

11

(4.20a)

zj

M

m

N

n

TMnm

TMjITM

y

TMjI

max max

jI

jI

eb

yncos

axm

cosby

ax

bn

byn

sena

xmcos

by

bCe

βγγ

γ

πππ

ππγωεβ

= =

+

+

+

−= ∑∑

1

1

1

1

01011

1

0 1 1

1

1

1

1

0101

1

0

221

2

22

11

1 1

1

11

(4.20b)

zj

M

m

N

n

TMnm

TMx

TMjI

max max

jI

jI

eb

yncos

axm

cosby

ax

bn

byn

sena

xmcos

by

bCh

βγγ

γ

πππ

ππγ

= =

+

+

+

= ∑∑

1

1

1

1

01011

1

0 1 1

1

1

1

1

0101

1

221

2

22

11

1 1

1

11

(4.20c)

DBD
PUC-Rio - Certificação Digital Nº 0124855/CA
Page 9: 4 Transição entre Guias de Onda Superquadráticos · (4.19b), e efetuando-se as diferenciações indicadas, obtém-se o seguinte conjunto de equações: j. z M m N n TM m n TM TM

Transição entre Guias de Onda Superquadráticos 50

zj

M

m

N

n

TMnm

TMy

TMjI

max max

jI

jI

eb

ynsen

axm

senby

ax

am

byn

sena

xmcos

ax

aCh

βγγ

γ

πππ

ππγ

= =

+

= ∑∑

1

1

1

1

01011

1

0 1 1

1

1

1

1

0101

1

221

2

22

11

1

11

(4.20d)

Para a região II, segue-se o mesmo procedimento realizado para a região I,

levando-se em consideração as dimensões do segundo guia, conforme ilustrado na

figura 4.2.

O cálculo dos cálculo dos elementos das matrizes [R], [Q] e [P] será

apresentado a seguir.

4.2.2. Cálculo dos Elementos da Matriz [R]

A integração descrita na equação (4.11) é realizada na região à esquerda da

descontinuidade, guia I, e envolve os campos elétrico e magnético da região I,

seção SI, conforme mostrado na figura 4.2.

- Modos TE

De acordo com a equação (4.11):

∫ ⋅×=IS

TEjI

TEjI

TEjj sdher rrr (4.21)

que pode ser reescrita como [22]:

∫=I

jI

S

TEjI

cTETEjITE

jj dsk

r 2

0

2

ψωµ

β (4.22)

Substituindo-se em (4.22) a expressão de TEjIψ dada em (4.16a), resulta:

DBD
PUC-Rio - Certificação Digital Nº 0124855/CA
Page 10: 4 Transição entre Guias de Onda Superquadráticos · (4.19b), e efetuando-se as diferenciações indicadas, obtém-se o seguinte conjunto de equações: j. z M m N n TM m n TM TM

Transição entre Guias de Onda Superquadráticos 51

( )

dxdyb

yncosb

yncos

axmsen

axmsenCC

kr

xf

y

a

xnm

TEnm

nm

TEnm

cTETEjITE

jj

I

jIjIjI

∫∑∑

=

=

=

0 1

2

1

1

0 1

2

1

1

0

2

22

22

01

11

22

22

11

ππ

ππωµ

β

(4.23)

onde ( ) ( )[ ] 11

10101 1 γγaxbxf I −= .

- Modos TM

∫∫ =⋅×=I

jI

I S

TMjI

cTMTMjI

S

TMjI

TMjI

TMjj ds

ksdher 2

0

2

ψωε

βrrr (4.24)

Substituindo-se em (4.24) a expressão de TMjIψ dada em (4.16b), resulta:

( )

dxdyb

ynsenb

ynsenby

ax

axmcos

axmcosCC

kr

xf

y

a

xnm

TMnm

nm

TMnm

cTMTMjITM

jj

I

jIjIjI

∫∑∑

=

=

+

=

0 1

2

1

1

2

0101

0 1

2

1

1

0

2

221

22

11

01

11

22

22

11

ππ

ππωµ

β

γγ

(4.25)

4.2.3. Cálculo dos Elementos da Matriz [Q]

De forma análoga a efetuada para o cálculo da matriz [R], com mudanças

apenas nas dimensões do guia, tem-se:

- Modos TE

( )

dxdyb

yncosb

yncos

axmsen

axmsenCC

kq

xf

y

a

xnm

TEnm

nm

TEnm

cTETEiIITE

ii

II

iIIiIIiII

∫∑∑

=

=

=

0 2

2

2

1

0 2

2

2

1

0

2

22

22

02

11

22

22

11

ππ

ππωµ

β

(4.26)

DBD
PUC-Rio - Certificação Digital Nº 0124855/CA
Page 11: 4 Transição entre Guias de Onda Superquadráticos · (4.19b), e efetuando-se as diferenciações indicadas, obtém-se o seguinte conjunto de equações: j. z M m N n TM m n TM TM

Transição entre Guias de Onda Superquadráticos 52

onde ( ) ( )[ ] 21

20202 1 γγaxbxf II −= , sendo γ2 o parâmetro γ relativo ao segundo

guia (seção SII).

- Modos TM

( )

dxdyb

ynsenb

ynsenby

ax

axmcos

axmcosCC

kq

xf

y

a

xnm

TMnm

nm

TMnm

cTMTMiIITM

ii

II

iIIiIIiII

∫∑∑

=

=

+

=

0 2

2

2

1

2

0202

0 2

2

2

1

0

2

221

22

22

02

11

22

22

11

ππ

ππωµ

β

γγ

(4.27)

4.2.4. Cálculo dos Elementos da Matriz [P]

A integração descrita na equação (4.8) é realizada na região à esquerda da

descontinuidade, guia I, e envolve o campo elétrico da região I, seção SI, e o

campo magnético da região II, seção SII, conforme mostrado na figura 4.2.

- Para modos TE no primeiro guia e modos TE no segundo guia

De acordo com a equação (4.8), os elementos pij são dados por:

∫ ⋅×=IS

TEiII

TEjI

TE/TEij sdhep rrr (4.28)

Para os modos TE, os campos elétrico e magnético transversais podem ser

expressos como [22]:

TETTTH ψ

ωµβ

∇−=0

r (4.29a)

TETTzTzT aHaE ψ

βωµ

∇×=×−=rrrr

0 (4.29b)

DBD
PUC-Rio - Certificação Digital Nº 0124855/CA
Page 12: 4 Transição entre Guias de Onda Superquadráticos · (4.19b), e efetuando-se as diferenciações indicadas, obtém-se o seguinte conjunto de equações: j. z M m N n TM m n TM TM

Transição entre Guias de Onda Superquadráticos 53

Substituindo-se (4.29a) e (4.29b) em (4.28), obtém-se:

( )[ ] dsaap z

S

TETT

TETTz

TEiIITE/TE

ij

I

iIIjI

rr⋅∇×∇×−= ∫ ψψ

ωµβ

0

(4.30)

Aplicando-se a propriedade vetorial ( ) ( )BACCABCBArrrrrrrrr

⋅−⋅=⋅× na

equação (4.30), resulta:

( ) ( )[ ] dsaaap z

S

TETTz

TETT

TETT

TETTz

TEiIITE/TE

ij

I

iIIjIiIIjI

rrr⋅∇⋅∇−∇⋅∇= ∫ ψψψψ

ωµβ

0

(4.31)

A segunda parcela da integral é nula, e a equação (4.31) é reduzida a seguinte

expressão:

dspI

iIIjI

S

TETT

TETT

TEiIITE/TE

ij ∫ ∇⋅∇−= ψψωµβ

0

(4.32)

Da primeira identidade de Green [23], tem-se:

∫∫∫ ∂

∂+∇−=∇⋅∇

C

TETTE

T

S

TET

TET

S

TETT

TETT ds

ndsds jI

iIIjIiIIiIIjI

ψψψψψψ 2 (4.33)

As condições de contorno do problema impõem que 0=∂

n

TETjI

ψ . Combinando-se as

equações (4.33) e (4.32), resulta:

∫ ∇−=I

jIiII

S

TET

TET

TEiITE/TE

ij dydxp ψψωµβ 2

0

(4.34)

mas 022 =+∇ TETcTE

TET jIjIjI

k ψψ e, em conseqüência, a equação (4.34) assume a

seguinte forma:

DBD
PUC-Rio - Certificação Digital Nº 0124855/CA
Page 13: 4 Transição entre Guias de Onda Superquadráticos · (4.19b), e efetuando-se as diferenciações indicadas, obtém-se o seguinte conjunto de equações: j. z M m N n TM m n TM TM

Transição entre Guias de Onda Superquadráticos 54

∫=I

jIiII

jI

S

TET

TET

cTETEiIITE/TE

ij dydxk

p ψψωµ

β

0

2

(4.35)

Substituindo-se a expressão do potencial dada em (4.16a) em (4.35),

levando-se em consideração as dimensões do guia, obtém-se:

( )

dxdyb

yncos

byn

cos

axm

sena

xmsenCC

kp

xf

y

a

xnm

TEnm

nm

TEnm

cTETEiIITE/TE

ij

I

iIIjIjI

∫∑∑

=

=

=

0 2

2

1

1

0 2

2

1

1

0

2

22

22

01

11

22

22

11

ππ

ππωµ

β

(4.36)

- Para modos TM no primeiro guia e modos TM no segundo guia

Segundo a equação (4.8), os elementos pij são dados por:

∫ ⋅×=IS

TMiII

TMjI

TM/TMij sdhep rrr (4.37)

O procedimento para determinação dos elementos TM/TMijp é análogo ao

realizado anteriormente para os elementos TE/TEijp . Dessa forma, obtém-se a

seguinte expressão:

∫=I

iIIjI

iII

S

TMT

TMT

cTMTMjITM/TM

ij dydxk

p ψψωε

β

0

2

(4.38)

Substituindo-se a expressão do potencial dada em (4.16b) em (4.38),

levando-se em consideração as dimensões dos guias, obtém-se:

( )

dxdyb

ynsen

byn

senby

ax

by

ax

axm

cosa

xmcosCC

kp

xf

y

a

xnm

TMnm

nm

TMnm

cTMTMjITM/TM

ij

I

iIIjIiII

∫∑∑

=

=

+

+

=

0 2

2

1

1

02020101

0 2

2

1

1

0

2

2211

22

2211

01

11

22

22

11

ππ

ππωε

β

γγγγ

(4.39)

DBD
PUC-Rio - Certificação Digital Nº 0124855/CA
Page 14: 4 Transição entre Guias de Onda Superquadráticos · (4.19b), e efetuando-se as diferenciações indicadas, obtém-se o seguinte conjunto de equações: j. z M m N n TM m n TM TM

Transição entre Guias de Onda Superquadráticos 55

- Para modos TE no primeiro guia e modos TM no segundo guia

De acordo com a equação (4.8), os elementos pij são dados por:

( )∫∫ −=⋅×=I

iIIjIiIIjI

I S

TMx

TEy

TMy

TEx

S

TMiII

TEjI

TM/TEij dydxhehesdhep rrr (4.40)

Substituindo-se em (4.40) as expressões das componentes transversais dos

campos elétrico e magnético dadas no conjunto de equações (4.18) e (4.20),

levando-se em consideração as dimensões dos guias, resulta:

( )

( )

( )

( )

+

+

+

+

=

∫∫

∫∫

∫∫

∫∫

∑∑

=

=

==

==

=

=

dxdyb

ynsen

byn

cosby

axm

cosa

xmcos

bam

_dxdyb

ynsen

byn

senby

axm

sena

xmsen

am

bn

dxdyb

ynsen

byn

senax

axm

sena

xmsen

am

bn

dxdyb

ynsen

byn

senax

axm

cosa

xmsen

abn

CCp

xf

y

a

x

xf

y

a

x

xf

y

a

x

xf

y

a

x

nm

TMnm

nm

TEnm

TM/TEij

I

I

I

I

iIIjI

0 2

2

1

1

1

020 2

2

1

1

02

2

1

1

0 2

2

1

1

020 2

2

1

1

2

2

1

1

0 2

2

1

1

0 022

2

1

1

2

2

1

1

0 2

2

1

1

1

020 2

2

1

1

02

2

1

1

22222

222222

221

2222

22222

201

201

01 2

201

11

22

22

11

ππππγπ

ππππππ

ππππππ

ππππγπ

γ

γ

γ

γ

( )

( )

∫∫

∫∫

==

==

dxdyb

yncos

byn

cosby

axm

cosa

xmcos

bn

am

dxdyb

yncos

byn

cosa

xmcos

axm

cosax

bn

am

xf

y

a

x

xf

y

a

x

I

I

0 2

2

1

1

020 2

2

1

1

2

2

1

1

0 2

2

1

1

0 2

2

1

1

022

2

1

1

222222

22221

22

201

01 2

ππππππ

ππππππ

γ

γ

(4.41)

- Para modos TM no primeiro guia e modos TE no segundo guia

0=TE/TMijp (4.42)

No cômputo dos elementos das matrizes [P], [Q] e [R], as integrais em

relação a y são efetuadas analiticamente e as integrais em x numericamente.

DBD
PUC-Rio - Certificação Digital Nº 0124855/CA
Page 15: 4 Transição entre Guias de Onda Superquadráticos · (4.19b), e efetuando-se as diferenciações indicadas, obtém-se o seguinte conjunto de equações: j. z M m N n TM m n TM TM

Transição entre Guias de Onda Superquadráticos 56

4.3. Resultados Numéricos

A formulação descrita anteriormente para o cálculo da matriz de

espalhamento de descontinuidades entre guias superquadráticos de diferentes

seções transversais foi implementada computacionalmente em linguagem

FORTRAN.

Como primeira aplicação, considerou-se uma descontinuidade entre dois

guias de onda circulares de raios 5 mm e 6 mm, respectivamente. A faixa de

freqüências considerada é de 20 a 30 GHz, garantindo que apenas o modo TE11 se

propaga. Considerou-se, em ambos os guias, o número de harmônicos das funções

trigonométricas, Mmax e Nmax, iguais a 10.

A convergência do método, em função do número de harmônicos das

funções de base (Mmax = Nmax), é mostrada na figura 4.3, onde estão indicados os

valores do módulo, em dB, e da fase, em graus, do coeficiente de reflexão na porta

I, para o modo fundamental, na freqüência de 20 GHz. Observa-se que 20 modos

são suficientes para assegurar a convergência dos resultados.

Os valores de |S11| e |S21|, para o modo fundamental, em função da

freqüência, são apresentados nas figuras 4.4 e 4.5, respectivamente. Nestas figuras

também estão indicados os valores calculados pela formulação dada em [21], que

aplica o método do casamento de modos diretamente a guias circulares. Observa-

se que, para valores do |S11| acima de –30dB, as discrepâncias entre os dois

métodos são menores do que 0,8 dB. As figuras 4.6 e 4.7 apresentam o

comportamento das fases de S11 e S21, respectivamente, para o modo

fundamental.

Como segundo exemplo, a transição entre um guia de onda retangular com

seção transversal 19,05 mm x 9,525 mm para um guia circular de raio 19,05 mm,

foi considerada. A convergência dos resultados, em função do número de modos

utilizados, para a freqüência de 9 GHz, é mostrada na tabela 4.1.

A figura 4.8 apresenta os valores |S11| e |S21|, para o modo fundamental, na

faixa de freqüências de 8 a 15 GHz. Considerou-se, para o primeiro guia, o

número de harmônicos das funções trigonométricas, Mmax e Nmax, iguais a 8, e para

o segundo guia, 16. Nota-se que os valores computados utilizando o método

proposto concordam satisfatoriamente com a solução apresentada em [18],

DBD
PUC-Rio - Certificação Digital Nº 0124855/CA
Page 16: 4 Transição entre Guias de Onda Superquadráticos · (4.19b), e efetuando-se as diferenciações indicadas, obtém-se o seguinte conjunto de equações: j. z M m N n TM m n TM TM

Transição entre Guias de Onda Superquadráticos 57

apontando discrepâncias menores do que 0,4 dB. As fases de S11 e S21, para o

modo fundamental, são mostradas na figura 4.9.

Como última aplicação, foi calculada a matriz de espalhamento de uma

descontinuidade entre um guia de retangular maior para um guia circular menor.

As dimensões transversais do guia de onda retangular são 2a = 22,86 mm e

2b = 10,16 mm. Foram feitas computações para três valores de raio do guia

circular: 2,54 mm, 3,81mm e 5,08 mm. Na faixa de freqüências utilizada, o único

modo propagante é o modo fundamental TE10 no guia retangular. Todos os modos

do guia circular são evanescentes. A convergência dos resultados, em função do

número de modos utilizados, para as freqüências de 8 GHz e 14 GHz, é mostrada

na tabela 4.2.

A susceptância normalizada de uma descontinuidade pode ser determinada a

partir de S2211:

11

11

221221

SSjB

+−

= (4.70)

onde B é a susceptância normalizada da descontinuidade e S2211 é o elemento de

S22 correspondente ao primeiro modo.

A figura 4.10 mostra a susceptância normalizada da descontinuidade, B,

para três valores de raio do guia circular, em função de a/λ, onde λ é o

comprimento de onda. Considerou-se, para o primeiro guia, o número de

harmônicos das funções trigonométricas, Mmax e Nmax, iguais a 4, e para o segundo

guia, 12. Na mesma figura são apresentados os valores obtidos de acordo com a

formulação dada em [17]. Verifica-se que as discrepâncias entre os dois métodos

são pequenas.

DBD
PUC-Rio - Certificação Digital Nº 0124855/CA
Page 17: 4 Transição entre Guias de Onda Superquadráticos · (4.19b), e efetuando-se as diferenciações indicadas, obtém-se o seguinte conjunto de equações: j. z M m N n TM m n TM TM

Transição entre Guias de Onda Superquadráticos 58

Figura 4.3 – Valores do módulo de S1111 em dB e da fase em graus, em função do

número de harmônicos das funções de base (Nmax), para uma descontinuidade entre

guias de onda circulares. Os raios do primeiro e segundo guia são 5 mm e 6 mm,

respectivamente.

Figura 4.4 – Valores do módulo de S1111 em dB, em função da freqüência em GHz, para

uma descontinuidade entre guias de onda circulares. Os raios do primeiro e segundo

guia são 5 mm e 6 mm, respectivamente.

20.0 22.0 24.0 26.0 28.0 30.0

Frequência (GHz)

-50.0

-40.0

-30.0

-20.0

|S11

| (d

B)

De acordo com [21]

Método Proposto

0 4 8 12 16 20 24 28 32

No. de Harmônicos (Nmax)

-21.0

-20.8

-20.6

-20.4

|S1

1|

(dB

)

-200

-100

0

100

200

Fase

de

S11 (

Gra

us)

|S11|

Fase de S11

DBD
PUC-Rio - Certificação Digital Nº 0124855/CA
Page 18: 4 Transição entre Guias de Onda Superquadráticos · (4.19b), e efetuando-se as diferenciações indicadas, obtém-se o seguinte conjunto de equações: j. z M m N n TM m n TM TM

Transição entre Guias de Onda Superquadráticos 59

Figura 4.5 – Valores do módulo de S2111 em dB, em função da freqüência em GHz, para

uma descontinuidade entre guias de onda circulares. Os raios do primeiro e segundo

guia são 5 mm e 6 mm, respectivamente.

Figura 4.6 – Valores da fase de S1111 em graus, em função da freqüência em GHz, para

uma descontinuidade entre guias de onda circulares. Os raios do primeiro e segundo

guia são 5 mm e 6 mm, respectivamente.

20.0 22.0 24.0 26.0 28.0 30.0

Frequência (GHz)

-0.04

-0.03

-0.02

-0.01

0.00

|S21

| (d

B)

De acordo com [21]

Método Proposto

20.0 22.0 24.0 26.0 28.0 30.0

Frequência (GHz)

-100

0

100

200

Fase

de

S11 (

Gra

us) De acordo com [21]

Método Proposto

DBD
PUC-Rio - Certificação Digital Nº 0124855/CA
Page 19: 4 Transição entre Guias de Onda Superquadráticos · (4.19b), e efetuando-se as diferenciações indicadas, obtém-se o seguinte conjunto de equações: j. z M m N n TM m n TM TM

Transição entre Guias de Onda Superquadráticos 60

Figura 4.7 – Valores da fase de S2111 em graus, em função da freqüência em GHz, para

uma descontinuidade entre guias de onda circulares. Os raios do primeiro e segundo

guia são 5 mm e 6 mm, respectivamente.

Número de Modos

Guia Retangular Guia Circular S1111 S2111

TE TM TE TM

6 4 41 42 -0,127-j0,678 0,570-j0,446

9 6 63 64 -0,137-j0,678 0,566-j0,448

12 9 87 90 -0,133-j0,678 0,568-j0,448

16 12 116 121 -0,138-j0,678 0,566-j0,449

20 16 151 156 -0,136-j0,677 0,567-j0,448

Tabela 4.1 – Valores de S1111 e S2111 para uma descontinuidade entre um guia de onda

retangular com seção transversal 19,05 mm x 9,525 mm e um guia circular de raio 19,05

mm, em função do número de modos utilizados, na freqüência de 9 GHz.

20.0 22.0 24.0 26.0 28.0 30.0

Frequência (GHz)

-4

-2

0

2

Fase

de

S21 (

Gra

us)

De acordo com [21]

Método Proposto

DBD
PUC-Rio - Certificação Digital Nº 0124855/CA
Page 20: 4 Transição entre Guias de Onda Superquadráticos · (4.19b), e efetuando-se as diferenciações indicadas, obtém-se o seguinte conjunto de equações: j. z M m N n TM m n TM TM

Transição entre Guias de Onda Superquadráticos 61

Figura 4.8 – Valores dos módulos de S1111 e de S2111 em dB, em função da freqüência

em GHz, para uma descontinuidade entre um guia de onda retangular com seção

transversal 19,05 mm x 9,525 mm e um guia circular de raio 19,05 mm.

Figura 4.9 – Valores das fases de S1111 e de S2111 em graus, em função da freqüência

em GHz, para uma descontinuidade entre um guia de onda retangular com seção

transversal 19,05 mm x 9,525 mm e um guia circular de raio 19,05 mm.

8.0 9.0 10.0 11.0 12.0 13.0 14.0 15.0

Frequência (GHz)

-50.0

-40.0

-30.0

-20.0

-10.0

0.0

|S1

1|

e |S

21|

(dB

)

|S11| - De acordo com [18]

|S11| - Método Proposto

|S21| - De acordo com [18]

|S21| - Método Proposto

8.0 9.0 10.0 11.0 12.0 13.0 14.0 15.0

Frequência (GHz)

-200

-100

0

100

200

Fase

(G

raus)

S11

S21

DBD
PUC-Rio - Certificação Digital Nº 0124855/CA
Page 21: 4 Transição entre Guias de Onda Superquadráticos · (4.19b), e efetuando-se as diferenciações indicadas, obtém-se o seguinte conjunto de equações: j. z M m N n TM m n TM TM

Transição entre Guias de Onda Superquadráticos 62

Número de Modos

Guia Circular Guia Retangular S2211 em 8 GHz S2211 em 14 GHz

TE TM TE TM

6 6 42 36 -1,000+j0,027 -0,997+j0,080

9 9 64 56 -1,000+j0,027 -0,997+j0,080

12 12 90 81 -1,000+j0,027 -0,997+j0,079

16 16 121 110 -1,000+j0,027 -0,997+j0,079

20 20 156 144 -1,000+j0,027 -0,997+j0,079

Tabela 4.2 – Valores de S2211 para uma descontinuidade entre um guia de onda circular

de raio 2,54 mm e um guia retangular com seção transversal 22,86 mm x 10,16 mm, em

função do número de modos utilizados, nas freqüências de 8 GHz e 14 GHz.

Figura 4.10 – Valores da susceptância de uma descontinuidade entre um guia de onda

retangular e um guia circular, em função de a/λ, para três valores de raios, r, do guia

circular. As dimensões transversais do guia retangular são: 2a = 22,86 mm e 2b = 10,16

mm.

0.6 0.8 1.0 1.2 1.4 1.6

/ λ

0.0

1.0

2.0

3.0

Log

|B|

a

r = 5.08 mm

r = 3.81 mm

r = 2.54 mm

De acordo com [17]

Método Proposto

DBD
PUC-Rio - Certificação Digital Nº 0124855/CA