22.11.2014 Views

Download full text - ELSA - Europa

Download full text - ELSA - Europa

Download full text - ELSA - Europa

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Multiple scales in frequency<br />

• Some large (CPU!) problems contain multiple frequency scales.<br />

• Impact raises high frequencies<br />

and non-linear effects locally.<br />

• However, the bulk structure has<br />

low-frequency vibratory<br />

behaviour (linear).<br />

• Fine FE discretization (to fit<br />

complex geometry) is too detailed<br />

to describe linear vibrations.<br />

• Use modal projection: replace FE unknowns (dofs) by less<br />

numerous ones, specific to the chosen frequency spectrum.<br />

129<br />

Multiple scales in frequency (2)<br />

• Domain decomposition allows to separate the frequency scales.<br />

• How can one suitably couple a linear, “modal” S/D with a nonlinear,<br />

“direct” one? One arrives to a form (S/D 1 is modal):<br />

⎡ ∆t ˆ<br />

∆t ˆ T ⎤ ⎡ ∆t<br />

ˆ ⎤<br />

⎢<br />

M1 0 − C1 1 1<br />

2 2 ⎥ ⎢<br />

Kα<br />

2 ⎥<br />

⎢ ⎥⎡<br />

α1<br />

⎤ ⎢ ⎥<br />

⎢ ∆t ∆t T<br />

∆t<br />

0 M2 − C ⎥⎢<br />

2<br />

U<br />

⎥<br />

2<br />

= ⎢ F ⎥<br />

2<br />

⎢ 2 2 ⎥⎢<br />

⎥ ⎢ 2 ⎥<br />

⎢<br />

t ˆ t<br />

⎥⎢⎣Λ<br />

⎥⎦<br />

⎢<br />

ˆ p n p n⎥<br />

⎢<br />

∆ ∆<br />

− C<br />

1 1 2 2<br />

1<br />

− C2<br />

0 ⎥ ⎢C α + C U<br />

⎥<br />

⎢⎣ 2 2<br />

⎥⎦ ⎢⎣ ⎥⎦<br />

• This coupling is totally independent<br />

from the nature of the base used for<br />

modal reduction (free choice!).<br />

• This is possible thanks to dual approach:<br />

interface dofs don’t need to remain “visible”<br />

after reduction.<br />

α1 generalized unknowns:<br />

⎧U<br />

1<br />

=Φ1α<br />

1<br />

⎪<br />

⎨U<br />

<br />

1<br />

=Φ <br />

1α<br />

1<br />

⎪<br />

U<br />

⎩ 1<br />

=Φ <br />

1α<br />

1<br />

linear behaviour:<br />

F = KU = KΦα<br />

1 1 1 1 1 1<br />

"projected" operators:<br />

ˆ T<br />

⎧M1 Φ1M1Φ1<br />

⎪ ˆ T<br />

⎨K1 Φ1K1Φ1<br />

⎪<br />

Cˆ<br />

⎪⎩ 1<br />

C1Φ<br />

1<br />

130<br />

65

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

Saved successfully!

Ooh no, something went wrong!