A.

ENONCÉ

Figure 1: A single story frame

Let us considerate the structure in figure 1 subjected to the loading history shown on figure 2
below.

A.I. Part 1

We suppose that the structure has a linear elastic behavior.

Figure 2: Force variation

Parameters Expression
i Rank of group
h ( i )=4.5+0.1 (−1 ) ( i−1 ) [ m ]
i
Height, h
l ( i ) =6.2+0.1 (−1 ) ( i−1 ) [ m ]
i
Span, l
Load, p p ( i )=37+i [ kN /m ]
ζ ( i )=48+2.5 (−1 ) ( i−1 ) [ % ]
i
Damping factor
Young’s modulus E=[ 200+5 ( i−1 ) ] GPa
 Inertia I =[ 9875+ 10 (−1 )i ( i−1 ) ] ×10 4 mm4
Time step Δ t=0.01 s
Time of simulation T =1 s
First method Linear acceleration method
Second method Runge-Kutta method of order 4
Third method* Finite difference method if i≡ 1 [ 7 ]
Wilson method ( θ=1.38 ) if i≡ 2 [ 7 ]
Duhamel integral if i≡ 3 [ 7 ]
Average acceleration method if i≡ 4 [ 7 ]
Runge-Kutta method of order 3 if i≡ 5 [ 7 ]
Central difference method if i≡ 6 [ 7 ]
Houbolt method i≡ 0 [ 7 ]

a) Give the theoretical basis, the algorithm and the flowchart of the first method. Use
this method to compute displacement y ( t ) , velocity ẏ ( t ) and acceleration ÿ ( t ) . Draw
y ( t ) , ẏ ( t ) and ÿ ( t ) separately;
b) Do the same for the second method;
c) Do the same for the third method;
d) Plot displacement of the three methods on the same figure. Compare the three
methods.

A.II. Part 2

We suppose that the structure has an elastoplastic behavior illustrates on the figure 3 below.

Parameters Expression
i Rank of group
Restoring force r ( i )=25.4 i+0.1 (−1 ) î ( i−1 ) [ kN ]

Figure 3: Elastoplastic behavior

a) Give the algorithm and the flowchart of the first method. Use this method to compute
displacement y ( t ) , velocity ẏ ( t ) and acceleration ÿ ( t ) . Draw y ( t ) , ẏ ( t ) and ÿ ( t )
separately ;
b) Do the same for the second method;
c) Do the same for the third method;
d) Plot displacement of the three methods on the same figure. Compare the three
methods.

A.III. Part 3
a) Plot displacements of the first method with the two types of behavior on the same
figure. Comment the figure ;
b) Plot velocities of the second method with the two types of behavior on the same
figure. Comment the figure ;
c) Plot accelerations of the third method with the two types of behavior on the same
figure. Comment the figure ;
B. Données de l’exercice
 B.I. Paramètre du système (à copier dans le fichier Excel)

Différents paramètres du système

Paramètres Expressions Valeurs
i Numéro du groupe 5

Hauteur h i 4,1
h ( i )=4 , 5+0 , 1 (−1 ) ( i−1 ) [m]

Portée L i
L ( i ) =6 , 2+ 0 ,1 (−1 ) ( i−1 ) [m] 5,8

Poids de la structure p ( i )=37+i ¿ 42

Module de Young E=( 200+5 ( i−1 ) ) GPa 220

Inertie i
I =(9875+10 (−1 ) (i−1 ))×10 mm
4 4 98350000

Pas de variation temporel ∆t[s] 0,1

Temps de simulation T[s] 1

Facteur d'amortissement ζ =48+2,5(−1)i(i−1 ¿[% ] 38

 Caractéristiques intrinsèques du système

Caractéristiques Expressions Valeurs

Masse m m= p∗L/g [kg] 24857,14286

Rigidité K (6 E∗I ) 7534539,545

K= ¿
(h)3

Pulsation propre W 17,41015967

W=
K
m √
[rad /s ]

C Facteur d'amortissement 13,23172135

 B.II. Expression de la force extérieure

Le chargement étant variable dans le temps, son expression l’est également. On a la

fonction escalier définie par (le temps est exprimé en secondes et la force en Kilonewton) :

 t ∈ [ 0 ; 0.2 ] : F ( t )=55 i √ 1−25(t−0.20)2 (elliptique)

 t ∈ [ 0.2; 0.3 ] : F ( t )=55 i (constante)
 t ∈ [ 0.3 ; 0.5 ] : F ( t )=−275 i∗t+137.5 i (linéaire)
 t ∈ [ 0.5 ; 0.6 ] : F ( t )=0 (constante)
 t ∈ [ 0.6 ; 0.7 ] : F ( t )=4400i ( t−0.6)2 (parabolique)
 t ∈ [ 0.7 ; 0.95 ] : F ( t )=−85.60 i∗t 3 +73.3592 i (cubique)

B.III. Première partie : La structure a un comportement élastique linéaire

 Modélisation du système

Les deux montants du système ayant la même hauteur et bi encastrés de part et d’autre aux
12 EI
différents nœuds, ils ont la même rigidité : soit K 1=K 2 = 3 ; La traverse est supposée de
h
rigidité infinie.

Les deux montants étant en parallèle, la rigidité équivalente du système est la somme des
6 EI
deux rigidités : soit K=K 1+ K 2= 3 . On obtient le modèle suivant :
h

Figure 4: Modélisation du système

  Equation du système : D’après le second principe de Newton sur la mécanique, on a
l’équation du système donnée par : m ÿ ( t ) +c ẏ ( t ) +ky ( t )=F(t)

B.I. Première méthode : Méthode de l’accélération linéaire

1.a. Base théorique

La méthode de l’accélération linéaire est utilisée pour résoudre les équations

différentielles de mouvement. Généralement utilisé lorsque la sollicitation externe est
quelconque.

F (t )
ÿ ( t ) +2 ζ ω n ẏ ( t )+ ωn2 y (t )= (1)
m

Elle fait partie de la famille de méthodes qui se basent sur une approximation des intégrales
suivantes :

t + Δt
ẏ ( t+ Δt )= ẏ ( t ) + ∫ ÿ ( t+ Δt ) dt (2)
t

t + Δt
y ( t+ Δt )= y ( t ) + ∫ ẏ ( t+ Δt ) dt (3)
t

Et qui postulent une expression analytique relativement simple pour ẍ de façon à

pouvoir réaliser ces intégrales explicitement. En l’occurrence, la méthode de l’accélération
linéaire (par pas de temps Δt ) suppose que :

ÿ ( t + Δt )− ÿ ( t )
y⃛ (t )= (4)
Δt

En utilisant le développement en série de fourrier on a :

2 3
(Δt ) ( Δt )
y ( t+ Δt )= y ( t ) + ẏ ( t ) Δt+ ÿ ( t ) + ⃛y ( t ) + … (5)
2 6

( Δt )2
ẏ ( t+ Δt )= ẏ ( t ) + ÿ ( t ) Δt+ ⃛y ( t ) +… (6)
2

En introduisant (4) dans (5) et (6) on a :

 2 2
(Δt ) ( Δt)
y ( t+ Δt )= y ( t ) + ẏ ( t ) Δt+ ÿ ( t ) + ÿ ( t+ Δt ) (7)
3 6

ÿ ( t ) Δt ÿ ( t+ Δt ) Δ t
ẏ ( t+ Δt )= ẏ ( t ) + + (8)
2 2

En introduisant ces deux relations dans l’équation du mouvement en t+ Δt on a :

2 2
ÿ ( t ) Δt ÿ ( t + Δt ) Δ t ( Δt ) ( Δt ) F ( t + Δt )
ÿ ( t+ Δt )+ 2 ζ ωn ( ẏ ( t )+ + )+ω n2( y ( t ) + ẏ ( t ) Δt + ÿ ( t ) + ÿ ( t ) )=
2 2 3 6 m

F ( t+ Δt )−By ( t )−C ẏ ( t )−D ÿ (t )

⟹ ÿ ( t+ Δt )=
A

Avec :

( Δt )2 2 2 ( Δt )2
A=(1+ζ ωn Δt +ω n2 ), B=ω n , C=(2ζ ωn +ω n Δt ) D=(ζ ω n Δt+ ωn2 )
6 3

1.b. Organigramme de la méthode

Les calculs se font pour une durée de simulation T fixé et un pas Δt . L’intervalle de
T
temps [ 0 , T ] est donc subdivisé en N= intervalles [ t n ,t n+1 ] avec t 0=0 , t 1=Δt , t 2=2 Δt …
Δt

Les calculs se font de façon itérative de la façon suivante :

 Conditions initiales

ÿ ( t 0 )=a

ẏ ( t 0 )=b

y ( t 0 )=c

 Itérations (calcul de ÿ ( t n+1 ), ẏ ( t n ), y ( t n ) )

F ( t n+1 )
−By ( t n )−C ẏ ( t n ) −D ÿ ( t n )
m
ÿ ( t n+1 )=
A

ÿ ( t n ) Δt ÿ ( t n+1 ) Δ t
ẏ ( t n+1 )= ẏ ( t n ) + +
2 2
2 2
(Δt ) ( Δt )
y ( t n+1 )= y ( t n ) + ẏ ( t n ) Δt+ ÿ ( t n ) + ÿ ( t n+1 )
3 6
 On obtient ainsi les valeurs des déplacements, vitesses et accélération pour tous les
points t n on peut donc à partir de ceux-ci construire la courbe de réponses du système

1.c.Algorithme

Figure 5: Algorithme permettant de transformer le vecteur force en vecteur accélération pour

la mise en œuvre de la méthode de l'accélération linéaire
Figure 6: Méthode de l'accélération linéaire
 1.d. Valeurs prises par l’accélération, la vitesse et le déplacement sur le temps
de simulation

i t F(t) y''(t) y'(t) Y(t)

0 0 0 0 0 0
1 0,01 85868,72 3,19471562 0,01597357 5,32453E-05
2 9 8
2 0,02 119869,7 3,97923364 0,05184332 0,000385792
2 7 5
3 0,03 144865,2 4,22116381 0,09284531 0,001107219
4 1 2
4 0,04 165000 4,13432173 0,13462274 0,002245283
4
5 0,05 181895,4 3,82169585 0,17440282 0,003793016
2 7
6 0,06 196389,2 3,35084914 0,21026555 0,005720282
8 2
7 0,07 208981,9 2,77295333 0,24088456 0,007980848
1 4 5
8 0,08 220000 2,12946242 0,26539664 0,010517617
9 4
9 0,09 229670,2 1,45487605 0,28331833 0,013266813
8 5 6
10 0,1 238156,9 0,77801934 0,29448281 0,016161459
9 8 3
11 0,11 245582,8 0,12270969 0,29898645 0,019134267
5 9 8
12 0,12 252041,6 -0,49183104 0,29714085 0,022120024
6 2
13 0,13 257606,1 -1,05065849 0,28942840 0,025057528
7 4
14 0,14 262333,2 -1,54280779 0,27646107 0,027891076
8 3
15 0,15 266267,6 -1,96098299 0,25894211 0,030571577
 1 9
16 0,16 269443,8 -2,30122179 0,23763109 0,033057278
7 5
17 0,17 271888,6 -2,56253509 0,21331231 0,035314173
5
18 0,18 273621,5 -2,74652618 0,18676700 0,037316103
5 4
19 0,19 274656,0 -2,85699779 0,15874938 0,039044605
3 4
20 0,2 275000 -2,8995563 0,12996661 0,04048854
4
21 0,21 275000 -2,86842536 0,10112670 0,041643747
6
22 0,22 275000 -2,76069759 0,07298109 0,042513388
1
23 0,23 275000 -2,58876561 0,04623377 0,04310803
5
24 0,24 275000 -2,36530534 0,02146342 0,043444654
25 0,25 275000 -2,10288601 -0,00087754 0,043545396
26 0,26 275000 -1,81363089 -0,02046012 0,043436298
27 0,27 275000 -1,50893201 -0,03707294 0,043146093
28 0,28 275000 -1,19922033 -0,0506137 0,042705079
29 0,29 275000 -0,89379091 -0,06107875 0,042144072
30 0,3 275000 -0,60068139 -0,06855112 0,04149348
31 0,31 261250 -0,83816445 -0,07574534 0,040773976
32 0,32 247500 -1,02309266 -0,08505163 0,039971533
33 0,33 233750 -1,15671315 -0,09595066 0,039067635
34 0,34 220000 -1,24155866 -0,10794202 0,038048878
35 0,35 206250 -1,28121924 -0,12055591 0,036906719
36 0,36 192500 -1,28011225 -0,13336257 0,035637118
37 0,37 178750 -1,24325729 -0,14597941 0,034240101
38 0,38 165000 -1,17606178 -0,15807601 0,032719263
39 0,39 151250 -1,08412197 -0,16937693 0,031081233
40 0,4 137500 -0,97304293 -0,17966275 0,029335109
41 0,41 123750 -0,84828036 -0,18876937 0,027491908
42 0,42 110000 -0,71500594 -0,1965858 0,025564022
43 0,43 96250 -0,5779971 -0,20305081 0,023564697
44 0,44 82500 -0,44155138 -0,20814856 0,021507563
45 0,45 68750 -0,30942463 -0,21190344 0,019406202
46 0,46 55000 -0,18479208 -0,21437452 0,017273774
47 0,47 41250 -0,07023033 -0,21564963 0,015122698
48 0,48 27500 0,03228162 -0,21583938 0,012964399
7
49 0,49 13750 0,1213447 -0,21507124 0,010809104
50 0,5 0 0,19610801 -0,21348398 0,008665705
51 0,51 0 0,76778955 -0,20866449 0,006550198
1
52 0,52 0 1,24799443 -0,19858557 0,004509946
3
53 0,53 0 1,63445543 -0,18417332 0,002592931
7
54 0,54 0 1,92781055 -0,16636199 0,00083781
6
55 0,55 0 2,13122960 -0,14606679 -
7 0,000726029
56 0,56 0 2,25001557 -0,12416057 -
2 0,002078156
57 0,57 0 2,29119491 -0,10145451 -
7 0,003206574
58 0,58 0 2,26310973 -0,07868299 -
6 0,004107028
59 0,59 0 2,17502292 -0,05649233 -0,00478217
2
60 0,6 0 2,03674583 -0,03543348 -
8 0,005240647
61 0,61 2200 1,94014640 -0,01554902 -
4 0,005494755
62 0,62 8800 1,96468274 0,00397512 -
 4 3 0,005552829
63 0,63 19800 2,09484416 0,02427275 -
7 7 0,005412674
64 0,64 35200 2,31405713 0,04631726 -
8 4 0,005061551
65 0,65 55000 2,60527991 0,07091394 -
1 9 0,004477821
66 0,66 79200 2,95153753 0,09869803 -
6 0,003632647
67 0,67 107800 3,33638949 0,13013767 -
1 0,002491676
68 0,68 140800 3,74432518 0,16554124 -0,00101668
8 5
69 0,69 178200 4,16108491 0,20506829 0,000832894
8 5
70 0,7 220000 4,57390652 0,24874325 0,003098512
9 3
71 0,71 213621,4 3,01553401 0,28669045 0,005788667
7 8 5
72 0,72 207061,8 1,55658568 0,30955105 0,008782032
3 6 4
73 0,73 200317,4 0,22565546 0,31846226 0,01193319
3 8
74 0,74 193385,7 -0,95570924 0,31481199 0,015109406
1
75 0,75 186264,0 -1,97279024 0,30016949 0,018192789
6 3
76 0,76 178949,9 -2,81727548 0,27621916 0,02108177
6 5
77 0,77 171440,8 -3,48670812 0,24469924 0,02369194
4 7
78 0,78 163734,1 -3,98383892 0,20734651 0,025956312
1 2
79 0,79 155827,2 -4,31591195 0,16584775 0,02782505
 2 7
80 0,8 147717,6 -4,49391165 0,12179863 0,029264766
9
81 0,81 139402,6 -4,53179722 0,07667009 0,030257425
8 5
82 0,82 130879,9 -4,44574737 0,03178237 0,03079897
1 2
83 0,83 122146,7 -4,25343515 -0,01171354 0,030897712
84 0,84 113200,5 -3,97334973 -0,05284747 0,030572573
85 0,85 104038,7 -3,62417826 -0,09083511 0,02985125
4
86 0,86 94658,84 -3,22425778 -0,12507729 0,028768356
9
87 0,87 85058,26 -2,79110408 -0,15515409 0,027363589
7
88 0,88 75234,42 -2,34102119 -0,18081472 0,025679994
6
89 0,89 65184,75 -1,88879263 -0,20196379 0,023762333
9
90 0,9 54906,7 -1,44745306 -0,21864502 0,021655611
91 0,91 44397,68 -1,02813692 -0,23102297 0,019403777
3
92 0,92 33655,14 -0,63999881 -0,23936365 0,01704861
2
93 0,93 22676,51 -0,29019916 -0,24401464 0,014628803
1
94 0,94 11459,22 0,01605247 -0,24538537 0,012179251
3 9
95 0,95 0 0,27537823 -0,24392822 0,009730522
8
96 0,96 0 0,92158845 -0,23794338 0,007315779
3
97 0,97 0 1,46289102 -0,22602099 0,004991446
2
 98 0,98 0 1,89694956 -0,20922178 0,002811615
6
99 0,99 0 2,22471980 -0,18861344 0,000819708
5
100 1 0 2,45001993 -0,16523974 -
2 0,000951436
 1.e.Courbes

deplacement Y(t) vitesse y'(t)

0.06
0.6

0.05
0.4
0.04

0.03 0.2

0.02 0
0 0.2 0.4 0.6 0.8 1 1.2
0.01
-0.2
0
0 0.2 0.4 0.6 0.8 1 1.2
-0.4
-0.01

-0.02 -0.6

accelerati on y"(t)
0.05

0.04

0.03

0.02

0.01

0
0 0.2 0.4 0.6 0.8 1 1.2

-0.01
 B.II. Deuxième méthode : Méthode de Runge Kutta d’ordre 4
2. a. Base théorique

Les méthodes de Runge-Kutta sont des méthodes d'analyse

numérique d'approximation de solutions d'équations différentielles. Elles ont été
nommées ainsi en l'honneur des mathématiciens Carl Runge et Martin Wilhelm
Kutta lesquels élaborèrent la méthode en 1901.

Ces méthodes reposent sur le principe de l'itération, c'est-à-dire qu'une

première estimation de la solution est utilisée pour calculer une seconde estimation,
plus précise, et ainsi de suite.

Il en existe plusieurs mais celle que nous traiterons est la méthode d’ordre 4
qui est une amélioration de la méthode de base. Elle possède une étape supplémentaire
dans le raffinement du calcul.

Les méthodes de Runge Kutta s’écrivent de manière générale

K 1=h¿

K 2=h¿

………… …….

K n=h ¿

y k +1= y k +h ¿

Les coefficients sont déterminés à fin que l’ordre soit le plus élevé possible. On note A la
matrice de coefficients (α i , j), ᴦ le vecteur des coefficients γi et Ө le vecteur coefficients de θi

Une représentation en forme de tableau des équation donne

γ1 γ2 ……… γn
θ1 α 1 ,1 α 1 ,2 …………. α 1 ,n
θ2 α 2 ,1 α 2 ,2 …………. α 2 ,n
………….. ………… ……………. ……...
θn α n ,1 α n ,2 ……………. α n ,n
 Avec cette représentation, la méthode Runge-Kutta explicite à deux points qui
est d’ordre 4 est représentée par le tableau suivant

1/6 1/3 1/3 1/6

0 0 0 0 0
1/2 1/2 0 0 0
1/2 0 1/2 0 0
1 0 0 1 0

2. b. Organigramme
1.a. Calculs initiaux
i. Calculez k, c, m.
ii. Initialiser les variables
v 0=u̇0
1 ˙
v˙0=
m
[ F ( 0 ) −c v 0−k u 0 ]

iii. Sélectionnez un pas de temps approprié ∆t.

1.b. Pour chaque pas de temps

Calcul au début de l'intervalle de temps.

t=t i
x=u=ui
v 1=v i
1
v̇ 1=v i
m
[ F ( t ) −c v i−k ui ]

i.Calcul au premier point médian de l'intervalle de temps

Δt
t=t i+
2
vi Δ t
x=u Δt =ui +
i+
2
2
v˙1 Δ t
v 2=v Δt =v i +
i+
2
2
 1
v̇ 2=v̇
i+
Δt =
m
[ F ( t ) −c v 2−kx ]
2

ii. Calcul au deuxième point médian de l'intervalle de temps

Δt
t=t i+
2
Δt
x=u Δt =ui +v 2
i+
2
2
Δt
v 3=v Δt =v i + v˙2
i+
2
2
1
v˙3=
m
[ F ( t )−c v 3−kx ]
iii.Calculer les déplacements à la fin de l'intervalle de temps
t=t i+ Δ t
x=ui + Δt =ui + v 3 Δ t
v 4=v i+ Δ t=v i+ v˙3 Δ t
1
v˙4=
m
[ F ( t )−c v 4−kx ]
iv.Calculer le déplacement et la vitesse à la fin de l'intervalle de temps
Δt
ui +1=ui +
6 1
[ v +2 v 2 +2 v 3 + v 4 ]
Δt
v i+1=v i+
6 1
[ v˙ + 2 v˙2+ 2 v˙3+ v˙4 ]
 2. c. Valeurs prises par le déplacement, la vitesse et l’accélération sur le
temps de simulation

F(t) i t y(t) y'(t) y''(t)

0 0 0,01 0 0 0
61,106336 1
6
85,868722 2 0,02 7,92327E-05 0,02103567 3,1521341
5 3 4
104,49095 3
4
119,86972 4 0,03 0,00046319 0,05704323 3,9271640
1 5 5 6
133,13380 5
3
144,86523 6 0,04 0,00123552 0,09780908 4,1592238
9 8 9 4
155,4115 7
165 8 0,05 0,00242108 0,13914224 4,0629762
3 8 5
173,78933 9
9
181,89540 10 0,06 0,00401107 0,17832382 3,7422876
3 7 3 9
189,40593 11
8
196,38928 12 0,07 0,00597400 0,21346900 3,2653533
2 5 1 3
202,89987 13
5
208,98190 14 0,08 0,00826260 0,24328333 2,6837532
7 2 4 9
214,67180 15
6
220 16 0,09 0,01081905 0,26693526 2,0391582
5 1 9
224,99218 17
7
229,67028 18 0,1 0,01357919 0,28397283 1,3661204
4 1 4
234,05311 19
2
238,15698 20 0,11 0,01647595 0,29425979 0,6933778
6 4 2 4
241,99609 21
4
245,58285 22 0,12 0,01944225 0,29792122 0,0445450
3 9 3 9
248,92817 23
1
252,04166 24 0,13 0,02241338 0,29529460 -
3 3 3 0,5614511
254,93182 25 8
3
257,60616 26 0,14 0,02532884 0,28688467 -
7 7 4 1,1100311
260,07135 27 3
4
262,33328 28 0,15 0,02813388 0,27332158 -
2 1 1,5906413
264,39716 29 1
2
266,26760 30 0,16 0,03078048 0,25532224 -
5 4 4 1,9964226
267,94866 31 9
269,44387 32 0,17 0,03322810 0,23365517 -
2 2 1 2,3238557
270,75631 33 1
9
271,88864 34 0,18 0,03544399 0,20910894 -
9 8 2,5723820
272,84310 35 3
4
273,62154 36 0,19 0,03740330 0,18246455 -
5 1 1 2,7440083
274,22547 37 2
2
274,65603 38 0,2 0,03908881 0,15447180 -
5 4 3 2,8429006
274,91404 39 5
9
275 40 0,21 0,04049059 0,12582973 -
7 2,8749786
275 41 2
275 42 0,22 0,04160549 0,09721561 -
3 6 2,8343046
275 43 6
275 44 0,23 0,04243758 0,06939463 -
6 7 2,7184042
275 45 2
275 46 0,24 0,04299836 0,04305638 -
3 6 2,5398827
275 47 7
275 48 0,25 0,04330556 0,018763 -
6 2,3115569
275 49 3
275 50 0,26 0,04338192 -0,00305096 -
2 2,0460653
275 51 7
275 52 0,27 0,04325388 -0,02207502 -1,7555329
 275 53
275 54 0,28 0,04295039 -0,03811618 -
6 1,4512908
275 55 3
275 56 0,29 0,04250181 -0,05108983 -
275 57 1,1436543
7
275 58 0,3 0,04193880 -0,06100874 -
8 0,8417564
275 59 7
275 60 0,31 0,04129151 -0,06797058 -
4 0,5534356
268,125 61 3
261,25 62 0,32 0,04057979 -0,07478556 -
2 0,8006900
254,375 63 9
247,5 64 0,33 0,03978844 -0,08380619 -
9 0,9946256
240,625 65 8
233,75 66 0,34 0,03889808 -0,09450366 -
2 1,1363581
226,875 67 8
220 68 0,35 0,03789449 -0,10636719 -
6 1,2283435
213,125 69 2
206,25 70 0,36 0,03676846 -0,1189163 -
2 1,2741420
199,375 71 7
192,5 72 0,37 0,03551536 -0,13171063 -
5 1,2781809
185,625 73 4
178,75 74 0,38 0,03413475 -0,14435745 -
6 1,2455212
171,875 75 4
165 76 0,39 0,03262985 -0,15651688 -
5 1,1816360
158,125 77 9
151,25 78 0,4 0,03100700 -0,16790498 -
5 1,0922046
144,375 79 4
137,5 80 0,41 0,02927509 -0,17829483 -
9 0,9829256
130,625 81 5
123,75 82 0,42 0,02744501 -0,18751591 -0,8593537
9
116,875 83
110 84 0,43 0,02552907 -0,19545193 -
9 0,7267597
103,125 85 2
96,25 86 0,44 0,02354049 -0,2020374 -
 6 0,5900169
89,375 87 1
82,5 88 0,45 0,02149290 -0,20725321 -
9 0,4535120
75,625 89 1
68,75 90 0,46 0,01939993 -0,21112138 -
5 0,3210813
61,875 91 5
55 92 0,47 0,01727479 -0,2136994 -
1 0,1959704
48,125 93 1
41,25 94 0,48 0,01512996 -0,21507414 -
5 0,0808151
34,375 95 6
27,5 96 0,49 0,01297695 -0,21535581 0,0223570
7 5
20,625 97
13,75 98 0,5 0,01082607 -0,2146719 0,1121095
1 1
6,875 99
0 100 0,51 0,00868627 -0,21316143 0,1875636
7
0
0 0,52 0,00657402 -0,20832877 0,7638701
9 8
0
0 0,53 0,00453738 -0,19819518 1,2471196
6 9
0
0 0,54 0,00262465 -0,18370462 1,6351595
5 2
0
0 0,55 0,00087464 -0,16580742 1,9288001
7 6
0
0 0,56 - -0,14543266 2,1314250
0,00068324 4
4
0
0 0,57 -0,00202871 -0,12346438 2,2485758
0 4
0 0,58 - -0,10072222 2,2875280
0,00314997 4
2
0
0 0,59 - -0,07794616 2,2568700
0,00404306 3
5
0
0 0,6 - -0,0557854 2,1660971
 0,00471097 4
3
0
4,3387E-27 0,61 - -0,03479129 2,0252303
0,00516268 8
9
0,55
2,2 0,62 - -0,01512876 1,9289856
0,00541150 7
2
4,95
8,8 0,63 - 0,00419735 1,9554231
0,00546639 5 4
3
13,75
19,8 0,64 - 0,02433526 2,0885895
0,00532485 4 2
2
26,95
35,2 0,65 - 0,04626787 2,3114861
0,00497370 2 6
5
44,55
55 0,66 - 0,07080544 2,6066799
0,00439080 2 4
9
66,55
79,2 0,67 - 0,09858405 2,9568508
0,00354679 7
1
92,95
107,8 0,68 - 0,13006914 3,3452685
0,00240677 9 7
3
123,75
140,8 0,69 - 0,16556354 3,7561929
0,00093204 3 7
6
158,95
178,2 0,7 0,00091836 0,20521896 4,1751973
2 3 4
198,55
220 0,71 0,00318624 0,24905049 4,5894141
3 4 2
216,83279
7
213,62146 0,72 0,00587958 0,28699773 3,0143117
5 1 5
210,36458
3
207,06183 0,73 0,00887522 0,30968217 1,5422440
3 3 2
203,71288
6
200,31742 0,74 0,01202620 0,31828407 0,2019941
9 4 7 7
196,87513
9
193,38569 0,75 0,01519867 0,31423521 -
5 8 0,9849152
189,84877 1
7
186,26406 0,76 0,01827407 0,29914719 -
3 3 6 2,0039715
182,63123 5
2
178,94996 0,77 0,02115054 0,27474359 -
5 9 1 2,8472144
175,21994 7
171,44083 0,78 0,02374379 0,24279719 -
6 8 2 3,5126491
167,61233 8
3
163,73411 0,79 0,02598724 0,20507373 -
7 4 4,0035636
159,80584 3
5
155,82722 0,8 0,02783173 0,16328275 -4,3277807
8 7
151,79791
1
147,7176 0,81 0,02924475 0,11903615 -
1 7 4,4968745
143,58596 9
5
139,40268 0,82 0,03020925 0,07381464 -
4 3 1 4,5253781
135,16743 5
9
130,87990 0,83 0,03072225 0,02894204 -
6 5 8 4,4300046
126,53976 2
7
122,1467 0,84 0,03079313 -0,0144327 -
9 4,2289040
117,70038 8
4
113,20049 0,85 0,03044185 -0,05534644 -
9 9 3,9409713
108,64672 6
4
104,03873 0,86 0,02969704 -0,09302607 -
7 6 3,5852186
99,376219 9
4
94,658848 0,87 0,02859412 -0,12688683 -
8 3,1802226
89,886304 4
9
85,058266 0,88 0,02717343 -0,15652556 -
9 4 2,7436521
80,174414 1
1
75,234425 0,89 0,02547850 -0,18170958 -
6 6 2,2918803
70,237980 6
7
65,184758 0,9 0,02355439 -0,20236219 -
7 1,8396819
60,074438 8
7
54,9067 0,91 0,02144619 -0,21854547 -
8 1,4000125
49,681221 2
8
44,397683 0,92 0,01919780 -0,23044137 -0,9838671
3 5
39,055763
8
33,655142 0,93 0,01685075 -0,23833165 -
4 0,6002119
28,195498 5
4
22,676511 0,94 0,01444334 -0,24257755 -
1 1 0,2559819
17,097859
6
11,459223 0,95 0,01200996 -0,24359974 0,0438639
2 1 8
5,7602810
9
0,0007125 0,96 0,00958056 -0,24185898 0,2962385
7 1
0
0 0,97 0,00718793 -0,23559847 0,9386124
4 5
0
0 0,98 0,00488826 -0,22344176 1,4748182
8 6
0
0 0,99 0,00273517 -0,20646382 1,9028028
 5 6
0
0 1 0,00077146 -0,18574299 2,2238580
4 9
0
- -0,16232995 2,4421439
0,00097072 2
2
0
-
0,00246950
1
4) Courbes décrites par le déplacement, la vitesse et l’accélération sur le temps de simulation

déplacement u
vitesse v
0.05
0.4

0.04
0.3

0.03 0.2

0.02 0.1

0
0.01 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 1.05

-0.1
0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 1.05
-0.2
-0.01
-0.3
 accélerati on a
6

0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 1.05

-2

-4

-6
B.III. Superposition des courbes des déplacements sur une même courbe et interprétation

0.05 dep (runge cuta) dep (acc lin)

0.045

0.04

0.035

0.03

0.025

0.02

0.015

0.01

0.005

0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 1.05
-0.005

-0.01
0.35
0.3
0.25
0.2
0.15
0.1
0.05
0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 1.05
-0.05
-0.1
-0.15
-0.2
-0.25
-0.3

vitesse (runge cuta) vitesse (acc lin)

 5.5
5
4.5
4
3.5
3
2.5
2
1.5
1
0.5
0
-0.5 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 1.05
-1
-1.5
-2
-2.5
-3
-3.5
-4
-4.5
-5

ACC (runge cuta) ACC(acc lin)

 B.IV. Partie 2 : La structure a un comportement élasto-plastique

1.a. Première méthode : Accélération linéaire

1.a. Base théorique

Lorsque la structure a un comportement élasto-plastique, l’équation du mouvement peut se

mettre sur la forme : F I ( t )+ F d ( t )+ F s (t)=F (t) où F I ( t ), F d ( t ) , F s ( t ) et F (t) représente la force
d’inertie, la force d’inertie, la force de rappel et le chargement extérieur. Cette relation étant
quelque soit le temps, pour un incrément de temps ∆t, on a la relation suivante :

F I ( t+ ∆ t ) + F d ( t+ ∆ t )+ F s ( t+ ∆ t)=F (t+ ∆ t) Ainsi on a :

m ∆ ÿ (t)

{
F I ( t +∆ t )−F I ( t )=m ∆ ÿ (t)
F ( t +∆ t )−F ( t )=C (t)∆ y (t) ˙
d d
F s ( t +∆ t )−F s ( t )=K ( t ) y (t)

˙ ∆ R ( t )=∆ P (t). A partir de celle

L’équation du mouvement devient : m ∆ ÿ ( t ) +C ∆ y (t)+
relation, on développe l’algorithme de calcul à partir de cette relation.

1.b. Organigramme
i. Pour la méthode de l’accélération de l’accélération linéaire, on a les
1 1
coefficients suivants : β= et γ =
6 2
˙
P0−C y 0−R 0
ii. Initialisation : ÿ 0=
∆t
iii. Sélectionner le pas de variation temporelle : ∆ t

iv. Calculer les constantes : a=

m
+
β∆t β
γC
, b=
m
2β
+ ∆t
γ
2β
−1 C( )
v. Calculer les variables à chaque étape et à chaque itération
~
∆ P i=∆ Pi +a ẏ i +b ÿ i
vi. Déterminer la rigidité tangente : k i
γ m
vii. Calculer la rigidité finale : ǩ i=k i + C+
β ∆t β ∆t
2
 ~
∆ Pi
viii. Déterminer ∆ y i =
ki

ix. Calculer ∆ ẏ i =
γ
β ∆t
γ
(
∆ y i− ẏ i +∆ t 1−
β
γ
)ÿ
2β i
1 1 1
x. ∆ ÿ i = ∆ y i− ẏ i− ÿ
β ∆t
2
β ∆t 2β i
xi. Ensuite itérer pour les prochaines valeurs par la relation :
 y i+ 1= y i˙+ ¿ ∆ ẏ i ¿
 y i+ 1= y i+ ∆ y i
¨
 y i+ 1=¿ y i+¨∆˙y i ¿
 1.c.Valeurs prises par le déplacement, la vitesse et l’accélération

i t F(t) fs y(t) y'(t) y''(t)

0 0 0 0 0 0 0
1 0,01 309127, 311,906 0,00010 0,03245 6,49150
4 9 8 8 2
2 0,02 431531 2262,57 0,00078 0,26790 5,77185
9 5 5 7
3 0,03 521514, 10862,5 0,00376 0,47454 4,59143
9 8 8 2
4 0,04 594000 25233,8 0,00875 0,63823 3,51313
1 3 3 8
5 0,05 654823, 44159,2 0,01531 0,76320 2,63410
4 6 8 7 3
6 0,06 707001, 66554,4 0,02308 0,85674 1,94071
4 6 4
7 0,07 752334, 66554,4 0,03175 0,92548 1,95986
9 3 3 4
8 0,08 792000 66554,4 0,04111 0,99691 1,81193
5 5 9
9 0,09 826813 66554,4 0,05118 1,06275 1,63611
3 9 7
10 0,1 857365, 66554,4 0,0619 1,12212 1,45879
1 1 9
11 0,11 884098, 66554,4 0,0732 1,17499 1,29002
3 7 3
12 0,12 907350 66554,4 0,08501 1,22171 1,13260
9 4 6
13 0,13 927382, 66554,4 0,09729 1,26269 0,98636
2 8 1 5
14 0,14 944399, 66554,4 0,10997 1,29833 0,85000
8 8 1 8
15 0,15 958563, 66554,4 0,12300 1,32899 0,72196
4 6 2
16 0,16 969997, 66554,4 0,13633 1,35497 0,60068
9 5 1 6
17 0,17 978799, 66554,4 0,14991 1,37651 0,48480
1 6 2 5
18 0,18 985037, 66554,4 0,16370 1,39380 0,37309
6 7 6 8
19 0,19 988761, 66554,4 0,17766 1,40700 0,26449
7 4 2 5
20 0,2 990000 66554,4 0,19174 1,41621 0,15803
7 7
21 0,21 990000 66554,4 0,20591 1,42163 0,07887
7 3 5
22 0,22 990000 66554,4 0,22013 1,42434 0,03936
7 6
23 0,23 990000 66554,4 0,23438 1,42569 0,01964
2 1 7
24 0,24 990000 66554,4 0,24864 1,42636 0,00980
6 6
25 0,25 990000 66554,4 0,26290 1,42670 0,00489
4 2 4
26 0,26 990000 66554,4 0,27717 1,42687 0,00244
1 3
27 0,27 990000 66554,4 0,29144 1,42695 0,00121
4 9
28 0,28 990000 66554,4 0,30571 1,42699 0,00060
6 8
29 0,29 990000 66554,4 0,31998 1,42701 0,00030
7 4
30 0,3 990000 66554,4 0,33425 1,42702 0,00015
7 2
31 0,31 940500 66554,4 0,34850 1,42183 -
 3 1,04045
32 0,32 891000 66554,4 0,36265 1,38091 -1,5598
4 9
33 0,33 841500 66554,4 0,37637 1,32218 -
1 4 1,81901
34 0,34 792000 66554,4 0,38948 1,25455 -
9 4 1,94837
35 0,35 742500 66554,4 0,40192 1,18248 -
4 3 2,01294
36 0,36 693000 66554,4 0,41363 1,10819 -
6 7 2,04516
37 0,37 643500 66554,4 0,42460 1,03280 -
3 4 2,06125
38 0,38 594000 66554,4 0,43481 0,95686 -
5 2,06927
39 0,39 544500 66554,4 0,44426 0,88064 -
8 2,07328
40 0,4 495000 66554,4 0,45295 0,80428 -
8 3 2,07528
41 0,41 445500 66554,4 0,46088 0,72785 -
4 7 2,07628
42 0,42 396000 66554,4 0,46804 0,65139 -
6 7 2,07678
43 0,43 346500 66554,4 0,47444 0,57491 -
4 9 2,07702
44 0,44 297000 66554,4 0,48007 0,49843 -
6 3 2,07715
45 0,45 247500 66554,4 0,48494 0,42194 -
4 3 2,07721
46 0,46 198000 66554,4 0,48904 0,34545 -
7 1 2,07724
47 0,47 148500 66554,4 0,49238 0,26895 -
5 8 2,07726
48 0,48 99000 66554,4 0,49495 0,19246 -
 8 4 2,07726
49 0,49 49500 66554,4 0,49676 0,11597 -
6 2,07727
50 0,5 0 66554,4 0,49781 0,03947 -
6 2,07727
51 0,51 0 66554,4 0,49810 - -
5 0,03182 1,03675
52 0,52 0 66554,4 0,49773 -0,0674 -
8 0,51743
53 0,53 0 66554,4 0,49703 - -
9 0,08516 0,25825
54 0,54 0 66554,4 0,49617 - -
5 0,09402 0,12889
55 0,55 0 66554,4 0,49522 - -
9 0,09844 0,06433
56 0,56 0 66554,4 0,49424 - -
1 0,10065 0,03211
57 0,57 0 66554,4 0,49323 - -
3 0,10175 0,01602
58 0,58 0 66554,4 0,49221 -0,1023 -0,008
5
59 0,59 0 66554,4 0,49119 - -
1 0,10258 0,00399
60 0,6 0 66554,4 0,49016 - -
6 0,10271 0,00199
61 0,61 7920 66554,4 0,48914 - 0,16548
1 0,10195 9
62 0,62 31680 66554,4 0,48813 - 0,58204
8 0,09377 5
63 0,63 71280 66554,4 0,48724 - 1,12291
2 0,06963 2
64 0,64 126720 66554,4 0,48661 - 1,72582
8 0,02527 1
65 0,65 198000 66554,4 0,48647 0,04145 2,35969
 3 3 6
66 0,66 285120 66554,4 0,48703 0,13159 3,00902
1 4 5
67 0,67 388080 66554,4 0,48852 0,24568 3,66606
6 5 7
68 0,68 506880 66554,4 0,4912 0,38399 4,32695
8
69 0,69 641520 66554,4 0,49529 0,54664 4,98977
4 2 1
70 0,7 792000 66554,4 0,50105 0,73370 5,65354
1 6 3
71 0,71 769037, 66554,4 0,50865 0,92532 2,33894
3 1 8
72 0,72 745422, 66554,4 0,51800 1,00311 0,67095
6 6 2 4
73 0,73 721142, 66554,4 0,52806 1,02358 -
7 1 7 0,17551
74 0,74 696188, 66554,4 0,53828 1,01494 -
5 1 0,61215
75 0,75 670550, 66554,4 0,54839 0,99123 -
6 1 7 0,84444
76 0,76 644219, 66554,4 0,55825 0,95948 -
9 4 8 0,97494
77 0,77 617187 66554,4 0,56779 0,92318 -
3 7 1,05484
78 0,78 589442, 66554,4 0,57696 0,88406 -
8 4 9 1,10966
79 0,79 560978 66554,4 0,58574 0,84299 -
2 4 1,15217
80 0,8 531783, 66554,4 0,59410 0,80038 -
4 7 3 1,18873
81 0,81 501849, 66554,4 0,60204 0,75644 -
7 3 1,22251
82 0,82 471167, 66554,4 0,60953 0,71125 -1,2551
 7 9 9
83 0,83 439728, 66554,4 0,61658 0,66488 -
1 1,28729
84 0,84 407521, 66554,4 0,62315 0,61731 -
8 6 5 1,31947
85 0,85 374539, 66554,4 0,62925 0,56856 -
5 5 4 1,35185
86 0,86 340771, 66554,4 0,63486 0,51862 -
9 4 1,38451
87 0,87 306209, 66554,4 0,63997 0,46747 -
8 2 1 1,41752
88 0,88 270843, 66554,4 0,64456 0,41510 -
9 7 5 1,45088
89 0,89 234665, 66554,4 0,64863 0,36150 -
1 6 8 1,48463
90 0,9 197664, 66554,4 0,65216 0,30666 -
1 7 7 1,51875
91 0,91 159831, 66554,4 0,65514 0,25056 -
7 8 7 1,55326
92 0,92 121158, 66554,4 0,65756 0,19319 -
5 6 5 1,58815
93 0,93 81635,4 66554,4 0,65940 0,13453 -
4 9 6 1,62344
94 0,94 41253,2 66554,4 0,66066 0,07457 -
2 5 1,65911
95 0,95 0 66554,4 0,66131 0,01329 -
4 9 1,69522
96 0,96 0 66554,4 0,66136 - -
6 0,04488 0,84607
97 0,97 0 66554,4 0,66087 - -
7 0,07392 0,42227
98 0,98 0 66554,4 0,66011 - -
8 0,08841 0,21075
99 0,99 0 66554,4 0,65922 - -
 4 0,09564 0,10518
100 1 0 66554,4 0,65826 - -0,0525
2 0,09925
 1.d. Tracer les courbes décrites par le déplacement, la vitesse et l’accélération

DEPLACEMENT Y(T)
0.7

0.6

0.5

0.4

0.3

0.2

0.1

0
0 0.2 0.4 0.6 0.8 1 1.2
 VITESSE y'(t)
1.6

1.4

1.2

0.8

0.6

0.4

0.2

0
0 0.2 0.4 0.6 0.8 1 1.2
-0.2
 ACCELERATION y''(t)
7

0
0 0.2 0.4 0.6 0.8 1 1.2
-1

-2

-3
 2. Deuxième méthode : Méthode de Runge Kutta d’ordre 4

2. a. Base théorique

Les méthodes de Runge-Kutta sont des méthodes d'analyse

numérique d'approximation de solutions d'équations différentielles. Elles ont été nommées
ainsi en l'honneur des mathématiciens Carl Runge et Martin Wilhelm Kutta lesquels
élaborèrent la méthode en 1901.

Ces méthodes reposent sur le principe de l'itération, c'est-à-dire qu'une première

estimation de la solution est utilisée pour calculer une seconde estimation, plus précise, et
ainsi de suite.

Il en existe plusieurs mais celle que nous traiterons est la méthode d’ordre 4 qui est
une amélioration de la méthode de base. Elle possède une étape supplémentaire dans le
raffinement du calcul.

Les méthodes de Runge Kutta s’écrivent de manière générale

K 1=h¿

K 2=h¿

………… …….
K n=h ¿
y k +1= y k +h ¿

Les coefficients sont déterminés à fin que l’ordre soit le plus élevé possible. On note A
la matrice de coefficients (α i , j), ᴦ le vecteur des coefficients γi et Ө le vecteur coefficients de
θi

Une représentation en forme de tableau des équation donne

γ1 γ2 ……… γn
θ1 α 1 ,1 α 1 ,2 …………. α 1 ,n
θ2 α 2 ,1 α 2 ,2 …………. α 2 ,n
………….. ………… ……………. ……...
θn α n ,1 α n ,2 ……………. α n ,n
 Avec cette représentation, la méthode Runge-Kutta explicite à deux points qui est
d’ordre 4 est représentée par le tableau suivant

1/6 1/3 1/3 1/6

0 0 0 0 0
1/2 1/2 0 0 0
1/2 0 1/2 0 0
1 0 0 1 0

2. b. Organigramme
Soit f s la force de rappel de notre structure
A l’étape i+1, (f ¿¿ s )i +1=(f ¿¿ s)i +k i ui ¿¿
Dans le domaine élastique k =k
Dans le domaine élastoplasique k =0 et fs=r(i)

 Calculs initiaux
iv. Calculez k, c, m.
v. Initialiser les variables
(f ¿¿ s )0=k u0 ¿
¿
v 0=u̇0
˙1
v˙0= ¿¿
m
vi. Sélectionnez un pas de temps approprié ∆t.

 Pour chaque pas de temps

v.Calcul au début de l'intervalle de temps.

t=t i
x=u=ui
¿
 v 1=v i
1
v̇ 1=v i ¿
m

vi.Calcul au premier point médian de l'intervalle de temps

Δt
t=t i+
2
vi Δ t
x=u Δt =ui +
i+
2
2
¿

v˙1 Δ t
v 2=v Δt =v i +
i+
2
2
1
v̇ 2=v̇ Δt = ¿
i+
2
m

vii. Calcul au deuxième point médian de l'intervalle de temps

Δt
t=t i+
2
Δt
x=u Δt =ui +v 2
i+
2
2
¿
Δt
v 3=v Δt =v i + v˙2
i+
2
2
1
v˙3= ¿
m
viii.Calculer les déplacements à la fin de l'intervalle de temps
t=t i+ Δ t
x=ui + Δt =ui + v 3 Δ t
¿
v 4=v i+ Δ t=v i+ v˙3 Δ t
 1
v˙4= ¿
m
ix.Calculer le déplacement et la vitesse à la fin de l'intervalle de temps
Δt
ui +1=ui +
6 1
[ v +2 v 2 +2 v 3 + v 4 ]
Δt
v i+1=v i+
6 1
[ v˙ + 2 v˙2+ 2 v˙3+ v˙4 ]
 2. d. Valeurs prises par le déplacement, la vitesse et l’accélération sur le
temps de simulation

F(t) i t y(t) y'(t) y''(t)

0 0 0,01 0 0 0
61,1063366 1
0,02103567
7,92327E-05 3,15213414
85,8687225 2 0,02 3
104,490954 3
0,05704323
0,000463195 3,92716406
119,869721 4 0,03 5
133,133803 5
0,09780908
0,001235528 4,15922384
144,865239 6 0,04 9
155,4115 7
0,13914224
0,002421083 4,06297625
165 8 0,05 8
173,789339 9
0,17832382
0,004011077 3,74228769
181,895403 10 0,06 3
189,405938 11
0,21346900
0,005974005 3,26535333
196,389282 12 0,07 1
202,899875 13
0,24328333
0,008262602 2,68375329
208,981907 14 0,08 4
214,671806 15
0,26693526
0,010819055 2,03915829
220 16 0,09 1
224,992187 17
0,28397283
0,013579194 0,38905615
229,67028 18 0,1 1
234,053112 19
0,28928628
0,016443228 0,66016922
238,156986 20 0,11 4
241,996094 21
0,29693152
0,01937267 0,85775127
245,582853 22 0,12 6
248,928171 23
0,30623864
0,022387381 0,9944393
252,041663 24 0,13 3
254,931823 25
0,31665183
0,025501117 1,08051417
257,606167 26 0,14 8
260,071354 27
262,33328 28 0,15 0,028722553 0,32770843 1,12438762
 2
264,397162 29
266,267605 30 0,16 0,032056136 0,33902207 1,13296613
267,94866 31
0,35026912
0,035502768 1,11192905
269,443872 32 0,17 7
270,756319 33
0,36117755
0,039060387 1,06594488
271,888649 34 0,18 6
272,843104 35
0,37151766
0,042724424 0,99884165
273,621545 36 0,19 4
274,225472 37
0,38109441
0,046488196 0,91374214
274,656035 38 0,2 6
274,914049 39
0,38974093
0,050343214 0,81317149
275 40 0,21 5
275 41
275 42 0,22 0,054279548 0,39735761 0,71238977
275 43
275 44 0,23 0,058287224 0,4040303 0,62409859
275 45
275 46 0,24 0,062357401 0,409876 0,54674992
275 47
0,41499720
0,066482333 0,47898759
275 48 0,25 3
275 49
0,41948370
0,070655233 0,41962349
275 50 0,26 2
275 51
0,42341415
0,074870156 0,36761677
275 52 0,27 9
275 53
0,42685748
0,079121895 0,32205559
275 54 0,28 9
275 55
0,42987406
0,083405886 0,28214111
275 56 0,29 4
275 57
0,43251677
0,087718132 0,2471735
275 58 0,3 5
275 59
0,43483195
0,092055131 0,21653965
275 60 0,31 7
268,125 61
0,43421235 -
0,096404902
261,25 62 0,32 2 0,32842282
254,375 63
 0,42848828 -
0,10072239
247,5 64 0,33 5 0,80584449
240,625 65
0,41829238 -
0,104959785
233,75 66 0,34 7 1,22409613
226,875 67
0,40417888 -
0,1090752
220 68 0,35 2 1,59051107
213,125 69
0,38663330 -
0,11303194
206,25 70 0,36 3 1,91151379
199,375 71
0,36608101 -
0,116797859
192,5 72 0,37 2 2,19273252
185,625 73
0,34289464 -
0,120344794
178,75 74 0,38 9 2,43909795
171,875 75
0,31740067 -
0,123648072
165 76 0,39 3 2,65492968
158,125 77
0,28988508 -
0,126686079
151,25 78 0,4 1 2,84401195
144,375 79
0,26059842 -
0,129439879
137,5 80 0,41 6 3,00966001
130,625 81
0,22976020 -
0,131892884
123,75 82 0,42 7 3,15477821
116,875 83
-
0,134030559 0,19756272
110 84 0,43 3,28191095
103,125 85
0,16417442 -
0,135840175
96,25 86 0,44 7 3,39328729
89,375 87
0,12974291 -
0,137310576
82,5 88 0,45 4 3,49086002
75,625 89
0,09439747 -
0,138431991
68,75 90 0,46 3 3,57633991
61,875 91
0,05825137 -
0,139195861
55 92 0,47 4 3,65122572
48,125 93
0,02140384 -
0,139594685
41,25 94 0,48 8 3,71683044
34,375 95
-
0,139621893 -0,01605817
27,5 96 0,49 3,77430432
 20,625 97
-
0,139271729 -0,05405853
13,75 98 0,5 3,82465509
6,875 99
-
0,138539152 -0,09253051
0 100 0,51 3,86876556
0
-
0,137428659 -0,1287678
0 0,52 3,38928384
0
-
0,135978743 -0,16051396
0 0,53 2,96922747
0
-
0,134231474 -0,18832561
0 0,54 2,60123146
0
-
0,132223703 -0,21269038
0 0,55 2,27884363
0
-
0,129987716 -0,23403546
0 0,56 1,99641146
0
-
0,127551798 -0,2527351
0 0,57 1,74898298
0
-
0,124940727 -0,26911718
0 0,58 1,53221994
0
-
0,122176212 -0,28346891
0 0,59 1,34232179
0
-
0,119277269 -0,29604194
0 0,6 1,17595896
0
-
0,116260559 -0,30705671
4,3387E-27 0,61 1,03021459
0,55
-
0,113141391 -0,31642077
2,2 0,62 0,81780628
4,95
-
0,109941603 -0,32294796
8,8 0,63 0,46592306
13,75
19,8 0,64 0,106696239 -0,32533188 0,00814904
26,95
35,2 0,65 0,103452581 -0,32242803 0,58926634
44,55
55 0,66 0,100268631 -0,31323376 1,26416201
66,55
 79,2 0,67 0,097211782 -0,29687067 2,02121349
92,95
107,8 0,68 0,094357655 -0,27256926 2,85023866
123,75
140,8 0,69 0,091789078 -0,23965536 3,74231734
158,95
178,2 0,7 0,089595194 -0,19753838 4,68963489
198,55
220 0,71 0,087870679 -0,14570092 5,68534521
216,832797
213,621465 0,72 0,086681707 -0,09367118 4,74029455
210,364583
207,06183 0,73 0,085967679 -0,05052804 3,90554319
203,712886
200,317429 0,74 0,085645006 -0,01523921 3,16728468
196,875139
0,01309863
0,085639761 2,51346334
193,385695 0,75 6
189,848777
0,03527589
0,085886474 1,9335176
186,264063 0,76 2
182,631232
0,05198402
0,086327076 1,41819487
178,949965 0,77 9
175,21994
0,06382785
0,086909965 0,9593893
171,440836 0,78 3
167,612333
0,07133624
0,087589203 0,54999964
163,73411 0,79 5
159,805845
0,07497157
0,0883238 0,18380476
155,82722 0,8 1
151,797911
0,07513792 -
0,089077089
147,7176 0,81 5 0,14464547
143,585965
0,07218835 -
0,089816187
139,402684 0,82 6 0,44012569
135,167439
0,06643118 -
0,090511511
130,879906 0,83 4 0,70681878
126,539767
0,05813555 -
0,091136362
122,1467 0,84 6 0,94838923
117,700384
0,04753629 -
0,091666555
113,200499 0,85 2 1,16804734
108,646724
 0,03483814 -
0,092080102
104,038737 0,86 1 1,36860557
99,3762194
0,02021951 -
0,092356926
94,6588488 0,87 1 1,55252778
89,8863049
0,00383573 -
0,092478617
85,0582669 0,88 2 1,72197248
80,1744141
-
0,092428214 -0,01417809
75,2344256 0,89 1,87883067
70,2379807
-
0,09219002 -0,03370454
65,1847587 0,9 2,02475898
60,0744387
-
0,091749428 -0,05464176
54,9067 0,91 2,16120871
49,6812218
-
0,091092783 -0,07690144
44,3976833 0,92 2,28945131
39,0557638
-
0,090207251 -0,10040721
33,6551424 0,93 2,41060065
28,1954984
-
0,08908071 -0,12509313
22,6765111 0,94 2,52563253
17,0978596
-
0,087701649 -0,15090239
11,4592232 0,95 2,63540184
5,76028109
-
0,086059085 -0,17778619
0,0007125 0,96 2,74065749
0
-
0,084150033 -0,20345713
0 0,97 2,40101534
0
-
0,082000531 -0,22594655
0 0,98 2,10344162
0
-
0,079640379 -0,24564871
0 0,99 1,84274818
0
0 1 0,077095683 -0,26290905 -1,6143642
0
-
0,074389317 -0,2780302
1,41428536
0
0,071541317
2. e. Courbes décrites par le déplacement, la vitesse et l’accélération sur le temps de simulation

Déplacement u

0.16

0.14

0.12

0.1

0.08

0.06

0.04

0.02

0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 1.05
 vitesse v
0.5

0.4

0.3

0.2

0.1

0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 1.05

-0.1

-0.2

-0.3

-0.4
 accélération a

0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 1.05

-2

-4

-6