Maria Bayard Dühring - Solid Mechanics
Maria Bayard Dühring - Solid Mechanics
Maria Bayard Dühring - Solid Mechanics
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.