Calcul Dans EXCEL

Universit
Paul
Sabatier
Master 1 Conception et Calcul de Structures.
Toulouse III
Mthodes de Calcul dans Excel.

Michel SUDRE
Mars 2007
page 1
Mthode des Diffrences Finies.
1/5
Soit la recherche dune fonction f sur un domaine D. Cette fonction est solution dune quation
diffrentielle et vrifie des conditions imposes sur le bord du domaine D.
Exemples:
1 f(x): fonction dune seule variable qui reprsente la flche dune poutre reposant sur 2 appuis
et soumise son poids propre caractris par le poids linique p.
y
A
D est le domaine [A,B]
EIz
AB=L
f(x) vrifie: EIz.d4f = -p
dx4
B
f(0)=0
avec les conditions aux frontires:
d2f(0)
=0
dx2
f(L)=0
d2f(L)
=0
dx2
2 f(x,y): fonction de 2 variables qui reprsente la flche dune plaque carre reposant sur son
bord extrieur et soumise son poids propre caractris par le poids surfacique p.
C z
D est le domaine carr ABCD
E,
D
B
x
f(x,y) vrifie: (f) = -12p3(1-2)

E.e
y
A
reprsente le Laplacien
avec les conditions aux frontires:

f=0
sur [DA] et [BC]
f
2 =0
x
2
f=0
et
sur [AB] et [CD]

2f
2 =0
y
page 2
2/5
Avec la technique des diffrences finies, plutt que de chercher f en tous points de D, on se
limite chercher f en un nombre fini de points qui sont les noeuds dun maillage.
Exemples:
1 poutre:
1
...
I-1
I+1
...
y
2 plaque carre:
h
par symtrie,
on ne considre
quun 1/4 de la plaque.
h est le pas du maillage.
Les oprateurs de drivation vont tre remplacs par des combinaisons des valeurs prises par la
fonction f aux noeuds du maillage.
Soit le cas dune fonction dune seule variable: f(x). On cre une partition de laxe x avec un
pas h. Le noeud dabscisse x0 est repr par lindice I. Les noeuds voisins dabscisses x0 - h et
x0 + h sont reprs par I-1 et I+1.
On utilise le dveloppement de Taylor:
2
2
3
3
f(x0- h) = f(x0) - h.( df )x0 + h .( d 2f )x0 - h .( d 3f )x0 + ......
dx
2 dx
6 dx
[a]
en limitant au 2terme, on obtient:

f(x0- h)
f(x0) - h.( df )x0
dx
do lexpression de la drive premire en arrire:
( df )I
dx
f(x0) - f(x0 - h) = f(I) - f(I-1)

h
h
En partant du dveloppement:
2
2
3
3
f(x0+ h) = f(x0) + h.( df )x0 + h .( d 2f )x0 + h .( d 3f )x0 + ......
dx
2 dx
6 dx
[b]
on obtient de manire identique lexpression de la drive premire en avant:

( df )I
dx
f(x0 + h) - f(x0) = f(I+1) - f(I)

h
h
page 3
3/5
2
n
Il est possible dexprimer les oprateurs drive seconde: ( d 2f )I puis les suivants ( d nf )I .
dx
dx
A titre dexercice, montrer que:

2
en arrire:
( d 2f )I
dx
en avant:
( d 2f )I
dx
f(I) - 2.f(I-1) + f(I-2)

h2
f(I+2) - 2.f(I+1) + f(I)
h2
Tous ces rsultats jusqu n=4 peuvent tre regroups dans les 2 tableaux ci dessous:
arrire
2
3
4
h( df ) h2( d 2f ) h3( d 3f ) h4( d 4f )
dx
dx
dx
dx
avant
2
3
4
h( df ) h2( d 2f ) h3( d 3f ) h4( d 4f )
dx
dx
dx
dx
f(I)
f(I)
-1
-1
f(I-1)
-1
-2
-3
-4
f(I+1)
-2
-4
f(I+2)
-3
-1
-4
f(I+3)
-4
f(I+4)
f(I-2)
f(I-3)
f(I-4)
Diffrence au centre.
Considrons les expressions de ( df )I en arrire et en avant.
dx
Dans les 2 cas, le premier terme nglig dans le dveloppement de Taylor est
h2 .( d2f ) .
2 dx2 I
2
Donc le premier terme nglig dans lexpression de ( df )I est h .( d 2f )I .
dx
2 dx
Reprenons les expressions [a] et [b] de la page 2/5 et construisons [b] - [a]:
3
3
f(x0+ h) - f(x0- h) = 2.h.( df )x0 + h .( d 3f )x0 + ......
dx
3 dx
df
On obtient une nouvelle forme approche de ( )I :
dx
f(x
+
h)
f(x0- h) = f(I+1) - f(I-1)
df
0
( )I
2h
dx
2h
2 d3f
h
Dans cette nouvelle expression, le premier terme nglig est
.(
) .
6 dx3 I
[b] - [a]
Cette nouvelle formule dite au centre est donc plus prcise.

A titre dexercice, montrer que:
au centre:
( d 2f )I
dx
f(I+1) - 2.f(I) + f(I-1)

h2
page 4
4/5
Les formules au centre jusqu n=4 peuvent tre regroupes dans le tableau ci dessous:
centre
2
3
4
h( df ) h2( d 2f ) h3( d 3f ) h4( d 4f )
dx
dx
dx
dx
f(I2 )
f(I-1)
-1/2
-4
1
-2
f(I)
f(I+1)
-1/2
1/2
f(I+2)
-1
-4
1/2
En diffrence au centre, les expressions sont plus prcises. Il faudra donc les utiliser prioritairement. Les oprateurs en arrire et en avant pourront nanmoins tre exploits au voisinage
des frontires.
Dans le cas de fonctions de 2 variables x et y avec un pas h gal en x et en y, on pourra utiliser

les schmas suivants:
y
1
1
x
-4
1
h2
-1/4
1/
4
-1/4
/4
2
2
f= 2f + 2f
1
h2
2f
xy
1/
2
-1/2
-2
-1/2
/2
/2
1
h3
-1/2
-1/2
-1
-1/2
/2
/2
1
h3
1
-8
-8
20
-8
-8
1
h4
1
1
(f)
x
3f
xy2
(f)
page 5
5/5
La rsolution est traite dans le tableur Excel. Chaque noeud du maillage correspond une
cellule du tableur.
Les relations Diffrences Finies sont entres dans les cellules comme des formules. Le calcul
doit tre de type itratif cause des rfrences circulaires.
Il faut donc cocher la case Itration dans le menu Outils>Options>Calcul.
Exemple 1 : avec N=7.

poutre:
1
h=L/6
p=100 N/m
EIz=1.e6 Nm2
L=1m
page 6
Problmes avec conditions initiales.
1/6
Soit la recherche dune fonction x(t) vrifiant:

dx(t)
= f(x,t)
dt
avec la condition initiale: x(0) = x0
Exemple:
Un solide de masse m se dplace vitesse v(t). Il est soumis une force propulsive FP et une
force rsistante FR.
v(t)
FP(t)
la force propulsive FP est une fonction connue du temps
FR(v)
la force rsistante FR est une fonction connue de la vitesse
On dsigne par v0 = v(0) la vitesse initiale.

Lapplication du PFD conduit
m. dv(t) = FP(t)
dt
FR(v)
do:
dv(t)
= f(v,t)
dt
avec la condition initiale: v(0) = v0
Plutt que de chercher la fonction tout instant, on ralise un maillage du temps avec un pas
de temps t et on se limite une recherche aux noeuds du temps:
x0
0
...
I-1
I+1
...
recherche de x(t)
recherche de xI
I=1,2....
page 7
2/6
Mthode dEULER:
Si on connait x linstant t alors on connait f(x,t) et donc la drive dx(t) .
dt
Remplaons cette drive par une formule de diffrences finies.
Deux schmas sont envisags:
- schma explicite:
on crit au noeud I du temps:
f(xI,tI) = fI = ( dx )I
dt
x(tI+ t) - x(tI) = xI+1 - xI

t
t
do lexpression directe de la valeur inconnue:

xI+1 =
xI + fI.t
Cette mthode, trs peu prcise, nest pas utilisable.
- schma implicite:
on crit au noeud I+1 du temps:
f(xI+1,tI+1) = fI+1 = ( dx )I+1
dt
x(tI+ t) - x(tI) = xI+1 - xI

t
t
do lquation suivante:
xI+1
- f(xI+1,tI+1).t
= xI
qui ne permet pas le calcul immdiat de la valeur inconnue xI+1 .
page 8
3/6
Mthodes de RUNGE-KUTTA:
- formule dordre 2:
On crit comme dans la formule explicite prcdente:
x I+1 =
xI + f(xI,tI).t
(1) point P
puis en injectant le rsultat prcdent dans la formule implicite:

x I+1 =
xI + f(xI+1,tI+1).t
(2) point P
et on adopte finalement comme solution pour x I+1:

x I+1 = 1 ( x I+1 + x I+1)
2
Interprtation graphique:
x
P
xI+1
k2
P
k1
xI
I+1
k1= f(xI,tI).t
k2= f(xI+1,tI+1).t
xI+1 = 1 (k1+k2)
temps
Dans le cas dune quation du second ordre: par exemple loscillateur 1 DDL
x(t)
c
m
k
2
m. d 2x + c. dx + k.x = F(t)
dt
dt
F(t)
[a]
dx (0) = v
0
dt
x (0) = x0
On peut se ramener des quations diffrentielles dordre 1 en remplaant lquation [a] par un
systme:
dx
v
dt =
x (0) = x0
dv 1 . ( F(t) - c.v - k.x )

=
dt m
v (0) = v0
page 9
4/6
On calcule alors dans lordre:
x I+1 =
xI + V(xI,tI).t
(1) point P
puis:
V I+1
VI
(FI - c.vI - k.xI) .t
(2) point Q
puis:
x I+1 =
xI + V(xI+1,tI+1).t
(3) point P
puis:
V I+1
VI
1
(FI+1 - c.vI+1 - k.xI+1).t
+ m
(4) point Q
et on adopte finalement pour x I+1 et V I+1:

x I+1 = 1 ( x I+1 + x I+1)
2
V I+1 =
1
( V I+1 + V I+1)
2
x
P
k2
xI+1
k1= VI.t
P
xI
I+1
k2= V I+1 .t
k1
xI+1 = 1 (k1+k2)
temps
l2
VI+1
Q
VI
I+1
l1
temps
1
l1= m. (FI - c.vI - k.xI).t
l2=
1.
(FI+1 - c.vI+1 - k.xI+1).t

vI+1 = 1 (l1+l2)
2
page 10
5/6
- formule dordre 4:
Il sagit de la mthode la plus prcise et la plus utilise pour la recherche de x(t) vrifiant:
dx(t)
= f(x,t)
dt
avec la condition initiale: x(0) = x0
On calcule alors dans lordre:
p1 =
f( xI, tI)
p2 =
f( x I +
(1) pente au dpart
puis:
p1 t
, tI+ 1 )
2
(2) pente au point P
puis:
p2 t
p3 =
f ( x I+
p4 =
f( xI+p3t , tI+1)
, tI+ 1 )
2
(3) pente au point P
puis:
(4) pente la fin
et on adopte finalement pour x (I+1):

x (I+1) = x (I) + 1 (p1 + 2p2 + 2p3 + p4) .t
6
x
pente p4
1
pente p= (p1 + 2p2 + 2p3 + p4)
6
xI+1
pente p3
P
P
xI
I+ 12
I
t/2
pente p1
I+1
temps
pente p2
t/2
page 11
6/6
Problme de loscillateur 1 DDL dans Excel

x(t)
c
m
k
2
m. d 2x + c. dx + k.x = F(t)
dt
dt
F(t)
[a]
dx (0)
= v0
dt
x (0) = x0
- m, k, c et le pas de temps [delta] sont paramtrs.

- chaque ligne correspond un noeud du temps I=0,1,2,3....
- la nouvelle valeur de x est calcule partir de la valeur situe au dessus : L(-1)C
augmente de 1 (k1+k2) :
2
(LC(2)+LC(3))/2
page 12
Equations non-linaires.
1/5
Recherche des racines dune quation non-linaire R(u)=0.

Exemple: Assemblage de 2 barres:
A
E,S,L
E,S,L
1
Relation effort-dplacement:
u
F=
2ES
L
L
2
L -2hu+u2
-1 (h-u)
Mthode de Newton:
Ecrivons lquation sous la forme R(u)=0 et calculons la drive R(u) =
dR(u)
.
du
Une solution x0 doit tre propose pour amorcer le processus itratif.

Lors du passage de litration I litration I+1 , une approximation au 1 ordre est utilise:
R(uI+1) = R(uI) + R(uI) . [uI+1-uI]
Si on suppose qu lissue de cette itration, la solution est atteinte, alors:
R(uI+1) = 0
do laccroissement de uI vers uI+1:
On a donc le processus suivant : u0
u 1 = u 0-
R(u0)
R(u0)
uI+1 = uI -
u 2 = u 1-
R(uI)
R(uI)
R(u1)
R(u1)
.........
jusqu convergence de la solution.

R(u)
u
u0
u1
u2
page 13
2/5
Retour sur lexemple prcdent:

F
Equation dquilibre du noeud A:
F + 2N
h-u
=0
L+
[1]
reprsente la variation (algbrique) de longueur.

Gomtriquement, on exprime :
= L2-2hu+u2 - L
[2]
N = k
[3]
N est li par:
De ces 3 relations, on dduit:
F
phnomne du Snap-Through
F = 2k.
L
L2-2hu+u2
-1 (h-u)
2h
u
Rsolution dans Excel:
R(u) = 2k.
L
2
L -2hu+u2
-1 (h-u) - F
solution de dpart
R(u)
page 14
3/5
Rsolution dun systme de N quations non-linaires.

Exemple [N=2] : Systme de 2 barres en grands dplacements:
u:
v
composantes de dplacement du noeud A.
E,S,L
A
1
2
E,S,L
Mthode de Newton-Raphson:
Les 2 quations dquilibre sont crites sous la forme:
R1(u,v) = 0
[1]
R2(u,v) = 0
[2]
Do les expressions des drives partielles:

R1
= R1,u
u
R2
= R2,u
u
R1
= R1,v
v
R2
= R2,v
v
Passage de litration I litration I+1:
{ } { } [
R1
R2
I+1
R1
R2
R1,u R1,v
R2,u R2,v
] { } { }
x
(la solution est atteinte)
Do le systme:
R1,u R1,v
] { } { }
x
R2,u R2,v
I
(matrice Jacobienne)
= -
R1
R2
qui donne laccroissement
{ } {} {}
u
de
u
.
vers
v
v
I
I+1
page 15
4/5
Retour sur lexemple: Systme de 2 barres en grands dplacements

F
Equations dquilibre du noeud A:
N1. L-u
N2. u
=0
N1.v
N2. L+v
=F
N1
N2
do:
N1= F.
N2= F.
u
L-u+v
-u +
v2
FL u
u2
+
=
ES L-u+v
2L
2L
v+
v2
FL L-u
u2
=
+
ES L-u+v
2L
2L
puis:
L-u
L-u+v
Par combinaison, les 2 quations prcdentes peuvent tre transformes en:

2uL2 - 3u2L + 2u3- v2L + 2v2u +2uvL = 0
[1]
R1(u,v) = 0
[1]
3
2vL2 + 3v2L + 2v3 + u2L + 2u2v -2uvL - 2 FL = 0
ES
[2]
R2(u,v) = 0
[2]
Do:
R1
= 2L2 - 6uL + 6u2 + 2v2 +2vL
u
R2
= - 2vL + 4uv +2uL
u
R1
= - 2vL + 4uv +2uL
v
R2
= 2L2 + 6vL + 6v2 + 2u2 -2uL
v
Dans Excel, il faut prvoir la zone dfinie ci-dessous:
matrice Jacobienne
La mthode de Cramer est utilise pour rsoudre:
K11
K12
K21
K22
] { } { }
x
= -
R1
R2
page 16
5/5
Exemple prcdent dans Excel:
page 17

Calcul Dans EXCEL

Transféré par

Droits d'auteur :

Formats disponibles

Calcul Dans EXCEL

Transféré par

Informations du document

Description originale:

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

Calcul Dans EXCEL

Transféré par

Droits d'auteur :

Formats disponibles

Universit

Master 1 Conception et Calcul de Structures.

Mthodes de Calcul dans Excel.

Mthode des Diffrences Finies.

D est le domaine [A,B]

f(x,y) vrifie: (f) = -12p3(1-2)

avec les conditions aux frontires:

sur [AB] et [CD]

Mthode des Diffrences Finies.

h est le pas du maillage.

en limitant au 2terme, on obtient:

f(x0) - f(x0 - h) = f(I) - f(I-1)

on obtient de manire identique lexpression de la drive premire en avant:

f(x0 + h) - f(x0) = f(I+1) - f(I)

Mthode des Diffrences Finies.

A titre dexercice, montrer que:

f(I) - 2.f(I-1) + f(I-2)

Cette nouvelle formule dite au centre est donc plus prcise.

f(I+1) - 2.f(I) + f(I-1)

Mthode des Diffrences Finies.

Dans le cas de fonctions de 2 variables x et y avec un pas h gal en x et en y, on pourra utiliser

Mthode des Diffrences Finies.

Exemple 1 : avec N=7.

Problmes avec conditions initiales.

Soit la recherche dune fonction x(t) vrifiant:

la force propulsive FP est une fonction connue du temps

la force rsistante FR est une fonction connue de la vitesse

On dsigne par v0 = v(0) la vitesse initiale.

Problmes avec conditions initiales.

x(tI+ t) - x(tI) = xI+1 - xI

do lexpression directe de la valeur inconnue:

Cette mthode, trs peu prcise, nest pas utilisable.

x(tI+ t) - x(tI) = xI+1 - xI

qui ne permet pas le calcul immdiat de la valeur inconnue xI+1 .

Problmes avec conditions initiales.

puis en injectant le rsultat prcdent dans la formule implicite:

et on adopte finalement comme solution pour x I+1:

dv 1 . ( F(t) - c.v - k.x )

Problmes avec conditions initiales.

On calcule alors dans lordre:

(FI - c.vI - k.xI) .t

et on adopte finalement pour x I+1 et V I+1:

(FI+1 - c.vI+1 - k.xI+1).t

Problmes avec conditions initiales.

(1) pente au dpart

(2) pente au point P

(3) pente au point P

et on adopte finalement pour x (I+1):

Problmes avec conditions initiales.

Problme de loscillateur 1 DDL dans Excel

- m, k, c et le pas de temps [delta] sont paramtrs.

Recherche des racines dune quation non-linaire R(u)=0.

Une solution x0 doit tre propose pour amorcer le processus itratif.

do laccroissement de uI vers uI+1:

On a donc le processus suivant : u0

jusqu convergence de la solution.

Retour sur lexemple prcdent:

Equation dquilibre du noeud A: