Dr: Belarbi M.
T Notions de fonctions d’interpolation
Chapitre III
NOTIONS DE FONCTIONS D’INTERPOLATION
SAVOIR CHOISIR
1-Introduction
L'analyse d'une structure en barres et poutres, par la méthode des déplacements, réalise la
compatibilité cinématique et l'équilibre, à l'aide de conditions formulées aux nœuds, entre
lesquels les barres et poutres constituent les éléments finis. Les propriétés de ces derniers
s'expriment par une relation matricielle liant les forces et déplacements nodaux, où intervient la
loi constitutive. Cette relation est exacte dans le cadre de la théorie des barres et poutres (Chap.
I).
Si l'on peut voir dans la méthode des éléments finis une extension de ces notions classiques, il
existe toutefois d'importantes différences par rapport aux poutres.
1.1 Interpolation
On ne connaît pas l'expression analytique des fonctions décrivant le champ des déplacements
dans les éléments finis pour un problème quelconque donné (Figure 3.1).
z,w(x,y)
1 déformé
4
2 w x
y
initial y
x
3
Figure 3.1 Elément rectangulaire de plaque fléchie à quatre nœuds : quelle déformée transversale
w(x,y) ?
On doit donc faire des hypothèses sur des fonctions: c'est le processus d'approximation ou
interpolation, dans lequel on choisit les fonctions d'interpolation. Dans les poutres, ces
fonctions sont connues de la mécanique des structures (Figure 3.2).
y,v(x)
déformé
v(x)
1 2
x EI
initial
Figure 3.2 Elément de poutre de Bernoulli à deux nœuds ; flèche v(x,y) connue: v cubique et M linéaire
(sans charge transversale entre les noeuds).
Université de Biskra, Département de Génie Civil 103
Dr: Belarbi M.T Notions de fonctions d’interpolation
1.2 Continuité cinématique
On doit s'efforcer d'assure la compatibilité cinématique des déplacements tout le long des
frontières séparant les éléments finis, et non pas aux seuls nœuds (Figure 3.3), sinon il en
résulterait des discontinuités des déplacements et des concentrations de contrainte inadmissible
sur ces frontières (figure 3.4).
initial
noeud
y
frontière
x
déformé
Figure 3.3 Elément fini d’élasticité plane, Figure 3.4 essai de traction discrétisé par
triangulaire, à six nœuds deux éléments à quatre nœuds,
avec discontinuité inadmissible
des déplacements.
1.3 Précision et convergence
L'interpolation peut elle se faire librement? Que vaut la solution approximative trouvée?
L'erreur est elle chiffrable? Dans quelles circonstances converge t on vers une solution exacte
(Figure 3.5, 3.6 et 3.7)? Ces importantes questions ont mis quelque temps à trouver des
réponses satisfaisantes, et certaines d'entre elles restent encore ouvertes aujourd'hui.
Figure 3.5 : Cylindre pincé avec diaphragmes rigides.
Université de Biskra, Département de Génie Civil 104
Dr: Belarbi M.T Notions de fonctions d’interpolation
Figure 3.6 : Cylindre pincé avec diaphragmes rigides. Discrétisation.
Figure 3.7 : Cylindre pincé avec diaphragmes rigides. Déplacement vertical en C.
Université de Biskra, Département de Génie Civil 105
Dr: Belarbi M.T Notions de fonctions d’interpolation
1.4 Construire un élément fini
La construction d'un élément fini implique de nombreux choix:
Forme géométrique;
Nœuds (nombre, emplacement);
Degrés de liberté (nature, nombre et type par noeud);
Fonction (s) d'inetrpolation.
Ces choix sont en outre soumis à certaines contraintes:
La construction d'un élément fini implique de nombreux choix:
Propres à la méthode des éléments finis (critère de convergence, Chap. II);
Liées à la modélisation, c'est-à-dire à la théorie choisie pour analyser la structure.
Ils sont essentiellement responsables de la qualité des résultats, dont ils soulignent -une fois de
plus- le caractère d'approximation, et constituent donc le problème clé de la méthode des
éléments fini.
Remarque:
C'est la fonction qui nous permette de décrire le déplacement et même les contraintes dans
le volume de l'élément en fonction des déplacements nodaux (ou contraintes nodaux).
les fonctions d'interpolation doivent être dérivables (calcul des déformations Chap. II) et,
dans une combinaison, linéairement indépendantes. Elles sont le plus souvent de forme
simple, presque toujours polynomiale; dans la suite, seule l'interpolation polynomiale est
retenue.
2- Champs interpolés
Soit :
Le champ de déplacement cinématiquement admissible choisi pour étudier un modèle
déplacement donné ; U(x,y,…) est le vecteur des composantes du déplacement, fonction d’un
certain nombre de paramètres on constantes , groupées dans le vecteur ;
x,y,… sont les variables indépendantes ; on écrit :
(3.1)
étant une matrice rectangulaire, dont chaque terme est appelé mode de
déplacement ; les paramètres définissant les intensités de ces modes ; le produit
représente une séquences de fonctions décrivant [U], appelées fonctions de déplacement.
Aucune de ces quantités n’a d’interprétation physique immédiate.
Les déplacements nodaux {di} de l’élément (inconnues) sont groupés dans le vecteur {d }; on
les choisit en nombre égal à celui les paramètres , qui vont être éliminés ; en effet, les {di}
peuvent être exprimés en fonction des à l’aide de (3.1), par introduction des coordonnées
des nœuds :
(3.2)
Université de Biskra, Département de Génie Civil 106
Dr: Belarbi M.T Notions de fonctions d’interpolation
La matrice [A] est carrée et ne contient que des nombres. on a :
(3.3)
Le champ de déplacement peut s’écrire :
(3.4)
Et par identification,
(3.5)
Les fonctions décrites par la matrice [N] se nomment fonctions d’interpolation (ou par fois
fonctions poids, fonction de forme) ; elles représentent le champ des déplacements en fonction
des valeurs nodales q de ces déplacements.
(3.6)
3- Élément, nœuds, déplacement:
Les six éléments décrits ci-dessous sont les plus simples et les plus connus dans l’analyse des
problèmes plans de l’élasticité ; les trois premiers ou les trois suivants peuvent être employés
simultanément dans une même discrétisation ; les trois éléments quadratiques, dont les
frontières se déforment en paraboles, sont évidemment meilleurs que les trois éléments
linéaires.
triangle linéaire :
Y
CHAMPS :
3
X
1 2
rectangle linéaire :
CHAMPS :
Y
4 3
X
Université de Biskra, Département de Génie Civil 107
1 2
Dr: Belarbi M.T Notions de fonctions d’interpolation
raidisseur linéaire :
X
1 2
triangle quadratique:
Y
5
6 4
X
1 2 3
rectangle quadratique:
6
7 5
8 4
X
1 2 3
raidisseur quadratique:
Université de Biskra, Département de Génie Civil 108
Dr: Belarbi M.T Notions de fonctions d’interpolation
X
1 2 3
4. Exemples d'applications
4.1- Exemple 1 : interpolation 1-D
U1 U2
X
L
1 2
Repère local (propre à l’élément)
Déplacements nodaux (U1 U2)
Pb : trouvez
2 nœuds (points) → interpolation linéaire
En écriture matriciel :
(*)
Université de Biskra, Département de Génie Civil 109
Dr: Belarbi M.T Notions de fonctions d’interpolation
Avec :
Ne : fonctions d’interpolation d’un élément
de : déplacements nodaux d’un élément
Où :
N1 : fonction d’interpolation propre du nœud 1
N2 : fonction d’interpolation propre du nœud 2
N1 N2
1 1 1
X X
1 2 1 2
Propriétés des fonctions d’interpolation :
Au nœud 1
a)-pour un nœud (i)
Aux autres nœuds (ici au nœud 2
b)- ce qui correspond d’ailleurs aux modes rigides de translation.
4.2- Exemple 2 : interpolation parabolique 1-D
U1 U2 U3
X
1-a 2 3
0 +a
Repère local (propre à l’élément)
Université de Biskra, Département de Génie Civil 110
Dr: Belarbi M.T Notions de fonctions d’interpolation
Déplacements nodaux (U1 U2 U3)
Pb : trouvez
3 points (nœuds) → parabolique
Avec :
Où :
Ou
Ou
Ou
Université de Biskra, Département de Génie Civil 111
Dr: Belarbi M.T Notions de fonctions d’interpolation
Figure. 3.8 Fonctions d'interpolation lagrangiennes; paraboliques.
Le déplacement à l’intérieur d’un élément (e) est donné par :
4.3- Exemple 3 : interpolation cubique 1-D
Même démarche que l'exemple 2, on obtient les fonctions d'interpolation suivantes:
Figure 3.9 Fonctions d'interpolation lagrangiennes; cubiques.
Remarque :
Ni l’épaisseur t de l’élément, ni les propriétés E (v) de son matériau (Hooke) n’interviennent à
ce niveau.
5. Propriétés des Ni - conformité
P1 : fonction de forme continues sur l’élément.
P2 : fonction de forme continues à l’interface entredeux éléments.
P3 : les fonctions de forme doivent être telles que l’on ne trouve pas de ‘’déformations’’ dans
l’élément lorsque les déplacements des nœuds résultent d’une translation ou une rotation
d’ensemble e l’élément (déplacement de corps rigide).
‘’Terme constant obligatoire’’.
P4 : les fonctions de forme doivent être telles que l’on puisse trouver éventuellement un champ
de déformation constant sur l’élément
‘’Terme linéaire obligatoire’’.
Université de Biskra, Département de Génie Civil 112
Dr: Belarbi M.T Notions de fonctions d’interpolation
*** Si les quatre propriétés sont vérifiées ‘’élément conforme’’.
Remarque :
Tout élément conforme est convergent.
Pour une interpolation polynomiale, le champ inconnu est aussi un polynôme le long d'une
frontière. Par conséquent, avec des polynômes complets, on assure la continuité C0 si, pour un
polynôme de degré p sur une frontière, on place p + 1 nœuds le long de cette frontière.
6. Eléments finis triangulaires
On construit quasi spontanément des éléments triangulaires à
nœuds (figure 3.10), correspondant à un polynôme complet de degré p en les variables
indépendantes x et y. avec continuité C0 du même degré sur chacune des trois frontières.
(a) (b) (c)
Figure 3.10 Eléments finis triangulaires avec continuité C0; (a) linéaire; (b) quadratique; (c) cubique.
Le polynôme d'interpolation s'écrit (composante u par exemple)
et ses termes apparaissent naturellement grâce au triangle de Pascal (figure 3.11).
Triangle Degré p Nbre de termes
1 constant 0 1
x y linéaire 1 3
x2 xy y2 quadratique 2 6
3 2 2
x xy xy y3 cubique 3 10
x4 x3y x2y2 xy3 y4 quartique 4 15
5 4 3 2 2 3 4
x xy xy xy xy y5 quintique 5 21
… … … … …
Figure 3.11 Triangle de Pascal.
Dans les éléments finis destinés aux calculs basés sur la convergence h (voir chap. II), on ne
dépasse pratiquement jamais le degré 3. La matrice [A] peut toujours être inversée sans
difficulté. On peut aussi dessiner aisément les fonctions d'interpolation (figure. 3.12) et par
suite en chercher les équations. On vérifie que leur somme vaut 1
Université de Biskra, Département de Génie Civil 113
Dr: Belarbi M.T Notions de fonctions d’interpolation
Figure 3.12 Fonctions d'interpolation du triangle quadratique.
7. Eléments finis rectangulaires
Plaçant l'origine des coordonnées au centre de l'élément rectangulaire, de dimension 2a x 2b, on
définit avantageusement les coordonnées naturelles
Qui valent sur les frontières du rectangle (figure. 3.13).
En partant de l'interpolation de Lagrange à une variable indépendante, on peut, pour les
domaines rectangulaires, construire deux familles de fonctions d'interpolation.
y,
4 3 (a,b)(1,1)
b
x,
1 a a 2
Figure 3.13 Elément rectangulaire et ses coordonnées cartésiennes (x, y) et naturelles (,).
Famille de Lagrange
Les nœuds sont placés aux intersections des lignes d'un quadrillage régulier (figure. 3.14) et
les fonctions d'interpolation peuvent s'obtenir par le produit des polynômes de Lagrange (3.10)
de chacune des coordonnées x et y, ou et . L'élément cubique est, déjà, rarement utilisé.
La figure 3.15 montre les fonctions d'interpolation typiques de l'élément quadratique. La
fonction du nœud central s'appelle une fonction bulle, pour des raisons évidentes; elle
s'associe à des degrés de liberté internes qui peuvent être condensés. En tout point, la somme
des fonctions vaut 1.
Université de Biskra, Département de Génie Civil 114
Dr: Belarbi M.T Notions de fonctions d’interpolation
Figure 3.15 Fonctions d'interpolation typiques de l'élément biquadratique.
4 3
(a)
1 2
9
(b)
1 5
(c) 13
1 2
Figure.3.14 Famille des rectangles de Lagrange et quelques fonctions d'interpolation: (a)
élément bilinéaire; (b) élément biquadratique; (c) élément bicubique.
L'interpolation paramétrique comprend les termes contenus dans un losange issu
du triangle de Pascal, comme l'indique la figure 3.16. Le polynôme est complet jusqu'au degré
p, puis incomplet jusqu'au degré 2p. Sur des lignes x (ou ) = cste, ou y (ou ) = cste,
interpolation est complète au degré p.
mode de
l’élément
1 biquadratique
2 2
3 2 2 2
4 3 2 2 3 2
4 3 2 2 3 4
4 2 3 3 2 4
4 3 3 4
4 4
Université de Biskra, Département de Génie Civil 115
Dr: Belarbi M.T Notions de fonctions d’interpolation
Figure 3.16 Termes de l'interpolation polynomiale des rectangles de Lagrange.
Famille de Serendip
Ces éléments sont construits sans nœuds internes (figure 3.17). Ils sont donc plus simples que
les précédents et, souvent, préférés. Les fonctions d'interpolation sont faciles à visualiser. En
tout point d'un élément, leur somme est encore unité. L'élément bilinéaire est identique à celui
de la famille de Lagrange.
(a) (b) (c)
(d)
Nœuds d'angle:
Nœuds médians:
( i = numéro du nœud; (i,i)= coordonnées du nœud i.)
Figure. 3.17 Famille des rectangles de Serendip: (a) élément bilinéaire ; (b) élément biquadratique ;
(c) élément bicubique; (d) fonctions d'interpolation de l'élément biquadratique.
L'interpolation paramétrique est encore extrait du triangle de Pascal (Figure 3.18). Tant que
p3, le polynôme, de degré p+1, est complet au degré p. quel que soit p, il demeure complet
au degré p sur les lignes x (ou ) = cste, ou y (ou ) = cste. en pratique, on ne dépasse guère
p=3.
mode de
l’élément
1
bicubique
2 2
3 2 2 2
4 3 3 2
4 4
Figure 3.18 Mode de l'interpolation polynomiale des rectangles de Serendip.
Université de Biskra, Département de Génie Civil 116
Dr: Belarbi M.T Notions de fonctions d’interpolation
Université de Biskra, Département de Génie Civil 117