Interaction sol-structure en milieu stochastique - Bibliothèque Ecole ...
Interaction sol-structure en milieu stochastique - Bibliothèque Ecole ...
Interaction sol-structure en milieu stochastique - Bibliothèque Ecole ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
INTERACTION SOL-STRUCrURE EN MILIEU STOCHASTIQUE<br />
de simulations de <strong>milieu</strong> aléatoire associé aux fluctuations de la rigidité.<br />
Ainsi, à partir des équations bilan et sans approximation dans l'expression analytique initiale<br />
des grandeurs OL et<br />
G (voir équations 2.13, 2.14), on va calculer les mom<strong>en</strong>ts statistiques<br />
du tassem<strong>en</strong>t et du dévers, au moy<strong>en</strong> d'un grand nombre de réalisations du <strong>milieu</strong> (1000<br />
réalisations de <strong>milieu</strong>).<br />
Pour ce faire, on utilisera les deux grandeurs statistiques données:<br />
- le rayon de corrélation lié à la variable -<br />
k0<br />
- l'écart type a associé à cette même variable.<br />
Dans une première étape, à l'aide de ces différ<strong>en</strong>tes données précitées ci-dessus, on génère des<br />
<strong>milieu</strong>x aléatoires (Schuëller, 1987) pour représ<strong>en</strong>ter les fluctuations de rigidité adim<strong>en</strong>sionnelles,<br />
notées f(x);<br />
ainsi<br />
Ek(x)<br />
N<br />
f(x) = k0<br />
= acos(kx+p)<br />
(2.46)<br />
où (p est une phase aléatoire indép<strong>en</strong>dante, uniformém<strong>en</strong>t distribuée <strong>en</strong>tre O et 2 ir, et k représ<strong>en</strong>te<br />
le nombre d'onde. Les coeffici<strong>en</strong>ts a sont connus à partir de la d<strong>en</strong>sité spectrale de puissance<br />
(DSP). Cette dernière est donnée à partir de la transformée de Fourrier de la fonction<br />
d'autocorrélation.<br />
Pour une fonction d'autocorrélation expon<strong>en</strong>tielle définie par:<br />
II<br />
R(f () = E[f(x + )f(x)] = a2e r<br />
(2.47)<br />
la DSP est alors donnée par:<br />
Sff(kfl)<br />
a2r0<br />
ir(l +(kr0)2)<br />
(2.48)<br />
avec k = nb.k. La longueur du pas ik s'évalue à partir de la représ<strong>en</strong>tation de la DSP <strong>en</strong> estimant<br />
la valeur maximum de k. De plus, il est nécessaire que l'aire de la DSP permette de retrouver<br />
l'écart type des fluctuations.<br />
68