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.

from which we can derive through (22.174) the eigenfrequencies<br />

ωn ≈ nπc � n<br />

1 + ε 2π2 �<br />

. (22.176)<br />

L<br />

2<br />

The eigenfrequencies are raised by the stiffness. The<br />

difference increases as n2 . The inharmonicity coefficient<br />

can be defined as<br />

i = ε n2π 2<br />

. (22.177)<br />

2<br />

22.4.3 Flexural Vibrations of Thin Plates<br />

The thin plate, or Kirchhoff–Love plate, described below<br />

is a generalization of the Euler–Bernoulli beam<br />

in 2-D [22.13]. The general case of orthotropic plates<br />

is treated here. The problem is solved in Cartesian<br />

coordinates, where w(x, y, t) denotes the transverse displacement.<br />

It is assumed that the coordinates coincide<br />

with the symmetry axes of the material. The assumption<br />

of orthotropy leads <strong>to</strong> the following relations between<br />

plane stresses σij and plane strains εkl<br />

⎛<br />

⎜<br />

⎝<br />

σxx<br />

σyy<br />

σxy<br />

⎛<br />

Ex νyx Ex<br />

⎞ ⎜<br />

⎜1<br />

− νxyνyx 1 − νxyνyx<br />

⎜<br />

⎟ ⎜<br />

⎠ = ⎜ νyx Ex Ey<br />

⎜<br />

⎜1<br />

− νxyνyx 1 − νxyνyx<br />

⎝<br />

0<br />

0<br />

⎞<br />

⎟<br />

⎠<br />

⎛ ⎞<br />

εxx<br />

⎜ ⎟<br />

× ⎝ ⎠ .<br />

0 0 2Gxy<br />

(22.178)<br />

εyy<br />

εxy<br />

These relations involve two elastic moduli Ex and<br />

Ey, two Poisson’s ratios νxy and one <strong>to</strong>rsional<br />

modulus Gxy. In addition, we have the property<br />

Exνyx = Eyνxy [22.14]. The bending and twisting moments<br />

are obtained by integration of the elementary<br />

moments over the thickness h of the plate:<br />

Mx =<br />

My =<br />

�h/2<br />

−h/2<br />

�h/2<br />

−h/2<br />

Mxy = Myx =<br />

zσxx dz ;<br />

zσyy dz ;<br />

�h/2<br />

−h/2<br />

zσxydz (22.179)<br />

Structural <strong>Acoustics</strong> and Vibrations 22.4 Bars, Plates and Shells 923<br />

which leads <strong>to</strong> the matrix relation between the moments<br />

and curvatures<br />

⎛<br />

∂<br />

⎛ ⎞ ⎛<br />

Mx D1<br />

⎜ ⎟ ⎜<br />

⎝ My ⎠ =−⎝D2/2<br />

Mxy 0<br />

D2/2<br />

D3<br />

0<br />

⎞ ⎜<br />

0 ⎜<br />

⎟ ⎜<br />

0 ⎠ × ⎜<br />

D4/2 ⎜<br />

⎝<br />

2w ∂x2 ∂2w ∂y2 ∂2 ⎞<br />

⎟ ,<br />

⎟<br />

w ⎠<br />

∂x∂y<br />

(22.180)<br />

where<br />

D1 =<br />

D2 =<br />

D3 =<br />

Exh 3<br />

12(1 − νxyνyx) ;<br />

Exνyxh3 Eyνxyh3 = , (22.181)<br />

6(1 − νxyνyx) 6(1 − νxyνyx)<br />

Eyh 3<br />

12(1 − νxyνyx) ;<br />

Gxyh3 D4 = . (22.182)<br />

3<br />

Again, the equation of motion can be obtained using<br />

Hamil<strong>to</strong>n’s principle [22.12]. Finally, given an external<br />

force density term f (x, y, t), we obtain the equation of<br />

motion for the plate:<br />

ρph ∂2w ∂t2 = ∂2Mx ∂x2 + ∂2My ∂y2 + 2 ∂2Mxy + f (x, y, t) .<br />

∂x∂y<br />

(22.183)<br />

Boundary Conditions. As for the beam, the boundary<br />

conditions follow from integration by parts of the elastic<br />

energy. The number of possible conditions depends on<br />

the selected geometry. For rectangular plates, for example,<br />

Leissa lists 21 possible boundary conditions [22.15].<br />

Along the principal directions, for example, the conditions<br />

of greatest practical interest are:<br />

1. Clamped edge: displacement w = 0 and rotation<br />

∂w<br />

∂x = 0;<br />

2. Simply supported edge: displacement w = 0and<br />

bending moment Mx = 0;<br />

3. Free edge: bending moment Mx = 0 and shear force<br />

∂Mx ∂Mxy<br />

∂x + 2 ∂y = 0.<br />

Part G 22.4

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

Saved successfully!

Ooh no, something went wrong!