A multiple monopole method based on the generalized multipole technique is presented for the calculation of band structures of two-dimensional mixed solid/fluid phononic crystals. In this method, the fields are expanded by using the fundamental solutions with multiple origins. Besides the sources used to expand the wave fields, an extra monopole source is introduced as the external excitation. By varying the frequency of the excitation, the eigenvalues can be localized as the extreme points of an appropriately chosen function. By sweeping the frequency range of interest and sweeping the boundary of the irreducible first Brillouin zone, the band structure of the phononic crystals can be obtained. The method can consider the fluid-solid interface conditions and the transverse wave mode in the solid component strictly. Some typical examples are illustrated to discuss the accuracy of the present method. multiply connected domains. The layer potential techniques, which are in the same spirit of the GMT, were developed by Ammari et al. [22][23][24] for studying the wave scattering from high-contrast microstructures, where the spectral properties exhibit a bandgap opening.Recently, the GMT has been developed to calculate the bandgaps of scalar waves in 2D phononic crystals by the present authors [25] based on the method of Reutskiy. Phononic crystals are artificial inhomogeneous composite materials constructed by scatterers periodically embedded in a homogeneous host material and have received a great deal of attention in the last decades [26,27]. The most interesting and important feature of phononic crystals is that such artificial composites can exhibit elastic/acoustic wave bandgaps, in which wave propagation and vibration are all forbidden regardless of the polarization and propagating direction of the elastic/acoustic waves. This behavior leads to many interesting physical properties and potential applications of phononic crystals such as acoustic filters, control of vibration isolation, noise suppression, and design of new transducers [28,29]. In order to determine the band structures and search for the so-called complete phononic bandgaps, an eigenvalue problem has to be solved. For this purpose, several numerical methods have been proposed to compute the band structures. However, some of them have difficulties in dealing with the mixed fluid-solid systems [30,31]. This is due to the fact that different wave modes propagate in solids and fluids, respectively. The mixed longitudinal and transverse elastic waves can propagate in solids, while only the purely longitudinal acoustic wave can propagate in fluids. So the solid-fluid interaction must be taken into account at the interface between the solid and the fluid, which makes it difficult to compute the band structures of phononic crystals composed of solid and fluid components. The widely used plane wave expansion method [26,27] and the recently developed wavelet method [32] were applied to calculate the band structures of mixed solid/fluid phononic crystals by ...