A 3D time-harmonic problem for an infinite elastic matrix with an arbitrarily located interacting rigid disk-shaped inclusion and a penny-shaped crack is analyzed by the boundary integral equation method. Perfect bonding between the matrix and the moving inclusion is assumed. The crack faces are subjected to time-harmonic loading. The boundary integral equations (BIEs) obtained are solved numerically by the implementation of regularization and discretization procedures. Numerical calculations are carried out for a crack under tensile loading of constant amplitude, where an interacting inclusion is perpendicular to the crack and has the same radius. Both the normal crack-opening-displacement and the mode-I stress intensity factor are investigated for different wave numbers and distances between the crack and the inclusion.