Abstract A modified incremental harmonic balance (IHB) method is used to determine periodic solutions of wave propagation in discrete, strongly nonlinear, periodic structures, and solutions are found to be in a two-dimensional hyperplane. A novel method based on the Hill’s method is developed to analyze stability and bifurcations of periodic solutions. A simplified model of wave propagation in a strongly nonlinear monatomic chain is examined in detail. The study reveals the amplitude-dependent property of nonlinear wave propagation in the structure and relationships among the frequency, the amplitude, the propagation constant, and the nonlinear stiffness. Numerous bifurcations are identified for the strongly nonlinear chain. Attenuation zones for wave propagation that are determined using an analysis of results from the modified IHB method and directly using the modified IHB method are in excellent agreement. Two frequency formulae for weakly and strongly nonlinear monatomic chains are obtained by a fitting method for results from the modified IHB method, and the one for a weakly nonlinear monatomic chain is consistent with the result from a perturbation method in the literature.