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: v2 · 2026-08-30 · History
Domain-specific #: 1726
Origin domain: structural dynamics
Subdomain: linear vibration analysis
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.

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 prime: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) — more general patterns this builds on

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 → Convolution → 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 — Unclustered & Miscellaneous (1565 abstractions)

Nearest neighbors

Computed from structural-signature embeddings · 2026-09-08