Untitled 2

Chapitre 6 : Approximation, Fonctions
d
d'Interpolation
Interpolation unidimensionn
onnelle de type Lagrange, Interpolation polynô
ômiale :
Fonctions de formes, Polynôm me de Lagrange, Polynôme d'Hermite, Trian ngle de
Pascal, Conditions de conformit
mité.
Fonctions d'interpolation
Les fonctions de forme ou fon

onctions d’interpolation sont les fonctions Ni qui
ui relient
les déplacements d’un point ququelconque intérieur à un élément aux n déplace acements
nodaux qi qui sont les degrés dee liberté dans le cas de l’approche cinématiquee : il y a
pour un élément autant de foncti
tions de forme que de degrés de liberté dans l’éllément.
Elles assurent le passage du p

problème continu au problème discret, la connaaissance
du déplacement en quelques nœ œuds discrets permettant de reconstruire le ch
hamp de
déplacement dans l’élément. Lee déplacement en un point quelconque de l’élém ment est
une combinaison linéaire des déplacements nodaux, dont les coefficients sont son les
valeurs des fonctions de forme een ce point.
Règles à respecter pour la cconstruction des fonction de forme:
1- Principe de compatibilité (de continuité):
La construction de la solutition approchée sur un élément fini doit essentieellement

satisfaire les exigences du prob oblème à résoudre et la géométrie de l'élémen nt. Pour
illustrer cette affirmation, considdérons les problèmes de barre et de poutre repré
résentés
respectivement sur les figures 7.47.4a et b. Sous l'effet de la force appliquée P, chaque
section transversale A de la barre
arre est soumise à une contrainte constante = P/A.P En
conséquence, la barre est sou soumise à une déformation constante = /E, /E où E
représente le module d'élasticitéé du matériau.
Dans un contexte unidimens nsionnel, la déformation normale est en fait donnée

comme une dérivée directe du d déplacement u(x) ; c'est-à-dire = du(x)/dx. Pu uisque la
déformation est constante sur tou
oute la barre, il s'ensuit que le déplacement u(x)
x) est une
fonction linéaire de x. En cons
onséquence, il est possible de construire une fonction
f
approchée u(x) pour le déplacemment en utilisant un polynôme linéaire:
u(x) = 1 + 2x
65
Les paramètres 1 et 2 sont identifiés à l'aide des deux valeurs nodales d'extrémité
U1 et U2. Le problème des barres est classé comme un problème C0. La solution
approchée doit être continue et sa dérivée doit exister.
Considérons maintenant le problème du poutre. Sous la charge appliquée

uniformément répartie, chaque section transversale de la poutre est soumise à un
déplacement vertical w(x) et à une rotation (x). A partir de la théorie des poutres
classique, la rotation (x) est obtenue comme la dérivée première de la déviation w(x);
c'est-à-dire (x) = dw(x)/dx. La pente (x) doit être continue, sinon la poutre
développerait des « plis » dans sa forme déviée. Par conséquent, si nous sommes sur
le point de construire une fonction approchée w(x) pour la flexion, alors la fonction
approchée et sa dérivée première doivent être continues. La dérivée seconde, qui
représente la courbure de la poutre, doit exister. Une fonction appropriée qui satisfait
à ces exigences serait:
w(x) = 1 + 2 x+ 3 x2 + 4 x3
Les quatre paramètres 1, 2, 3 et 4 peuvent être identifiés en utilisant les deux

valeurs nodales d'extrémité pour la flèche, w1, w2, et les deux valeurs d'extrémité pour
la pente, 1 et 2. Le problème du poutre est classé comme un problème C1. La
solution approchée et sa dérivée première doivent être continues, la dérivée seconde
doit exister.
En général, le principe de compatibilité peut être formulé comme suit :
• Pour un problème de classe C0 (continuité C0), la solution approchée doit être

continue à travers la frontière des éléments mais pas nécessairement ses dérivées.
• Pour un problème de classe C1 (continuité C1), la solution approchée et ses

dérivées du premier ordre doivent être continues à travers la frontière des éléments
mais pas nécessairement ses dérivées du second ordre.
• Pour un problème de classe Cn (continuité Cn), la solution approchée et ses

dérivées d'ordre (n−1) doivent être continues à travers la frontière des éléments mais
pas nécessairement ses dérivées d'ordre n.
2- Principe de la complétude
Encore une fois, considérons le problème des barres de la figure 7.4a. Si la force
appliquée P est différente de zéro, alors le déplacement u(x) a une valeur finie
différente de zéro en tout point x appartenant à la barre sauf en x = 0, où un
déplacement égal à zéro est imposé (condition aux limites) . Si nous choisissons de
discrétiser la barre avec un élément linéaire à deux nœuds, alors la fonction approchée
adoptée donnée dans l'équation (7.20) fera un choix approprié car si la taille des
éléments diminue jusqu'à zéro, c'est-à-dire limx 0 u( x) = 1, qui est une constante
représentant la valeur réelle du déplacement en ce point. Cependant, si la fonction
d'essai ne contenait pas de terme constant, limx 0 u(x) sera égal à zéro, ce qui ne
66
représente pas le cas réel. Dee plus, le terme constant est nécessaire pourr que la
fonction approchée puisse repré résenter un mouvement de corps rigide. Dans ce cas,
tous les points doivent avoirr le même déplacement u(x) = . De plus, us, on a
du(x)/dx= 2, ce qui représente le cas réel de la barre à déformation constantte. Cela
conduit à la définition du princcipe de complétude, qui peut être énoncé comm mme suit.
Lorsque la taille de l'élément see réduit à zéro, la fonction approchée doit être capable
de représenter :
• Pour un problème de claasse C0 (continuité C0), une valeur constantte de la

fonction exacte ainsi que des vaaleurs constantes de ses dérivées du premier ord
dre.
• Pour un problème de claasse C1 (continuité C1), une valeur constantte de la

fonction exacte ainsi que des valeurs constantes de ses dérivées du premier
er et du
second ordre.
• Pour un problème de claasse Cn (continuité Cn), une valeur constantte de la

fonction exacte ainsi que des vaaleurs constantes de ses dérivées jusqu'au nième ordre.
o
Ces conditions, telles qu qu'énoncées par les principes de compatib bilité et

d'exhaustivité, sont suffisantess pour garantir que la solution par élémen nts finis
converge vers la solution exactee. Heureusement, de nos jours, nous n'avons pass besoin
d'observer ces principes à chaaque fois que nous résolvons un problème avec la
méthode des éléments finis. Tous les éléments communs qui sont utilisés dans la
pratique ont été développés et vérifiés selon ces principes, et plus encore. . . Leurs
formulations géométriques et analytiques sont fournies dans les biblio othèques
d'éléments de la plupart des lo ogiciels d'analyse par éléments finis. Cependan nt, il ne
suffit jamais de rappeler que less so
solutions obtenues avec la méthode des élémen nts finis
ne sont que des approximations
ons dde la solution exacte. Par conséquent, il est intére
éressant
de comprendre ces principes afiafin d'évaluer la précision ou de faire un diagnos
gnostic d'un
modèle d'éléments finis.
Les fonctions de forme de L

Lagrange:
La plus simple est l'emploi d

des polynômes de Lagrange. Le polynôme de Laagrange
d’ordre i passe exactement par 1 au point xi et par 0 sur tous les autres pointts xj. On
peut donc l’utiliser comme foncction de forme :
Pour représenter un champ du premier degré le long d’un bord, ce qui suppose
suppos
donc deux inconnues, il faut deeux connecteurs indépendants donc deux nœudsuds : un à
chaque extrémité de l’arête. Pouour une barre du premier degré de longueur L et de
caractéristiques constantes,
s, on p
peut écrire directement : x1 =0 et x2 = L.
67
On retrouve les deux fonctions
ons d’interpolation précédemment calculées.
Par application de la formulee donnant l’expression des polynômes de Lagraange, on

détermine les trois fonctions d’’interpolation associées aux trois degrés de lib
berté de
cet élément de barre du second
ond d
degré :
binaison linéaire des trois monômes 1, x et x2, mais

Toute parabole est une comb m est
également une combinaison linéa
éaire unique des trois fonctions ci-dessus.
Les fonctions de forme d'H

'Hermite:
D'autres solutions peuvent exxister pour les fonctions de forme. On cite ici un seul
exemple les éléments finis d'Her
ermite qui ont la particularité d'avoir deux foncti
tions de
base associées à chaque nœud. Dans cette version, la valeur de la solution estt ajustée
avec la première fonction alorrs que la deuxième permet d'ajuster la valeu ur de la
dérivée. Ce type de fonctions de base peut avoir un intérêt pour la résolu ution de
certaines équations aux dérivéess partielles (par exemple l'équation des poutress ou des
plaques), même si elle nécessitee d'avoir deux fois plus de fonctions pour un maillage
m
donné.
Si on prend un exemple dee la poutre classique, la rotation est la dérivée

ée de la
flèche:
( )
= ( )
Si l'élément pris en compte eest l'élément poutre de deux nœuds avec deux
ux degrés
de liberté par nœuds, w et , on p
peut écrire l'approximation de w comme suit:
w(x) = 1 + 2 x+ 3 x2 + 4 x3
la rotation est la dérivée de w(x)), donc elle s'écrit comme suit:

( )= 2 + 2 3 x + 3 4 x2
Les conditions aux limites de la pou
poutre sont:
w(0)=w1, w(L)=w2
(0) = , ( )=
68
Après résolution de ces équations, les valeurs des coefficients a0, a1, a2 et a3
seront trouvé:
Puis, les fonctions de formes seront trouvées comme suit:
3x 2 2 x3 2 x 2 x3 3x 2 2 x 3 x 2 x3
N1 = 1 − + N 2 = x − + N 3 = − N 4 = − +
L2 L3 L L2 L2 L3 L L2
Avec:
w = N1w1 + N 2θ1 + N 3 w2 + N 4θ 2
69
Références
1- Abdelghani SEGHIR, Cours Méthode des Éléments Finis, (Polycopié

présenté à l'université Abderrahmane Mira – Bejaia).
2- Dimitrios G. Pavlou, Essentials of the Finite Element Method (2005, Elsevier

Inc).
3- Gouri Dhatt, Gilbert Touzot, Emmanuel Lefrançois- Finite element method.
(2012, ISTE Ltd and John Wiley & Sons, Inc).
4- J.N. Reddy - An Introduction to the Finite Element Method, 3rd Edition

(2005, McGraw-Hill Education (ISE Editions).
5- Jean-Charles Craveur, Modélisation des éléments finis,(2008, Dunod).
6- LOUHIBI M. Z, la méthode des éléments finis, (université Djillali Liabes -

Sidi Bel Abbes).
7- Michel Cazenave, Méthode des éléments finis Approche pratique en
mécanique des structures, (2010, Dunod).
8- Ouinas Djamel, Application de la méthode des éléments finis, Cours et

exercices corrigés, (2em édition, 2014, OPU).
9- Vincent Manet, Méthode des éléments finis, (2013- Vincent Manet).
70

Untitled 2

Transféré par

Droits d'auteur :

Formats disponibles

Untitled 2

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Untitled 2

Transféré par

Droits d'auteur :

Formats disponibles

Chapitre 6 : Approximation, Fonctions

Les fonctions de forme ou fon

Elles assurent le passage du p

Règles à respecter pour la cconstruction des fonction de forme:

1- Principe de compatibilité (de continuité):

La construction de la solutition approchée sur un élément fini doit essentieellement

Dans un contexte unidimens nsionnel, la déformation normale est en fait donnée

Considérons maintenant le problème du poutre. Sous la charge appliquée

Les quatre paramètres 1, 2, 3 et 4 peuvent être identifiés en utilisant les deux

En général, le principe de compatibilité peut être formulé comme suit :

• Pour un problème de classe C0 (continuité C0), la solution approchée doit être

• Pour un problème de classe C1 (continuité C1), la solution approchée et ses

• Pour un problème de classe Cn (continuité Cn), la solution approchée et ses

• Pour un problème de claasse C0 (continuité C0), une valeur constantte de la

• Pour un problème de claasse C1 (continuité C1), une valeur constantte de la

• Pour un problème de claasse Cn (continuité Cn), une valeur constantte de la

Ces conditions, telles qu qu'énoncées par les principes de compatib bilité et

Les fonctions de forme de L

La plus simple est l'emploi d

Par application de la formulee donnant l’expression des polynômes de Lagraange, on

binaison linéaire des trois monômes 1, x et x2, mais

Les fonctions de forme d'H

Si on prend un exemple dee la poutre classique, la rotation est la dérivée

la rotation est la dérivée de w(x)), donc elle s'écrit comme suit:

Puis, les fonctions de formes seront trouvées comme suit:

1- Abdelghani SEGHIR, Cours Méthode des Éléments Finis, (Polycopié

2- Dimitrios G. Pavlou, Essentials of the Finite Element Method (2005, Elsevier

4- J.N. Reddy - An Introduction to the Finite Element Method, 3rd Edition

6- LOUHIBI M. Z, la méthode des éléments finis, (université Djillali Liabes -

8- Ouinas Djamel, Application de la méthode des éléments finis, Cours et

9- Vincent Manet, Méthode des éléments finis, (2013- Vincent Manet).

Vous aimerez peut-être aussi