In this paper, the recursive T-matrix formalism is developed to determine the dynamic Green's function for SH-wave in the multi-dipping layers with arbitrary shape embedded in a half-space. In each layer, the wave field generated by a harmonic line load can be separated into two parts: the source term and the complementary part. The source term is exactly the Green's function for SH-wave in two-dimensional half-space, and the complementary part which causes the SH-wave to satisfy the boundary conditions at the interface is determined by the wave function expansion method. Using the polar coordinates, the basis functions are constructed by Bessel typed cylindrical functions, and their orthogonality conditions at the corresponding interfaces can be derived by means of Betti's third identity. Applying the conditions of continuity at the interface, the recursive T-matrix formalism is developed for determining the expansion coefficients of the complementary part from the associated source dependent constants.
