We develop analytical approaches and analyze their range of applicability to oscillators with two nonlinear springs in parallel and series connections. Specifically, we focus our study on systems with hardening springs and cubic nonlinearities. In both cases, three dimensionless parameters govern the oscillator namely lambda = k(2)/k(1) and epsilon(1,2) = epsilon(1,2)A(2). Here, k(1,2) and epsilon(1,2) are the linear stiffness and the nonlinearity coefficient of both springs, respectively, and A is the amplitude of the position of the mass, in the parallel case, or the deflection of the spring connected to the mass (k(2), epsilon(2)), in the series case. It is found that, in parallel configuration, for lambda > 0 and 0 < epsilon(1,2) <= 1 the analytical solution gives an excellent approach to the exact solution found numerically. However, in series connection the numerical simulations show that the solution of the oscillator becomes much more complex than in parallel connection, and the analytical approaches work excellently in the ranges 0 < lambda <= 1, 0 < epsilon(1,2) <= 0.1, and, 0 < lambda <= 0.1, 0 < epsilon(1,2) <= 1. (C) 2015 Elsevier Ltd. All rights reserved.