Due to the unique magnetic, mechanical and thermal properties, magnetic nanoparticles (MNPs) have comprehensive applications as the contrast and therapeutic agents in biomedical imaging and magnetic hyperthermia. The linear and nonlinear magnetoacoustic responses determined by the magnetic properties of MNPs have attracted more and more attention in biomedical engineering. By considering the relaxation time of MNPs, we derive the formulae of second harmonic magnetoacoustic responses (2H-MARs) for a cylindrical MNP solution model based on the mechanical oscillations of MNPs in magnetoacoustic tomography with magnetic induction (MAT-MI). It is proved that only the second harmonic magnetoacoustic oscillations can be generated by MNPs under an alternating magnetic excitation. The acoustic pressure of the 2H-MAR is proportional to the square of the magnetic field intensity and exhibits a linear increase with the concentration of MNPs. Numerical simulations of the 2H-MAR are confirmed by the experimental measurements for various magnetic field intensities and solution concentrations using a laser vibrometer. The favorable results demonstrate the feasibility of the harmonic measurements without the fundamental interference of the electromagnetic excitation, and suggest a new harmonic imaging strategy of MAT-MI for MNPs with enhanced spatial resolution and improved signal-to-noise ratio in biomedical applications.