Skip to main content

Lateral vibration of a simply supported beam

Short answer

The Euler-Bernoulli equation is EI y'''' + rho A y'' = 0; for simply supported ends, sin(beta l)=0 gives omega_n = n^2 pi^2/l^2 sqrt(EI/(rho A)). With I=4.1667e-6 m^4 and A=0.005 m^2, the first three frequencies are 25.43 Hz, 101.72 Hz, and 228.88 Hz.

(a) Derive the equation for the lateral (flexural) vibration of a uniform beam and obtain the natural frequencies of a simply supported beam. (b) A simply supported steel beam of length 3 m has a rectangular section 50 mm wide and 100 mm deep (vibrating in the 100 mm deep direction). Find the first three natural frequencies. E=200E = 200 GPa, ρ=7850\rho = 7850 kg/m³.

(a) Euler-Bernoulli beam equation

Consider an element dxdx of a beam with bending stiffness EIEI, mass per unit length ρA\rho A, deflection y(x,t)y(x,t), shear force VV and bending moment MM.

     V+dV
 M ->|====|<- M+dM
     V

Vertical equilibrium:   V−(V+dV)=ρA dx ∂2y∂t2\;V - (V + dV) = \rho A\,dx\,\dfrac{\partial^2 y}{\partial t^2}, that is −∂V∂x=ρA∂2y∂t2-\dfrac{\partial V}{\partial x} = \rho A\dfrac{\partial^2 y}{\partial t^2}.

Moments (neglecting rotary inertia and shear deformation): V=∂M∂xV = \dfrac{\partial M}{\partial x} and M=EI∂2y∂x2M = EI\dfrac{\partial^2 y}{\partial x^2}. Hence

EI∂4y∂x4+ρA∂2y∂t2=0EI\frac{\partial^4 y}{\partial x^4} + \rho A\frac{\partial^2 y}{\partial t^2} = 0

Put y=Y(x)sin⁡ωty = Y(x)\sin\omega t and β4=ρAω2EI\beta^4 = \dfrac{\rho A\omega^2}{EI}:

d4Ydx4−β4Y=0  ⇒  Y=C1cos⁡βx+C2sin⁡βx+C3cosh⁡βx+C4sinh⁡βx\frac{d^4Y}{dx^4} - \beta^4Y = 0 \;\Rightarrow\; Y = C_1\cos\beta x + C_2\sin\beta x + C_3\cosh\beta x + C_4\sinh\beta x

Simply supported ends: Y(0)=Y′′(0)=0⇒C1=C3=0Y(0) = Y''(0) = 0 \Rightarrow C_1 = C_3 = 0. At x=lx = l, Y(l)=Y′′(l)=0Y(l) = Y''(l) = 0 gives C4=0C_4 = 0 and

sin⁡βl=0  ⇒  βnl=nπ\sin\beta l = 0 \;\Rightarrow\; \beta_nl = n\pi ωn=(nπ)2EIρAl4=n2π2l2EIρA,Yn=sin⁡nπxl\omega_n = (n\pi)^2\sqrt{\frac{EI}{\rho Al^4}} = \frac{n^2\pi^2}{l^2}\sqrt{\frac{EI}{\rho A}}, \qquad Y_n = \sin\frac{n\pi x}{l}

(b) Numerical

I=bh312=0.05×0.1312=4.1667×10−6 m4,A=bh=0.005 m2I = \frac{bh^3}{12} = \frac{0.05 \times 0.1^3}{12} = 4.1667\times10^{-6}\ \text{m}^4, \qquad A = bh = 0.005\ \text{m}^2 EIρA=200×109×4.1667×10−67850×0.005=145.710 m2/s\sqrt{\frac{EI}{\rho A}} = \sqrt{\frac{200\times10^9 \times 4.1667\times10^{-6}}{7850 \times 0.005}} = 145.710\ \text{m}^2/\text{s} ω1=π232×145.710=159.79 rad/s,f1=25.43 Hz\omega_1 = \frac{\pi^2}{3^2}\times 145.710 = 159.79\ \text{rad/s}, \qquad f_1 = 25.43\ \text{Hz} f2=4f1=101.72 Hz,f3=9f1=228.88 Hzf_2 = 4f_1 = 101.72\ \text{Hz}, \qquad f_3 = 9f_1 = 228.88\ \text{Hz}

Answer: f1=25.43f_1 = 25.43 Hz, f2=101.72f_2 = 101.72 Hz, f3=228.88f_3 = 228.88 Hz.

✓ Answer checked by Claude Opus 5.5 (high effort)

All Theory of Mechanism and Machine II chapter 10 questions →