AbstractA time domain boundary element method (BEM) for evaluating stresses in an axisymmetric soil mass undergoing consolidation has been developed. Previous BEM work on axisymmetric poroelasticity for boundary displacements and pore pressures is extended to permit the computation of stresses at both boundary and interior points. The stress formulation preserves the surface‐only discretization.The boundary displacement integral equation is progressively differentiated to obtain the related stress and strain integral equations. Explicit expressions for the steady‐state axisymmetric fundamental solutions are derived in this process. The transient components of the integrands are obtained directly from the transformation of the three‐dimensional kernels into a cylindrical system. Numerical implementation of these integral equations is carried out within a general purpose BEM computer code and several illustrative examples are presented to validate the method.