Section 9.6 Modeling with Homogeneous Linear Equations
Linear homogeneous constant coefficient (LHCC) equations model systems where forces or rates depend linearly on the state and its derivatives, with no external driving. These equations appear in mechanical vibrations, electrical circuits, and many other applications where natural oscillations occur.
In this section, we’ll model a damped mass-spring system and explore the three distinct types of motion that can result: overdamping, critical damping, and underdamping.
Subsection The Problem: Suspension System Design
Automobile suspension systems use springs and shock absorbers (dampers) to provide a smooth ride. When a car hits a bump, the suspension compresses and then returns to equilibrium. The question is: how quickly should it return, and should it oscillate?
Different damping levels produce qualitatively different behaviors, each appropriate for different applications.
Subsection Assumptions
-
The spring obeys Hooke’s law: \(F_s = -kx\)
-
Damping is proportional to velocity: \(F_d = -c\frac{dx}{dt}\)
-
The mass is constant
-
No external forces act on the system
-
Motion is vertical (one-dimensional)
Subsection Building the Model
Let \(x(t)\) represent the displacement from equilibrium. Newton’s second law gives:
\begin{equation*}
m\frac{d^2x}{dt^2} = -kx - c\frac{dx}{dt}
\end{equation*}
Rearranging to standard form:
\begin{equation*}
\frac{d^2x}{dt^2} + \frac{c}{m}\frac{dx}{dt} + \frac{k}{m}x = 0
\end{equation*}
This is a second-order linear homogeneous equation with constant coefficients.
Parameters:
-
\(m\text{:}\) mass (kg)
-
\(k\text{:}\) spring constant (N/m)
-
\(c\text{:}\) damping coefficient (N·s/m or kg/s)
-
\(\omega_0 = \sqrt{k/m}\text{:}\) natural frequency (rad/s)
-
\(\zeta = \frac{c}{2\sqrt{km}}\text{:}\) damping ratio (dimensionless)
The characteristic equation is:
\begin{equation*}
r^2 + \frac{c}{m}r + \frac{k}{m} = 0
\end{equation*}
The discriminant determines the type of motion:
-
Overdamped (\(\zeta > 1\)): Two distinct real roots, no oscillation
-
Critically damped (\(\zeta = 1\)): Repeated real root, fastest return without overshoot
-
Underdamped (\(0 < \zeta < 1\)): Complex roots, damped oscillation
Subsection Physical Activity: Spring-Mass-Damper System
Objective: Build physical systems demonstrating different damping regimes.
Materials:
-
Springs (slinky or coil springs)
-
Masses (washers, weights)
-
Damping materials (cardboard, foam, or water)
-
Stopwatch
-
Ruler or motion sensor
-
Optional: video recording device
Setup Variations:
-
Underdamped: Mass-spring in air (minimal damping)
-
Increased Damping: Attach cardboard fins or foam to increase air resistance
-
Heavy Damping: Mass-spring system moving through water or oil
Procedure:
-
For each setup, pull the mass down a fixed distance and release.
-
Record the motion (video or position measurements over time).
-
Count oscillations (if any) and measure the period.
-
Measure the decay time (time to reduce amplitude by half).
-
Note whether the system overshoots equilibrium.
Analysis:
-
Plot displacement versus time for each damping level.
-
Classify each system as underdamped, critically damped, or overdamped.
-
For underdamped cases, measure the damped frequency and compare with the natural frequency.
-
Estimate the damping ratio \(\zeta\) from the logarithmic decrement method (if applicable).
Subsection Computational Activity
Use the exact solutions to generate plots:
Case 1: Underdamped (\(\zeta = 0.2\))
\begin{equation*}
x(t) = e^{-\zeta\omega_0 t}(A\cos(\omega_d t) + B\sin(\omega_d t))
\end{equation*}
where \(\omega_d = \omega_0\sqrt{1-\zeta^2}\)
Case 2: Critically Damped (\(\zeta = 1\))
\begin{equation*}
x(t) = e^{-\omega_0 t}(A + Bt)
\end{equation*}
Case 3: Overdamped (\(\zeta = 2\))
\begin{equation*}
x(t) = Ae^{r_1 t} + Be^{r_2 t}
\end{equation*}
where \(r_1, r_2\) are negative real roots.
-
Solve for constants \(A\) and \(B\)
-
Plot \(x(t)\) for \(0 \leq t \leq 20\)
-
Compare the three cases on one graph
-
Which returns to equilibrium fastest? Which overshoots?
Subsection Analytical Questions
-
Verify that \(x(t) = e^{-\zeta\omega_0 t}\cos(\omega_d t)\) satisfies the differential equation for underdamped motion.
-
Show that critical damping occurs when \(c = 2\sqrt{km}\text{.}\)
-
For underdamped motion, derive the relationship between consecutive peak amplitudes (logarithmic decrement).
-
Why is critical damping preferred for car suspensions and door closers?
-
In what applications would you want underdamping? Overdamping?
-
How does doubling the mass affect the natural frequency? The damping ratio?
Subsection Design Challenge
You are designing a suspension system for a 1000 kg vehicle. The spring constant is \(k = 40000\) N/m.
-
Calculate the natural frequency \(\omega_0\text{.}\)
-
For a comfortable ride, you want \(\zeta = 0.7\) (underdamped but not too oscillatory). What damping coefficient \(c\) is needed?
-
If the car hits a bump causing an initial displacement of 0.1 m, how long until the displacement is less than 0.01 m?
-
Compare this design with a critically damped system. Which returns to equilibrium faster?
Subsection Final Report
Prepare a detailed report (4-5 pages) including:
-
Derivation of the damped harmonic oscillator equation from Newton’s laws.
-
Solution of the characteristic equation for all three damping cases.
-
Physical experiment description with photos/diagrams and data.
-
Classification of your experimental systems by damping type.
-
Graphs comparing the three damping regimes (both experimental and theoretical).
-
Parameter estimation for your physical systems.
-
Complete solution to the design challenge with justification.
-
Discussion of real-world applications where each damping type is preferred.
-
Analysis of how parameter changes affect system behavior.
Linear homogeneous constant coefficient equations provide exact, predictable solutions for a wide class of physical systems. The damped harmonic oscillator is fundamental in engineering design, from vehicle suspensions to building earthquake resistance, from electrical filters to quantum mechanics. Understanding how the characteristic equation’s roots determine system behavior—whether oscillatory, exponential, or critical—is essential for controlling and optimizing these systems. The mathematical classification (overdamped, critically damped, underdamped) directly translates to observable physical phenomena, demonstrating the power of differential equations in engineering and science.
You have attempted of activities on this page.
