République Algérienne Démocratique et Populaire

Ministère de l’enseignement supérieur et de la recherche scientifique

Université 8Mai 1945 – Guelma
Faculté des sciences et de la Technologie
Département Electronique et Télécommunications

Mémoire
Présenté pour obtenir le diplôme de Magister en Télécommunications

OPTION: Télécommunications
Par:
BOUSSAHA SOMIA

Thème

Contribution à l'étude et à la modélisation des antennes

imprimées rectangulaires
Applications aux systèmes de communication
Devant le Jury :

Président Pr. TEBBIKH Hicham [Link]

Rapporteur Dr. BOUALLEG Ahcene [Link]
Examinateurs Pr. BOUMAZA Mohamed Seghir [Link]
Dr. REDADAA SALAH [Link]
Invité Dr. TABBA Mohamed Tahar [Link]

SESSION 2012
 ‫ﺑﺴﻢ اﷲ اﻟﺮﺣﻤﻦ اﻟﺮﺣﯿﻢ واﻟﺤﻤﺪ ﷲ رب اﻟﻌﺎﻟﻤﯿﻦ واﻟﺼﻼة واﻟﺴﻼم‬
.‫ﻋﻠﻰ ﻧﺒﯿﻨﺎ ﻣﺤﻤﺪ ﺻﻠﻰ اﷲ ﻋﻠﯿﮫ و ﺳﻠﻢ‬

A travers ce modeste travail, je tiens à remercier vivement notre

encadreur Dr .Boualleg Ahcene pour ses conseils précieux et pour
toutes les commodités et aisances qu’il j’ai apportées durant notre
étude et réalisation de ce projet.
Je remercie Monsieur Mohamed Tebbikh Hicham, Professeur à
l’Université de Guelma d’avoir accepté de présider le jury.
J'adresse mes remerciements à [Link] Mohamed Seghir et le
Dr. Redadaa Salah et aussi Dr. Tabba Mohamed Tahar pour avoir
accepté de juger ce travail.
Sans omettre bien sur de remercier profondément à tous ceux qui ont
contribué de prés ou de loin à la réalisation de ce travail en
particulier Mme .Medjaldi Assia responsable du laboratoire de
circuit imprimée pour leur aide.
Je tiens aussi à remercier les gens du laboratoire de
recherche, (LT) pour leur aide et leurs
judicieux conseils.
 Résumé

La miniaturisation des antennes a suscité ces dernières années un très grand intérêt au prés des chercheurs
dans le monde des télécommunications. La course à cette miniaturisation est motivée par des exigences
d’intégration des antennes dans l’architecture des terminaux mobiles, pour réduire au maximum leur
encombrement. Cette structure convient particulièrement bien aux applications de pointes, comme les
communications mobiles cellulaires ou la réception de signaux de télédiffusion satellite. L’une des applications
actuelles qui nécessite aussi l’utilisation d’antennes de petites dimensions est celle des interconnections sans
fils.
Parmi les différentes techniques de miniaturisation des antennes, nous avons exploré dans ce travail une
étude générale sur les antennes imprimée .Ce type d'antenne planaire s'utilise, soit sous la forme d'antenne à
élément unique, soit sous la forme de réseaux où plusieurs antennes sont groupées sur un même substrat.
Nous avons évalué les performances de l’antenne imprimée rectangulaire par le HFSS pour les différant
méthodes d’alimentation. Enfin, nous avons expliqué le phénomène de couplage, qu’est un critère important
qu’il faut donc considérer lors du calcul des caractéristiques d’un réseau.
 Mots clés:
Antenne imprimée, réseau d’antennes, diagramme de rayonnement, synthèse, optimisation, Méthodes
d’analyses, modèle simulé, simulation par le HFSS.

Abstrat

The miniaturization of antennas has attracted in recent years great interest to researchers in the fields of
telecommunications. The race for miniaturization is driven by the requirements of integration of antennas in
the architecture of mobile devices, to minimize their size. This structure is particularly suitable for advanced
applications such as mobile phones or reception of satellite broadcast signals. One of the current applications
which also require the use of small antennas is that of the son without interconnections.
Among the various techniques of miniaturization of antennas, we explored in this paper a general study of
printed antennas. This type of planar antenna used either as single antenna element, either in the form of
networks where multiple antennas are grouped on the same substrate.
We have modeled the performance of the rectangular patch antenna by HFSS for differing feeding methods.
Finally, we have explained the phenomenon of coupling, which is an important criterion that must be assumed
when calculating the characteristics a network.
 Key words:
Printed antenna, antenna array, radiation pattern, synthesis, optimization, analysis methods, simulated model,
the HFSS simulation.

‫ﻣﻠﺨﺺ‬

‫ اﻟﺴﺒﺎق ﻋﻠﻰ‬.‫ﻗﺪ اﺟﺘﺬب اﻟﺘﺼﻐﯿﺮ ﻣﻦ اﻟﮭﻮاﺋﯿﺎت ﻓﻲ اﻟﺴﻨﻮات اﻷﺧﯿﺮة اھﺘﻤﺎﻣﺎ ﻛﺒﯿﺮا ﻟﻠﺒﺎﺣﺜﯿﻦ ﻓﻲ ﻣﺠﺎﻻت اﻻﺗﺼﺎﻻت اﻟﺴﻠﻜﯿﺔ واﻟﻼﺳﻠﻜﯿﺔ‬
‫ ھﺬا اﻟﮭﯿﻜﻞ ھﻮ ﻣﻨﺎﺳﺒﺔ ﺧﺎﺻﺔ‬.‫ ﻟﺘﻘﻠﯿﻞ ﺑﺼﻤﺘﮭﺎ‬،‫اﻟﺘﺼﻐﯿﺮ اﻟﺘﻲ ﺗﻘﻮم ﻋﻠﻰ ﺿﺮورة اﻟﺘﻜﺎﻣﻞ ﻣﻦ اﻟﮭﻮاﺋﯿﺎت ﻓﻲ اﻟﮭﻨﺪﺳﺔ اﻟﻜﮭﺮﺑﺎﺋﯿﺔ ﻣﻦ اﻷﺟﮭﺰة اﻟﻨﻘﺎﻟﺔ‬
‫ واﺣﺪة ﻣﻦ اﻟﺘﻄﺒﯿﻘﺎت اﻟﺤﺎﻟﯿﺔ اﻟﺘﻲ ﺗﺘﻄﻠﺐ أﯾﻀﺎ اﺳﺘﺨﺪام‬.‫ ﻣﺜﻞ اﻟﮭﻮاﺗﻒ اﻟﻤﺤﻤﻮﻟﺔ أو اﺳﺘﻘﺒﺎل إﺷﺎرات اﻟﺒﺚ ﻋﺒﺮ اﻷﻗﻤﺎر اﻟﺼﻨﺎﻋﯿﺔ‬،‫ﻟﻠﺘﻄﺒﯿﻘﺎت ﻣﻊ ﻧﻘﺎط‬
‫ھﻮاﺋﯿﺎت ﺻﻐﯿﺮة ھﻮ ان اﻻﺑﻦ ﺑﺪون اﻟﺘﺮاﺑﻂ‬
‫ وﯾﺴﺘﺨﺪم ھﺬا اﻟﻨﻮع ﻣﻦ‬.‫ واﻛﺘﺸﻔﻨﺎ ﻓﻲ ھﺬه اﻟﺪراﺳﺔ دراﺳﺔ ﻋﺎﻣﺔ ﻣﻦ اﻟﮭﻮاﺋﯿﺎت اﻟﻤﻄﺒﻮﻋﺔ‬،‫ﻣﻦ ﺑﯿﻦ ﺗﻘﻨﯿﺎت ﻣﺨﺘﻠﻔﺔ ﻣﻦ اﻟﺘﺼﻐﯿﺮ ﻣﻦ اﻟﮭﻮاﺋﯿﺎت‬
.‫ إﻣﺎ ﻓﻲ ﺷﻜﻞ ﺷﺒﻜﺎت ﺣﯿﺚ ﯾﺘﻢ ﺗﺠﻤﯿﻊ ﻋﺪة ھﻮاﺋﯿﺎت ﻋﻠﻰ اﻟﺮﻛﯿﺰة ﻧﻔﺴﮫ‬،‫ اﻣﺎ اﻟﻌﻨﺼﺮ ھﻮاﺋﻲ واﺣﺪ‬،‫اﻟﮭﻮاﺋﯿﺎت اﻟﻤﺴﻄﺤﺔ‬
‫ واﻟﺬي ھﻮ ﻣﻌﯿﺎر ھﺎم‬،‫ ﺷﺮﺣﻨﺎ ﻇﺎھﺮة اﻗﺘﺮان‬،‫ أﺧﯿﺮا‬.‫ ﻻﺧﺘﻼف أﺳﺎﻟﯿﺐ اﻟﺘﻐﺬﯾﺔ‬HFSS ‫ﻧﺤﻦ ﻋﻠﻰ ﻏﺮار أداء اﻟﮭﻮاﺋﻲ رﻗﻌﺔ ﻣﺴﺘﻄﯿﻠﺔ ﻣﻦ ﻗﺒﻞ‬
.‫اﻟﺬي ﯾﺠﺐ ان ﺗﻀﻄﻠﻊ ﻋﻨﺪ ﺣﺴﺎب ﺧﺼﺎﺋﺺ اﻟﺸﺒﻜﺔ‬

:‫اﻟﻜﻠﻤﺎت اﻻﻓﺘﺘﺎ ﺣﯿﺔ‬

HFSS ‫ ﻣﺤﺎﻛﺎة‬،‫ وﺗﺤﻠﯿﻞ ﻃﺮق وﻧﻤﻮذج اﻟﻤﺤﺎﻛﺎة‬،‫ واﻟﺘﺤﺴﯿﻦ‬،‫ واﻟﺘﻮﻟﯿﻒ‬،‫ وﻧﻤﻂ إﺷﻌﺎع‬،‫ وﻣﺠﻤﻮﻋﺔ اﻟﮭﻮاﺋﻲ‬،‫اﻟﻤﻄﺒﻮﻋﺔ اﻟﮭﻮاﺋﻲ‬
  Page de titre
 Remerciements
 Résumé
 Sommaire
 Introduction Générale………………………………………………………………...1
Chapitre I : Généralités sur les Antennes Imprimées

I.1-Introduction …………………………………………………………………………………………......2
I.2- Description des antennes à éléments imprimés …………………………………………..2
I.2.1-Plan de masse …………………………………………………………………………………………2
I.2.2-Substrat ……………………………………………………………………………………………......2
I.2.3-Éléments rayonnants ………………………………………………………………………………….3
I.2.4- Alimentation des éléments rayonnants ………………………………………………………………5
1-Alimentation par ligne micro ruban …………………………………………………………………...5
 Alimentation par ligne micro ruban à travers un bord rayonnant………………………………6
 Alimentation par ligne micro ruban à travers un bord non rayonnant……………………….....6
 Antenne patch rectangulaire couplée capacitivement à une ligne micro ruban…………………6
2-Alimentation par une sonde coaxiale …………………………………………………………………..6
3-Alimentation couplée par fente (ouverture) …………………………………………………………...7
4-Alimentation couplée par proximité……………………………………………………………………7
 Comparaison entre ces différentes techniques d’alimentations………………………………..8
5-Alimentation par guide d’ondes coplanaires………………………………………………………….....8
6-Alimentation par fibres optiques …………………………………………………………......................9
I.3-Caractéristiques des antennes imprimées …………………………………………………9
I.3.1- Représentation en quadripôles………………………………………………………….......................9
I.3.2- Coefficient de réflexion …………………………………………………………................................9
I.3.3- Impédance d’entrée de l’antenne …………………………………………………………................10
I.3.4-Fréquence de résonnance f r …………………………………………………………........................10
I.3.5-Bande passante ………………………………………………………….............................................10
I.3.6-Taux d’onde stationnaire …………………………………………………………..............................11
I.3.7- Directivité et gain…………………………………………………………........................................11
I.3.8- Angle d’ouverture a - 3 dB ………………………………………………………….........................12
I.4- Alimentation des éléments rayonnants pour polarisations orthogonales ou Circulaires.12
I.4.1-Polarisations orthogonales…………………………………………………………...........................12
I.4.2-Polarisations circulaires…………………………………………………………...............................13
I.5-Mécanisme de rayonnement…………………………………………………………......................13
I.6-Rayonnements d’un élément rectangulaire ……………………………………………………...14
I.6.1-Principe du rayonnement de cet élément ……………………………………………………............14
I.6.2-Calcul du champ rayonné par une fente ……………………………………………………..............15
I.6.3-Calcul du champ rayonné par l’élément rectangulaire ……………………………………………….17
I.7-Choix de la technologie des antennes patchs…………………………………………………….18
I.8-Applications des antennes imprimées……………………………………………………………..18
I.9-Conclusion ……………………………………………………………………………………………19

Bibliographie ……………………………………………………………………………………………..20
Chapitre II : Méthodes d’Analyse des antennes imprimées

II.1-Introduction…………………………………………………………………………………………...21
II.2- Différentes Méthodes d'analyse ………………………………………………………………….21
 II.2.1- Méthodes analytiques ……………………………………………………………………………...21
II.2.1.1-Modèle de la ligne de transmission…………………………………………………………….22
1-Introduction………………………………………………………………………………………...22
2- Modèle et caractéristiques………………………………………………………………………...22
3- calcul de l’impédance d’entrée……………………………………………………………………25
4- Calcul de la fréquence de résonance………………………………………………………………26
II.2.1.2-Modèle de cavité………………………………………………………………………………..26
 Calcul du champ rayonné ……………………………………………………………………...28
Rayonnement des courants électriques (  r  1 ) ……………………………………….28
a) Diagramme de rayonnement dans le plan E (   0 ) …………………………..28
b) Diagramme de rayonnement dans le plan H      …………………………29
 
 2 
Rayonnement par application du théorème d’équivalence …………………………….30
a) Le diagramme de rayonnement dans le plan E ………………………………...31
b) Le diagramme de rayonnement dans le plan H………………………………...32
II.2.1.3-Comparaison des deux modèles analytiques…………………………......................................33
II.2.2- Méthodes numériques exactes (full-wave) …………………………...............................................34
II.2.2.1--méthode des moments …………………………........................................................................34
 Méthode de Green…………………………...............................................................................36
II.2.2.2- Matrices de ligne de transmission TML…………………………..............................................37
II.2.2.3-La méthode des différences finies…………………………........................................................38
1-Principe de la méthode …………………………..............................................................................38
2-Point de départ : les équations de Maxwell…………………………................................................38
3-Principe des différences finies centrées………………………….....................................................38
4-La discrétisation des équations de Maxwell …………………………..............................................39
4.1-La discrétisation spatiale…………………………...................................................................39
4.2- Construction du maillage…………………………..................................................................41
4.3- La discrétisation temporelle………………………….............................................................41
4.4- Détermination du critère de stabilité…………………………...............................................42
4.5- Dispersion numérique …………………………......................................................................42
II.2.2.4-La méthode des éléments finis …………………………...........................................................43
 Théorie de la FEM …………………………...........................................................................43
II.3-Conclusion …………………………...................................................................................................45
Bibliographie …………………………..................................................................................................46
Chapitre III : Modélisation et simulation (par HFSS) d’une antenne patch en
utilisant les différentes techniques d’alimentation

III.1-Introduction…………………………...............................................................................................48
III.2- Description du logiciel HFSS ………………………….............................................................48
III.3-procédure de conception de l’antenne rectangulaire…………................................................49
III.4-Techniques d’alimentation d’une antenne imprimée rectangulaire en utilisant HFSS .50
III.4.1-Alimentation par ligne micro ruban …………................................................................................50
III.4.1.1-Présentation de l’antenne étudiée …………............................................................................50
III.4.1.2- Influence des paramètres géométriques sur l’antenne………….............................................56
III.[Link] -l’influence de la permittivité relative du substrat ..........................................................56
III.[Link] - l’influence de hauteur du substrat................................................................................57

III.[Link]- l’influence de largeur du patch.......................................................................................58

III.[Link]- l’influence de longueur du patch.................................................................................59
III.[Link]- l’influence de l’épaisseur du patch ...............................................................................59
III.[Link]-effet de conductivité du patch......................................................................................60
III.4.2-Alimentation par ligne micro ruban avec un transformateur d’impédance...................................64
III.4.3-Alimentation par ligne enterrée ....................................................................................................91
III.4.4-l’alimentation par câble coaxiale ................................................................................................. .67
III.5-Antennes pour application radiocommunication…..................................................................70
III.6-Conclusion…......................................................................................................................................72
Bibliographie ….....................................................................................................................................73
Chapitre IV : Modélisation et simulation d’un réseau d’antennes patch
en utilisant HFSS
IV.1- Introduction…...................................................................................................................................74
IV.2-Définition du couplage dans le plan E et le plan H…..............................................................................74
 Couplage vertical ou couplage dans le plan H..........................................................................74
 Couplage horizontal ou couplage dans le plan E......................................................................75
IV.3-Contribution à la modélisation du couplage..............................................................................75
IV.3.1-L’équation matricielle....................................................................................................................75
IV.3.2-Modélisation par la méthode analytique .........................................................................................75
IV. 3. 3-Modélisation par la méthode numérique.....................................................................................77
 Convergence des résultats numériques..................................................................................79
 Ports d’impédance.................................................................................................................79
IV.4-Les réseaux d’antennes................................................................................................................80
IV.4.1-Alimentation des réseaux d’antennes imprimées........................................................................82
1-Réseaux linéaires à alimentation série...........................................................................................83
2-Réseaux linéaires à alimentation parallèle....................................................................................83
3-Réseaux bidimensionnels..............................................................................................................84
IV.5- Modélisation d’un réseau de deux patchs par le HFSS...................................................85
IV.5.1-Réseau linéaire de deux patchs à alimentation série..................................................................85
 Couplage dans le plan H..........................................................................................................85
 Couplage dans le plan E..........................................................................................................87
IV.5.2-Réseau linéaire de deux patchs à alimentation parallèle..........................................................89
IV.6- Modélisation d’un réseau à quatre patchs alimenté en parallèle .......................................91
IV.7- Modélisation d’un réseau bidimensionnel par le HFSS................................................94
IV.7.1- Réseau bidimensionnel à quatre patchs .................................................................................94
IV.7.2- Réseau bidimensionnel à 8 patchs......................................................................................95
IV.8-Conclusion..............................................................................................................................97
Bibliographie .................................................................................................................................98
Chapitre V : Réalisations

V.1- Introduction...................................................................................................................................99
V.2- Les étape de réalisation..............................................................................................................99
V.3- exemples de réalisation et simulation......................................................................................99
Premier exemple : Patch alimenté par un transformateur d’impédance..........................................99
 Plan horizontal……………………………………………………………………………...100
 Plan verticale……………………………………………………………………………......101
Deuxiéme exemple : Patch alimenté par ligne microruban directe………………………………..102
 Plan horizontal……………………………………………………………………………..103
 Plan verticale……………………………………………………………………………….104
Troisiéme exemple : Réseau à quatre patchs alimentée en paralléle……………………………..105
V.5-conclusion…………………………………………………………………………………………..107
 Conclusion Générale…………………………………………………………………....108
 Annexes…………………………………………………………………………109
 Liste des figures
 Liste des tableaux
 Introduction Générale

Introduction Générale
Dans un système de radiocommunication, l’antenne est un composant à part entière qui nécessite une
étude particulière. Tout en cherchant à optimiser les performances radioélectriques d’une antenne, on doit
l’adapter aux applications les plus récentes (communications : spatiales, mobiles, sans fil ou encore antennes
pour automobiles) et répondre aux exigences d’intégration dans l’architecture des terminaux. L’antenne doit
améliorer les performances radio globales d’un terminal. Il est donc nécessaire de concevoir des antennes
bien adaptées, avec un bon gain, une bonne pureté de polarisation, linéaire ou circulaire sur des bandes des
fréquences souvent très larges pour les antennes planaires, omnidirectionnelles ou au contraire très
directives, pilotables en fréquence ou en polarisation etc…. Mais on doit chercher aussi des antennes à faible
coût et les plus petites possibles (on dispose de peu de place sur un terminal mobile) au plus près desquelles,
on doit avoir la possibilité d’intégrer facilement des composants devant assurer d’autres fonctions
électroniques, et à partir de la, les chercheurs ont commencé à orienter leurs travaux vers la conception
d’antennes pouvant satisfaire les nécessités et contraintes précédemment cités, ces antennes sont
communément appelées antennes « patch » ou « antennes imprimées ».

L’antenne patch est conçue pour satisfaire ces besoins, c’est un conducteur métallique de forme
particulière placé sur un substrat terminé par un plan de masse, son caractère en miniature offre la possibilité
de l’intégrer facilement dans les systèmes d’émission réception.

L’objectif de notre travail est de faire une étude des antennes imprimées qui ont des applications aux
systèmes radiocommunications. De ce fait, notre mémoire est réparti en cinq chapitres :

 Le premier chapitre de ce manuscrit est consacré à une étude générale des antennes imprimées,
ainsi que leurs caractéristiques physiques, électriques et de rayonnement. Nous avons présenté
ensuite les différentes techniques d’alimentation et le mécanisme de rayonnement d’un élément
rectangulaire.

 Le deuxième chapitre énumère quelques méthodes d’analyses: des méthodes analytiques simples
(modèle de la ligne de transmission et le modèle de la cavité), et d’autre méthodes numériques
(méthode des moments, méthode des différences finies, méthode des éléments finis Matrices de ligne
de transmission TML).

 Le troisième chapitre est consacré à la modélisation et la simulation d’une antenne imprimée de

forme rectangulaire en utilisant les différentes techniques d’alimentation (ligne directe, ligne avec
transformateur, ligne enterré et câble coaxial). Les résultats de simulations obtenus sont effectués par
le logiciel HFSS (High Frequency Structure Simulator).

 Le quatrième chapitre est réservé à l’étude de l’effet du couplage sur les performances de
l’antenne (tel que : le gain, bande passante et coefficient de réflexion) d’un réseau linéaire pour les
deux type d’alimentation série et parallèle et d’un réseau bidimensionnel.

 Enfin, Le cinquième chapitre est dédié à la réalisation et la caractérisation de deux antennes

imprimées l’une alimentée par ligne directe et l’autre à travers un transformateur d’impédance. En
plus nous avons réalisé un modèle d’un réseau d’antennes linéaire à quatre patchs alimentés en
parallèle, la réalisation est faite au niveau du laboratoire de Télécommunications au sein du
département Electronique et Télécommunication

1
 Chapitre I Généralités sur les Antennes Imprimées
I.1-Introduction :
Les antennes sont devenues les éléments de base pour une liaison radioélectrique. Les technologies
modernes s’orientent vers la miniaturisation de ces antennes tout en essayant de garder les meilleures
performances.
Dans ces dernières années, divers types d’antennes imprimées ont été développés, et différentes
techniques d’alimentation sont mises en œuvre pour améliorer les performances de ce type d’antennes dans
plusieurs cas d’applications. La technique des circuits imprimés a révolutionné les domaines de
l’électronique, et plus récemment celui des hyperfréquences, où elle se prête particulièrement bien à la
réalisation de circuits et d’antennes.
L’antenne imprimée (patch) est conçue pour satisfaire à ces besoins, son caractère en miniature offre la
possibilité de l’intégrer facilement dans les systèmes d’émission - réception.
Dans ce chapitre nous présentons la description et le mécanisme de fonctionnement, les avantages, les
inconvénients des antennes patchs, ainsi que les différents types d'alimentation.

I.2- Description des antennes à éléments imprimés :

Une antenne à éléments rayonnants imprimés est constituée d’un plan de masse et d’un substrat
diélectrique, dont la surface porte un ou plusieurs éléments métallisés (figureI.1), appelés en anglais «
patch » (« pavé » ou « motif » en français) [I.1] [I.2] [I.3].

figure I.1 : Schéma de principe d’une antenne à éléments rayonnants imprimés

Ces éléments rayonnants sont de différentes formes et leur alimentation peut être effectuée par divers
procédés permettant d’obtenir un diagramme en polarisation linéaire ou circulaire [I.2] [I.4].
Les paramètres physiques et géométriques liés à cette structure sont :
 La permittivité relative de diélectrique (  r ).
 La tangente des pertes ( tag ) dans ce même substrat, avec dominance des pertes par effet joule.
 L’épaisseur du diélectrique (elle doit rester faible par rapport à la longueur d’onde à transmettre).
 Les dimensions de l’élément rayonnant (L, W).

I.2.1-Plan de masse : [I.6]

Doit être théoriquement de dimension infinie mais en pratique celui-ci est plutôt, pour des raisons
d’encombrement, de l’ordre de 3 ou 4 longueur d’onde. Il arrive cependant que celui-ci soit très réduit, ce
qui bien évidement modifie les caractéristiques de l’antenne initiale, parfois de maniéré très prononcée.

I.2.2-Substrat : [I.6] [I.7]

Les caractéristiques du substrat conditionnent également le comportement de l’antenne. On choisit en
général un diélectrique de faible hauteur vis-à-vis de la longueur d’onde  d’utilisation (quelques  /10 ),
idéalement la permittivité relative du substrat doit être faible (  r  3 ) de façon à permettre le rayonnement,
en évitant le confinement des champs dans la cavité comprise entre l’élément imprimé et le plan de masse.

2
Chapitre I Généralités sur les Antennes Imprimées
Les matériaux les plus couramment utilisés sont des composites à base de téflon
( 2   r  3 et tg   10 ), du polypropylène (  r  2,18 et tg   3.10 4 ) ainsi que des mousses synthétiques
3

contenant beaucoup de minuscules poches d’air (  r  1, 03 et tg   103 )

Le substrat est caractérisé par sa permittivité relative  r , sa hauteur h et son facteur de perte tan  , il
doit répondre à certaines exigences :
 Une conservation de la forme originale
 La stabilité en température (pour soudure) et aux conditions atmosphériques
 L’uniformité de la permittivité relative et de hauteur (sur de grande plaque)
 faible pertes diélectrique tan   103 pour une meilleure efficacité
 Une résistance aux produits chimiques
 Une faible rugosité
Substrat mince avec des constantes diélectrique élevées sont souhaitables pour les circuits micro-ondes,
car ils nécessitent étroitement liée champ pour minimiser le rayonnement indésirable et le couplage

I.2.3-Éléments rayonnants :
Les métallisations sont en très bons conducteurs tels que le cuivre, l’argent ou l’or (de conductivités
respectives   5,8.107 ; 6,2.107 et 4,1.107 S / m ).
Afin de simplifier l'analyse et l'estimation des performances, le patch a généralement une forme
rectangulaire, circulaire, triangulaire, elliptique ou autre forme connue comme indiqué dans la figure (I.2)
Les divers types d’éléments rayonnants les plus simples sont [I.1] [I.4] [I.7]:

figure I.2 : Divers types d’éléments rayonnants.

Il existe une grande variété de forme que peut prendre une antenne plaquée ou imprimée. Il y a tout
d’abord le dipôle imprimé, qui attire par sa petitesse quelque application. Cependant son exploitation reste
limitée à cause de sa bande passante très étroite. Les plus souvent rencontre, sons les antennes rectangulaires
ou carres, circulaires et elliptiques. Elles sont simples a réaliser, et ont été étudiés et modélisées par plusieurs
auteurs [I.7] [I.8].
On trouve également d’autres formes plus complexes et plus difficiles à analyser. Elle résulte souvent
de la combinaison de deux formes simples et sont utilisées dans certaines applications particulières [I.7].
Le choix de la forme d’antenne dépendra du type d’application désiré et des paramètres qu’on cherche à
optimiser (bande passante, gain et efficacité, impédance d’entrée, niveaux des lobs secondaires, …etc.)[I.7].
Parmi toutes les formes l’élément rectangulaire est sans contexte le plus facile à appréhender pour la
compréhension des mécanismes de rayonnements des antennes patch. Mais tout d’abord se pose la question
d’alimenter une telle structure [I.7].

3
Chapitre I Généralités sur les Antennes Imprimées

figure I.3 : Présentation d'une antenne patch a élément rectangulaire

Pour un patch rectangulaire, sa longueur L est généralement 0.33330  L  0.50 , ou 0 est la
longueur d'onde dans l'espace libre. Le patch est conçue de façon qu'il soit très mince ( t  0 , ou t est
l'épaisseur du patch) La taille h du substrat diélectrique est habituellement 0.0030  h  0.050 [I.7] [I.10]
[I.8] [I.9].

Pour augmenter la bande passante, on peut aussi utiliser des éléments rayonnants étagés (figure I.4) :
[I.2]

figure I.4 : éléments rayonnants étagés

 Un premier élément E1 sur un substrat de permittivité  1 , de fréquence centrale f 1 [I.2] [I.10].
 Un second élément E 2 sur un substrat de permittivité  2 , de fréquence centrale f 2 , au-dessus du plan
de masse de la structure [I.2].
Si les fréquences f 1 et f 2 sont relativement voisines. Le couplage électromagnétique entre les deux
éléments superposés permet un élargissement de la courbe de réponse fréquentielle de l’ensemble [I.2] [I.5].
C’est L’élément rayonnant le plus utilisé est celui de forme rectangulaire, dont l’axe de symétrie passe
par le point d’excitation [I.2].

Ses dimensions sont la largeur w qui doit être inférieure à d  0  r pour éviter les modes d’ordre
supérieur et sa longueur L , très légèrement inférieur à la demi-longueur d’onde dans le diélectrique de
permittivité équivalente  e [I.2] [I.5] [I.9].

L  0,49e  0,49 0
e
(I .1)
 e se calcule conformément aux formules situées au-dessous [I.11] :
1 1  h
1
 w 
2

 e    r  1    r  1  1  12   0.04  1    pour : w / h  1 (I .2)

2 2  w  h  

1
1 1  h
 e    r  1    r  1  1  12  pour : w / h  1 (I .3)
2 2  w

4
Chapitre I Généralités sur les Antennes Imprimées
Vu la faible hauteur h des substrats par rapport à la largeur w de la métallisation, on se trouve dans le
cas ou w h  1 , si bien que  e est voisin de  r pour les valeurs de  r utilisées [I.2] [I.4].
Il existe aussi des éléments rectangulaires dont la longueur L  e 4 , l’extrémité opposée à
l’alimentation est court-circuitée par des petits barreaux relient le bord de l’élément au plan de masse. On
utilise également des éléments de longueur L  e , avec une alimentation centrale [I.2].

I.2.4- Alimentation des éléments rayonnants :

L’alimentation joue un rôle très important lors de la conception des antennes imprimées. Elle peut
modifier les caractéristiques de rayonnement de celles-ci. L’alimentation pourra être soit centrée soit
arbitraire par rapport à la géométrie de l’élément. Par sa position, elle détermine l’adaptation de l’antenne à
la ligne d’alimentation. Il existe une large variété de techniques d’alimentation permettant ainsi de fournir à
l’antenne, l’énergie qui lui est nécessaire pour rayonner [I.7].

1-Alimentation par ligne micro ruban :

Dans cette méthode d’alimentation, une ligne microruban (généralement de 50 ) est connecte
directement au bord du patch rayonnant. La longueur de la bande conductrice est plus petite par rapport au
patch. Ce genre d’alimentation a l’avantage que l’alimentation peut être gravée sur le même substrat pour
fournir une structure planaire .C’est la technique la moins coûteuse et la plus simple a réalisé [I.7] [I.12]
[I.13].

figure I.5 : Alimentation par ligne Micro ruban

L’alimentation peut se faire par connexion directe à une ligne microbande dont le point de jonction est
sur l’axe de symétrie de l’élément (figure 6.a) ou décalé par rapport à cet axe de symétrie (figure6.b) si cela
permet une meilleure adaptation d’impédance, dans ce but, l’alimentation axiale avec une encoche donne de
bons résultats (figure 6.c) [I.2] [I.4][I.14].

figure I.6: Alimentation par une ligne microbande : a. axiale, b. décalée, c. axiale avec encoche.
L’avantage de ce type d'alimentation est sûrement la moins coûteux et la plus aisée est qu’elle peut être
gravée sur la même face de l’antenne, elle est facile à fabriquer, et simple à adapter à la résonance. C’est
l’une des techniques les plus utilisées dans les antennes imprimées, Les principaux inconvénients de cette
structure simple sont [I.15]:
- les couplages parasites dû la réalisation de la ligne de transmission et du résonateur sur le même substrat,
et qui diminuent les performances de l’antenne au-delà de quelques GHz
- la nécessité d’ajouter un transformateur d’impédance pour obtenir l’adaptation entre la ligne
d’alimentation et le résonateur.
On peut distinguer trois types d’alimentation [I.3] :
5
Chapitre I Généralités sur les Antennes Imprimées
 Alimentation par ligne micro ruban à travers un bord rayonnant:
Dans ce cas la ligne d'alimentation entre en connexion directement sur le côté rayonnant du patch
figure I.7.a, cette technique est l'une des méthodes les plus utilisées. Il existe une variation de l'impédance
d'entrée car le point d'alimentation est déplacé le long de ce côté.

figure I.7. a : Le patch alimenté par ligne micro ruban à travers un bord rayonnant

 Alimentation par ligne micro ruban à travers un bord non rayonnant:

La ligne d'alimentation entre en contact avec le patch sur le bord non rayonnant figure I.7.b, dans ce cas
les calculs sont plus complexes, car les courants de la ligne d'alimentation, sont orthogonaux à ceux du patch

figure I.7.b : Le patch alimenté par ligne micro ruban à travers un bord non rayonnant

 Antenne patch rectangulaire couplée capacitivement à une ligne micro ruban:

Dans ce cas l’élément rayonnant est couplé capacitivement à la ligne micro ruban injectée dans le
substrat figure I.7.c. L’utilisation de substrats différents entre ligne d’excitation et antenne peut améliorer la
largeur de la bande passante de l’antenne.

figure I.7.c : Le patch couplé capacitivement à ligne micro ruban

2-Alimentation par une sonde coaxiale :

Cette méthode d’alimentation a beaucoup été étudiée dans la littérature. Elle à l’avantage de positionner
l’alimentation sous l’élément rayonnant et à utiliser une ligne non rayonnant (coaxial), ce qui diminue
fortement de rayonnement parasite [I.1] [I.14], et qu’elle peut être appliquée à n'importe quel endroit choisi
à l’intérieur du patch, avec une facilité de fabrication [I.12].
Le conducteur central du câble coaxial est soudé sur l’élément rayonnant après avoir traversé le plan de
masse et le diélectrique. Le conducteur extérieur est quant à lui soudé au plan de masse [I.13] [I.15].
Cependant, cette méthode présente des inconvénients au niveau du diagramme de rayonnement. En
effet, la connexion génère un pic de courant localisé au niveau de l’élément rayonnant qui peut induire une
dissymétrie dans le diagramme de rayonnement. De plus, des pertes apparaissent avec le perçage du plan de
masse, du diélectrique ainsi que de l’élément plaqué [I.1] [I.9] [I.12].

6
Chapitre I Généralités sur les Antennes Imprimées

figure I.8 : Alimentation par connexion directe à une ligne coaxiale

3-Alimentation couplée par fente (ouverture) :

Dans ce type d'alimentation, le patch de rayonnement et la ligne d'alimentation du micro ruban sont
séparés en plan de masse comme indique sur la figure [Link] conjonction entre le patch et la ligne
d'alimentation est faite par une ouverture ou une fente dans le plan de masse .Généralement, les matériaux à
permittivité élevée sont employés pour le substrat inférieur par contre les matériaux à faible constante
diélectrique sont utilisés pour le substrat supérieur dans le but d’optimiser le rayonnement du patch
[I.12][I.13] [I.15][I.14].
Ce type d'alimentation est difficile à concevoir à cause des couches multiples, qui augmentent
l'épaisseur d'antenne. Cependant elle offre l’élargissement de bande passante [I.12].
L’inconvénient de telle technique et sa complexité de fabrication. Cette complexité réside dans : l’ouverture
dans le plans de masse, La réalisation de deux couche diélectriques [I.7] [I.10].
Elle reste tout de même utilisée pour certaines applications, car elle permet une amélioration de la bande
passante qui peut aller jusqu'à 25% [I.7]

figure I.9 : Alimentation par couplage électromagnétique à une ligne microbande à travers une fente
dans le plan de masse.

4-Alimentation couplée par proximité :

Cette technique d’alimentation est connue également sous le nom de couplage électromagnétique. Deux
substrats diélectriques sont employés tels que la ligne d’alimentation est située entre les deux substrats et le
patch est imprimé sur le substrat supérieur (figure I.10) [I.7] [I.12] [I.14].
L’avantage principal de cette technique est qu’elle élimine le rayonnement parasite due à l’alimentation
et fournit une largeur de bande très élevée (plus que 13%), en raison de l’augmentation globale de
l’épaisseur de l’antenne microbande. Cette technique fournit également des choix entre deux milieux
diélectriques différents, un pour le patch et un pour la ligne d’alimentation pour optimiser les performances
de l’antenne [I.7] [I.12].
Parmi les inconvénients de cette méthode d'alimentation nous citons la difficulté de fabrication à cause
des deux couches diélectriques nécessitant un alignement approprié, ainsi que la difficulté d’intégration de
dispositifs actifs [I.12].

7
Chapitre I Généralités sur les Antennes Imprimées

figure I.10 : Alimentation par couplage électromagnétique à une ligne microbande

 Comparaison entre ces différentes techniques d’alimentations : [I.9] [I.16]

Tableau I.1: Comparaison entre les différentes techniques d’alimentations

5-Alimentation par guide d’ondes coplanaires :

Le plan de masse et la ligne d’alimentation se trouvent sur la même face du substrat. A cette partie est
rajoutée un élément rayonnant imprimé sur un second substrat déposé au-dessus à hauteur des fentes
rayonnantes. Le pavé rayonnant est alors excité par couplage électromagnétique [I.1].

figure I.11 : Alimentation par guide d’ondes coplanaires

8
Chapitre I Généralités sur les Antennes Imprimées
6-Alimentation par fibres optiques :
Ce type de liaison présente plusieurs avantages, comme une large bande passante, une forte liaison, une
grande immunité au bruit, une taille et un poids réduits. Le concept est le suivant : une diode laser est
modulée à l’aide d’une fibre optique. Sur le substrat de l’antenne une photodiode démodule le signal
lumineux qui sera ensuite couplé à l’antenne imprimée par un circuit d’adaptation [I.7].

figure I.12: Alimentation par fibre optique

I.3-Caractéristiques des antennes imprimées :

Dans les communications sans fil, chaque application met en relief certaines caractéristiques des
antennes. D'une manière générale, une antenne utilisée dans un type d’application ne peut pas l'être dans
d’autres.
Une antenne peut être caractérisée par :

I.3.1- Représentation en quadripôles : [I.7]

Une antenne peut être représentée comme un quadripôle (Figure I.4), défini par les paramètres S :

figure I.13 : Représentation en quadripôle

A1 , A2 , B1 , B2 Sont des ondes de puissance.
On a les relations suivantes :
B1  S11. A1  S12 . A2
(I.4)
B2  S 21. A1  S22 . A2
 Le coefficient S11 correspond à la réflexion en entrée des quadripôles lorsque A2  0 .
 Le coefficient S12 représente la transmission de la puissance entrant en sortie vers l’entrée
lorsque A1  0 .
 Le coefficient S 21 est le gain du quadripôle lorsque A2  0 .
 Le coefficient S 22 est la réflexion en sortie du quadripôle lorsque A1  0 .
I.3.2- Coefficient de réflexion : [I.6] [I.7]
La représentation du coefficient de réflexion S11 (ou  ) en fonction de la fréquence permet de
déterminer la fréquence de résonance de l’antenne pour des dimensions données.

9
Chapitre I Généralités sur les Antennes Imprimées
Z e  Z c onde réflechie
  (I.5)
Z e  Z c onde incidente
Ze  Zc
S11 (dB)  10 log (I.6)
Ze  Zc
Z e : Impédance d’entrée de l’antenne
Z c : Impédance caractéristique de la ligne qui alimenté l’antenne ( Z c  50 ).
Une antenne est un circuit résonnant et l’impédance est ramenée peut être égale à l’impédance
caractéristique à cette fréquence de résonance.
Une antenne résonne si la partie imaginaire de son impédance d'entrée est nulle (Pratiquement on prend
le minimum de la partie imaginaire),
  0  Re Z e   Z c (I.7)

Et dans laquelle la directivité de l'antenne est maximal, et le gain est plus élevé, c'est à dire le
rayonnement est plus directif, et c'est la définition pour déterminer la fréquence de résonance.

I.3.3- Impédance d’entrée de l’antenne : [I.6] [I.16]

En ce qui concerne l’émission et conformément à la théorie des circuits linéaires, l’antenne peut être
représentée par un dipôle d’impédance d’entrée complexe Z e  f   Re  f   jX e  f  « chargeant » la sortie
du circuit émetteur, modélisé par une source d’impédance interne Z s  f   Rs  f   jX s  f  . Précisons que
f est la fréquence de résonance de l’antenne.

figure I.14 : schéma équivalent du dispositif à l’émission

D’une part, la partie réelle de l’impédance d’entrée Re  f  caractérise la puissance dissipée par
l’antenne sous forme de rayonnement et de pertes diverses dans sa structure (effet Joule, onde de surface,
etc.) et d’autre part, la partie imaginaire X e  f  , d’interprétation plus délicate est liée à la puissance réactive
concentrée au voisinage immédiat de l’antenne et au déphasage.
I.3.4-Fréquence de résonnance f r : [I.6]
C’est la fréquence ou l’antenne présente un minimum d’amplitude du coefficient de réflexion. Elle est
donnée par cette relation :
c 1
fr  (I.8)
2  r L  2L
Elle est en fonction de permittivité et les dimensions de l’antenne patch.

I.3.5-Bande passante : [I.6] [I.7]

La bande passante en pourcent est donnée par la relation suivante :
f  fi
B  %   100  s (I.9)
fc
f s : la fréquence limite supérieur
f i : La fréquence limite inferieur
10
Chapitre I Généralités sur les Antennes Imprimées
f c : La fréquence centrale d’utilisation pour laquelle l’antenne est conçue.

La bande passante donnée pour une certaine valeur de rapport d’onde stationnaire (ROS), S à l’entrée de
f S 1
l’antenne est donnée par :  (I.10)
fr Q S

I.3.6-Taux d’onde stationnaire : [I.6] [I.7]

Il est donné par la relation suivante en fonction du coefficient de réflexion :
1 
TOS  (I.11)
1 
Le coefficient de réflexion en fonction du TOS :
TOS  1

TOS  1
(I.12)
I.3.7- Directivité et gain : [I.2]
La directivité maximale est le produit par 4 du rapport puissance PM rayonnée par unité d’angle
solide dans la direction de rayonnement maximal, à la puissance totale rayonnée Pr .
La direction de rayonnement maximal est    2 ,    2 ; la puissance PM est :
1  1 1
PM   E  H   R 2  E2  R 2
2 2 120 (I.13)
V w
Avec: E  2 0
R 0 (I.14)
2
1 w
PM  V02  
60  0  (I.15)
1
La puissance totale rayonnée par une fente étant V02 G , la puissance totale rayonnée par l’élément
2
2
imprimé est : V0 G où G est la conductance d’une fente.
La directivité de l’élément rectangulaire est donc :
2
P 1 w
D  4 M   
Pr 15G  0  (I.16)
Nous avons vu que G  G f  Gm , si par une approximation au premier ordre, on néglige Gm devant G f ,
nous avons :

2
1 w
Df   
15G f  0  (I.17)
2
1 1200 w w
 Si w 0  1, Df     8 
15 w  0   0  (I.18)
2 2
1  0 
w
 Si w 0  1, Df     6
90 
 0 
15  w  (I.19)
En tenant compte de Gm , la directivité D est donnée par :
2
1 w
D  
15G f  Gm   0  (I.20)

11
Chapitre I Généralités sur les Antennes Imprimées
Df
Ou : D
1  g12 (I.21)
Gm
g12  est appelée la conductance mutuelle normalisée.
Gf
Le gain est déduit de la directivité connaissant le rendement  de l’antenne G0  D ; Le gain
augmente avec la hauteur h du substrat et si  r diminue.
I.3.8- Angle d’ouverture a - 3 dB : [I.2]
Les angles d’ouverture à- 3 dB, notés 2 3 , sont calculés à partir des formules donnant les champs
rayonnés dans les plans H et E . Moyennant quelques approximations et après calculs, Bahlet Bhartia
Ont obtenu les formules :
1 2
  w  
2 3 H  2 cos  21 
1
 

   0  (I.22)
1 2

2 3 E 1
 
 2 cos 1  3k 02 L2  k 02 h 2  
7  (I.23)
I.4- Alimentation des éléments rayonnants pour polarisations orthogonales ou
Circulaires :
I.4.1-Polarisations orthogonales: [I.2]

figure I.15: Orientation du plan E et des courants surfaciques sur l’élément rayonnant
a. selon la position .b. du point d’excitation

Les courants surfaciques qui sont produits sur l’élément métallique sont dirigés selon l’axe de symétrie
passant par le point d’excitation (figure I.15.a et b) : ils sont donc à polarisation uniforme et il en résulte que
les champs E et H rayonnés pour une excitation donnée sont à polarisations rectilignes.
Si les excitations sont effectuées sur deux axes de symétrie orthogonaux, les courants et les champs produits
sont à polarisations orthogonales, c’est-à-dire que le plan E lié à une excitation correspond au plan H de
l’autre, et réciproquement

figure I.16 : excitation permettant d’obtenir des polarisations orthogonales pour un élément : a-carré, b-
circulaire, c- à grilles orthogonales
Si ces excitations sont effectuées simultanément (figure I.16.a et 16.b) en des points homologues situés
sur de axes de symétrie d’un élément rayonnant carré ou circulaire, le rayonnement de l'antenne se fera
selon deux polarisations orthogonales.
Pour obtenir chacune des deux polarisations orthogonales avec une haute pureté de polarisation, de
chercheurs du CENT et de l’université de Nice ont proposé se constituer les éléments métalliques carrés par
deux grilles de barres orthogonales (figure I.16.c), chacune des grilles étant excitée de manière indépendante

12
Chapitre I Généralités sur les Antennes Imprimées
par des fentes à couplage électromagnétique orthogonales. Dans ces condition, le niveau de polarisation
croisée est inférieur à -30db sur une bande de 50℅ et de l’ordre de -40 dB sur une bande plus étroite
(quelques pour-cent).

I.4.2-Polarisations circulaires: [I.2] [I.8]

figure I.17 : création d’un déphasage de 90° entre les points d’alimentation :
pour des lignes dont les longueurs différent de  / 4 d’après Zurcher et Gardiol b- par un coupleur à 3dB

Pour avoir une polarisation circulaire, il faut que les deux excitations orthogonales dans l’espace soient,
de plus, déphasées de 90° dans le temps. Ceci s’obtient soit en utilisant deux lignes d’alimentation dont les
longueurs différent de  / 4 (figureI.17.a), soit en alimentant les deux points d’excitation à partir d’un
coupleur à 3dB (figureI.17.b) dont on sait qu’il fournit des ondes d’amplitudes identiques, mais déphasées
de 90˚.Pour un déphasage de + 90°la polarisation est circulaire droite (P.C.D), tandis qu’elle est circulaire
gauche (P.C.G) pour un déphasage de – 90°

figure I.18: Point d’alimentation pour obtenir une polarisation circulaire avec un élément
a. presque carré, b. carré avec une fente en diagonale, c. carré à coins tronqués.
Polarisation circulaire peut être aussi obtenue en alimentant un élément métallique rectangulaire,
presque carré (figure I.18.a) à partir d’un point situé sur une diagonale, ou encore avec des éléments carrés
alimentés en un seul point d’un de leurs axes se symétrie (figure I.18.b et I.18.c) si le carré est rendu
asymétrique, par exemple par une petite fente découpe selon une diagonale du carré ou en tronquant deux
coins opposés pour tous les cas représentés à la figure I.16, l’élément métallique entre en résonance du fait
de sa faible asymétrie, sur deux modes orthogonaux excités à partir du seul point d’alimentation .en
contrepartie de sa simplicité, cette méthodes ne fournit la polarisation circulaire que dans une bande de
fréquences très étroite et dans des conditions d’alimentation assez difficiles à maîtriser
I.5-Mécanisme de rayonnement: [I.16]
Pour mieux comprendre le mécanisme de rayonnement des antennes imprimées, la littérature
scientifique s’appuie généralement sur la théorie de l’optique.

figure I.19: Mécanisme de rayonnement

13
Chapitre I Généralités sur les Antennes Imprimées
L’élément rayonnant est en fait une multitude de sources ponctuelles similaires a celle qui se trouve au
point a. ces sources isotropiques vont rayonnées dans toute les directions.
Une partie de celui-ci est réfléchie par le plan de masse, puis par le conducteur supérieur et ainsi de
suite. Certains de ces rayons aboutissent sur l’arrête du conducteur (point « b ») qui les diffracte et par
conséquent le champ se décompose en deux composantes normale et tangentielle par rapport au plan de
masse. On peut donc envisager trois régions :

Région A :
C’est la plus dense en rayons. Dans cette région le champ électromagnétique s’accumule plus
fortement. Cette propriété est d’ailleurs utilisée pour la propagation le long des lignes micro rubans. Lorsque
les fréquences sont basses ; le champ électromagnétique est concentré dans cette zone. Donc il y a
propagation sans rayonnement.

Région B :
Région de l’espace libre au-dessus du substrat dans laquelle il y a dispersion libre du champ. La
contribution des sources ponctuelles est faible. La partie la plus grande du rayonnement provient des rayons
à l’extrémité du conducteur supérieur (point « b »), le reste du champ demeure à l’intérieur du substrat et
sera considéré comme perte dans le diélectrique. Nous obtenant ainsi une antenne dont le rendement en
rayonnement est très limitée.

Région C :
Certains rayons atteignent les surfaces de séparation (air-diélectrique) avec une incidence rasante et
restent piégés à l’intérieur du diélectrique. Il peut arriver que certain de ces rayons soit difracté en atteignant
le bord (Point « c ») du substrat, engendrant ainsi, un rayonnement parasite. Cette dernière est à l’ origine du
couplage inter éléments, lorsque les antennes imprimées sont mise en réseau. Au de la d’une certaine
fréquence, ce phénomène affecte considérablement le rayonnement de l’antenne. Ce si, détériore la qualité
du diagramme de rayonnement, particulièrement au niveau du lobe secondaire

I.6-Rayonnements d’un élément rectangulaire :

I.6.1-Principe du rayonnement de cet élément : [I.2] [I.4]

figure I.20 : Lignes de champ entre les bords d’un élément imprimé et le plan de masse

Un élément imprimé ABCD sur substrat de diélectrique et plan de masse rayonne comme les 4 fentes
qui le délimitent :
 Fentes AB et DC parallèles à Oz .
 Fentes AD et BC parallèles à Ox .
Les lignes de champ électrique qui relient l’élément imprimé au plan de masse sont indiquées sur
(figureI.20), pour un élément de longueur e 2 alimenté en S à l’une de ses extrémités. Il en résulte que les
fentes AB et DC , qui sont distantes de e 2 , présentent des lignes de champ en opposition de phase. Pour
les fentes AD et BC , il y a un renversement de phase au passage par x  0 , si bien que, d’une moitié à
l’autre de chacune de ces fentes, les lignes de champ sont en opposition de phase.

14
Chapitre I Généralités sur les Antennes Imprimées

figure I.21 : Lignes de champ entourant un élément imprimé : a. composantes horizontales, b.

composantes verticales.
Les (figures I.21.a et b) représentent les composantes horizontales et verticales de ces lignes de champ.
Si l’on raisonne par rapport aux plans principaux xOy (plan E) et yOz (plan H), on voit bien que la
contribution des fentes AD et BC est nulle puisqu’il y a toujours deux composantes situées en des points
symétries par rapport à ces plans, qui sont en opposition de phase. Il est de même pour la contribution des
composantes verticales des champs dans les fentes AB et DC .par contre, les composantes horizontales des
champs dans ces fentes créent des rayonnements qui sont en phase dans le plan yOz et qui interfèrent dans le
plan xOy .
Finalement, le rayonnement de l’élément imprimé étudié peut être modélisé par celui de deux fentes
parallèles à Oz , distantes de L, de dimensions w selon Oz et h (égale à l’épaisseur du substrat). Ces deux

fentes sont soumises à un champ E uniforme, polarisé parallèlement à l’axe Ox .

I.6.2-Calcul du champ rayonné par une fente : [I.2]

figure I.22 : Modélisation d’une fente excitée par un champ électrique à polarisation uniforme

La fente, excitée par un champ électrique à polarisation uniforme selon Ox , peut être modélisée par un
courant surfacique de type magnétique dirigé selon Oz .
m     
I s V / m   n  E   y  E x  E z (I.24)

Comme la largeur h de la fente est très faible devant  et que E est constant selon cette largeur, nous
pouvons donc assimiler le rayonnement de la fente à celui d’un courant magnétique linéique dont l’intensité
m m  
est : I V   I s h  Eh z  V0 z (I.25)
V0 Représente la tension entre les deux côtés de la fente.
L’équivalence que nous venons d’établir est importante car nous connaissons le champ rayonné par un
doublet électrique de longueur dz et nous pouvons en déduire, par le principe de dualité, le champ rayonné
par un doublet magnétique de longueur dz .à partir de là, nous saurons calculer le champ rayonné par le
m 
courant magnétique I  V0 z qui modélise la fente.

15
Chapitre I Généralités sur les Antennes Imprimées

   
 
figure I.23 : a-coordonnées et composantes  ,  pour les champs E  et H  rayonnés par un doublet
électrique, b- coordonnées et notations pour calculer le rayonnement d’un courant magnétique de
longueur W
Le champ électromagnétique rayonné, en zone lointaine, par un doublet électrique de hauteur dz ,
parcouru par un courant constant I 0  A à deux composantes.
120 1
E  j I 0 dz sin  e  jk0 R (I.26)
2 0 R
2
Avec : K 0  (I.27)
0
E
H   (I.28)
120
0
Ou : 120  (I.29)
0
Le champ électromagnétique rayonné par un doublet magnétique parcouru par un courant V0 (volts) a
deux composantes qui sont déduites des deux précédentes par le principe de dualité :
1 1 1
H  j V0 dz sin  e  jK0 R (I.30)
2 120 0 R
1
E  H .120  j V0 dz sin  e jK0 R (I.31)
20 R
m 
Calculons maintenant le champ électrique rayonné par le courant magnétique I  V0 z de longueur W
.prenons un point variable Q tel que QP  r :
W / 2
V0 e  jK0 r
20 W/ 2
E  j sin  Q dz (I.32)
r
Lorsque le point P est en zone lointaine de l’élément rayonnant, on peut considérer que QP // OP et
faire les approximations suivantes :
 Dans le terme de phase : r  R  z cos 
 Dans le terme d’amplitude : r  R
Dans ces conditions :
W / 2
V0
sin  e  e jK0 Z cos dz
 jK 0 R
E  j (I.33)
20 W / 2
Le calcul de l’intégrale est classique :
KW 
W /2 sin  0 cos  
2
I   e jK0 RZ cos dZ  2   (I.34)
W / 2
K 0 cos 
On en déduit :

16
Chapitre I Généralités sur les Antennes Imprimées
 W 
sin  cos  
V
E  j 0 sin   0  e jK0 R (I.35)
2 R cos 
En réalité, compte tenu de la présence du plan de masse, le champ rayonné doit être multiple par deux,
conformément au principe des images. Nous avons donc, en module :
 W 
sin  cos  
V
E  0 sin   0  (I.36)
R cos 

I.6.3-Calcul du champ rayonné par l’élément rectangulaire : [I.2]

figure I.24 : Modèle équivalent utilisé pour calculer le rayonnement d’un élément rectangulaire.
a).représentation dans l’espace, b).schéma dans le plan xOy .

Le rayonnement de cet élément est équivalent à celui des deux fentes représentées sur la figure I.24.a.
Pour l’étudier, nous allons prendre l’origine des axes dans une position symétrique par rapport aux deux
fentes.
 Dans le plan H (plan yOz correspondant à    2 ) :
Les champs rayonnés par les deux fentes s’additionnent en phase et le diagramme de rayonnement est
le même que celui d’une seule fente.
Le champ rayonné est donc :
 w 
sin  cos  
 V  0  x
E  2 0 sin 
R cos  (I.37)
 Dans le plan E (plan xOy correspondant à    2 ) :
le champ rayonné par chacune des fentes doit être calculé d’après la relation
w 
sin  cos  
V
E  0 sin   0  pour :  
R cos  2 (I.38)
sin x w
Pour faire ce calcul, nous devons mettre le dernier facteur sous la forme avec x  cos  afin
x 0
sin x
de pouvoir utiliser la propriété  1 quand x  0 .
x
Dans ces conditions :
  V w
E     0
 2  R 0 (I.39)

17
Chapitre I Généralités sur les Antennes Imprimées
Pour calculer le champ rayonné par les deux fentes (figure. I.24.b), nous devons tenir compte de ce
2 L
que la fente située à l’abscisse x   L 2 rayonne avec un retard de phase  cos  tandis que celle qui
0 2
2 L
est située à l’abscisse x  L 2 rayonne avec une avance de phase cos  .
0 2
Le champ total rayonné est donc :
 j   L cos  

0 R

V w  i   L cos
E  0 e  0   e  0  
(I.40)
 Vw  L 
D’où : E   2 0 cos cos  
0 R  0  (I.41)

 Dans le plan à φ =cte :

Passant par l’axe Oz , le champ rayonné dans une direction  ,   est donné, en module, par :
 w 
sin cos  
 V  0  cos L sin  cos  
E   2 0 sin   
R cos   0  (I.42)
I.7-Choix de la technologie des antennes patchs :
Actuellement, les antennes patchs deviennent de plus en plus importantes en raison de la grande variété
des formes géométriques qu'elles peuvent prendre et qui les rendent applicables aux différentes situations
d’intégration [I.3] [I.7].
Parmi les avantages de cette technologie d'antennes, on peut citer : leur volume réduit,
la conformabilité et la possibilité d'intégrer les circuits micro-ondes au niveau des antennes. Ajoutons que la
simplicité de leurs structures fait que ce type d’antenne est adapté au faible coût de fabrication. Ceci est une
propriété clé qui explique l’introduction des antennes patchs dans les applications de la communication
mobile [I.3] [I.7].
 D’une façon générale, les avantages des antennes imprimées sont : [I.15] [I.17] [I.18]
1 .Faible poids.
2 .Limitation du risque de panne.
3. Faible coût de fabrication, production en masse possible.
4 .Polarisation linéaire et circulaire.
5 .Antennes multi bandes, multi polarisations possibles.
6 .Simples, Robustes.
7 .non encombrants.
 Cependant, elle présente les inconvénients suivants : [I.15] [I.17] [I.18]
1. Bande passante étroite, souvent associée avec les problèmes de tolérances (géométriques et physiques).
2. Généralement faible gain (6 dB).
3. La plupart des antennes rayonnent uniquement dans le demi-plan.
4. Supportent uniquement des faibles puissances (quelques dizaines de watts).
5. Rayonnement parasite de l’alimentation.
6. Circuit d’alimentation parfois très complexe.
Donc, la conception des antennes doit répondre à des compromis en termes de performances et de
complexité de réalisation.

I.8-Applications des antennes imprimées : [I.1] [I.18]

Les demandes des antennes imprimées sont au loin très nombreuses, de l’électronique grand public aux
systèmes technologiques de pointe. Les antennes imprimées sont utilisées dans :
1. La télémétrie des missiles.
2. Les contrôles et les commandes à distance.
3. Les radars.
18
Chapitre I Généralités sur les Antennes Imprimées
4. Les Réseaux mobiles (GSM 900, DCS 1800, UMTS,…).
5. Réseaux sans fils : (WLAN, BLUETOOTH, WIFI, ...).
6 .Les communications spatiales (GPS).
Plusieurs systèmes complexes comme les radars de poursuite ont été miniaturisés à partir de la taille
réduite des antennes planaires.
Quelques domaines d’applications des antennes imprimées et leurs largeurs de bandes :

I.9-Conclusion :
La technologie microbande en générale, et les antennes patch en particulier connaissent un succès
croissant aux prés des industriels et des professionnels de télécommunication, qu’elle soit spatiale ou
terrestre. Il est nécessaire de bien définir l’intégralité des paramètres des antennes, afin de prévoir son
comportement avant même la réalisation, et de s’assurer qu’elle se conformera aux exigences des systèmes.
Nous avons essayé à travers ce chapitre de présenter de façon concise, l’essentiel de la théorie des
antennes imprimées dont aurait besoin tout ingénieur ou chercheur s’intéressant à ce type d’antenne.
On a cité quelques applications de ce type d’antennes, l’étude du fonctionnement d’une antenne
imprimée, on s’est intéressé par la suite aux plus importantes techniques utilisées pour son alimentation,
mécanisme de rayonnement ainsi qu’aux différentes caractéristiques de son rayonnement.

19
Chapitre I Généralités sur les Antennes Imprimées

Bibliographie

[I.1] [Link]/base.../antennes.../antennes-imprimées-e3310/
[I.2] Paul [Link] « Micro ondes [Link] passifs, propagation, antennes. Cours et exercice » Dunod,
paris 1997.
[I.3] Lotfi Djouane « Analyse rigoureuse d’une plaque microbande rectangulaire sur plan de masse ave
ouverture rectangulaire par la méthode des moments et les dyades hybrides de Green » mémoire de
magister en électronique de l’université de Batna.
[I.4] Mouhouche Faiza « Analysis of small microstrip patch antennas for mobile communication » mémoire
de magister en électronique de l’université de Boumerdes.
[I.5] [Link]
[I.6] Barkat Ouarda« Etude et synthèse des caractéristiques de réseaux d’antennes imprimées
supraconductrices dans la bande millimétrique » Thése de doctorat de l’université Mentouri
Constantine décembre 2009.
[I.7] Abla Behloul « Contribution à l'étude des antennes microstrip supraconducteur de forme triangulaire»
mémoire de magister en électronique de l’université Mentouri Constantine Avril 2010.
[I.8]« Antenna theory analysis and design » third edition .Constantine, A Balanis ajohn wiley a
[Link].
[I.9]Vivekananda Lanka Subramanian « Pattern analysis of the rectangular microstrip patch antenna » final
master degree thesis 30 ECTS, thesis No.: 4/2009 MSc. electrical engineering – communication &
signal processing January 2009.
[I.10] [Link]/aliouneba/pdf/[Link]
[I.11] Paul F .Combes« Micro ondes1 .lignes, guides et cavités, cours et exercices » Dunod, paris 1997.
[I.12] Chouti Leila« Contribution à l’étude d’antennes imprimés rectangulaires double bande et multi bandes
tenant compte de l’effet de couplage » Mémoire pour obtenir le diplôme de magister janvier 2009,
université Mentouri-Constantine.
[I.13] Hafedh ben Ibrahim Gaha «Analyse et conception des antennes fractales applications à la
télécommunication large bande » Thèse de doctorat juillet 2007.
[I.14] Addaci Rafik «Evaluation de la fréquence de résonance complexe et de la bande passante d’une
antenne microbande rectangulaire par la méthode neurospectrale » mémoire pour l’obtention du
diplôme de magister en électronique 16/10/2006, Université de Constantine.
[I.15] Hatem Rmili «Étude, réalisation et caractérisation d’une antenne plaquée en polyanilline
fonctionnant a 10 GHz » Thése pour obtenir le garde de docteur de l’université de Bordeaux I,
novembre 2004.
[I.16] Abdellatif Hafiane « Etude du couplage d’antennes imprimées par la méthode des différences finies et
extraction du modèle électrique» Mémoire pour le garde de docteur de l’université de Marne la
Vallee.Mars2003
[I.17] David Bidou «Contribution a l’étude d’antenne Fort Gain .Application Dans Le Domaine
Millimétrique» Thèse pour obtenir le garde de docteur de l’université de l’Imoges Mai 2003.
[I.18] Jerome Colinas «Etude et amélioration des performances d’une méthode de mesure en champ proche :
la méthode des courants magnétiques équivalents» Mémoire présenté en vue de l’obtention du
maîtrises et sciences appliquées (M-Sc-A) Août 2001, université de Montréal.

20
Chapitre II Méthodes d’Analyse des antennes imprimées

II.1-Introduction :
Etant donné la complexité de la modélisation des antennes imprimées et la grande variété des formes de
l'élément rayonnant (rectangulaire, circulaire, triangulaire…etc.), plusieurs méthodes basées sur des modèles
approchés ont été mises au point.
Les antennes en technologie micro ruban ont généralement un rayonnement bidimensionnel sur un
substrat diélectrique et peuvent donc être classées comme composants planaires. Les objectifs de l'analyse
d'antenne sont de prévoir les caractéristiques de rayonnement telles que des modèles de rayonnement, le
gain, et la polarisation aussi bien que des caractéristiques des champs proches, telles que l'impédance
d'entrée, la largeur de bande, le couplage, et l'efficacité de l’antenne.
La présence de l'inhomogénéité diélectrique rend l'analyse des antennes micro ruban plus complexe. De
plus, les conditions de frontière non homogènes, les caractéristiques étroites de bande de fréquence, la
grande variété d'alimentation et de la forme des plaques compliquées donne l’avantage à ce type d’analyse.
Parmi les méthodes analytiques principales: le modèle de ligne de transmission et le modèle de cavité ont été
utilisés. Ces méthodes maintiennent la simplicité. Les modèles analytiques emploient des suppositions
simplifiées, mais ils offrent généralement des solutions simples.

II.2- Différentes Méthodes d'analyse :

Deux types de méthodes permettent d’analyser les antennes imprimées [II.1]:
Les méthodes analytiques sont issues d’une approche physique du phénomène ou des hypothèses
simplificatrices permettent d’obtenir des équations assez simples mais en contrepartie, on obtient des
résultats approchés. Ces méthodes sont généralement basées sur les courants magnétiques équivalents le
long des bords du patch.
Les méthodes numériques exactes (full-wave) donnent des résultats plus précis au détriment de la perte de
la signification du phénomène physique. Elles nécessitent l’utilisation d’algorithmes numériques puissants et
un temps de calcul plus élevé que les méthodes analytiques. Elles sont généralement basées sur l’obtention
des distributions de courants électriques sur le patch et sur le plan de masse.
Nous pouvons rajouter que lorsque les antennes à analyser possèdent des géométries simples, elles
peuvent être facilement modélisées par des méthodes analytiques mais lorsque leur forme devient complexe
et non intuitive, il faut alors utiliser des méthodes numériques.

II.2.1- Méthodes analytiques : [II.1]

Deux méthodes analytiques permettent de décrire simplement le fonctionnement des antennes
imprimées résonnantes :
— le modèle de la ligne de transmission (valable notamment pour l’élément rectangulaire).
— le modèle de la cavité électromagnétique (valable notamment pour tout autre élément rayonnant de
géométrie quelconque simple).
Ces modèles permettent une étude simplifiée avec des formules directement utilisables et des résultats
satisfaisants. Cette approche fut celle utilisée par les premiers chercheurs à s’intéresser aux antennes
imprimées dans les années 1970.
Ces modèles ne sont pas tous équivalents. Certains permettent de décrire plus facilement le
fonctionnement de l’antenne en termes d’impédance (modèle ligne de transmission,) ou bien plus facilement
en termes de rayonnement (modèle de la cavité et des fentes rayonnantes).
Globalement, le modèle ligne de transmission est le plus simple mais aussi le moins précis, le modèle de
la cavité, plus précis, permet également une bonne approche physique du fonctionnement

21
Chapitre II Méthodes d’Analyse des antennes imprimées
II.2.1.1-Modèle de la ligne de transmission :
1-Introduction :
L’antenne est considérée comme une ligne de transmission dont les deux extrémités sont des ouvertures
(fentes) rayonnantes (figure II.1), ce modèle est complété en tenant compte des pertes dans le diélectrique,
dans les conducteurs et du couplage entre les deux ouvertures [II.1] [II.2] [II.3].
Pour arriver au modèle final, il faut tenir compte de certains phénomènes tels que [II.2]:
 L’alimentation.
 Le rayonnement.
 La permittivité effective du substrat.
 Les pertes diélectriques et de conductions.
 Les dimensions effectives de l'élément.

figure II.1 : Modélisation de l’antenne patch rectangulaire par deux fentes couplées
2- Modèle et caractéristiques :
Dans une première étape l'élément rayonnant est modélisé par une ligne de transmission de longueur L
chargé à ses deux extrémités par une admittance qui caractérise les discontinuités dues aux bords de
l’élément (figure II.2) [II.1] [II.2].

figure II.2 : schéma équivalent de l’antenne

 Le rayonnement : est le facteur le plus important .il est modéliser par grandeurs G et B, obtenues à
partir des formules suivantes [II.2] [II.4] [II.5]:
W  1 2
G 1   K0h   (II .1)
1200  24 
W h 1
B 1  0.636 ln  K 0 h   Pour  (II .2)
1200 0 10
Où 0 est la longueur d’onde dans le vide.
2
et K 0  : le nombre d’onde.
0
W : la largeur de l’élément rayonnant.
 L’alimentation : l’excitation est ramenée a l’antenne par une sonde coaxial, les mesures de
l’impédance s’effectuent au niveau de la sonde [II.2].
 Permittivité effective : du fait que la propagation se fait a la fois dans le diélectrique et dans l’air on
introduit une constante diélectrique effective  e relative à la discontinuité a l’interface air-
diélectrique.
Pour W  h  e : est donnée par James-Hall [II.2] [II.6] :

22
Chapitre II Méthodes d’Analyse des antennes imprimées

t
  1  r  1  10h 
0.5  r  1
e  r  h
1    (II .3)
2 2  W  W
4.6
h
 e : Permittivité effective.
 r : Permittivité relative.

h: la hauteur du substrat.
W : La largeur du substrat.

t : l’épaisseur de l’élément rayon

Ou, par Hammer stade [II.1] [II.2] [II.4] [II.5] [II.7] [II.8]:
0.5
  1  r  1  12h 
e  r  1   (II .4)
2 2  W 
Ainsi, les champs rayonnés par l’antenne micro ruban modélisée par le modèle de la ligne de transmission
sont définis par les relations [II.2]:
 Plan E :
 K 0W 
hWK0 E0e  Jk0 r
 sin sin  
2
E  j  cos  K W  (II .5)
2 r  0
sin  
 2 
 Plan H :
 KW 
sin  0 sin   
hWK 0 E0e  Jk0 r   2   cos  K 0 Le sin   (II .6)
E  2 j cos   
2 r  K 0W
sin    2 
 2 
 Phénomène de dispersion : pour les fréquences élevées, la permittivité du diélectrique varie avec la
fréquence. Cette variation est donnée par [II.2] [II.6]:
 
e  f   r  r e (II .7)
 f 
1 G  
 f 
 p
Zc
Avec : f p  (II .8)
2 0 h
Z c : Impédance caractéristique.

G : facteur empirique.

La figure ci-dessous donne la variation de  e en fonction de la fréquence :

 e quand f  0 
Avec :  e  f     (II .9)
 r quand f   

23
Chapitre II Méthodes d’Analyse des antennes imprimées

figure II.3: variation de  e en fonction de la fréquence

 Pertes : Les signaux qui se propagent sur une ligne micro ruban connaissent une atténuation causée
par la conductivité finie du conducteur et du plan de masse, et par les pertes dans le diélectrique.
Cette atténuation est exprimée par le coefficient d'atténuation qui à deux composantes, l'une
correspondant aux pertes dans le conducteur  c , et l'autre correspondant aux pertes dans le
diélectrique  d . Les expressions pour ces coefficients sont [II.6]:
f  r   e  f  1 
 d  27.3. tg   (II .10)
c  e  f    r 1 
tg : tangente de perte.
Pour W/h ≥1 :
2
R  32  We / h   h  1.25  2. .W  
 c  1.38. s   1  1 ln    (II .11)
Z c .h  32  We / h   We    t  
 . f .0
Avec : Rs  (II .12)

 : La conductivité des conducteurs.
We : La largeur effective donnée par [II.5]:
1.25.t   2.h  
We  W   1  ln   pour : W / h  1/ 2 (II .13)
   t 
1.25.t   4. .W  
We  W   1  ln    pour : W / h  1/ 2 (II .14)
   t 
 Longueur effective [II.2] [II.4] [II.8] [II.9]:
Le  L  2L (II .15)
Les deux extrémités de l’élément se comportent comme des continuités de type circuit ouvert
équivalent à deux capacités, chaque capacité peut être remplacée par une extension équivalente L
- Selon Lier [II.1] [II.2][II.7] [II.10]:
L  0.412.h. e
  0.3 . W / h  0.264  (II .16)
  e  0.258  . W / h  0.8 
-Selon Garg [II.2] [II.11]:
h  W / h  0.366    1  W 
L     0.28   0.274  ln   20518    (II .17)
2.  W / h  0.556   r   h 
Le modèle final de l’antenne est représente par la figure (II. 4)

24
Chapitre II Méthodes d’Analyse des antennes imprimées

figure II.4 : schéma équivalent final de l’antenne

3- calcul de l’impédance d’entrée :

Le modèle de l’antenne est une ligne de transmission chargé à ces deux extrémités par une résistance de
rayonnement R, cette ligne est caractérisée par son impédance Z c et ses pertes diélectriques et métalliques.
L’impédance caractéristique est donnée par [II.2] [II.4] [II.5]:
1
120.  We W 
Zc  .  1.393  0.667 ln  e  1.444   (II .18)
e  h  h 
L’impédance d’entrée est calculée en ramenant les deux résistances R au niveau du plan de la sonde
d’une distance L1 et L2. Le schéma équivalent devient [II.2]:

Z i : L’impédance ramenée au niveau de la sonde

Z c  R  Z c th  .Li  
Zi  (II .19)
Z c  [Link]  .Li 
Avec :     j  : Constante de propagation.
   d   c : Coefficient d’affaiblissement.
2. . f
 k   e  f  : Nombre d’onde.
c
On pose :
Z 1  [Link]    j   .Li 
yi  c  (II .20)
Zi z  th    j   .Li 
ye  y1  y2
On a donc :

1  [Link] .Li  z  th .Li  1   tg  .Li  

2
2
Re  ye     Ry (II .21)
i 1  z  th .Li    tg  .Li  . 1  [Link] .Li 
2 2 2

Im  ye  

tg  .Li .  z  th .Li   1  [Link] .Li 
2 2
 I (II .22)
 z  th .Li    tg  .Li  . 1  [Link] .Li 
2 2 y

25
Chapitre II Méthodes d’Analyse des antennes imprimées

 R 
z  Z 
 c

Avec :  L1  x  L  (II .23)
 L  L  2 L  L 
 2 1

 

L’impédance d’entré est donné par :

 R .Z  I .Z  
Z e   2y 0 2  j  X L  2y 0 2   (II .24)
 Ry  I y  
 Ry  I y 

4- Calcul de la fréquence de résonance: [II.2]

La fréquence de résonnasse est définie comme étant celle pour laquelle la partie réelle de l’impédance
d’entrée est maximal. Pour trouver la fréquence de résonance on cherche le MAX  Re  Z e  
II.2.1.2-Modèle de cavité :
Dans l’analyse des antennes imprimées par le modèle de la cavité, l’antenne est modélisée par une
cavité résonnante, dans laquelle on est capable de mettre en évidence les modes résonnants qui peuvent
s’installer. La nature des modes dépend de la forme du patch [II.12] [II.13] [II.14]. Chaque mode conduit à
une distribution de courant sur le patch, et le diagramme de rayonnement peut être calculé en faisant
rayonner ces courants. La nature des modes dépend de la forme du patch. Nous allons développer un modèle
pour un patch rectangulaire [II.13] [II.28]

figure II.5: Modèle de la cavité.

Les champs qui se trouvent à l’intérieur du substrat diélectrique sont ceux qui se trouvent dans une
cavité limitée par deux murs électriques horizontaux et quatre murs magnétiques latéraux [II.12] [II.16]
[II.11] [II.28].

figure II.6 : Représentation des murs électriques et magnétiques dans le modèle de la cavité
Afin de permettre une mise en équation des modes stationnaires, la cavité est placée dans un repère

 o, x, y, z  et présente les dimensions suivantes : a suivant l’axe des x, b suivant l’axe des y, h suivant l’axe
des z [II.15] [II.28].

26
Chapitre II Méthodes d’Analyse des antennes imprimées

figure II.7 : Représentation de la cavité dans un repère (O, x, y, z)

Les hypothèses qui vont être utilisées pour déterminer les modes qui s’installent dans cette cavité sont
les suivantes [II.13] [II.28].
- Les murs magnétiques dans les plans y = 0, y = b, x = 0, x = a sont idéaux

- La hauteur h du patch au-dessus du plan de masse est très inférieure  correspondant à la fréquence de
fonctionnement de l’antenne

Nous désignerons le champ électrique interne à la cavité porté par la direction OZ par : [II.15]
[II.16] [II.17] [II.28]. (Voir Annexe -A.1-)

 m   n 
E z  x , y   A cos  x  cos  y (II .25)
 a   b 

Les composantes de champ magnétiques s’obtiennent [II.15] [II.28]:

n  m   n 
Bx  x, y   A cos  x  sin  y (II .26)
jb  a   b 

m  m   n 
B y  x, y    A sin  x  cos  y (II .27)
j a  a   b 

La fréquence de fonctionnement de l’antenne [II.15] [II.17] [II.28] (voir Annexe -A .1-):

2 2
1  m    n 
f      (II .28)
2   a   b 

Nous proposons de donner une représentation de ces éléments dans le cas du mode fondamental (1, 0,
0) En remplaçant les valeurs de m = 1 et n = 0 dans l’expression des champs (II.25) (II.26) (II.27), et en
choisissant arbitrairement la constante A = 1, on obtient :

 
E z  x , y   cos  x  (II .29)
a 

H x x , y   0 (II .30)

  
H y x , y   j sin  x  (II .31)
wa a 

Les courants à la surface du patch sont donnés [II.13] [II.28]:

     
J x , y   j sin  x U x (II .32)
w a a 

27
Chapitre II Méthodes d’Analyse des antennes imprimées
 Calcul du champ rayonné :
Il existe plusieurs manières de faire un calcul analytique du champ rayonné dans la zone de
rayonnement lointain. Nous allons présenter dans cet partie deux méthodes, une très simple, mais valide
uniquement lorsque la permittivité relative du diélectrique est égale à 1 (le diélectrique est assimilé à l’air),
une autre un peu plus complexe, mais applicable avec n’importe quelle permittivité relative [II.13] [II.28].

Rayonnement des courants électriques (  r  1 ) :

Le modèle de la cavité a permis de déterminer l’expression des courants électriques qui prennent
naissance sur le patch
Dans un premier temps, nous calculons le champ lointain rayonné en espace libre au point P repéré
par ses coordonnées sphériques  r , ,   [II.13] [II.17] [II.28].

    
jk
  
E  p    R  J  x, y U U e jKOM .ds
4 s
(II .33)

S: représente la surface du patch sur laquelle on trouve les courants

0
  120    : est l’impédance d’onde du vide
0

  R   e  jKR / R (II .34)


U : est un vecteur unitaire dans la direction de propagation
 
K  KU : est le vecteur d’onde

figureII.8 : Représentation du référentiel et des coordonnées pour le calcul du champ lointain rayonné
par les courants.
Le champ total rayonné par l’antenne sur son plan de masse (voir Annexe -A .2-):

  
sin kh cos  ejVb 1 k e 1
jWa
 
E ptotal C
2j jV k2 W2
 
cos cos U  sin U   (II .35)

On représente souvent le diagramme de rayonnement dans deux plans principaux pour lesquels la
polarisation de l’onde est parfaitement déterminée.

a) Diagramme de rayonnement dans le plan E (   0 ) :

On a:   0  V  0 etW  k sin  : d’où l’expression simplifiée du champ rayonné dans le plan E, dans
laquelle la constante C a été modifiée pour prendre en compte tous les termes constants qui sont apparus :

28
Chapitre II Méthodes d’Analyse des antennes imprimées

  e jka sin  1 
E  P  planE  c sin  kh cos  
 cos  
2
cos  U    (II .36)

Le diagramme de rayonnement prend donc une expression particulièrement simple dans ce cas particulier

  sin  
E    cos   (II .37)
 2 

0.9

0.8

0.7

0.6
E()

0.5

0.4

0.3

0.2

0.1

0
-80 -60 -40 -20 0 20 40 60 80
theta(deg)

figure II.9 : Diagramme de rayonnement dans le plan E en champ

   sin   
2

E dB    10  LOG  cos    (II .38)

  2  

-5

-10

-15

-20
E()en dB

-25

-30

-35

-40

-45

-50
-80 -60 -40 -20 0 20 40 60 80
theta(deg)

figure II.10 : Diagramme de rayonnement dans le plan E en puissance

On peut noter que, sur le mode fondamental, le diagramme de rayonnement dans le plan E est
indépendant de la dimension a.

b) Diagramme de rayonnement dans le plan H    :

  
 2 

On a:    V  K sin  etW  0 : d’où l’expression simplifiée du champ rayonné dans le plan H,
2
dans laquelle la constante C a été modifiée une nouvelle fois pour prendre en compte tous les termes
constants qui sont apparus :
 e jkb sin   1 
E  P  planH  C sin  kh cos  
sin 
U   (II .39)

29
Chapitre II Méthodes d’Analyse des antennes imprimées
On en déduit l’expression du diagramme de rayonnement en champ :
  s in  b 
s in 
2

a 
(II .40)

E    cos
 s in  b
2 a

1
b/a=1
0.9 b/a=2
b/a=3
0.8
b/a=5
b/a=7
0.7

E(\théta) 0.6

0.5

0.4

0.3

0.2

0.1

0
-80 -60 -40 -20 0 20 40 60 80
théta(deg)

figure II.11: Diagramme de rayonnement dans le plan H en champ

   sin  b  
2

 sin   
  2 a  (II .41)
E dB     LOG    sin  b
10 cos 
 
 2 a 
 

0
b/a=1
b/a=2
-5 b/a=3
b/a=5
b/a=7
-10
E()en dB

-15

-20

-25

-30
-80 -60 -40 -20 0 20 40 60 80
theta(deg)

figure II.12 : Diagramme de rayonnement dans le plan H en puissance

Contrairement au plan E, le diagramme de rayonnement dans le plan H sur le mode fondamental
dépend de la dimension b, ou plus précisément du rapport b/a : plus cette dimension est grande, plus le
diagramme est directif. Au-delà de b = 2a, il y a apparition de lobes secondaires.
Rayonnement par application du théorème d’équivalence : [II.8] [II.28]

figureII.13 : Courants magnétiques équivalents qui entourent le patch rectangulaire

30
Chapitre II Méthodes d’Analyse des antennes imprimées
En posant AB = DC = a et AD = BC = b
  x 
-sur la portion AB: J M  E 0 cos( )U x (II .42)
a
 
- Sur la portion BC : J M   E 0U y (II .43)
  x 
- Sur la portion DC : J M   E 0 cos( )U x (II .44)
a
 
- Sur la portion AD : J M   E 0U y (II .45)
Le calcul du champ rayonné va s’effectuer à l’aide de l’intégrale de rayonnement généralisée au courant
magnétique :
       
E  p 
jK
4
 
  r   J M  U e j K OM dI  C  JM  U e j K
I I
  (II .46)

 
e j K OM  e  j Wx Vy
(II .47)
Enfin, le champ total rayonné s’obtient en sommant toutes les contributions on obtient (voir Annexe-A.3-)
  Vb   Wa     Vb   Wa  
 

sin  
 Wa   2   Vb 
cos  
2 
  
 
sin  
 Wa   2   Vb 
cos  
2    (II .48)
E  p  b cos   cos W sin    sin U  b cos   sin W sin    cos  cosU
  2  Vb  2  2  
2
   2  Vb  2 W2    a 
W    
 2 a   2  2 

a) Le diagramme de rayonnement dans le plan E :

  sin  
E    cos   (II .49)
 2  
 r 

   sin   
E dB    20  LOG  cos   (II .50)
 2  
  r  

L’évolution des diagrammes de rayonnement lorsque  r varie

0
r=2
-1
r=5
-2
r=7
r=9
-3

-4
E()en dB

-5

-6

-7

-8

-9

-10
-80 -60 -40 -20 0 20 40 60 80
theta(deg)

figureII.14 : diagrammes de rayonnement en puissance avec une variation de  r dans le plan E

Plus la permittivité est élevée, plus la dimension de l’antenne est petite en terme de longueur d’onde, et
plus le diagramme de rayonnement est large. Jusqu’à obtenir un diagramme quasi omni directionnel pour de
fortes valeurs de  r
31
Chapitre II Méthodes d’Analyse des antennes imprimées
b) Le diagramme de rayonnement dans le plan H:
L’expression générale du diagramme de rayonnement dans le plan H, en fonction du rapport des
dimensions du patch b / a :

  sin  b 
sin 
 2  a 
E     r  cos  (II .51)
 sin  b
2 r a

  sin  b 
sin 
 2  a 
E dB    20  LOG  sin  r b  cos (II .52)

2 r a

Les figures suivantes représentent l’évolution des digrammes, paramétrée par le rapport b/a, et par la
permittivité relative  r

0
b/a=1
-5 b/a=2
b/a=3
-10
b/a=5
b/a=7
-15

-20
E()en dB

-25

-30

-35

-40

-45

-50
-80 -60 -40 -20 0 20 40 60 80
theta(deg)

figure II.15 : Le diagramme de rayonnement dans le plan H avec variation de b/a lorsque  r  2
0
b/a=1
-5 b/a=2
b/a=3
-10
b/a=5
b/a=7
-15

-20
E()en dB

-25

-30

-35

-40

-45

-50
-80 -60 -40 -20 0 20 40 60 80
theta(deg)

figure II.16 : Le diagramme de rayonnement dans le plan H avec variation de b/a lorsque  r  5

32
Chapitre II Méthodes d’Analyse des antennes imprimées

0
b/a=1
-5 b/a=2
b/a=3
-10
b/a=5
b/a=7
-15

-20

E()en dB
-25

-30

-35

-40

-45

-50
-80 -60 -40 -20 0 20 40 60 80
theta(deg)

figure II.17 : Le diagramme de rayonnement dans le plan H avec variation de b/a lorsque  r  7
0
b/a=1
-5 b/a=2
b/a=3
-10
b/a=5
b/a=7
-15

-20
E()en dB

-25

-30

-35

-40

-45

-50
-80 -60 -40 -20 0 20 40 60 80
theta(deg)

figureII.18 : Le diagramme de rayonnement dans le plan H avec variation de b/a lorsque  r  9

Pour une valeur de permittivité donnée, plus le rapport b/a n’augmente, le diagramme de
rayonnement est plus directif. Au-delà qu’il y a apparition de lobes secondaires. Et aussi plus la valeur de
permittivité augmente, le diagramme de rayonnement est large jusqu’au obtenir un diagramme quasi omni
directionnel et les lobes secondaires est diminués aussi.
II.2.1.3-Comparaison des deux modèles analytiques :
Nous avons abordé les deux principaux modèles analytiques des antennes microruban, à savoir le
modèle de ligne de transmission et le modèle de cavité. La comparaison des deux modèles est donnée au
tableau suivant [II.2]:

Tableau II.1: Comparaison des deux modèles analytiques

33
Chapitre II Méthodes d’Analyse des antennes imprimées

II.2.2- Méthodes numériques exactes (full-wave) :

Il existe principalement quatre familles de méthodes numériques pour analyser les antennes imprimées.
Elles sont issues des méthodes utilisées en électromagnétisme : [II.18] [II.19] [II.20]
— la méthode des moments (MoM).
— la méthode des éléments finis (FEM).
— la méthode des différences finies (FDTD).
— la méthode de la matrice des lignes de transmissions (TLM).
Ces méthodes permettent d’obtenir des résultats très précis au détriment de temps de calculs longs car il
s’agit de discrétiser le volume ou la surface à analyser mais aussi le temps selon la méthode utilisée, ce qui
conduit à un très grand nombre d’inconnues. Les temps de calculs peuvent cependant être réduits si l’on
possède une puissance de calcul conséquente [II.1].

II.2.2.1--méthode des moments :

C’est une méthode d'analyse basée sur la résolution des équations intégrales des courants.
Dans le domaine cartésien, l’élément rayonnant est découpé en cellules rectangulaire supportant chacune un
courant de forme triangulaire (méthode MPIE: "Mixed Potential Integral Equation"). Cette discrétisation
permet de modéliser des formes complexes, à angle droit (dipôle replié, crois, encoche,…) [II.21].
Dans le domaine spectral, les fonctions inconnues sont développées sur un ensemble de fonctions de base
généralement globales. Cela permet d'obtenir un grand nombre d'équations algébriques que l'on peut intégrer
numériquement. Elles permettent d'arriver rapidement au résultat mais nécessitent une étude numérique
délicate (choix automatique des fonctions, intégrations,...) [II.21].
Une des méthodes, qui fournie l'analyse full-wave des antennes microbandes, est méthode de moments
(method of Moments, MOM). Dans cette méthode, les courants de surface sont employés pour modéliser le
patch et les courants de polarisation de volume sont employés pour modéliser les champs dans le substrat
diélectrique. Il a été montré par Newman [II.19] [II.3] et Ulyathan où une équation intégrale est obtenue
pour ces courants inconnus et l'usage de la méthode de moments, cette équation intégrale est convertie en
équations matricielle, ce qui peut alors être résolu par des diverses techniques algébriques pour fournir le
résultat.
Une vue générale de la méthode des moments décrite par Harrington [II.21] [II.22] et Balanis est
donnée ci-dessous [II.3]. La méthode des moments est une technique numérique qui permet de résoudre
efficacement le système d'équations intégrales en le transformant en un système matriciel résolu par
calculateur [II.22].
Elle est basée sur le critère de nullité d'une fonctionnelle constituée à partir d'une intégrale des résidus,
générée par la différence entre la solution approximative (fonction d'essai) et la solution exacte, pondérée par
des fonctions de poids (fonctions de test) [II.17].
La fonction d'essai est exprimée sous forme de série de fonctions de base connues dont les coefficients
de pondération sont déterminés en résolvant le système linéaire [II.17].
La forme de base de l'équation à résoudre par la méthode de moment est [II.22] [II.3] [II.23]:
F g  h (II .53)
Où F est un opérateur linéaire connu, g est une fonction inconnue, et h est la source ou fonction d'excitation.
Le but ici est de trouver g lorsque F et h sont connus. La fonction inconnue g peut être développée comme
une combinaison linéaire de N termes pour donner [II.22] [II.3]:
N
g   an g n  a1 g1  a2 g 2  .......  aN g N (II .54)
n 1

34
Chapitre II Méthodes d’Analyse des antennes imprimées

Où an est un constant inconnu et g n est une fonction connue appelée fonction de base. Substituant l'équation
(II.53) dans (II.54) et en utilisant la propriété de linéarités du l’opérateur F nous pouvons écrire [II.22]
[II.3]:
N

 a f  g  h
n 1
n n (II .55)

Les fonctions de base g n doivent être sélectionnées d'une telle manière, que chaque f  g n  Peut être
calculé. Les constantes inconnues an ne peuvent pas être déterminées directement parce qu’il y a N inconnus
et seulement une équation. Une des méthodes utiliser pour trouver ces constantes est la méthode des
résiduels pondérés [II.2] [II.3]. Dans cette méthode, un ensemble de solutions initiales est établi avec une
ou plusieurs variables. Les résiduels sont une mesure de la différence entre la solution initiale et la solution
vraie. Les variables sont choisis d'une manière à garantir un meilleur ajustement des fonctions initiale basées
sur la minimisation des résiduels. Ceci est effectué en définissant un ensemble de N fonction pondérée (ou
examinant) wm   w1 , w2 ,....., wN dans le domaine du l’opérateur F, prenant le produit de ces fonctions,
l'équation (II .57) devient [II.22] [II.3]:
N

a
n 1
n wm , f  g n   wm , h (II .56)

Où m = 1,2,…, N
En écrivant sous forme matricielle, nous obtenons [II.22] [II.21] [II.23]:
 f mn  an    hm  (II .57)
 w1 , f  g1  w1 , f  g 2  ....   w1 , h 
   a1 
a   
 w2 , f  g1  w2 , f  g 2  ....   w2 , h 
   2
 
 f mn   .  n  a3   hm    w3 , h 
 a   (II .58)
.   
.  . 
   
 w , f  g1  w , f  g 2  .... a  w , h
 N 
 
 N N
 N

Les constantes inconnues peuvent maintenant résolue en utilisant des techniques Algébriques telles que
la décomposition LU ou l’élimination gaussienne. On doit se rappeler que la fonction calculé doit être choisi
convenablement de sorte que les éléments wm soient non seulement linéairement indépendants mais
également réduisent au minimum les calculs requis pour évaluer le produit. Un choix possible des fonctions
d’essais calculé est de considérer ces dernières identiques aux fonctions de bases, c à d, wm  g n ceci est
connu sous le nom de méthode du Galerkin [II.19] comme décrit par Kantorovich et Akilov [II.3] [II.23].
Ainsi, la méthode de moment a été brièvement expliquée [II.3] [II.24].
 
 
Le but de cette méthode est de déterminer les courants de surface inconnus J s et M s .présent sur l’élément
rayonnant, le plan de masse et les murs magnétiques.
D’après L’équation de Richmond se présente sous la forme suivante :
       
J
 s t s t
.E  M .H dS   
 J a .Et  M a .H t d (II .59)
 
S V

Où S est la surface de l’élément rayonnant, V est le volume entourant la source d’alimentation, Et et H t


sont les champs « test » relatifs à une source électrique test quelconque placée sur l’élément rayonnant, J a et

M a sont les densités de courant électrique et magnétique d’excitation.

On considère que les plans conducteurs sont parfaits, par conséquent la densité de courant magnétique M s
sur ces plans est nulle. L’équation de Richmond se réduit à :

35
Chapitre II Méthodes d’Analyse des antennes imprimées
   
 J .E dS    J .E d
s t a t (II .60)

S V

La méthode des moments consiste à décomposer la densité de courant surfacique inconnue J s , en la

  
 
projetant sur la base J1 , J 2 ,...., J N , de sorte à pouvoir écrire [II.17] :
 
N
Js   In J n (II .61)
n 1
Si on considère le cas de N sources test à l’intérieur de la surface S, chacune d’elles va créer un

champ Em qui doit satisfaire l’équation de réaction de Richmond et qui s’écrit sous cette forme :
 N     
 
S  n1  I n J n  . E m dS    J a .Em d (II .62)
V

Ou encore :
N      
 I n  n
 J . E m dS n     J a .Em d (II .63)
n 1  Sn  V

Sn est la surface du niéme mode de décomposition.

 
Les fonctions J n n = 1,…, N, sont alors appelées Fonctions de base, et les champs Em , m = 1,……, N, sont
appelés Fonctions de test.
On peut formuler cette dernière équation de la manière suivante :
N

I Z n mn  Vm (II .64)
 
n 1

Avec : Vm   J a .Em d (II .65)

 
V

Et : Z mn    J n .Em dS n (II .66)

Sn

Les N sources forment ainsi un système de N équations linéaires, qui s’exprime sous la forme matricielle
suivante :
 Z  I   V  (II .67)
Dont les inconnus sont les coefficients I n du vecteur  I  , avec n=1,2,…, N.
La difficulté réside dans le choix des fonctions de base et des fonctions de test, ainsi que le choix de la
valeur de N.
Cette méthode est utilisée pour la détermination rigoureuse des grandeurs physiques suivantes :
 Les courants superficiels présents sur l’élément rayonnant, nécessaire pour l’évaluation des
champs rayonnés par celui-ci.
 L’impédance d’entrée de l’antenne.
 La fréquence de résonance.
Cette méthode se place dans le domaine temporel aussi bien que dans le domaine spectral, dans le quel
l’étude s’effectue dans l’espace transformé de Fourier. C’est une technique qui est très adaptée aux
structures multicouches constituées d’éléments rayonnants de formes simples. Un cas particulier de cette
méthode consiste à prendre les fonctions de bases égales aux fonctions test, ce cas particulier est appelé
méthode de Galerkin [II.9].
 Méthode de Green: [II.9]
Cette méthode est utilisée en général en association avec la méthode des moments, ou avec celle de
Galerkin.
En premier lieu, on cherche à déterminer la fonction de Green pour une courante unité localisée sur le

plan de l’antenne micro ruban. Donc, pour déterminer les champs électriques et magnétiques E et

36
Chapitre II Méthodes d’Analyse des antennes imprimées
 
H rayonnés à grande distance par l’antenne il suffit de connaître le vecteur potentiel A . Ce potentiel doit
cependant vérifier l’équation d’onde qui dirige son comportement.
La méthode de Green exprime ce vecteur potentiel sous forme d’intégrale d’une fonction de Green GA

calculé par la densité de courant inconnue J s Ceci se formule mathématiquement par :
 
A  x, y, z    J S  x ' , y ' , z ' .GA  x ' , y ' , z '  dx ' dy ' (II .68)
S
Ainsi, le champ électromagnétique rayonné est déterminé pour n’importe quel point à l’intérieur du
diélectrique. On obtient ainsi une équation intégro-différentielle dont l’inconnue est la densité de courant
superficielle, qui est résolue en forçant le champ électrique tangentiel total à être nul sur la surface de
l’antenne micro ruban.
Pour se faire, la méthode des moments ou de Galerkin est utilisée en choisissant les fonctions de base
et de test adéquates, afin de réduire le problème à une équation matricielle linéaire dont l’inconnue est la
distribution du courant sur l’antenne. La connaissance de cette dernière permet de déterminer les autres
propriétés de l’antenne microruban.
Cette méthode ne demande pas une très grande précision dans le calcul numérique des éléments de la
matrice du système. La convergence vers la solution dépend uniquement du choix des fonctions de base et
de test. Cette technique d’analyse est parfaitement applicable pour les antennes de forme circulaire,
triangulaire, elliptique,…etc. On peut facilement l’introduire dans la résolution de problèmes plus complexes
comme l’étude du phénomène de couplage entre élément adjacents du réseau et l’analyse des antennes
multicouches.

II.2.2.2- Matrices de ligne de transmission TML: [II.17]

Cette méthode numérique, la matrice des lignes de transmission ou Transmission Line Matrix en
anglais (T.L.M.), appliquées aux équations de Maxwell s’appuie sur les similitudes qui existent d’une part,
entre les équations de Maxwell vérifiées par les champs électromagnétiques d’une structure et d’autre part
les équations de Kirchhoff vérifiées par les courants et tensions d’un réseau de lignes de transmission.
La modélisation par la méthode des matrices de lignes de transmission implique deux étapes de base :
 Remplacement du problème électromagnétique par le réseau équivalent dérivant de l’analogie entre
le champ et les réseaux électrique
 Résolution du réseau équivalent par des méthodes itératives
Huygens, dans son traité de la lumière a défini le front d’onde comme la somme d’une infinité de
sources de rayonnements secondaires qui produisent, en trois dimensions, des ondes lettes sphériques. Leur
enveloppe constitue le front d’ondes. La méthode TLM discrétise les équations des télégraphistes appliquées
de façon locale. Grâce à la TLM, on appelle sources de rayonnement, les noeuds du domaine de calcul
maillés.
On réalise la mise en équation du problème en considérant le réseau maillé comme une série
d’intersections orthogonales de lignes de transmissions. Une cellule est reliée avec ses voisines par une
matrice S. Son calcul est itéré dans le temps par incrémentation, ce qui permet de suivre la propagation d’un
signal dans le réseau.
De plus, des améliorations sont depuis régulièrement apportées, telle que le maillage à pas variable,
l’extension aux milieux anisotropes, etc. On utilise cette méthode pour simuler des phénomènes de
propagation d’onde dans le domaine temporel.
Soulignons que la modélisation des matériaux composant une structure par un ensemble de
capacitances, d’inductances, et de résistances, réduit sa représentation réaliste car l’électromagnétique
utilisera davantage une méthode utilisant les grandeurs  ,  ,  et  .
Enfin, la méthode TLM ne présente pas de problème de convergence ni de stabilité. Elle est cependant
exigeante en place mémoire et en temps de calcul pour deux raisons principales :
- L’étude de zones où les champs électromagnétiques présentent de fortes discontinuités nécessite un
maillage plus fin, ce qui entraîne l’utilisation d’espace mémoire relativement important.

37
Chapitre II Méthodes d’Analyse des antennes imprimées
- Afin de limiter le phénomène dû à une réponse temporelle tronquée, on doit recueillir un grand nombre
d’impulsions sur une longue durée. Ceci nécessite un nombre d’itérations important entraînant une
augmentation du temps de calcul.
II.2.2.3-La méthode des différences finies :
1-Principe de la méthode :
La méthode F.D.T.D (Finite Difference Time Domain) ou la méthode des différences finies dans le
domaine temporel est une approche numérique permettant la résolution des équations différentielles dans le
domaine temps [II.19] [II.21].
Cette approche a pour but de mettre au point des méthodes de calcul permettant d’analyser la
propagation d’ondes hyperfréquences dans des lignes planaires. Nous avons choisi la simulation dans le
domaine temporel, selon la méthode des différences finies, car elle a l’avantage de permettre une
caractérisation dans une large bande de fréquence d’une structure en une seule simulation [II.19] [II.17]. En
appliquant la transformée de Fourier aux signaux temporels obtenus, il nous est possible de déterminer les
caractéristiques du guide, notamment les paramètres de propagation, l’impédance caractéristique en fonction
de la fréquence [II.17].
L’application de cette méthode aux équations de Maxwell dans l’espace libre a été introduite pour la
première fois par Yee en 1966[II.26] [II.24], la méthode consiste à approcher les dérivées ponctuelles
spatiales et temporelles qui apparaissent dans les équations de Maxwell par des différences finies centrées
[II.22] [II.25]. Le schéma proposé par Yee permet de surmonter la difficulté due à la dépendance simultanée
 
des champs électrique E et champs magnétique H entre eux, en effet, on obtient un schéma numérique
explicite permettant le calcul du champ électromagnétique dans tout le volume d’étude en fonction du temps.
Les composantes des champs électrique et magnétique sont décalées d’un demi-pas spatial et calculées de
manière alternative et itérative, respectivement à des multiples pairs et impairs du demi-pas temporel [II.17].
2-Point de départ : les équations de Maxwell:
Le point de départ de la méthode des différences finies est la forme différentielle des équations de
Maxwell
(II .71) et (II .72). Lorsque l’on considère un milieu linéaire, homogène, isotrope et sans perte, ces équations
ont la forme suivante [II.17] [II.22] [II.24] :

H 1 
   E (II .69)
t 

E 1 
   H (II .70)
t 
Lorsque ces équations sont projetées suivant les directions x, y et z, on obtient les équations suivantes [II.17]
[II.21] [II.26] :
     
Hx 1   Ez  E y   Ex 1   H z  H y 
    (II .71)     (II .72)
t   y z  t   y z 
     
H y 1   Ex  Ez   Ey 1   H x  H z 
    (II .73)     (II .74)
t   z x  t   z x 
     
Hz 1   E y  Ex   Ez 1   H y  H x 
    (II .75)     (II .76)
t   x y  t   x y 
Pour implémenter ces équations dans le programme, on doit donc les discrétiser à partir de l’approximation
de la dérivée centrée.

3-Principe des différences finies centrées : [II.17]

Prenons f(x), une fonction continue et dérivable, il est possible d’obtenir une approximation pour la
dérivée au point x0 , en se basant sur le développement en série de Taylor :

38
Chapitre II Méthodes d’Analyse des antennes imprimées
le développement en série de Taylor de la différenciation à droite est donné par :
 h h ' h 2 '' h3 '''
f  x0    f  x0   f  x0   f  x0   f  x0   ... (II .77)
 2 2 8 48
 h h h 2 '' h3
à gauche est donnée par : f  x0    f  x0   f '  x0   f  x0   (II .78)
 2 2 8 48
Maintenant on soustrait ces deux équations et on divise le résultat par h, ce qui nous donne le résultat suivant
 h  h
f  x0    f  x0  
2 2
f '  x0        h2  (II .79)
h
  h2  : représente l’erreur d’ordre 2 commise, sera négligé par la suite.
On appelle approximation centrée cette approximation de la dérivée. Les résultats qu’elle offre sont
plus précis en comparaison avec ceux donnés par d’autres types d’approximations dites droits ou gauches,
dont les formules sont décrites respectivement ci-dessous :
f  x0  h   f  x0 
f '  x0      h (II.80)
h
f  x0   f  x0  h 
f '  x0     h (II .81)
h
On remarque que le terme   h  qui est du premier ordre est moins précis en comparaison avec celui du
deuxième ordre de la dérivée centrée. Par conséquent, on utilisera l’approximation centrée dans notre étude
pour discrétiser les dérivées partielles, spatiales et temporelles présentes dans les équations de Maxwell.
On calcule la dérivée centrée d’une fonction (figure II.19) au centre d’un intervalle en se basant sur les
valeurs de la fonction aux extrémités.

figure II.19 : Point d’évaluation du calcul de la dérivée centrée

Par conséquent, les six équations de Maxwell seront exprimées sur chaque intervalle à partir de cette
approximation. On peut noter que les dérivées spatiales sont liées aux dérivés temporels. Elles feront
chacune l’objet d’une discrétisation particulière, soit par rapport à l’espace, soit par rapport au temps.

4-La discrétisation des équations de Maxwell :

4.1-La discrétisation spatiale: [II.17]

Abordons maintenant la façon dont discrétise les six équations de Maxwell. Pour cela, étudions l’une
des six équations de Maxwell :
  
H y 1   Ex  Ez 
   
t   z x 
Dans le calcul de H y , on fait intervenir la dérivée partielle de Ex par rapport à z, donc d’après la
définition de la dérivé centrée, le point où l’on calcule H y doit se trouver au milieu d’un segment parallèle à
l’axe oz , ayant comme extrémités deux points où Ez est connu. De même, le calcule de H y fait intervenir la

39
Chapitre II Méthodes d’Analyse des antennes imprimées

dérivée partielle de Ez par rapport à x . Donc le point ou l’on calcule H y se trouve également au milieu d’un
segment parallèle à ox ayant pour extrémités deux points où Ex est connu. En conséquence, H y représenté
sur la figure II.20, doit se trouver au milieu des deux points Ex et des deux points Ez
Notons que les valeurs du champ électrique et du champ magnétique seront calculées en des points
différents du maillage. Elles seront respectivement appelées noeuds électriques et noeuds magnétique.

figure II.20 : Circulation du champ E autour de H

figure II.21 : Disposition des noeuds électriques et des noeuds magnétiques dans le plan xoy
On vient de voir que les noeuds H y doivent se trouver entre deux noeuds Ex et entre deux noeuds Ez ,
de même les noeuds H x doivent se trouver entre deux noeuds Ez et entre deux nœuds E y (équations (II .71)
- (II .73)). Selon les équations (II .74)-(II .76), on constate également que le noeud E y , représenté sur la
figure II.22, doit se trouver entre deux noeuds H x et entre deux noeuds H z . Les noeuds Ex doivent se
trouver entre deux noeuds H y et entre deux noeuds H z . Les noeuds Ez doivent se trouver entre deux noeuds
H x et entre deux noeuds H y .

figure II.22 : Circulation du Champ H autour de E

L’arrangement des nœuds électriques et magnétiques doit donc respecter toutes ces conditions, et
conduit au schéma de la maille de Yee représenté par la figure II.23 suivantes [II.24] [II.27]

figure II.23 : Cellule de YEE

Les parallélépipèdes ou mailles élémentaires constituent le volume de calcul. Afin de le représenter
selon le schéma précédemment décrit, on doit construire un maillage pour la structure à étudier. Précisons
que dans le volume de calcul, sont toujours présents un nœud magnétique entre quatre noeuds électriques et
un noeud électrique entre quatre nœuds magnétiques. Ainsi la dérivée centrée est utilisée pour toutes les
dérivées spatiales présentes dans les équations de Maxwell. Pour représenter le volume de calcul, il est
nécessaire de construire un maillage [II.24].

40
Chapitre II Méthodes d’Analyse des antennes imprimées
4.2- Construction du maillage: [II.17]
Une discrétisation spatio-temporelle est nécessaire pour résoudre les équations de Maxwell (II .78).La
discrétisation spatiale s’effectue dans un volume nécessairement fini. Dans le cas d’un maillage régulier, les
dérivées spatiales sont évaluées dans les trois directions Ox , Oy, Oz avec des incréments constants dx, dy, dz ,
appelés pas spatiaux. Ces derniers sont choisis par l’utilisateur et dépendent de la plus petite longueur d’onde
présente dans la bande de fréquence d’analyse et de la géométrie de la structure à étudier. Le volume de
calcul est donc un parallélépipède comme le montre la figure II.24. Il est composé de nx, ny, nz cellules (ou
mailles) élémentaires de taille dx, dy, dz . On va associer trois noeuds électriques et trois nœuds magnétiques
pour chaque cellule élémentaire. Les valeurs du champ en ces noeuds seront notées
Ex  i, j , k  , E y  i, j , k  , Ez  i, j , k  , H x  i, j , k  , H y  i, j , k  , H z  i, j , k  et sont représentées ci-dessous dans la
maille de Yee figure II.24 où les entiers i, j, k représentent les indices de la cellule dans le maillage et
varient respectivement de 1 à nx 1 à ny et 1 à nz .

figure II.24 : Extraction d’une cellule élémentaire

La construction d’un maillage irrégulier est possible. Dans ce cas, les pas de discrétisation varient selon
les directions Ox , Oy et Oz . Par conséquent, on décrira le maillage à partir de trois tableaux représentant les
valeurs de ces pas.
 dx  i   i  1 nx (II .82)
 dy  j   j  1 ny (II .83)
 dz  k   k  1 nz (II .84)
Notons que dans ce cas, on doit prendre des précautions en raison des noeuds électriques qui ne seront
plus exactement entre deux noeuds magnétiques. Quant aux noeuds magnétiques, on les trouve toujours au
milieu des noeuds électriques, comme on peut le constater sur figure II.25. Pour ne pas s’éloigner de
l’hypothèse de la dérivée centrée, la variation entre les dimensions de deux mailles voisines ne doit pas
dépasser 20 %.

figure II.25 : Disposition des noeuds électriques et des nœuds magnétiques dans le plan oxy , pour le cas
d’un maillage irrégulier
Après s’être intéressé aux problèmes de discrétisation spatiale, nous allons maintenant aborder la
résolution des équations de Maxwell. Par conséquent, une discrétisation temporelle s’impose.

4.3- La discrétisation temporelle: [II.17]

Prenons (II .74) une des six équations de Maxwell :
  
H y 1   Ex  Ez 
   
t   z x 
Elle fait apparaître dans le membre de gauche la dérivée du champ magnétique par rapport au temps,
tandis que le membre de droite est considéré à une instante t. Si on prend en considération le principe de la

41
Chapitre II Méthodes d’Analyse des antennes imprimées
dérivée centrée, on en déduit que le membre de droite (le champ électrique) doit être calculé entre deux
instants successifs où on calcule le membre de gauche (le champ magnétique).
En tenant compte des six équations de Maxwell, on arrive à la conclusion que le champ électrique et le
champ magnétique ne doivent pas être calculés aux mêmes instants, mais à des instants décalés.
Pour le cas d’une discrétisation temporelle uniforme, avec un pas d’échantillonnage dt , le champ électrique
sera calculé pour des multiples pairs de dt / 2 , et le champ magnétique pour les multiples impairs de dt .,
Comme le montre la figure II.26 :

figure II.26: Calcul de H à l’instant ndt et calcul de E à l’instant (n+0.5) dt

4.4- Détermination du critère de stabilité : [II.17]
Le critère de stabilité numérique répond à deux impératifs, l’un spatial et l’autre, numérique. En effet,
les valeurs du champ en tout point du volume de calcul dépendent des valeurs obtenues aux instants
précédents. Ce schéma fonctionne sous la contrainte d’un critère de stabilité reliant le pas temporel dt aux
pas spatiaux dx, dy et dz
En raison de ce critère de stabilité, l’échantillonnage numérique doit être suffisamment fin pour pouvoir
suivre l’évolution temporelle du champ électromagnétique.
Par conséquent, afin qu’il n’y ait pas d’instabilités numériques, le pas de temps pour le calcul doit être plus
petit que le pas de temps correspondant à la propagation de l’onde sur une maille. On représente la condition
d’échantillonnage par la relation suivante [II.22] [II.26] [II.27]:
1
dt  (II .85)
1 1 1
c  
dx 2 dy 2 dz 2
Où c est la vitesse de la lumière, dx, dy, dz sont les pas de discrétisation spatiale.
Précisons que dans le cas où ce critère de stabilité n’est pas respecté, l’algorithme sera divergent. En
pratique, on choisit t assez proche de sa valeur maximum pour réduire autant que possible le nombre
d’itérations, donc la durée du calcul informatique.

4.5- Dispersion numérique : [II.17] [II.21]

Lorsqu’un signal électromagnétique se propage dans un domaine de calcul maillé par les différences
finies, il subit des transformations (distorsion, atténuation) dues, entre autres, aux effets dispersifs du
maillage. Ces effets sont dus à la discrétisation qui donne une représentation approchée des signaux, mais
aussi et surtout, aux précisions des formulations utilisées pour approcher les dérivés partiels. En d’autres
termes, cette dispersion dépend, d’une part, de la taille de la cellule dx, dy, dz par rapport à la plus petite
longueur d’onde présente dans le spectre d’analyse, et d’autre part, de l’ordre de l’erreur commise lors de
l’évaluation des dérivées partielles qui figurent dans les équations de Maxwell.
Pour quantifier la dispersion numérique de l’algorithme F.D.T.D, Taflove a fait une étude sur les deux
paramètres physiques : la vitesse de phase V et la vitesse de groupe Vg qui sont définies respectivement
comme suit :
  d  (II .86)
V  et V g   
k  dk 
  rad / s  : Pulsation du signal
A la différence de la vitesse de phase qui est la vitesse d’un front d’onde monochromatique, la vitesse de
groupe est la vitesse d’une impulsion riche en fréquences. Elle tient compte à la fois de la dispersion sur le
vecteur d’onde (anisotropie due au maillage), et de la dispersion fréquentielle.

42
Chapitre II Méthodes d’Analyse des antennes imprimées

En résumé, en choisissant une valeur d’incrément spatial inférieur à la valeur  /10 , il paraît juste de
dire que le phénomène de dispersion est négligeable.

II.2.2.4-La méthode des éléments finis :

 Théorie de la FEM :
La méthode des éléments finis est conceptualisée par A. Hrennikoff et R. Courant dans les années
1940 pour résoudre des problèmes de mécanique de structures. Quelques années plus tard, elle est introduite
en électromagnétisme pour résoudre l’équation de Maxwell. C’est un outil mathématique qui résout de
manière discrète les équations aux dérivées partielles [II.6] [II.29]. De manière générale, l’équation porte
sur une fonction spécifique définie sur un domaine et comporte des conditions aux bords permettant
d’assurer existence et unicité de la solution [II.29]. Les avantages de cette méthode proviennent de sa
capacité à s’adapter à des structures de formes géométriques relativement complexes [II.24] [II.29].
Cette méthode s’applique aux dispositifs micro ondes de formes quelconques [II.19]. La méthode des
éléments finis est une méthode rigoureuse, mais nécessite des temps de calcul importants ainsi qu’une
grande ressource de mémoire du calculateur. Il existe de nombreux ouvrages traitant du sujet [II.24].
Après la définition de la structure, on va dans un premier temps discrétiser le domaine de calcul en
morceaux homogènes. L’ensemble de ce grand nombre de petites régions est appelé le maillage [II.6]
[II.19] [II.29]. En 2D on préfère de discrétiser en triangles pour les éléments surfaciques et en 3D en
tétraèdres pour les volumes, en raison de la grande flexibilité et donc leur facilite de s’adapter aux structures
plus complexes [II.1] [II.2] [II.6]. L’avantage de la méthode des éléments finis est lié au fait que la forme
tétraédrique et la variation des dimensions des cellules élémentaires caractérisant le volume discrétisé,
confère au maillage une très grande souplesse [II.1] [II.19]. Cette méthode permet l’étude de structures
géométriques complexes et elle est donc parfaitement adaptée pour traiter les problèmes d’interconnexions
lors de la conception des modules hyperfréquences [II.6]. Par contre, elle nécessite de gros moyens
informatiques pour l’étude de structure complexes [II.2] [II.19].
Les éléments peuvent être très petits lorsque la géométrie est détaillée et plus grands ailleurs. Ces
milieux sont caractérisés par leur permittivité  r , leur perméabilité r ou éventuellement leur conductivité 
Ces grandeurs peuvent être réelles, complexes et/ou tensorielles .Le but de FEM est de calculer le champ
aux coins des éléments [II.6].
La figure suivante montre un exemple de maillage du plan transverse d’une carte électronique disposant
de deux pistes parallèles [II.6] [II.24]:

figure II.27: Maillage triangulaire adaptif en FEM

Dans chaque élément, une variation simple (souvent linéaire) du champ est assurée. Les coins des
éléments sont appelés des noeuds et le but de l’analyse est de déterminer la valeur des champs en chacun des
noeuds du maillage [II.24].
La résolution de l’équation aux dérivées partielles sur chaque noeud du réseau des éléments finis donne
une solution approximée par des fonctions d’interpolation on cherche une solution approchée suffisamment
fiable.

43
Chapitre II Méthodes d’Analyse des antennes imprimées
L’équation aux dérivées partielles ne peut pas être résolue directement dans le maillage [II.6]. Elle doit
pour cela être écrite sous une forme variationnelle. Cette forme contient les informations de l’équation aux
dérivées partielles et les conditions aux limites [II.6] [II.29].
La plupart des méthodes éléments finis sont des techniques dites variationnelles. Ces méthodes
fonctionnent en minimisant ou en maximisant une expression qui est connue pour être stationnaire autour de
la valeur exacte. Généralement, les méthodes éléments finis résolvent les champs inconnus en minimisant
une fonction d’énergie [II.6] [II.24].
En 3D, pour les problèmes harmoniques en temps, on peut écrire [II.6] :
H2 E 2 J .E
F     dV (II .87)
V
2 2 2 j
Les deux premiers termes représentent l’énergie électrique et l’énergie magnétique stockée dans le
volume, le troisième terme est l’énergie dissipée (ou appliquée) par les courants de conduction. Dans cette
expression on peut écrire H en fonction de E et J est la source donc connue [II.6].
Parallèlement, la solution est approximée sur un élément par des fonctions d’interpolation  i . La même
base de fonctions d’interpolation est utilisée pour tous les éléments qui sont très souvent une famille de
fonctions polynômes [II.6] [II.29]. Pour simplifier, les fonctions d’interpolation d’un élément sont nulles au-
delà de cet élément [II.6] [II.29]. La somme de ces fonctions d’interpolation pondérées par des
coefficients  i nous donne la solution sur un élément [II.6] [II.29] :
N
E   i i (II .88)
i 1

Ou N est le nombre de noeuds par élément,  i la valeur inconnue du champ aux nœuds et αi la fonction
d’interpolation de chaque nœud [II.6] [II.29].
Maintenant on peut réunir toutes les étapes précédentes pour obtenir une solution approximative. A ce
sujet, on va remplacer équation (II.88) dans équation (II.87).L’expression obtenue, on va minimiser à
chaque  i . De cette façon, on obtient un système linéaire d’équations dont les inconnues sont les
coefficients  i . Souvent un élément interagit seulement avec les éléments de son voisinage, il en résulte que
la matrice du système soit très souvent creuse. On utilise des méthodes spécifiques pour résoudre ce système
d’équations. La solution finale du champ est continue, approximée et connue sur l’ensemble du domaine
[II.6] [II.29].

figure II.28: Les différentes ´étapes de la méthode FEM

Un des grands avantages de FEM est sa capacité à traiter les structures complexes grâce au maillage
adaptable. Au contraire, la FDTD (finite différence time Domain méthode) est limitée à traiter les formes
rectangulaires. Autrement, au niveau de la programmation la FDTD est beaucoup plus facile à implémenter,
il suffit de discrétiser les équations d’onde [II.6].
Dans les années 1980 et 1990, beaucoup de laboratoires académiques et privés ont développé leurs
propres codes de calcul mais aujourd’hui, on trouve de très nombreux logiciels commerciaux performants,
les uns étant responsables des autres. Cependant l’emploi de ces logiciels ne dispense pas d’une bonne

44
Chapitre II Méthodes d’Analyse des antennes imprimées
connaissance de l’électromagnétisme et en particulier de la méthode d’analyse utilisée. L’ingénieur novice
doit notamment faire l’effort de l’appréhender.
Aujourd’hui, les efforts de recherche se portent sur le fait de combiner ces différentes méthodes pour
constituer des méthodes dites hybrides et multi échelles et ainsi obtenir de meilleurs résultats en termes de
vitesse de calcul et de précision : MoM-FDTD, MoM-FEM, multipôle GMT. Le tableau II.2 résume les
principales caractéristiques des méthodes détaillées dans ce paragraphe [II.1].

Tableau II.2 : Méthodes numériques pour l’analyse des antennes imprimées

II.3-Conclusion :
Dans ce chapitre nous avons présenté de façon générale les différentes méthodes et modèles pour
l’analyse d’une antenne imprimée, il y a des modèles analytiques simples tel que Modèle de la ligne de
transmission et modèle de cavité, et d’autre méthodes (techniques) numériques rigoureuses tel que : la
méthode des moments, la méthode des matrices des lignes de transmission, la FDTD et FEM qui son les
plus répandues. Ces technique basent principalement sur la discrétisation du volume à étudier ou la surface,
ses résultats son très précis.

45
Chapitre II Méthodes d’Analyse des antennes imprimées

Bibliographie
[II.1] [Link]/base.../antennes.../antennes-imprimées-e3310/
[II.2] Abla Behloul « Contribution à l'étude des antennes microstrip supraconducteur de forme triangulaire»

mémoire de magister en électronique de l’université Mentouri Constantine Avril 2010.

[II.3] Lotfi Djouane « Analyse rigoureuse d’une plaque microbande rectangulaire sur plan de masse ave
ouverture rectangulaire par la méthode des moments et les dyades hybrides de Green» mémoire de
magister.
[II.4] Vinod K. Singh et Zakir Ali« Dual band U- shaped microstrip antenna for wireless communication»
international journal of engineering science and technology ISSN: 0975-5462.
[II.5] Hatem Rmili « Étude, réalisation et caractérisation d’une antenne plaquée en polyanilline
fonctionnant A 10 GHz » Thèse pour obtenir le garde de docteur de l’université de Bordeaux I,
novembre 2004.
[II.6] [Link]/download/mth_32.pdf.

[II.7] Vivekananda Lanka Subramanian « Pattern analysis of ’the rectangular microstrip patch antenna »
final master degree thesis 30 ECTS, Thesis No.: 4/2009 MSc. electrical engineering –
Communication & Signal Processing January, 2009.
[II.8] [Link]/~trharris/files/text/[Link]
[II.9] Jerome Colinas« Etude et amélioration des performances d’une méthode de mesure en champ proche :
la méthode des courants magnétiques equivalents» Memoire présenté en vue de l’obtention du
maîtrises et sciences appliquées (M-Sc-A) Août 2001.
[II.10] Michel Clénet« Design of UHF circularly polarized patch antenna as a feed for a 9.1 metre parabolic
reflector» Defence R&D Canada -Ottawa June 2004.
[II.11] Mohammad A.A « Cavity model analysis of rectangular microstrip antenna operating in TM03 mode».

[II.12] [Link] « Analyse du rayonnement d’une antenne imprimée en forme d’anneau à l’aide de modèle
équivalents ».
[II.13] Vaudon Patrick-Master de recherche télécommunication hautes fréquences et optiques «XII.-Les
antennes imprimée» IRCOM- Université de Limoges.
[II.14] Moncef Kadi « Antennes miniatures sur des substrats de hautes permittivités ».

[II.15] Paul [Link] «[Link] passifs, propagation, antennes. Cours et exercice» Dunod,
paris 1997.
[II.16] «Antenna theory analysis and design» Constantine, A Balanis ajohn wiley a [Link].
[II.17] Abdellatif Hafiane «Etude du couplage d’antennes imprimées par la méthode des différences finies et
extraction du modèle électrique» Mémoire pour le garde de docteur de l’université de Marne la
Vallee.Mars2003.
[II.18] Mouhouche Faiza « Analysis of small microstrip patch antennas for mobile communication »
mémoire de magister en électronique de l’université de Boumerdes.
[II.19] Chouti Leila «Contribution à l’étude d’antennes imprimés rectangulaires double bande et multi
bandes tenant compte de l’effet de couplage »Mémoire pour obtenir le diplôme de magister janvier
2009.

46
Chapitre II Méthodes d’Analyse des antennes imprimées
[II.20] Hafedh ben Ibrahim Gaha «Analyse et conception des antennes fractales applications aux
télécommunications large bande» Thèse de doctorat juillet 2007.
[II.21] Michel Ney «Simulation électromagnétique -Outils de conception» Traité électronique : [E1030]
08/09/2008.
[II.22] Rodrigo Ponce Heredia« Contribution à la modélisation du rayonnement d’antennes conformées
pour des applications C.E.M. en aéronautique» Thèse pour obtenir le grade de docteur de
l’université de l’Umoges juillet 2008.
[II.23] Mustapha Himeur « Modélisation numérique pour La compatibilité électromagnétique de circuits
d’électronique de puissance» thèse pour obtenir le grade de docteur décembre 2002.
[II.24] Sylvain Demarty« Contribution à l’étude électromagnétique théorique et expérimentale des cartes de
circuit imprimé »Thèse pour obtenir le grade de Docteur de L’Université de Limoges, août 2006.
[II.25] David Bidou «contribution a l’étude d’antenne fort gain .Application dans Le domaine
millimétrique» Thése pour obtenir le garde de docteur de l’université de l’Imoges Mai 2003.
[II.26] Christos Kalialakis «Finite différence time domain analysis of microstrip antenna-circuit modules»
thèses de doctorat July 1999, université de Birmingham.
[II.27] Carlos Pereira «Etude avancée des canaux de transmission radio en contexte MIMO :
environnements complexes et couplage inter-antennes très large bande» thèse pour l'obtention du
grade de docteur de l'université de Poitiers, Aout 2006.

[II.28] Somia Boussaha et Mekhancha Sarra «Analyse du diagramme de rayonnement quasi

omnidirectionnel d’une antenne patch en utilisant le modèle de cavité» International conférence
systems and processing information May 15-17, 2011, Guelma.
[II.29] Guillaume Lunet « Radome actif utilisant des matériaux et structures A propriétés
électromagnétiques contrôlées» thèse pour l'obtention du grade de docteur de l'université de
Bordeaux, octobre 2009.

47
 Chapitre III Modélisation et simulation (par HFSS) d’une antenne
patch en utilisant les différentes techniques d’alimentation

III.1-Introduction :
La simulation est une phase indispensable et nécessaire qui nous donne une approximation des résultats
comme elle est réelle, pour choisir des bons et meilleurs paramètres avant la réalisation. Donc pour réaliser
un objet (une antenne imprimée, circuit, composant, carte imprimée…), nous devons certainement passer par
trois étapes nécessaires, elles sont comme suit :

 La modélisation: par des programmes ou des formulations.

 La simulation par des logiciels convenables.
 La réalisation (caractérisation par analyseur de réseau vectoriel,………)

Nous avons simulé les antennes imprimées à l’aide du logiciel Ansoft-HFSS.

III.2- Description du logiciel HFSS :

HFSS (High Frequency Structure Simulator) est un des outils d’analyse électromagnétique 3D basé sur
la méthode des éléments finis (FEM) parmi les plus aboutis. C’est un logiciel développé par la société
américaine Ansoft. Il présente une interface graphique avancée et de nombreuses fonctionnalités, il a
bénéficié de nombreuses optimisations pour accélérer les calculs. Il se révèle relativement bien adapté pour
décrire les cartes électroniques de manière précise mais nécessite des temps de calculs importants [III.1].
HFSS est un logiciel qui calcule le comportement électromagnétique d’une structure, pour analyser ce
comportement en détails, le logiciel met à notre disposition des outils d’interprétation post-traitement. Il
effectue une modélisation électromagnétique par résolution des équations de Maxwell à l’aide de la
méthode des éléments finis [III.2]. Le principe de la méthode utilisée consiste à diviser l’espace d’étude en
un grand nombre des petites régions (tétraèdres), puis à calculer localement le champ électromagnétique
dans chaque élément. En d’autre terme, les étapes de la FEM sont :
 Discrétisation en éléments finie de la région de résolution
 Discrétisation des équations pour chaque élément
 Assemblage de tous les éléments dans la région de résolution
 Résolution du système d’équation
HFSS utilise une méthode d’interpolation combinée avec un processus itératif dans lequel un maillage est
créé et automatiquement redéfini dans les régions critiques. Le simulateur généré une solution basée sur le
maillage initial prédéfini. Ensuite, il affine le maillage dans les régions ou il existe une haute densité
d’erreurs, et généré une nouvelle solution [III.3].
Le calcul des paramètres ( S11 par exemple) de la structure hyperfréquence suit les étapes suivantes [III.3]:
 Division de la structure en un nombre fini d’éléments.
 Excitation de chaque port de la structure avec une onde se propagent le long d’une structure guide
d’onde uniforme ou d’une ligne de transmission qui possède la même section que le port.
 Calcul de la configuration total du champ électromagnétique à l’intérieur de la structure

48
 Chapitre III Modélisation et simulation (par HFSS) d’une antenne
patch en utilisant les différentes techniques d’alimentation

figure III.1 : Division de la structure en tétraèdre

On peut ajouter aussi un balayage de fréquences (frequency sweep), de cette façon HFSS résout la
structure pour plusieurs fréquences mais toujours à partir du même maillage, le maillage de la fréquence
principale. Afin d’obtenir des bons résultats pour chaque fréquence il faut que la fréquence principale soit la
plus élevée. [III.4]
Un projet HFSS est un dossier qui contient un (ou plusieurs) modèle(s) appelé(s) design, chaque
modèle contient une structure géométrique, ses conditions aux limites et les matériaux utilisés, ainsi que les
solutions de champs électromagnétiques et les interprétations post-traitement. Tous ces composants
peuvent être atteints par la fenêtre Project Manager dans l'environnement HFSS [III.2].

III.3-Procédure de conception d’une antenne rectangulaire :

Sur la base de la formulation mathématique simplifiée, la procédure de modélisation est décrite, ce qui
conduit à des conceptions pratiques d’antenne microruban rectangulaire. Cette procédure suppose que les
informations spécifiées incluent le diélectrique du substrat  r , la fréquence de résonance f r , et la hauteur du
substrat h.

1) On a spécifie:  r , f r et h, et on a calculé les dimensions du patch W et L

Tel que : 0.0030  h  0.050 (III.1)

Calcul W :

c 2 (III.2)
W
2 fr r 1

 On détermine la permittivité effective de l’antenne utilisée :

1

  1  r  1  12h  2
e  r  1   (III.3)
2 2  W 

 On détermine L :

L  0.412.h.
  e  0.3 . W / h  0.264  (III.4)
  e  0.258  . W / h  0.8 
Calcul du L :

49
 Chapitre III Modélisation et simulation (par HFSS) d’une antenne
patch en utilisant les différentes techniques d’alimentation
c
L  2L (III.5)
2 fr  e

Ces étapes sont nécessaires pour constituer l’antenne imprimée rectangulaire.

2) On peut utiliser aussi Microstrip Patch Antenna Calculator [III.5] pour déterminer les dimensions
du patch W et L (voir Annexe -A.4-).
3) Lorsque l’antenne imprimée est alimentée par une ligne micro ruban, on peut utiliser Microstrip
Line Calculator pour déterminer les dimensions W et L de cette ligne [III.6] (voir Annexe- A.5-).
Les principaux inconvénients de l’alimentation par ligne micro ruban sont [III.3]:
- les couplages parasites dû à la réalisation de la ligne de transmission et du résonateur (patch) sur le même
substrat, et qui diminuent les performances de l’antenne au-delà de quelques GHz

- la nécessité d’ajouter un transformateur d’impédance de longueur pour obtenir l’adaptation entre la
4
ligne d’alimentation et le résonateur (patch) à 50 .

III.4-Techniques d’alimentation d’une antenne imprimée rectangulaire en utilisant

HFSS :
III.4.1-Alimentation par ligne micro ruban :
III.4.1.1-Présentation de l’antenne étudiée :
Le modèle d'antenne patch utilisé pour la simulation numérique avec Ansoft HFSS est représenté sur la
figureIII.2. L'antenne patch est conçue pour les 9 GHz sur un substrat d’alumine dont la permittivité relative
est 9.6 et d’épaisseur h de 0.635 mm. La largeur de la ligne microruban, permettant d’avoir une impédance
d’entrée de 50 Ω.

figureIII.2 : Antenne imprimé alimenté par ligne micro ruban

figure III.3 : Schéma de l’antenne, W : largeur ; L : longueur ; e : épaisseur du patch ; h : épaisseur du

substrat ; Wi : largeur de la ligne d’alimentation.

50
 Chapitre III Modélisation et simulation (par HFSS) d’une antenne
patch en utilisant les différentes techniques d’alimentation
Les tableaux III.1 et III.2 récapitulent les propriétés de l’antenne simulée

Matériaux Paramètres

Plan de masse : cuivre 5.8 107 S / m (conductivité)

Substrat : alumine  r  9.6 (permittivité relative) et Tan  0.001 (tangente des pertes)

Patch : cuivre 5.8 107 S / m (conductivité)

Tableau III.1 : propriétés des matériaux

Eléments aire (mm  mm) épaisseur (mm)

Substrat 15  15 h=0.635

Plan de masse 15  15 0.05

Patch (W  L) 7.23  5.23 0.1

Ligne micro ruban (W l  Ll ) 0.632  4.12 0.1

Tableau III.2: Paramètres de l’antenne.

Les résultats de la simulation sont représentés sur les figures suivantes :
Ansoft LLC S11 en dB HFSSDesign1 ANSOFT
0.00
Curve Info
dB(S(1,1))
Setup1 : Sw eep1

-5.00

-10.00
dB(S(1,1))

-15.00

-20.00

-25.00

-30.00
7.00 7.50 8.00 8.50 9.00 9.50 10.00
Freq [GHz]

figure III.4 : Variation du coefficient de réflexion en fonction de la fréquence

-La variation du coefficient de réflexion en fonction de la fréquence montre une résonance approximative à
la fréquence de f r  8.45GHz

-Le coefficient de réflexion S11  28dB

-La bande passante pour S11  10 dB s’étende de 8.091 à 8.904 GHz, en pourcent égale B p  %   9.6%

51
 Chapitre III Modélisation et simulation (par HFSS) d’une antenne
patch en utilisant les différentes techniques d’alimentation
Ansoft LLC Re(S11),Im(S11) HFSSDesign1 ANSOFT
0.25
Curve Info
im(S(1,1))
Setup1 : Sw eep1
re(S(1,1))
Setup1 : Sw eep1

0.00
Y1

-0.25

-0.50

-0.75
7.00 7.50 8.00 8.50 9.00 9.50 10.00
Freq [GHz]

figure III.5 : Représentation de la partie réelle et imaginaire du coefficient de réflexion en

fonction de la fréquence
Re ( S11 ) max= 8.87 GHz
Im ( S11 ) min=7 GHz, Im ( S11 ) max= 8.38 GHz
Cette figure représente la variation de la partie réelle et imaginaire du coefficient de réflexion normalisé
  R  jX
r   /  max
  r.  max : Cette relation permet de lire les valeurs de la partie réelle et imaginaire du coefficient de
réflexion à partir de l’abaque de Smith
 : C’est le coefficient de réflexion de forme complexe
 : C’est le module du coefficient de réflexion
r : Le coefficient de réflexion normalisé
D’après la figure III.4 :  max  28  28
Ansoft LLC S11 HFSSDesign1 ANSOFT

100 90 80 Curve Info

110 1.00 70 S(1,1)
120 60 Setup1 : Sw eep1

130 0.50 2.00 50

140 40

150 30

160 0.20 5.00 20

170 10

0.00 0.20 0.50 1.00 2.00 5.00

180 0.00 0

-170 -10

-160 -0.20 -5.00 -20

-150 -30

-140 -40
-0.50 -2.00
-130 -50
-120 -60
-110 -1.00 -70
-100 -90 -80

figure III.6 : Représentation du coefficient de réflexion sur l’abaque de Smith

Nous pouvons déterminer à partir l’abaque de Smith les paramètres suivants : le coefficient de réflexion
complexe et le ROS en n’importe quel point de cette courbe son connaitre la fréquence.

52
 Chapitre III Modélisation et simulation (par HFSS) d’une antenne
patch en utilisant les différentes techniques d’alimentation
Ansoft LLC Re(Z),Im(Z) HFSSDesign1 ANSOFT
80.00
Curve Info
im(Z(1,1))
Setup1 : Sw eep1
re(Z(1,1))
60.00 Setup1 : Sw eep1

40.00

20.00
Y1

0.00

-20.00

-40.00

-60.00
7.00 7.50 8.00 8.50 9.00 9.50 10.00
Freq [GHz]

figure III.7 : Représentation de la partie réelle et imaginaire de l’impédance en fonction de la

fréquence
A la fréquence de résonance la partie imaginaire est presque nulle par contre la partie réelle est environ de
50.84 
Ansoft LLC ROS en dB HFSSDesign1 ANSOFT
15.00
Curve Info
dB(VSWR(1))
Setup1 : Sw eep1

12.50

10.00
dB(VSWR(1))

7.50

5.00

2.50

0.00
7.00 7.50 8.00 8.50 9.00 9.50 10.00
Freq [GHz]

figure III.8 : Représentation du ROS en dB en fonction de la fréquence

A la fréquence de résonnance le ROS est minimal il est d’environ de 0.70 dB, d’autre terme quand le
coefficient de réflexion est proche de 1 (0 dB) le ROS est élevé et vis versa.

figure III.9 : Représentation 3-D du Gain total en dB

53
 Chapitre III Modélisation et simulation (par HFSS) d’une antenne
patch en utilisant les différentes techniques d’alimentation
La représentation en coordonnées polaires du gain en dB nous montre bien son variation en fonction
des angles d’Euler  et  .
Ansoft LLC Plan theta HFSSDesign1 ANSOFT

0 Curve Info
dB(GainTotal)
-30 30 Setup1 : Sw eep1
Freq='8.45631068GHz' Phi='0deg'
-6.60

-8.20
-60 60
-9.80

-11.40

-90 90

-120 120

-150 150

-180

figureIII.10 : Digramme de rayonnement de Gain total en dB (plan  )

Cette figure représente la variation du gain en dB en fonction de l’angle d’élévation  pour ce

cas   0 , la valeur maximale du gain est dans la direction   0 , il est égal à -6.93dB et il diminue si 
varie jusqu’à une valeur minimale pour    90 il est égal à-12.41 dB
Ansoft LLC Plan phi HFSSDesign1 ANSOFT

0 Curve Inf o
dB(GainTotal)
-30 30 Setup1 : Sw eep1
Freq='8.45631068GHz' Theta='90deg'
-12.28

-12.56
-60 60
-12.84

-13.12

-90 90

-120 120

-150 150

-180

figureIII.11 : Digramme de rayonnement en Gain en dB (plan  )

Cette figure représente la variation du gain en dB en fonction de l’angle d’azimute  , pour ce

cas   90 , la valeur maximale du gain est dans la direction   135 , il est égal à -12.08 dB et la valeur
minimale est pour   280 , il est égal à -13.26dB

figure III.12 : Représentation rectangulaire 3-D du Gain total en dB

54
 Chapitre III Modélisation et simulation (par HFSS) d’une antenne
patch en utilisant les différentes techniques d’alimentation
Aussi, la figure III.12 montre bien la variation du gain en fonction de  et  , on peut donc lire à partir de
n’importe quel point la valeur du gain

figureIII.13 : Distribution du courant surfacique sur le patch

La figure III.13 nous donne la répartition de la densité de courant surfacique, à une phase donnée et
montre la concentration du champ électromagnétique sur les bords du patch qui est parallèle à oy et au bout
de la ligne d’alimentation. Ce résultat confirme la théorie des antennes plaquées qui prévoit un rayonnement
par les bords du résonateur et la résonance de la structure excitée par la ligne microruban située. On peut
aussi distinguer les deux extrémités qui sont parallèle à ox à une distribution plus dense par rapport à toute la
surface du patch.

figureIII.14 : Distribution du champ électrique E sur le patch

La distribution du champ électrique sur le patch est maximale aux extrémités de ce dernier suivant oy
et minimale au son milieu. Cette remarque explique le rayonnement des deux fentes aux deux bords du
patch. Par contre la distribution du champ est maximale sur la ligne d’alimentation.

55
 Chapitre III Modélisation et simulation (par HFSS) d’une antenne
patch en utilisant les différentes techniques d’alimentation

figureIII.15 : Distribution du champ électrique E sur le substrat

La distribution du champ E sur le substrat est maximale aux dessous des deux fentes situés suivant oy
et minimale au son milieu. Est une distribution maximale sur la ligne.

figureIII.16 : Distribution du champ électrique E sur le substrat et le patch

Concernant la distribution du champ électrique sur le substrat et le patch nous pouvons dire que le
rayonnement du champ E dépasse la longueur du patch. Donc, ce rayonnement explique bien la longueur
effective et la permittivité effective que nous avons les inclut dans nos calculs.
III.4.1.2- Influence des paramètres géométriques sur l’antenne :
Dans cette partie, on va étudier les influences de plusieurs paramètres sur le rayonnement de l’antenne
tels que: la hauteur du substrat, la permittivité relative du substrat, la largeur du patch, la longueur du patch,
l’épaisseur du patch et la conductivité du patch. Ces paramètres géométriques influent directement sur les
performances de l’antenne à savoir: la bande passante, la fréquence de résonance et le niveau du coefficient
de réflexion.
III.[Link] -l’influence de la permittivité relative du substrat :
A cause de son influence sur la plupart des performances de l’antenne, la permittivité relative du
substrat doit être choisie pour avoir les propriétés radioélectriques désirées. La figure III.17 montre les
variations du coefficient de réflexion en fonction de la fréquence de travail. Dans notre simulation on a
choisi trois matériaux de permittivités relatives différentes: alumine (  r  9.6 ) (c’est le modèle de
référence), AI-N (  r  8.8 ), marble (  r  8.3 )

56
 Chapitre III Modélisation et simulation (par HFSS) d’une antenne
patch en utilisant les différentes techniques d’alimentation
Ansoft LLC différents valeurs de permittivité HFSSDesign1 ANSOFT
0.00
Curve Inf o
dB(S(1,1))
Setup1 : Sw eep1
dB(S(1,1))_1
Imported
dB(S(1,1))_2
Imported

-12.50
Y1

-25.00
permittivité=9.6
permittivité=8.8

permittivité=8.3

-37.50
7.00 7.50 8.00 8.50 9.00 9.50 10.00
Freq [GHz]

figureIII.17 : Variation du coefficient de réflexion en fonction de la fréquence

pour différentes valeurs de permittivité relative

 r  9.6  r  8.8  r  8.3

Fréquence de résonance GHz 8.45 8.79 9.02

Coefficient de réflexion S11 en dB -28 -31.37 -35.94

la perte diélectrique Tan 0.001 0 0

Etendu de la bande passante 8.091 -8.904 8.401-8.298 8.610-9.564

bande passante(%) 9.6 10.2 10.5

Tableau III.3 : Interprétation du résultat de la simulation pour différentes valeurs de  r

À noter qu’une diminution dans la valeur de  r produit une augmentation: de la fréquence de
résonnance, la bande passante et le niveau du coefficient de réflexion.
III.[Link] - l’influence de la hauteur du substrat :
Le substrat utilisé dans l’antenne microruban est caractérisé par son épaisseur. A cause de son influence
sur la plupart des propriétés de l’antenne, l’épaisseur du substrat doit être choisie pour avoir les propriétés
radioélectriques désirées. La figure III.18 montre les variations de S11 (coefficient de réflexion) en fonction
de la fréquence de travail. Dans notre simulation on a choisi trois valeur du substrat : h=0.7mm, h=0.635mm
(c’est le modèle de référence), h=0.55mm.
Ansoft LLC différent valeurs de h HFSSDesign1 ANSOFT
0.00
Curve Info
dB(S(1,1))
Setup1 : Sw eep1
dB(S(1,1))_1
Imported
-5.00
dB(S(1,1))_2
Imported

-10.00
Y1

-15.00
h=0.55mm

-20.00

-25.00 h=0.635mm
h=0.7mm

-30.00
7.00 7.50 8.00 8.50 9.00 9.50 10.00
Freq [GHz]

figureIII.18 : Variation du coefficient de réflexion en fonction de la fréquence

pour différent valeur de hauteur h
57
 Chapitre III Modélisation et simulation (par HFSS) d’une antenne
patch en utilisant les différentes techniques d’alimentation
0.7 mm 0.635 mm 0.55

Fréquence de résonance GHz 8.27 8.45 8.64

Coefficient de réflexion S11 en dB -30.41 -28 -15.89

Etendu de la bande passante 7.711-9.138 8.091 -8.904 8.394-8.917

bande passante(%) 17.25 9.6 6.05

Tableau III.4 : Interprétation du résultat de la simulation pour différentes valeurs de h

À noter qu’une diminution dans la valeur de h produit une augmentation dans la fréquence de
résonnance, une diminution de la bande passante et une diminution du niveau de coefficient de réflexion.
L’utilisation de minces substrats est conseillée pour les circuits micro-ondes parce qu’elle minimise les
ondes de surfaces, les radiations non désirées et le volume de l’antenne.

III.[Link]- l’influence de la largeur du patch :

Le patch utilisé dans l’antenne microruban est caractérisé par son largeur. A cause de son influence sur
la plupart des propriétés d’antenne, la largeur du patch doit être choisie pour avoir les propriétés
radioélectriques désirées. La figure III.19 montre les variations de S11 en fonction de la fréquence de travail
en prenant largeur du patch comme paramétré à varier. . Dans notre simulation on a choisi trois valeur de la
largeur du patch : W=8mm, W=7.23mm (c’est le modèle de référence), W=6mm.
Ansoft LLC différent valeurs de W HFSSDesign1 ANSOFT
0.00

-5.00

-10.00

-15.00

W=6mm
-20.00
Y1

-25.00
W=7.23mm

-30.00 Curve Inf o

dB(S(1,1))
Setup1 : Sw eep1
-35.00 dB(S(1,1))_1
Imported
W=8mm
dB(S(1,1))_2
Imported
-40.00

-45.00
7.00 7.50 8.00 8.50 9.00 9.50 10.00
Freq [GHz]

figureIII.19 : Variation du coefficient de réflexion en fonction de la fréquence

pour différent valeur de la largeur W

W=8mm W=7.23mm W=6mm

Fréquence de résonance GHz 8.43 8.45 8.50

Coefficient de réflexion S11 en dB -41.17 -28 -18.65

Etendu de la bande passante 8.091-8.870 8.091-8.904 8.125-8.969

bande passante(%) 9.24 9.6 9.91

Tableau III.5 : Interprétation du résultat de la simulation pour différentes valeurs de W

58
 Chapitre III Modélisation et simulation (par HFSS) d’une antenne
patch en utilisant les différentes techniques d’alimentation
On observe que lorsque la largeur du patch démunie la fréquence de résonnance augmentation, la
bande passante augmente et le niveau du coefficient de réflexion diminue.
III.[Link]- l’influence de la longueur du patch :
De même façon que la largeur, nous avons choisi trois valeur de la longueur du patch : L=5.23mm
c’est le modèle de référence, L=6mm, L=7mm. Les résultats de simulation sont donnés sur la figure
suivante :
Ansoft LLC différents valeurs de L HFSSDesign1 ANSOFT
0.00
Curve Info
dB(S(1,1))
Setup1 : Sw eep1
dB(S(1,1))_1
Imported
-5.00
dB(S(1,1))_2
Imported

-10.00
Y1

-15.00

L=7mm

-20.00
L=6mm

-25.00

L=5.23mm

-30.00
5.00 5.50 6.00 6.50 7.00 7.50 8.00 8.50 9.00
Freq [GHz]

figureIII.20 : Variation du coefficient de réflexion en fonction de la fréquence

pour différentes valeurs de la longueur L

L=5.23mm L=6mm L=7mm

Fréquence de résonance GHz 8.45 7.43 6.42

Coefficient de réflexion S11 en dB -28 -23.48 -20.83

Etendu de la bande passante 8.091-8.904 7.129-7.798 6.172-6.714

bande passante(%) 9.6 9 8.44

Tableau III.6 : Interprétation du résultat de la simulation pour différentes valeurs de L

On observe que l’augmentation de la longueur du patch, produit une diminution dans la fréquence de
résonnance, une diminution de la bande passante et une diminution du niveau de coefficient de réflexion.

III.[Link]- l’influence de l’épaisseur du patch :

Nous avons choisi trois valeur de l’épaisseur du patch : e=0.2mm, e=0.1mm c’est le modèle de
référence, e=0.05mm. Les résultats de simulation sont donnés sur la figure suivante :
Ansoft LLC différent valeurs de e HFSSDesign1 ANSOFT
0.00
Curve Info
dB(S(1,1))
Setup1 : Sw eep1
-5.00
dB(S(1,1))_1
Imported
dB(S(1,1))_2
-10.00 Imported

-15.00

-20.00
e=0.05mm
Y1

-25.00 e=0.1mm

-30.00

-35.00

-40.00 e=0.2mm

-45.00
7.00 7.50 8.00 8.50 9.00 9.50 10.00
Freq [GHz]

figureIII.21 : Variation du coefficient de réflexion en fonction de la fréquence

pour différent valeur de l’épaisseur e
59
 Chapitre III Modélisation et simulation (par HFSS) d’une antenne
patch en utilisant les différentes techniques d’alimentation
e=0.2mm e=0.1mm e=0.05mm

Fréquence de résonance GHz 8.56 8.45 8.43

Coefficient de réflexion S11 en dB -42.08 -28 -21.27

Etendu de la bande passante 7.824-9.890 8.091-8.904 8.126-8.810

bande passante(%) 24.13 9.6 8.11

Tableau III.7 : Interprétation du résultat de la simulation pour différentes valeurs de e

On observe que la diminution de l’épaisseur du patch, produit une diminution dans la fréquence de
résonnance, une diminution de la bande passante et une diminution du niveau de coefficient de réflexion.

III.[Link]-Effet de la conductivité du patch :

La théorie des antennes imprimées prévoit que l’énergie électromagnétique est rayonnée seulement à
travers les bords de la ‘cavité’ délimitée par le patch et le plan de masse. Ce qui signifie que le patch
(résonateur) doit être suffisamment conducteur pour éviter le rayonnement à travers sa surface. D’où l’intérêt
d’étudier l’effet de la conductivité du patch sur les performances de l’antenne. La figure III.22, présente la
variation du coefficient de réflexion S11 à l’entrée de l’antenne en fonction de la fréquence, pour des patchs
d’épaisseur e = 0.1mm, de permittivité  r  1 et de conductivités variables. Nous avons choisi trois
matériaux (conducteurs) du patch de permittivité égal à 1 et avec des conductivités différentes tels que:
  5.8  107 (cuivre) c’est le modèle de référence,   3.8  107 (aluminium),   1.5  107 (brass). Les
résultats de simulation sont donnés sur la figure suivante :

Ansoft LLC différent valeurs de conductivité HFSSDesign1 ANSOFT

0.00

-5.00

-10.00

conductivité=58000000
Curve Inf o
Y1

-15.00
dB(S(1,1))
Setup1 : Sw eep1
dB(S(1,1))_1
Imported
-20.00 dB(S(1,1))_2
Imported

-25.00 conductivité=38000000
conductivité=15000000

-30.00
7.00 7.50 8.00 8.50 9.00 9.50 10.00
Freq [GHz]

figure. III.22 : Variation du coefficient de réflexion en fonction de la fréquence

pour patch à conductivité différente

60
 Chapitre III Modélisation et simulation (par HFSS) d’une antenne
patch en utilisant les différentes techniques d’alimentation

  5.8  107   3.8  107   1.5  107

Fréquence de résonance GHz 8.45 8.46 8.46

Coefficient de réflexion S11 en dB -28 -25.99 -26.40

Etendu de la bande passante 8.091-8.904 8.089-8.911 8.085-8.911

bande passante(%) 9.6 9.7 9.7

Tableau III.8 : Interprétation du résultat de la simulation pour différentes valeurs de conductivité

On observe que la diminution de la conductivité du patch, produit une petite augmentation au niveau
du coefficient de réflexion, mais la fréquence de résonnance et la bande passante ne change pas.
Donc, la conductivité influe uniquement sur le niveau du coefficient de réflexion. Le meilleur S11 est atteint
pour une conductivité comparable à celle des métaux.

III.4.2-Alimentation par ligne micro ruban avec un transformateur d’impédance :


Le modèle de l’antenne à ligne micro ruban avec un transformateur de longueur (voir Annexe A.4)
4
réalisé par HFSS est donné par la figure III.23. Les paramètres de simulation de la structure sont les même
que ceux fixés lors la fabrication de l’antenne alimenté par ligne micro ruban. Ils sont donnés aux tableaux
III.9 et III.10 qui récapitulent les propriétés de l’antenne simulée.

Matériaux Paramètres

Plan de masse : cuivre 5.8 107 S / m

Substrat : alumine  r  9.6 et Tan  0.001

Patch : cuivre 5.8 107 S / m

Tableau III.9 : Propriétés des matériaux

Eléments aire (mm  mm) épaisseur (mm)

Substrat 15  15 h=0.635

Plan de masse 15  15 0.05

Patch (W  L) 7.23  5.23 0.1

Ligne micro ruban (W l  Ll ) 0.632  4.12 0.1

Transformateur d’impédance (W t  Lt ) 0.332  3.23 0.1

Tableau III.10: Paramètres de l’antenne.

61
 Chapitre III Modélisation et simulation (par HFSS) d’une antenne
patch en utilisant les différentes techniques d’alimentation
Le schéma de l’antenne à simuler selon HFSS est présenté dans la figure III.23:

figureIII.23 : Antenne imprimé alimenté par un transformateur d’impédance selon HFSS

Le résultat de la simulation est représenté sur les figures suivantes :
Ansoft LLC S11 en dB HFSSDesign1 ANSOFT
-2.50
Curve Info
dB(S(1,1))
Setup1 : Sw eep1

-5.00

-7.50

-10.00
dB(S(1,1))

-12.50

-15.00

-17.50

-20.00
7.00 7.50 8.00 8.50 9.00 9.50 10.00
Freq [GHz]

figure III.24 : Variation de coefficient de réflexion en fonction de la fréquence

- La fréquence de résonance f r  8.33GHz

-Le coefficient de réflexion S11  19.19 dB

-La bande passante BP  667 MHz ( BP  %   8.01% ).

Ansoft LLC Re (Z),Im (Z) HFSSDesign1 ANSOFT
120.00
Curve Info
im(Z(1,1))
Setup1 : Sw eep1
re(Z(1,1))
Setup1 : Sw eep1
100.00

80.00
Y1

60.00

40.00

20.00

0.00
7.00 7.50 8.00 8.50 9.00 9.50 10.00
Freq [GHz]

figure III.25 : Représentation de partie réelle et imaginaire de l’impédance en fonction de la

fréquence
62
 Chapitre III Modélisation et simulation (par HFSS) d’une antenne
patch en utilisant les différentes techniques d’alimentation
A la fréquence de résonance la partie imaginaire est nulle par contre la partie réelle est environ de 48.68
Ansoft LLC ROS en dB HFSSDesign1 ANSOFT
15.00
Curve Inf o
dB(VSWR(1))
Setup1 : Sw eep1

12.50

10.00
dB(VSWR(1))

7.50

5.00

2.50

0.00
7.00 7.50 8.00 8.50 9.00 9.50 10.00
Freq [GHz]

figure III.26 : Représentation du ROS en dB en fonction de la fréquence

La valeur du ROS est minimale à la fréquence de résonance f r  8.33GHz il est d’environ de 1.89 dB

figure III.27 : Représentation 3-D du Gain total en dB

Ansoft LLC S11 en dB HFSSDesign1 ANSOFT

0 Curve Inf o
dB(GainTotal)
-30 30 Setup1 : Sw eep1
Freq='8.330097087GHz' Phi='0deg'
-7.60

-9.20
-60 60
-10.80

-12.40

-90 90

-120 120

-150 150

-180

figure III.28 : Diagramme du Gain en dB (plan  )

Le gain est maximal pour   0 et égal à -6.55dB et il est minimal pour   90 et égal à -13.48dB.

63
 Chapitre III Modélisation et simulation (par HFSS) d’une antenne
patch en utilisant les différentes techniques d’alimentation
Ansoft LLC Plan phi HFSSDesign1 ANSOFT

0 Curve Info
dB(GainTotal)
-30 30 Setup1 : Sw eep1
Freq='8.330097087GHz' Theta='90deg'
-11.60

-12.20
-60 60
-12.80

-13.40

-90 90

-120 120

-150 150

-180

figure III.29 : Diagramme du Gain en dB (plan  )

Le gain est maximal pour   270 et égal à -11.21dB et il est minimal pour   150 et égal à -13.85dB

III.4.3-Alimentation par ligne enterrée :

L'antenne à ligne enterrée réalisée est une structure à deux substrats, avec une ligne microruban plaquée
sur le substrat inférieur qui se termine en un circuit ouvert sous le patch imprimé sur le substrat supérieur. Ce
couplage à proximité permet d'améliorer la bande passante et de réduire les rayonnements parasites.
Le modèle de l’antenne à ligne enterrée réalisé par HFSS est donné dans la figure III.30. Les paramètres de
simulation de la structure sont les même que ceux fixés lors la fabrication de l’antenne alimenté par ligne
micro ruban. Ils sont donnés aux tableaux III.1 et III.2.

figureIII.30 : Antenne imprimé alimenté par ligne enterrée selon HFSS.

figure III.31 : Schéma de l’antenne à ligne enterrée.

64
 Chapitre III Modélisation et simulation (par HFSS) d’une antenne
patch en utilisant les différentes techniques d’alimentation
Le largueur de la ligne microruban est Wl  0.632mm , permet d’avoir une impédance de 50 à
l’entrée de l’antenne. Une bonne adaptation est obtenue pour un microruban de longueur L  6.735mm .
Les deux substrats en alumine sont modélisés par deux couches diélectriques identiques d’épaisseur
0.635mm , la partie réelle de la permittivité est 9,6 et l’angle de perte  est égal à 0,001. Le plan de masse
est modélisé par une mince couche en cuivre d’épaisseur 0.05mm et de conductivité 5,8.107 S / m .
De même, le ruban d’alimentation est modélisé par une couche de cuivre d’épaisseur 0.1mm . La surface du
substrat ainsi que du plan de masse est de 15 15 mm2 .
D’autre part, les mêmes conditions de simulation numérique (port d’excitation, conditions aux limites
absorbantes, capacité de calcul…) utilisées pour l’antenne à alimentation directe (voir [Link].2), ont été
appliquées pour simuler l’antenne à ligne enterrée.
Le résultat de la simulation est représenté sur les figures suivantes :
Ansoft LLC S11 en dB HFSSDesign1 ANSOFT
0.00
Curve Info
dB(S(1,1))
Setup1 : Sw eep1

-5.00

-10.00

-15.00
dB(S(1,1))

-20.00

-25.00

-30.00

-35.00
6.00 7.00 8.00 9.00 10.00 11.00
Freq [GHz]

figure III.32 : Variation de coefficient de réflexion en fonction de la fréquence

- La fréquence de résonance f r  7.87GHz

-Le coefficient de réflexion S11  34.69dB

-La bande passante BP  590MHz ( BP  %   7.49% ).

Ansoft LLC Re(Z),Im(Z) HFSSDesign1 ANSOFT
300.00
Curve Info
im(Z(1,1))
Setup1 : Sw eep1
250.00
re(Z(1,1))
Setup1 : Sw eep1

200.00

150.00

100.00
Y1

50.00

0.00

-50.00

-100.00

-150.00
6.00 7.00 8.00 9.00 10.00 11.00
Freq [GHz]

figure III.33 : Représentation de partie réelle et imaginaire de l’impédance en fonction de la fréquence

A la fréquence de résonance la partie imaginaire est nulle par contre la partie réelle est environ de 50.10

65
 Chapitre III Modélisation et simulation (par HFSS) d’une antenne
patch en utilisant les différentes techniques d’alimentation
Ansoft LLC ROS en dB HFSSDesign1 ANSOFT
25.00
Curve Info
dB(VSWR(1))
Setup1 : Sw eep1

20.00

15.00
dB(VSWR(1))

10.00

5.00

0.00
6.00 7.00 8.00 9.00 10.00 11.00
Freq [GHz]

figure III.34 : Représentation du ROS en dB en fonction de la fréquence

La valeur du ROS est minimale à la fréquence de résonance f r  7.87GHz il est d’environ de 0.32 dB

figure III.35 : Représentation 3-D du Gain total en dB

Ansoft LLC Plan theta HFSSDesign1 ANSOFT

0 Curve Info
dB(GainTotal)
-30 30 Setup1 : Sw eep1
Freq='7.877022654GHz' Phi='0deg'
-2.00

-4.00
-60 60
-6.00

-8.00

-90 90

-120 120

-150 150

-180

figureIII.36 : Digramme de rayonnement en Gain en dB (plan  )

Le gain est maximal pour   0 et égal à -0.87dB et il est minimal pour   90 et égal à -7.99dB.

66
 Chapitre III Modélisation et simulation (par HFSS) d’une antenne
patch en utilisant les différentes techniques d’alimentation
Ansoft LLC Plan phi HFSSDesign1 ANSOFT

0 Curve Info
dB(GainTotal)
-30 30 Setup1 : Sw eep1
Freq='7.877022654GHz' Theta='90deg'
-7.10

-7.45
-60 60
-7.80

-8.15

-90 90

-120 120

-150 150

-180

figureIII.37 : Digramme de rayonnement en Gain en dB (plan  )

Le gain est maximal pour   90 et égal à -6.8dB et il est minimal pour   220 et égal à -8.33dB

III.4.4-Alimentation par câble coaxiale :

Le modèle d'antenne imprimé utilisé pour la simulation numérique dans Ansoft HFSS est représenté sur
la figureIII.38. L'antenne patch est conçue pour le 9GHz .Les paramètres de simulation de la structure sont
le même que ceux fixés lors la fabrication de l’antenne alimenté par ligne micro ruban. Ils sont donnés aux
tableaux III.11 et III.12 qui récapitulent les propriétés de l’antenne simulée.

Matériaux Paramètres

Plan de masse : cuivre 5.8 107 S / m

Substrat : alumine  r  9.6 et Tan  0.001

Patch : cuivre 5.8 107 S / m

Tableau III.11 : Propriétés des matériaux

Eléments aire (mm  mm) épaisseur (mm)

Substrat 15  15 h=0.635

Plan de masse 15  15 0.05

Patch (W  L) 7.23  5.23 0.1

Feed Rayon=1.6 0.5

Pin Rayon=0.7 0.685

Tableau III.12: Paramètres de l’antenne.

Le schéma de l’antenne à simuler selon HFSS est présenté dans la figure III.38:

67
 Chapitre III Modélisation et simulation (par HFSS) d’une antenne
patch en utilisant les différentes techniques d’alimentation

figureIII.38 : Antenne imprimé alimenté par câble coaxiale selon HFSS

Le résultat de la simulation est représenté sur les figures suivantes :
Ansoft LLC S11 en dB HFSSDesign1 ANSOFT
0.00
Curve Info
dB(S(1,1))
Setup1 : Sw eep1
-5.00

-10.00

-15.00

-20.00
dB(S(1,1))

-25.00

-30.00

-35.00

-40.00

-45.00
5.00 5.25 5.50 5.75 6.00 6.25 6.50 6.75 7.00
Freq [GHz]

figure III.39 : Variation de coefficient de réflexion en fonction de la fréquence

-- La fréquence de résonance f r  6.26GHz

-Le coefficient de réflexion S11  40.94 dB

-La bande passante BP  246MHz ( BP  %   3.92% ).

Ansoft LLC Re(Z),Im(Z) HFSSDesign1 ANSOFT
62.50
Curve Info
im(Z(1,1))
Setup1 : Sw eep1
re(Z(1,1))
Setup1 : Sw eep1
50.00

37.50
Y1

25.00

12.50

0.00

-12.50
5.00 5.25 5.50 5.75 6.00 6.25 6.50 6.75 7.00
Freq [GHz]

figure III.40 : Représentation de partie réelle et imaginaire de l’impédance en fonction de la fréquence

68
 Chapitre III Modélisation et simulation (par HFSS) d’une antenne
patch en utilisant les différentes techniques d’alimentation
A la fréquence de résonance la partie imaginaire est nulle par contre la partie réelle est environ de 50.55
Ansoft LLC ROS en dB HFSSDesign1 ANSOFT
35.00
Curve Info
dB(VSWR(1))
Setup1 : Sw eep1

30.00

25.00
dB(VSWR(1))

20.00

15.00

10.00

5.00

0.00
5.00 5.25 5.50 5.75 6.00 6.25 6.50 6.75 7.00
Freq [GHz]

figure III.41 : Représentation du ROS en dB en fonction de la fréquence

La valeur du ROS est minimale à la fréquence de résonance f r  6.26GHz il est d’environ de 0.15 dB

figure III.42 : Représentation 3-D du Gain total en dB

Ansoft LLC Plan theta HFSSDesign1 ANSOFT

0 Curve Inf o
dB(GainTotal)
-30 30 Setup1 : Sw eep1
Freq='6.26GHz' Phi='0deg'
-6.40

-7.80
-60 60
-9.20

-10.60

-90 90

-120 120

-150 150

-180

figure III.43 : Diagramme du Gain en dB (plan  )

Le gain est maximal pour   5 et égal à -5.15dB et il est minimal pour   90 et égal à -11.29dB.

69
 Chapitre III Modélisation et simulation (par HFSS) d’une antenne
patch en utilisant les différentes techniques d’alimentation
Ansoft LLC Plan phi HFSSDesign1 ANSOFT

0 Curve Inf o
dB(GainTotal)
-30 30 Setup1 : Sw eep1
Freq='6.26GHz' Theta='90deg'
-11.00

-11.20
-60 60
-11.40

-11.60

-90 90

-120 120

-150 150

-180

figure III.44 : Diagramme du Gain en dB (plan  )

Le gain est maximal pour   0 et égal à -10.85dB et il est minimal pour   105 et égal à -11.74dB

III.5-Antennes pour application radiocommunication :

Dans cette partie, nous avons simulé des antennes imprimées à élément rayonnant rectangulaire pour
des différentes applications. Les paramètres de simulation de la structure sont fixés lors de la fabrication de
l’antenne alimentée par ligne micro ruban. Les dimensions du patch (W et L) et celle de la ligne micro
ruban ont été optimisées à l’aide du logiciel Ansoft- HFSS pour avoir un coefficient de réflexion minimum
à l’entrée de l’antenne au voisinage de la fréquence demandée.
Le tableau ci-dessous montre quelques exemples d’application de l’antenne imprimée rectangulaires tels
que :
 Bande passante du GSM 900 est [Rx : 935-960MHz et Tx : 890-915].
 Bande passante du GSM 1800 est [Tx : 1710-1785 MHz et Rx : 1805-1880MHz]
 Bande passante de l’UMTS est [1.885-2.025GHz].
 Bande passante du WiFi est [2.412, 2.484 GHz].

Antenne GSM 1800 W  L  4.07mm , h  0.015mm ,  r  81

Antenne UMTS W  L  5.07mm , h  0.015mm ,  r  81

70
 Chapitre III Modélisation et simulation (par HFSS) d’une antenne
patch en utilisant les différentes techniques d’alimentation

Antenne WiFi W  4.07, L  5.07mm , h  0.015mm ,  r  81

Fréquence de Bande Gain

résonance passante dB
variation de coefficient de réflexion en fonction de la fréquence GHz

Ansoft LLC S11 en dB HFSSDesign1 ANSOFT

0.00
Curve Info
dB(S(1,1))
Setup1 : Sw eep1
-2.00

-4.00
f r  1.76GHz [1.71-1.81] -4.6
GSM -6.00
dB(S(1,1))

-8.00

1800 -10.00

-12.00

-14.00
0.50 1.00 1.50 2.00 2.50 3.00
Freq [GHz]

Ansoft LLC S11 en dB HFSSDesign1 ANSOFT

0.00
Curve Info
dB(S(1,1))
-1.25 Setup1 : Sw eep1

UMTS -3.25 f r  1.92GHz [1.91-1.93] -5.1

dB(S(1,1))

-5.25

-7.25

-9.25

-11.25
0.50 1.00 1.50 2.00 2.50 3.00
Freq [GHz]

Ansoft LLC S11 en dB HFSSDesign1 ANSOFT

0.00
Curve Info
dB(S(1,1))
Setup1 : Sw eep1
-2.00

WiFi -4.00
f r  2.43GHz [2.40-2.46] -4.49
dB(S(1,1))

-6.00

-8.00

-10.00

-12.00
0.50 1.00 1.50 2.00 2.50 3.00
Freq [GHz]

Tableau III.13: Exemples d’application de l’antenne imprimée.

71
 Chapitre III Modélisation et simulation (par HFSS) d’une antenne
patch en utilisant les différentes techniques d’alimentation

III.6-Conclusion:
A travers ce chapitre, on a essayé de présenter en premier lieu la simulation d’une antenne imprimée
rectangulaire en basant sur l’alimentation directe par ligne micro ruban. Ensuite, on a étudié l’influence des
différents paramètres géométriques de l’antenne imprimée sur les caractéristiques de rayonnement de
l’antenne a savoir: le coefficient de reflation, la fréquence de résonance, la bande passante,………..
En deuxième lieu, on a simulé le même modèle précédent mais avec d’autres types d’alimentations
comme : ligne micro ruban avec un transformateur, ligne enterré et par câble coaxiale.
Le résultat obtenu montre que l’alimentation par ligne microruban directe ayant une large bande
passante et un grand gain par apport aux autres types d’alimentation.

72
 Chapitre III Modélisation et simulation (par HFSS) d’une antenne
patch en utilisant les différentes techniques d’alimentation

Bibliographie
[III.1] Sylvain Demarty « Contribution à l’étude électromagnétique théorique et expérimentale des cartes de
circuit imprimé » Thèse pour obtenir le grade de Docteur de l’université de limoges le 31 août 2006

[III.2] Hamid Meghdadi, Jean-Charles Poux « Simulation 3D HFSS de filtres » Novembre 2005 – Mai 2006

[III.3]Hatem Rmili « Étude, réalisation et caractérisation d’une antenne plaquée en polyanilline fonctionnant
à 10 GHz » Thèse pour obtenir le garde de docteur de l’université de Bordeaux I, novembre 2004

[III.4] [Link]

[III.5] [Link]

[III.6][Link]
Operation=Synthesize&Wa=&W_units_list=Wmm&La=&L_units_list=Lmm

73
 Chapitre IV Modélisation et simulation d’un réseau d’antennes patch
en utilisant HFSS
IV.1- Introduction : [IV.1] [IV.2]
Le couplage mutuel est dû aux effets simultanés du rayonnement en espace libre et à la propagation des
ondes de surface. C’est un critère important qu’il faut donc considérer lors du calcul des caractéristiques des
réseaux. Le calcul théorique du couplage mutuel a été proposé par Kraus 1950 et Balanis 1982, il dépend
du type d’antenne et de la distance entre les éléments. Ainsi Jedlicka et Carver ont étudié
expérimentalement l’effet de couplage entre les antennes plaquées pour des géométries rectangulaires et
circulaires. Différentes méthodes ont été présentées pour calculer le coefficient de couplage entre les
antennes rectangulaires micro ruban. Ces méthodes ont été proposées par différents auteurs tels que
Sindoris, Pénard, Pozar, etc.

IV.2-Définition du couplage dans le plan E et le plan H : [IV.1] [IV.2] [IV.3]

Les diagrammes de rayonnement sont généralement représentés dans deux plans orthogonaux "plan E
et plan H", par rapport à la direction principale.

 Plan E : lieu des points de l'espace où le champ électrique rayonné est contenu dans ce plan.
 Plan H : lieu des points de l'espace où le champ magnétique rayonné est contenu dans ce
plan.
On distingue deux types de couplage voir figure IV.1:

figure IV. 1 : Couplage entre les antennes patch dans le plan E et H

 Couplage vertical ou couplage dans le plan H : c'est-à-dire le couplage entre deux éléments dans le
même substrat, le long de la direction de y avec une séparation d'accouplement Sh, les longueurs L des
patchs (dans la direction de x) sont les mêmes. La variation du champ électrique a lieu sur les côtés nommés
1, 2, 3,4, des patchs qui correspondent aux fentes rayonnantes voir la figure IV.2.

figure IV.2 : Configuration plan H

Les paramètres géométriques de l'antenne sont :
W: largeur du patch
L: longueur du patch
l : longueur de la ligne d'alimentation
D : distance qui sépare les deux patchs

74
 Chapitre IV Modélisation et simulation d’un réseau d’antennes patch
en utilisant HFSS
 Couplage horizontal ou couplage dans le plan E : c'est-à-dire le couplage entre deux éléments dans le
même substrat, le long de la direction de x avec une séparation d'accouplement Se, toutes les largeurs W des
patchs (dans la direction de y) ont les mêmes tailles. La variation du champ électrique a lieu sur les côtés
nommés 5, 6, 7,8 des patchs qui correspondent aux fentes rayonnantes voir la figure IV.3.

figure IV.3 : Configuration plan E

IV.3-Contribution à la modélisation du couplage :

Cette partie consiste donc à déterminer les paramètres électriques dans le cas d'un réseau de deux patchs
couplés.

IV.3.1-L’équation matricielle : [IV.2]

Le courant qui circule dans chaque antenne, induit des courants dans toutes les autres antennes
alimentées ou non. On peut modéliser le couplage mutuel entre les antennes par des paramètres tels que
l’impédance mutuelle, l’admittance mutuelle, le coefficient de transmission, etc.
En tenant compte de l’effet de couplage, nous aurons besoin de calculer le courant de chaque élément à
partir des tensions appliquées et de la matrice impédance du réseau.
Dans le cas pratique, les éléments rayonnants sont alimentés par des courants de même amplitude et de
même phase, les courants de surface créés sur chaque élément sont obtenus à partir de l’équation
matricielle suivante :

V    Z . I  (IV.1)
Dans laquelle :

 Z  : est la matrice d’impédance.

I : élément du vecteur courant  I  .

V : élément du vecteur excitation V  .

Les courants induits I ont été représentés à partir des tensions appliquées V à chaque élément. Dans les
cas pratiques, la matrice V  est généralement connue et nous pouvons alors déduire le courant à l'entrée de
chaque élément en connaissant la matrice d’impédance, donc on peut écrire :

 I   V . Z 
1

(IV.2)

IV.3.2-Modélisation par la méthode analytique : [IV.1] [IV.2]

On considère un réseau linaire formé de deux antennes disposées chacune sur l’un des axes à une
distance S comme l’indique la figure IV.4. L’antenne utilisée dans cette étude à une forme rectangulaire,
qui est simple à étudier due à la simplicité des formes des équations modélisant le rayonnement.

75
 Chapitre IV Modélisation et simulation d’un réseau d’antennes patch
en utilisant HFSS

figure IV.4 : Schémas de deux patchs

Le couplage mutuel entre deux antennes est modélisé en utilisant les paramètres Z, en admettant que les
modes sont quasi TEM. En effet, on associe aux deux antennes couplées un quadripôle décrit par sa matrice
d’impédance  Z  qui est défini de la façon suivante :

V 1  Z 11I 1  Z 12 I 2
V 2  Z 21I 1  Z 22 I 2 (IV.3)

Z 11 et Z 22 représentent respectivement l'impédance d'entrée de l'antenne 1 en présence de l'antenne 2 et

l'impédance d'entrée de l'antenne 2 en présence de l'antenne 1. Z 21 et Z 12 sont les impédances mutuelles et
traduisent couplage entre les deux antennes. L'impédance mutuelle Z 21 ( Z 12 ) entre les deux antennes
s'exprime en fonction de :

1
Z 12  
I 1I 2  E 2  J 1dxdy
(IV.4)

I1 et I2 sont les intensités des courants d'entrée qui alimentent respectivement l'antenne 1 et l'antenne 2. J1 est
la distribution du courant dans la première antenne. E2 est l'intensité du champ électrique produit le long de
la seconde antenne. I1 et I2 sont les amplitudes de courant dans le premier et le second patch.
Pour deux antennes patch, le schéma électrique équivalent qui prend en compte le couplage entre les patchs
est montré sur la Figure IV.5. Les éléments qui caractérisent le couplage sont Cc2 et Cc1.
Le modèle que nous proposons est constitué de deux circuits électriques équivalents d’antennes patchs.
Chaque circuit électrique équivalent d’antennes est constitué d’un circuit RLC en parallèle (une résistance
R1, une capacité C1 et une self L1) représente la première antenne, en parallèle avec (une résistance R2, une
capacité C2 et une self L2) représente la deuxième antenne, les deux antennes sont identiques (R1=R2=R,
C1=C2=C et L1=L2=L), les éléments qui caractérisent le couplage sontCc2 disposé horizontalement et Cc1
disposé verticalement.

figure IV.5 : Circuit électrique équivalent de deux patchs

En utilisant les transformations suivantes [I] = [Y]. [V], on assimile notre circuit à une cellule en Π avec :

Y 11 Y 1  i  C c 2
(IV.5)

76
 Chapitre IV Modélisation et simulation d’un réseau d’antennes patch
en utilisant HFSS
Y 12  i  C c 2
(IV.6)

Y 21 Y 12
(IV.7)

Y 22 Y 11
(IV.8)

1 1
Où : Y 1    i  C  C c 1 
iL  R (IV.9)

D’autre part, on applique les relations de passage pour exprimer les paramètres Z ij et S ij en fonction de
l’admittance:

 Z 11 Z 12  1 Y 22 Y 12 
Z Z   Y Y  (IV.10)
 21 22  Y  21 11 

On peut déterminer la matrice [S] à partir de la matrice [Y] du quadripôle:

2 y 12
S 12 
1  y 11 1  y 22   y 12 y 21 (IV.11)

S 11 
1  y 11 1  y 22   y 12 y 21
1  y 11 1  y 22   y 12 y 21 (IV.12)

2 y 21
S 21 
1  y 11 1  y 22   y 12 y 21 (IV.13)

S 22 
1  y 11 1  y 22   y 12 y 21
1  y 11 1  y 22   y 12 y 21 (IV.14)

Avec:

Y  yY 0
(IV.15)

1
Y0 
Z0 (IV.16)

Z 0 : correspond à l'impédance caractéristique de chaque accès

IV. 3. 3-Modélisation par la méthode numérique :

Dans ce paragraphe la disposition des deux antennes imprimées périodiques est déjà présentée sur la
Figure IV.4. La seule modification est que la première antenne est alimentée par une sonde coaxiale à la
position  x p , y p  et la deuxième par une ligne micro ruban Figure IV.6. Nous allons appliquer la méthode
des moments procédure de Galerkin pour déterminer les paramètres caractéristiques de couplage.

77
 Chapitre IV Modélisation et simulation d’un réseau d’antennes patch
en utilisant HFSS

figure IV. 6 : Exemple de couplage entre deux éléments d’une antenne micro ruban
En se basant sur la solution de la méthode de Galerkin qui consiste à développer la distribution du
courant sur le patch en série de fonction de base : [IV.4]
 N 
J s   In J n
n 1
(IV.17)

Où :

J : La densité surfacique de courant sur le patch.
 s
J n : La nime fonction de base selon x ou y.
I n : Amplitude du courant inconnue.
En résolvant l’équation :

V    Z . I  (IV.18)

Nous déterminons les I n .

on peut écrire :

Z mn    Em J n ds
(IV.19)
Sn


Vm    Em J i dv
(IV.20)
vi

Avec :

Em : Le champ électrique dû au courant J m du patch considéré.
J i : Le courant délivré par la source imprimée du patch.
 

  G  K , K .F  J  .F  J  .F  J  .F  J  dK dK
 
Z xx

mn xx x y x xn y xn x ym y ym x y
(IV.21)
 
Et les éléments du vecteur tension sont :
 
Vm    Q  K , K .  F . J  .F  J  exp  i.K x
 
v x y x ym y ym x p  i.K y y p  dK x dK y
(IV.22)
où:
iZ 0 K s K x   r  1 sin  K1d   iK x K1Tm sin  K1d 
2

Qv  
4 K 0  rTeTm (IV.23)
Z
x p  0 tan  r K 0 d
r
  (IV.24)

78
 Chapitre IV Modélisation et simulation d’un réseau d’antennes patch
en utilisant HFSS
0
Z0 
0 (IV.25)
 2 
K 0 : Nombre d’onde  K 0 
  
 Convergence des résultats numériques :
Les intégrales des éléments matriciels de l'impédance et du vecteur de tension sont évaluées
numériquement, et pour la convergence des résultats numériques on fait la conversion en coordonnées
polaires.
K x   cos 
K y   sin  (IV.26)

Les doubles intégrations infinies selon  K x , K y  sont changés aux coordonnées polaires double intégration
suivant   ,  
   2


 
dK x dK y     d d 
0 0
(IV.27)

Ce changement de variables, à pour avantage de transformer la double intégrale infinie en une seule intégrale
infinie et une autre finie, dans ce cas l’effet des ondes de surface sont inclus dans les calculs .
Le chemin d’intégration requis pour atteindre la convergence numérique est :
  0 /2
 0 (IV.28)

Pour réduire au minimum le calcul numérique, on choisit la fonction de symétrie suivante :

Fx   K x   Fx  K x  et Fy   K y   Fy  K y 
(IV.29)
Ce qui permet d’écrire les relations (IV.21) et (IV.22) comme suit :
 /2 
Z xx
4  G K , K  R  Fx  J xn  Fx  J ym    Re  Fy  J xn  Fy  J y   d  d 
mn xx x y e     (IV.30)
0 0

 /2 

  Q  K , K  Im  F  J  exp  i.K x  * R  Fy  J ym  exp  i.K x x p    d  d

  

Vm  4.i
0 0
v x y x ym x p e
(IV.31)

Les équations (IV.30) et (IV.31) sont les formes finales utilisées pour le calcul des éléments de la matrice
impédance et le vecteur tension.
La formule (IV.31) correspond à l’alimentation par une sonde coaxiale, cependant lorsque l’excitation est
We
faite par une ligne micro ruban, il faut la multiplier par le facteur
H
Si W est la largeur de la ligne d’alimentation, We est la largeur effective :
   0.30   W  0.262 H 
We  W  0.412  e  H
  e  0.258   W  0.813H  (IV.32)
Où  e est la permittivité effective de la ligne :
1   r  r  1  10 H 
e   1  
2 2  W  (IV.33)
 Ports d’impédance :
Considérons les deux éléments de l’antenne micro ruban de la géométrie de la Figure IV. 4. La relation
entre les ports de tension et les courants est :
79
 Chapitre IV Modélisation et simulation d’un réseau d’antennes patch
en utilisant HFSS
V P    Z P  .  I P 
(IV.34)
Où l’indice P est utilisé pour différentier entre les deux ports et aussi pour distinguer entre cette équation et
l’équation (III.20).
La matrice  Z P  est toujours une matrice 2× 2 (pour une géométrie à 2 ports).
Les ports d’impédance  Z P  nous permettent de déterminer l’impédance d’entrée et l’impédance de
couplage mutuelle.
Z11P est l’impédance d’entrée de l’élément 1 avec l’élément 2 en circuit ouvert, s’écrit sous la forme suivante :
N
Z11P   I nVn 
1

n 1
(IV.35)
Où I n sont les amplitudes des modes de courant trouvées à partir de (IV.18), et Vn1 sont les tensions dues à
la source 1.
L’impédance de transfert du circuit ouvert entre les ports 1 et 2 peut être écrite sous la forme :
  2    2
Z 21P  v
  E .J i dv
 Ii  (IV.36)
2

  2   2 
Où E est le champ électrique total causé par le courant de source J i du port 2, L’utilisation de l’équation
(IV.20) nous donne :
N
Z12P  Z 21P   I nVn
2

n 1
(IV.37)
IV.4-Les réseaux d’antennes :
Dans cette partie, nous allons traiter des méthodes de mise en réseau des éléments imprimés. Cette
application a connu un développement important depuis une vingtaine d’années car les éléments imprimés se
prêtent très bien à la réalisation de réseaux de plusieurs dizaines, ou plusieurs d’éléments, cela pour diverses
raisons [IV.5] [IV.6] :
-les faibles dimensions des éléments et leur simplicité de réalisation.
-leur aptitude à être imprimés sur des surfaces planes ou incurvées, facilement intégrables, par exemple, au
fuselage d’un avion ou d’une fusée.
-la possibilité de disposer les éléments rayonnants et le dispositif d’alimentation de part et d’autre d’un plan
de masse commun, en une structure très compacte, légère et peu coûteuse.
La conception et la réalisation de réseaux à éléments imprimés est un sujet très vaste, encore en pleine
évolution.
Les réseaux sont constitués d’un ensemble d’antennes élémentaires réunies par un distributeur
d’alimentation (figure IV.7) [IV.7] [IV.8].

figure IV.7: Réseau d’antenne schématisé avec son réseau de distribution.

80
 Chapitre IV Modélisation et simulation d’un réseau d’antennes patch
en utilisant HFSS
L’antenne élémentaire peut être un dipôle, une antenne imprimée (ou « patch »), un monopole, une
fente, un cornet …, Le réseau de distribution permet de pondérer soit en amplitude soit en phase chaque sous
élément du réseau ce qui permet une grande maniabilité au niveau des diagrammes de rayonnement [IV.7].

figure IV.8: Exemples de réseaux d'antennes a 16 éléments

Dans le cas d’un réseau imprimé, le gain maximal est de l’ordre de 28 dB qui avec une optimisation
peut atteindre les 30 dB. La bande de fonctionnement du réseau par rapport à l’élément seul est quant à elle
réduite et ne permet d’obtenir que 5% [IV.7].
Les antennes élémentaires fonctionnent à la même fréquence et sont disposées périodiquement selon une ou
plusieurs directions.
Cette mise en réseau des sources crée des rayonnements parasites (figure IV.9) qui sont dus au
débordement des champs magnétiques, à l'épaisseur de substrat utilisé, à l'orientation et à la distance entre
les sources [IV.9] [IV.10].

figure IV.9: Mécanisme des rayonnements parasites

Le groupement en réseau le plus simple est obtenu avec des sources identiques qui se déduisent les unes
des autres par translation pour former des réseaux linéaires et plans [IV.7] [IV.9] [IV.10].

figure IV.10: Architecture de base des réseaux d'antennes : a- Réseaux linéaires b- Réseaux plans
 Pour le réseau linéaire, on cherche à conformer le diagramme de rayonnement seulement dans le plan
contenant les sources. Lors d'une modification du diagramme de rayonnement sur l'ensemble de
l'hémisphère, les sources élémentaires doivent être disposées suivant le réseau bidimensionnel [IV.7]
[IV.11] [IV.10] [IV.8].
 Dans une antenne réseau, l'énergie est distribuée entre les diverses sources rayonnantes selon une loi
donnée : les caractéristiques de rayonnement du système dépendent à la fois du diagramme de
rayonnement de l'élément de base, des coefficients de réflexion en amplitude et en phase sur chaque
source et de la distance entre éléments [IV.7] [IV.11] [IV.10].
81
 Chapitre IV Modélisation et simulation d’un réseau d’antennes patch
en utilisant HFSS
La géométrie du réseau est déterminée grâce aux étapes suivantes : [IV.7] [IV.10]
 Le gain, la taille maximale de l'antenne, le nombre d'éléments rayonnants nécessaires.
 La distance entre deux éléments rayonnants successifs peut varier de 0,6 à 0,75 0 pour que la
directivité du réseau soit maximale. Si les éléments sont trop rapprochés, le phénomène de couplage qui se
crée fait diminuer la valeur maximale de la directivité. S'ils sont trop éloignés, des lobes de réseau
apparaissent.
 Ensuite, à partir de la théorie des lignes, il faut élaborer le circuit des lignes micro rubans alimentant
les différents éléments rayonnants. L'objectif de cette étape est : il faut réaliser l'adaptation 50 Ω de l'antenne
au niveau du point d'excitation.
 La dernière étape consiste à faire l'analyse du réseau et des lignes d'alimentation.
 Les réseaux imprimés présentent comme avantages : [IV.7] [IV.10]
- un faible encombrement du fait de l’épaisseur de l’élément répétitif face à la longueur d’onde (de l’ordre de

0 / 20 à 0 /100 ).

-Ils sont devenus plus faible coût du fait de leur banalisation dans le domaine des télécommunications ce qui
n’empêche pas des conceptions de circuits très fiables et très performantes que ce soit en gravure chimique
ou mécanique.
-Le fait que l’on puisse conformer le diagramme de rayonnement à souhait en fonction de l’application
envisagée est un avantage indiscutable des réseaux d’antennes imprimées par rapport aux autres types
d’antennes.
-La possibilité de générer de la polarisation circulaire est également très intéressant pour un grand nombre
d’applications telles que les communications, la navigation, les systèmes radar…
-leur aptitude à être imprimés sur des surfaces planes ou incurvées, facilement intégrables, par exemple, au
fuselage d’un avion ou d’une fusée.
 Toutefois, ce système présente des inconvénients : [IV.7] [IV.10]
-Le premier est la faible bande d’utilisation de 5% qui peut tout de même atteindre un maximum de 10%.
- De plus, le couplage entre les antennes élémentaires et les pertes générées par le réseau d’alimentation
réduisent considérablement le gain maximum à 30 dB. Au-delà, l’augmentation du gain est compensée par la
perte dans les lignes de transmission.
-Aux fréquences supérieures à 40 GHz, des problèmes de tolérances de réalisations et de connectique ne
permettent pas de comparer les résultats théoriques et expérimentaux. Outre ce problème à hautes fréquences,
le rayonnement parasite des différents composants passifs tels que les lignes de transmission, les coudes,
apparaît aux basses fréquences, qui est du au débordement des champs magnétiques, à l'épaisseur de substrat
utilisé, à l'orientation et à la distance entre les sources.

IV.4.1-Alimentation des réseaux d’antennes imprimées: [IV.10]

Le réseau d’alimentation permet d’exciter les éléments rayonnants d’une structure donnée de réseau
d’antennes, il réalise aussi la pondération complexe d’alimentation des sources. La précision suivant laquelle
est réalisée cette pondération influe directement sur les performances et les caractéristiques du réseau
d’antennes.
En effet, on est confronté à des problèmes d’encombrement, lorsque le nombre des sources devient
important, et cela peut même engendrer le problème de couplage entre les éléments rayonnants.
Dans le cas d’un réseau à grand nombre d’éléments, l’alimentation par lignes coaxiales devient très coûteuse
par ce que chaque source rayonnante doit être alimentée séparément par une sonde coaxiale.

82
 Chapitre IV Modélisation et simulation d’un réseau d’antennes patch
en utilisant HFSS
Il existe deux géométries d’alimentation :
 L’alimentation en série.
 L’alimentation en parallèle.

1-Réseaux linéaires à alimentation série: [IV.13] [IV.11]

figure IV.11: Réseau linéaire d’éléments rectangulaires à alimentation série

Ces réseaux constitués (figure IV.11) par des éléments rayonnants rectangulaires, reliés en chaîne les
uns aux autres par des tronçons de ligne microbande. Deux points homologues de deux éléments rayonnants
successifs sont séparés par une distance d . Comme la conductance d’un élément et, par conséquent, sa
puissance rayonnée est proportionnelle à sa largeur w , les largeurs des éléments diminuent, symétriquement,
en allant du centre O vers le bord pour avoir une loi d’éclairement du réseau anodisée. L’alimentation peut
être à partir d’un point tel que E par exemple, afin que les fentes qui bordent les éléments rayonnent en phase
ou avec le même gradient de phase. Elle peut aussi être effectuée à partir d’une extrémité A du réseau. Si la
fréquence est telle que d  e  0  e , tous les éléments rayonnent en phase et l’on a un réseau à
rayonnement transversal tandis que si d  e , il y a un gradient de phase et c’est un réseau à rayonnement
oblique.
Pour les guides d’ondes à fentes rayonnantes on peut avoir un fonctionnement en mode résonnant si les
extrémités libres sont en court-circuit ou en circuit ouvert . Ce fonctionnement, en ondes stationnaires, est
particulièrement indiqué si la fréquence est telle que : d  e .On peut aussi avoir un fonctionnement en
mode non résonnant si les extrémités libres sont fermées sur une charge adaptée. Ce fonctionnement, en
ondes progressives permet d’utiliser l’antenne aussi bien en rayonnement transversal d  e  qu’en
rayonnement oblique d  e  .

2-Réseaux linéaires à alimentation parallèle : [IV.5] [IV.6] [IV.8]

figure IV.12: Alimentation parallèle de 23éléments rayonnants par 3 étages de diviseurs de puissance

2 n éléments rayonnants peuvent être alimentés en parallèle (figure IV.12) pour n étages de 2 m 1 diviseurs
de puissance à chaque étage m . Afin que cette alimentation soit totalement adaptée, il est judicieux d’utiliser
des diviseurs de type Wilkinson.
Cette alimentation présente des avantages intéressants.
1) elle est à large bande puisque étant complètement symétrique à partir de l’entrée, les distances
parcourues jusqu’à chacun des éléments rayonnants sont les mêmes.

83
 Chapitre IV Modélisation et simulation d’un réseau d’antennes patch
en utilisant HFSS
2) Il est possible de disposer, sur les voies d’accès à ces éléments, des déphaseurs et des atténuateurs
(ou des amplificateurs) afin de faire varier leur phase et leur amplitude, pour conformer localement
ou globalement le diagramme de rayonnement.
3) Enfin, ce réseau d’alimentation peut être placé d’un coté du plan de masse, alors que les éléments
rayonnants sont placés de l’autre coté ; il suffit, pour cela, de se servir d’un dispositif d’excitation de
ces éléments
Par ailleurs, Munson a montré que les éléments rayonnants pouvaient être jointifs, de manière à former
une bande rayonnante alimentée par un dispositif parallèle (figure IV.13), sous la condition N S  W .

figure IV.13: Bande rayonnante à éléments jointifs à alimentation parallèle

N S est le nombre d’éléments. W  W e avec e  0  e . Cette structure est apte à être courbée pour
être enroulée en anneau autour d’un corps cylindrique, par exemple celui d’un missile ou d’une fusée.

3-Réseaux bidimensionnels : [IV.5] [IV.6]

Les réseaux linéaires permettent de conformer le diagramme de rayonnement dans un plan passant par
les sources du réseau. Afin d’obtenir un effet de réseau dans tout l’espace et, en particulier, dans les plans E

et H , il est nécessaire d’avoir recours à des réseaux bidimensionnels. Ils peuvent être obtenus en mode non
résonant à partir de réseaux linéaires à alimentation série (figure IV.14.a) ou parallèle (figure IV.14.b)

figure IV.14: Réseau bidimensionnel formé de réseaux linéaires, a- alimentation série, b- alimentation
parallèle
Dans le premier cas, il y a un gradient de phase aussi bien entre les sources disposées horizontalement
qu’entre celles disposées verticalement, ce qui donne un rayonnement oblique dans les deux plans.
L’utilisation des 4 voies d’entrée A, B, C et D permet, si nécessaire, d’obtenir 4 faisceaux rayonnés
simultanément.
Dans le second cas, les 4 réseaux linéaires sont alimentés en phase et l’inclinaison du faisceau rayonné
ne peut être obtenue que dans seul plan parallèle aux réseaux. Pour avoir un rayonnement oblique dans le
plan orthogonal. Il faudrait disposer des déphaseurs sur chacune des voies afin de créer un gradient de phase.
Avec la double alimentation parallèle, schématisée sur la figure IV.15, on produire un effet de réseau
dans deux plans orthogonaux avec, de plus, des polarisations orthogonales, voire circulaires, si les
excitations des éléments sont conformes à ce que nous avons décrit.

84
 Chapitre IV Modélisation et simulation d’un réseau d’antennes patch
en utilisant HFSS

figure IV.15: Réseau bidimensionnel à double alimentation parallèle

IV.5- Modélisation d’un réseau de deux patchs par HFSS :

IV.5.1-Réseau linéaire de deux patchs à alimentation série:
Ce réseau est constitué par deux éléments rayonnants rectangulaires identiques, reliés en chaîne les uns
aux autres par des tronçons de ligne microbande de 50 . Deux points homologues de deux éléments
rayonnants successifs sont séparés par une distance d  13.12mm (tel que d  e  0  e ), les éléments
rayonnants sont alimentés en phase, on utilise alors le patch du modèle de la ligne micro ruban étudié au
chapitre 3.
Les tableaux IV.1 et IV.2 récapitulent les propriétés de l’antenne simulée

Matériaux Paramètres

Plan de masse : cuivre 5.8 107 S / m

Substrat : alumine  r  9.6 et Tan  0.001

Patch : cuivre 5.8 107 S / m

Tableau IV.1 : propriétés des matériaux

Eléments aire (mm  mm) épaisseur (mm)

Substrat 30  30 h=0.635

Plan de masse 30  30 0.05

Patch (W  L) 7.23  5.23 0.1

Ligne micro ruban à 50 (W  L) 0.632  5.89 0.1

Tableau IV.2: Paramètres de l’antenne

 Couplage dans le plan H :
C’est le couplage entre deux éléments dans le même substrat le long de la direction y. La simulation
numérique avec Ansoft HFSS du modèle précédent est représentée sur la figure IV.16.

85
 Chapitre IV Modélisation et simulation d’un réseau d’antennes patch
en utilisant HFSS

figure IV.16 : Réseau linéaire de deux patches alimentés en série par ligne micro ruban (planH)
Le résultat de la simulation est représenté sur les figures suivantes :
Ansoft LLC S11 en dB HFSSDesign1 ANSOFT
-3.75
Curve Info
dB(S(1,1))
Setup1 : Sw eep1

-5.00

-6.25

-7.50
dB(S(1,1))

-8.75

-10.00

-11.25

-12.50
5.50 5.60 5.70 5.80 5.90 6.00 6.10 6.20
Freq [GHz]

figure IV.17: Variation du coefficient de réflexion en fonction de la fréquence

- La fréquence de résonance f r  5.76GHz

-Le coefficient de réflexion S 11  12.03dB

-La bande passante BP  137 MHz ( BP  %   2.37% ).

La représentation 3D et en coordonnées polaire du gain nous montre bien son variation en fonction des
deux angles  et 

figure IV.18 : Représentation 3-D du Gain total en dB

86
 Chapitre IV Modélisation et simulation d’un réseau d’antennes patch
en utilisant HFSS
Ansoft LLC Plan theta HFSSDesign1 ANSOFT

0 Curve Info
dB(GainTotal)
-30 30 Setup1 : Sw eep1
Freq='5.760517799GHz' Phi='0deg'
-14.60

-16.20
-60 60
-17.80

-19.40

-90 90

-120 120

-150 150

-180

figure IV.19 : digramme de rayonnement du Gain total en dB (plan  )

Le gain est maximal pour   0 et égal à -13.96dB et il est minimal pour   90 et égal à -20.64dB.
Ansoft LLC Plan phi HFSSDesign1 ANSOFT

0 Curve Info
dB(GainTotal)
-30 30 Setup1 : Sw eep1
Freq='5.760517799GHz' Theta='90deg'
-14.80

-17.60
-60 60
-20.40

-23.20

-90 90

-120 120

-150 150

-180

figure IV.20 : digramme de rayonnement en Gain en dB (plan  )

Le gain est maximal pour   270 et égal à -12.42dB et il est minimal pour   165 et égal à -24.6dB

 Couplage dans le plan E :

C’est le couplage entre deux éléments déposés sur le même substrat le long de la direction x. La
simulation numérique est représentée sur la figure IV.21.

figure IV.21 : Réseau linéaires de deux patches alimentés en série par ligne micro ruban (planE)
Le résultat de la simulation est représenté sur les figures suivantes :

87
 Chapitre IV Modélisation et simulation d’un réseau d’antennes patch
en utilisant HFSS
Ansoft LLC S11 en dB HFSSDesign1 ANSOFT
0.00
Curve Info
dB(S(1,1))
Setup1 : Sw eep1
-2.50

-5.00

-7.50

dB(S(1,1))
-10.00

-12.50

-15.00

-17.50

-20.00
5.50 5.75 6.00 6.25 6.50
Freq [GHz]

figure IV.22: variation du coefficient de réflexion en fonction de la fréquence

- La fréquence de résonance f r  6.05GHz

-Le coefficient de réflexion S11  19.14dB

-La bande passante BP  115MHz ( BP  %   1.9% ).

figure IV.23: représentation 3-D du Gain total en dB

Ansoft LLC Plan theta HFSSDesign1 ANSOFT

0 Curve Info
dB(GainTotal)
-30 30 Setup1 : Sw eep1
Freq='6.050161812GHz' Phi='0deg'
-12.40

-14.80
-60 60
-17.20

-19.60

-90 90

-120 120

-150 150

-180

figure IV.24 : digramme de rayonnement de Gain total en dB (plan  )

Le gain est maximal pour   50 et égal à -10.14dB et il est minimal pour   15 et égal à -21.36dB.

88
 Chapitre IV Modélisation et simulation d’un réseau d’antennes patch
en utilisant HFSS
Ansoft LLC Plan phi HFSSDesign1 ANSOFT

0 Curve Info
dB(GainTotal)
-30 30 Setup1 : Sw eep1
Freq='6.050161812GHz' Theta='90deg'
-15.20

-18.40
-60 60
-21.60

-24.80

-90 90

-120 120

-150 150

-180

figure IV.25:.digramme de rayonnement en Gain en dB (plan  )

Le gain est maximal pour   175 et égal à -13.63dB et il est minimal pour   70 et égal à -26.59dB

IV.5.2-Réseau linéaire de deux patchs à alimentation parallèle :

Deux éléments rayonnants identiques peuvent être alimentés en parallèle (figure IV.26). Afin que cette
alimentation soit totalement adaptée, il est judicieux d’utiliser des diviseurs de type Wilkinson [IV.14]. Les
éléments rayonnants sont en phase.
Les tableaux IV.3 et IV.4 récapitulent les propriétés de l’antenne simulée

Matériaux Paramètres

Plan de masse : cuivre 5.8 107 S / m

Substrat : alumine  r  9.6 et Tan  0.001

Patch : cuivre 5.8 107 S / m

Tableau IV.3 : propriétés des matériaux

Eléments Aire (mm  mm) épaisseur (mm)

Substrat 20  20 h=0.635

Plan de masse 20  20 0.05

Patch (W  L ) 7.23  5.23 0.1

Ligne micro ruban à 50 (W  L ) 0.632  4.12 0.1

Ligne micro ruban à 70.7 0.92  3.32 0.1

(W  L )

Ligne micro ruban à 100 0.09  11.885 0.1

(W  L )

Tableau IV.4: Paramètres de l’antenne.

89
 Chapitre IV Modélisation et simulation d’un réseau d’antennes patch
en utilisant HFSS

figure IV.26 : Réseau linéaire de deux patches alimentés en parallèle par diviseur de Wilkinson
Le résultat de la simulation est représenté sur les figures suivantes :
Ansoft LLC S11 en dB HFSSDesign1 ANSOFT
-7.50
Curve Info
dB(S(1,1))
Setup1 : Sw eep1

-8.00

-8.50
dB(S(1,1))

-9.00

-9.50

-10.00

-10.50
5.00 5.25 5.50 5.75 6.00 6.25 6.50
Freq [GHz]

figure IV.27: variation de coefficient de réflexion en fonction de la fréquence

- La fréquence de résonance f r  5.6GHz

-Le coefficient de réflexion S11  10.46dB

-La bande passante BP  396 MHz ( BP  %   7.07% ).

figure IV.28: représentation 3-D du Gain total en dB

90
 Chapitre IV Modélisation et simulation d’un réseau d’antennes patch
en utilisant HFSS
Ansoft LLC Plan theta HFSSDesign1 ANSOFT

0 Curve Info
dB(GainTotal)
-30 30 Setup1 : Sw eep1
Freq='5.601941748GHz' Phi='0deg'
-8.00

-10.00
-60 60
-12.00

-14.00

-90 90

-120 120

-150 150

-180

figure IV.29 : digramme de rayonnement du Gain total en dB (plan  )

Le gain est maximal pour   0 et égal à -5.64dB et il est minimal pour   90 et égal à -15.92dB.
Ansoft LLC Plan phi HFSSDesign1 ANSOFT

0 Curve Inf o
dB(GainTotal)
-30 30 Setup1 : Sw eep1
Freq='5.601941748GHz' Theta='90deg'
-14.10

-14.70
-60 60
-15.30

-15.90

-90 90

-120 120

-150 150

-180

figure IV.30:.digramme de rayonnement en Gain en dB (plan  )

Le gain est maximal pour   255 et égal à -13.9dB et il est minimal pour   20 et égal à -16.25dB

IV.6- Modélisation d’un réseau à quatre patchs alimentés en parallèle:

Quatre éléments rayonnants identiques peuvent être alimentés en parallèle en utilisant les diviseurs de
type Wilkinson (figure IV.31). Deux points homologues de deux éléments rayonnants successifs sont séparés
par une distance d  13.12mm . Les éléments rayonnants sont en phase, on utilise le même modèle précédant.

figure IV.31 : Réseau linéaires de quatre patches alimentés en parallèle par diviseur de Wilkinson
Les tableaux IV.5 et IV.6 récapitulent les propriétés du réseau simulé

91
 Chapitre IV Modélisation et simulation d’un réseau d’antennes patch
en utilisant HFSS
Matériaux Paramètres

Plan de masse : cuivre 5.8 107 S / m

Substrat : alumine  r  9.6 et Tan  0.001

Patch : cuivre 5.8 107 S / m

Tableaux IV.5 : propriétés des matériaux

Eléments Aire (mm  mm) Epaisseur (mm)

Substrat 50  30 h=0.635

Plan de masse 50  30 0.05

Patch (W  L ) 7.23  5.23 0.1

Ligne micro ruban à 50 (W  L ) 0.632  4.12 0.1

Ligne micro ruban à 70.7 0.92  3.32 0.1

(W  L )

Ligne micro ruban à 100 0.09  3.43 0.1

(W  L )

Ligne micro ruban à 141.4 0.018  13.21 0.1

(W  L )

Ligne micro ruban à 100 0.09  18.715 0.1

(W  L )

Tableau IV.6: Paramètres de l’antenne.

Les résultats de la simulation sont représentés sur les figures suivantes :
Ansoft LLC S11 en dB HFSSDesign1 ANSOFT
-2.00
Curve Info
dB(S(1,1))
Setup1 : Sw eep1
-4.00

-6.00

-8.00
dB(S(1,1))

-10.00

-12.00

-14.00

-16.00

-18.00
5.50 5.75 6.00 6.25 6.50 6.75 7.00
Freq [GHz]

figure IV.32: variation de coefficient de réflexion en fonction de la fréquence

92
 Chapitre IV Modélisation et simulation d’un réseau d’antennes patch
en utilisant HFSS
-La fréquence de résonance f r  5.87GHz

-Le coefficient de réflexion S11  16.92dB

-La bande passante BP  158MHz  BP  %   2.69% 

figure IV.33: représentation 3-D du Gain total en dB

Ansoft LLC Plan theta HFSSDesign1 ANSOFT

0 Curve Inf o
dB(GainTotal)
-30 30 Setup1 : Sw eep1
Freq='5.873786408GHz' Phi='0deg'
-1.00

-7.00
-60 60
-13.00

-19.00

-90 90

-120 120

-150 150

-180

figure IV.34 : digramme de rayonnement de Gain total en dB (plan  )

Le gain est maximal pour   0 et égal à 1.36dB et il est minimal pour   55 et égal à -22.43dB
Ansoft LLC Plan phi HFSSDesign1 ANSOFT

0 Curve Info
dB(GainTotal)
-30 30 Setup1 : Sw eep1
Freq='5.873786408GHz' Theta='90deg'
-10.00

-15.00
-60 60
-20.00

-25.00

-90 90

-120 120

-150 150

-180

figure IV.35:.digramme de rayonnement en Gain en dB (plan  )

Le gain est maximal pour   90 et égal à -5.29dB et il est minimal pour   35 et égal à -28.05dB

93
 Chapitre IV Modélisation et simulation d’un réseau d’antennes patch
en utilisant HFSS
IV.6- Modélisation d’un réseau bidimensionnel par HFSS :
IV.6.1- Réseau bidimensionnel à quatre patchs :
Le réseau bidimensionnel peut être obtenu à partir de réseaux linéaires à alimentation parallèle. Les 4
patchs identiques sont alimentés en phase et l’inclinaison du faisceau rayonné ne peut être obtenue que dans
un seul plan parallèle au réseau. Deux points homologues de deux éléments rayonnants successifs sont
séparés par une distance d  13.12mm .

figure IV.36 : Réseau bidimensionnel à quatre patchs

Ansoft LLC S11 en dB HFSSDesign1 ANSOFT
0.00
Curve Info
dB(S(1,1))
Setup1 : Sw eep1

-5.00

-10.00
dB(S(1,1))

-15.00

-20.00

-25.00
1.00 1.20 1.40 1.60 1.80 2.00
Freq [GHz]

figure IV.37: variation de coefficient de réflexion en fonction de la fréquence

-La fréquence de résonance f r  1.29GHz

-Le coefficient de réflexion S11  24.90dB

-La bande passante BP  112 MHz  BP  %   8.68%  .

La représentation 3D en coordonnées polaire du gain en dB nous montre bien son variation en fonction
de deux angles  et 

94
 Chapitre IV Modélisation et simulation d’un réseau d’antennes patch
en utilisant HFSS

figure IV.38: représentation 3-D du Gain total en dB

Ansoft LLC Plan theta HFSSDesign1 ANSOFT

0 Curve Info
dB(GainTotal)
-30 30 Setup1 : Sw eep1
Freq='1.291262136GHz' Phi='0deg'
-16.76

-16.92
-60 60
-17.08

-17.24

-90 90

-120 120

-150 150

-180

figure IV.39 : digramme de rayonnement de Gain total en dB (plan  )

Le gain est maximal pour   5 et égal à 1.36dB et il est minimal pour   90 et égal à -17.3dB
Ansoft LLC Plan phi HFSSDesign1 ANSOFT

0 Curve Info
dB(GainTotal)
-30 30 Setup1 : Sw eep1
Freq='1.291262136GHz' Theta='90deg'
-19.00

-23.00
-60 60
-27.00

-31.00

-90 90

-120 120

-150 150

-180

figure IV.40:.digramme de rayonnement en Gain en dB (plan  )

Le gain est maximal pour   175 et égal à -5.29dB et il est minimal pour   265 et égal à -32.77dB

IV.6.1- Réseau bidimensionnel à 8 patchs :

Le réseau bidimensionnel à 8 patchs de 9GHz peut être obtenu à partir de réseaux linéaires à
alimentation parallèle.

95
 Chapitre IV Modélisation et simulation d’un réseau d’antennes patch
en utilisant HFSS

figure IV.41: Réseau bidimensionnel à 8 patchs

Le résultat de la simulation est représenté sur les figures suivantes :
Ansoft LLC S11 en dB HFSSDesign1 ANSOFT
0.00
Curve Info
dB(S(1,1))
Setup1 : Sw eep1

-10.00

-20.00
dB(S(1,1))

-30.00

-40.00

-50.00
1.00 1.20 1.40 1.60 1.80 2.00
Freq [GHz]

figure IV.42: variation de coefficient de réflexion en fonction de la fréquence

-La fréquence de résonance f r  1.31GHz

-Le coefficient de réflexion S11  48.36dB

-La bande passante BP  140 MHz  BP  %   10.68%  .

figure IV.43: représentation 3-D du Gain total en dB

96
 Chapitre IV Modélisation et simulation d’un réseau d’antennes patch
en utilisant HFSS
Ansoft LLC Plan theta HFSSDesign1 ANSOFT

0 Curve Info
dB(GainTotal)
-30 30 Setup1 : Sw eep1
Freq='1.310679612GHz' Phi='0deg'
-10.52

-10.84
-60 60
-11.16

-11.48

-90 90

-120 120

-150 150

-180

figure IV.44 : digramme de rayonnement de Gain total en dB (plan  )

Le gain est maximal pour   0 et égal à -10.35dB et il est minimal pour   90 et égal à -11.7dB.
Ansoft LLC Plan phi HFSSDesign1 ANSOFT

0 Curve Inf o
dB(GainTotal)
-30 30 Setup1 : Sw eep1
Freq='1.310679612GHz' Theta='90deg'
-15.00

-20.00
-60 60
-25.00

-30.00

-90 90

-120 120

-150 150

-180

figure IV.45:.digramme de rayonnement en Gain en dB (plan  )

Le gain est maximal pour   175 et égal à -11.66dB et il est minimal pour   270 et égal à -30dB

IV.8-Conclusion :
Dans ce chapitre, on a basé sur deux paramètres principaux tels que : le type d’alimentation (série ou
parallèle) et le type de réseau (linéaire ou bidimensionnel) pour avoir l’influence de ces deux paramètres sur
les caractéristiques de rayonnement d’un réseau d’antennes à savoir la fréquence de résonance le cœfficient
de réflexion, la bande passante et le gain.
Pour le cas d’alimentation en série de deux patchs, on a remarqué que le couplage dans le plan H ayant
une bande passante plus large que le couplage dans plan E, mais un gain faible.
Pour le cas d’alimentation en parallèle de deux patchs, on a trouvé que le couplage ayant une bande
passante large mais a une grande atténuation et faible gain. Par contre le couplage entre quatre patchs ayant
une bande passante d’environ 3% et fort gain et faible atténuation.
Donc, pour une large bande passante et meilleur gain en utilise un réseau à plusieurs patchs alimenté en
parallèle.
Par contre le réseau bidimensionnel ayant une large bande passante quand l’élément de réseau
augmente, et aussi une faible atténuation et fort gain.

97
 Chapitre IV Modélisation et simulation d’un réseau d’antennes patch
en utilisant HFSS
Bibliographie
[IV.1] Abdellatif Hafiane « Etude du couplage d’antennes imprimées par la méthode des différences finies et
extraction du modèle électrique » Mémoire pour le garde de docteur de l’université de Marne la
Vallee.Mars2003.
[IV.2] Chouti Leila « Contribution à l’étude d’antennes imprimés rectangulaires double bande et multi
bandes tenant compte de l’effet de couplage » Mémoire pour obtenir le diplôme de magister
janvier 2009.
[IV.3] Tomader Mazri et Fatima Riouch « Improved performance of a microstrip antenna array using a tree
structure patch fed by electromagnetic coupling» ICONS 2011: The Sixth International Conference
on Systems.

[IV.4] Benachour Leila « Modélisation des structures planaires de formes régulières implantées sur substrat
d’une anisotropie tensorielle’’ Mémoire présenté pour l’obtention du diplôme de magistère en
électronique 2008.

[IV.5] Paul [Link] « Micro ondes [Link] passifs, propagation, antennes. Cours et exercice »

[IV.6] « Handbook of microstrip antennas, volume 2 » Jim R. James, institution of electrical engineers
[IV.7] Barkat Ouarda« Etude et synthèses des caractéristique de réseaux d’antennes imprimées
supraconductrices dans la bande millimétrique » Thése de doctorat de l’université Mentouri
Constantine décembre 2009
[IV.8] « Antenna Theory Analysis and Design» third edition .Constantine, A Balanis Ajohn Wiley a
[Link].
[IV.9] « Application of antenna arrays to mobile communications, part I: performance improvement,
feasibility, and system considerations » IEEE, Vol. 85, No. 7, July 1997

[IV.10] David Bidou «Contribution a l’étude d’antenne fort gain .application dans le domaine
millimétrique» Thése pour obtenir le garde de docteur de l’université de l’Imoges Mai 2003.
[IV.11] Chouti Leila « Contribution à l’étude d’antennes imprimés rectangulaires double bande et multi
bandes tenant compte de l’effet de couplage » Mémoire pour obtenir le diplôme de magister
janvier 2009.

[IV.12] Hatem Rmili « Étude, réalisation et caractérisation d’une antenne plaquée en polyanilline
fonctionnant à 10 GHz » Thèse pour obtenir le garde de docteur de l’université de Bordeaux I,
novembre 2004.
[IV.13] [Link]/base.../antennes.../antennes-imprimées-e3310/
[IV.14] Daniel Courivaud – Sigtel - Groupe ESIEE « Micro-ondes » Paris 2002.

98
Chapitre V Réalisations

V.1- Introduction :
Dans cette partie, nous avons réalisé des modèles d’antennes imprimées à élément rectangulaire. La
réalisation est faite au niveau du laboratoire de Télécommunications de notre département Electronique et
Télécommunications en basant sur la technologie des circuits imprimés. La caractérisation des modèles
réalisés est effectuée au niveau du laboratoire cité précédemment où on a utilisé le banc d’essai d’antenne
qui est géré par le logiciel Labsoft.

V.2- Les étapes de réalisation :

En tenant compte des étapes de conception citées dans le chapitre précédent, nous avons réalisé trois
types d’antennes imprimées, en utilisant comme substrat l ‘Epoxy ‘ ayant une permittivité relative  r  4.8
et d’épaisseur 1.6 mm
La réalisation d’antennes imprimée s’effectuée selon les étapes suivantes :
-1ère étape : On a tracé la forme d’antenne simulée .
- 2ème étape : Insolation : On couvre les traçages par un encre noir, puis on les attaque par l’ultra-violet .
- 3ème étape : Révélation : On met la plaquette dans une cuvette contenant NaOH .
- 4ème étape : Gravure : la gravure est s’effectue dans une machine contenant PerChlorure de Fer, elle va
attaquer la partie non protégée de l’ultra-violet et elle l’enlève .

V.3- Exemples de réalisation et simulation :

Après modélisation et simulation, les antennes sont réalisées suivant les procédures de conception d’une
antenne imprimée rectangulaire. Avec des petites approximations on a obtenu la fréquence de résonnance à
8.5GHz pour les deux premiers exemples et 5.66GHz pour le troisième exemple. Les modèles réalisés sont
caractérisés au niveau du laboratoire de Télécommunications.

 Premier exemple : Patch alimenté par un transformateur d’impédance :

figure V.1 : patch alimenté par ligne avec transformateur résonant à 8.5GHz

figure V.2 : Le modèle de patch alimenté par ligne avec transformateur par le HFSS
99
Chapitre V Réalisations
Le tableau V.1 et V.2 qui récapitulent les propriétés de l’antenne simulée.

Matériaux Paramètres

Plan de masse : cuivre 5.8 107 S / m

Substrat : Bakelite (Epoxy)  r  4.8 et Tan  0.002

Patch : cuivre 5.8 107 S / m

Tableau V.1 : Propriétés des matériaux

Eléments aire (mm  mm) épaisseur (mm)

Substrat 20  20 h=1.6

Plan de masse 20  20 0.05

Patch (W  L) 9.36  7.1 0.05

Ligne micro ruban (W l  Ll ) 2.86  7.39 0.05

Transformateur d’impédance (W t  Lt ) 0.935  4.62 0.05

Tableau V.2: Paramètres de l’antenne.

Ensuite, nous avons simulé l’antenne avec les dimensions obtenues, avec des petites approximations
nous avons obtenu les résultats suivants .
Ansoft LLC S11 en dB HFSSDesign1 ANSOFT
-5.00

Curve Info
dB(S(1,1))
-10.00 Setup1 : Sw eep1

-15.00
dB(S(1,1))

-20.00

-25.00

-30.00

-35.00
7.50 8.00 8.50 9.00 9.50 10.00
Freq [GHz]

figure V.3 : Variation du coefficient de réflexion en fonction de la fréquence.

- La fréquence de résonance f r  8.5GHz

-Le coefficient de réflexion S11  33.65dB

-La bande passante BP  959MHz ( BP  %   11.28% ).

 Plan horizontal :

100
Chapitre V Réalisations
Ansoft LLC Plan theta HFSSDesign1 ANSOFT

0 Curve Info
dB(GainTotal)
-30 30 Setup1 : Sw eep1
Freq='8.503236246GHz' Phi='0deg'
-2.00

-4.00
-60 60
-6.00

-8.00

-90 90

-120 120

-150 150

-180

figure V.4: Digramme de rayonnement de Gain total en dB (plan  )

Le gain est maximal pour   0 et égal à -0.11dB et il est minimal pour   90 et égal à -9.11dB.

figure V.5: Digramme de rayonnement mesuré de Gain total (plan  )

 Plan verticale:
Ansoft LLC Plan phi HFSSDesign1 ANSOFT

0 Curve Inf o
dB(GainTotal)
-30 30 Setup1 : Sw eep1
Freq='8.503236246GHz' Theta='90deg'
-6.80

-7.60
-60 60
-8.40

-9.20

-90 90

-120 120

-150 150

-180

figure V.6: Digramme de rayonnement en Gain en dB (plan  )

Le gain est maximal pour   235 et égal à -6.17dB et il est minimal pour   140 et égal à -9.97dB

figure V.7: Digramme de rayonnement mesuré de Gain total (plan  )

101
Chapitre V Réalisations

figure V.8 : Représentation 3-D du Gain total en dB

La représentation 3D en coordonnées polaire du gain en dB nous montre bien son variation en fonction
de deux angles  et  .

 Deuxième exemple: Patch alimenté par ligne microruban directe :

figure V.9 : patch alimenté par ligne à 8.5GHz

figure V.10 : Le modèle de patch alimenté par ligne a réalisé

Les tableaux V. 3 et V.4 récapitulent les propriétés de l’antenne simulée

Matériaux Paramètres

Plan de masse : cuivre 5.8 107 S / m

Substrat : Bakelite (Epoxy)  r  4.8 et Tan  0.002

Patch : cuivre 5.8 107 S / m

Tableau V. 3 : propriétés des matériaux

102
Chapitre V Réalisations
Eléments aire (mm  mm) épaisseur (mm)

Substrat 20  20 h=1.6

Plan de masse 20  20 0.05

Patch (W  L) 10.36  4.36 0.05

Ligne micro ruban (W l  Ll ) 2.86  7.39 0.05

Tableau V. 4 : Paramètres de l’antenne

Ensuite, nous avons simulé l’antenne avec les dimensions obtenues, avec des petites approximations nous
avons obtenu les résultats suivants .
Ansoft LLC S11 en dB HFSSDesign1 ANSOFT
-7.50
Curve Info
dB(S(1,1))
Setup1 : Sw eep1
-10.00

-12.50

-15.00
dB(S(1,1))

-17.50

-20.00

-22.50

-25.00

-27.50
7.50 8.00 8.50 9.00 9.50 10.00
Freq [GHz]

figure V.11 : Variation du coefficient de réflexion en fonction de la fréquence.

- La fréquence de résonance f r  8.5GHz

-Le coefficient de réflexion S11  25.92dB

-La bande passante BP  1928MHz ( BP  %   22.68% ).

 Plan horizontal :
Ansoft LLC Plan theta HFSSDesign1 ANSOFT

0 Curve Inf o
dB(GainTotal)
-30 30 Setup1 : Sw eep1
Freq='8.503236246GHz' Phi='0deg'
-12.00

-14.00
-60 60
-16.00

-18.00

-90 90

-120 120

-150 150

-180

figure V.12: Digramme de rayonnement de Gain total en dB (plan  )

Le gain est maximal pour   0 et égal à -11.17dB et il est minimal pour   90 et égal à -19.14dB.

103
Chapitre V Réalisations

figure V.13: Digramme de rayonnement mesurée de Gain total (plan  )

 Plan verticale:
Ansoft LLC Plan phi HFSSDesign1 ANSOFT

0 Curve Info
dB(GainTotal)
-30 30 Setup1 : Sw eep1
Freq='8.503236246GHz' Theta='90deg'
-12.40

-14.80
-60 60
-17.20

-19.60

-90 90

-120 120

-150 150

-180

figure V.14: Digramme de rayonnement en Gain en dB (plan  )

Le gain est maximal pour   70 et égal à -11.79dB et il est minimal pour   345 et égal à -21.23dB

figure V.15: Digramme de rayonnement mesurée de Gain total (plan  )

104
Chapitre V Réalisations

figure V.16 : Représentation 3-D du Gain total en dB

 Toisième exemple : Réseau à quatre patchs alimenté en paralléle :

figure V.17 : Réseau à quatre patchs alimenté en paralléle

figure V.18 : Le modèle réalisé

Les tableaux V.5 et V.6 récapitulent les propriétés de l’antenne simulée

Matériaux Paramètres

Plan de masse : cuivre 5.8 107 S / m

Substrat : Bakelite (Epoxy)  r  4.8 et Tan  0.002

Patch : cuivre 5.8 107 S / m

Tableau V.5 : propriétés des matériaux

105
Chapitre V Réalisations
Eléments Aire (mm  mm) Epaisseur (mm)

Substrat 50  40 h=1.6

Plan de masse 50  40 0.05

Patch (W  L ) 10.36  7.39 0.05

Ligne micro ruban à 50 (W  L ) 2.86  7.39 0.05

Ligne micro ruban à 70.7 1.49  4.77 0.05

(W  L )

Ligne micro ruban à 100 0.637  4.63 0.05

(W  L )

Ligne micro ruban à 141.4 0.196  13.757 0.05

(W  L )

Ligne micro ruban à 100 0.637  18.221 0.05

(W  L )

Tableau V.6: Paramètres de l’antenne.

Ensuite, nous avons simulé l’antenne avec les dimensions obtenues, avec des petites approximations nous
avons obtenu les résultats suivants .
Ansoft LLC S11 en dB HFSSDesign1 ANSOFT
-5.00
Curve Info
dB(S(1,1))
Setup1 : Sw eep1
-7.50

-10.00

-12.50
dB(S(1,1))

-15.00

-17.50

-20.00

-22.50

-25.00
5.00 5.25 5.50 5.75 6.00 6.25 6.50 6.75 7.00
Freq [GHz]

figure V.19 : Variation du coefficient de réflexion en fonction de la fréquence.

- La fréquence de résonance f r  5.66GHz

-Le coefficient de réflexion S11  24.67 dB

-La bande passante BP  516MHz ( BP  %   9.11% ).

106
Chapitre V Réalisations
Ansoft LLC Plan theta HFSSDesign1 ANSOFT

0 Curve Inf o
dB(GainTotal)
-30 30 Setup1 : Sw eep1
Freq='5.666666667GHz' Phi='0deg'
0.00

-5.00
-60 60
-10.00

-15.00

-90 90

-120 120

-150 150

-180

figure V.20: Digramme de rayonnement de Gain total en dB (plan  )

Le gain est maximal pour   0 et égal à 2.24dB et il est minimal pour   65 et égal à -19.19dB.
Ansoft LLC Plan phi HFSSDesign1 ANSOFT

0 Curve Info
dB(GainTotal)
-30 30 Setup1 : Sw eep1
Freq='5.666666667GHz' Theta='90deg'
-5.00

-10.00
-60 60
-15.00

-20.00

-90 90

-120 120

-150 150

-180

figure V.21: Digramme de rayonnement en Gain en dB (plan  )

Le gain est maximal pour   90 et égal à -3.15dB et il est minimal pour   160 et égal à -23.52dB

figure V.22 : Représentation 3-D du Gain total en dB

V.5-conclusion :
Après avoir simulé les différents types d’antenne patch rectangulaire disposé en réseau ou comme un
seul élément, on a passé à l’étape de réalisation et caractérisation. Les antennes patchs réalisées sont les plus
utilisées dans le domaine de télécommunications sans fil telles que l’antenne patch rectangulaire à
alimentation directe par ligne micro ruban ou à travers un transformateur d’impédance. .
A travers cette partie de réalisation, nous avons vu les outils de la technologie des circuits imprimés.
Enfin, nous pouvons dire que les resultats trouvés par la réalisation d’une antenne patch qui rayonne à
8.5GHz se conforme bien avec les resultats trouvés par simulation. Les petites différences entre les
prévisions théoriques et les mesures trouvées sont dues essentiellement à l’incertitude de dimensionnement
des antennes imprimées lors de la réalisation.

107
 Conclusion Générale

Conclusion Générale
Ces dernières années, l'augmentation du nombre de standards de communications (GSM 900, DCS
1800, UMTS, WLAN, BLUETOOTH...), a fortement encouragé les chercheurs à trouver des modèles
d'antennes patch de plus en plus complexes faisant généralement appel à des structures plus compliqué
(court-circuitées, éléments avec fentes rayonnantes) qui permettent d’obtenir les caractéristiques de
rayonnement désirables, une large bande passante et une antenne de plus petite taille possible.

En premier lieu nous avons commencé par l’essentiel de la théorie des antennes imprimées, ainsi que
les différentes méthodes d’analyse électromagnétiques existantes. Ensuite nous présentons la modélisation et
la simulation d’une antenne imprimée rectangulaire à alimentation directe par ligne microruban en basant
sur le logiciel HFSS (High Frequency Structure Simulator). En plus, nous avons évoqué l’influence de
quelques paramètres caractéristiques sur les performances de l’antenne. On à aussi simulé cette même
antenne et avec le même logiciel en changeant uniquement la méthode d’alimentation (câble coaxiale, ligne
micro ruban avec un transformateur et avec une ligne enterré).
Le résultat obtenu montre que l’alimentation par ligne microruban directe ayant une large bande
passante par apport aux autres types d’alimentation.

En fin, une étude d’un réseau d’antennes imprimées est basée sur le couplage entre les éléments
rayonnants. Dans ce contexte, nous avons modélisé et simulé différents type de réseau d’antennes (réseau : à
deux éléments, à quatre éléments et huit éléments) en basant sur les deux types d’alimentations série et
parallèle
Pour atteindre notre objectif, des simulations numériques avec le logiciel HFSS sont réalisées. Ce
logiciel nous a permis de nous rapprocher du meilleur cas pratique. A travers ses simulations nous avons
étudié le comportement physique, électrique et de rayonnement de différentes antennes plaquées (la
fréquence de résonance, la largeur de bande, diagramme de rayonnement du gain,…) en fonction des
paramètres qui intervient dans la conception (propriétés des matériaux, paramètres de l’antenne,…).

Les principaux résultats obtenus sont :

 Une mince couche du diélectrique à permittivité donnée, cause une sévère dégradation dans la bande
passante des antennes microbandes.

 La variation des paramètres géométriques de l’antenne (permittivité et hauteur du substrat, largeur,

longueur et la conductivité du patch), influe sur la fréquence de résonance, le coefficient de réflexion,
la bande passante et le gain.

 L’alimentation par ligne microruban directe ayant une large bande passante et fort gain est simple à
réaliser par apport aux autres types d’alimentation.

 Le couplage dans le plan E est plus élevé que dans le plan H.

 Pour une large bande passante et meilleur gain en utilisant un réseau à plusieurs patchs alimentés en
parallèle.

 Le réseau bidimensionnel ayant une large bande passante lorsque le nombre d’éléments rayonnants
augmente, et aussi amène à une faible atténuation et fort gain.

 D’après une comparaison entre les antennes simulées et les antennes réalisées, nous pouvons dire que
les résultats pratiques sont très proches aux résultats simulés.

108
 Annexes
A.1- Calcul du champ électrique à l’intérieur du cavité : [A.1]
Le champ électrique à l’intérieur du cavité ne peut être orienté que suivant oz , et doit être constant sur
une ligne de champ entre les deux plans parfaitement conducteur (car les lignes de champ sont parallèles).
Nous désignerons donc le champ électrique à l’intérieur du cavité, porté par la direction oz par :
 
E  E z  x, y  U z (1)
De l’équation de MAXWELL :
  
  E   j B (2)
Dans laquelle B représente le champ magnétique interne, on déduit :
  
 
 x   0    j Bx 
      

 y    0    j B y  (3)
   Ez    j B 
   z 

 
 z 
Soit encore :
 1  Ez  Ez    
B  U  U y  B U  B U (4)
j  y
x x x y y
x 
On note qu’il n’existe pas de composante de champ magnétique suivant oz , et donc que les lignes de
champ magnétiques sont parallèles au plan z = 0.
L’équation de propagation en champ électrique est donnée par :
 2  
 E  k2 E  0 (5)
Elle sera simplifie car il n’existe qu’une seule composante de champ suivant oz et cette composante n’a
pas de dépendance en z :
 2 E z  x, y   2 E z  x, y 
  k 2 Ez  x, y   0 (6)
x 2
y 2

La recherche d’une solution à variables séparées conduit à poser :

E z  x, y   f  x  .g  y  (7)
k 2  kx2  k y 2 (8)
L’équation (5) se ramène à la résolution de deux équations indépendantes :
2 f  x
 kx 2 f  x   0 (9)
x 2

 g  y
2

 ky2 g  y   0 (10)
y 2
Dont les solutions évidentes sont :
f  x   [Link]  k x .x   [Link]  k x .x  (11)
g  x   [Link]  k y . y   [Link]  k y . y  (12)
Où A, B, C, D, est 4 constantes qui dépendent des conditions aux limites.
Le champ électrique Ez se présente donc sous la forme :
Ez  x, y    [Link]  k x .x   [Link]  k x .x   . [Link]  k y . y   [Link]  k y . y   (13)
Les conditions aux limites sur les surfaces parfaitement conductrices ont déjà été utilisées pour préciser
la direction du champ électrique dans la cavité.
Les conditions aux limites sur les murs magnétiques parfaits imposent un champ magnétique tangentiel
nul sur ces parois, soit :
Bx  0 en y  0 et y  b (14)

109
 Annexes
By  0 en x  0 et x  a (15)
De (4) et (14) on déduit:
Ez  x, y 
Bx  0 0 (16)
y
C.k y .sin  k y . y   D.k y .cos  k y . y   0 en y  0 et y  b (17)
en y  0  D  0 (18)
en y  b  k y .b  n.  k y  n. / b (19)
n étant un entier à priori quelconque.
De (4) et (15) on déduit :
Ez  x, y 
By  0 0 (20)
x
 A.K x .sin  K x .x   B.K x .cos  K x .x   0 en x  0 et x  b (21)
en x  0  B  0 (22)
en x  b  k x .a  m.  k x  m. / a (23)
m étant un entier à priori quelconque.
En utilisant les expressions de k x et k y obtenus ci-dessus, et en tenant compte de la nullité des
constantes B et D, le champ électrique (13) s’écrit sous la forme suivante, à une constante multiplicative près
qui est posée arbitrairement égale à A : [A.1] [A.2]:
 m   n 
E z  x , y   A cos  x  cos  y (24)
 a   b 
 la fréquence de résonance : [A.1] [A.2]
La relation (8) permet de calculer la fréquence de résonance de chacun de ces modes. Réécrite avec les
expressions de k x et k y données en (19) et (24), elle permet de prévoir la fréquence de fonctionnement de
l’antenne :
2 2
 m   n 
k 
2
   (25)
 a   b 
k 2   2 . .  4. 2 . f 2 . . (26)
D’où on déduit:
2 2
1  m    n 
f      (27)
2   a   b 

A.2- Calcul du champ rayonné pour (  r  1 ) : [A.1] [A.2]

Dans un premier temps, nous calculons le champ lointain rayonné en espace libre au point P repéré
par ses coordonnées polaires  r , ,  

    
jk
  
E  p    R  J  x, y U U e jKOM .ds
4 s
(28)

Les courants à la surface du patch sont donnés:

110
 Annexes
     
J x , y   j sin  x U x (29)
w a a 
En substituant les courants (29) dans l’intégrale (28), et en regroupant les termes constants sous une
constante C, on obtient :
           

E  p   C  sin  x   e j k OM .ds. U x  U  U
 a 
 (30)
s 
Le vecteur de polarisation de l’onde rayonnée est donné par le double produit vectoriel :
 sin  cos   1   1   0 
   
        
U x  U  U   cos  sin     0     0     cos  sin    (31)
  sin    0    0   sin  

M est un point courant de la surface, de coordonnées  x, y, h  où est situé le courant. Le produit scalaire
de (30) s’évalue de la manière suivante :
 k sin  cos   x 
    
k .OM   k sin  sin   y   k  x sin  cos   y sin  sin   h cos   (32)
 k cos   h 
  

En utilisant ces deux derniers résultats, l’intégrale (I .30) s’écrit :

 b a
   j k x sin cos  
E  p  C e j k h cos
0 e j k y sin  sin 
dy 0  a 
sin  x  e  dx  cos  cos  U   sin  U   (33)

Afin d’alléger l’écriture, on peut poser :

V  k sin  sin  (34)

W  k sin  cos  (35)
Les intégrales présentes dans (33) deviennent alors:

b y b
e jV y e jV b  1
e dy  
jV y
Pour V différent de 0 (36)
0
jV y 0
jV
b

e dy  b Pour V  0
jV y
(37)
0

   /a
a

e
jWx
sin  x  dx   e jWa  1 Pour W différent de k   / a (38)
 / a   W
2
0 a  2

 
a
a
e sin  x  dx  j Pour W  k   / a
jWx
(39)
0 a  2
D’où l’écriture globale du champ rayonné en espace libre par les courants présents sur le patch :
 jWa
 jVb
  e  1  
E  p   Ce j k h cos e 1 a
jV    2
 cos cos  U   sin  U   (40)
  W
2

a
Cette écriture peut encore se simplifier, en rappelant que k  2 /  , et que sur le mode fondamental, et
pour un diélectrique tel que  r  1 , la dimension de l’antenne a est égal à  / 2 .On en déduit que dans ce cas
particulier k   / a .
Nous devons maintenant prendre en compte l’influence du plan de masse en considérant le
rayonnement des courants images qui sont opposés au courants du patch, et situés à la cote –h, d’où leur
expression déduite de (40):
111
 Annexes
  1 k  e  1  
jVb jWa

E  p image  Ce j k h cos

jV
e
k W
2 2  
 cos  cos  U  sin  U  (41)

En sommant (I .40) et (I .41) on obtient le champ total rayonné par l’antenne sur son plan de masse :
 
sin khcos ejVb 1 k e 1
jWa
  
E ptotal C
2j jV k2 W2

cos cos U sinU  (42)

A.3- Rayonnement par application du théorème d’équivalence : [A.1]

L’application du théorème d’équivalence se fait en respectant les étapes suivantes :

 Etape 1 : définition de la surface de HUYGENS

La surface de HUYGENS est constituée par un plan théoriquement infini, situé juste au-dessus du patch
et du diélectrique.

 Etape 2 : Définition des courants équivalents

Les courants électriques et magnétiques équivalents sont définis à partir des champs présents sur la
surface de HUYGENS par les relations :
     
JE  n H J M  n  E (43)
Sur le patch, le champ électrique tangentiel étant nul, il n’existe que des courants électriques
équivalents.
D’après le théorème d’équivalence, le rayonnement de ces courants équivalents conduit à un champ
électromagnétique inchangé au-dessus de la surface de HUYGENS, et à un champ électromagnétique nul en
dessous de cette surface.

 Etape 3 : On remplit la région où le champ est nul par un conducteur électrique parfait

Puisque le champ électromagnétique est nul dans le demi espace inférieur, on peut, à condition de
prendre en compte les modifications que cela implique à remplir cet espace par un matériau parfaitement
conducteur.

 Etape 4 : On applique le théorème des images par rapport à la surface de HUYGENS

Chaque courant électrique tangentiel au-dessus de la surface a une image opposée en dessous : on en
déduit qu’à la limite, les deux courants s’annulent lorsqu’ils sont tous les deux au niveau de la surface. Par
conséquent, le rayonnement des courants électriques est nul.
Chaque courant magnétique tangentiel au-dessus de la surface a une image de même sens en dessous :
on en déduit qu’à la limite, les deux courants s’ajoutent lorsqu’ils sont tous les deux au niveau de la surface.
112
 Annexes
Par conséquent, le rayonnement des courants magnétiques est doublé par la présence du matériau
parfaitement conducteur

 Etape 5 : Simplification
Les courants magnétiques équivalents n’ont une valeur significative qu’au voisinage immédiat du patch
: on montre que leur décroissance est exponentielle au fur et à mesure que l’on s’éloigne du bord. On peut
donc considérer avec une bonne approximation que le champ rayonné par l’antenne peut être déduit des
courants magnétiques localisés autour du patch :

Moyennant quelques approximations, l’évaluation des courants magnétiques équivalents peut se faire
de la manière suivante (figure -1-) : sur les côtés du patch (la figure représente un côté parallèle à oy ), on
définit un parcours (C) = MNOP :

figure -1- : coupe dans un plan xoz

La relation d’Ampère sur un contour fermé (C) = MNOP s’écrirait en termes de courants électriques :
 
 H .d I  I e
MNOP
(44)

Généralisée aux courants magnétiques, elle devient :

 
 E.d I  I m
MNOP
(45)

La circulation du vecteur champ électrique peut se décomposer :

 
 E.d I   E x dx   E z dz   E x dx   E z dz  I m
i e e e

MNOP MN NO OP PM (46)
1  2  3  4
où les indices i et e désignent respectivement les champs intérieurs et extérieurs à la cavité.

Si on suppose le champ électrique confiné dans la cavité, les intégrales (1), (2), (3), sont nulles.
D’où l’écriture vectorielle de la densité linéique de courant magnétique qui entoure le patch, sur la portion
AD ou BC :
 
J M   Ez i U y (47)
On peut représenter une ligne de densité courant magnétique qui entoure le patch (figure -2-)

113
 Annexes

figure-2- : Courants magnétiques équivalents qui entourent le patch rectangulaire


En posant AB = DC = a et AD = BC = b les courants magnétiques J M sont déduits du champ
interne Ez , et s’expriment de la manière suivante, sur le mode fondamental : [A.1] [A.2]
  x 
-sur la portion AB: J M  E 0 cos( )U x (48)
a
 
- Sur la portion BC : J M   E 0U y (49)
  x 
- Sur la portion DC : J M   E 0 cos( )U x (50)
a
 
- Sur la portion AD : J M   E 0U y (51)
Le calcul du champ rayonné va s’effectuer à l’aide de l’intégrale de rayonnement généralisée au courant
magnétique :

       

E  p 
jK
4
 
  r   J M  U e j K OM dI  C  JM  U e j K
I I
  (52)
 
e j K OM  e  
j Wx Vy

Dans la suite du calcul, nous poserons la constante C  j k   r  / 4

Nous posons également dans un souci de simplification d’écriture :
W  k sin  cos  et V  k sin  sin  (53)
D’où l’expression déjà rencontrée de l’exponentielle présente dans (60) :
 
e j K OM  e 
j W x V y 
(54)
Le calcul du champ total rayonné s’obtient en évaluant les contributions de chaque ligne de courant
magnétique équivalent.

 Champ lointain rayonné par les courants magnétiques sur la portion AD :

 y b   j k .OM

E AD  P    E0C  
U y U e
y 0

dy (55)
Soit, en tenant compte du fait que x = 0 entre A et D :
  
 
y b
E AD  P    E0C  e jVy dy cos U   cos  sin U  (56)
y 0

 e jVb  1  
E AD  P    E0C
jV

cos U   cos  sin U   (57)

Vb
 jVb sin  
E AD  P   bE0Ce 2
Vb 
2 cos U   cos  sin U 
 (58)

2
Champ lointain rayonné par les courants magnétiques sur la portion BC : [A.1]
Par rapport au cas précédent, les courants sont de même sens. Puisque le point M se déplace sur la
 
j Wx Vy 
droite x = a, on a : e j K OM  e
114
 Annexes
et donc :
Vb
 jVb sin  
E BC  P   bE0C e jWa e 2 
Vb
2 cos U   cos  sin U 
 (59)

2
Champ lointain rayonné par les portions AD et BC :
Il s’obtient en sommant les contributions d’AD (58) et BC (59) :
Vb
 jVb sin  
E AD , BC  P   bE0C 1  e jWa  e 2
Vb 
2 cos U   cos  sin U 
 (60)

2
Vb
 j Wa Vb  sin  
E AD , BC  P   2bE0C e 2 cos 
 Wa 

 2  Vb
 2 cos U   cos  sin U 
 (61)

2
Champ lointain rayonné par la portion AB :
Sur la portion AB, le courant est donné par l’expression (48). D’après l’intégrale de rayonnement
généralisée au courant magnétique (52) le champ lointain s’obtient en calculant :
   x    j  k .OM 
 

 
xa
E AB  P   C  E0 cos   U x  U e dx (62)
x 0
 a 
Puisque y = 0 sur cette portion, on obtient après développement du produit vectoriel :
   x  jWx  
 
xa
E AB  P   CE0  cos   e dx sin  U   cos  cos  U  (63)
x0
 a 
Wa
 jWa cos  
E AB  P   CE0 jWe 2 2
 
2
sin  U   cos  cos  U  (64)
W  
2

a
Champ lointain rayonné par la portion DC :
Il se déduit du champ lointain rayonné par la portion AB, en notant que le sens des courants est inversé,
et que l’intégrale se calcule sur une portion de droite telle que y = b :
 jWa
cos Wa / 2   
E DC  P   CE0 e jVb jWe 2
W 2   / a 
 2
sin  U 
  cos  cos U  (65)

Champ lointain rayonné par les portions AB et DC :

Il s’obtient en sommant les contributions de AD (I .64) et BC (I .65) :
 jWa
cos Wa / 2   
E AB , DC  P   CE0  e jVb  1 jWe 2
W 2   / a 
2
sin U   cos  cos  U   (66)

 j Wa Vb 
cos Wa / 2   
Vb
E AB , DC  P   2CE0 sin We 2
2 W 2   / a 
2
sin  U  cos  cos  U   (67)

Enfin, le champ total rayonné s’obtient en sommant toutes les contributions, soit les relations (61) et
(67) .En omettant les termes constants qui sont normalisés dans le diagramme de rayonnement, on obtient :

  sinVb/2 cosWa/2   sinVb/2 cosWa/2  

E p bcosWa/2 cosWsinVb/2 2 sin  U bcosWa/2 sinWsinVb/2 cos cosU (68)
 Vb/2 W /a
2 
  Vb/2 W2 a/2 
A.4-Détermination des paramètres du patch :
On a utilisé Microstrip Patch Antenna Calculator [A.3] pour déterminer les dimensions W et L du
patch

115
 Annexes

Donc, la longueur de patch est : 5.23mm et la largeur égale à : 7.23mm

A.5-Détermination des paramètres de la ligne micro ruban :
1. Largueur de la ligne :
Lorsque l’antenne imprimée est alimentée par ligne micro ruban, le largueur de la ligne est donnée
par la Méthode de Wheeler qui permet de déterminer les caractéristiques physiques de la ligne à partir
de l’impédance caractéristique de cette ligne [A.4]:

W W 8 eA
 Si 2 :  2A (69)
h h e 2

W W  r 1  0.61  2
 Si 2 :   ln  B  1  0.39     B  1ln  2 B  1  (70)
h h  r  r  

 Zc  r  1  0.11 
 A   2   r  1   0.23  
 Z0  r  1  r 
Avec :  (71)
B   Z0
 2  r Zc


 On peut utiliser aussi Microstrip Line Calculator [A.5] pour déterminer largueur de la ligne.

2. Longueur de la ligne :
Pour déterminer la longueur de la ligne, il faut d’abord déterminer la longueur électrique de la ligne
en degrés tel que : Longueur électrique=longueur physique / Vf
 Longueur physique : [A.6]
1 1
On à utilisé une onde de fréquence 9GHz à une période N de : 0.11ns ( N   ).pendant cette
f 9.109
durée elle parcourt 0.033m dans l’air, à la vitesse de lumière qui de 3  108 m / s .

116
 Annexes
1
Dans une ligne dont le coefficient de vélocité Vf =0.32 (avec : V f  ), la vitesse de propagation est
r
d’environ : 0.96  108 m / s ce qui signifie qu’en une période de 0.11ns elle aura parcouru 0.0105m
(Lphys=0.0105m).
Donc, une ligne de longueur physique 0.0105m et dont le coefficient de vélocité est de 0.32 à une longueur
électrique de 0.032m.
 La longueur électrique en dégrée est déterminée par la méthode de conversion Degrés
électrique de/vers Mètres [A.7].

 Enfin, on à utilisé Microstrip Line Calculator pour déterminer la longueur de la ligne [A.5]

117
 Annexes

Donc, la longueur de ligne est : 4.12mm et la largeur égale à : 0.632mm

3. Longueur de transformateur d’impédance:

On à utilisé un transformateur d’impédance de longueur pour obtenir l’adaptation entre la ligne
4
d’alimentation et le résonateur (patch) à 50 .
c
f
Avec :  
e
Tel que, la formule de  e a été donnée par Hammerstad [A.4] :
1
1 1  h 2
 Pour les bandes telles que W / h  1 :  e    r  1    r  1  1  12 
2 2  W

 1
2
1 1  h  2  W
 Pour les bandes telles que W / h  1 :  e    r  1    r  1 1  12   0.04 1   
2 2  W  h 
 

118
 Annexes
Bibliographie
[A.1] Vaudon Patrick-Master de recherche télécommunication hautes fréquences et optiques « XII.-Les
antennes imprimée » IRCOM- Université de Limoges.

[A.2] Somia Boussaha et Mekhancha Sarra « Analyse du diagramme de rayonnement quasi

omnidirectionnel d’une antenne patch en utilisant le modèle de cavité » International conférence
systems and processing information May 15-17, 2011, Guelma

[A.3] [Link]
[A.4] Paul F .Combes « Micro ondes1 .lignes, guides et cavités, cours et exercices »

[A.5] [Link]
50&EL=113.47& Operation=Synthesize&Wa=&W_units_list=Wmm&La=&L_units_list=Lmm

[A.6] [Link]

[A.7] [Link]
&u=[Link]
KQ-1BgoOQ

119
 Chapitre I : Généralités sur les Antennes Imprimées
figure I.1 : Schéma de principe d’une antenne à éléments rayonnants imprimés………………………………2
figure I.2 : Divers types d’éléments rayonnants. ……………………………………………………………….3
figure I.3 : Présentation d'une antenne patch a élément rectangulaire………………………………………......4
figure I.4 : éléments rayonnants étagés………………………………………....................................................4
figure I.5 : Alimentation par ligne Micro ruban……………………………………….......................................5
figure I.6: Alimentation par une ligne microbande : a. axiale, b. décalée, c. axiale avec encoche ………….....5
figure I.7. a : Le patch alimenté par ligne micro ruban à travers un bord rayonnant…………………………...6
figure I.7.b : Le patch alimenté par ligne micro ruban à travers un bord non rayonnant…………………….....6
figure I.7.c : Le patch couplé capacitivement à ligne micro ruban……………………………………………..6
figure I.8 : Alimentation par connexion directe à une ligne coaxiale…………………………………………...7
figure I.9 : Alimentation par couplage électromagnétique à une ligne microbande à travers une fente dans le
plan de masse………………………………………………………………………………………………….....7
figure I.10 : Alimentation par couplage électromagnétique à une ligne microbande….......................................8
figure I.11 : Alimentation par guide d’ondes coplanaires…………………………………………………........8
figure I.12: Alimentation par fibre optique…………………………………………………..............................9
figure I.13 : Représentation en quadripôle…………………………………………………...............................9
figure I.14 : schéma équivalent du dispositif à l’émission…………………………………………………….10
figure I.15: Orientation du plan E et des courants surfaciques sur l’élément rayonnant
a. selon la position .b. du point d’excitation………………………………………………………………….12
figure I.16 : excitation permettant d’obtenir des polarisations orthogonales pour un élément : a-carré, b-
circulaire, c- à grilles orthogonales…………………………………………………………………………......12
figure I.17 : création d’un déphasage de 90° entre les points d’alimentation : pour des lignes dont les
longueurs différent de  / 4 d’après Zurcher et Gardiol b- par un coupleur à 3dB…………………………......13
figure I.18: Point d’alimentation pour obtenir une polarisation circulaire avec un élément a. presque carré, b.
carré avec une fente en diagonale, c. carré à coins tronqués………………………….......................................13
figure I.19: Mécanisme de rayonnement…………………………....................................................................13
figure I.20 : Lignes de champ entre les bords d’un élément imprimé et le plan de masse…………………….14
figure I.21 : Lignes de champ entourant un élément imprimé : a. composantes horizontales, b. composantes
verticales…………………………......................................................................................................................15
figure I.22 : Modélisation d’une fente excitée par un champ électrique à polarisation uniforme…………….15
   
 
figure I.23 : a-coordonnées et composantes  ,  pour les champs E  et H  rayonnés par un doublet
électrique, b- coordonnées et notations pour calculer le rayonnement d’un courant magnétique de longueur
W…………………………..................................................................................................................................16
figure I.24 : Modèle équivalent utilisé pour calculer le rayonnement d’un élément rectangulaire
a).représentation dans l’espace, b).schéma dans le plan xOy …………………………......................................17
Chapitre II : Méthodes d’Analyse des antennes imprimées
figure II.1 : Modélisation de l’antenne patch rectangulaire par deux fentes couplées………………………...22
figure II.2 : schéma équivalent de l’antenne………………………………………………………………......22
figure II.3: variation de  e en fonction de la fréquence……………………………………………………….24
figure II.4 : schéma équivalent final de l’antenne…………………………………………………………......25
figure II.5: Modèle de la cavité…………………………………………………………..................................26
figure II.6 : Représentation des murs électriques et magnétiques dans le modèle de la cavité……………......26
figure II.7 : Représentation de la cavité dans un repère (O, x, y, z) ………………………………………......27
figureII.8 : Représentation du référentiel et des coordonnées pour le calcul du champ lointain rayonné par
les courants………………………………………..............................................................................................28
figure II.9 : Diagramme de rayonnement dans le plan E en champ……………………………………….......29
figure II.10 : Diagramme de rayonnement dans le plan E en puissance……………………………………...29
figure II.11: Diagramme de rayonnement dans le plan H en champ………………………………………......30
figure II.12 : Diagramme de rayonnement dans le plan H en puissance………………………………………30
figureII.13 : Courants magnétiques équivalents qui entourent le patch rectangulaire……………………......30
figureII.14 : diagrammes de rayonnement en puissance avec une variation de  r dans le plan E…………...31
figure II.15 : Le diagramme de rayonnement dans le plan H avec variation de b/a lorsque  r  2 ………......32
figure II.16 : Le diagramme de rayonnement dans le plan H avec variation de b/a lorsque  r  5 ….............32
figure II.17 : Le diagramme de rayonnement dans le plan H avec variation de b/a lorsque  r  7 …..............33
figureII.18 : Le diagramme de rayonnement dans le plan H avec variation de b/a lorsque  r  9 …...............33
figure II.19 : Point d’évaluation du calcul de la dérivée centrée…....................................................................39
figure II.20 : Circulation du champ E autour de H….........................................................................................40
figure II.21 : Disposition des nœuds électriques et des nœuds magnétiques dans le plan xoy …......................40
figure II.22 : Circulation du Champ H autour de E…........................................................................................40
figure II.23 : Cellule de YEE…..........................................................................................................................40
figure II.24 : Extraction d’une cellule élémentaire….........................................................................................41
figure II.25 : Disposition des nœuds électriques et des nœuds magnétiques dans le plan oxy , pour le cas
d’un maillage irrégulier…...................................................................................................................................41
figure II.26: Calcul de H à l’instant ndt et calcul de E à l’instant (n+0.5) dt….................................................42
figure II.27: Maillage triangulaire adaptif en FEM…........................................................................................43
figure II.28: Les différentes ´étapes de la méthode FEM…...............................................................................44
Chapitre III : Modélisation et simulation (par HFSS) d’une antenne patch en utilisant les
différentes techniques d’alimentation
figure III.1 : Division de la structure en tétraèdres…........................................................................................49
figureIII.2 : Antenne imprimé alimenté par ligne micro ruban…......................................................................50
figure III.3 : Schéma de l’antenne, W : largeur ; L : longueur ; e : épaisseur du patch ; h : épaisseur du substrat
; Wi : largeur de la ligne d’alimentation…..........................................................................................................50
figure III.4 : Variation de coefficient de réflexion en fonction de la fréquence…............................................51
figure III.5 : Représentation de partie réelle et imaginaire du coefficient de réflexion en fonction de la
fréquence…..........................................................................................................................................................52
figure III.6 : Représentation du coefficient de réflexion sur l’abaque de Smith...............................................52
figure III.7 : Représentation de partie réelle et imaginaire de l’impédance en fonction de la fréquence..........53
figure III.8 : Représentation du ROS en dB en fonction de la fréquence.........................................................53
figure III.9 : Représentation 3-D du Gain total en dB........................................................................................53
figureIII.10 : Digramme de rayonnement de Gain total en dB (plan  ).........................................................54
figureIII.11 : Digramme de rayonnement en Gain en dB (plan  )..................................................................54
figure III.12 : Représentation rectangulaire 3-D du Gain total en dB................................................................54
figureIII.13 : Distribution du courant surfacique sur le patch...........................................................................55
figureIII.14 : Distribution du champ électrique E sur le patch.........................................................................55
figureIII.15 : Distribution du champ électrique E sur le substrat.....................................................................56
figureIII.16 : Distribution du champ électrique E sur le substrat et le patch.....................................................56
figureIII.17 : Variation du coefficient de réflexion en fonction de la fréquence pour différent valeur de
permittivité relative .............................................................................................................................................57
figureIII.18 : Variation du coefficient de réflexion en fonction de la fréquence pour différent valeur de
hauteur h.............................................................................................................................................................57
figureIII.19 : Variation du coefficient de réflexion en fonction de la fréquence pour différent valeur de la
largeur W............................................................................................................................................................58
figureIII.20 : Variation du coefficient de réflexion en fonction de la fréquence pour différent valeur de la
longueur L...........................................................................................................................................................59
figureIII.21 : Variation du coefficient de réflexion en fonction de la fréquence pour différent valeur de
l’épaisseur e........................................................................................................................................................59
figure. III.22 : Variation du coefficient de réflexion en fonction de la fréquence pour patch à conductivité
différente..............................................................................................................................................................60
figureIII.23 : Antenne imprimé alimenté par un transformateur d’impédance selon HFSS.............................62
figure III.24 : Variation de coefficient de réflexion en fonction de la fréquence..............................................62
figure III.25 : Représentation de partie réelle et imaginaire de l’impédance en fonction de la fréquence........62
figure III.26 : Représentation du ROS en dB en fonction de la fréquence.......................................................63
figure III.27 : Représentation 3-D du Gain total en dB......................................................................................63
figure III.28 : Diagramme du Gain en dB (plan  )...........................................................................................63
figure III.29 : Diagramme du Gain en dB (plan  )...........................................................................................64
figureIII.30 : Antenne imprimé alimenté par ligne enterrée selon HFSS.........................................................64
figure III.31 : Schéma de l’antenne à ligne enterrée..........................................................................................64
figure III.32 : Variation de coefficient de réflexion en fonction de la fréquence..............................................65
figure III.33 : Représentation de partie réelle et imaginaire de l’impédance en fonction de la fréquence........65
figure III.34 : Représentation du ROS en dB en fonction de la fréquence.......................................................66
figure III.35 : Représentation 3-D du Gain total en dB......................................................................................66
figureIII.36 : Digramme de rayonnement en Gain en dB (plan  )..................................................................66
figureIII.37 : Digramme de rayonnement en Gain en dB (plan  )..................................................................67
figureIII.38 : Antenne imprimé alimenté par câble coaxiale selon HFSS.........................................................68
figure III.39 : Variation de coefficient de réflexion en fonction de la fréquence..............................................68
figure III.40 : Représentation de partie réelle et imaginaire de l’impédance en fonction de la fréquence........68
figure III.41 : Représentation du ROS en dB en fonction de la fréquence.......................................................69
figure III.42 : Représentation 3-D du Gain total en dB......................................................................................69
figure III.43 : Diagramme du Gain en dB (plan  )...........................................................................................69
figure III.44 : Diagramme du Gain en dB (plan  )...........................................................................................70
Chapitre IV : Modélisation et simulation d’un réseau d’antennes patch en utilisant
HFSS
figure IV. 1 : Couplage entre les antennes patch dans le plan E et H.................................................................74
figure IV.2 : Configuration plan H.....................................................................................................................74
figure IV.3 : Configuration plan E......................................................................................................................75
figure IV.4 : Schémas de deux patchs................................................................................................................76
figure IV.5 : Circuit électrique équivalent de deux patchs.................................................................................76
figure IV. 6 : Exemple de couplage entre deux éléments d’une antenne micro ruban.......................................78
figure IV.7: Réseau d’antenne schématisé avec son réseau de distribution.......................................................80
figure IV.8: Exemples de réseaux d'antennes a 16 éléments..............................................................................81
figure IV.9: Mécanisme des rayonnements parasites.........................................................................................81
figure IV.10: Architecture de base des réseaux d'antennes : a- Réseaux linéaires b- Réseaux plans.............81
figure IV.11: Réseau linéaire d’éléments rectangulaires à alimentation série...................................................83
figure IV.12: Alimentation parallèle de 23éléments rayonnants par 3 étages de diviseurs de puissance...........83
figure IV.13: Bande rayonnante à éléments jointifs à alimentation parallèle....................................................84
figure IV.14: Réseau bidimensionnel formé de réseaux linéaires, a- alimentation série, b- alimentation
parallèle................................................................................................................................................................84
figure IV.15: Réseau bidimensionnel à double alimentation parallèle...............................................................85
figure IV.16 : Réseau linéaires de deux patches alimentés en série par ligne micro ruban (planH)..................86
figure IV.17: Variation de coefficient de réflexion en fonction de la fréquence................................................86
figure IV.18 : Représentation 3-D du Gain total en dB......................................................................................86
figure IV.19 : digramme de rayonnement de Gain total en dB ( plan  )........................................................87
figure IV.20 : digramme de rayonnement en Gain en dB (plan  )..................................................................87
figure IV.21 : Réseau linéaires de deux patches alimentés en série par ligne micro ruban (planE)..................87
figure IV.22: variation de coefficient de réflexion en fonction de la fréquence................................................88
figure IV.23: représentation 3-D du Gain total en dB........................................................................................88
figure IV.24 : digramme de rayonnement de Gain total en dB (plan  ).........................................................88
figure IV.25:digramme de rayonnement en Gain en dB (plan  )....................................................................89
figure IV.26 : Réseau linéaires de deux patches alimentés en parallèle par diviseur de Wilkinson..................90
figure IV.27: variation de coefficient de réflexion en fonction de la fréquence................................................90
figure IV.28: représentation 3-D du Gain total en dB........................................................................................90
figure IV.29 : digramme de rayonnement de Gain total en dB (plan  ).........................................................91
figure IV.30:digramme de rayonnement en Gain en dB (plan  )....................................................................91
figure IV.31 : Réseau linéaires de quatre patches alimentés en parallèle par diviseur de Wilkinson...............91
figure IV.32: variation de coefficient de réflexion en fonction de la fréquence................................................92
figure IV.33: représentation 3-D du Gain total en dB........................................................................................93
figure IV.34 : digramme de rayonnement de Gain total en dB (plan  ).........................................................93
figure IV.35:digramme de rayonnement en Gain en dB (plan  )....................................................................93
figure IV.36 : Réseau bidimensionnel à quatre patchs......................................................................................94
figure IV.37: variation de coefficient de réflexion en fonction de la fréquence................................................94
figure IV.38: représentation 3-D du Gain total en dB.......................................................................................95
figure IV.39 : digramme de rayonnement de Gain total en dB (plan  ).........................................................95
figure IV.40:digramme de rayonnement en Gain en dB (plan  )...................................................................95
figure IV.41 : Réseau bidimensionnel à 8 patchs...............................................................................................96

figure IV.42: variation de coefficient de réflexion en fonction de la fréquence................................................96

figure IV.43: représentation 3-D du Gain total en dB.......................................................................................96
figure IV.44 : digramme de rayonnement de Gain total en dB (plan  ).........................................................97
figure IV.45:digramme de rayonnement en Gain en dB (plan  )....................................................................97
Chapitre V : Réalisations
figure V.1 : patch alimenté par ligne avec transformateur à 8.5GHz................................................................99
figure V.2 : Le modèle de patch alimenté par ligne avec transformateur par le HFSS.....................................99
figure V.3 : Variation du coefficient de réflexion en fonction de la fréquence..............................................143
figure V.4: Digramme de rayonnement de Gain total en dB (plan  )..........................................................101
figure V.5: Digramme de rayonnement mesurée de Gain total (plan  )......................................................101
figure V.6: Digramme de rayonnement en Gain en dB (plan  )...................................................................101
figure V.7: Digramme de rayonnement mesurée de Gain total (plan  )......................................................101
figure V.8 : Représentation 3-D du Gain total en dB.......................................................................................102
figure V.9 : patch alimenté par ligne à 8.5GHz...............................................................................................102
figure V.10 : Le modèle de patch alimenté par ligne a réalisé........................................................................102
figure V.11 : Variation du coefficient de réflexion en fonction de la fréquence............................................103
figure V.12: Digramme de rayonnement de Gain total en dB (plan  ).........................................................103
figure V.13: Digramme de rayonnement mesurée de Gain total (plan  ).....................................................104
figure V.14: Digramme de rayonnement en Gain en dB (plan  )..................................................................104
figure V.15: Digramme de rayonnement mesurée de Gain total (plan  ).....................................................104
figure V.16 : Représentation 3-D du Gain total en dB.....................................................................................105
figure V.17 : Réseau à quatre patchs alimentée en paralléle...........................................................................105
figure V.18 : Le modèle de patch a réalisé.......................................................................................................105
figure V.19 : Variation du coefficient de réflexion en fonction de la fréquence.............................................106
figure V.20: Digramme de rayonnement de Gain total en dB (plan  ).........................................................107
figure V.21: Digramme de rayonnement en Gain en dB (plan  )..................................................................107
figure V.22 : Représentation 3-D du Gain total en dB.....................................................................................107
Tableau I.1: Comparaison entre les différentes techniques d’alimentations……………………………….....8
Tableau II.1: Comparaison des deux modèles analytiques……………………………….............................33
Tableau II.2 : Méthodes numériques pour l’analyse des antennes imprimées……………………………...45
Tableau III.1 : propriétés des matériaux…………………………….............................................................51
Tableau III.2: Paramètres de l’antenne……………………………...............................................................51
Tableau III.3 : Interprétation du résultat de la simulation pour différentes valeurs de  r ………………….57
Tableau III.4 : Interprétation du résultat de la simulation pour différentes valeurs de h…………………...58
Tableau III.5 : Interprétation du résultat de la simulation pour différentes valeurs de W…………………..58
Tableau III.6 : Interprétation du résultat de la simulation pour différentes valeurs de L…………………...59
Tableau III.7 : Interprétation du résultat de la simulation pour différentes valeurs de e……………………60
Tableau III.8 : Interprétation du résultat de la simulation pour différentes valeurs de conductivité………..61
Tableau III.9 : Propriétés des matériaux………………………………………………………………….....61
Tableau III.10: Paramètres de l’antenne………………………………………………………………….....61
Tableau III.11: Propriétés des matériaux……………………………………………………………………67
Tableau III.12: Paramètres de l’antenne…………………………………………………………………….67
Tableau III.13: Exemple d’application de l’antenne imprimée……………………………………………..71
Tableau IV.1 : propriétés des matériaux………………………………………………………………….....85
Tableau IV.2: Paramètres de l’antenne……………………………………………………………………...85
Tableau IV.3 : propriétés des matériaux………………………………………………………………….....89
Tableau IV.4: Paramètres de l’antenne………………………………………………………………….......89
Tableau IV.5 : propriétés des matériaux………………………………………………………………….....92
Tableau IV.6: Paramètres de l’antenne……………………………………………………………………...92
Tableau V.1 : Propriétés des matériaux………………………………………………………………….....100
Tableau V.2: Paramètres de l’antenne………………………………………………………………….......100
Tableau V.3 : propriétés des matériaux…………………………………………………………………….102
Tableau V. 4 : Paramètres de l’antenne…………………………………………………………………….103
Tableaux V.5 : propriétés des matériaux…………………………………………………………………...152
Tableau V.6: Paramètres de l’antenne………………………………………………………………….......152