07.01.2015 Views

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 ...

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!