ZQ Feng

IUP-Evry 1
Licence GMP - ME23
IUP, Université d’Evry – Val d’Essonne

MEF
Méthode des Eléments Finis

Zhi-Qiang FENG (CM, TD, TP) [Link]

Jean Michel CROS (TD, TP) [Link]
Gregory TURBELIN (TP) [Link]

ZQ Feng
IUP-Evry 2

Objectifs

u Apprendre la méthode des éléments finis (MEF)

u Maîtriser les concepts de base de la modélisation
numérique
u Être capable de résoudre des problèmes mécaniques
et physiques
u Apprendre à utiliser le code industriel ANSYS

Admissibilité : 3/5 examen + 2/5 TP ≥ 10/20

ZQ Feng
IUP-Evry Enseignement en calcul des 3

structures à l’IUP
ME23
ME23 :: Méthode
Méthodedes
des éléments
éléments finis
finis

CS31
CS31 :: Calcul
Calcul des
des structures
structures

DESS
DESS
Codes
Codes de
de calcul
calcul s
M SSSS cttuurrees
Mééca D
canniq DEESSSS DDEE sttrruuc
s
iquue
ennoon e ddeess
nlliinn iqiquu e
ééaair a
na m
m
iree DDyyn

ZQ Feng
IUP-Evry 4

Contenu de l’enseignement
u introduction
u méthodes matricielles pour résoudre les problèmes
discrets : système ressorts, réseaux électriques,
réseaux hydrauliques
u étapes logiques du calcul par éléments finis
u étapes pratiques du calcul par EF
u problèmes d'équilibre ou de valeurs aux limites
u approximation et discrétisation
u principe d’assemblage des matrices élémentaires
u introduction des conditions aux limites
ZQ Feng
IUP-Evry 5

Contenu de l’enseignement (suite)

u formulation intégrale
u méthode de résidus pondérés (Galerkin, méthode de
collocation …)
u éléments finis 1D, 2D et 3D
u systèmes de coordonnées
u éléments de référence et fonctions de forme
u transformation des dérivations et des intégrales
u intégration numérique (méthode de Gauss)
u équation aux dérivées partielles décrivant les
phénomènes physiques rencontré par l’ingénieur

ZQ Feng
IUP-Evry 6

Contenu de l’enseignement (suite)

u passage du continu au discret

u conditions aux limites portant sur la dérivée
u application en thermique (transfert de chaleur par
conduction et convection)
u application en électromagnétique
ZQ Feng
IUP-Evry 7

Méthode g énérale d'analyse num érique

Système physique Lois de la physique

Formulation des
Sciences de l’ingénieur
équations
Équations aux dérivées partielles

Transformation des Méthode des

équations Formulation intégrale
résidus pondérés

Système d’équations algébriques

Résolution Méthode de Gauss
numérique Méthode de Newton
Solution approchée (numérique)

ZQ Feng
IUP-Evry 8

Connaissances requises et à acquérir

u L'analyse numérique fait appel aux connaissances

dans les trois domaines suivants :
u Sciences de l'ingénieur pour établir les équations
u Méthodes numériques pour transformer ces équations
u Programmation et informatique pour exécuter
efficacement les calculs sur l'ordinateur.
ZQ Feng
IUP-Evry 9

Bibliographies

u Livres :
u The Finite Element Method
O.C. Zienkiewicz and R.L. Taylor
u Une présentation de la méthode des éléments finis
G. Dhatt et G. Touzot
u Modélisation des structures par éléments finis
J.L. Batoz et G. Dhatt

ZQ Feng
IUP-Evry 10

Bibliographies (suite)

u Revues :
u Int. J. for Numerical Methods in Engineering
u Computers & Structures
u Revue Européenne des Eléments Finis
ZQ Feng
IUP-Evry 11

Codes de calcul

u NASTRAN, ANSYS, ABAQUS, ADINA, SAP,

MARC, I-DEAS, DYNA3D, PAM/System,
RADIOSS, SAMCEF, …

u MEF/MOSAIC, SYSTUS, CESA, CASTEM,

ASTOR, ASTER, ZEBULON, MODULEF, …

u à l’IUP : ANSYS, I-DEAS, CASTEM, FER/Solid,

FER/Contact, ...

ZQ Feng
IUP-Evry 12

Domaines d'application de la MEF

u analyse linéaire (statique et dynamique)

u analyse non linéaire (grands déplacements, grandes
déformations, contact et frottement, flambage, ...)
u mise en forme des matériaux
u thermique (en régime permanent et transitoire, ...)
u mécanique des fluides
u électromagnétisme
u dynamique rapide (choc, impact, crash)
u optimisation des structures
ZQ Feng
IUP-Evry 13

Secteurs d'utilisation de la MEF

u génie mécanique
u génie civile
u transport
u aéronautique
u espace
u nucléaire
u énergétique
u militaire
u ....

ZQ Feng
IUP-Evry 14

Exemple: but ée de choc

DEMO
DEMO
• Grandes déformations
• Matériau: caoutchouc
• Mooney-Rivlin: hyperélasticité
incompressible
• Auto-contact multi-zones
ZQ Feng
IUP-Evry 15

Exemple: Projection à plasma

DEMO
DEMO

Ü Ü

Ü Ü

ZQ Feng
IUP-Evry Problèmes d'équilibre 16

(Système discret)
u Pour un système discret (système de ressorts,
réseaux électriques, réseaux hydrauliques,...), les
équations de comportement peuvent en général
s'écrire sous la forme matricielle suivante :
[K] {U} = {F}
[K] - matrice caractérisant le système
{U} - variables inconnues du problème
DDL – Degré De Liberté
DOF – Degree Of Freedom
{F} - sollicitations connues (second membre)
ZQ Feng
IUP-Evry 17

Signification (1D)

Problème K U F

Ressort raideur k déplacement force

Barre T/C rigidité EA/L déplacement force
Thermique conductivité température flux de
αA/L chaleur
Electricité résistivité 1/R potentiel intensité
Ecoulement perméabilité charge débit
αA/L hydraulique

ZQ Feng
IUP-Evry 18

Principes de la MEF

u La MEF est basée sur une idée

simple : subdiviser (discrétiser)
une forme complexe en un grand
nombre de sous-domaines
élémentaires de forme géométrique
simple (éléments finis)
interconnectés en des points
appelés nœuds.
ZQ Feng
IUP-Evry 19

Principes de la MEF (suite)

u Nous considérons le comportement mécanique de chaque

élément séparément, puis nous assemblons ces éléments
de telle façon que l’équilibre des forces et la compatibilité
des déplacements soient satisfaits en chaque nœud.
u La MEF utilise des approximations simples des variables
inconnues dans chaque élément pour transformer les
équations aux dérivées partielles en équations algébriques.
u Les nœuds et les éléments n’ont pas forcement de
signification physique particulière, mais sont basés sur des
considérations de précision de l’approximation.

ZQ Feng
IUP-Evry
Étapes logiques du calcul par 20

éléments finis
u Définir les nœuds et les éléments (Créer le maillage)
u Pour chaque élément, établir la matrice de rigidité élémentaire [ke]
reliant les degrés de libertés (déplacements) nodaux {ue} et les forces
{fe} appliquées aux nœuds : [ke] {ue} = {fe}
u Assembler les matrices et les vecteurs élémentaires en un système
global [K] {U} = {F} de manière à satisfaire les conditions
d’équilibre aux nœuds
u Modifier le système global en tenant compte des conditions aux
limites
u Résoudre le système [K] {U} = {F} et obtenir les déplacements {U}
aux nœuds
u Calculer les gradients (flux de chaleur, déformations et contraintes)
dans les éléments et les réactions aux nœuds sur lesquels les
conditions aux limites sont imposées.
ZQ Feng
IUP-Evry Assemblage 21

{ke}, {fe} ==> [K], {F}

u construction de la matrice étendue [Ke] et du vecteur
étendu {Fe } de chaque élément

0 0 0
 
 
[K ]
e
= 0

[k ]
e
0

 
0 0 0 

u addition des matrices et des vecteurs étendus

[K ] =  K 1  + K 2  + ... = ∑ [K ]
e
e

Exemple

ZQ Feng
IUP-Evry 22

Introduction des conditions aux limites

u méthode du terme diagonal dominant

u méthode du terme unité sur la diagonal
u méthode de suppression des équations

k1 k2 k3

CL : u1=u1
1 4 2 3

f1 f4 f2 f3

Exemple
ZQ Feng
IUP-Evry 23

Etapes pratiques du calcul par EF

Organisation générale d’un code de calcul

Pré-processeur

Solveur

Post-processeur

ZQ Feng
IUP-Evry 24

Pré-processeur

u Choisir le type d’éléments

u Entrer les propriétés géométriques
u Entrer les paramètres physiques
u Créer le modèle géométrique
u Créer le maillage : définir les nœuds et les éléments
u Appliquer les sollicitations
u Imposer les conditions aux limites
ZQ Feng
IUP-Evry 25

Solveur

u Choisir le type d’analyse (statique, dynamique,…)

u Construire la matrice et le vecteur élémentaire [k e], {fe}
u Assembler [ke] et {fe} dans [K] et {F}
u Prendre en compte les conditions aux limites
u Résoudre le système d’équations [K] {U} = {F}
u Calculer les variations additionnelles (gradients, réactions,
σ, ε…)

ZQ Feng
IUP-Evry 26

Post-processeur

u Présenter les résultats de façon intelligible et synthétique :

u sous forme numérique
u sous forme graphique
u Effectuer des fonctions complémentaires : combinaisons,
interprétations, interpolations, animation, …

Demo ANSYS
ZQ Feng
IUP-Evry Problèmes d'équilibre 27

(système continu)
u Pour un système continu, prenons comme exemple
le problème thermique :
∂2T ∂ 2T
K +K + Q = 0 dans V
∂x 2 ∂y 2
V ∂T
K = h ( T − Tf ) sur Sf
Q ∂n
Sf T=?
Comment ???
h, Tf
[K] {U} = {F}

ZQ Feng
IUP-Evry 28

Système continu

u D'une manière générale, le comportement d'un

système continu est décrit par les équations aux
dérivées partielles (EDP): S
L(u) + f v = 0 sur un domaine V V
C(u) = fs sur la frontière S de V
L, C- opérateurs différentiels caractérisant
le système,
u - fonctions inconnues,
fv, fs - fonctions connues appelées sollicitations.
ZQ Feng
IUP-Evry 29

Quelques notions

u un système discret est linéaire si les termes de [K] et

{F} sont des constantes indépendantes de {U}

u un système continu est linéaire si les expressions

L(u) et C(u) sont linéaires en u et ses dérivées. De
plus fv et f s sont indépendantes de u. Nous pouvons
alors écrire :
{L(u)} = [L]{u}
{C(u)} = [C]{u}
u un système d'équations différentielles est d'ordre m
s'il fait intervenir des dérivées de u jusqu'à l'ordre m

ZQ Feng
IUP-Evry 30

Quelques notions

u un opérateur différentiel L est dit homogène si :

L(u=0) = 0
u un système d'équations linéaires différentielles
[L]{u} + {fv} = 0
est dit homogène si : {fv} = 0
u les conditions aux limites
[C]{u} = {fs}
sont dites homogènes si : {fs} = 0
ZQ Feng
IUP-Evry 31

Quelques notions
u un système différentiel linéaire est dit auto-adjoint
ou symétrique si :

∫ <v> [L] {u} dV = ∫ <u> [L] {v} dV

Pour toutes les fcts u et v qui satisfont :
[C]{u} = [C]{v} = 0
u un système différentiel linéaire est dit positif si :

∫ <u> [L] {u} dV ≥ 0

Pour toutes les fcts u qui satisfont [C]{u} = 0

ZQ Feng
IUP-Evry 32

Approximation

u Un modèle mathématique d'un système physique fait

intervenir plusieurs variables ou fonctions, dites
exactes uex(x) : températures, déplacements,
potentiels, vitesse, etc. Celles-ci sont représentées
par des fonctions "approchées" u(x) telles que la
différence :
e(x) = u(x) - uex(x)
soit "petite" (de l'ordre de grandeur de la précision
désirée).
ZQ Feng
IUP-Evry Construction d’une fonction approchée 33

(Étape 1 )
u Choisir un ensemble fini de fonctions dépendant de
n paramètres ai : u(x, a 1, a2, ..., a n)
u(x) = P1 (x) a1 + P2 (x) a2 + ... + Pn(x) an = <P> {a}

u Les fonctions sont souvent choisies de manière à

être faciles à évaluer, à intégrer ou dériver
explicitement.
u Polynômes: u(x) = a1 + a2 x + ... + anx n-1
u fonctions trigonométriques:

u(x) = a 1sin(πx) + a2sin(2πx) + ... + ansin(nπx)

ZQ Feng
IUP-Evry Construction d’une fonction approchée 34

(Étape 2 )
u Déterminer les paramètres a1, a2, ..., an en faisant
coïncider uex(x) et u(x) en n points x1, x2, ..., xn,
c'est-à-dire en annulant e(x) en ces n points.

L'approximation peut fournir

u une solution approchée en tout point x d'une fonction
difficile à évaluer ou connue seulement en certains
points.
u une solution approchée d'une équation différentielle ou
aux dérivées partielles.
Exemple
ZQ Feng
IUP-Evry 35

Approximation nodale

u En général, les paramètres a1, a2, ..., an n'ont pas de sens

physique. Cependant nous pouvons choisir comme
paramètres ai les valeurs de la fonction u ex(x) en n points
appelés NOEUDS de coordonnées x1, x2 , ..., xn. Imposons
de plus que la fonction approchée u(x) coïncide avec la
fonction exacte uex(x) en ces nœuds.
u(x 1) = uex( x1 ) = u1
u(x 2) = uex( x2 ) = u2
......
u(x n) = uex( xn ) = un

ZQ Feng
IUP-Evry 36

Approximation nodale

u La fonction approchée (approximation globale)

u(x) = P1(x) a1 + P 2(x) a2 + ... + Pn (x) an = <P(x)> {a}
s'écrit alors (approximation nodale) :
u(x) = N 1(x) u1 + N2(x) u 2 + ... + Nn (x) un = <N(x)>{ue}

u Définitions :
{a} - paramètres généraux de l'approximation
{ue} - variables nodales de l'approximation
<P(x)> - fonctions de base de l'approximation
<N(x)> - fonctions d'interpolation (fonctions de forme)
ZQ Feng
IUP-Evry 37

Approximation nodale

u L'approximation nodale possède deux propriétés

fondamentales :
1) comme u(xi) = ui, les fonctions Ni vérifient :

Ni(xj) = δij (symbole de Kronecker)

2) l'erreur d'approximation s'annule en tous les
nœuds x i
e(xi) = 0
Exemple

ZQ Feng
IUP-Evry 38

Approximation nodale

u La méthode d'approximation nodale d'une fonction

d'une variable s'étend directement à l'approximation
de plusieurs fonctions de plusieurs variables.
u EXEMPLE : déplacements d'une structure 3D
u(x,y,z) = <N(x,y,z)> {un}
v(x,y,z) = <N(x,y,z)> {vn}
w(x,y,z) = <N(x,y,z)> {wn}
ZQ Feng
39
IUP-Evry
Discr étisation
Approximation nodale par sous-domaines
sous-domaines

u La construction d'une fonction approchée u(x) est difficile lorsque le

nombre de nœuds et donc de paramètres inconnus ui devient
important. Le problème se complique encore si le domaine V a une
forme complexe et si la fonction u(x) doit satisfaire des conditions
aux limites sur la frontière de V.
u La méthode d'approximation nodale par sous-domaines simplifie la
construction de u(x). Elle consiste à :
1) identifier un ensemble de sous-domaines Ve du domaine V.
2) définir une fonction approchée ue(x) différente sur chaque sous-
domaine par la méthode d'approximation nodale. Chaque fonction
ue(x) peut dépendre des variables nodales d'autres sous-domaines
comme c'est le cas dans l'approximation de type "Spline".

ZQ Feng
IUP-Evry Discr étisation 40

Approximation nodale par éléments finis

u La méthode d'approximation nodale par éléments finis est
une méthode particulière d'approximation nodale par sous-
domaines qui présente les particularités suivantes :
u L'approximation nodale sur chaque sous-domaine Ve ne fait
intervenir que les variables nodales attachées à des nœuds situés
sur Ve et sur sa frontière.
u Les fonctions approchées ue(x) sur chaque sous-domaine Ve sont
construites de manière à être continues sur Ve et elles satisfont des
conditions de continuité entre les différents sous-domaines. Les
sous-domaines Ve sont appelés des ELEMENTS connectés par
des NOEUDS.
ZQ Feng
IUP-Evry Discr étisation 41

Règles de partition du domaine en éléments

u Deux éléments distincts ne peuvent avoir en commun que
des nœuds situés sur leur frontières, si elle existe.

1D 2D 3D

u L'ensemble de tous les éléments doit constituer un domaine

aussi proche que possible du domaine donné.
u Le recouvrement de deux éléments et les trous entre
éléments sont inadmissibles.

ZQ Feng
IUP-Evry
Discr étisation 42

Erreur de discr étisation g éométrique

u Lorsque la frontière du domaine est constituée par des

courbes ou des surfaces plus complexes que celles qui
définissent les frontières des éléments, une erreur est
inévitable. Cette erreur est appelée "erreur de discrétisation
géométrique". Elle peut être réduite
u en diminuant la taille des éléments
u en utilisant des éléments à frontières plus complexes.
ZQ Feng
IUP-Evry Discr étisation 43

Formes d' éléments classiques

u 1D
u linéaire ou quadratique : L2, L3
u 2D
u éléments triangulaires: T3, T6
u éléments quadrilatéraux: Q4, Q8
u 3D
u éléments tétraédriques: TE4, TE10
u éléments hexaédriques: H8, H20
u éléments prismatiques: P6, P15

ZQ Feng
IUP-Evry 44

Formulation intégrale
Pour résoudre le système d'équations différentielles avec des
conditions aux limites :
S
L(u) + fv = 0 sur un domaine V
C(u) = fs sur la frontière S de V V

la meilleure solution est la solution analytique.

Mais, en pratique, la méthode analytique est inapplicable pour
des problèmes réels:
u domaine irrégulier
u multi-matériaux
u matériaux anisotropes
u équations non linéaires, .....
ZQ Feng
IUP-Evry 45

Méthodes numériques

u Trouver des solutions approchées aux points

discrets.
Les méthodes numériques sont classées en 3
catégories :
u Méthode des différences finies
u cette méthode marche plutôt pour des domaines rectangulaires
u il est très difficile d'écrire un code général pour cette méthode

u Méthode variationnelle
u Méthode des résidus pondérés
Base mathématique de la MEF

ZQ Feng
IUP-Evry 46

Méthode variationnelle
Pour résoudre l'équation différentielle suivante :
ìï d 2f
ïïD + Q = 0 sur 0 < x < H
ï dx 2
í
ïï
ïï f( x = 0) = f0 ; f( x =H ) = fH
î

la méthode variationnelle consiste à trouver une fonction test

φ = φ(x) telle que :
æ
H D æd f ö
2 ö
P = ò 0 ççç çç ÷÷÷ + Qf÷÷÷dx ® min
è
çè 2 dx ø ø÷

Cette méthode est inapplicable dans les cas où l'équation

différentielle contient des termes en dérivée première.
ZQ Feng
IUP-Evry 47

Méthode des résidus pondérés

u Supposons que u = h(x) est une solution approchée

de l'équation différentielle.
résidu : R(u) = L(u) + fv ≠ 0
u La méthode des résidus pondérés consiste à
rechercher des fonctions u qui annulent la
forme intégrale globale :
W = ∫ <ψ> {R(u)} dV = ∫ <ψ> {L(u) + fv} dV = 0
fonction de pondération (fonction test)
Le nombre de ψi(x) est égal au nombre de paramètres inconnus

ZQ Feng
IUP-Evry 48

Forme int égrale faible

u L'intégration par parties de la forme intégrale

globale fournit des formes intégrales dites faibles
qui présentent les avantages suivants :
u L'ordre maximum des dérivées de u qui apparaissent dans
la forme intégrale globale diminue. Les conditions de
dérivabilité de u sont donc moins fortes.
u Certaines des conditions aux limites peuvent être prises
en compte dans la forme intégrale faible.
ZQ Feng
IUP-Evry 49

Problème continu à deux dimensions

∂ 2 u ∂ 2u
u Équation de Poisson : + + fv = 0 sur V
∂x 2 ∂y 2

u Condition aux limites sur u (dite condition de Dirichlet) :

u = us sur Su
∂u
u Condition aux limites sur ∂n : (condition de flux)
∂u
∂n + α u = fs sur Sf
u si α = 0, cette condition est dite de Neuman
u si α ≠ 0, cette condition est dite de Cauchy

ZQ Feng
IUP-Evry 50

Forme intégrale globale

 ∂ 2u ∂ 2 u 
Wg =
∫ ψ (x , y ) 2 + 2 + f v dV = 0
 ∂x ∂y 
u u est dérivable deux fois et doit satisfaire toutes les
conditions aux limites sur Su et Sf. Les fonctions
ψ(x,y) ne sont soumises à aucune condition.
ZQ Feng
IUP-Evry 51

Forme int égrale faible

∫ ∫ ∫
 ∂ψ ∂u ∂ψ ∂u  ∂u ∂u
Wf = −  + − ψ f v dV + ψ dS + ψ dS
V  ∂x ∂x ∂y ∂y  Su ∂n Sf ∂n

u Les fonctions ψ(x,y) et u(x,y) doivent être une fois

dérivables. Les termes de contour peuvent être
déterminés avec les conditions aux limites sur Su et
Sf.
∂u
u sur Sf : ∂n = fs - α u
u sur Su, on peut imposer ψ = 0.

ZQ Feng
IUP-Evry 52

Forme int égrale faible

∫ ∫
 ∂ψ ∂u ∂ψ ∂u 
Wf = −  + − ψ f v dV + ψ (fs − α u )dS
V  ∂x ∂x ∂y ∂y  Sf

où u et ψ doivent satisfaire les conditions aux limites :

ψ = 0; u = u s sur Su
ZQ Feng
53

Choix de ψi(x)
IUP-Evry

Le choix du type de ψi(x) conduit à différentes

méthodes:

méthode de collocation
méthode de sous-domaines
méthode des moindres carrés

méthode de Galerkin

ZQ Feng
IUP-Evry 54

Méthode de collocation
u La fonction ψi(x) est la distribution de Dirac δ(xi) au point
xi, dit point de collocation. La forme intégrale s'écrit :

∫δ(xi) R(x, u) dV = R(xi,u) = 0

Autrement dit, R(xi, u) = 0 sur des points spécifiques. En
pratique, cette méthode est peu utilisée car elle est difficile à
mettre en oeuvre avec une approximation par éléments finis.
De plus elle conduit à un système d'équations non
symétrique. Par contre, elle a l'avantage d'éviter l'intégration
sur le volume.
ZQ Feng
IUP-Evry 55

Méthode de sous-domaines

ψi(x) = 1 sur certaines zones.

Autrement dit,

∫ R(u) dV = 0 sur des zones spécifiques.

Cette méthode est peu utilisée car le choix des sous-
domaines est difficile. De plus, elle nécessite des
intégrations sur le volume.

ZQ Feng
IUP-Evry 56

Méthode des moindres carr és

ψi(x) = R(x)
Cette méthode consiste à minimiser l'expression
Er = ∫ R(x)2 dV
par rapport aux paramètres a1 , a2 , ... an.
Les conditions de stationnarité sont : W = δEr = 0
Cette relation est équivalente aux n équations algébriques :

Wi(a) = ∫ L(Pi) (L(<P>{a})+fv ) dV = 0 i=1,2,...,n

Cette méthode est peu utilisée car elle ne permet pas l'intégration par
parties, et impose donc des conditions plus strictes sur l'approximation
de u que la méthode de Galerkin. Par contre, elle conduit à un système
symétrique et défini-positif quelle que soit l'équation différentielle.
ZQ Feng
IUP-Evry 57

Méthode de Galerkin
Les fonction ψ(x) sont constituées par l'ensemble des
variations δu des fonction u :
ψ = δu = <P> {δa} = <δa> {P} ∀ {δa}

W = ∫ δu (L(u) + fv) dV = 0

W = <δa> ∫{P}(L(<P>{a}) + fv ) dV = 0 ∀ {δa}

Comme W doit s'annuler pour tout {δa}, la relation précédente
est équivalente aux n équations algébriques :

Wi(a) = ∫ P i(L(<P>{a}) + fv) dV = 0 i = 1,..,n

Ce système est symétrique si L est auto-adjoint.
Exemple

ZQ Feng
IUP-Evry 58

Élément 1D
u Règles du maillage φ Transfert de chaleur
par conduction et
u la densité du maillage peut convection
être variable en fonction de
la variation des inconnues
u placer un nœud où D, Q
changent brutalement
u placer un nœud à l'endroit e
où on veut connaître la i j
solution
ìï d 2f
ïï D + Q = 0 sur 0 < x < H
ïï dx 2
ï
í
ïï
ïï f( x =0) = f0 ; D df = h( ff - fH )
ïïî dx ( x =H )