Applications of Higher-Order DEs

Learning Objectives

  • Understand the mathematical models for Spring-Mass Systems.
  • Classify the damping regimes (overdamped, critically damped, underdamped).
  • Analyze the small-angle approximation for simple pendulums.
  • Relate electrical RLC circuits to mechanical damped systems.
  • Formulate and solve the simple beam deflection differential equation.

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. An external driving force F(t)F(t) can be applied to the mass.

Dividing the general equation by mm gives the standard form, which uses the natural frequency Ο‰n\omega_n and the damping ratio ΞΆ\zeta.

  • 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=c2βˆ’4mkD = 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)).

Newton's Second Law for Spring-Mass Systems

The general differential equation for a driven, damped harmonic oscillator.

md2xdt2+cdxdt+kx=F(t)m \frac{d^2x}{dt^2} + c \frac{dx}{dt} + kx = F(t)

Variables

SymbolDescriptionUnit
mmMasskg
ccDamping coefficientNΒ·s/m
kkSpring constantN/m
xxDisplacement from equilibriumm
F(t)F(t)External driving forceN
ttTimes

Standard Form of Spring-Mass Equation

The normalized differential equation used to analyze system behavior.

d2xdt2+2ΞΆΟ‰ndxdt+Ο‰n2x=f(t)\frac{d^2x}{dt^2} + 2\zeta\omega_n \frac{dx}{dt} + \omega_n^2x = f(t)

Variables

SymbolDescriptionUnit
ΞΆ\zetaDamping ratio, \zeta = c / (2\sqrt{km})dimensionless
Ο‰n\omega_nNatural frequency, \omega_n = \sqrt{k/m}rad/s
f(t)f(t)Normalized external force, f(t) = F(t)/mN/kg

Spring-Mass Simulation

Interactive Simulation

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

Interact with the RLC circuit simulation below to see higher-order ODE applications.

RLC Circuit Model

Ld2qdt2+Rdqdt+1Cq=0L \frac{d^2q}{dt^2} + R \frac{dq}{dt} + \frac{1}{C} q = 0
1.0 H
1.0 Ξ©\Omega
2.0 F⁻¹
Charge q(t)
Time t = 0.0s

The Simple Pendulum

A simple pendulum of mass mm attached to a string of length LL swinging under gravity gg is governed by a non-linear differential equation.

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 where 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}.

Non-Linear Pendulum Equation

The exact differential equation governing a simple pendulum.

d2θdt2+gLsin⁑(θ)=0\frac{d^2\theta}{dt^2} + \frac{g}{L}\sin(\theta) = 0

Variables

SymbolDescriptionUnit
ΞΈ\thetaAngular displacementrad
ggAcceleration due to gravitym/sΒ²
LLLength of the pendulum stringm
ttTimes

Linearized Pendulum Equation

The small-angle approximation for a simple pendulum, resulting in simple harmonic motion.

d2ΞΈdt2+gLΞΈ=0\frac{d^2\theta}{dt^2} + \frac{g}{L}\theta = 0

Variables

SymbolDescriptionUnit
ΞΈ\thetaAngular displacement (small angles only)rad
ggAcceleration due to gravitym/sΒ²
LLLength of the pendulum stringm
ttTimes

RLC Circuits and Mechanical Analogies

An electrical circuit with a resistor (RR), inductor (LL), and capacitor (CC) connected in series is modeled by Kirchhoff's Voltage Law. The charge on the capacitor q(t)q(t) and the current i(t)i(t) exhibit behavior directly analogous to a damped mechanical system.

Analogy between Mechanical and Electrical Systems:
  • Mass m↔m \leftrightarrow Inductance LL
  • Damping c↔c \leftrightarrow Resistance RR
  • Spring Constant k↔k \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)

Kirchhoff's Voltage Law for RLC Circuits

The second-order differential equation for the charge q(t) on the capacitor.

Ld2qdt2+Rdqdt+1Cq=E(t)L \frac{d^2q}{dt^2} + R \frac{dq}{dt} + \frac{1}{C}q = E(t)

Variables

SymbolDescriptionUnit
LLInductanceH
RRResistanceΞ©\Omega
CCCapacitanceF
qqCharge on the capacitorC
E(t)E(t)Electromotive force (voltage)V
ttTimes

Current Equation for RLC Circuits

The differential equation for current i(t), found by differentiating the charge equation.

Ld2idt2+Rdidt+1Ci=dEdtL \frac{d^2i}{dt^2} + R \frac{di}{dt} + \frac{1}{C}i = \frac{dE}{dt}

Variables

SymbolDescriptionUnit
LLInductanceH
RRResistanceΞ©\Omega
CCCapacitanceF
iiCurrent, i(t) = dq/dtA
EEElectromotive force (voltage)V
ttTimes

Simple Beam Deflection

Elastic Beam Theory

The vertical deflection y(x)y(x) of an elastic beam under a distributed transverse load w(x)w(x) is governed by a fourth-order linear differential equation. 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).

Simple Beam Deflection Equation

The 4th-order differential equation governing beam deflection.

EId4ydx4=w(x)EI \frac{d^4y}{dx^4} = w(x)

Variables

SymbolDescriptionUnit
EEYoung's ModulusPa
IIArea moment of inertiam4m^4
EIEIFlexural rigidityNβ‹…m2NΒ·m^2
yyVertical deflectionm
xxDistance along the beamm
w(x)w(x)Distributed transverse loadN/m
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 (L↔m,R↔c,1/C↔kL \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).