ODE: Periodic Forcing and Laplace Transforms

In engineering, systems from vibrating bridges to AC circuits are often driven by periodic external forces. Predicting their response is essential for stability, efficiency, and safety. The Laplace transform is a mathematical technique that converts differential equations into algebraic ones, making it exceptionally powerful for analyzing systems under periodic forcing, where the input function repeats over time.

Foundations: Periodic Forcing and Laplace Transform Essentials

A periodic function $f(t)$ satisfies $f(t+T)=f(t)$ for all $t$ and some fixed period $T>0$. Common examples include sine waves and square waves, which model alternating voltages or repetitive mechanical loads. The Laplace transform, defined as $\mathcal{L}\{f(t)\}=F(s)={\int }_0^{\infty }{e}^{-st}f(t)dt$, is your primary tool for solving linear ordinary differential equations (ODEs) with constant coefficients. It converts derivatives into algebraic terms involving the complex frequency variable $s$, incorporating initial conditions seamlessly. When facing an exam question on forced vibrations or circuits, your first step should always be to identify if the forcing function is periodic and recall that the Laplace transform method systematically handles the resulting ODE.

The Laplace Transform Formula for Periodic Functions

Instead of laboriously integrating over infinity for a periodic function, a specialized formula exists. If $f(t)$ is periodic with period $T$, its Laplace transform is given by: $$\mathcal{L}\{f(t)\}=\frac{1}{1-e^{-sT}}{\int }_0^T{e}^{-st}f(t)dt$$ This formula derives from expressing the infinite integral as a sum of integrals over successive periods. You only need to integrate over one period, from $0$ to $T$, and then divide by the factor $(1-e^{-sT})$. For instance, consider a simple periodic square wave where $f(t)=1$ for $0<t<a$ and $f(t)=0$ for $a<t<T$ within one period. Applying the formula: $$\mathcal{L}\{f(t)\}=\frac{1}{1-e^{-sT}}{\int }_0^a{e}^{-st}(1)dt=\frac{1}{1-e^{-sT}}\cdot \frac{1-e^{-sa}}{s}$$ A common exam trap is misapplying this formula to non-periodic functions or incorrectly identifying the period $T$. Always verify periodicity before using it.

Step-by-Step Solution of ODEs with Periodic Forcing

Consider a standard second-order system, like a spring-mass-damper with periodic external force $F(t)$: $$m{y}^{\prime \prime }(t)+c{y}^{\prime }(t)+ky(t)=F(t)$$ with initial conditions $y(0)$ and $y^{\prime }(0)$. The method proceeds in four clear steps:

1. Take the Laplace transform of both sides. Using linearity and derivative properties: $\mathcal{L}\{y^{\prime \prime }\}=s^{2}Y(s)-sy(0)-y^{\prime }(0)$ and $\mathcal{L}\{y^{\prime }\}=sY(s)-y(0)$. The right-hand side requires transforming $F(t)$. If $F(t)$ is periodic, use the formula from the previous section.
2. Solve algebraically for $Y(s)$. This yields an expression like $Y(s)=\frac{\mathcal{L}\{F(t)\}+(ms+c)y(0)+m{y}^{\prime }(0)}{m{s}^2+cs+k}$.
3. Perform partial fraction decomposition. Break $Y(s)$ into simpler terms whose inverse transforms are known. This is where the response splits naturally into components.
4. Take the inverse Laplace transform to obtain the time-domain solution $y(t)$.

For example, let $m=1,c=2,k=2$, with zero initial conditions, and a sinusoidal force $F(t)=\sin (t)$, which is periodic. Its Laplace transform is $\mathcal{L}\{\sin (t)\}=\frac{1}{s^2+1}$. The transformed equation becomes $(s^2+2s+2)Y(s)=\frac{1}{s^2+1}$. Solving: $$Y(s)=\frac{1}{(s^2+1)(s^2+2s+2)}$$ After partial fractions, the inverse transform gives the full solution. In exam settings, show every algebraic step clearly; a frequent error is making sign mistakes during partial fraction expansion.

Analyzing Steady-State and Transient Responses

The complete solution $y(t)$ from the Laplace method contains two distinct parts: the transient response and the steady-state response. The transient response consists of terms that decay to zero as $t\to \infty$, typically arising from the system's homogeneous solution and depending on initial conditions. The steady-state response is the persistent, long-term behavior that mirrors the frequency of the periodic forcing function. In the $Y(s)$ expression, terms with denominators containing the system's characteristic polynomial (like $s^2+2s+2$) often lead to transient terms after inversion, while terms stemming from the transform of the forcing function (like $1/(s^2+1)$) yield the steady-state. For the sinusoidal forcing example above, the steady-state is a sinusoidal oscillation at the driving frequency. Recognizing this separation is crucial for engineering design, where the transient might indicate startup stresses, and the steady-state determines operational performance.

Resonance and Practical Applications in Vibrating Systems

Resonance occurs when the frequency of the periodic forcing function matches or nears the natural frequency of the system, leading to dramatically amplified oscillations. For the ODE $m{y}^{\prime \prime }+c{y}^{\prime }+ky=F_0\sin (\omega t)$, the natural frequency is ${\omega }_n=\sqrt{k/m}$. The steady-state amplitude peaks when the driving frequency $\omega$ is close to ${\omega }_n$, especially with light damping (small $c$). In the Laplace domain, this condition manifests as a denominator term approaching zero, making the corresponding coefficient in the partial fraction expansion very large. For instance, if the characteristic polynomial has roots close to $\pm i\omega$, resonance effects dominate.

This has direct applications to vibrating systems with periodic excitation, such as aircraft wings buffeted by wind, rotating machinery with imbalance, or buildings during earthquakes. Engineers use Laplace transform analysis to predict resonant frequencies and design systems to avoid them, perhaps by adding dampers or tuning mass and stiffness. In exam problems, you might be asked to find the driving frequency that maximizes amplitude or to interpret a given transfer function's poles in the s-plane. Always check for potential resonance when the forcing frequency is variable.

Common Pitfalls

1. Incorrectly applying the periodic transform formula: Using the formula $\frac{1}{1-e^{-sT}}{\int }_0^T{e}^{-st}f(t)dt$ for a function that is not truly periodic over all $t>0$ is a critical error. Correction: Confirm $f(t+T)=f(t)$ for all $t$ before proceeding. Sketch the function if unsure.

1. Neglecting initial conditions: When taking the Laplace transform of derivatives, it's easy to forget the terms involving $y(0)$ and $y^{\prime }(0)$, especially under exam pressure. Correction: Always write out the derivative properties explicitly: $\mathcal{L}\{y^{\prime }\}=sY(s)-y(0)$ and $\mathcal{L}\{y^{\prime \prime }\}=s^{2}Y(s)-sy(0)-y^{\prime }(0)$.

1. Confusing transient and steady-state components: Students sometimes misidentify which part of the solution persists. Correction: After inverse transforming, identify terms with exponential decay (e.g., $e^{-at}$ where $a>0$) as transient. Terms that are sinusoidal or constant constitute the steady-state.

1. Overlooking resonance conditions in damped systems: For lightly damped systems, resonance occurs near, but not exactly at, the undamped natural frequency. Correction: The peak amplitude for the system $m{y}^{\prime \prime }+c{y}^{\prime }+ky=F_0\sin (\omega t)$ occurs at ${\omega }_r={\omega }_n\sqrt{1-2{\zeta }^2}$, where $\zeta =c/(2\sqrt{mk})$ is the damping ratio. Remember to calculate this if damping is present.

Summary

• The Laplace transform of a periodic function with period $T$ is efficiently computed using $\mathcal{L}\{f(t)\}=\frac{1}{1-e^{-sT}}{\int }_0^T{e}^{-st}f(t)dt$, requiring integration over only one period.
• Solving ODEs with periodic forcing involves transforming the entire equation, solving algebraically for $Y(s)$, performing partial fraction decomposition, and inverting back to the time domain.
• The total solution separates into a transient response that decays over time and depends on initial conditions, and a steady-state response that persists and matches the forcing frequency.
• Resonance—a large amplitude response—occurs when the forcing frequency nears the system's natural frequency, a critical consideration in designing vibrating systems like machinery or structures.
• This Laplace transform framework provides a systematic, powerful method for analyzing and predicting the behavior of engineering systems under periodic excitation, from conceptual understanding to quantitative design.