République Tunisienne

Ministère de la défense nationale

Académie militaire Fondouk Jedid

Rapport du projet de fin d’année

Département de l’électromécanique

Titre du projet:
Comportement des
poutres en matériaux
FGM flexoélectriques.

Élaboré par:

Gandouzi Wassel

Encadré par:

Mr. Sami Abdelkahlek

Academic year: 2022/2023

Promotion:
Mohamed Moncef Bey (2019/2024)
Dédicace
“ Ce projet de fin d’année est dédié à mes chers parents qui
m’ont toujours poussé et motivé dans mes études, Sans
eux, nous n’aurons certainement pas faire d’études
longues,
Ce projet de fin d’année représente donc l’aboutissement
du soutien et des encouragements qu’ils m’ont prodigués
tout au long de mes études. Qu’ils en soient remerciés par
cette trop modeste dédicace,
C’est un moment de plaisir de dédier cet œuvre, à nos
frères, en signe d’amour, de reconnaissance et de gratitude
pour le dévouement et les sacrifices dont vous avez fait
toujours preuve à notre égard. Merci. ”
- Wassel

I
Remerciements
Je tiens particulièrement à remercier Allah le
tout puissant, ce mémoire n’aurait jamais été
réalisé sans sa bénédiction.
Je tiens à remercier ma famille de m’avoir
soutenu, je ne serai jamais assez reconnaissant
envers mes parents qui ont toujours tout mis en
œuvre pour qu’on s’épanouisse dans tout ce
qu’on entreprend.
J’adresse mes remerciements à mon encadrant
Monsieur Sami Abdelkhalek, de m’avoir
accordé toute sa confiance et de m’avoir guidée
tout au long de cette PFA, pour son aide
consistante, ses conseils judicieux, et pour ses
remarques précises.
Enfin, je tiens à remercier toute personne m’a
aidé de près ou de loin durant mon travail.

II
Table des matières

1 Introduction générale .............................................................................. 1

2 Etat de l’art et histoire ............................................................................. 3

2.1 Histoire et évolutions ............................................................................. 4

2.2 Les matériaux FGM ................................................................................ 5
2.3 Matériaux FGM-Flexo............................................................................. 6
2.4 Types des matériaux FGM-Flexoélectriques ........................................... 7
2.4.1 Matériaux solides............................................................................. 7
2.4.2 Les cristaux liquides ......................................................................... 8
2.4.3 Les polymères .................................................................................. 8
2.5 MAN....................................................................................................... 9
2.5.1 MAN (définition) .............................................................................. 9
2.5.2 Le MAN pour le calcul des comportements des matériaux FGM-
flexoeléctrique........................................................................................... 10
3 Formulation mathématique des poutres flexoélectriques .......................13

1.1 Choix de la forme de la structure ....................................................... 14

1.2 Mise en équation du système ............................................................. 14
1.3 La procédure automatique de calcul des termes MAN ..................... 17
1.4 Validation de la MAN pour les poutres élastiques avec non-linéarité
géométrique ................................................................................................ 21
4 Test et validation de l'analyse des comportements des matériaux FGM
(MAN automatique) ......................................................................................23

1.1 Validation de FGM par rapport les matériaux homogènes (Force

uniquement) ................................................................................................ 24

III
 1.2 Étude de cas d'un matériau FGM composé de deux matériaux
différents (force uniquement) .................................................................... 26
1.3 Validation de FGM par rapport les matériaux homogènes (voltage
uniquement) ................................................................................................ 28
1.4 Étude de cas d'un matériau FGM composé de deux matériaux
différents (voltage uniquement) ................................................................ 29
4.1 Effet flexoélectrique dans les nano-poutres FGM ................................. 32
5 Conclusion et perspectives......................................................................33

6 Bibliographie ..........................................................................................36

IV
Table des figures

Figure 1 Poutre à extrémités fixes soumise à une force verticale. ................... 14

Figure 2 Déflexion de la micro-poutre à extrémité fixe : (a) Forme de la
déflexion, (b) Déflexion maximale en fonction de la charge. ............................ 21
Figure 3 Flambage de la poutre élastique du Tableau 2. .................................. 22
Figure 4 Déflexion maximale en fonction de la charge. ................................... 25
Figure 5 Déflexion maximale en fonction de la charge(cas flexo/cas FGM). .... 25
Figure 6 Déflexion maximale en fonction de la charge (cas matériau homogène
/cas FGM). ........................................................................................................ 27
Figure 7 Le mode de déflexion sous tension de la nano—poutre ( cas matériau
homogène /cas FGM). ...................................................................................... 27
Figure 8 Déflexion maximale en fonction de la tension éléctrique. .................. 28
Figure 9 Déflexion maximale en fonction de la tension éléctrique (cas de
flexo/cas FGM) ................................................................................................. 29
Figure 10 Déflexion maximale en fonction de la charge (cas1er matériau
homogène /cas FGM/cas 2eme matériau homogène ). .................................... 31
Figure 11 Le mode de déflexion sous tension de la nano—poutre ( cas 2em
matériau homogène /cas FGM). ....................................................................... 31

V
Liste des tables

Tableau 1 Les propriétés de la micro-poutre élastique utilisées pour la

comparaison avec le modèle basé sur la méthode DQM de Baroudi et al. […]
sont les suivantes : ........................................................................................... 21
Tableau 2 Propriétés de la poutre élastique utilisées pour la validation des
tests de flambage : ........................................................................................... 22
Tableau 3 Les caractéristiques de la poutre nanométrique flexoélectrique
étudiée par Yan et al. [...] ................................................................................. 24
Tableau 4 Les caractéristiques de la poutre nanométrique flexoélectrique
(deuxième matériau) ........................................................................................ 26
Tableau 5 Les caractéristiques de la micro-poutre flexoélectrique (1er matériau)
......................................................................................................................... 30
Tableau 6 Les caractéristiques de la micro-poutre flexoélectrique (2eme
matériau) ......................................................................................................... 30

VI
1 Introduction générale

Dans ce projet, nous présentons une technique novatrice qui facilite le

développement d'un modèle numérique pour les poutres FGM-
flexoélectriques, prenant en compte la non-linéarité géométrique grâce à la
Méthode Asymptotique Numérique (MAN). Les poutres FGM-flexoélectriques
sont des structures composites constituées de matériaux à gradient de fonction
(FGM) qui présentent des propriétés flexoélectriques.

La non-linéarité géométrique est un aspect essentiel à considérer dans l'analyse

des structures flexibles, car elle se manifeste lorsque les déformations sont
suffisamment importantes pour ne pas être négligées. Dans le cas des poutres
FGM-flexoélectriques, les effets de la non-linéarité géométrique peuvent avoir
une influence significative sur leur comportement mécanique et électrique.

La Méthode Asymptotique Numérique est une approche mathématique qui

permet de résoudre des problèmes complexes en décomposant le système en
plusieurs sous-systèmes plus simples, en utilisant des approximations
asymptotiques. Cette méthode offre une alternative efficace à la méthode des
éléments finis traditionnelle, en réduisant la complexité du modèle et en
améliorant la précision des résultats.

Dans ce rapport, nous nous appuyons sur la théorie d'Euler-Bernoulli pour

modéliser le comportement des poutres FGM-flexoélectriques. La théorie
d'Euler-Bernoulli est largement utilisée pour l'analyse des poutres, en
considérant des hypothèses simplificatrices telles que la petite déformation et
la linéarité des matériaux. Cependant, ces hypothèses peuvent être limitantes

1
dans le cas des poutres FGM-flexoélectriques en raison de leur comportement
non-linéaire.

Afin de valider notre modèle, une comparaison avec des travaux antérieurs est
réalisée. Cette validation permet de démontrer la fiabilité de notre approche et
l'adéquation de notre modèle aux résultats expérimentaux et aux modèles
analytiques existants.

Le reste du rapport est organisé comme suit : dans la section suivante, nous
détaillons la méthodologie utilisée pour développer le modèle numérique en
utilisant la Méthode Asymptotique Numérique et la théorie d'Euler-Bernoulli.
Ensuite, nous présentons les résultats de la validation et discutons de leur
pertinence.

2
Chapitre 1

2 Etat de l’art et histoire

3
2.1 Histoire et évolutions

Les matériaux FGM, ou matériaux à gradient de fonction, sont des matériaux

composites conçus avec des variations progressives de composition, de
structure et de propriétés à travers leur épaisseur. Ces matériaux sont conçus
pour combiner les avantages de différents matériaux et offrir des performances
améliorées dans diverses applications. L'histoire et l'évolution des matériaux
FGM remontent aux années 1980, et depuis lors, ils ont suscité un intérêt
croissant dans le domaine de la science des matériaux.

L'idée fondamentale des matériaux FGM est de créer une transition graduelle
entre deux matériaux différents plutôt que de les superposer brusquement.
Cela permet de minimiser les contraintes et les tensions causées par les
différences de dilatation thermique et de propriétés mécaniques entre les
matériaux. En utilisant cette approche, les matériaux FGM peuvent présenter
des combinaisons de propriétés uniques qui ne sont pas atteignables avec les
matériaux traditionnels.

L'une des premières applications des matériaux FGM a été dans le domaine des
revêtements thermiques pour les turbines à gaz. Les revêtements FGM
permettent une transition graduelle entre le matériau de base et le
revêtement, améliorant ainsi l'adhérence et la résistance aux contraintes
thermiques. Cela a conduit à une meilleure performance et une durée de vie
accrue des composants de turbine.

Au fil du temps, les matériaux FGM ont été explorés dans de nombreux autres
domaines, tels que les revêtements protecteurs, les matériaux structurels, les
matériaux de friction, les matériaux biomédicaux, etc. Les chercheurs ont
étudié différentes combinaisons de matériaux, y compris les céramiques, les
métaux et les polymères, pour créer des matériaux FGM adaptés à chaque
application spécifique.

Les méthodes de fabrication des matériaux FGM ont également évolué au fil du
temps. Les techniques couramment utilisées comprennent la métallurgie des
poudres, la pulvérisation thermique, la méthode sol-gel, la fabrication additive
et d'autres processus de fabrication avancés. Ces techniques permettent de

4
contrôler précisément les variations de composition et de structure à travers
l'épaisseur du matériau.

L'évolution récente des matériaux FGM comprend des avancées telles que
l'utilisation de nanomatériaux pour créer des gradients de taille de grain, la
conception assistée par ordinateur pour optimiser les profils de gradients et
l'exploration de nouvelles combinaisons de matériaux pour obtenir des
propriétés encore meilleures.

En résumé, les matériaux FGM sont des matériaux composites conçus avec des
variations progressives de composition et de propriétés à travers leur épaisseur.
Leur histoire remonte aux années 1980, et depuis lors, ils ont été largement
explorés et utilisés dans diverses applications. Les avancées dans la
compréhension des matériaux et les techniques de fabrication ont permis de
développer des matériaux FGM aux propriétés uniques, offrant ainsi de
nouvelles possibilités dans le domaine de la science des matériaux.

2.2 Les matériaux FGM

Les matériaux FGM (Fonctionally Graded Materials) sont des matériaux qui
présentent une transition graduelle des propriétés physiques et chimiques à
travers leur épaisseur. Contrairement aux matériaux homogènes, où les
propriétés sont uniformes dans tout le matériau, les matériaux FGM
permettent de créer des gradients de propriétés sur une échelle
macroscopique.

La transition graduelle des propriétés dans les matériaux FGM peut être
obtenue de différentes manières. Une méthode couramment utilisée consiste à
mélanger deux matériaux de base avec des propriétés différentes, tels que la
céramique et le métal, en variant progressivement les proportions des
matériaux le long de l'épaisseur du matériau.

Les matériaux FGM présentent plusieurs avantages. Ils peuvent être conçus
pour répondre à des exigences spécifiques en matière de propriétés, telles que
la résistance, la conductivité thermique, la conductivité électrique, etc., dans
différentes parties d'un matériau. Par exemple, un matériau FGM peut être
5
conçu pour avoir une résistance élevée d'un côté et une ductilité accrue de
l'autre côté, offrant ainsi une combinaison de propriétés mécaniques
souhaitables.

Les applications des matériaux FGM sont nombreuses. Ils sont utilisés dans des
domaines tels que l'aérospatiale, l'ingénierie des matériaux, la fabrication
d'outils et d'objets de haute technologie. Par exemple, les aubes de turbine
dans les moteurs d'avion peuvent être fabriquées en utilisant des matériaux
FGM pour améliorer la résistance et la durabilité.

En résumé, les matériaux FGM sont des matériaux qui présentent une
transition graduelle des propriétés physiques et chimiques à travers leur
épaisseur. Ils offrent des avantages en termes de conception et de combinaison
de propriétés, ce qui les rend utiles dans de nombreux domaines d'application.

2.3 Matériaux FGM-Flexo

Les matériaux FGM flexoélectriques sont une classe spécifique de matériaux

fonctionnellement gradués qui combinent les propriétés des matériaux FGM
avec celles des matériaux flexoélectriques. Le flexoélectricité est un
phénomène qui se produit dans certains matériaux non centrosymétriques, où
une polarisation électrique est générée en réponse à un gradient de contrainte
mécanique.

Dans les matériaux FGM flexoélectriques, on cherche à créer une transition

graduelle des propriétés flexoélectriques à travers l'épaisseur du matériau. Cela
signifie que la quantité de flexoélectricité générée varie progressivement dans
le matériau. Par exemple, on peut concevoir un matériau FGM flexoélectrique
dans lequel la flexoélectricité est plus élevée d'un côté et diminue
progressivement de l'autre côté.

Les propriétés flexoélectriques des matériaux FGM flexoélectriques peuvent

être utilisées dans diverses applications. Par exemple, ils peuvent être utilisés
pour concevoir des capteurs de pression flexoélectriques qui convertissent les

6
gradients de contrainte en signaux électriques. Ces capteurs peuvent être
utilisés dans les systèmes de surveillance et de contrôle de la pression, tels que
les dispositifs médicaux, les robots souples, etc.

Les matériaux FGM flexoélectriques peuvent également être utilisés dans les
actionneurs, qui sont des dispositifs capables de convertir l'énergie électrique
en énergie mécanique. Lorsqu'un champ électrique est appliqué à un matériau
FGM flexoélectrique, il génère un gradient de contrainte mécanique qui peut
être utilisé pour produire des mouvements et des déformations contrôlés.

En résumé, les matériaux FGM flexoélectriques sont des matériaux

fonctionnellement gradués qui combinent les propriétés des matériaux FGM
avec celles des matériaux flexoélectriques. Ils présentent une transition
graduelle des propriétés flexoélectriques à travers leur épaisseur et peuvent
être utilisés dans des applications telles que les capteurs de pression et les
actionneurs.

2.4 Types des matériaux FGM-Flexoélectriques

Par rapport à la piézoélectricité, qui n'existe que dans les 20 groupes ponctuels
des systèmes qui sont non centrosymétriques, la flexoélectricité a l'avantage
d'exister dans les 32 groupes ponctuels cristallins.

2.4.1 Matériaux solides

En raison de leurs applications potentielles, la classe des céramiques
piézoélectriques a suscité un grand intérêt car elles peuvent être fabriquées
pour des applications spécifiques en fonction de leurs propriétés particulières.
Le titanate de baryum BaTiO3 présente un énorme coefficient flexoélectrique
ainsi qu'un fort champ diélectrique, ce qui en fait un bon exemple pour les
applications ultérieures de la flexoélectricité. C'est pourquoi de nombreux
chercheurs le choisissent comme matériau d'échantillon pour leurs études .
Maranganti et Sharma ont rapporté que les trois composantes indépendantes
du tenseur flexoélectrique BaTiO3 : µ11, µ12 et µ14 sont respectivement de
l'ordre de 0,15 nC/m, -5,5 nC/m et -1,9 nC/m. Dans d'autres études , les auteurs
ont utilisé la méthode de nano-indentation pour estimer le coefficient
flexoélectrique transversal et ont constaté qu'il était d'environ 4 µC/m. Fu et al.
sont les premiers à déterminer expérimentalement le coefficient
7
flexoélectrique longitudinal µ11 à la fois pour le titanate de baryum BaTiO3 et
le titanate de strontium SrTiO3.

D'autres piézocéramiques ont fait l'objet de mesures expérimentales , telles

que le niobate de magnésium plomb, le titanate de baryum-strontium et le
titanate de plomb-zirconate, évaluant leurs coefficients flexoélectriques de
l'ordre de 10^-6 C/m. Le groupe de Sharma a également prédit les coefficients
flexoélectriques pour les cristaux non ferroélectriques (GaAs, GaP, ZnS, NaCl et
KCl) et a trouvé des coefficients flexoélectriques dans la plage de 10^-12 à 10^-
11 C/m, tandis que la plage pour les ferroélectriques est de 10^-9 à 10^-8 C/m.

2.4.2 Les cristaux liquides

En réalité, la flexoélectricité a été découverte pour la première fois dans les
cristaux liquides par Meyer, ce qui a été interprété par la relation entre les
distorsions de courbure et la polarisation électrique. L'étude théorique et
expérimentale de la flexoélectricité dans les cristaux liquides (CL) est un défi .
En particulier, l'intérêt a été porté sur les cristaux liquides nématiques, par
exemple, Takahashi et al. ont proposé une nouvelle méthode pour mesurer les
coefficients flexoélectriques des nématiques (CL). Nous pouvons également
citer le travail de Kotov et al. , qui a démontré que pour une orientation plane,
l'effet flexoélectrique entraîne un nouveau type de transition d'orientation.
Harden et al. ont utilisé une méthode directe via le courant électrique produit
par la déformation périodique de flexion de la surface de confinement des CL et
ont trouvé un coefficient flexoélectrique géant (35 × 10^-9 C/m) dans un CL
nématique à noyau courbé.

2.4.3 Les polymères

La preuve de la flexoélectricité dans les polymères est fournie par les travaux de
Marvan et Havránek , où un modèle simple de mécanisme d'orientation des
dipôles contrôlé par un volume libre anisotrope a été utilisé et le gradient de
déformation a été expérimentalement réalisé dans une pyramide tronquée
découpée à partir d'élastomères soumis à une pression axiale. L'ordre de
grandeur du coefficient flexoélectrique pour les élastomères étudiés était
évalué à 10^-11 - 10^-10 C/m. Plus tard, dans ses travaux, Baskaran et al. ont
montré la flexoélectricité géante dans les films de polymère en polyfluorure de
vinylidène (PVDF) en phase α. Dans un autre travail, Baskaran et al. ont tenté
de calculer le coefficient flexoélectrique transversal effectif des films de PVDF
8
sur la base d'un essai de flexion. La méthode consiste à coller le film de
polymère avec une électrode de mise à la terre et une électrode de pointe sur
une poutre en porte-à-faux pour former une structure unimorphe. Chu et
Salem ont récemment étudié la flexoélectricité dans plusieurs polymères
thermoplastiques et thermodurcissables (dont le PVDF, le PET orienté, le
polyéthylène et l'époxy) en adoptant l'essai de flexion en porte-à-faux avec une
électrode supérieure entièrement recouverte. Les coefficients flexoélectriques
de ces polymères sont de l'ordre de 10^-9 à 10^-8 C/m.

2.5 MAN
2.5.1 MAN (définition)
La méthode asymptotique numérique est une approche mathématique qui vise
à estimer de manière précise et efficace les solutions de problèmes complexes
en se basant sur les comportements asymptotiques des variables ou des
paramètres du système. Cette méthode combine à la fois des techniques
numériques et des analyses asymptotiques afin d'obtenir des estimations
fiables des solutions.

La méthode asymptotique numérique s'appuie sur les principes suivants :

Identification des régimes asymptotiques : Cette étape consiste à déterminer

les régimes limites du système, où les variables ou les paramètres du problème
deviennent très grands ou très petits. Ces régimes limites sont souvent
caractérisés par des propriétés spécifiques du système, ce qui permet de
simplifier les équations ou d'ignorer certains termes.

Décomposition en termes asymptotiques : Les équations du système sont

décomposées en une série de termes asymptotiques. Chaque terme est associé
à une puissance appropriée des variables ou des paramètres asymptotiques.
Cette décomposition permet de séparer les contributions dominantes des
contributions négligeables dans les équations.

Combinaison des termes asymptotiques : Les termes asymptotiques obtenus

dans l'étape précédente sont combinés de manière à obtenir une
approximation des solutions du problème. Cette combinaison peut impliquer
des techniques numériques pour résoudre les équations, ainsi que des analyses
asymptotiques pour évaluer les contributions des différents termes.

9
En utilisant ces principes, la méthode asymptotique numérique permet
d'obtenir des estimations précises et efficaces des solutions de problèmes
complexes. Elle exploite les comportements asymptotiques du système pour
simplifier les équations et identifier les contributions les plus importantes, ce
qui facilite la résolution numérique et améliore la compréhension du problème.

Résolution des termes asymptotiques : Chaque terme de la série asymptotique

est résolu individuellement en utilisant des méthodes numériques adaptées à
sa nature. Cela peut comprendre des techniques de résolution analytique, des
approximations polynomiales, des méthodes de perturbation ou d'autres
méthodes spécifiques à chaque terme.

Combinaison des solutions : Les solutions des termes asymptotiques sont

combinées pour obtenir une approximation globale de la solution du problème.
Cette combinaison peut être réalisée à l'aide de techniques de sommation ou
d'assemblage des termes asymptotiques.

2.5.2 Le MAN pour le calcul des comportements des matériaux FGM-

flexoeléctrique

[Link] Avantage de la MAN

Calcul automatique du pas

À chaque étape, la méthode asymptotique numérique (MAN) cherche à

résoudre le problème en construisant soigneusement des séries de Taylor pour
calculer les inconnues. Nous abordons un problème qui est présent depuis le
début et nous essayons de l'aborder avec des critères précis. Nous avons
l'impression de nous rapprocher d'une solution, même si celle-ci reste encore
un peu floue. Trouver une solution non linéaire n'est pas une tâche facile. Cela
nécessite de reconstruire une série de courbes polynomiales. Cependant,
même si nous pouvons en faire une approximation, cela ne signifie pas que
nous avons décrypté le code. Travailler avec des fractions rationnelles n'est pas
seulement une théorie, c'est une réalité grâce aux approximations de Pade
(PADE). Ces approximations des équations aux dérivées partielles (EDP) sont
organisées en une série d'ordres qui peuvent être utilisés dans de nombreux

10
calculs complexes. L'utilisateur a la responsabilité de décider à l'avance de
l'ordre à utiliser, et plusieurs auteurs ont souligné qu'une sélection appropriée
de ce type de séquence peut être très bénéfique. Selon Boutyour et al., le
nombre optimal d'activités semble se situer entre quinze et vingt. Cette plage
est assez pratique, si vous me le demandez. L'amplitude du pas correspond
précisément au rayon de convergence de cette série spécifique. Il est bien
connu que cette distance varie en fonction des différentes séquences
numériques.

Calcul direct de la solution à chaque pas

À chaque itération, les méthodes de prédiction-correction (PC) consistent à

effectuer une série de prédictions linéaires du problème en estimant le champ
de solution, puis à calculer l'erreur correspondante (résidu) et à corriger par la
suite. Une itération de prédiction-correction comprend donc une prédiction
suivie d'une correction. Dans un code basé sur la méthode des éléments finis
(EF), une correction implique de multiplier le résidu par l'inverse de la matrice
tangente, également appelée matrice jacobienne. Cette matrice tangente
représente la pente du problème non linéaire pour chaque valeur des
inconnues considérée. Les inversions de matrices tangentes rendent les calculs
coûteux en termes de temps CPU lorsqu'il s'agit de méthodes itératives. En
revanche, avec la méthode asymptotique numérique (MAN), un calcul direct de
la solution sur un pas est avantageux, évitant ainsi les inversions de matrices et
les coûts associés.

[Link] Inconvénient de la MAN

Le problème tangent utilisé dans les algorithmes itératifs pour la prédiction
n'est autre que le premier terme de la série de Taylor utilisée dans la méthode
asymptotique numérique (MAN). C'est pourquoi on considère que la MAN est
une généralisation des méthodes de prédiction-correction. Les problèmes
d'ordres supérieurs qui apparaissent dans le développement en série sont dus à
la non-linéarité du problème. Ils doivent généralement être déterminés
analytiquement et implémentés dans les codes de calcul. Bien que la MAN offre
de nombreux avantages pour la modélisation des problèmes non linéaires, son
inconvénient principal réside dans la complexité du développement de ces
termes supplémentaires. Par exemple, dans le cas du calcul de poutres
élastiques soumises à de grands déplacements (non-linéarité géométrique), le
11
développement analytique des termes d'ordres supérieurs peut être fastidieux.
Cet inconvénient est encore plus prononcé lors de la résolution de problèmes
impliquant des couplages multiphysiques, tels que les structures
flexoélectriques. C'est pourquoi l'objectif de ce projet est l'automatisation de la
méthode asymptotique numérique (MAN).

Dans ce rapport, nous exposons une technique qui facilite le développement

d'un modèle numérique pour les poutres FGM-flexoélectriques prenant en
compte la non-linéarité géométrique grâce à la MAN. Le travail est basé sur la
méthode des éléments finis pour la théorie d'Euler-Bernoulli. Une validation
avec des travaux antérieurs est effectuée pour démontrer la fiabilité du modèle.
De plus, nous étudions le comportement des matériaux FGM.

12
Chapitre 2

3 Formulation mathématique des poutres

flexoélectriques

13
1.1 Choix de la forme de la structure

Dans ce rapport, nous nous concentrons particulièrement sur les poutres Euler-
Bernoulli à extrémités fixes, initialement droites, soumises à une force verticale
à mi-longueur (voir Figure 1) et avec des sections transversales rectangulaires.
Elles représentent le problème que nous essayons de présenter dans ce travail,
et une extension à d'autres types de charges et de conditions aux limites est
possible.

Figure 1 Poutre à extrémités fixes soumise à une force verticale.

1.2 Mise en équation du système

L'énergie potentielle de la structure nano flexoélectrique est donnée par

l'équation suivante [Baroudi]:

1
H =  WL −   +   P d 
 2
1 (1)
=  W L − ,2i + ,i Pi d  i = x, y , z
 2

14
Où  est la permittivité,  est le potentiel électrique et P est la densité de
champ de polarisation. W L est la densité d'énergie interne qui est :

1 1 1 1
WL = P  a  P +  : c :  +  : d  P + P  f :  +   : g :  
2 2 2 2
1 1 1 (2)
= aij Pi Pj + cijkl  ij  kl + dijk  ij Pk + f ijkl  ij , k Pl + g ijklmn  ij , k  lm , n
2 2 2

Où  est le tenseur de contrainte, a est le tenseur de susceptibilité

diélectrique réciproque d'ordre deux, c est le tenseur de constantes élastiques
d'ordre quatre, d est le tenseur de constantes piézoélectriques d'ordre trois, f
est le tenseur de propriétés flexoélectriques d'ordre quatre et g représente
l'effet élastique non local du gradient de déformation. L'hypothèse de la poutre
conduit à une expression réduite de la variation de l'enthalpie électrique
[Baroudi]:

 H =  T11 xx + T 111 xx , x + T 113 xx , x d  +



 ( P −  ) 
 x ,x ,x d +
(3)
 ( P −  ) 
 z ,z ,z d

Où Tij est le tenseur de contrainte et T ijk est le tenseur de contrainte du

gradient de déformation, définis respectivement par :

15
 W L
Tij =
 ij
W L
T ijk = (4)
 ij , k

On prend c1111 = c11 , d111 = d11 , d113 = d13 , f1111 = f11 , f1133 = f14 et g111111 = g113113 = g11 , les

relations (2) et (4) donne:

T11 = c11 xx + d11 Px + d13 Pz −  H +  Wnc = 0

T 111 = f11 Px + g11 xx , x
(5)
T 113 = f14 Pz + g11 xx , x

Grâce à l'hypothèse de la poutre, seules les composantes Px et Pz du vecteur de

polarisation P sont non nulles, et elles sont exprimées comme suit [Sourour]:

( , x + d11 xx + f11 xx , x )
1
Px = −
a11

( , z + d13  xx + f14  xx , z )
1
Pz = − (6)
a33

16
L'équilibre électromécanique de la structure est régi par l'équation suivante :

Où Wnc est le travail des charges externes. Soit :

w 2
 xx = u  − z  +
2
 = w
  1 
 ( x, z ) =  ( x ) cos 
(7)
z  +  + z hv
h  2 

L'assemblage et le reste des calculs se font à l'aide de machines.

1.3 La procédure automatique de calcul des termes

MAN

MAN nécessite le développement de U , N et  en fonction de la longueur de

pas a comme suit :

n
U = U0 +  a pU p
p =1
n
N = N0 +  a p N p
p =1 (9)
n
 = 0 +  a p  p
p =1

17
où n est l'ordre de la série. Ces expressions asymptotiques sont injectées dans
l'équation (11) et nous obtenons des systèmes pour chaque ordre p comme
suit :
Order 1:

 EAu1 u dx + EAw0 w1 u dx = 0

 L L

 L EIw1 wdx + L ( N 0 + EAw0 ) w1 wdx + L EAw0 u1 wdx = 1 F  w ( xF )
2 (10)

Order p  1:

 1 p −1
 L ( EAu 
p + EAw 0 p )
 w   u dx = − L 2 
EA wiw p −i  u dx
 i =1

 L EIwp w dx + L ( N 0 + EAw0 ) wp  w dx + L EAw0u p  w dx =
2

 (11)
p −1
  1 
  p F  w ( xF ) −    N i + EA w0 wi  wp −i  w  dx
 i =1  2 
L

L'interpolation d'Hermite (cubique) est utilisée pour les déplacements u et w .

Ainsi, la discrétisation des éléments finis de u et w est écrite comme la relation
suivante :

u = H T uˆ e
(12)
w = H T wˆ e
18
où uˆ e and wˆ e sont les vecteurs des inconnues nodales de chaque élément de
poutre, où uˆ e contient les déplacements nodaux û et les extensions û  , tandis
que wˆ e contient les flèches nodales ŵ et les rotations ŵ . H T est la transposée
du vecteur de fonction de forme. Ainsi, la forme discrète des termes des
systèmes (12) et (13) pour chaque élément e est écrite comme suit :

Order 1:

 EAu1 u dx + EAw0 w1 u dx =  uˆ eT K uue uˆ1e +  uˆ eT K uw

  wˆ 1e
e

 le le

 l EIw1 wdx + l ( N 0 + EAw0 ) w1 wdx + l EAw0u1 wdx =
2

 e e e (13)
  wˆ K Eww wˆ 1 +  wˆ K Gww
eT e e eT e
wˆ 1e +  wˆ eT K wu
e
uˆ1e


Order p  1:

 EAu p  u dx + EAw0 wp  u dx =  uˆ eT K uue uˆ ep +  uˆ eT K uw

 L L
e
wˆ ep
 1 p −1
− L EA  wiwp −i  u dx =  uˆ Fnle
eT up

 2 i =1

 L EIwp w dx + L EA ( N 0 + DAw0 ) wp  wdx + (14)
2


  EAw0 u p  wdx =  wˆ eT K Eww wˆ ep +  wˆ eT K Gww wˆ ep +  wˆ eT K wu uˆ ep
e e e

 L
p −1
  1 
− L   N i + DA 2 w0 wi  wp −i  w dx =  wˆ Fnle
eT wp

 i =1  

Où :

19
 2
K uue = EA H H T dr
le le
2
e
K uw = EA w0 H H T dr
le le
2
e
K wu = EA w0 H H T dr
le le
8
= EI  H H T dr
e
K Eww
le 2
le (15)

e
K Gww = EA
le
2
le
( N 0 + EAw02 ) H H T dr

e
K ww = K Eww
e
+ K Gww
e

La variable r est la coordonnée locale, le est la longueur de l'élément de poutre

et K est la rigidité géométrique élémentaire. L'assemblage (représenté par
e
Gww

l'opérateur U ( • ) ) fournit les systèmes globaux comme suit :

e

Ordre 1 :

 K uu uˆ1 + K uw wˆ 1 = 0
 (16)
 K ww wˆ 1 + K wu uˆ1 = 1 F
w

Order p  1:

 K uu uˆ p + K uw wˆ p = 0

(17)

 K ww wˆ p + K wu uˆ p =  p F + Fnl
w wp


Où :

K uu = U K uu
e
, K ww = U K ww
e
, K uw = K wu
T
= U K uw
e
, Fnlwp = U Fnle
wp
(18)
e e e e

Le vecteur Fw contient les charges nodales selon l'axe y .

20
1.4 Validation de la MAN pour les poutres élastiques
avec non-linéarité géométrique

Le modèle avec MAN automatique dans le contexte de la non-linéarité

géométrique est testé pour la poutre qui présente les caractéristiques
détaillées dans le Tableau 1. Les résultats sont proches de ceux fournis par une
méthode de quadrature différentielle (DQM) basée sur un modèle développé
par Baroudi et al. (voir Figure 2). Nous avons également constaté un accord
avec les résultats du développement de l'ANM classique.
Tableau 1 Les propriétés de la micro-poutre élastique utilisées pour la comparaison avec le modèle basé sur la méthode
DQM de Baroudi et al. sont les suivantes :

Longueur L (m) 400 × 10−6

Largeur b (m) 20 × 10−6
Épaisseur h (m) 20 × 10−6
Module de Young E (GPa) 131
Nombre d'éléments 40

Figure 2 Déflexion de la micro-poutre à extrémité fixe : (a) Forme de la déflexion, (b) Déflexion maximale en fonction de la
charge.

21
Le modèle est également testé pour l'instabilité de flambage des poutres
élastiques. Nous simulons le cas d'une force de compression axiale pour
laquelle on peut calculer analytiquement la force critique de flambage Fc . La
courbe de post-flambage pour une poutre à extrémité fixe avec les propriétés
géométriques et matérielles décrites dans le Tableau 2 est présentée dans la
Figure 3. Elle montre un bon accord entre nos résultats numériques et la valeur
analytique de la force critique Fc .
Tableau 2 Propriétés de la poutre élastique utilisées pour la validation des tests de flambage :

Longueur L (m) 500 × 10−3

Largeur b (m) 20 × 10−3

Épaisseur h (m) 1 × 10−3

Module de Young E (GPa) 200

Nombre d'éléments 40

Avec Fc analytique = 52.637 N .

Figure 3 Flambage de la poutre élastique du Tableau 2.

22
Chapitre 3

4 Test et validation de l'analyse des

comportements des matériaux FGM (MAN
automatique)

23
1.1 Validation de FGM par rapport les matériaux

homogènes (Force uniquement)

Nous considérons la poutre nanométrique à extrémité fixe avec les

caractéristiques présentées dans le Tableau 3, soumise à une force verticale
F = 300nN . La Figure 4 présente les réponses pour un seul matériau que
nous avons validé dans le Chapitre 2.

Tableau 3 Les caractéristiques de la poutre nanométrique flexoélectrique étudiée par Yan et al. [...]

Longueur L (m) 400 × 10−9

Largeur b (m) 20 × 10−9
Épaisseur h (m) 20 × 10−9
Module de Young E = c11 (GPa) 131
Constantes diélectriques a11 (Vm/C) 0.8971 × 108
a33 (Vm/C) 0.788 × 108
Constantes piézoélectriques d11 (V/m) 0
d31 (V/m) 1.87 × 108
Constantes flexoélectriques f11 (V) 0
f14 (V) 5
Échelle de longueur élastique l0 (m) 0
non-locale
Perméabilité  (C2 V-1 m-1) 1.115 × 10−8
Nombre d'éléments 40

24
 Figure 4 Déflexion maximale en fonction de la charge.

Maintenant, nous utilisons le calcul pour les matériaux FGM composés de deux
matériaux supposés différents. Cependant, dans ce cas, nous utilisons deux
matériaux ayant des caractéristiques similaires, et nous devons constater que
les deux résultats sont très proches. Cela est illustré dans la Figure 5.

Figure 5 Déflexion maximale en fonction de la charge(cas flexo/cas FGM).

25
1.2 Étude de cas d'un matériau FGM composé de deux

matériaux différents (force uniquement)

Dans cette étude, nous nous intéressons à un matériau FGM (Matériau

Gradient Fonctionnel) constitué de deux matériaux distincts. L'objectif est
d'analyser le comportement mécanique de ce matériau composite et d'évaluer
ses performances par rapport à des applications spécifiques.

Les caractéristiques du premier matériau sont répertoriées dans le Tableau 3.

Les caractéristiques du deuxième matériau sont répertoriées dans le Tableau 4

Tableau 4 Les caractéristiques de la poutre nanométrique flexoélectrique (deuxième matériau)

Longueur L (m) 400 × 10−9

Largeur b (m) 20 × 10−9
Épaisseur h (m) 20 × 10−9
Module de Young E = c11 (GPa) 3.7
Constantes diélectriques a11 (Vm/C) 1.38 × 1010
a33 (Vm/C) 1.38 × 1010
Constantes piézoélectriques d11 (V/m) 0
d31 (V/m) −1.02 × 109
Constantes flexoélectriques f11 (V) 0
f14 (V) 179.13
Échelle de longueur élastique l0 (m) 5.32 × 10−7
non-locale
Perméabilité  (C2 V-1 m-1) 8.142 × 10−11
Nombre d'éléments 40

26
Le résultat du matériau FGM composé de ces deux matériaux est obtenu dans
la Figure 6, et nous le comparons au résultat du premier matériau homogène.

Figure 6 Déflexion maximale en fonction de la charge (cas matériau homogène /cas FGM).

Figure 7 Le mode de déflexion sous tension de la nano—poutre ( cas matériau homogène /cas FGM).

27
1.3 Validation de FGM par rapport les matériaux

homogènes (voltage uniquement)

Nous considérons la poutre nanométrique à extrémité fixe avec les

caractéristiques présentées dans le Tableau 3, à une tension électrique v = 10V
.

La Figure 7 présente les réponses pour un seul matériau que nous avons validé
dans le paragraphe 3.2.

Figure 8 Déflexion maximale en fonction de la tension éléctrique.

Maintenant, nous utilisons le calcul pour les matériaux FGM composés de deux
matériaux supposés différents. Cependant, dans ce cas, nous utilisons deux
matériaux ayant des caractéristiques similaires, et nous devons constater que
les deux résultats sont très proches. Cela est illustré dans la Figure9.

28
 Figure 9 Déflexion maximale en fonction de la tension éléctrique (cas de flexo/cas FGM)

1.4 Étude de cas d'un matériau FGM composé de deux

matériaux différents (voltage uniquement)

Dans cette étude, nous nous intéressons à un matériau FGM (Matériau

Gradient Fonctionnel) constitué de deux matériaux distincts. L'objectif est
d'analyser le comportement mécanique de ce matériau composite et d'évaluer
ses performances par rapport à des applications spécifiques.

Les caractéristiques du premier matériau sont répertoriées dans le Tableau 5.

Les caractéristiques du deuxième matériau sont répertoriées dans le Tableau 6.

Le résultat du matériau FGM composé de ces deux matériaux est obtenu dans
la Figure 10, et nous le comparons au résultat du premier matériau 2eme
matériau homogène.

29
Tableau 5 Les caractéristiques de la micro-poutre flexoélectrique (1er matériau)

Longueur L (m) 100 × 10−6

Largeur b (m) 10 × 10−6
Épaisseur h (m) 1 × 10−6
Module de Young E = c11 (GPa) 131
Constantes diélectriques a11 (Vm/C) 0.8971 × 108
a33 (Vm/C) 0.788 × 108
Constantes piézoélectriques d11 (V/m) 0
d31 (V/m) 1.87 × 108
Constantes flexoélectriques f11 (V) 0
f14 (V) 5
Échelle de longueur élastique l0 (m) 0
non-locale
Perméabilité  (C2 V-1 m-1) 1.115 × 10−8
Nombre d'éléments 40
Tableau 6 Les caractéristiques de la micro-poutre flexoélectrique (2eme matériau)

Langueur L (m) 100 × 10−6

Largeur b (m) 10 × 10−6
Epaisseur h (m) 1 × 10−6
Module de Young E = c11 (GPa) 3.7
Constantes diélectriques a11 (Vm/C) 1.38 × 1010
a33 (Vm/C) 1.38 × 1010
Constantes d11 (V/m) 0
piézoélectriques d31 (V/m) −1.02 × 109
Constantes f11 (V) −220.47
flexoélectriques f14 (V) −179.13
Échelle de longueur l0 (m) 11.2 × 10−9
élastique non-locale
Perméabilité  (C2 V-1 m-1) 8.14 × 10−11
Nombre d'éléments 40

30
Figure 10 Déflexion maximale en fonction de la charge (cas1er matériau homogène /cas FGM/cas 2eme matériau
homogène ).

Figure 11 Le mode de déflexion sous tension de la nano—poutre ( cas 2em matériau homogène /cas FGM).

31
4.1 Effet flexoélectrique dans les nano-poutres FGM

Le nano-poutre FGM est composé de deux matériaux différents, avec une

distribution connue le long des axes x et z. On suppose que la distribution du
matériau est constante le long de l'axe y.

La longueur de la poutre est L, l'épaisseur est h et la largeur est b.

Les propriétés du matériau varient de telle sorte que, pour une propriété du

matériau  ( x, z ) la relation entre ses valeurs pour les deux constituants du

FGM, notées 1 et  2 , est donnée par :
nz
 x z−z
nx

 ( x, z ) = 1 + ( 2 − 1 )    
  
L h  (19)

où nx et nz sont des paramètres de matériau d'exposant de loi de puissance

le long des directions x et z , respectivement. z est la position de fibre

neutre par rapport à la face inférieure du poutre.

  x 
nx
x 
nx

h ( nz + 1)  c1  nz + 2 − 2    + 2c2   
   L    L  
  
z ( x) =
  x 
nx
x 
nx

2 ( nz + 2 )  c1  nz + 1 −    + c2    (20)
   L    L  
 

32
5 Conclusion et perspectives

33
Conclusion
En conclusion, cet article présente une technique novatrice pour le
développement d'un modèle numérique prenant en compte la non-linéarité
géométrique des poutres FGM-flexoélectriques. En utilisant la Méthode
Asymptotique Numérique associée à la méthode des éléments finis pour la
théorie d'Euler-Bernoulli, les auteurs parviennent à créer un modèle fiable et
précis.

L'article aborde également l'étude du comportement des matériaux FGM, en se

concentrant sur les matériaux solides, les cristaux liquides et les polymères.
Une validation rigoureuse du modèle est réalisée en comparant les résultats
obtenus avec des travaux antérieurs, ce qui confirme sa fiabilité.

De plus, l'analyse des comportements des matériaux FGM est approfondie à

travers des tests et des validations spécifiques. L'effet flexoélectrique dans les
nano-poutres FGM est étudié, ce qui permet de mieux comprendre les
propriétés électromécaniques de ces matériaux.

En résumé, cet article contribue à l'avancement de la recherche sur les poutres

FGM-flexoélectriques en proposant une méthode de développement de
modèle numérique prenant en compte la non-linéarité géométrique. Les
résultats obtenus et la validation réalisée démontrent la pertinence et la
fiabilité de cette approche. Ces avancées permettent d'approfondir notre
compréhension du comportement des matériaux FGM et ouvrent de nouvelles
perspectives pour leur utilisation dans des applications pratiques.

Perspectives

Ce projet ouvre plusieurs perspectives intéressantes pour la poursuite des

recherches dans le domaine des poutres FGM-flexoélectriques. Voici quelques-
unes des perspectives suggérées par les résultats et les conclusions de l'étude :

• Développement de modèles plus complets : Bien que cet article présente

une technique prometteuse pour modéliser les poutres FGM-
34
 flexoélectriques prenant en compte la non-linéarité géométrique, il est
possible d'explorer des modèles plus complexes et sophistiqués. Par
exemple, l'inclusion d'autres effets physiques tels que la
thermoélectricité ou les propriétés électromagnétiques pourrait
permettre une modélisation plus précise du comportement des poutres
FGM-flexoélectriques.
• Optimisation des propriétés des matériaux FGM-flexoélectriques : Les
matériaux FGM offrent une grande flexibilité dans la conception de
structures ayant des propriétés spécifiques. En utilisant le modèle
numérique développé dans cet article, il serait possible d'optimiser les
propriétés des matériaux FGM-flexoélectriques pour des applications
particulières. Par exemple, en ajustant la distribution des matériaux
constitutifs, il serait possible d'optimiser la réponse électromécanique
des poutres pour des performances maximales dans des domaines tels
que la génération d'énergie, les capteurs ou les actionneurs.
• Étude de comportement dynamique : L'analyse présentée dans cet
article se concentre principalement sur le comportement statique des
poutres FGM-flexoélectriques. Cependant, l'étude du comportement
dynamique, y compris les vibrations et les réponses dynamiques sous
différentes charges et excitations, est d'une importance capitale pour de
nombreuses applications pratiques. Les futures recherches pourraient
donc se concentrer sur l'extension de la méthode numérique proposée
pour inclure des analyses dynamiques.
• Expérimentation et validation : Bien que ce projet présente une
validation du modèle proposé en le comparant à des travaux antérieurs, il
serait également important de réaliser des expériences physiques pour
valider les résultats obtenus. L'expérimentation permettrait de confirmer
la validité du modèle et d'évaluer son applicabilité dans des situations
réelles.

35
6 Bibliographie

[Baroudi] STATIC AN DYNAMIC BEHAVIORS OF FLEXOELECTRIC PIEZOELECTRIC

NANOSTRUCTURES By Sourour BAROUDI Ph.D. Thesis submitted in partial
fulfillment of the requirements for the Degree o of Doctor in Applied Mechanics

36
Résumé
Ce projet est réalisé dans le cadre de mon projet de fin d’année pour l’année
académique 2022/2023 pour mes études d’ingenieurie en éléctromécanique à
l’acdemie militaire de fondouk jdid.

Ce projet présente une nouvelle technique innovante pour faciliter le

développement de modèles numériques pour les poutres FGM-
flexoélectriques.

La technique proposée repose sur l'utilisation de la méthode asymptotique

numérique combinée à la méthode des éléments finis pour la théorie d'Euler-
Bernoulli. La méthode asymptotique numérique simplifie les équations
complexes résultant de la non-linéarité géométrique en utilisant des
développements en séries asymptotiques pour obtenir des expressions
analytiques approximatives qui préservent les caractéristiques essentielles du
problème.

En utilisant la méthode des éléments finis, largement utilisée pour résoudre les
équations différentielles régissant le comportement des structures, cette
approche permet de modéliser avec précision le comportement non-linéaire
des poutres FGM-flexoélectriques en prenant en compte les effets de l'énergie
électrique et mécanique.

Mots clefs : Poutres FGM-flexoélectriques , MAN , Méthode des éléments

finis , Validation du modèle , Effets de l'énergie électrique et mécanique

37