Antenassss

ii
AGRADECIMENTOS
Ao Prof. Dr. Wilson Arnaldo Artuzi Júnior por seus ensinamentos e orientações.
Ao Prof. Dr. José Ricardo Descardeci por seus ensinamentos e palavras de
incentivo.
A todos os alunos do curso de mestrado que me apoiaram e forneceram
valiosas informações. Em especial, aos amigos: Marcos Vinícius de Paula Pinto,
Michelle Foltran Miranda e Cícero Hildenberg.
A minha mãe e minha esposa pelo amor, apoio e palavras de incentivo
dedicados durante esse período tão conturbado.
Aos meus amigos Maxwell Almeida de Freitas e Marcelo Vieira Feliciano pelo
apoio, incentivo e presença nos momentos mais importantes.
Aos professores e funcionários do Departamento de Engenharia Elétrica da
Universidade Federal do Paraná.
iii
SUMÁRIO
AGRADECIMENTOS ................................................................................................... II
RESUMO ...................................................................................................................... V
ABSTRACT ................................................................................................................. VI
LISTA DE FIGURAS .................................................................................................. VII
LISTA DE TABELAS .................................................................................................. IX
SIMBOLOGIA............................................................................................................... X
CAPÍTULO 1 ................................................................................................................ 1
1 INTRODUÇÃO .................................................................................................. 1
CAPÍTULO 2 ................................................................................................................ 4
2 MÉTODO FDTD ................................................................................................ 4
2.1 INTRODUÇÃO................................................................................................... 4
2.2 O MÉTODO FDTD............................................................................................. 4
2.2.1 Grandezas eletromagnéticas e formas diferenciais .................................... 6
2.2.2 Discretização das Equações de Maxwell ...................................................... 8
2.2.3 Equação Global da Onda .............................................................................. 11
2.2.4 Modelagem de superfícies............................................................................ 12
2.2.5 Modelagem de fios ........................................................................................ 12
2.2.6 Discretização temporal ................................................................................. 13
2.3 MODELAGEM COMPUTACIONAL................................................................. 13
2.3.1 Parâmetros de Discretização ....................................................................... 14
2.3.2 Malha Hexaédrica .......................................................................................... 15
2.3.3 Malha Tetraédrica .......................................................................................... 16
CAPÍTULO 3 .............................................................................................................. 20
3 ANTENAS DE MICROFITA RETANGULAR E CIRCULAR .......................... 20
3.1 INTRODUÇÃO................................................................................................. 20
3.2 ANTENA DE MICROFITA RETANGULAR ..................................................... 21
3.3 ANTENA DE MICROFITA CIRCULAR............................................................ 22
3.4 MÉTODO DE ALIMENTAÇÃO ........................................................................ 22
3.5 IMPEDÂNCIA DE ENTRADA .......................................................................... 23
3.5.1 Equacionamento............................................................................................ 23
3.5.2 Diâmetro do condutor interno do cabo coaxial ......................................... 24
3.5.3 Perda de retorno ............................................................................................ 26
3.6 EFEITO DE BORDA ........................................................................................ 26
3.7 DIAGRAMA DE RADIAÇÃO............................................................................ 28
3.7.1 Equacionamento............................................................................................ 29
iv
CAPÍTULO 4 .............................................................................................................. 32
4 RESULTADOS PARA A ANTENA DE MICROFITA RETANGULAR ........... 32
4.1 MONTAGEM DA ESTRUTURA ...................................................................... 32
4.2 ALIMENTAÇÃO ............................................................................................... 33
4.3 ESPAÇO COMPUTACIONAL ......................................................................... 34
4.3.1 Discretização do espaço computacional em malha hexaédrica .............. 35
4.3.2 Discretização do espaço computacional em malha tetraédrica............... 35
4.4 CARACTERÍSTICAS DOS MATERIAIS UTILIZADOS ................................... 36
4.5 RESULTADOS ................................................................................................ 36
4.6 ANÁLISE DOS RESULTADOS ....................................................................... 44
4.6.1 Efeito da indutância de alimentação ........................................................... 44
4.6.2 Perda de retorno ............................................................................................ 44
4.6.3 Análise do efeito de borda............................................................................ 45
4.6.4 Diagramas de Radiação ................................................................................ 45
4.6.5 Comparação dos recursos computacionais gastos nas simulações...... 46
CAPÍTULO 5 .............................................................................................................. 47
5 RESULTADOS PARA A ANTENA DE MICROFITA CIRCULAR.................. 47
5.1 MONTAGEM DA ESTRUTURA ...................................................................... 47
5.2 ALIMENTAÇÃO ............................................................................................... 48
5.3 ESPAÇO COMPUTACIONAL ......................................................................... 48
5.4 CARACTERÍSTICAS DOS MATERIAIS UTILIZADOS ................................... 49
5.5 RESULTADOS ................................................................................................ 50
5.6 ANÁLISE DOS RESULTADOS ....................................................................... 56
5.6.1 Efeito da indutância de alimentação ........................................................... 56
5.6.2 Perda de retorno ............................................................................................ 56
5.6.3 Análise do efeito de borda............................................................................ 56
5.6.4 Diagrama de radiação ................................................................................... 57
CAPÍTULO 6 .............................................................................................................. 58
6 CONCLUSÃO.................................................................................................. 58
6.1 TRABALHOS FUTUROS................................................................................. 60
ANEXO I ..................................................................................................................... 61
ANEXO II .................................................................................................................... 63
ANEXO III ................................................................................................................... 65
ANEXO IV................................................................................................................... 71
REFERÊNCIAS BIBLIOGRÁFICAS.......................................................................... 74
v
RESUMO
O objetivo principal deste trabalho é a utilização do método FDTD (Diferenças

Finitas no Domínio do Tempo), com discretização do espaço computacional em
malha tetraédrica, a fim de se verificar a viabilidade da aplicação do método em
antenas de microfita com formato não retangular.
Inicialmente, serão analisados dois métodos de discretização: hexaédrico e
tetraédrico, empregados no estudo de uma antena de microfita quadrada. Após as
simulações, a validade dos resultados será verificada através da comparação com
os valores extraídos experimentalmente.
Como última parte do trabalho, será feita uma simulação com uma antena de
microfita em formato circular, discretizada apenas através de uma malha tetraédrica,
comparando os resultados com os valores experimentais.
Palavras chave: FDTD, antena de microfita.

vi
ABSTRACT
The main objective of this work is the use of the FDTD (Finite Difference Time
Domain) method, using tetrahedral discretization of the computational space, in order
to analyze microstrip antennas of non-rectangular shapes.
In the beginning, two discretization methods will be analyzed: hexahedral and
tetrahedral, to be applied in the study of a square microstrip antenna. After
simulations, the confiability of the results will be evaluated by means of comparisons
with experimental values.
Finally, a circular microstrip antenna will be simulated by using only tetrahedral
discretization and the results will be compared with the experimental values.
Keywords: FDTD, microstrip antenna.

vii
LISTA DE FIGURAS
FIGURA 1.1- ANTENA DE MICROFITA COM ALIMENTAÇÃO COAXIAL................. 1

FIGURA 2.1- ESTRUTURA DA CÉLULA TETRAÉDRICA. ......................................... 5
FIGURA 2.2- SUBDIVISÕES DA CÉLULA PARA ANÁLISE DE VALORES (A)
ELÉTRICOS E (B) MAGNÉTICOS. OS VOLUMES HACHURADOS INDICAM AS
SUBSEÇÕES (A) V23 E (B) V3. .................................................................................... 9
FIGURA 2.3- ANTENA DE MICROFITA QUADRADA EM MALHA HEXAÉDRICA. . 16
FIGURA 2.4- ANTENA DE MICROFITA QUADRADA EM MALHA TETRAÉDRICA.17
FIGURA 2.5- QUALIDADE DOS ELEMENTOS TETRAÉDRICOS QUE COMPÕEM A
MALHA. ...................................................................................................................... 18
FIGURA 2.6- QUALIDADE DOS ELEMENTOS TRIANGULARES QUE COMPÕEM A
MALHA. ...................................................................................................................... 19
FIGURA 3.1- ANTENA DE MICROFITA RETANGULAR COM ALIMENTAÇÃO
COAXIAL. ................................................................................................................... 21
FIGURA 3.2-ANTENA DE MICROFITA CIRCULAR COM ALIMENTAÇÃO COAXIAL.
.................................................................................................................................... 22
FIGURA 3.3-CIRCUITO EQUIVALENTE PARA O CÁLCULO DA IMPEDÂNCIA DE
ENTRADA................................................................................................................... 23
FIGURA 3.4- VARIAÇÃO TÍPICA DA RESISTÊNCIA E DA REATÂNCIA EM
FUNÇÃO DA FREQÜÊNCIA...................................................................................... 25
FIGURA 3.5- DIMENSÕES FÍSICA E EFETIVA DE UMA ANTENA DE MICROFITA
RETANGULAR. .......................................................................................................... 27
FIGURA 3.6- SISTEMA DE COORDENADAS PARA OBTENÇÃO DO DIAGRAMA
DE RADIAÇÃO. .......................................................................................................... 29
FIGURA 4.1- (A) PONTO DE ALIMENTAÇÃO (P) DA ANTENA DE MICROFITA. (B)
CONEXÃO DO CABO COAXIAL. .............................................................................. 32
FIGURA 4.2 PULSO GAUSSIANO DE EXCITAÇÃO DA FONTE.......................... 33
FIGURA 4.3- MODELO COMPUTACIONAL DA ESTRUTURA DA ANTENA DE
MICROFITA QUADRADA. ......................................................................................... 35
FIGURA 4.4- CAMPO ELÉTRICO NA FONTE EM MALHA HEXAÉDRICA. ............ 37
FIGURA 4.5- CAMPO ELÉTRICO NA FONTE EM MALHA TETRAÉDRICA. .......... 37
viii
FIGURA 4.6- RESISTÊNCIA DA ANTENA DE MICROFITA QUADRADA. .............. 38

FIGURA 4.7- REATÂNCIA DA ANTENA DE MICROFITA QUADRADA................... 39
FIGURA 4.8- PERDA DE RETORNO DA ANTENA QUADRADA SEM O AJUSTE DA
INDUTÂNCIA DE ALIMENTAÇÃO............................................................................. 40
FIGURA 4.9- PERDA DE RETORNO DA ANTENA QUADRADA COM O AJUSTE
DA INDUTÂNCIA DE ALIMENTAÇÃO....................................................................... 41
FIGURA 4.10- PLANO = 0 PARA CONSTRUÇÃO DO DIAGRAMA DE
RADIAÇÃO DA ANTENA QUADRADA...................................................................... 42
FIGURA 4.11- DIAGRAMA DE RADIAÇÃO DA ANTENA QUADRADA PARA =0 .
.................................................................................................................................... 42
RADIAÇÃO DA ANTENA QUADRADA...................................................................... 43
FIGURA 4.13- DIAGRAMA DE RADIAÇÃO DA ANTENA QUADRADA PARA =
45 ............................................................................................................................... 43
FIGURA 5.1- (A) PONTO DE ALIMENTAÇÃO (P) DA ANTENA DE MICROFITA
CIRCULAR. (B) CONEXÃO DO CABO COAXIAL..................................................... 47
FIGURA 5.2- MODELO COMPUTACIONAL DA ESTRUTURA DA ANTENA DE
MICROFITA CIRCULAR ............................................................................................ 48
FIGURA 5.3- CAMPO ELÉTRICO NA FONTE DA ANTENA CIRCULAR................. 50
FIGURA 5.4- RESISTÊNCIA E REATÂNCIA NA ANTENA DE MICROFITA
CIRCULAR. ................................................................................................................ 51
FIGURA 5.5- PERDA DE RETORNO DA ANTENA CIRCULAR SEM O AJUSTE DA
FIGURA 5.6- PERDA DE RETORNO DA ANTENA CIRCULAR COM O AJUSTE DA
RADIAÇÃO DA ANTENA CIRCULAR........................................................................ 54
FIGURA 5.8- DIAGRAMA DE RADIAÇÃO DA ANTENA CIRCULAR PARA = 30 .55
FIGURA 6.1- ANTENA DE MICROFITA SOBRE SUPERFÍCIE CÔNICA. ............... 60
ix
LISTA DE TABELAS
TABELA 1- MATERIAIS UTILIZADOS NA ANTENA QUADRADA ........................... 36

TABELA 2-EFEITO DA INDUTÂNCIA DO CONDUTOR DE ALIMENTAÇÃO NA
ANTENA QUADRADA................................................................................................ 44
TABELA 3 -COMPARAÇÃO DO EFEITO DE BORDA .............................................. 45
TABELA 4- COMPARAÇÃO DOS MÉTODOS DE DISCRETIZAÇÃO PARA UMA
ANTENA DE MICROFITA QUADRADA..................................................................... 46
TABELA 5- MATERIAIS UTILIZADOS NA ANTENA CIRCULAR ............................. 49
TABELA 6- EFEITO DA INDUTÂNCIA DO CONDUTOR DE ALIMENTAÇÃO NA
ANTENA CIRCULAR.................................................................................................. 56
x
SIMBOLOGIA
E Campo elétrico
H Campo magnético
D Densidade de fluxo elétrico
B Densidade de fluxo magnético
J Densidade de corrente elétrica
Permeabilidade magnética
Permissividade elétrica
Condutividade elétrica
Densidade de carga
Relutividade magnética
[dp,ij] Operador primário de derivada exterior discreta
[dkl,q] Operador secundário de derivada exterior discreta
JS Fonte de corrente de excitação
Vij Fração de volume corresponde à aresta ij
Vp Fração de volume corresponde à face p
Wij Função de base
Wp Função de base
Operador Hodge
dLij Elemento diferencial de superfície secundária
dSp Elemento diferencial de comprimento secundário
Wql Função de base
Wq Função de base
[Eij] Matriz de coeficientes numéricos de campo elétrico
[Bp] Matriz de coeficientes numéricos de campo magnético
[Jij] Matriz de coeficientes numéricos de correntes
[ kl,ij] Matriz de constantes numéricas
Ne Número de arestas da célula poliédrica
Nf Número de faces da célula poliédrica
LMED Tamanho médio da aresta
RMIN Raio do círculo inscrito na menor célula
xi
c Velocidade da luz no vácuo

r Permissividade elétrica relativa
fmax Freqüência máxima de excitação
t Espessura da antena de microfita
0 Comprimento de onda no espaço livre
h Espessura do dielétrico
L Comprimento da antena de microfita retangular
W Largura da antena de microfita retangular
r Raio da antena de microfita circular
reff Raio da antena de microfita circular considerando o efeito de borda
Za Impedância da antena
Ra Resistência da antena
Xa Reatância da antena
Freqüência angular
I Transformada de Fourier da corrente de excitação da antena
V Transformada de Fourier do campo elétrico na fonte
R Resistência da fonte
Xf Reatância de alimentação
Impedância intrínseca do meio
k Constante de propagação da onda
0 Permeabilidade elétrica no vácuo
0 Permissividade magnética no vácuo
Comprimento médio da aresta da fonte
d diâmetro do condutor coaxial
dest diâmetro estimado do condutor coaxial
deff diâmetro efetivo do condutor coaxial
XMIN Valor mínimo da curva de reatância
XMAX Valor máximo da curva de reatância
Lf Indutância do condutor coaxial
Zf Impedância do equipamento de medição
reff Constante dielétrica efetiva
L Variação no comprimento da antena devido ao efeito de borda
xii
f0 Freqüência de ressonância
A Vetor potencial magnético
Js Vetor densidade de corrente superficial
J Vetor densidade de corrente na direção
J Vetor densidade de corrente na direção
xi coordenada x do centro da célula superficial

yi coordenada y do centro da célula superficial
Si Área da célula superficial da antena
N Número de células em que foi dividida a superfície da antena
AF Fator de correção para plano terra infinito
i Corrente de excitação da antena
I0 Valor de pico da corrente de excitação da antena
T Duração da fonte
t Passo da simulação
Rs Resistência superficial
Rl Resistência linear
r Variação do raio da antena circular devido ao efeito de borda
ABC Condição de contorno de absorção (Absorbing boundary condition)
PEC Condutor elétrico perfeito (Perfect electric conductor)
PMC Condutor magnético perfeito (Perfect magnetic conductor)
1
CAPÍTULO 1
Introdução Geral
1 INTRODUÇÃO
Uma antena de microfita consiste de uma fina camada metálica, que atua como
elemento radiador, separada de seu plano terra por uma camada composta de um
substrato dielétrico. Sua alimentação pode ser feita de várias maneiras, sendo a
alimentação via cabo coaxial aplicada às antenas utilizadas neste trabalho.
Antena metálica
Condutor central do
cabo coaxial
Substrato dielétrico
Plano terra
Blindagem do cabo
coaxial
Figura 1.1- Antena de microfita com alimentação coaxial.

2
O conceito de Antena de Microfita surgiu há aproximadamente cinqüenta anos

atrás com o trabalho de Deschamps [3] nos Estados Unidos e de Gutton e Baissinot
[2] na França. Nenhum outro trabalho de mesma importância foi apresentado até o
início dos anos 70, quando Byron [7], em 1970, descreveu um radiador a partir de
uma fita condutora separada de um plano terra através de um substrato dielétrico.
Dois anos mais tarde, em 1972, dados sobre elementos de microfita de
formas circulares e retangulares foram apresentados por Howell [4], mas foi Munson
[8], que em 1973, patenteou pela primeira vez um elemento de microfita, visando
aplicações em foguetes. Weinschel [9], em 1973, desenvolveu várias geometrias de
elementos de microfita para uso em arranjos operando na banda-S de freqüência,
montados sobre estruturas cilíndricas.
Sanford [10], em 1974, mostrou em seu trabalho que o elemento de microfita
poderia ser utilizado em um arranjo conformal de elementos, desenvolvido para
operação na banda-L de freqüência, utilizado na comunicação de uma aeronave KC-
135 com um satélite ATS-6.
A introdução de elementos conformais para utilização em veículos
aeroespaciais, e a óbvia possibilidade de implementação de arranjos de elementos
de microfita, de baixo peso, foram responsáveis pelo grande interesse na expansão
desta tecnologia.
Sistemas de antenas de microfita são empregados em grande número nas
comunicações sem fio, pois eles se adaptam bem a qualquer superfície, são de
pequena dimensão e podem fornecer igual, ou melhor, performance elétrica do que
as arquiteturas das antenas clássicas [2]. Além disso, antenas integradas
diretamente nos equipamentos de comunicação, tendem a ser mais confiáveis e
possuem um custo de instalação muito menor [2]. Desta forma, é extremamente
importante o desenvolvimento de estudos a respeito desta tecnologia, tendo em
vista sua grande importância para as comunicações sem fio, aplicações militares e
na área de telefonia.
O método FDTD vem sendo aplicado em larga escala para o estudo das
antenas de microfita [1], tendo em vista que nos últimos anos houve uma enorme
evolução dos computadores, máquinas estas que tornam possível a resolução de
3
milhões de cálculos em frações de segundos, tornando possível a implementação do

método FDTD.
Neste trabalho vamos analisar os principais parâmetros (impedância de
entrada, diagrama de radiação e freqüência de ressonância) de uma antena de
microfita circular e outra retangular, utilizando-se do método FDTD. Conjuntamente
com essas análises, estudaremos também a relação dos resultados com o método
de discretização utilizado para cada estrutura, fazendo comparativos em termos dos
parâmetros supracitados e dos esforços computacionais (tempo gasto nas
simulações e memória computacional utilizada) aplicados na obtenção dos
resultados. Ainda, será apresentado um novo método de discretização, utilizando
malha tetraédrica, visando a análise de antenas conformais de formato não
retangular.
No capítulo 2 será apresentada uma breve introdução ao método FDTD,
seguido do desenvolvimento do novo método de discretização para antenas de
microfita de formato não retangular. No capítulo 3 serão descritos os parâmetros
analisados das antenas de microfita com seus respectivos equacionamentos. No
capítulo 4 e 5 serão apresentados e analisados os resultados obtidos para as
antenas de microfita retangular e circular, respectivamente. A conclusão deste
trabalho será apresentada no capítulo 6.
4
CAPÍTULO 2
Método FDTD e Formas Diferenciais
2 MÉTODO FDTD
2.1 INTRODUÇÃO
A origem do método FDTD se deu em 1966. Neste ano, Kane S. Yee
desenvolveu um trabalho que chamou de Numerical Solution of Initial Boundary
Value Problems Involving Maxwell s Equations in Isotropic Media [10]. O trabalho de
Yee marcou o início do estudo do método FDTD. O termo FDTD foi utilizado pela
primeira vez por Allen Taflove em 1980 [12].
Apesar da sua simplicidade, o método FDTD teve uma lenta evolução por
necessitar de computadores com um elevado desempenho, pois requer muita
mémoria e velocidade de processamento para sua execução. Contudo, a evolução
da microinformática foi fortíssima nas últimas décadas, fazendo com que o
desempenho das estações de trabalho de alguns anos atrás tenha sido superado
pelos atuais microcomputadores domésticos.
2.2 O MÉTODO FDTD

Nesta dissertação será utilizada uma formulação diferente daquela
convencionalmente adotada por Taflove [12], a fim de introduzir a discretização em
malha tetraédrica. Esta formulação foi inicialmente apresentada em [11] e será
explanada em detalhe ao longo deste capítulo.
O tratamento matemático para a aplicação do método será focado,
inicialmente, nos elementos de malha e a solução geral do problema será alcançada
posteriormente por meio de uma superposição de elementos. Chamaremos de célula
a região de meio homogêneo compreendida entre determinado número de faces,
5
contendo uma malha elementar primária (arestas que compõem as faces) e outra
secundária (arestas interiores que partem de um ponto no interior da célula até um
ponto localizado em cada uma das suas faces).
Por simplificação, a descrição do método se baseará em uma célula
tetraédrica, conforme mostrado na fig. 2.1, podendo ser ampliada a qualquer célula
poliédrica convexa.
3
2
1
P
3
Figura 2.1- Estrutura da célula tetraédrica.
Na figura 2.1, os pontos 1,2,3 e 4 compõem a malha primária (linhas sólidas) e

os pontos 0, 1 , 2 , 3 e 4 referem-se à malha secundária (linhas pontilhadas). As
superfícies hachuradas indicam a face primária 3 (azul) e face secundária 23
(laranja). Os pontos 1,2,3 e 4 são fornecidos por meio de softwares geradores de
malha e formam a malha primária. A reta formada por dois destes pontos é chamada
de aresta primária, será denominada pelos dois índices dos pontos que a forma.
Três destas arestas formam a face primária, que será denominada pelo índice do
ponto que não faz parte de nenhuma destas arestas. A malha secundária consiste
no ponto 0 no interior da célula (baricentro do poliedro) e por um ponto no baricentro
de cada face 1 , 2 , 3 e 4 . A reta formada pelo ponto 0 e pelo baricentro da face
6
primária será denominada aresta secundária, cujo índice será dado pelo índice do
ponto da face na qual está ligada. Duas destas arestas secundárias definem uma
face secundária que terá como índice os pontos que não fazem parte desta face. O
ponto P é o baricentro da aresta.
2.2.1 Grandezas eletromagnéticas e formas diferenciais

A partir das leis de Maxwell na forma integral, podemos determinar o grau da
forma diferencial que representará os vários valores de campo. Na notação vetorial,
temos as seguintes grandezas eletromagnéticas:
B ds E dl , (2.1)
tS L
J ds D ds H dl , (2.2)
S tS L
D ds dv , (2.3)
S V
B ds 0, (2.4)
S
onde E é o vetor campo elétrico, B é o vetor fluxo magnético, H é o vetor campo

magnético, J é o vetor densidade de fluxo elétrico e D é o vetor densidade de fluxo
elétrico. S é a superfície delimitada pelo caminho L, V é o volume delimitado pela
superfície S e é a densidade volumétrica de carga. A intensidade de campo
elétrico é integrada sobre um caminho, então ele se torna uma forma diferencial de
grau 1. A intensidade de campo magnético, como também é integrada sobre um
caminho, se torna, também, uma forma diferencial de grau 1. As densidades de fluxo
elétrico e magnético são integradas em uma superfície, tornando-se formas
diferenciais de grau 2. As fontes são densidades de corrente elétrica, que são
formas diferenciais de grau 2. A densidade de carga volumétrica é integrada em um
volume, sendo uma forma diferencial de grau 3. Existe ainda o potencial elétrico ,
que é uma forma diferencial de grau zero.
7
Sendo E E1 x E2 y E3 z uma grandeza vetorial dada em v/m, e
dl dx x dy y dz z , um elemento vetorial diferencial de comprimento, dado em

metros, tal que :
E dl E1dx E2 dy E3dz E , (2.5)
onde E é uma forma diferencial de grau 1, dada em volts, pronta para ser integrada
ao longo de uma linha uma vez que inclui a diferencial de comprimento.
De forma semelhante, B B1 x B2 y B3 z é uma grandeza vetorial dada em
Wb/m2, e ds dzdy x dxdz y dxdy z , um elemento vetorial diferencial de superfície,

dado em m2, tal que:
B ds B1dydz B2 dxdz B3dydx B, (2.6)
onde B é uma forma diferencial de grau 2, dada em Wb, pronta para ser integrada
em uma superfície uma vez que inclui a diferencial de área.
Analogamente para as demais grandezas eletromagnéticas, usando as
formas diferenciais anteriormente definidas, as leis de Maxwell podem ser escritas
como[16]:
B E , (2.7)
tS L
J D H, (2.8)
S tS L
D , (2.9)
S V
B 0 . (2.10)
S
As formas diferenciais são contrações do produto escalar de uma grandeza

pelo seu respectivo elemento diferencial de integração.
Entretanto, para o desenvolvimento a seguir, será necessária a utilização das
equações de Maxwell em sua forma pontual, como se segue:
B dE , (2.11)
t
D J dH , (2.12)
t
8
onde
d dx dy dz . (2.13)
x y z
2.2.2 Discretização das Equações de Maxwell

As equações (2.11) e (2.12) podem ser discretizadas em relação a uma célula
e representadas por equações matriciais pela integração de formas-1 nas arestas e
formas-2 nas faces, obtendo-se:
Bp d p ,ij Eij , (2.14)

t
Dkl J kl d kl , q Hq , (2.15)
t
onde
Bp B e E ij E , (2.16)
p ij
são escalares que representam o fluxo magnético na face primária p e a diferença

de potencial elétrico ao longo da aresta primária ij, e
Dkl D , Hq H e J kl J, (2.17)
kl q kl
são escalares que representam o fluxo elétrico na face secundária kl, a diferença de
potencial magnético ao longo da aresta secundária q e a densidade de corrente
elétrica na face kl, respectivamente. d p ,ij e d kl ,q são matrizes cujos elementos
são +1, -1 e 0, os quais realizam a somatória das quantidades associadas às arestas

T
de uma face, respeitando suas orientações, sendo que d kl ,q d p ,ij [11], onde o
índice T indica a matriz transposta.

Relações métricas e constitutivas estão contidas nas seguintes equações
D E , H B e J E JS , (2.18)
sendo , e , respectivamente, permissividade elétrica, relutividade magnética e

condutividade elétrica que caracterizam o meio homogêneo no interior da célula e JS
é a fonte de corrente (excitação). O símbolo representa o operador Hodge, que
9
converte uma forma de grau a em uma forma de grau 3 a , preservando sua

ortogonalidade.
Para conseguir a forma discreta das relações constitutivas em (2.18), é
necessária a realização de algumas aproximações. No método FETD (Elementos
finitos no domínio do tempo), B e E dentro da célula são normalmente expandidos
em funções de base linear, porém, neste trabalho usaremos funções de base
constante para expandir E e H no interior da célula[11].
2 1 2 1
3 3
4 4
(a ) (b )
Figura 2.2- Subdivisões da célula para análise de valores (a) elétricos e (b) magnéticos.
Os volumes hachurados indicam as subseções (a) V23 e (b) V3.
Para análise de valores elétricos subdividiremos a célula em 6 subseções

volumétricas, cada uma contendo uma aresta primária e uma face secundária,
chamadas de Vij, conforme mostrado na figura 2.2(a). Para análise dos valores
magnéticos a célula será subdividida em 4 subseções volumétricas, cada uma
contendo uma face primária e uma aresta secundária, chamadas de Vp, conforme
figura 2.2(b).
Assim as grandezas E e B podem ser expandidas como
E EijWij e B B pW p , (2.19)
ij p
dLij / Lij dS p / S p
onde Wij e Wp ,
0, fora de Vij 0, fora de Vp
sendo dLij é um elemento diferencial de comprimento de aresta primária e dSp é um

elemento diferencial de superfície primária.
10
Substituindo (2.19) em (2.18) temos que

dS p / S p , em V p
H Bp Wp , Wp , (2.20)
p 0, fora de V p
dLij / Lij em Vij

D Eij Wij , Wij , (2.21)
ij 0, fora de Vij
dLij / Lij , em Vij

J Eij Wij J ij Wij , Wij , (2.22)
ij ij 0, demais volumes
onde
dLij é um elemento diferencial de superfície secundária ortogonal a dLij.
dSp é um elemento diferencial de comprimento secundário ortogonal a dSp.
Então, substituindo (2.19) a (2.22) nas Leis de Maxwell (2.11) e (2.12) temos,
respectivamente:
Eij dWij B pW p , (2.23)

ij t p
Bp d Wp Eij Wij J ij Wij Eij Wij . (2.24)

p t ij ij ij
Valendo-se do método de Galerkin [5], multiplicaremos ambos os membros da

Lei de Faraday na equação (2.23) por ( Wq), integrando-se sobre o volume do
elemento. Faremos o mesmo com a Lei de Ampère-Maxwell na equação (2.24),
utilizando Wql, obtendo-se
d q ,ij Eij q, p Bp , (2.25)

t
d kl , p B p kl ,ij Eij kl ,ij Eij kl ,ij J ij , (2.26)

t
onde
Eij, Bp são os coeficientes numéricos a serem determinados.
Jij são coeficientes numéricos conhecidos.
e
q, p WqW p , (2.27)
11
kl ,ij Wkl Wij , (2.28)
kl Wkl Wij , (2.29)
kl ,ij Wkl Wij , (2.30)
d q ,ij Wq dWij , (2.31)
d kl , p Wkl d W p (2.32)
são constantes conhecidas em função dos parâmetros geométricos e constitutivos

do elemento.
A expressão (2.25) é a equação discretizada da Lei de Faraday, e associa o
somatório das diferenças de potencial ao longo das arestas de uma face com o fluxo
magnético através da referida face.
A expressão (2.26) é a equação discretizada da Lei de Ampère-Maxwell, e
associa a somatória dos fluxos magnéticos através das faces comuns a uma aresta
com a diferença de potencial ao longo da referida aresta.
2.2.3 Equação Global da Onda
Para se obter a equação global da onda multiplicaremos a equação (2.25) por

1
q, p e derivaremos a equação (2.26) em relação ao tempo. Após estas
operações, substituiremos a equação (2.25) modificada após a multiplicação por

1
q, p em (2.26), também modificada após a derivação em relação ao tempo,
obtendo-se
kl , ij Eij kl ,ij Eij kl ,ij E ij kl ,ij J ij , (2.33)

t2 t t
onde
1
kl ,ij d kl , p q, p d q ,ij .
Podendo-se reescrever a equação (2.33)

2
E E E J , (2.34)
t2 t t
12
que é a Equação Global da Onda, na qual desejamos obter os valores de E nas

arestas do poliedros que discretizam o espaço computacional como resposta ao
fluxo de corrente de excitação J.
Considerando a subdivisão baricêntrica da célula, dividindo o seu volume V
em partes iguais, obtemos as seguintes constantes volumétricas:
V V V V
kl ,ij kl ,ij , kl ,ij kl ,ij , kl ,ij kl ,ij e q, p q, p ,
N e Lij 2 N e Lij 2 Ne Lij 2 2
Sp N f
(2.35)
onde kl,ij e q,p são as funções delta de Kronecker, Ne é o número de arestas, Nf é o
número de faces do poliedro no qual a malha foi discretizada e Lij a aresta do
poliedro.
2.2.4 Modelagem de superfícies

Analogamente às constantes volumétricas, considerando a subdivisão
baricêntrica da célula triangular ou quadrada de área S, temos que:
S
kl ,ij kl ,ij , (2.36)
Rs N e Lij 2
onde Rs é a resistência superficial do material, utilizando

PEC RS = 10-6 (RS 0)
10
PMC RS = 10 (RS )
ABC Rs , onde e são características eletromagnéticas do meio
adjacente à superfície
Ar RS = 377
Cobre RS = 10-3
Teflon RS = 236
2.2.5 Modelagem de fios

Para a modelagem de fios utilizaremos a seguinte equação:
Lij
ij ,ij , (2.37)
Rl
13
onde Rl é a resistência linear do material, utilizando os mesmos valores adotados

para RS nos materiais citados no subitem 2.2.4.
2.2.6 Discretização temporal

Para a discretização temporal utilizaremos o método da diferença finita
centrada, que fará a interação dos campos elétrico e magnético no tempo, partindo
da seguinte expressão:
t t
f (t ) f (t )
f (t ) 2 2 , (2.38)
t t
podemos, então, fazer a seguinte simplificação:
t=n t E(t) = En , (2.39)
onde t é o passo no tempo e n é um índice que representa o número de passos.
Com isso, podemos fazer
E(t+ t) = En+1 , (2.40)
E(t- t) = En-1 . (2.41)
Substituindo (2.39), (2.40) e (2.41) em (2.34) temos que:
1 1 2 1 1
2
[ ] [ ] [E n 1] 2
[ ] [ ] [E n ] 2
[ ] [ ] [E n 1] [ ][ J n ] ,
t 2 t t t 2 t t
(2.42)
que é a equação da onda discretizada no tempo e no espaço.
2.3 MODELAGEM COMPUTACIONAL

A modelagem computacional das estruturas foi realizada através da construção
das antenas de microfita em um software, onde foram respeitadas as dimensões
reais do problema. Para a implementação do método FDTD é necessário efetuar-se
a discretização da estrutura considerada em pequenos poliedros, em cujas arestas
serão obtidos valores de campo elétrico. As dimensões desses poliedros
dependerão dos valores dos parâmetros de discretização que devem obedecer aos
limites de dispersão e estabilidade numérica do método FDTD [17].
Após os cálculos dos parâmetros de discretização, deve-se evitar a utilização
de valores de arestas muito pequenos, pois isso, dependendo da capacidade de
processamento do computador que se está utilizando, pode inviabilizar a solução,
14
pois demandará grande quantidade de memória RAM e muito espaço no disco

rígido.
Obtém-se, também, a partir desta discretização, as coordenadas dos vértices
que compõem o poliedro, tendo-se, então, parâmetros suficientes para o cálculo de
comprimento das arestas e áreas das faces, elementos estes, essenciais para a
implementação do método FDTD.
2.3.1 Parâmetros de Discretização

Para a divisão do espaço computacional em poliedros, devem ser respeitados
alguns critérios de dispersão e estabilidade numérica do método [17]. Para
estabelecer esses critérios deve-se, em primeiro lugar, determinar o limite máximo
de freqüência que será excitado na simulação.
O princípio básico a ser seguido para o dimensionamento da célula é que a
célula na qual a estrutura foi discretizada deve ser bem menor que o menor
comprimento de onda excitado na simulação. Quanto menor for a célula maior será a
precisão da resposta. O que se deve ter em mente é que devemos respeitar o limite
de amostragem de Nyquist, onde, = 2 x, com x sendo o comprimento médio das
arestas. O fato é que, como não podemos determinar precisamente o comprimento
de onda e a discretização não é uniforme, adotam-se dimensões de células bem
menores do que o comprimento de onda, fazendo x /10 [17].
Por isso, foi verificado a partir de resultados numéricos obtidos [17] que, para
se respeitar o limite de dispersão numérica, deve ser observada a seguinte equação:
c
LMED , (2.43)
10 f max r
onde, através da equação (2.43), definindo o limite máximo de freqüência de

excitação fmax, podemos obter o comprimento médio das arestas, LMED. Onde c é a
velocidade da luz no vácuo e r é a permissividade elétrica relativa do material
dielétrico de maior permissividade aplicado na estrutura a ser simulada.
15
Para definir-se o passo no tempo devemos atender a condição de Courant [17],

que, para uma célula cúbica, pode ser obtido através da seguinte expressão:
Lcub
t (2.44)
c 3
onde Lcub é o comprimento da aresta.
Para entendermos melhor a base para a condição de Courant,
consideremos uma onda plana se propagando em uma malha discretizada em
FDTD. Em um passo no tempo t, qualquer ponto dessa onda não deve passar por
mais de uma célula, pois durante um passo no tempo a onda somente poderá se
propagar do interior de uma célula até o limite interno de suas arestas.
Com as dimensões médias da célula definidas a partir da equação (2.43),
obtém-se o passo de tempo máximo para a simulação utilizando a seguinte
expressão:
2 Rmin
t , (2.45)
3c
onde, Rmin é o raio do círculo inscrito na menor célula oriunda da discretização do
espaço computacional [11]. O valor de Rmin dependerá da malha que será gerada
por software. A expressão (2.45) é obtida empiricamente com base na condição de
Courant [17], definida em (2.44).
2.3.2 Malha Hexaédrica

Para a discretização em malha hexaédrica, teremos que segmentar as
arestas das estruturas nas direções x, y e z, de forma que todo o espaço
computacional tridimensional seja formado por pequenos hexaedros e as superfícies
discretizadas em quadriláteros, conforme se vê na figura 2.3.
Para a escolha do número de partes em que serão divididos os lados da
estrutura, deve-se respeitar o limite de dispersão numérica do método, observando a
equação (2.43).
Pode-se observar na figura 2.3 que a estrutura foi discretizada em hexaedros
com arestas menores perto das bordas da antena de microfita, resultando em uma
concentração maior de células em alguns pontos da estrutura. Isso foi feito para se
evitar variações bruscas nos valores de campo quando se muda de um meio
16
material para outro, tendo em vista que nessa situação ocorrerá uma mudança na
velocidade da onda [17].
Figura 2.3- Antena de microfita quadrada em malha hexaédrica.
2.3.3 Malha Tetraédrica

Na discretização em malha tetraédrica, dividiremos todo o espaço
computacional em poliedros, utilizando como valor de aresta, o resultado da análise
dos parâmetros de discretização. Para esse tipo de discretização a única entrada
será o valor da aresta escolhida a partir da equação (2.43), não sendo necessária a
seleção dos lados que serão particionados, como no caso da malha hexaédrica. Na
discretização em malha tetraédrica é difícil estimar o valor das arestas em
determinados pontos da estrutura, uma vez que esta estimativa dependerá
principalmente da geometria do problema. Pode-se observar na figura 2.4 uma
estrutura discretizada em malha tetraédrica. Neste tipo de discretização as
superfícies são compostas por triângulos.
17
Figura 2.4- Antena de microfita quadrada em malha tetraédrica.
Pela mesma razão que no caso da malha hexaédrica, algumas regiões foram
discretizadas com arestas menores causando uma concentração de células em
alguns pontos da estrutura. No caso da malha tetraédrica, esse artifício evita,
também, deformações nas arestas geradas pelo método de discretização,
melhorando a qualidade da malha gerada, conforme pode ser observado nas figuras
2.5 e 2.6.
18
A figura 2.5 é obtida através do software gerador da malha, e indica a

qualidade da malha discretizada em função do ângulo sólido formado pelas faces do
tetraedro. O eixo das ordenadas apresenta o número de elementos tetraédricos em
que foi discretizado o espaço computacional e o eixo das abscissas apresenta os
ângulos em graus.
Figura 2.5- Qualidade dos elementos tetraédricos que compõem a malha.

19
A figura 2.6 é obtida através do software gerador da malha, e indica a

qualidade da malha discretizada em função do ângulo plano formado pelas arestas
das faces. O eixo das ordenadas apresenta o número de elementos triangulares em
que foi discretizado o espaço computacional e o eixo das abscissas apresenta os
ângulos em graus.
Figura 2.6- Qualidade dos elementos triangulares que compõem a malha.

20
CAPÍTULO 3
Parâmetros Analisados nas Antenas de Microfita
3 ANTENAS DE MICROFITA RETANGULAR E CIRCULAR
3.1 INTRODUÇÃO
Apesar das antenas de microfita terem surgido a mais de meio século [3],
esses dispositivos apenas começaram a ganhar alguma atenção nos anos 70.
Antenas de microfita, como a mostrada na figura 3.1, consistem de uma superfície
metálica bem fina (t 0, onde 0 é o comprimento de onda no espaço livre na
freqüência de operação e t a espessura da antena metálica), assentada sobre um
substrato dielétrico de pequena espessura (h 0, geralmente 0,003 0 h
0,05 0) sobre um plano terra. A antena de microfita é desenvolvida para que sua
irradiação máxima seja normal a sua superfície. Isto pode ser alcançado através de
uma escolha apropriada do modo de excitação (configuração do campo) embaixo da
antena.
21
substrato
dielétrico
Antena
metálica W
h
z
y
Alimentação plano terra
x coaxial
Figura 3.1- Antena de microfita retangular com alimentação coaxial.
Freqüentemente as antenas de microfita são denominadas, também, como

antenas patch. O elemento irradiante da antena de microfita pode ser quadrado,
retangular, em forma de fita (dipolo), circular, elíptico e triangular, entre outros.
3.2 ANTENA DE MICROFITA RETANGULAR

Para uma antena de microfita retangular, como a vista na figura 3.1, a
dimensão L situa-se , geralmente, na faixa entre 0 /3 L 0 / 2 . O substrato
dielétrico que separa a antena do plano terra possui constante dielétrica, r, com
valor situado na faixa entre 2, 2 r 12 . Os substratos utilizados nos dispositivos
em que se deseja uma boa performance como antena, são mais espessos e
possuem constantes dielétricas situadas próximas ao limite inferior da faixa, pois,
dessa forma, eles fornecem uma melhor eficiência e maior largura de banda [5]. Um
material comumente utilizado como substrato dielétrico em antenas de microfita é o
PTFE (politetrafluoretileno), material semelhante ao teflon.
Substratos menos espessos com constantes dielétricas maiores dentro da faixa
de permissividades elétricas relativas são mais desejáveis para circuitos em
microondas, pois dessa forma, irão gerar menos campos nas bordas e menos efeitos
de acoplamento, possuindo dimensões bem pequenas [5].
22
3.3 ANTENA DE MICROFITA CIRCULAR

A montagem de uma antena de microfita circular segue as mesmas regras da
retangular, entretanto, o seu raio r é a principal dimensão a ser determinada.
centro
ponto de
alimentação
substrato dielétrico antena
r
h
plano terra
conector coaxial
Figura 3.2-Antena de microfita circular com alimentação coaxial.
3.4 MÉTODO DE ALIMENTAÇÃO

Existem várias configurações de alimentação para antenas de microfita. As
quatro mais utilizadas são através de uma linha de transmissão, por acoplamento de
abertura, por acoplamento de proximidade e via cabo coaxial [5].
Neste trabalho a configuração empregada para se alimentar as antenas é via
cabo coaxial, conforme se pode observar na figura 3.2.
23
A alimentação via cabo coaxial, onde o condutor interno do cabo é ligado à

antena e o condutor externo fica acoplado ao plano terra, é comumente usada para
se alimentar antenas de microfita.
3.5 IMPEDÂNCIA DE ENTRADA

Denomina-se impedância de entrada de uma antena a impedância que esta
apresenta à linha de transmissão ou à estrutura de acoplamento através da qual ela
é alimentada [13]. Seu conhecimento é de fundamental importância, pois a eficiência
com que se efetua a transferência de energia do transmissor para a antena, ou da
antena para o receptor acha-se diretamente ligada à impedância da antena.
3.5.1 Equacionamento
Reatância indutiva do
condutor de alimentação
Resistência da fonte
Xf j Lf
Z
Impedância da antena
I
V R
Za Ra jX a
Figura 3.3-Circuito equivalente para o cálculo da impedância de entrada.
O circuito equivalente mostrado na figura 3.3 pode ser analisado através da

seguinte equação:
V RZ R ( Ra jXa j Lf )
, (3.1)
I R Z R Ra jXa j Lf
onde I é a transforma de Fourier da corrente de excitação da antena, que é um pulso

Gaussiano de corrente, V é a transformada de Fourier da integral do campo elétrico
24
ao longo da fonte de excitação, Z é a associação em série da impedância da antena

(Za) com a reatância indutiva da alimentação (Xf), R é a resistência linear atribuída à
fonte com a finalidade de fazer decair rapidamente para zero a tensão gerada e Lf é
a indutância oriunda do condutor interno do cabo coaxial, gerando uma reatância de
alimentação, X f .
3.5.2 Diâmetro do condutor interno do cabo coaxial

Para se estimar o diâmetro do condutor interno do cabo coaxial, devemos
observar a seguinte equação [5]:
Xf 2
0,577
kh
4e
d , (3.2)
k
onde k é a constante de propagação da onda, é a impedância
intrínseca do meio , h a espessura do substrato dielétrico e d é o diâmetro do

condutor interno do cabo coaxial.
Na figura 3.4, pode-se observar que a resistência máxima(Rmax) irá ocorrer na
freqüência de ressonância, que vale 1,225 GHz. X f é a reatância indutiva inserida
pelo condutor central do cabo coaxial.

25
40
pico de resistência na
Resistência, R ressonância (Rmax)
Resistência, reatância (ohms) Xmax
20
Xf
Reatância, X
Xmin
-20
1,200 1,225 1,250
Freqüência (GHz)
Figura 3.4- Variação típica da resistência e da reatância em função da freqüência.
Para a obtenção da reatância de alimentação X f a partir do gráfico da figura
3.4, utilizaremos a seguinte equação:
X min X max
Xf , (3.3)
2
onde Xmin e Xmax são, respectivamente, os valores mínimo e máximo da curva de

reatância, conforme visto no gráfico da figura 3.4.
Para X f em função da indutância de alimentação temos:
Xf j Lf . (3.4)
26
3.5.3 Perda de retorno

A perda por reflexão na conexão do cabo coaxial de alimentação com a antena
pode ser calculada a partir da seguinte equação:
Z Zf
RL 20 Log , (3.5)
Z Zf
onde Z é a impedância de entrada da antena, dada pela equação (3.1) e Z f é a
impedância do equipamento de medição, à qual será atribuído o valor de 50 , que é

o valor padrão dos equipamentos de medição comumente utilizados.
3.6 EFEITO DE BORDA

Como as dimensões da antena de microfita são finitas em comprimento e
largura, os campos elétricos nas bordas sofrem um espalhamento para fora dos
limites laterais da antena. Esse fenômeno é conhecido como efeito de borda [5],
como mostrado na figura 3.5(b). L representa uma dimensão elétrica, isto é, na
realidade a antena de microfita se comportará, eletricamente, como se suas
dimensões fossem maiores do que realmente são, pois suas linhas de campo
elétrico se espalham para fora dos limites das bordas da antena metálica, surgindo
uma constante dielétrica efetiva reff.
27
(a) vista superior
substrato dielétrico antena metálica

linhas de campo
h r
(b) vista lateral plano terra
Figura 3.5- Dimensões física e efetiva de uma antena de microfita retangular.
A intensidade do efeito de borda é função das dimensões (L e W) da antena e

da altura do substrato dielétrico (h), observados na figura 3.5(a) e (b). Desde que
L/h 1 [5], o efeito de borda é bastante reduzido; entretanto deve ser levado em
consideração pois ele modificará a freqüência de ressonância da antena. O mesmo
se aplica à dimensão W.
O resultado deste efeito de borda se faz sentir na freqüência de ressonância da
antena, quanto maior o número de linhas de campo concentradas no dielétrico,
maior será a freqüência de ressonância. Isto quer dizer que, à medida que o efeito
de borda aumenta, a freqüência de operação da antena diminui [5].
Utilizando o gráfico da figura 3.4 e a equação 3.6, podemos estimar o valor de
L para uma antena de microfita circular a partir da freqüência de operação
(ressonância) da antena, utilizando a seguinte equação [5]:
1 1
L 2L , (3.6)
2 f0 reff 0 0
28
onde f0 é a freqüência de ressonância da antena extraída do gráfico da resistência
da antena, L é o incremento no comprimento da antena devido ao efeito de borda e
0e 0 são, respectivamente, a permeabilidade magnética e constante dielétrica do
vácuo.
Para obtenção de reff, temos que [5]:

1/ 2
r 1 r 1 h
reff 1 12 , (3.7)
2 2 W
onde reff é a constante dielétrica efetiva do material dielétrico devido ao efeito de

borda, r é a constante dielétrica física do material utilizado, h é a espessura do
dielétrico e W a largura da antena.
Para a antena de microfita circular utilizaremos a seguinte equação [18] para
estimarmos o seu raio efetivo:
1/ 2
2h r
reff r 1 ln 1, 7726 , (3.8)
r r 2h
onde reff é o raio efetivo da antena circular, considerando seu efeito de borda e h é a
espessura do dielétrico.
3.7 DIAGRAMA DE RADIAÇÃO

Uma antena apresenta direções preferenciais de irradiação, que são direções
segundo as quais, para uma mesma potência total irradiada é possível registrar
diferentes quantidades de campo elétrico, a uma mesma distância r do centro de
coordenadas.
Os diagramas de radiação são representações gráficas das propriedades de
radiação de uma antena em função de coordenadas espaciais. Usualmente, estes
diagramas tridimensionais são calculados, medidos e registrados em uma série de
diagramas bidimensionais. De uma forma ainda mais simples, o desempenho da
antena será descrito em termos de diagramas bidimensionais segundo alguns
planos principais [13].
29
Na maioria das vezes os diagramas de radiação são determinados na região

de campo distante e representados em função de coordenadas direcionais [5].
3.7.1 Equacionamento
A figura 3.6 mostra o sistema de coordenadas utilizado para o equacionamento
do diagrama de radiação da antena. O desenvolvimento será feito em termos de
campo distante, considerando a fonte de campo elétrico fora da origem. Na figura
observa-se o ponto P, onde será analisado o campo elétrico; nessa figura R é a
distância do elemento radiador até o ponto P, r é a distância da origem até P, r é a
distância da origem até o elemento radiador, é o ângulo entre a projeção de r no
plano xy e o eixo x, é o ângulo r e o eixo x, é o ângulo entre r e o eixo z e éo
ângulo entre r e o eixo z.
z
P
y
r'
Figura 3.6- Sistema de coordenadas para obtenção do diagrama de radiação.
Para o desenvolvimento matemático das expressões que irão gerar os

diagramas de radiação das antenas de microfita estudadas neste trabalho,
utilizaremos as seguintes equações:
30
jkR
e
A JS ds ' , (3.9)
4 s' R
E j A. (3.10)
A equação (3.9) define uma função auxiliar para o cálculo do campo elétrico
radiado pela antena, o vetor potencial magnético A [5]. Para o cálculo de A ,
utilizaremos o vetor densidade de corrente elétrica na superfície da antena J S e a
constante de propagação da onda k, onde k 0 0 .
A equação (3.10) resume a relação entre o vetor potencial magnético e o

campo elétrico E [5] para análise do campo distante da antena de microfita.
Sabendo-se que (x,y,z) são as coordenadas cartesianas do ponto de análise
do campo elétrico e (x ,y ,z ) são as coordenadas cartesianas da fonte de geradora
do campo elétrico, observando a figura 3.6 , temos que:
R (x x' ) 2 (y y' ) 2 (z z' ) 2 . (3.11)
Sabendo-se que em coordenadas esféricas

x = r sen cos , (3.12)
y = r sen sin , (3.13)
z = r cos . (3.14)
Desenvolvendo a equação (3.11) em coordenadas esféricas, observando na
figura 3.6 que ' e considerando r r (campo distante), temos que:

2
r'
R r 1 sen cos( ') . (3.15)
r
Substituindo a equação (3.15) na equação (3.9), temos
e jkr
A 0
J S e jkr 'sen cos( ')
ds ' . (3.16)
4 r s'
Sabendo que J s é um vetor de corrente elétrica atuante na superfície da
antena de microfita, podemos fazer:

JS Jxx Jyy. (3.17)
31
Em coordenadas esféricas
J J x cos cos J y sen cos , (3.18)
J J x sen J y cos . (3.19)
Sabendo-se que as fontes elementares (células superficiais) da antena

encontram-se no plano xy, sendo essas fontes elementares polígonos, tendo como
coordenadas de seu centro xi e yi e área Si, considerando a distância do centro
desses polígonos ao centro de coordenadas igual a ri, podemos, então, substituir
(3.18) e (3.19) em (3.16), obtendo-se as expressões de A e A .
N jkr
0e
A Si ( J xi cos cos J yi sen cos )e jkri sin cos( i )
, (3.20)
i 1 4 r
N jkr
0e
A Si ( J xi sen J yi cos )e jkri sin cos( i)
. (3.21)
i 1 4 r
Considerando que J xi e J yi são constantes no interior de cada uma de suas
células e substituindo as equações (3.20) e (3.21) em (3.10), aplicando algumas

propriedades trigonométricas no resultado obtido, e ainda considerando o plano terra
como sendo infinito, temos as expressões das componentes de campo elétrico nas
direções e :
N jkr
j 0e
E AF Si ( J xi cos cos J yi sen cos )e jk sen ( xi cos yi sen )
, (3.22)
i 1 4 r
N jkr
j 0e
E AF Si ( J xi sen J yi cos )e jk sen ( xi cos yi sen )
, (3.23)
i 1 4 r
AF 2 j sen(kh cos ) , (3.24)
onde N é o numero de células em que foi dividida a superfície da antena de

microfita, AF é o fator de correção para se realizar o cálculo considerando o efeito do
plano terra infinito e sendo h a espessura do substrato dielétrico.
32
CAPÍTULO 4
Resultados
4 RESULTADOS PARA A ANTENA DE MICROFITA RETANGULAR
4.1 MONTAGEM DA ESTRUTURA

A fim de comparar os resultados obtidos foi utilizado como referência o artigo
de Alain Reineix e Bernard Jecko [2], onde é apresentada uma antena de microfita
quadrada com 10mm de lado. Neste dispositivo, o substrato dielétrico que separa a
antena do plano terra possui uma constante dielétrica igual a 2,33 e uma espessura
de 1,57mm, conforme pode ser observado na figura 4.1(a) e (b).
10mm
10 mm (a)
P
2.5mm
2.5mm
1,57 mm
r 2, 33 (b)
plano terra
conector coaxial
Figura 4.1- (a) Ponto de alimentação (P) da antena de microfita. (b) Conexão do cabo
coaxial.
33
4.2 ALIMENTAÇÃO
A alimentação da antena de microfita quadrada é feita através de cabo coaxial
conforme mostrado na figura 4.1(b). Para modelar este tipo de alimentação, utilizou-
se um segmento de reta, com resistência linear igual a 200 ligando o plano da
antena ao plano terra. Esse valor de resistência produz um caimento mais rápido
para zero da resposta de tensão na fonte. No segmento de reta que representa a
fonte é aplicado um pulso Gaussiano de corrente de excitação i, com amplitude igual
a 1 A e duração de 0,1 ns, conforme se pode observar na figura 4.2.
Figura 4.2 Pulso Gaussiano de excitação da fonte.
O gráfico do pulso Gaussiano de excitação é obtido através da seguinte

equação:
2
n t 1
8
T 2
i I0 e (4.1)
34
onde I0 = 1A é o valor de pico do pulso, n é o número de passos no tempo, t é o

valor do passo no tempo, conforme equação (2.45) e T é a duração da fonte obtida
da seguinte equação:
1
T (4.2)
f max
sendo f max o limite superior da faixa de freqüências que se pretende excitar na
antena de microfita. Neste trabalho o valor de f max é igual a 10GHz.
4.3 ESPAÇO COMPUTACIONAL

Para o modelamento do espaço computacional utilizamos uma camada de ar
com altura igual a 20mm, partindo-se da superfície da antena até o limite superior do
referido espaço, e uma distância lateral das bordas da antena ao limite lateral do
espaço computacional de 12mm, obtendo-se um volume total de
3
34,0x34,0x21,57mm . A adoção desses valores foi obtida através da montagem de
várias estruturas com diferentes dimensões. Os resultados eram analisados a cada
simulação e comparados com os valores experimentais. As dimensões escolhidas
foram as da estrutura, a partir da qual, não ocorreram modificações nos resultados
obtidos.
A região compreendida entre a antena de microfita e o plano terra foi
preenchida por um substrato dielétrico de espessura igual a 1,57mm. As paredes
externas do espaço computacional são constituídas por superfícies absorventes,
conforme mostrado na figura 4.3. A escolha desses materiais para constituir os
limites externos do espaço computacional visa suprimir, ao máximo, a reflexão da
onda, o que ocasionaria medições errôneas dos parâmetros da antena. A antena e o
plano terra são constituídos de cobre com resistência superficial igual a 10-3 .
35
Figura 4.3- Modelo computacional da estrutura da antena de microfita quadrada.
4.3.1 Discretização do espaço computacional em malha hexaédrica

Para a discretização da estrutura da antena de microfita quadrada em uma
malha hexaédrica estruturada, implementou-se os seguintes passos:
dividiu-se o segmento de reta que representa a fonte de alimentação em 4
partes;
cada lado do quadrado que compõe a antena de microfita foi segmentado em
20 partes;
a altura da camada de ar foi segmentada em 20 partes;
a distância entre os lados do quadrado e o limite do espaço computacional foi
dividida em 12 partes;
Para a antena de microfita quadrada foram adotados os seguintes parâmetros
de discretização:
LMED = 1,8 mm, conforme equação (2.43), t = 6,164.10-5 ns, conforme
equação (2.45) e T = 0,1 ns, conforme equação (4.2).
4.3.2 Discretização do espaço computacional em malha tetraédrica

Para a discretização em malha tetraédrica devemos apenas incluir o valor
médio da aresta , LMED , obtido através da equação (2.43) e deixar que o software
efetue a discretização do espaço computacional. Para a fonte incluiu-se um valor
36
médio de aresta menor, de modo que a mesma ficasse discretizada em 4

segmentos. O valor de LMED adotado para a estrutura foi de 1,5mm e para a fonte foi
de 0,39mm. O passo no tempo é t = 3,09.10-4 ns, conforme equação (2.45) e a
duração da fonte é T=0,1ns, conforme equação (4.2).
4.4 CARACTERÍSTICAS DOS MATERIAIS UTILIZADOS

Os materiais utilizados na simulação da antena de microfita quadrada têm suas
principais características conforme a tabela abaixo:
TABELA 1- Materiais utilizados na antena quadrada
Material Característica Valor

Ar (constante dielétrica) 1,0
r
Substrato dielétrico (constante dielétrica) 2,33

r
Cobre da antena Rs (resistência superficial) 1,0 m
Resistor da fonte Rl (resistência linear) 127,38 /mm
4.5 RESULTADOS
Os resultados foram obtidos através do software MatLab, gerando grande
volume de processamento, criando-se a necessidade de uma CPU com uma boa
capacidade de processamento.
As figuras 4.4 e 4.5 mostram a curva do campo elétrico na fonte da antena de
microfita quadrada com seu espaço computacional discretizado através de malhas
hexaédrica e tetraédrica, respectivamente.
37
Figura 4.4- Campo elétrico na fonte em malha hexaédrica.
Figura 4.5- Campo elétrico na fonte em malha tetraédrica.

38
Para os dois métodos de discretização os gráficos de resistência e reatância

são mostrados nas figuras 4.6 e 4.7, respectivamente e comparados com o resultado
experimental extraído do artigo de Reineix e Jecko [2]. As curvas de resistência e de
reatância são obtidas plotando-se a parte real e imaginária de Z, respectivamente,
na equação (3.1).
Figura 4.6- Resistência da antena de microfita quadrada.

39
Figura 4.7- Reatância da antena de microfita quadrada.
Nas simulações o diâmetro do condutor interno do cabo coaxial foi considerado

infinitesimal, entretanto surgirá um valor de Lf (indutância de alimentação) que
depende da malha gerada próxima ao ponto de alimentação, pois esta malha
causará um diâmetro efetivo para o condutor, originando a indutância Lf.
O diâmetro efetivo para o condutor de alimentação pode ser estimado pela
seguinte expressão [19]:
d eff 0, 27 , (4.3)
onde deff é o diâmetro efetivo do condutor de alimentação e é o comprimento

médio da aresta gerada na região próxima ao ponto de alimentação.
Com isso, a maneira utilizada para se compensar a diferença nos diâmetros
simulado e experimental foi subtraindo um valor de reatância indutiva de Z, esses
valores de reatância corresponderão a um valor de indutância de ajuste.
40
O gráfico de perda de retorno em decibéis (dB) para os dois métodos de

discretização é mostrado na figura 4.8. A linha cheia representa a malha tetraédrica,
a pontilhada a malha hexaédrica e os círculos representam os valores obtidos
experimentalmente.
Figura 4.8- Perda de retorno da antena quadrada sem o ajuste da indutância de

alimentação.
Para geração da curva da figura 4.9, foi feito um ajuste aproximado da

indutância de alimentação, subtraindo um valor de 0,5 nH para a discretização
tetraédrica e 0,25 nH no caso da discretização hexaédrica. Nessa figura linha cheia
representa a malha tetraédrica, a pontilhada a malha hexaédrica e os círculos
representam os valores obtidos experimentalmente.
41
Figura 4.9- Perda de retorno da antena quadrada com o ajuste da indutância de

alimentação.
Para uma melhor visualização e interpretação dos diagramas de radiação

obtidos na simulação da antena de microfita quadrada, utilizaremos como auxílio as
figuras 4.10 e 4.12, em que serão apresentados os planos no espaço em
coordenadas cartesianas onde serão traçados os respectivos diagramas de
radiação. A direção =0 está sobre o eixo z, em seu sentido positivo.
42
Ponto de
alimentação
Figura 4.10- Plano = 0 para construção do diagrama de radiação da antena quadrada.
Na figura 4.11 é mostrado o diagrama de radiação da antena de microfita

quadrada em função de , para = 0 . As curvas geradas para discretização
tetraédrica ( ), hexaédrica ( ) e experimental ( ) são apresentadas. As três
curvas praticamente coincidem.
Figura 4.11- Diagrama de radiação da antena quadrada para =0 .

43
Ponto de
alimentação
45
y
Figura 4.12- Plano = 45 para construção do diagrama de radiação da antena quadrada.
Na figura 4.13 é mostrado o diagrama de radiação da antena de microfita

quadrada em função de , para = 45 . As curvas geradas para discretização
tetraédrica ( ), hexaédrica ( ) e experimental ( ) são apresentadas. Neste caso,
as três curvas também coincidem.
Figura 4.13- Diagrama de radiação da antena quadrada para = 45 .

44
4.6 ANÁLISE DOS RESULTADOS
4.6.1 Efeito da indutância de alimentação

Os argumentos para a obtenção do diâmetro estimado do condutor interno do
cabo coaxial são extraídos do gráfico da figura 4.7, utilizando as equações (3.2) e
(3.3), considerando o meio material como sendo o substrato dielétrico aplicado na
estrutura. Para a obtenção do diâmetro efetivo aplica-se a equação (4.3), utilizando
como tamanho da aresta o comprimento médio dos segmentos que compõem a
fonte. Gera-se, então, a seguinte tabela:
TABELA 2-Efeito da indutância do condutor de alimentação na antena quadrada
Estrutura Diâmetro estimado Diâmetro efetivo

(mm) (mm)
Simulação com discretização hexaédrica 0,41 0,11
Simulação com discretização tetraédrica 0,09 0,11
Resultado experimental 1,06 1,27*
* para o caso experimental considerou-se o valor do diâmetro real.
A tabela 2 compara os diâmetros estimado e efetivo dos condutores de

alimentação nos dois métodos de discretização com o resultado experimental [2].
A estrutura experimental apresentou um diâmetro estimado de 1,06mm, valor
este bem próximo ao real de 1,27mm [2]. Entretanto, as outras estruturas
apresentaram alguma diferença em relação ao experimental, mostrando a existência
de uma grande dependência do diâmetro do condutor com a malha gerada.
A diferença nos diâmetros dos condutores acarreta valores diferentes de Lf e,
conseqüentemente, valores de Xf diferentes, causando modificações nas curvas de
resposta da impedância de entrada e perda de retorno.

As perdas de retorno nos dois métodos de discretização foram mostradas nas
figuras 4.8 (sem correção de Lf) e 4.9 (com correção de Lf), fazendo a comparação
com os valores experimentais [2]. A melhor aproximação entre as curvas aconteceu
no gráfico da figura 4.9, onde foi realizada uma subtração da reatância indutiva do
condutor coaxial, utilizando uma indutância de ajuste. Isso aconteceu porque, após a
45
subtração dos valores de indutância de ajuste, obteve-se uma correção aproximada

nos valores dos diâmetros dos condutores internos do cabo coaxial em cada
simulação.
4.6.3 Análise do efeito de borda

Os efeitos do campo elétrico nas bordas da antena retangular foram analisados
através da variável L estudada na seção 3.6. Os valores de L para os dois
métodos de discretização foram comparados com os valores experimentais [2]. Para
a obtenção dos valores de L foram utilizados como argumentos os valores da
freqüência de ressonância extraídos do gráfico da figura 4.6 e as dimensões da
antena. Com isso, utilizamos as equações (3.6) e (3.7) para gerar a seguinte tabela:
TABELA 3 -Comparação do efeito de borda
Estrutura L (mm)
Simulação com discretização hexaédrica 1,2
Simulação com discretização tetraédrica 1,3
Estrutura experimental 1,0
Observando os resultados contidos na tabela 3, pode-se constatar que o efeito

do campo elétrico nas bordas da antena resultou em valores semelhantes para as
duas formas de discretização, entretanto apresentou uma diferença aproximada de
0,2 a 0,3mm da estrutura experimental. Essa diferença é muito pequena, haja vista
que a antena é um quadrado de 10mm de lado, fazendo com que o erro seja em
torno de 3% da dimensão da estrutura. Com isso, a tabela 3 traduz as diferenças
encontradas no valor da freqüência de ressonância em termos das dimensões da
antena de microfita.
4.6.4 Diagramas de Radiação

Os diagramas de radiação visíveis nos gráficos das figuras 4.11 e 4.13
mostraram coerência entre os dois métodos de discretização e o resultado
experimental [2]. Percebeu-se através dos resultados mostrados nos gráficos que os
diagramas de radiação foram pouco sensíveis às mudanças no método de
discretização, isto é, os resultados foram satisfatórios tanto em malha tetraédrica
46
como em malha hexaédrica. Ainda, os resultados mostram que as aproximações

utilizadas para o cálculo do diagrama de radiação são perfeitamente aceitáveis.
4.6.5 Comparação dos recursos computacionais gastos nas simulações

Como foi citado no início da seção 4.5, para as simulações realizadas neste
trabalho foram necessários computadores com grande poder de processamento.
Quantizou-se, então, esse processamento através da tabela que se segue:
TABELA 4- Comparação dos métodos de discretização para uma antena de microfita

quadrada.
Parâmetro Hexaédrica Tetraédrica Experimental

Freqüência de ressonância 8,4 GHz 8,3 GHz 8,75GHz
Resistência de entrada da antena 205 220 190
Células primárias geradas 46464 76875 _____
Arestas primárias geradas 147600 94349 _____
Memória utilizada 200 Mb 96 Mb _____
Tempo de processamento 3 horas 21 minutos _____
A partir da tabela 4, pode-se observar que apesar de atingir resultados

coerentes e próximos ao valor experimental [2], uma discretização hexaédrica utiliza
muito mais os recursos da CPU e possui um tempo de simulação bastante superior,
apesar de gerar um menor número de células na discretização. Isso pode ser
explicado pelo fato de que na discretização hexaédrica são geradas um número bem
maior de arestas do que na tetraédrica e, além disso, utilizou-se células muito
pequenas nas bordas do elemento radiador.
47
CAPÍTULO 5
Resultados
5 RESULTADOS PARA A ANTENA DE MICROFITA CIRCULAR
5.1 MONTAGEM DA ESTRUTURA

A fim de comparar os resultados obtidos foi utilizado como referência o artigo
de M. Zweki et al [1], onde foi utilizada uma antena de microfita circular com 10mm
de raio. Neste dispositivo, o substrato dielétrico que separa a antena do plano terra
possui uma constante dielétrica igual a 2,55 e uma espessura de 1,524mm,
conforme pode ser observado na figura 5.1.
4.8mm (a)
P
5.2mm
1,524mm (b)
r 2,55
plano terra
conector coaxial
Figura 5.1- (a) Ponto de alimentação (P) da antena de microfita circular. (b) Conexão do
cabo coaxial.
48
5.2 ALIMENTAÇÃO
A alimentação da antena de microfita circular é feita de maneira praticamente
idêntica a antena retangular (ver seção 4.2). Com isso, utilizaremos f max =10GHz
para a obtenção do pulso gaussiano, tal qual o da figura 4.2.

A única alteração feita na alimentação da antena circular foi no valor da
resistência linear do segmento de reta que liga a antena ao plano terra. Nesta
modelagem utilizaremos uma resistência linear igual a 50 para aquele segmento,
pois esse valor de resistência gerou um rápido decaimento para zero da tensão na
fonte.
5.3 ESPAÇO COMPUTACIONAL

Para o modelamento do espaço computacional utilizamos uma camada de ar
com altura igual a 40mm, partindo-se da superfície da antena até o limite superior do
referido espaço, e uma distância lateral das bordas da antena ao limite lateral do
espaço computacional de 20mm, obtendo-se um volume total de 30x60x41,524mm3.
A metodologia utilizada para a obtenção dessas dimensões é idêntica à aplicada a
antena de microfita quadrada vista na seção 4.3. A região compreendida entre a
antena de microfita e o plano terra foi preenchida por um substrato dielétrico de
espessura igual a 1,524mm.
Figura 5.2- Modelo computacional da estrutura da antena de microfita circular

49
Aproveitando-se da simetria da estrutura, pode-se dividir o espaço

computacional em duas partes idênticas, conforme a figura 5.2, reduzindo-se, assim,
o tempo de processamento gasto na simulação. O plano de simetria é constituído de
PMC (perfect magnetic conductor) e as demais superfícies externas do espaço
computacional são constituídas por materiais absorventes, conforme mostrado na
figura 5.2. A antena e o plano terra são constituídos de cobre com resistência
superficial de 10-3 .
Na discretização da estrutura da antena circular utilizou-se a malha tetraédrica,
onde foi adotado o valor de LMED igual a 1,6mm, conforme equação (2.43). A fonte
foi discretizada em 5 segmentos, para isso aplicou-se LMED = 0,30mm em sua
estrutura. O passo é t = 3,385.10-4 ns, conforme equação (2.45) e a duração da
fonte é T=0,1 ns, conforme equação (4.2), respeitando-se, assim, os limites de
dispersão e estabilidade numérica do método vistos na subseção 2.3.1.
5.4 CARACTERÍSTICAS DOS MATERIAIS UTILIZADOS

Os materiais utilizados na simulação da antena de microfita circular têm suas
principais características conforme a tabela abaixo:
TABELA 5- Materiais utilizados na antena circular
Material Característica Valor

Ar (constante dielétrica) 1,0
r
Substrato dielétrico (constante dielétrica) 2,55

r
Cobre da antena Rs (resistência superficial) 1,0 m
Resistor da fonte Rl (resistência linear) 111,548 /mm
PMC Rs (resistência superficial) 1010

50
5.5 RESULTADOS
No gráfico da figura 5.3 é mostrada a forma de onda do campo elétrico na fonte
que alimenta a antena de microfita circular.
Figura 5.3- Campo elétrico na fonte da antena circular.

51
A figura 5.4 mostra a resistência e a reatância da antena de microfita circular

em função de sua freqüência. As curvas de resistência e de reatância são obtidas
plotando-se a parte real e imaginária de Z, respectivamente, a partir da equação
(3.1).
Figura 5.4- Resistência e reatância na antena de microfita circular.

52
A figura 5.5 mostra, para uma antena de microfita circular, a perda de retorno
em decibéis (dB) em função da freqüência sem o ajuste da indutância gerada pelo
condutor de alimentação. O ajuste de Lf será feito da mesma forma utilizada na
antena de microfita quadrada. A linha contínua representa a resposta simulada e os
círculos representam o resultado experimental obtido do artigo de M. Zweki [1].
Figura 5.5- Perda de retorno da antena circular sem o ajuste da indutância de

alimentação.
53
A figura 5.6 mostra a perda de retorno em decibéis (dB) em função da

freqüência com o ajuste da indutância gerada pelo condutor de alimentação. Para
geração da curva da figura 5.6, foi feito um ajuste na indutância de alimentação no
valor de 1,2 nH, o que melhorou os resultados em relação aos valores experimentais
[1].
Figura 5.6- Perda de retorno da antena circular com o ajuste da indutância de

alimentação.
54
Para uma melhor visualização e interpretação do diagrama de radiação obtidos

na simulação da antena de microfita quadrada, utilizaremos como auxílio a figura
5.7, onde é apresentado o plano no espaço em coordenadas cartesianas onde será
traçado o diagrama de radiação. A direção =0 está sobre o eixo z, em seu sentido
positivo.
30
Figura 5.7- Plano = 30 para construção do diagrama de radiação da antena circular.

55
Na figura 5.8 é mostrado o diagrama de radiação da antena circular em função

de , para = 30 . A linha cheia representa o diagrama da antena discretizada em
malha tetraédrica; a linha pontilhada representa o diagrama de radiação obtido
experimentalmente [1]. A diferença existente entre as duas curvas próximo a =
180 , deve-se ao fato do experimento ter plano terra finito.
Figura 5.8- Diagrama de radiação da antena circular para = 30 .

56
5.6 ANÁLISE DOS RESULTADOS
5.6.1 Efeito da indutância de alimentação

Os argumentos utilizados para o cálculo do diâmetro estimado do condutor
coaxial na simulação são extraídos do gráfico da figura 5.4, aplicando as equações
(3.2) e (3.3). Para o cálculo do diâmetro efetivo na simulação foi aplicada a equação
(4.3).
Tabela 6- Efeito da indutância do condutor de alimentação na antena circular
Estrutura Diâmetro estimado (mm) Diâmetro efetivo (mm)

Simulação 0,35 0,08
Resultado experimental 0,41*
* Para o caso experimental considerou-se o valor do diâmetro real.

Observou-se que o diâmetro estimado para o condutor central do cabo coaxial
na simulação correspondeu a um valor próximo ao que foi utilizado na estrutura
experimental [1]. Entretanto, houve uma diferença em relação ao diâmetro efetivo, o
que motivou o ajuste da indutância nessa estrutura e evidenciou a dependência
entre o valor do diâmetro e a malha gerada.

As perdas de retorno foram mostradas nas figuras 5.5 (sem o ajuste de Lf) e
5.6 (com o ajuste de Lf), fazendo a comparação com os valores experimentais [1]. A
melhor aproximação entre as curvas aconteceu no gráfico da figura 5.6, onde foi
realizada uma subtração da reatância indutiva do condutor coaxial, conforme visto
na seção 5.6.2. Observa-se na figura 5.3, que o valor de tensão gerada na fonte
decai para zero rapidamente devido à resistência linear de 50 inserida na
modelagem da alimentação.
5.6.3 Análise do efeito de borda

O efeito do campo elétrico na borda da antena circular foi analisado através da
variável reff estudada na seção 3.6. Para a obtenção do valor de reff foram utilizados
como argumentos os valores de h, r e r, aplicados na equação (3.8). Com isso
obteve-se o valor de reff = 10,8 mm, obtendo-se, assim um r, aproximado, de
0,8mm. Percebe-se que a variação do raio foi muito pequena em relação à dimensão
57
da antena circular, fazendo com que o valor da freqüência de operação da antena

não variasse muito em relação ao experimental.
5.6.4 Diagrama de radiação

O diagrama de radiação simulado para uma antena de microfita circular é
mostrado no gráfico da figura 5.8, onde é comparado com a diagrama obtido
experimentalmente. No caso da antena de microfita circular percebe-se que as
curvas se ajustam em praticamente todos os pontos de seu traçado, porém, próxima
à região de =180 , as curvas apresentam algumas diferenças em seu traçado
devido ao fato de o plano terra ser finito para o caso experimental e infinito na
simulação computacional.
58
CAPÍTULO 6
Conclusões e Trabalhos futuros
6 CONCLUSÃO
Duas maneiras de discretização foram vistas neste trabalho, sendo que a
discretização em malha hexaédrica mostrou melhores resultados para estruturas
retangulares, porém exigindo muito mais recursos da CPU. Uma nova maneira de
discretização, a tetraédrica, foi introduzida para ser aplicada a estruturas de formato
não retangular. A efetividade desse tipo de discretização foi comprovada através de
sua aplicação a uma estrutura retangular, fazendo a comparação de resultados com
outro método de discretização já conhecido e com resultados experimentais.
Posteriormente, já constatada sua efetividade, a discretização em malha tetraédrica
foi aplicada a uma estrutura circular.
A escolha da maneira de discretização vai depender da precisão desejada e
dos computadores disponíveis para a realização processamento da simulação. Uma
má escolha pode ocasionar em horas de processamento e até mesmo, na
impossibilidade de realização da simulação com um resultado confiável.
Observou-se também que, para a simples construção de diagramas de
radiação, a maneira de discretização escolhida pouco influenciará no resultado. Isso
foi observado quando se analisou a antena de microfita quadrada.
Analisando-se a antena de microfita quadrada, pode-se observar pequenas
discrepâncias nos resultados quando comparados aos valores experimentais. Isso
se deve ao fato de que o efeito de campo elétrico nas bordas da antena e do
diâmetro do condutor central do cabo coaxial influenciarem no resultado, alterando
sua impedância de entrada e freqüência de ressonância.
59
A antena circular, ao ser discretizada em malha tetraédrica, mostrou bons

resultados, tendo em vista que sua geometria favorece a aplicação desse método de
discretização. Entretanto, apesar dos resultados obtidos terem sido melhores do que
os obtidos na antena quadrada, foram, da mesma forma, influenciados pelo efeito de
borda e pela indutância de alimentação.
Utilizando-se o método FDTD com as duas formas de discretização, pode-se,
ainda, simular computacionalmente o comportamento elétrico de uma antena de
microfita antes de implementá-la fisicamente. Entretanto, para que essa simulação
forneça resultados confiáveis, far-se-á necessário um modelamento especial,
utilizando o método FDTD [19], para a fonte que alimentará a antena de microfita.
Com isso, foi mostrado neste trabalho que o método FDTD com discretização
da estrutura em malha tetraédrica é perfeitamente viável para a análise de
parâmetros de antenas de microfita com formato retangular ou circular,
principalmente para este último, onde a geometria da estrutura se adapta bem às
células geradas pela discretização.
60
6.1 TRABALHOS FUTUROS

Como trabalho futuro, sugere-se o estudo de antenas conformais [14], tendo
em vista a grande importância de alocarmos antenas de microfitas sobre superfícies
não planas, como por exemplo, o teto de um automóvel ou partes frontais de
aeronaves.
Figura 6.1- Antena de microfita sobre superfície cônica.
Para se estudar esse modelo de estrutura far-se-á necessária a mudança do

sistema de coordenadas utilizadas neste trabalho. Recomenda-se o método de
discretização através de uma malha tetraédrica, tendo em vista a forma curva da
estrutura.
61
ANEXO I
Algoritmo utilizado para o cálculo da impedância de entrada da antena de
microfita retangular
load minhavar;
%resistencia da fonte (ohms)
R=200;
L=1.57/4;
%calculo do V
V=(Ez(:,1)-Ez(N,1)+Ez(:,2)-Ez(N,2)+Ez(:,3)-Ez(N,3)+Ez(:,4)-Ez(N,4))*L;
%soma as tres colunas de Ez
%calculo de I
I=pulse(1:N,dt,T)';
%fft de I
fftI=(fft(I,2^20));
%fft do V
fftV=(fft(V,2^20));
%calculo Za
Za=((fftV.*R)./(fftI.*R-fftV));
df=(1/dt)/2^20; %passo da freq

f=[0:2^20-1]*df; %range da freq
%figure(1);
%plot(f,abs(fftV))
%axis([0 50 0 70])
figure(3);
plot(f,real(Za),'r',f,imag(Za),'b')
62
axis([0 15 -100 220])

grid on
%figure(4);
%plot(f,imag(Za))
%axis([0 30 -30 30])
%grid on
%Resultados experimentais Paper Jecko [2]

Fj=[7 7.25 7.5 7.75 8 8.25 8.5 8.75 9 9.25 9.5 9.75 10 10.25 10.5 10.75];
Rj=[9 8 9 19 28 61 134 190 128 67 40 29.5 21.5 20 18.5 18.6];
Xj=[48 52 65 79 91 116 127 20 -55 -45 -30 -15 -5 0 6 15];
%Comparacao com os resultados experimentais

figure(5)
plot(f,real(Za),'r',Fj,Rj,'b+')
axis([7 12 -10 220])
grid on
figure(6)
plot(f,imag(Za),'r',Fj,Xj,'b+')
axis([7 12 -100 220])
grid on
63
ANEXO II
Algoritmo utilizado para o cálculo da impedância de entrada da antena de
microfita circular
plot(ft,real(Zac)/2,'r',ft,imag(Zac)/2,'b')
axis([0 12 -100 250]) clear all;
close all;
load varIE;
%resistencia da fonte (ohms)
R=50;
L=1.524/5;
%calculo do V
V=(Ez(:,1)-Ez(N,1)+Ez(:,2)-Ez(N,2)+Ez(:,3)-Ez(N,3)+Ez(:,4)-Ez(N,4)+Ez(:,5)-
Ez(N,5))*L; %soma as tres colunas de Ez
%calculo de I
I=pulse(1:N,dt,T)';
%fft de I
fftI=(fft(I,2^20));
%fft do V
fftV=(fft(V,2^20));
df=(1/dt)/2^20; %passo da freq
ft=[0:2^20-1]*df; %range da freq
%calculo Za
Za=((fftV.*R)./(fftI.*R-fftV));
%calculo de Za corrigido
Zac=((fftV.*R)./(fftI.*R-fftV))-j*2*pi*ft'*(1.3);
%figure(1);
%plot(f,abs(fftV))
%axis([0 50 0 70])
Zf=50; % imp da fonte de medicao
fg=[2:0.1:12];% freq extraida do graf do paper - experimental
64
%Loss em dB do graf do paper experimental [1]

Ldbmed=[0 0 0 0 -0.15 -0.20 -0.21 -0.22 -0.23 -0.24 -0.25 -0.26 -0.27 -0.28 -
0.29 -0.30 -0.31 -0.32 -0.33 -0.34 -0.36 -0.44 -0.56 -0.65 -0.72 -1.08 -1.44 -1.81 -2.60
-3.70 -5.10 -7.16 -8.70 -8.20 -6.45 -4.86 -3.60 -2.90 -1.40 -1.28 -1.15 -1.13 -1.10 -
1.08 -1.07 -1.07 -1.05 -1.05 -1.04 -1.04 -1.06 -1.06 -1.06 -1.08 -1.10 -1.12 -1.15 -1.35
-1.50 -1.65 -1.85 -2.52 -3.25 -4.65 -6.80 -8.98 -8.95 -6.45 -5.00 -3.98 -3.30 -3.25 -
3.23 -3.22 -3.50 -3.70 -3.96 -6.72 -7.12 -8.70 -11.72 -13.30 -16.80 -16.08 -12.88 -
10.10 -8.06 -6.44 -5.38 -4.68 -3.95 -3.40 -3.00 -2.95 -2.70 -3.00 -3.50 -3.84 -4.32 -
4.10 -3.78];
% Loss calculado na simulacao
LdB=20*log10(abs((Za-Zf)./(Za+Zf)));
%valor corrigido
LdBc=20*log10(abs((Zac-Zf)./(Zac+Zf)));
figure(3);
grid on
figure(4);
plot(fg,Ldbmed,ft,LdB,'r');
axis([0 12 -50 10])
grid on
figure(5);
plot(fg,Ldbmed,ft,LdBc,'r');
axis([0 12 -50 10])
grid on
65
ANEXO III
Algoritmo utilizado para o cálculo do diagrama de radiação da antena de
microfita retangular
% BLOCO DO CALCULO HEXAEDRICO
clc
close all
%load s2d;carrega a matriz dos nos dos elementos da malha
%load P;carrega a matriz das coordenadas dos nos do projeto
load mvarh;
[n,c]=size(s2d);%n=numero de linhas de s2d(duracao sim/passo)
Sir=zeros(n,1);%matriz das areas dos retangulos
Crx=zeros(n,1);%matriz das coordenadas(x) dos centros dos elementos da
malha
Cry=zeros(n,1);%matriz das coordenadas(y) dos centros dos elementos da
malha
teta=[-pi/2:pi/180:pi/2];%faz a variacao de teta de 0 a 90 graus em rad
[m,l]=size(teta);
%N=fix(TT/dt);
input('Frequencia de ressonancia Hexaedrica(em GHz): ');
fo=fo*10^9;
input('Angulo Fi em graus: ');
Fi=Fi*pi/180;
w=2*pi*fo;
k=w*(muo*epo)^(1/2);
%altura do dieletrico
h=1.57;
for(i=1:n)
%obtem da matriz s2d o numero do no correspondente de cada elemento da
%malha
ni1=s2d(i,1);
ni2=s2d(i,2);
ni3=s2d(i,3);
ni4=s2d(i,4);
%obtem da matriz P as coordenadas(x,y,z) do no correspondente de cada

elemento da
%malha e joga na variavel Pn1
Pni1=P(ni1,:);
Pni2=P(ni2,:);
Pni3=P(ni3,:);
66
Pni4=P(ni4,:);
%obtem as coordenadas x,y apenas, de cada no e lanca na matrix Mi

xi1=Pni1(1,1);
yi1=Pni1(1,2);
xi2=Pni2(1,1);
yi2=Pni2(1,2);
xi3=Pni3(1,1);
yi3=Pni3(1,2);
xi4=Pni4(1,1);
yi4=Pni4(1,2);
Mi=[xi1 yi1 1;xi2 yi2 1;xi3 yi3 1];
%calculo da area do retangulo
Sir(i,1)=abs(det(Mi));
%calculo das coordenadas xi/yi do centros dos retang da malha
Crx(i,1)=[(xi1+xi2+xi3+xi4)/4];
Cry(i,1)=[(yi1+yi2+yi3+yi4)/4];
end
%criar matriz auxiliar para obtencao de Sir, Cx, Cy completos

mat_aux=ones(1,l);
Sircomp=Sir*mat_aux;
Crxcomp=Crx*mat_aux;
Crycomp=Cry*mat_aux;
%criar matriz auxiliar para obtencao de teta completo
mat_aux2=ones(n,1);
tetacomp=mat_aux2*teta;
%calculo DFT de Ex e Ey
%e=exp(-j*2*pi*fo*[1:N]*dt);
Exfo=Ex;
Eyfo=Ey;
%criar matriz de Exi e Eyi
Exi=Exfo'*mat_aux;
Eyi=Eyfo'*mat_aux;
67
% Fator de Array para considerar efeito do plano terra
AF=2*j*sin(k*h.*cos(tetacomp));
%Calculo Eteta e Efi
Eteta=AF.*Sircomp.*(Exi*cos(Fi).*cos(tetacomp)+Eyi*sin(Fi).*cos(tetacomp)).*e
xp(j*k*sin(tetacomp).*(Crxcomp*cos(Fi)+Crycomp*sin(Fi)));
EFi=AF.*Sircomp.*(-
Exi*sin(Fi)+Eyi*cos(Fi)).*exp(j*k*sin(tetacomp).*(Crxcomp*cos(Fi)+Crycomp*sin(Fi)));
%somatorio de Eteta e EFi
Eteta=abs(sum(Eteta));
EFi=abs(sum(EFi));
%Definindo o Plano E
PLE=sqrt(Eteta.^2+EFi.^2);
%Definindo o Plano H
%PLH=sqrt(Eteta.^2+EFi.^2);
%Normalizando Eteta e EFi
%EtetaN=Eteta/max(Eteta);
%EFiN=EFi/max(Eteta);
% BLOCO DO CALCULO TETRAEDRICO
clc
clear all;
close all
load mvart;
malha
malha
teta=[-pi/2:pi/180:pi/2];%faz a variacao de teta de 0 a 90 graus em rad
68
[m,l]=size(teta);
%N=fix(TT/dt);
input('Frequencia de ressonancia do GID Tetraedrica(em GHz): ');
fo=fo*10^9;
input('Angulo Fi em graus: ');
Fi=Fi*pi/180;
w=2*pi*fo;
h=1.57e-3;
for(i=1:n)
%malha
ni1=s2d(i,1);
ni2=s2d(i,2);
ni3=s2d(i,3);
%ni4=s2d(i,4);

elemento da
Pni1=P(ni1,:);
Pni2=P(ni2,:);
Pni3=P(ni3,:);
%Pni4=P(ni4,:);

xi1=Pni1(1,1);
yi1=Pni1(1,2);
xi2=Pni2(1,1);
yi2=Pni2(1,2);
xi3=Pni3(1,1);
yi3=Pni3(1,2);
%xi4=Pni4(1,1);
%yi4=Pni4(1,2);
%calculo da area do triangulo
Sir(i,1)=abs(det(Mi))/2;
Crx(i,1)=[(xi1+xi2+xi3)/4];
69
Cry(i,1)=[(yi1+yi2+yi3)/4];
end

mat_aux=ones(1,l);
mat_aux2=ones(n,1);
Exfo=Ex;
Eyfo=Ey;
Exi=Exfo'*mat_aux;
Eyi=Eyfo'*mat_aux;
EFi=AF.*Sircomp.*(-
Eteta=abs(sum(Eteta));
EFi=abs(sum(EFi));
70
PLEt=sqrt(Eteta.^2+EFi.^2);
%PLHt=sqrt(Eteta.^2+EFi.^2);
%Normalizando PLE e PLH
PLNt0=PLEt/max(PLEt);
%PLHNt0=PLHt/max(PLEt);
71
ANEXO IV
Algoritmo utilizado no cálculo do diagrama de radiação da antena de
microfita circular
clc
clear all
close all
load varDR;
malha
malha
%teta=[-pi/2:pi/180:pi/2];%faz a variacao de teta de 0 a 90 graus em rad
teta=[-pi/2:pi/180:pi/2];
[m,l]=size(teta);
fo=input('Frequencia de ressonancia (em GHz): ');
fo=fo*10^9;
Fi=input('Angulo Fi em graus: ');
Fi=Fi*pi/180;
w=2*pi*fo;
h=1.524;
for(i=1:n)
%malha
ni1=s2d(i,1);
ni2=s2d(i,2);
ni3=s2d(i,3);
%ni4=s2d(i,4);

elemento da
Pni1=P(ni1,:);
Pni2=P(ni2,:);
Pni3=P(ni3,:);
%Pni4=P(ni4,:);
72

xi1=Pni1(1,1);
yi1=Pni1(1,2);
xi2=Pni2(1,1);
yi2=Pni2(1,2);
xi3=Pni3(1,1);
yi3=Pni3(1,2);
%xi4=Pni4(1,1);
%yi4=Pni4(1,2);
%calculo da area do triangulo
Sir(i,1)=abs(det(Mi))/2;
Crx(i,1)=[(xi1+xi2+xi3)/4];
Cry(i,1)=[(yi1+yi2+yi3)/4];
end

mat_aux=ones(1,l);
mat_aux2=ones(n,1);
Exfo=Ex;
Eyfo=Ey;
Exi=Exfo'*mat_aux;
Eyi=Eyfo'*mat_aux;

73
EFi=AF.*Sircomp.*(-
Eteta=sum(Eteta);
EFi=sum(EFi);
PLEt=sqrt(Eteta.^2+EFi.^2);
PLHt=sqrt(Eteta.^2+EFi.^2);
%Normalizando PLE e PLH
PLENt=PLEt/max(PLEt);
PLHNt=PLHt/max(PLEt);
74
REFERÊNCIAS BIBLIOGRÁFICAS
[1] M.ZWEKI, R.A.; ABD-ALHAMEED, M.A.; MANGOUD, P.S.: Broadband

Analysis of Finite Microstrip Patch Antenna Structure using FDTD. 11th
International Conference on Antennas and propagat., UMIST, Manchester,U.K.,
pp.17-20, April 2001.
[2] REINEIX, A.; JECKO,B.: Analysis of Microstrip Patch Antennas Using

Finite Difference Time Domain Method. IEEE Trans. on Antennas and Propagat.,
vol.37,nº 11, pp. 1361-1369, November 1989.
[3] G. A. DESCHAMPS: Microstrip Microwave Antennas, 3rd USAF Symp. On

Antennas, 1953.
[4] J. Q. HOWELL: Microstrip Antennas, in Dig. Int. Symp. Antennas

Propagat. Soc., Williamsburg, VA, pp. 177-180, 1972.
[5] C. A. BALANIS, Antenna Theory: Analysis and Design, Second edition.

New York: John Wiley & Sons, Inc, 1997.
[6] C. A. BALANIS, Advanced Engineering Electromagnetics, Second

edition. New York: John Wiley & Sons, Inc, 1989.
[7] H. GUTTON E G. BAISSINOT, Flat Aerial Ultra High Frequency, French

Patent nº 703113, 1955.
[8] R. E. MUNSON, Single Slot Cavity Antennas Assembly, U.S. Patent nº

3713162, Jan. 23, 1973.
[9] H. D. WEINSHEL, Progress Report on Development of Microstrip

Cylindrical Arrays for Sounding Rockets, Phisic and Sci. Lab., New Mexico State
Univ., Las Cruces, 1973.
[10] YEE, K. S.: Numerical Solution of Initial Boundary Value Problems

Involving Maxwell s Equations in Isotropic Media. IEEE Trans. Antennas
Propagat., vol. 14, pp. 302-307, May 1966.
[11] WILSON A. ARTUZI JR.: Finite Element Time Domain Method Using
Piecewise Constants Basis Functions. SBMO/IEEE International Microwave and
Optoelectronics Conference, vol.2, pp. 1029-1032, Foz do Iguaçu, 2003.
[12] A. TAFLOVE, Computational Electrodynamics: The Finite Difference

Time-Domain Method, Artech House, Boston, 1995.
[13] LUÍS CLÁUDIO ESTEVES, Antenas: Teoria Básica e Aplicações, São

Paulo: McGraw-Hill do Brasil, 1980.
75
[14] ] W. Y. TAM; A. K. Y. LAI; K. M. LUK.: Input Impedance of Spherical

Microstrip Antenna. IEE Proc.- Microw. Antenas e Propag., vol.142,nº 3, pp. 285-
288, June 1995.
[15] P. MYTHILI; A. DAS.: Simple Approach to Determine Resonant

Frequencies of Microstrip Antennas. IEE Proc.- Microw. Antenas e Propag.,
vol.145,nº 2, pp. 159-162, April 1998.
[16] K. F. WARNICK; R. H. SELFRIDGE; D. V. ARNOLD.: Teaching

Electromagnetic Field Theory Using Differential Forms. IEEE Trans. on
Education, vol. 40, pp. 53-68, February 1997.
[17] KARL S. KUNZ; RAYMOND J. LUEBBERS, The Finite Difference Time

Domain Method for Eletromagnetics, New York: CRC Press LLC, 1993.
[18] A. A. KISHK; LOTFOLLAH SHAFAI.: The effect of various parameters of

circular microstrip antennas on their radiation efficiency and the mode
excitation. IEEE Trans. on Antennas and Propagat., vol.34, nº 8, pp. 969-976,
August 1989.
[19] SO-ICHI WATANABE; M. TAKI.: An Improved FDTD Model for the

Feeding Gap of Thin-Wire Antenna. IEEE Microw. and Guided Wave Letters, vol.
8, nº 4, April 1998.
This document was created with Win2PDF available at http://www.daneprairie.com.
The unregistered version of Win2PDF is for evaluation or non-commercial use only.

Antenassss

Enviado por

Direitos autorais:

Formatos disponíveis

Antenassss

Enviado por

Dados do documento

Título original

Direitos autorais

Formatos disponíveis

Compartilhar este documento

Compartilhar ou incorporar documento

Opções de compartilhamento

Você considera este documento útil?

Este conteúdo é inapropriado?

Direitos autorais:

Formatos disponíveis

Antenassss

Enviado por

Direitos autorais:

Formatos disponíveis

ii

O objetivo principal deste trabalho é a utilização do método FDTD (Diferenças

Palavras chave: FDTD, antena de microfita.

Keywords: FDTD, microstrip antenna.

FIGURA 1.1- ANTENA DE MICROFITA COM ALIMENTAÇÃO COAXIAL................. 1

FIGURA 4.6- RESISTÊNCIA DA ANTENA DE MICROFITA QUADRADA. .............. 38

TABELA 1- MATERIAIS UTILIZADOS NA ANTENA QUADRADA ........................... 36

c Velocidade da luz no vácuo

J Vetor densidade de corrente na direção

J Vetor densidade de corrente na direção

xi coordenada x do centro da célula superficial

Figura 1.1- Antena de microfita com alimentação coaxial.

O conceito de Antena de Microfita surgiu há aproximadamente cinqüenta anos

milhões de cálculos em frações de segundos, tornando possível a implementação do

2.2 O MÉTODO FDTD

Figura 2.1- Estrutura da célula tetraédrica.

Na figura 2.1, os pontos 1,2,3 e 4 compõem a malha primária (linhas sólidas) e

2.2.1 Grandezas eletromagnéticas e formas diferenciais

onde E é o vetor campo elétrico, B é o vetor fluxo magnético, H é o vetor campo

Sendo E E1 x E2 y E3 z uma grandeza vetorial dada em v/m, e

dl dx x dy y dz z , um elemento vetorial diferencial de comprimento, dado em

Wb/m2, e ds dzdy x dxdz y dxdy z , um elemento vetorial diferencial de superfície,

As formas diferenciais são contrações do produto escalar de uma grandeza

2.2.2 Discretização das Equações de Maxwell

Bp d p ,ij Eij , (2.14)

são escalares que representam o fluxo magnético na face primária p e a diferença

são +1, -1 e 0, os quais realizam a somatória das quantidades associadas às arestas

índice T indica a matriz transposta.

sendo , e , respectivamente, permissividade elétrica, relutividade magnética e

converte uma forma de grau a em uma forma de grau 3 a , preservando sua

Para análise de valores elétricos subdividiremos a célula em 6 subseções

sendo dLij é um elemento diferencial de comprimento de aresta primária e dSp é um

Substituindo (2.19) em (2.18) temos que

dLij / Lij em Vij

dLij / Lij , em Vij

Eij dWij B pW p , (2.23)

Bp d Wp Eij Wij J ij Wij Eij Wij . (2.24)

Valendo-se do método de Galerkin [5], multiplicaremos ambos os membros da

d q ,ij Eij q, p Bp , (2.25)

d kl , p B p kl ,ij Eij kl ,ij Eij kl ,ij J ij , (2.26)

kl ,ij Wkl Wij , (2.28)

kl Wkl Wij , (2.29)

kl ,ij Wkl Wij , (2.30)

d q ,ij Wq dWij , (2.31)

são constantes conhecidas em função dos parâmetros geométricos e constitutivos

2.2.3 Equação Global da Onda

Para se obter a equação global da onda multiplicaremos a equação (2.25) por

operações, substituiremos a equação (2.25) modificada após a multiplicação por

kl , ij Eij kl ,ij Eij kl ,ij E ij kl ,ij J ij , (2.33)

Podendo-se reescrever a equação (2.33)

que é a Equação Global da Onda, na qual desejamos obter os valores de E nas

2.2.4 Modelagem de superfícies

onde Rs é a resistência superficial do material, utilizando

ABC Rs , onde e são características eletromagnéticas do meio

2.2.5 Modelagem de fios

onde Rl é a resistência linear do material, utilizando os mesmos valores adotados

2.2.6 Discretização temporal

2.3 MODELAGEM COMPUTACIONAL

pois demandará grande quantidade de memória RAM e muito espaço no disco