<p>We solve the problem of finding the two-dimensional dynamic stress state for an elastic isotropic body containing a system of arbitrarily arranged thin rigid inclusions. The inclusions are subjected to the action of normal and shear forces and moments varying as harmonic functions of time. The original problem is reduced to a system of integral equations for the jumps of stresses on the surfaces of inclusions. To solve this system, we propose to use an iterative method, which enables us to avoid the necessity of solving high-dimensional systems of integral equations. The presented examples demonstrate the convergence and stability of the proposed method also in the case of arbitrary systems of inclusions of complex configurations.</p>