Duhamel's Integral¶

Compute the zero-initial-condition response of a linear time-invariant vibration system to an arbitrary load history by causally convolving that load with the system's impulse response.

Version: v3 · 2026-09-19 · History
Domain-specific #: 1726
Domain group: Applied Sciences & Engineering
Origin domain: Engineering & Design (beyond software)
Subdomains: Structural Dynamics, Linear Vibration Analysis → Engineering & Design (beyond software)
Aliases: Duhamel integral, Duhamel convolution integral, Response integral

Core Idea¶

Duhamel's integral is the causal time-domain response formula that turns an arbitrary load history into the forced response of a linear time-invariant vibration system. It treats the load as a continuum of infinitesimal impulses. Each impulse excites the system's known impulse response; linear superposition adds all delayed contributions that can have reached the observation time.

For a quiescent single-degree-of-freedom system

\[ m\ddot{x}(t)+c\dot{x}(t)+kx(t)=p(t), \]

with zero initial displacement and velocity, the response is.

How would you explain it like I'm…

Adding Up Tiny Taps

Think of pushing a swing. A long push is really lots of tiny taps, one after another. Each tap makes the swing wiggle in the same known way, and you add up the wiggles from all the taps so far to know where the swing is now. That adding-up is Duhamel's integral.

Many Little Pushes Add Up

Duhamel's integral is a way engineers figure out how something that vibrates, like a building, bridge or spring, moves when it's pushed by a changing force. The idea is to chop the force into lots of tiny quick taps. Engineers already know how the object responds to one tap: it jiggles and slowly settles down. The total motion at any moment is found by adding up the jiggles from all the taps that have already happened, with older taps having had more time to fade. Taps that haven't happened yet don't count.

Impulse-Response Convolution

Duhamel's integral gives the response of a linear vibrating system, like a mass on a spring with damping, to any force that changes over time. It treats the force history as a continuous series of tiny impulses. Each impulse causes the system's known impulse response, how it moves after a single sharp kick, delayed to start at that impulse's time. Because the system is linear, you can add all those delayed responses, and only impulses before the current time count, which makes it causal. For a system starting at rest, the displacement is x(t) = integral from 0 to t of p(tau) h(t - tau) d tau, where p is the force and h the impulse response. For an underdamped spring-mass system, h is a decaying sine wave.

Duhamel's integral is the causal time-domain formula for the forced response of a linear time-invariant vibration system to an arbitrary load history. The load is treated as a continuum of infinitesimal impulses p(tau) d tau; each excites the system's impulse response h, shifted to start at tau, and linear superposition sums every delayed contribution that can have reached the observation time t. For a single-degree-of-freedom system m x'' + c x' + k x = p(t) starting from zero displacement and velocity, x(t) = integral from 0 to t of p(tau) h(t - tau) d tau, where h is the causal displacement response to a unit force impulse. In the underdamped viscous case, h(s) = (1 / (m omega_d)) e^(-zeta omega_n s) sin(omega_d s) for s >= 0, with omega_n = sqrt(k/m), zeta = c / (2 m omega_n) and omega_d = omega_n sqrt(1 - zeta^2). Substituting gives an explicit integral for the response to any forcing. The concept is not merely integrating something: it is the engineering realization of causal convolution in which the kernel is the system's impulse response and the input is the dynamic excitation.

Scope of Application¶

The formula is standard in structural dynamics, mechanical vibration, earthquake engineering, aerospace dynamics, acoustics, and linear controls. It calculates response to pulses, blast and impact histories, wind or machinery loads, measured ground acceleration, arbitrary laboratory inputs, and other nonperiodic excitations.

For base excitation, the equation is usually written in relative coordinates and the effective load is proportional to ground acceleration. For multi-degree-of-freedom systems, modal decomposition can reduce a classically damped linear system to independent modal equations; each modal coordinate receives its own Duhamel integral and the physical response is reconstructed from the modes. Matrix impulse-response formulations provide an equivalent direct statement.

Clarity¶

The abstraction separates system memory from input history. The kernel \(h\) contains the mass–damping–stiffness dynamics; \(p\) contains the event. Once \(h\) is known, changing the load does not require solving the differential equation symbolically from scratch. Conversely, changing the system changes the kernel even when the input remains fixed.

Manages Complexity¶

An arbitrary load history is infinite-dimensional, and direct differential-equation reasoning can appear to require a new solution for every waveform. Duhamel's integral compresses all system-specific dynamics into one reusable kernel. The event then enters only as weights in a causal accumulation.

This decomposition makes sensitivity legible. Increasing damping changes the kernel's decay; changing stiffness shifts its oscillation frequency; changing mass scales the response.

Abstract Reasoning¶

The formula licenses several reusable inferences:

Impulse-response sufficiency: for an LTI system, knowing \(h\) determines the zero-state response to every admissible load.
Shift rule: delaying the load delays the response by the same amount because the kernel depends on \(t-\tau\), not absolute time.
Scaling rule: multiplying the load by a constant multiplies response by that constant.
Addition rule: response to \(p_1+p_2\) equals the sum of separate responses.
Memory diagnosis: the support and decay of \(h\) state how long past inputs matter.
Resonance reasoning: a load whose time structure aligns repeatedly with the oscillatory kernel can accumulate a large response.
Initial-condition decomposition: total response equals zero-input free response plus the Duhamel zero-state response.
Modal extension: when linear equations decouple into modes, apply the same convolution mode by mode and recombine.
Failure diagnosis: amplitude-dependent frequency, path-dependent hysteresis, stiffness change, or time-varying coefficients invalidate one fixed kernel.

Knowledge Transfer¶

Reuse across linear structural, mechanical, acoustic, electrical, and control systems is literal: each retains a time-indexed input, fixed causal impulse response, superposition, and convolutional output. Only the physical interpretation of input, output, and kernel changes.

The wider transferable structure belongs to Convolution, which covers filtering, probability, imaging, and distributed-lag models. Calling every convolution a Duhamel integral would erase the engineering commitments. Conversely, recognizing a Duhamel integral as convolution allows frequency-domain multiplication, deconvolution diagnostics, kernel identification, and fast-convolution algorithms to transfer into structural dynamics without confusing them with the named response formula.

Relationships to Other Abstractions¶

Current abstraction Duhamel's Integral Domain-specific

Parents (1) — outgoing typed relations

Duhamel's Integral is a kind of Convolution Prime

The minimal prospective placement is a strict specialization of prime:convolution.

Hierarchy path (1) — routes to 1 parentless root

Duhamel's Integral → is a kind of Convolution → is a kind of Function (Mapping)

Neighborhood in Abstraction Space¶

Duhamel's Integral sits in a sparse region of the domain-specific corpus (89^th percentile for distinctiveness): few abstractions share its structure, so a faithful description tends to retrieve it precisely.

Family — Structural Vibration & Seismic Analysis (13 abstractions)

Nearest neighbors

Computed from structural-signature embeddings · 2026-10-25