We use the incremental harmonic balance (IHB) method to analyse the dynamic stability problem of a nonlinear multiple-nanobeam system (MNBS) within the framework of Eringen's nonlocal elasticity theory. The nonlinear dynamic system under consideration includes MNBS embedded in a viscoelastic medium as clamped chain system, where every nanobeam in the system is subjected to time-dependent axial loads. By assuming the von Karman type of geometric nonlinearity, a system of m nonlinear partial differential equations of motion is derived based on the Euler-Bernoulli beam theory and D' Alembert's principle. All nanobeams in MNBS are considered with simply supported boundary conditions. Semi-analytical solutions for time response functions of the nonlinear MNBS are obtained by using the single-mode Galerkin discretization and IHB method, which are then validated by using the numerical integration method. Moreover, Floquet theory is employed to determine the sta- method, we obtain an incremental relationship with the frequency and amplitude of time-varying axial load, which defines stability boundaries. Numerical examples show the effects of different physical and material parameters such as the nonlocal parameter, stiffness of viscoelastic medium and number of nanobeams on Floquet multipliers, instability regions and nonlinear amplitude-frequency response curves of MNBS. The presented results can be useful as a first step in the study and design of complex micro/nanoelectromechanical systems.