SUMMARYThe stress model of the hybrid‐Trefftz finite element is formulated for the analysis of elastodynamic problems defined on unsaturated porous media. The supporting mathematical model is the theory of mixtures with interfaces and considers the full coupling between the solid, fluid and gas phases, including the effect of seepage acceleration. Hybrid‐Trefftz stress elements use the free‐field regular solutions of the homogeneous Navier (or Beltrami) equation to construct the approximation of the generalized stresses in the domain of the element. The influence of non‐homogeneous terms in the Navier equation is modelled using solutions of the corresponding static problem. The resulting elements are highly convergent under p‐refinement and robust to both low and high excitation frequencies, as the trial functions embody relevant physical information on the modelled phenomenon. Copyright © 2013 John Wiley & Sons, Ltd.