This paper investigates the mass-sensing capability of an array of a few identical electrostatically actuated microbeams, as a first step toward the implementation of arrays of thousands of such resonant sensors. A reducedorder model is considered and Taylor series are used to simplify the nonlinear electrostatic force. Then the Harmonic Balance Method associated with the Asymptotic Numerical Method, as well as time-integration or averaging methods are applied to this model and their results are compared. In this paper, two-and three-beam arrays are studied. The predicted responses exhibit complex branches of solutions with additional loops due to the influence of adjacent beams. Moreover, depending on the applied voltages, the solutions with and without added mass exhibit large differences in amplitude which can be used for detection. For symmetric configurations, the symmetry breaking induced by an added mass is exploited to improve mass sensing.