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:
Stabilité et flambage des
poutres flexoélectriques

Élaboré par:
Gatri Malek

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. ”
- Malek

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 ............................................................................. 4

2.1 Histoire et évolutions ............................................................................. 5

2.2 Les poutres flexoelectriques .................................................................. 6

2.3 L'effet flexoélectrique ............................................................................ 7

2.4 Types des matériaux Flexoélectriques .................................................... 8

2.5 Stabilité des poutres flexoélectriques .................................................... 8

2.6 MAN....................................................................................................... 9

2.6.1 MAN (définition) .............................................................................. 9

2.6.2 Le MAN pour le calcul des stabilités des poutres flexoeléctrique ... 10

2.6.3 L’utilisation de la MAN dans ce rapport.......................................... 12

3 Formulation mathématique des poutres flexoélectriques .......................13

3.1 Choix de la forme de la structure ......................................................... 14

3.2 Mise en équation du système .............................................................. 14

3.3 La procédure de calcul des termes MAN .............................................. 19

3.4 Validation de la MAN pour les poutres élastiques avec non-linéarité

géométrique ................................................................................................. 23

4 Stabilité des poutres flexoélectriques .....................................................25

4.1 MAN pour les poutres flexoélectriques avec non-linéarité géométrique

…………………………………………………………………………………………………………..26

III
 4.1.1 Tests préliminaires ......................................................................... 26

4.1.2 Modélisation du flambage des poutres flexoélectriques actionneurs

sous tension .............................................................................................. 27

4.1.3 Modélisation de la variation de tension dans les capteurs à poutre

flexoélectrique soumis à une force externe ............................................... 29

5 Conclusion et perspectives......................................................................32

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. ............................ 23
Figure 3 Flambage de la poutre élastique du Tableau 2. .................................. 24
Figure 4 Comparaison avec les résultats de Baroudi et al. [...]. ....................... 27
Figure 5 Flambage de la poutre flexoélectrique PVDF sous tension................ 28
Figure 6 Le mode de flambage sous tension de la micro-poutre piezoélectrique
en PVDF (Tableau 4). ........................................................................................ 28
Figure 7 Flambage de la poutre flexoélectrique en BaTiO3 sous tension.......... 29
Figure 8 Le comportement non linéaire de la nano-poutre flexoélectrique en
PVDF (Tableau 5) .............................................................................................. 30
Figure 9 Tension en fonction de la charge pour la nano-poutre flexoélectrique
en PVDF du Tableau 5 ...................................................................................... 31
Figure 10 Tension en fonction de la charge pour la nano-poutre flexoélectrique
en BaTiO3 du Tableau 3 ................................................................................... 31

V
Liste des tableaux

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. […] :23
Tableau 2 Propriétés de la poutre élastique utilisées pour la validation des
tests de flambage : ........................................................................................... 24
Tableau 3 Les caractéristiques de la poutre nanométrique flexoélectrique
étudiée par Yan et al. [...] ................................................................................. 26
Tableau 4 Les caractéristiques de la micro-poutre flexoélectrique en PVDF. . 27
Tableau 5 Les caractéristiques des nano-poutres flexoélectriques en PVDF ... 30

VI
1 Introduction générale

Dans cet article, nous présentons une technique novatrice qui facilite le
développement de modèles numériques pour les poutres flexoélectriques en
prenant en compte la non-linéarité géométrique. Notre approche 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.

Les poutres flexoélectriques sont des structures composées de matériaux

piezoélectriques, capables de convertir l'énergie mécanique en énergie
électrique et vice versa. La non-linéarité géométrique est un aspect crucial à
prendre en compte dans la modélisation de ces poutres, car elle affecte leur
comportement mécanique et électrique.

La méthode asymptotique numérique est une approche mathématique

puissante qui permet de simplifier les équations complexes résultant de la non-
linéarité géométrique. Elle se base sur des développements en séries
asymptotiques pour obtenir des expressions analytiques approximatives, tout
en préservant les caractéristiques essentielles du problème.

Nous mettons en œuvre cette technique dans le cadre de la méthode des

éléments finis, qui est une méthode numérique couramment utilisée pour
résoudre les équations différentielles régissant le comportement des
structures. Notre approche permet de modéliser avec précision le
comportement non-linéaire des poutres flexoélectriques, en prenant en
compte les effets de l'énergie électrique et mécanique.

Pour valider notre modèle, nous le comparons à des travaux antérieurs réalisés
dans ce domaine. Cette validation démontre la fiabilité et l'efficacité de notre
approche pour la modélisation des poutres flexoélectriques prenant en compte
la non-linéarité géométrique.

En outre, nous étudions le comportement électrique des poutres

flexoélectriques et mettons en évidence un phénomène important : la charge
électrique peut provoquer un flambage de la poutre. Cette découverte souligne
1
l'importance de considérer les interactions entre les aspects électriques et
mécaniques dans la conception et l'analyse des poutres flexoélectriques.

En conclusion, notre technique basée sur la méthode asymptotique numérique

combinée à la méthode des éléments finis offre une approche précise et
efficace pour modéliser les poutres flexoélectriques en prenant en compte la
non-linéarité géométrique. Notre étude met en évidence le comportement
électrique de ces structures et souligne l'importance de considérer les effets
électromécaniques dans leur analyse. Ces résultats peuvent contribuer à
améliorer la conception et la compréhension des poutres flexoélectriques,
ouvrant ainsi de nouvelles perspectives dans différents domaines d'application.

Le premier chapitre “État de l’art et histoire ” dans lequel nous présenterons

les principaux aspects couverts par cet article, en commençant par retracer
l'histoire et les évolutions de la flexoélectricité. Ensuite, nous nous pencherons
sur les poutres flexoélectriques elles-mêmes, en explorant leur structure et
leurs caractéristiques essentielles. L'effet flexoélectrique, qui est au cœur de ce
phénomène, sera également discuté en détail. Nous examinerons ensuite les
différents types de matériaux flexoélectriques, en mettant en évidence leurs
propriétés et leurs applications potentielles. L'un des aspects importants des
poutres flexoélectriques est leur stabilité, et nous discuterons des facteurs qui
influencent cette stabilité et des défis associés. Enfin, nous aborderons le
concept de MAN automatique (Méthode Asymptotique Numérique) et son rôle
dans le développement de modèles numériques pour les poutres
flexoélectriques.

Le deuxième chapitre “Formulation mathématique des poutres

flexoélectriques“ dans lequel Nous commençons par discuter du choix de la
forme de la structure, qui est une étape cruciale dans la modélisation des
poutres flexoélectriques . Ensuite, nous abordons la mise en équation du
système, c'est-à-dire comment décrire mathématiquement le comportement
d'une poutre flexoélectrique. Nous examinons les équations qui régissent les
interactions entre les déformations mécaniques et les charges électriques, en
prenant en compte la non-linéarité géométrique. Nous présentons également
la procédure automatique de calcul des termes de la Méthode Asymptotique
Numérique (MAN). Cette méthode, utilisée pour simplifier les équations
complexes, joue un rôle essentiel dans le développement d'un modèle
2
numérique précis pour les poutres flexoélectriques. Enfin, nous nous
intéressons à la validation de la MAN automatique pour les poutres élastiques
avec non-linéarité géométrique. Nous examinons des études antérieures et
comparons les résultats obtenus avec ceux de notre modèle, afin de démontrer
la fiabilité et l'efficacité de notre approche.

Le troisième chapitre “Stabilité des poutres flexoélectriques“ nous présentons

les tests préliminaires réalisés dans le cadre de la MAN automatique. Ces tests
servent à évaluer la fiabilité et l'efficacité de la MAN pour modéliser les poutres
flexoélectriques avec non-linéarité géométrique. Nous abordons également la
modélisation du flambage des poutres flexoélectriques actionneurs sous
tension. Le flambage est un phénomène critique qui peut compromettre la
stabilité et les performances des poutres flexoélectriques, et nous examinons
comment la MAN automatique permet de prendre en compte cet aspect dans
la modélisation. Enfin, nous discutons de la modélisation de la variation de
tension dans les capteurs à poutre flexoélectrique soumis à une force externe.
Cette modélisation est essentielle pour comprendre comment les capteurs
réagissent aux forces externes et comment la tension induite peut influencer la
stabilité de la poutre flexoélectrique.

3
Chapitre 1

2 Etat de l’art et histoire

4
2.1 Histoire et évolutions

La flexoélectricité est un phénomène qui se produit dans certains matériaux

lorsqu'ils sont soumis à une déformation mécanique, telle qu'une flexion ou
une torsion. Ce phénomène est lié à l'interaction entre la déformation
mécanique et les charges électriques présentes dans le matériau, ce qui
entraîne l'apparition d'une polarisation électrique.

L'histoire de la flexoélectricité remonte aux travaux du physicien français Marie

Curie [1], qui a découvert en 1880 l'effet piézoélectrique, un phénomène
similaire dans lequel une polarisation électrique se produit en réponse à une
contrainte mécanique. Au fil des années, des chercheurs ont exploré davantage
les propriétés piézoélectriques des matériaux, mais ce n'est qu'à partir des
années 1990 que le concept de flexoélectricité a commencé à être étudié en
profondeur.

L'évolution de la flexoélectricité a été étroitement liée aux avancées dans la

compréhension des matériaux et des phénomènes physiques qui la sous-
tendent. Des recherches ont été menées pour identifier les matériaux
présentant des propriétés flexoélectriques significatives, ainsi que pour
développer des techniques de mesure et de caractérisation de ces propriétés.

Plus récemment, des progrès ont été réalisés dans la modélisation et la

simulation de la flexoélectricité, ce qui a permis une meilleure compréhension
des mécanismes à l'œuvre et ouvert la voie au développement de nouvelles
applications. Par exemple, la flexoélectricité a été étudiée dans le contexte de la
conception de dispositifs microélectromécaniques (MEMS) et de capteurs
ultrasensibles.

En outre, la flexoélectricité a également suscité un intérêt croissant dans le

domaine des matériaux fonctionnels et des nanotechnologies. Les chercheurs
explorent les possibilités d'utiliser des matériaux flexoélectriques pour
développer des dispositifs électroniques plus efficaces et plus compacts, ainsi
que pour créer des systèmes de conversion d'énergie et des capteurs innovants.

En résumé, l'histoire de la flexoélectricité remonte aux travaux pionniers de

Marie Curie sur le phénomène piézoélectrique. Au fil du temps, la
5
compréhension de la flexoélectricité s'est approfondie grâce à des progrès dans
la modélisation, la caractérisation des matériaux et les applications potentielles.
Aujourd'hui, la flexoélectricité continue d'être un domaine de recherche actif,
offrant des opportunités pour l'innovation dans différents domaines
technologiques.

2.2 Les poutres flexoelectriques

Les poutres flexoélectriques sont des structures qui possèdent des propriétés
flexoélectriques, c'est-à-dire qu'elles peuvent générer une charge électrique
lorsqu'elles sont soumises à une contrainte mécanique, et vice versa. Ces
matériaux sont souvent utilisés dans le domaine de l'électronique et de la
robotique pour créer des capteurs et des actionneurs sensibles à la flexion.

L'effet flexoélectrique est similaire à l'effet piézoélectrique, où une charge

électrique est générée en réponse à une contrainte mécanique. Cependant,
dans le cas du flexoélectrique, la charge électrique est produite en raison d'un
gradient de contrainte à travers la section de la poutre, plutôt que d'une
contrainte uniforme.

Les poutres flexoélectriques peuvent être utilisées pour mesurer des forces, des
pressions ou des déformations mécaniques. Lorsqu'elles sont pliées ou
déformées, elles génèrent une charge électrique proportionnelle à la contrainte
appliquée. Cette charge électrique peut ensuite être mesurée et utilisée pour
détecter ou contrôler des mouvements, des déformations ou d'autres
paramètres physiques.

En raison de leurs propriétés flexoélectriques, les poutres flexoélectriques

présentent plusieurs avantages, tels qu'une grande sensibilité, une réponse
rapide et une faible consommation d'énergie. Elles peuvent être fabriquées à
partir de différents matériaux, tels que des polymères, des céramiques ou des
composites, en fonction des exigences spécifiques de l'application.
6
En résumé, les poutres flexoélectriques sont des structures qui convertissent les
contraintes mécaniques en signaux électriques grâce à l'effet flexoélectrique.
Elles trouvent des applications dans des domaines tels que l'électronique, la
robotique, les capteurs de force et les actionneurs.

2.3 L'effet flexoélectrique

La flexoélectricité fait référence au couplage entre la polarisation et les

gradients de contrainte, plutôt qu'entre la polarisation et la contrainte comme
c'est le cas pour la piézoélectricité. L'origine du nom provient du mot latin
"flexus" qui signifie plier, et cela est interprété par le fait qu'un gradient de
contrainte apparaît normalement dans bent plates[2]. Mashkevich [3] et
Tolpygo [4] ont été les premiers à prouver l'existence de l'effet flexoélectrique
théoriquement depuis 1957. En 1964, Kogan [5] a décrit cet effet de manière
phénoménologique en interprétant le contexte du couplage électron-phonon
dans les cristaux centro-symétriques, pour lesquels le couplage flexoélectrique
a une influence importante. Un autre cadre phénoménologique décrivant l'effet
a été élaboré par Mindlin [6].

Étant donné que de grands gradients de contrainte sont facilement générés à

l'échelle nanométrique, la flexoélectricité est considérée comme non
négligeable pour les dispositifs basés sur le nanomètre. Ainsi, étant donné que
cet effet était censé être faible à grande échelle, l'intérêt pour ce phénomène
n'a pas été très important au cours des dernières décennies. Récemment, avec
le développement de la nanotechnologie, la flexoélectricité a suscité de plus en
plus l'attention de la communauté scientifique et de recherche.

Pour mieux comprendre le phénomène de flexoélectricité, on peut se référer à

plusieurs articles tels que [2], où une revue complète de la théorie de la
flexoélectricité est développée, ainsi qu'un résumé concis plus récent traité par
Maranganti et al.[7]. De plus, Zubko et al.[8] ont passé en revue dans leur
article les fondements de l'effet flexoélectrique dans les solides et ont discuté
de sa présence dans de nombreux systèmes à l'échelle nanométrique.
7
2.4 Types des matériaux Flexoélectriques

• Les cristaux liquides

• Les membranes biologiques
• Les polymères

2.5 Stabilité des poutres flexoélectriques

La stabilité des poutres flexoélectriques dépend de plusieurs facteurs,

notamment la géométrie de la poutre, les propriétés mécaniques du matériau
flexoélectrique et les charges appliquées. Voici quelques considérations
importantes liées à la stabilité des poutres flexoélectriques :

Flambement : Le flambement est un phénomène de déformation élastique

instable qui se produit lorsque la poutre est soumise à une charge de
compression axiale. La flexoélectricité peut influencer le comportement de
flambement en générant une charge électrique lors de la déformation. Pour
prévenir le flambement, il est nécessaire de prendre en compte les propriétés
mécaniques du matériau et de concevoir la poutre avec des dimensions
appropriées.

Instabilités dynamiques : Lorsque les poutres flexoélectriques sont soumises à

des charges dynamiques, des instabilités peuvent se produire en raison de la
rétroaction entre la déformation mécanique et la génération de charge
électrique. Cela peut entraîner des vibrations indésirables ou des
comportements instables de la poutre. La conception et la modélisation
précises du système sont essentielles pour éviter de telles instabilités.

Effets thermiques : Les variations de température peuvent affecter les

propriétés mécaniques du matériau flexoélectrique, ce qui peut influencer la
stabilité de la poutre. Des contraintes thermiques peuvent être induites lors de
l'activation du matériau flexoélectrique, ce qui peut entraîner des déformations
non désirées ou des contraintes excessives. Une gestion appropriée des effets

8
thermiques est donc nécessaire pour assurer la stabilité des poutres
flexoélectriques.

Effets électriques : Les charges électriques générées par l'effet flexoélectrique

peuvent également influencer la stabilité de la poutre. Des phénomènes tels
que l'effet de contrainte, l'effet piézoélectrique et l'effet électrostrictif peuvent
interagir avec la déformation mécanique, ce qui peut affecter la réponse
globale de la poutre. La modélisation précise du couplage électromécanique est
essentielle pour prédire et gérer ces effets électriques.

2.6 MAN
2.6.1 MAN (définition)

La méthode asymptotique numérique est une approche mathématique utilisée

pour approximer les solutions de problèmes complexes en exploitant les
comportements asymptotiques des variables ou des paramètres du système.
Elle combine des techniques numériques et des analyses asymptotiques pour
obtenir des estimations précises et efficaces des solutions.

La méthode asymptotique numérique repose sur les principes suivants :

Identification des régimes asymptotiques : On détermine 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 associés à des caractéristiques
particulières du système et peuvent permettre de simplifier les équations ou de
négliger certains termes.

Décomposition en termes asymptotiques : On décompose les équations du

système en une série de termes asymptotiques, où 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.

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

est résolu individuellement à l'aide de méthodes numériques adaptées à la

9
nature du terme. Cela peut inclure 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 effectuée en utilisant des techniques de
sommation ou d'assemblage des termes asymptotiques.

La méthode asymptotique numérique est largement utilisée dans divers

domaines scientifiques et techniques, tels que la physique, les mathématiques
appliquées, l'ingénierie, la mécanique des fluides, etc. Elle permet de traiter des
problèmes complexes et d'obtenir des résultats précis tout en réduisant la
complexité des calculs. Cependant, il est important de noter que l'utilisation de
cette méthode nécessite une compréhension approfondie des propriétés
asymptotiques du système étudié et une validation rigoureuse des
approximations effectuées.

2.6.2 Le MAN pour le calcul des stabilités des poutres

flexoeléctrique

La Méthode Numérique Asymptotique (MAN) présente une efficacité notable

pour résoudre une grande variété de problèmes non linéaires, tels que la
dynamique structurelle [9], la dynamique des fluides [10], le flambage et la
stabilité structurels . La MAN offre également la possibilité de résoudre
plusieurs problèmes industriels, notamment la modélisation des défauts de
planéité dans le laminage des bandes . Elle consiste à approximer la solution
non linéaire sous forme de série polynomiale ou de série PADE [11], ce qui
présente les principaux avantages suivants :

Avantage de la MAN

Calcul automatique du pas

À chaque étape, le MAN trouve une solution en construisant soigneusement

des séries de Taylor pour calculer les inconnues. On s'attaque à un problème
10
qui existe depuis le début, et on essaie d'y arriver avec des critères précis. On a
l'impression de se rapprocher d'une solution, même si c'est encore un peu flou.
Trouver une réponse non linéaire n'est pas un jeu d'enfant. Cela nécessite de
reconstituer une série de courbes polynomiales. Pourtant, même si vous
pouvez en faire une approximation, cela ne signifie pas que vous avez déchiffré
le code. La notion de travailler avec des fractions rationnelles n'est pas
seulement une théorie ; c'est une réalité grâce aux 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 calculs complexes.
L'utilisateur a la responsabilité de décider au préalable d'une commande, et
plusieurs auteurs ont mis en évidence qu'une sélection appropriée de ce type
de séquence peut être très bénéfique. Le nombre idéal d'activités, selon
Boutyour et al., semble se situer entre quinze et vingt. Cette gamme est assez
pratique, si vous me demandez. L'amplitude du pas est précisément le rayon de
convergence de cette série particulière. C'est un fait connu que cette distance
varie selon les différentes séquences numériques.

Calcul direct de la solution à chaque pas

Sur chaque pas, les méthodes PC consistent à effectuer une succession de

prédictions linéaires du problème (estimation du (des) champ(s) solution(s)),
calculer le résidu correspondant (erreur) et corriger par la suite. Une prédiction-
correction correspond à une itération et sur chaque pas un calcul normal peut
effectuer entre 3 et 5 itérations. Dans un code (EF), une correction revient à
multiplier le résidu par l’inverse de la matrice dite tangente qui, par analogie,
représente la pente du problème non linéaire à chaque valeur considérée des
inconnues. Les inversions des matrices tangentes rendent les calculs couteux
(en temps CPU) avec les méthodes itératives, alors qu’un calcul MAN donne
l’avantage d’un calcul direct de la solution sur un pas.

Inconvénient de la MAN

Le problème tangent utilisé pour la prédiction dans les algorithmes itératifs

n’est rien autre que le premier terme de la série de Taylor de la MAN. C’est pour
cette raison que l’on considère que la MAN est une généralisation des
méthodes PC. Les problèmes à ordres supérieurs engendrés par le
développement en série apparaissent à cause de la non-linéarité. Ils sont

11
généralement à déterminer analytiquement et à être implémentés dans les
codes de calculs. Bien que la MAN offre plusieurs 3 avantages pour la
modélisation des problèmes non linéaires, elle présente néanmoins une
lourdeur dans le développement de ces termes, ce qui est l’inconvénient
principal de la MAN. Rien que dans un calcul de poutres élastiques en grands
déplacements (non-linéarité géométrique), le développement analytique des
termes d’ordres supérieurs est fastidieux. Cet inconvénient est de plus en plus
accentué lors de la résolution d’un problème avec des couplages multi
physiques (structures flexoélectiques par exemple) d’où l’objectif de ce projet :
automatisation de la MAN.

2.6.3 L’utilisation de la MAN dans ce rapport

Dans cet article, nous exposons une technique qui facilite le développement
d'un modèle numérique pour les poutres 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 électrique de la poutre et montrons que la
charge électrique peut provoquer un flambage.

12
Chapitre 2

3 Formulation mathématique des poutres

flexoélectriques

13
3.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.

3.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 :

W L
Tij =
 ij
W L (4)
T ijk =
 ij , k
15
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 = −
a33
(6)

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

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

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

Où u est le déplacement longitudinal, w est la flèche dans l'axe neutre de la

 
poutre et  est la courbure. Les notations ( g) et ( g) représentent
respectivement le premier et le second gradient des champs de déplacement
par rapport à la direction longitudinale x . v est la tension à travers l'électrode et
 ( x ) est l'expression de la distribution longitudinale du potentiel électrique. La

relation entre  and w  dans le système (8) est considérée comme une
contrainte du problème pour laquelle la méthode de pénalité est utilisée. En
combinant les équations (3), (5), (6), (7) et (8), nous obtenons l'équation
suivante :

17
 L
 d2 d2   w2  f11d11
−   c11 − 11 − 13  A  u  + − A ( u  + ww ) +
0
 a11 a33   2  a 11

f14 d13 d 2A d A   w2 

Aw − 11   ( x ) − 13 v    u  + −
a33 a11  a33 h   2 
 f11d11  w2   f112 
 
A u +  +  g11 −  A ( u  + ww ) −
 11
a  2   a11 

f11 2 A 
  ( x )   ( u  + ww ) +
a11  
  d112 d132   f142   d11 f11
  11
c − − +
  11
I g −  A w − I  +

  a11 a33   a33    a 11

d13 f14  w2  d13 2 A f A 

A  u + −  ( t ) + 14 v   w −
a33  2  a33  a33 h 

 f11d11  f112  
 I w +  11
g −  I     −
 a11  a11  
 2 d11  w2  2 f11
 A  u  + − A ( u  + ww ) −
  a11  2   a11
1  1  
A +     ( x )    ( x ) −
2  a11  
 2 d13 2 A  1  
 Aw −  +    ( x )   ( x ) −
  a33 2 h 2  a33  
 A d13  w2  A f14 A 1  
 
u + + w − 2  +   v v +
 h a33  2  h a33 h  a33  
 (  − w )  (  − w ) dx +  Wnc = 0
(8)

Où  est le paramètre de pénalité et g11 est lié à l'échelle de longueur élastique

non locale l0 comme suit :

18
 g11
l0 =
c11 (9)

La résolution de l'équation (9) avec la méthode MAN est délicate car elle
présente un grand nombre de variables. En effet, le risque d'erreurs dans les
développements mathématiques

et la programmation est critique. Par conséquent, une procédure générique qui

rend le calcul des termes MAN automatique est nécessaire.

3.3 La procédure 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 (11)
n
 = 0 +  a p  p
p =1

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:
19
 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 (12)

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

 (13)
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
(14)
w = H T wˆ e

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 :
20
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 (15)
  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

  wˆ ep
e

 L L
 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 +
2

 (16)
  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ù :

21
 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 (17)

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

 (18)
 K ww wˆ 1 + K wu uˆ1 = 1 F

w

Order p  1:

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

Où :

K uu = U K uue , K ww = U K ww
e
, K uw = K wu
T
= U K uw
e
, Fnlwp = U Fnlewp (11)
e e e e

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

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

Le modèle avec MAN 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. [13]
(voir Figure 2). Nous avons également constaté un accord avec les résultats du
développement de la MAN 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. [13] :

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.

23
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.

24
Chapitre 3

4 Stabilité des poutres flexoélectriques

25
4.1 MAN pour les poutres flexoélectriques avec non-
linéarité géométrique

4.1.1 Tests préliminaires

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 = −1nN et à une tension électrique v = −0.1V . La Figure 4 montre que les
réponses proposées par Baroudi et al. [16] sont proches de nos résultats avec et
sans flexoélectricité (uniquement l'élasticité).
Tableau 3 Les caractéristiques de la poutre nanométrique flexoélectrique étudiée par Yan et al. [13]

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

26
 Figure 4 Comparaison avec les résultats de Baroudi et al. [14].

4.1.2 Modélisation du flambage des poutres

flexoélectriques actionneurs sous tension

Une poutre flexoélectrique est un actionneur lorsque sa déformation est

contrôlée par une tension imposée. Nous considérons une poutre
flexoélectrique en PVDF avec une extrémité fixe soumise à une tension variable
v sans charge mécanique appliquée ( F = 0N ). Les caractéristiques
matérielles et géométriques sont présentées dans le Tableau 4.
Tableau 4 Les caractéristiques de la micro-poutre flexoélectrique en PVDF.

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

27
élastique non-locale
Perméabilité  (C2 V-1 m-1) 8.14 × 10−11
Nombre d'éléments 40

La Figure 5-a montre la variation de la déflexion maximale wmax lorsque la

tension v augmente. On observe que la poutre fléchit et qu'un point de
bifurcation est observé pour une tension v proche de 16 V, avec le mode
présenté dans la Figure 6. Dans ce cas, v est négatif et avec une tension
positive, le fléchissement ne se produit pas. En revanche, pour la poutre en
BaTio3 de la Table 3 dans des conditions similaires, le fléchissement est généré
par une tension positive, comme illustré dans la Figure 7. Cette différence est
liée au signe des paramètres piézoélectriques aij et f ij .

Figure 5 Flambage de la poutre flexoélectrique PVDF sous tension.

Figure 6 Le mode de flambage sous tension de la micro-poutre piezoélectrique en PVDF (Tableau 4).

28
 Figure 7 Flambage de la poutre flexoélectrique en BaTiO3 sous tension.

4.1.3 Modélisation de la variation de tension dans les

capteurs à poutre flexoélectrique soumis à une force
externe

Une poutre flexoélectrique est utilisée comme capteur lorsqu'aucune tension

v n'est imposée. Cependant, un potentiel électrique est révélé lorsque la
poutre est déformée sous une charge mécanique. Une nano-poutre
flexoélectrique en PVDF présentée dans le Tableau 5 est considérée dans une
telle situation. Elle est soumise à une force verticale variable F sans tension
v imposée, et son comportement non linéaire est tracé dans la Figure 8. La
déformation génère un potentiel électrique qui correspond à l'évolution de la
tension exposée dans la Figure 9. Des conditions similaires sont utilisées pour la
nano-poutre en BaTiO3 du Tableau 3. Le comportement est comparable, et
nous traçons l'évolution de la tension en fonction de la force F dans la Figure
10.

29
Tableau 5 Les caractéristiques des nano-poutres flexoélectriques en PVDF

Longueur L (m) 100 × 10−9

Largeur b (m) 10 × 10−9
Epaisseur h (m) 1 × 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 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

Dans la figure 8 w est entre 0 et -2.75 nm

Figure 8 Le comportement non linéaire de la nano-poutre flexoélectrique en PVDF (Tableau 5)

30
 Figure 9 Tension en fonction de la charge pour la nano-poutre flexoélectrique en PVDF (Tableau 5)

Figure 10 Tension en fonction de la charge pour la nano-poutre flexoélectrique en BaTiO3 du (Tableau 3)

31
5 Conclusion et perspectives

32
Conclusion générale

En conclusion, ce projet a abordé plusieurs aspects clés liés à la formulation

mathématique et à la stabilité des poutres flexoélectriques, en prenant en
compte la non-linéarité géométrique. Nous avons exploré la méthode
automatique de calcul de la Méthode Asymptotique Numérique (MAN) et son
rôle dans la modélisation précise des poutres flexoélectriques.

Nous avons souligné l'importance du choix de la forme de la structure et de la

mise en équation du système pour décrire mathématiquement le
comportement des poutres flexoélectriques. La MAN automatique s'est avérée
être une approche efficace pour simplifier les équations complexes résultant de
la non-linéarité géométrique, permettant ainsi le développement de modèles
numériques fiables.

Nous avons examiné les tests préliminaires qui ont permis de valider la MAN
automatique pour les poutres flexoélectriques avec non-linéarité géométrique.
Ces tests ont démontré la robustesse et l'efficacité de la méthode, renforçant
ainsi sa crédibilité pour l'étude du comportement des poutres flexoélectriques.

De plus, nous avons discuté de la modélisation du flambage des poutres

flexoélectriques actionneurs sous tension, mettant en évidence les effets
critiques de ce phénomène sur la stabilité des poutres. Nous avons également
exploré la modélisation de la variation de tension dans les capteurs à poutre
flexoélectrique soumis à une force externe, soulignant l'importance de
comprendre l'influence de la tension induite sur la stabilité de la poutre.

En résumé, ce projet fournit une approche complète et approfondie de la

formulation mathématique et de la stabilité des poutres flexoélectriques avec
non-linéarité géométrique. Les méthodes présentées offrent des outils précieux
pour la modélisation et l'analyse des poutres flexoélectriques, contribuant ainsi
à une meilleure compréhension de leur comportement mécanique et
électrique.

Ces résultats peuvent avoir des implications significatives dans divers domaines
d'application, tels que la conception de dispositifs flexoélectriques avancés, les

33
capteurs intelligents et les actionneurs précis. Ils ouvrent également la voie à de
nouvelles perspectives de recherche pour améliorer les performances et
l'efficacité des poutres flexoélectriques.

En fin de compte, les avancées présentées dans ce projet permettent de

promouvoir le développement de nouvelles technologies basées sur les poutres
flexoélectriques, ouvrant ainsi la voie à des applications innovantes et
impactantes dans divers domaines de l'ingénierie et de la science des
matériaux.

Perspectives

Cet article sur la formulation mathématique et la stabilité des poutres

flexoélectriques offre des perspectives intéressantes pour des recherches
futures. Voici quelques domaines clés qui pourraient bénéficier de
développements ultérieurs :

• Modèles avancés de non-linéarité géométrique : Bien que cet article ait

pris en compte la non-linéarité géométrique dans la modélisation des
poutres flexoélectriques, des modèles plus complexes et plus précis
pourraient être développés. Des recherches supplémentaires pourraient
se concentrer sur l'exploration de formulations mathématiques avancées
pour capturer des effets géométriques plus complexes, tels que la torsion
et la courbure plus prononcées.
• Analyse multiphysique : Les poutres flexoélectriques combinent des
comportements mécaniques et électriques, et une perspective
prometteuse consisterait à développer des modèles d'analyse
multiphysique plus complets. Cela permettrait de mieux comprendre
l'interaction entre les phénomènes électromécaniques et de prédire avec
précision les performances des poutres flexoélectriques dans des
conditions réelles.
• Optimisation des structures : Les méthodes de conception et
d'optimisation des poutres flexoélectriques pourraient être étudiées plus

34
 en profondeur. Des approches basées sur des algorithmes d'optimisation
pourraient être développées pour maximiser les performances des
poutres, en prenant en compte des critères tels que la stabilité, la
sensibilité électrique et la résistance mécanique.
• Applications émergentes : Les poutres flexoélectriques offrent un large
éventail d'applications potentielles, allant des capteurs et actionneurs
aux dispositifs de conversion d'énergie. Des recherches futures
pourraient se concentrer sur l'exploration de nouvelles applications
émergentes, telles que l'utilisation des poutres flexoélectriques dans les
domaines de la robotique souple, de la microélectronique et de la
médecine.

En intégrant ces perspectives de recherche dans les développements futurs, il

est possible de faire progresser davantage notre compréhension des poutres
flexoélectriques, de stimuler l'innovation technologique et d'explorer de
nouvelles frontières dans ce domaine prometteur.

35
6 Bibliographie

[1] Marie Curie, physicienne et chimiste polonaise naturalisée française, a

révolutionné la science en découvrant la radioactivité et en devenant la
première femme à remporter deux prix Nobel.

[2] P. V. Yudin et A. K. Tagantsev. Fundamentals of flexoelectricity in solids.

Nanotechnology, 24:432001, 2013.

[3] VS Mashkevich. Electrical, optical, and elastic properties of diamond-type

crystals ii. lattice vibrations with calculation of atomic dipole moments. Soviet
Physics JETP, 5 (4), 1957.

[4] KB Tolpygo. Long wavelength oscillations of diamond-type crystals including

long range forces. Soviet Physics-Solid State, 4(7):1297–1305, 1963.

[5] Sh M Kogan. Piezoelectric effect during inhomogeneous deformation and

acoustic scattering of carriers in crystals. Soviet Physics-Solid State, 5(10):2069–
2070, 1964.

[6] Raymond David Mindlin. Polarization gradient in elastic dielectrics.

International Journal of Solids and Structures, 4(6):637–642, 1968.

[7] R. Maranganti, N.D. Sharma, and P. Sharma. Electromechanical coupling in

nonpiezoelectric materials due to nanoscale nonlocal size effects: Green’s
function solutions and embedded inclusions. Physical Review B, 74(1):014110,
2006.

[8] Pavlo Zubko, Gustau Catalan, and Alexander K. Tagantsev. Flexoelectric effect
in solids. Annual Review of Materials Research, 43:387–421, 2013.

[9] La dynamique structurelle étudie les mouvements et les réponses des

structures aux charges dynamiques, fournissant des informations essentielles
pour leur conception et leur résistance.

[10] La dynamique des fluides concerne l'étude du mouvement et du

comportement des fluides, tels que les liquides et les gaz, ainsi que des forces
qui les influencent.

36
[11] La série de Padé est une méthode mathématique permettant d'approximer
une fonction complexe par un quotient de polynômes, offrant une meilleure
précision que la série de Taylor traditionnelle.

[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

37
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 flexoélectriques en tenant compte de la non-linéarité géométrique.

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 flexoélectriques en prenant en compte les effets de l'énergie
électrique et mécanique.

Le projet met également en évidence un phénomène important : la charge

électrique peut provoquer un flambage de la poutre flexoélectrique. Cela
souligne l'importance de considérer les interactions entre les aspects
électriques et mécaniques dans la conception et l'analyse des poutres
flexoélectriques.

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

Validation du modèle , Stabilété , Flambage

38