Aimed at meshing the ferromagnetic objects with both surface elements and volume elements, a magnetic integral equation method was proposed to solve the induced magnetic field of complex-geometry carriers. Numerical results of the shell sphere and solid sphere model obtained by the proposed method can match its analytical solution well, and the ARE (average relative error) is less than 0.27%; the induced magnetic field of the ship model can be modeled well with few hybrid elements, which demonstrates the advantage of the proposed method. Based on the proposed method, the AREs between numerical and experimental results of the simple submarineequipment mockup model and the submarine mockup model are both less than 2%. Magnetism deviations caused by the equipment model are displayed by the comparison between computational results of the simple submarineequipment mockup model and the simple submarine mockup model. Thus carrier′s induced magnetic field can be calculated by the proposed method effectively, which is practically valuable.