28.02.2013 Views

Introduction to Acoustics

Introduction to Acoustics

Introduction to Acoustics

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.

notes the value of this eigenmode at one particular point<br />

of the structure. For the structure vibrating in air, the<br />

displacement becomes<br />

ξ = φ10q1 + φ20q2 . (22.245)<br />

From (22.226) one can derive a first-order approximation<br />

for the bar displacement<br />

�<br />

� �<br />

�<br />

sC21<br />

ξ = φ10<br />

+ φ20 q20 + q10 .<br />

q10 + sC12<br />

q20<br />

D1<br />

D2<br />

(22.246)<br />

Operating Deflexion Shapes. Equation (22.246) can be<br />

rewritten<br />

�<br />

ξ = q10 φ10 + sC21<br />

�<br />

φ20<br />

D2<br />

�<br />

�<br />

sC12<br />

+ q20 φ20 + φ10<br />

D1<br />

with qi0 = φi0<br />

F . (22.247)<br />

mi Di<br />

In experiments on structures, sinusoidal excitation is<br />

often used. Imagine that we apply a sudden harmonic<br />

force F(t) = H(t)sinωt at time t = 0<strong>to</strong>thestructure,<br />

with excitation location and frequency such that q20 is<br />

negligible compared <strong>to</strong> q10. H(t) is the Heaviside function.<br />

In this case, the spatial pattern of the structure is<br />

given by<br />

sC21<br />

φ1 = φ10 + φ20 . (22.248)<br />

D2<br />

Because of the time dependence of the second term<br />

(through the Laplace variable s), the spatial shape<br />

evolves with time. Here, we can use the Laplace limit<br />

theorem, which states that the value of φ(t) as time tends<br />

<strong>to</strong> infinity is given by the product of sφ(s) ass tends <strong>to</strong><br />

zero. Since s2C12/D2 tends <strong>to</strong> zero as s tends <strong>to</strong> zero,<br />

the second term on the right-hand side of (22.248) vanishes<br />

after some time. Calculating the inverse Laplace<br />

transform shows that this decay time is of the order of<br />

the magnitude of the decay time of the second structural<br />

mode. After this transient regime, the spatial shape is<br />

nearly equal <strong>to</strong> the spatial shape in vacuo φ10.Inthemore<br />

general case, (22.247)showsthatFexcites both q10 and<br />

q20. After a transient regime, the bar displacement then<br />

finally converges <strong>to</strong>:<br />

� �<br />

φ10β1 φ20β2<br />

ξ(ω, x) = + . (22.249)<br />

D1(iω) D2(iω)<br />

The quantity between in square brackets is called the<br />

operating deflexion shape (ODS) ofthestructureat<br />

Structural <strong>Acoustics</strong> and Vibrations 22.5 Structural–Acoustic Coupling 931<br />

frequency ω. Since it is very difficult in practice <strong>to</strong> excite<br />

a single qi0, theODS describes the multi-mode<br />

shapes that are observed for sinusoidal excitation of<br />

structures.<br />

State-Space Formulation<br />

The transfer function formulation is convenient if the<br />

system is initially at rest, and for time-invariant systems.<br />

For other applications, such as sound control, it is<br />

useful <strong>to</strong> express the results in terms of state-space variables<br />

[22.26–29]. In most cases, the mechanical state of<br />

the system is given by the position and the velocity of<br />

the DOFs. Since the unloaded structure is described by<br />

its eigenmodes φn, all useful information about the state<br />

of the system is contained in the modal participation fac<strong>to</strong>rs<br />

qn and in their first derivatives with time ˙qn. This<br />

allows us <strong>to</strong> rewrite the equations for a 2-DOF coupled<br />

system as follows<br />

d<br />

dt<br />

where<br />

⎛<br />

⎜<br />

⎝<br />

X1<br />

X2<br />

X3<br />

X4<br />

⎞<br />

⎟<br />

⎠ =<br />

⎛<br />

0 0 1 0<br />

⎜ 0 0 0 1<br />

⎜<br />

⎝−ω2<br />

1 0 −2ζ1ω1 C12<br />

0 −ω2 ⎞<br />

⎟<br />

⎠<br />

2<br />

C21 −2ζ2ω2<br />

⎛ ⎞<br />

X1<br />

⎜ ⎟<br />

⎜X2⎟<br />

× ⎜ ⎟<br />

⎝ ⎠ +<br />

⎛ ⎞<br />

0<br />

⎜ ⎟<br />

⎜ 0 ⎟<br />

⎜ ⎟ F , (22.250)<br />

⎝ ⎠<br />

X3<br />

X4<br />

β1<br />

β2<br />

X1 = q1 ; X2 = q2 ; X3 = ˙q1 ; X4 = ˙q2 . (22.251)<br />

Equation (22.250) can then be formulated equivalently<br />

as<br />

˙X = AX + BF , (22.252)<br />

where F is the input and X is the state vec<strong>to</strong>r. The output<br />

Y depends on the investigated mechanical problem. If<br />

we decide, for example, <strong>to</strong> investigate the displacement,<br />

then we can write for the output<br />

� �t Y =<br />

X = Ɣ X . (22.253)<br />

φ1 φ2 00<br />

Equations (22.252)and(22.253) are general expressions<br />

for a linear system expressed in terms of state-space<br />

variables.<br />

Remark 1. The representation presented in (22.250) is<br />

not unique. Selecting, for example<br />

X1 = q1; X2 = ˙q1; X3 = q2; X4 = ˙q2 (22.254)<br />

leads <strong>to</strong> different values for A and B.<br />

Part G 22.5

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

Saved successfully!

Ooh no, something went wrong!