In previous work, we have performed amplitude expansions of the continuum equations for the Grinfeld instability and carried them to high orders. Nevertheless, the approach turned out to be restricted to relatively small amplitudes. In this article, we use a variational approach in terms of multi-cycloid curves instead. Besides its higher precision at given order, the method has the advantages of giving a transparent physical meaning to the appearance of cusp singularities and of not being restricted to interfaces representable as single-valued functions. Using a single cycloid as ansatz function, the entire calculation can be performed analytically, which gives a good qualitative overview of the system. Taking into account several but few cycloid modes, we obtain remarkably good quantitative agreement with previous numerical calculations. With a few more modes taken into consideration, we improve on the accuracy of those calculations. Our approach extends them to situations involving gravity effects. Results on the shape of steady-state solutions are presented at both large stresses and amplitudes. In addition, their stability is investigated.