Applications of Higher-Order DEs

Higher-order differential equations, especially second-order linear DEs, are fundamental in modeling oscillatory systems, electrical circuits, and structural mechanics.

Spring-Mass Systems

A mass mm attached to a spring with spring constant kk and damping coefficient cc is modeled by Newton's Second Law.
md2xdt2+cdxdt+kx=F(t)m \frac{d^2x}{dt^2} + c \frac{dx}{dt} + kx = F(t)
where F(t)F(t) is an external driving force. Dividing by mm gives the standard form:
d2xdt2+2ζωndxdt+ωn2x=f(t)\frac{d^2x}{dt^2} + 2\zeta\omega_n \frac{dx}{dt} + \omega_n^2x = f(t)
where ωn=k/m\omega_n = \sqrt{k/m} is the natural frequency, and ζ=c2km\zeta = \frac{c}{2\sqrt{km}} is the damping ratio.
  • Free Undamped Motion (c=0,F(t)=0c=0, F(t)=0): Simple Harmonic Motion (x+ωn2x=0x'' + \omega_n^2 x = 0). The solution is x(t)=C1cos(ωnt)+C2sin(ωnt)x(t) = C_1\cos(\omega_nt) + C_2\sin(\omega_nt).
  • Free Damped Motion (c>0,F(t)=0c > 0, F(t)=0): The auxiliary equation is mr2+cr+k=0m r^2 + cr + k = 0. The roots determine the type of damping based on the discriminant D=c24mkD = c^2 - 4mk.
    • Overdamped (D>0D > 0): Two real distinct roots. System slowly returns to equilibrium without oscillating. x(t)=c1er1t+c2er2tx(t) = c_1e^{r_1t} + c_2e^{r_2t}.
    • Critically Damped (D=0D = 0): One repeated real root. System returns to equilibrium as fast as possible without oscillating. x(t)=c1ert+c2tertx(t) = c_1e^{rt} + c_2te^{rt}.
    • Underdamped (D<0D < 0): Complex conjugate roots (α±βi\alpha \pm \beta i). System oscillates with exponentially decreasing amplitude. x(t)=eαt(c1cos(βt)+c2sin(βt))x(t) = e^{\alpha t}(c_1\cos(\beta t) + c_2\sin(\beta t)).
  • Forced Motion (F(t)0F(t) \neq 0): Add the particular solution xp(t)x_p(t). Resonance occurs when the driving frequency matches the natural frequency (ω=ωn\omega = \omega_n) in an undamped system, causing amplitude to grow infinitely over time (xp(t)tsin(ωnt)x_p(t) \propto t\sin(\omega_nt)).

Spring-Mass Simulation

Explore how mass, damping, and spring stiffness affect the motion of the system.

Spring-Mass-Damper System

Parameters

1 kg
0.5 Ns/m
10 N/m
1 m
System State
Underdamped
Crit. Damping = 6.32
m
Eq
Loading chart...

The Simple Pendulum

A simple pendulum of mass mm attached to a string of length LL swinging under gravity gg is governed by:
d2θdt2+gLsin(θ)=0\frac{d^2\theta}{dt^2} + \frac{g}{L}\sin(\theta) = 0
This is a non-linear differential equation. However, for small angles (θ0\theta \approx 0), we can use the Taylor series approximation sin(θ)θ\sin(\theta) \approx \theta.
The linearized equation becomes a simple harmonic oscillator:
d2θdt2+gLθ=0\frac{d^2\theta}{dt^2} + \frac{g}{L}\theta = 0
The natural frequency is ωn=g/L\omega_n = \sqrt{g/L}, and the period is T=2πL/gT = 2\pi\sqrt{L/g}.

RLC Circuits

An electrical circuit with a resistor (RR), inductor (LL), and capacitor (CC) connected in series is modeled by Kirchhoff's Voltage Law:
Ld2qdt2+Rdqdt+1Cq=E(t)L \frac{d^2q}{dt^2} + R \frac{dq}{dt} + \frac{1}{C}q = E(t)
where q(t)q(t) is the charge on the capacitor and E(t)E(t) is the electromotive force (voltage). Differentiating with respect to tt gives the equation for current i(t)i(t):
Ld2idt2+Rdidt+1Ci=dEdtL \frac{d^2i}{dt^2} + R \frac{di}{dt} + \frac{1}{C}i = \frac{dE}{dt}Analogy between Mechanical and Electrical Systems:
  • Mass mm \leftrightarrow Inductance LL
  • Damping cc \leftrightarrow Resistance RR
  • Spring Constant kk \leftrightarrow Inverse Capacitance 1/C1/C
  • Displacement x(t)x(t) \leftrightarrow Charge q(t)q(t)
  • Velocity v(t)v(t) \leftrightarrow Current i(t)i(t)

Simple Beam Deflection

The vertical deflection y(x)y(x) of an elastic beam under a distributed transverse load w(x)w(x) is governed by the fourth-order linear DE:
EId4ydx4=w(x)EI \frac{d^4y}{dx^4} = w(x)
where EE is Young's Modulus, II is the area moment of inertia, and EIEI is the flexural rigidity. This is solved by integrating four times and applying boundary conditions at the supports (e.g., fixed ends have y=0,y=0y=0, y'=0; pinned ends have y=0,y=0y=0, y''=0).
Key Takeaways
  • Spring-Mass Systems: Governed by mx+cx+kx=F(t)mx'' + cx' + kx = F(t). Damping ratio determines if motion is oscillatory.
  • Damping Regimes: Overdamped (c2>4mkc^2 > 4mk), Critically Damped (c2=4mkc^2 = 4mk), Underdamped (c2<4mkc^2 < 4mk).
  • Simple Pendulum: Modeled by θ+(g/L)θ=0\theta'' + (g/L)\theta = 0 for small angles.
  • RLC Circuits: Perfectly analogous to damped mechanical systems (Lm,Rc,1/CkL \leftrightarrow m, R \leftrightarrow c, 1/C \leftrightarrow k). The homogeneous solutions mirror underdamped, overdamped, and critically damped behavior.
  • Beam Deflection: Requires solving a 4th-order DE EIy(4)=w(x)EI y^{(4)} = w(x) with appropriate boundary conditions based on support types (fixed, pinned, free).