Shaoyun Wang
← All notes
Mechanics · 2024-10-14

Vibration of Continuous System

PDF

1 Questions

How to make equivalence from real structures?

2 Vibration mode

2.1 Longitudinal vibration

Area: AA, Length: LL.

Governing equation:

ρutt=Euxx\rho u_{tt} = E u_{xx}

Frequency of lowest mode under fixed-fixed, free-free condition:

ω=EρπL\omega = \sqrt{\frac{E}{\rho}} \frac{\pi}{L}

Frequency of lowest mode under fixed-free condition:

ω=Eρπ2L\omega = \sqrt{\frac{E}{\rho}} \frac{\pi}{2L}

Rao, S. S. (2019). Vibration of continuous systems. John Wiley & Sons.

2.2 Flexural vibration

The governing equation is

∂2∂x2(EI∂2w∂x2)+ρA∂2w∂t2=f(x,t)\frac{\partial^{2}}{\partial x^{2}}\left(E I \frac{\partial^{2} w}{\partial x^{2}}\right)+\rho A \frac{\partial^{2} w}{\partial t^{2}}=f(x, t)

where

I=Iy=∬Az2dAI=I_{y}=\iint_{A} z^{2} d A

Lowest mode for beam simply supported at both ends

ω=(β1l)2(EIρAl4)1/2=π2(EIρAL4)1/2.\omega=\left(\beta_{1} l\right)^{2}\left(\frac{E I}{\rho A l^{4}}\right)^{1 / 2}=\pi^{2}\left(\frac{E I}{\rho A L^{4}}\right)^{1 / 2}.

2.3 Torsional Vibration of Shafts

The governing equation is

∂∂x(GJ∂θ(x,t)∂x)+mt(x,t)=I0∂2θ(x,t)∂t2\frac{\partial}{\partial x}\left(G J \frac{\partial \theta(x, t)}{\partial x}\right)+m_{t}(x, t)=I_{0} \frac{\partial^{2} \theta(x, t)}{\partial t^{2}}

where I0=ρJI_0 = \rho J, J=∫Ar2dAJ = \int_A r^2 dA. Uniform and without attached mass mtm_t, we have

G∂2θ(x,t)∂x2=ρ∂2θ(x,t)∂t2G \frac{\partial^2 \theta(x, t)}{\partial x^2}=\rho \frac{\partial^{2} \theta(x, t)}{\partial t^{2}}

Frequency of lowest mode under fixed-fixed, free-free condition:

ω=GρπL\omega = \sqrt{\frac{G}{\rho}} \frac{\pi}{L}

Frequency of lowest mode under fixed-free condition:

ω=Gρπ2L\omega = \sqrt{\frac{G}{\rho}} \frac{\pi}{2L}

where

2.4 Shear vibration

Lowest frequency for shear mode is

ω=π2bμρ\omega=\frac{\pi}{2 b} \sqrt{\frac{\mu}{\rho}}

where bb is half of the thickness dd. Or

ω=π4dμρ\omega=\frac{\pi}{4 d} \sqrt{\frac{\mu}{\rho}}

where dd is the thickness.

2.5 Bending vibration of the plate

Governing equation is

D∇4w+ρhw¨−f=0D \nabla^{4} w+\rho h \ddot{w}-f=0

in which DD represents the flexural rigidity of the plate:

D=Eh312(1−ν2)D=\frac{E h^{3}}{12\left(1-\nu^{2}\right)}

Solution for a simply supported rectangular plate is

ωmn=λmn2(Dρh)1/2=π2[(ma)2+(nb)2](Dρh)1/2\omega_{m n}=\lambda_{m n}^{2}\left(\frac{D}{\rho h}\right)^{1 / 2}=\pi^{2}\left[\left(\frac{m}{a}\right)^{2}+\left(\frac{n}{b}\right)^{2}\right]\left(\frac{D}{\rho h}\right)^{1 / 2}

The lowest mode (if a>ba>b) is

ω=π2(1a)2(Dρh)1/2\omega=\pi^{2}\left(\frac{1}{a}\right)^{2}\left(\frac{D}{\rho h}\right)^{1 / 2}

Solution for circular plate is

ωmn=λmn2(Dρh)1/2\omega_{m n}=\lambda_{m n}^{2}\left(\frac{D}{\rho h}\right)^{1 / 2}

where λmn\lambda_{mn} is the solution of the following equation

Im(λa)Jm−1(λa)−Jm(λa)Im−1(λa)=0,m=0,1,2,…I_{m}(\lambda a) J_{m-1}(\lambda a)-J_{m}(\lambda a) I_{m-1}(\lambda a)=0, \quad m=0,1,2, \ldots

The lowest root is λ01a=3.196\lambda_{01} a=3.196. Therefore, the frequency of lowest mode is

ω=(3.196a)2(Dρh)1/2\omega=\left(\frac{3.196}{a}\right)^{2}\left(\frac{D}{\rho h}\right)^{1 / 2}