Engineering Optimization: Step-by-Step Derivations of Vibration Systems

Abstract

Vibration analysis forms the foundation of mechanical system design and optimization. This article derives the key relationships governing oscillatory motion in spring-mass systems, develops the concept of equivalent spring constants for complex structural elements, and demonstrates how energy methods simplify optimization problems. Through worked examples, we show how engineers can reduce complex geometries to tractable models and predict system behavior under dynamic loading.

Background

^[Vibration] is the oscillatory motion of a system around an equilibrium position. In engineering practice, controlling vibration is essential: excessive vibrations cause fatigue failure in structures, wear in machinery, and performance degradation in precision instruments. To optimize designs, engineers must predict how systems respond to dynamic loads.

The ^{[spring-mass model]} provides the simplest framework for understanding vibrations. A mass attached to a spring obeys Hooke's Law:

$F = -kx$

where $F$ is the restoring force, $k$ is the spring constant, and $x$ is displacement from equilibrium ^{[[spring-mass-model]}]. This linear relationship enables analytical solutions and serves as a foundation for more complex systems.

^[Stiffness] quantifies a material's resistance to deformation. Formally, stiffness is the ratio of applied force to resulting displacement:

$k = \frac{F}{x}$

^[[stiffness]]. In vibration problems, stiffness directly determines how quickly a system oscillates—stiffer systems vibrate faster. For engineers, understanding stiffness is critical because it governs natural frequencies, which in turn determine whether external excitations will cause resonance.

Key Results

Mechanical Energy and the Spring-Mass System

^{[Mechanical energy]} in a vibrating system comprises potential and kinetic components:

$E_{\text{mechanical}} = E_{\text{potential}} + E_{\text{kinetic}}$

For a spring-mass system, these are:

$PE = \frac{1}{2}kx^2, \quad KE = \frac{1}{2}mv^2$

^{[[mechanical-energy]}]. In an ideal, undamped system, total mechanical energy is conserved. This energy exchange is the mechanism of oscillation: when the spring is maximally compressed, all energy is potential and velocity is zero; at equilibrium, all energy is kinetic and velocity is maximum.

Energy conservation provides an elegant path to deriving the equation of motion. At any instant:

$E_{\text{total}} = \frac{1}{2}kx^2 + \frac{1}{2}mv^2 = \text{constant}$

Differentiating with respect to time:

$kx\frac{dx}{dt} + mv\frac{dv}{dt} = 0$

Since $v = \frac{dx}{dt}$ and $a = \frac{dv}{dt}$:

$kx \cdot v + m \cdot v \cdot a = 0$

Dividing by $v$ (when $v \neq 0$):

$kx + ma = 0$

$m\frac{d^2x}{dt^2} + kx = 0$

This is the fundamental equation of motion for an undamped spring-mass system. Its solution is harmonic oscillation with natural frequency:

$\omega_n = \sqrt{\frac{k}{m}}$

Equivalent Spring Constants for Structural Elements

Real engineering structures are not simple springs. Beams, rods, and shafts deform under load, and their stiffness can be modeled as equivalent springs. ^{[Equivalent massless spring constants]} allow engineers to replace complex geometries with simple spring models.

The equivalent spring constant depends on the structural configuration and boundary conditions ^{[[equivalent-massless-spring-constants]}]:

• Torsional spring: $k_t = \frac{EI}{L}$
• Rod in axial deformation: $k_a = \frac{EA}{L}$
• Shaft in torsion: $k_s = \frac{GJ}{L}$
• Helical spring: $k_h = \frac{Gd^4}{64nR^3}$
• Cantilever beam (tip load): $k_c = \frac{3EI}{L^3}$
• Pinned-pinned beam (midspan load): $k_{pp} = \frac{48EI}{L^3}$
• Clamped-clamped beam (midspan load): $k_{cc} = \frac{192EI}{L^3}$

Here, $E$ is Young's modulus, $G$ is shear modulus, $I$ is the second moment of inertia, $J$ is the polar moment, $A$ is cross-sectional area, and $L$ is length.

The physical interpretation is straightforward: stiffness is proportional to material stiffness ($E$ or $G$) and geometric resistance to deformation ($I$, $A$, or $J$), and inversely proportional to length. Longer elements are more compliant; thicker elements are stiffer.

Combining Spring Constants

When multiple springs act together, their combined stiffness depends on their arrangement. For springs in series (one after another), the compliances add:

$\frac{1}{k_{eq}} = \frac{1}{k_1} + \frac{1}{k_2} + \cdots$

For springs in parallel (side by side), the stiffnesses add directly:

$k_{eq} = k_1 + k_2 + \cdots$

These rules follow from force and displacement compatibility. In series, the same force passes through each spring, but displacements add. In parallel, the same displacement occurs across each spring, but forces add.

Worked Examples

Example 1: Natural Frequency of a Cantilever Beam with a Mass

A cantilever beam of length $L = 1$ m, with Young's modulus $E = 200$ GPa and second moment of inertia $I = 1 \times 10^{-8}$ m$^4$, supports a point mass $m = 10$ kg at its tip.

Step 1: Calculate the equivalent spring constant using the cantilever formula ^{[[equivalent-massless-spring-constants]}]:

$k_c = \frac{3EI}{L^3} = \frac{3 \times 200 \times 10^9 \times 1 \times 10^{-8}}{1^3} = 6 \times 10^4 \text{ N/m}$

Step 2: Calculate the natural frequency:

$\omega_n = \sqrt{\frac{k}{m}} = \sqrt{\frac{6 \times 10^4}{10}} = \sqrt{6000} \approx 77.5 \text{ rad/s}$

In Hz: $f_n = \frac{\omega_n}{2\pi} \approx 12.3$ Hz

Step 3: Interpret the result. The beam oscillates at approximately 12.3 cycles per second. If external vibrations near this frequency are present, resonance will occur, amplifying motion and potentially causing failure.

Example 2: Maximum Displacement from Energy Conservation

A spring-mass system with $k = 1000$ N/m and $m = 5$ kg is displaced $x_0 = 0.1$ m from equilibrium and released from rest.

Step 1: Calculate initial mechanical energy:

$E_{\text{total}} = \frac{1}{2}kx_0^2 = \frac{1}{2} \times 1000 \times (0.1)^2 = 5 \text{ J}$

Step 2: At maximum displacement (where velocity is zero), all energy is potential:

$E_{\text{total}} = \frac{1}{2}kx_{\max}^2$

$5 = \frac{1}{2} \times 1000 \times x_{\max}^2$

$x_{\max} = 0.1 \text{ m}$

Step 3: At equilibrium (where displacement is zero), all energy is kinetic:

$E_{\text{total}} = \frac{1}{2}mv_{\max}^2$

$5 = \frac{1}{2} \times 5 \times v_{\max}^2$

$v_{\max} = \sqrt{2} \approx 1.41 \text{ m/s}$

The system oscillates between $\pm 0.1$ m displacement with a maximum velocity of 1.41 m/s at the equilibrium position.

References

• ^[vibration]
• ^{[spring-mass-model]}
• ^{[mechanical-energy]}
• ^[stiffness]
• ^{[equivalent-massless-spring-constants]}
• ^{[mechanical-energy]}

AI Disclosure

This article was drafted with AI assistance from class notes organized in a Zettelkasten system. All mathematical derivations and technical claims are grounded in cited notes from an engineering optimization course. The worked examples were generated to illustrate the concepts but follow standard textbook approaches. The author retains responsibility for accuracy and interpretation.