28.08.2013 Views

Maria Bayard Dühring - Solid Mechanics

Maria Bayard Dühring - Solid Mechanics

Maria Bayard Dühring - Solid Mechanics

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

z [m]<br />

3<br />

2.5<br />

2<br />

3<br />

y [m]<br />

2<br />

1<br />

by Helmholtz equation for a homogeneous medium. The material properties for air in the room are ra and ka<br />

and it is assumed that there is no damping effect in the air. The material on the boundaries is inhomogeneous<br />

and the optimized design is a distribution of the usual reflecting material described by r2 and k2 and an<br />

absorbing material with the properties r1 ¼ 3:04 kg m 3 and k1 ¼ 7:90 10 5 Nm 2 . The absorbing material<br />

has an absorption coefficient equal to 0.1 which could be realized in practice by a cork sheet with the thickness<br />

of a few millimeters. It is then convenient to use the new variables<br />

8<br />

8<br />

~r ¼ r<br />

¼<br />

ra ><<br />

>:<br />

r1 ra r2 ra absorbing;<br />

reflecting;<br />

~k ¼ k<br />

ka<br />

><<br />

¼<br />

>:<br />

With this rescaling the acoustic model for the problem takes the form<br />

n r^p ¼<br />

ARTICLE IN PRESS<br />

M.B. <strong>Dühring</strong> et al. / Journal of Sound and Vibration 317 (2008) 557–575 567<br />

2<br />

1<br />

0 0x [m]<br />

k1<br />

ka<br />

k2<br />

ka<br />

absorbing;<br />

reflecting:<br />

(19)<br />

r 2 ^p þ ~o 2 ^p ¼ 0 Helmholtz equation, (20)<br />

pffiffiffiffiffiffiffiffiffi<br />

i ~o kara ^p; b:c: for surface with impedance Z (21)<br />

ZðrÞ<br />

n r^p ¼ i ~o ffiffiffiffiffiffiffiffiffi p<br />

karaU; b:c: for pulsating surface. (22)<br />

The inhomogeneities pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi<br />

on the walls are described in the boundary condition (21) by the impedance boundary<br />

ZðrÞ ¼raca ~kðrÞ~rðrÞ.<br />

The impedance boundary condition is only strictly valid for plane waves of normal<br />

incidence but is used for simplicity in this work. The material interpolation functions for ~rðxÞ and ~kðxÞ must<br />

now satisfy the requirements<br />

8<br />

r1 >< ra ~rðxÞ ¼ r2 >:<br />

x ¼ 0;<br />

x ¼ 1;<br />

8<br />

k1<br />

>< ;<br />

ka<br />

~kðxÞ ¼ k2<br />

>: ;<br />

x ¼ 0;<br />

x ¼ 1;<br />

(23)<br />

r a<br />

3<br />

4<br />

30 40 50<br />

frequency, f [Hz]<br />

60 70<br />

and again interpolation functions in the inverse material properties are used<br />

~rðxÞ 1 ¼ r1 ra 1<br />

þ x<br />

r2 ra 1<br />

r1 ra !<br />

1<br />

, (24)<br />

objective function, Φ [dB]<br />

130<br />

120<br />

110<br />

100<br />

90<br />

80<br />

70<br />

60<br />

50<br />

40<br />

ka<br />

initial guess<br />

optimized design<br />

frequency borders<br />

Fig. 8. Results of the optimization for the frequency interval [40.5;45.5] Hz with four target frequencies: (a) the optimized design and (b)<br />

the frequency response for the initial design and the optimized design.

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

Saved successfully!

Ooh no, something went wrong!