<p>This study addresses the anti-plane steady-state SH wave scattering problem in circular laminated structures with circular hole defects. An analytical method integrating multipolar coordinate transformation, Graf’s addition theorem, and wave function expansion is proposed, constructing for the first time an asymmetric Green’s function solution under the synergistic effect of bi-material media and defects. This approach overcomes the theoretical limitations of traditional one-phase media models in characterizing interfacial wave coupling effects. Based on the governing Helmholtz equations, a full-domain wavefield coupling model for bi-material media interfaces and defect boundaries is established using complex-plane coordinate transformation and multipolar coordinate expansion techniques. Integral transforms and Graf’s formula are introduced to rigorously satisfy the interfacial stress-displacement continuity conditions and the stress-free boundary conditions near the circular hole defects. The Green’s function solution set with explicit physical significance is ultimately obtained by solving the coefficient equations. A predictive model for the circumferential dynamic stress concentration factor (DSCF) is derived from the analytical solution, revealing the multiscale modulation mechanisms of shear modulus ratio, defect eccentricity, and incident wavenumber. Numerical validation shows that the absolute error between the maximum DSCF predicted by this method and the degenerated result is 0.74 across a wide frequency range. The research provides a high-precision analytical tool for SH wave scattering analysis in defective laminated structures, offering theoretical value and engineering potential for composite damage assessment and aerospace structural health monitoring.</p>