In this paper, an accurate series solution for the longitudinal vibration analysis of elastically coupled nanorods system is established, in which artificial springs are introduced to simulate such general coupling and boundary conditions. Energy formulation is derived for the description of axial dynamics of multiple coupled nanorods based on Eringen nonlocal elasticity. For each nanorod component, its longitudinal vibration displacement function is invariantly assumed as the superposition of Fourier series and boundary smoothed supplementary polynomials, with the aim to make all the spatial differential sufficiently continuous across each rod. All the unknown coefficients are determined in conjunction with Rayleigh-Ritz procedure to solve a standard eigenvalue matrix. Numerical results are presented for the particular cases of twonanorod and three-nanorod systems to demonstrate the correctness and effectiveness of the proposed model. Excellent agreements can be repeatedly observed in comparison with those from other approaches available in literature. Based on the model established, influence of elastic boundary and coupling conditions on the longitudinal modal characteristics of multiple nanorods system is investigated and addressed. This work can provide an efficient energy-based modeling framework for longitudinal vibration analysis of elastically connected nanorod system with complicated boundary conditions.