Computationally-useful methods of estimating the critical time step for linear triangular elements and for linear quadrilateral elements are given. Irregular nodal-point arrangements, position-dependent properties, and a variety of boundary conditions can be accommodated. The effects of boundary conditions and element shape on the critical time step are discussed. Numerical examples are presented to illustrate the effect of various boundary conditions and for comparison to the finite-difference method.