Engineering Optimization: Foundational Theorems in Mechanical Vibration

Abstract

Mechanical vibration analysis underpins the design and optimization of engineered systems across civil, mechanical, and aerospace domains. This article synthesizes core theoretical results governing vibrating systems: the spring-mass model, mechanical energy conservation, and equivalent stiffness formulations. We establish the mathematical foundations necessary for predicting system response and demonstrate how these principles enable systematic design optimization to minimize unwanted vibrations and ensure structural safety.

Background

^[Vibration] is the oscillatory motion of a system relative to equilibrium, characterized by displacement, velocity, and acceleration varying over time. In engineering practice, vibrations present both challenges and design constraints: excessive vibrations cause fatigue failure in structures and accelerated wear in machinery, while controlled vibrations are exploited in applications ranging from seismic isolation to precision manufacturing.

The fundamental challenge in optimization is predicting how a given system will vibrate under specified conditions, then modifying design parameters to achieve desired dynamic behavior. This requires a rigorous mathematical framework connecting physical properties (mass, stiffness, damping) to system response.

The ^{[spring-mass model]} provides the essential starting point. This idealized system—a point mass attached to a linear elastic element—captures the essential physics of countless real systems. By Hooke's Law, the restoring force is proportional to displacement:

$F = -kx$

where $k$ is the spring constant and $x$ is displacement from equilibrium ^{[[spring-mass-model]}]. This linear relationship is the foundation upon which all subsequent analysis rests.

Key Results

Energy Conservation in Vibrating Systems

A central theorem in vibration analysis concerns the exchange of energy forms. ^{[Mechanical energy]} in a vibrating system is conserved (in the absence of damping) and comprises two components:

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

For a spring-mass system, these are expressed as:

$PE = \frac{1}{2}kx^2$

$KE = \frac{1}{2}mv^2$

where $m$ is mass and $v$ is velocity ^{[[mechanical-energy]}].

Theorem 1 (Energy Conservation): In an undamped spring-mass system, the total mechanical energy remains constant throughout the oscillation cycle. At maximum displacement, all energy is potential; at the equilibrium position, all energy is kinetic.

This result is not merely descriptive—it provides a powerful tool for optimization. By controlling the stiffness $k$ and mass $m$, engineers directly control the energy storage capacity and oscillation frequency of the system. Systems designed with appropriate stiffness-to-mass ratios can be tuned to avoid resonance with external forcing frequencies.

Stiffness and System Response

^[Stiffness] is defined as the ratio of applied force to resulting displacement:

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

^[[stiffness]] is not merely a material property but a system property that depends on geometry, boundary conditions, and material composition. This distinction is critical for optimization: engineers can modify stiffness through geometric design without changing material.

Theorem 2 (Stiffness Determines Natural Frequency): The natural frequency of a spring-mass system is proportional to the square root of the stiffness-to-mass ratio. Increasing stiffness raises the natural frequency; increasing mass lowers it. This relationship is the basis for frequency-tuning in design.

Equivalent Spring Constants for Structural Elements

Real engineering systems are rarely simple point-mass-spring configurations. Beams, shafts, and composite structures must be reduced to equivalent spring constants for analysis. ^{[Equivalent massless spring constants]} provide formulas for common structural geometries:

For a cantilever beam with a concentrated load at the free end:

$k_c = \frac{3EI}{L^3}$

For a pinned-pinned beam with a load at midspan:

$k_{pp} = \frac{48EI}{L^3}$

For a clamped-clamped beam with a load at midspan:

$k_{cc} = \frac{192EI}{L^3}$

For axial deformation of a rod:

$k_a = \frac{EA}{L}$

For torsional deformation of a shaft:

$k_s = \frac{GJ}{L}$

where $E$ is Young's modulus, $G$ is the shear modulus, $I$ is the second moment of area, $J$ is the polar moment of inertia, $A$ is cross-sectional area, and $L$ is length ^{[[equivalent-massless-spring-constants]}].

Theorem 3 (Geometric Stiffness Scaling): For beam-like structures, stiffness scales inversely with the cube of length and directly with the second moment of area. This cubic dependence on length is the basis for the design principle that longer, more slender structures are significantly more flexible—a fact exploited in seismic design and vibration isolation.

Composite Systems and Series/Parallel Combinations

When multiple elastic elements are combined, the overall system stiffness depends on their arrangement. For springs in series, the reciprocal stiffnesses add; for springs in parallel, the stiffnesses add directly. This principle extends to complex structures through the ^{[equivalent spring constant]} framework, which allows reduction of multi-element systems to a single effective stiffness ^{[[equivalent-spring-constant]}].

Worked Examples

Example 1: Cantilever Beam Optimization

Consider a cantilever beam of length $L = 2$ m, with Young's modulus $E = 200$ GPa and second moment of area $I = 1 \times 10^{-5}$ m$^4$. A mass $m = 100$ kg is attached at the free end.

Using Theorem 3, the equivalent stiffness is:

$k_c = \frac{3 \times 200 \times 10^9 \times 1 \times 10^{-5}}{2^3} = \frac{6 \times 10^6}{8} = 7.5 \times 10^5 \text{ N/m}$

The natural frequency is:

$f_n = \frac{1}{2\pi}\sqrt{\frac{k}{m}} = \frac{1}{2\pi}\sqrt{\frac{7.5 \times 10^5}{100}} = \frac{1}{2\pi}\sqrt{7500} \approx 8.7 \text{ Hz}$

If external forcing occurs at 10 Hz, this system is near resonance and will experience large amplitudes. To increase the natural frequency to 12 Hz, we need:

$k_{\text{new}} = (2\pi \times 12)^2 \times 100 \approx 5.68 \times 10^6 \text{ N/m}$

This can be achieved by increasing $I$ (using a stiffer cross-section) or reducing $L$ (shortening the beam). The trade-off between these options depends on functional constraints.

Example 2: Energy-Based Design Check

For the same cantilever system, verify energy conservation. At maximum displacement $x_{\max} = 0.01$ m:

$PE_{\max} = \frac{1}{2} \times 7.5 \times 10^5 \times (0.01)^2 = 37.5 \text{ J}$

At the equilibrium position, all energy is kinetic:

$KE_{\max} = \frac{1}{2} \times 100 \times v_{\max}^2 = 37.5 \text{ J}$

$v_{\max} = \sqrt{\frac{2 \times 37.5}{100}} = 0.866 \text{ m/s}$

This consistency confirms the system model and provides a basis for predicting response to initial conditions or external forcing.

References

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

AI Disclosure

This article was drafted with AI assistance from class notes (Zettelkasten). All mathematical claims and theorems are grounded in cited notes. The worked examples and some explanatory passages were generated by AI to clarify the source material. The author reviewed all content for technical accuracy and relevance to the engineering optimization course context.