Abstract : Two approaches for sound field prediction of a spherical loudspeaker array operation are presented and compared to experimental measurements. In the first, real membrane modes from experimental modal analysis are used as input for BEM vibroacoustic simulations using LMS Virtual.Lab® software. In the second, an analytical solution based on the spherical harmonic expansion of an idealized velocity pattern over the spherical array is used. Results are presented in terms of sound power and directivity patterns, showing that the former has the same trend in all comparisons for the low frequency range, and that the analytical solution cannot be used for the high frequency range since it does not include the effect of the flexible membrane modes. Directivity patterns, however, show a good degree of similarity in all cases, and are not strongly affected by the flexible membrane modes. The differences found in the results and the unrealistic amplifications in the sound power curves from the simulations are caused mainly by three factors: acoustic cavity modes of the array, neglecting the acoustic coupling among loudspeakers for the operating condition and utilization of the same set of FRFs for all loudspeakers in the array. In a general way, it can be said that the flexible membrane modeling improves considerably the prediction of radiated sound power in the high frequency range, which cannot be obtained by the analytical model commonly used in the analysis of spherical loudspeaker arrays.