Abstract: A precise algorithm is proposed to evaluate the dynamic stiffness matrix of the rigid strip foundation embedded in layered transversely isotropic soil. The space transformation method is used to obtain the dynamic displacement response of arbitrary nodes at the surface or in the layered transversely isotropic soil. Then the dynamic flexibility matrix of the nodes in the excavation scope is established. Finally the flexible volume method is introduced to evaluate the dynamic stiffness matrix of the rigid strip foundation. The proposed method solves the dynamic flexibility coefficients of the layered soil in the frequency-wave number domain by using the precise integration method, and there is no limit to the number and thickness of layers. In addition, the computation is based on the matrix algebra with small dimension, and the computation is always stable with high efficiency. Numerical examples demonstrate the accuracy and the wide applicability of the proposed method to the layered transversely isotropic soil.