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.

the mass matrix is taken <strong>to</strong> be constant with time. In this<br />

case, (22.46) becomes<br />

with<br />

⎧<br />

M ¨X + KX = 0 (22.47)<br />

⎪⎨<br />

⎪⎩<br />

Mkj = M jk = ∂2 Ek<br />

∂ ˙qk∂ ˙q j<br />

Kkj = K jk = ∂2 V<br />

∂qk∂q j<br />

22.2.2 Eigenmodes<br />

and Eigenfrequencies<br />

. (22.48)<br />

Natural modes (or eigenmodes) are the sinusoidal (or<br />

harmonic) solutions of (22.47) with frequency ω in the<br />

absence of a driving term. As a consequence, the natural<br />

modes are solutions of the equation<br />

�<br />

− ω 2 �<br />

M + K X = 0 (22.49)<br />

with roots (or eigenvalue, or eigenfrequency) ωn solutions<br />

of the characteristic equation<br />

�<br />

det − ω 2 �<br />

M + K = 0 . (22.50)<br />

The eigenvec<strong>to</strong>r Φn associated <strong>to</strong> each eigenfrequency<br />

ωn is given by<br />

�<br />

− ω 2 �<br />

nM + K Φn = 0 . (22.51)<br />

Φn is defined with an arbitrary multiplicative constant.<br />

The natural modes are defined by the set of eigenvalues<br />

ωn and associated eigenvec<strong>to</strong>rs Φn. Spectral analysis<br />

theory shows that the Φn form an M and K orthogonal<br />

basis set, so that<br />

t ΦmMΦn = 0<br />

and<br />

t<br />

ΦmKΦn = 0 for m �= n (22.52)<br />

The orthogonality properties of the eigenmodes mean<br />

that the inertial (stiffness) forces developed in a given<br />

mode do not affect the motion of the other modes. The<br />

modes are mechanically independent. As a consequence,<br />

it is possible <strong>to</strong> expand any motion on<strong>to</strong> the eigenmodes.<br />

Given a force distribution F, the motion of the system<br />

is governed by the equation<br />

M ¨X + KX = F . (22.53)<br />

Structural <strong>Acoustics</strong> and Vibrations 22.2 Discrete Systems 909<br />

The modal projection is written in the form:<br />

X = �<br />

Φnqn(t) . (22.54)<br />

n<br />

The functions qn(t)in(22.54) are the modal participation<br />

fac<strong>to</strong>rs. Inserting (22.54) in<strong>to</strong>(22.53) and, after taking<br />

the scalar product of both sides of the equation with an<br />

eigenfunction Φm,wefind<br />

〈Φn, MΦn〉 ¨qn +〈Φn, KΦn〉qn =〈Φn, F〉 (22.55)<br />

where the notation 〈A, B〉 denotes the scalar product<br />

t A.B between vec<strong>to</strong>rs the A and B. Equation (22.55)<br />

shows that the generalized displacements are uncoupled.<br />

Each qn is governed by a single-DOF oscilla<strong>to</strong>r<br />

differential equation. The quantity<br />

mn =〈Φn, MΦn〉 (22.56)<br />

is the modal mass associated with the mode n. These<br />

coefficients are defined with an arbitrary multiplicative<br />

constant. Similarly,<br />

κn =〈Φn, KΦn〉 (22.57)<br />

is the modal stiffness, related <strong>to</strong> the modal mass through<br />

the relationship<br />

κn = mnω 2 n . (22.58)<br />

Finally, the quantity<br />

fn =〈Φn, F〉 (22.59)<br />

is the projection of the nonconservative forces on<strong>to</strong><br />

mode n. Each independent oscilla<strong>to</strong>r equation can then<br />

be rewritten as<br />

¨qn + ω 2 n qn = fn<br />

. (22.60)<br />

mn<br />

Remark. Since the eigenvec<strong>to</strong>rs are defined with arbitrary<br />

multiplicative constants Cn,(22.56) shows that the<br />

modal mass is proportional <strong>to</strong> C2 n . In addition, (22.58)<br />

shows that fn is also proportional <strong>to</strong> Cn. Therefore,<br />

through (22.60), qn is proportional <strong>to</strong> C−1 n and from<br />

(22.54), X is independent of Cn.<br />

Normal Modes. The normal modes Ψn correspond <strong>to</strong> the<br />

case where the arbitrary constant is such that the modal<br />

masses become unity. In this case, we can write<br />

Ψn = Φn<br />

√mn . (22.61)<br />

Part G 22.2

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

Saved successfully!

Ooh no, something went wrong!