  Amplitude expansions are used to determine steady states of a semi-infinite
solid subject to the Grinfeld instability in systems with a fixed (wave)length.
We present two methods to obtain high-order weakly nonlinear results. Using the
system size as a control parameter, we circumvent the problem that there is no
instability threshold for an extended system in the absence of gravity. This
way, the case without gravity becomes accessible to a weakly nonlinear
treatment. The dependence of the branch structure of solution space on the
level of gravity (or density difference) is exhibited. In the zero-gravity
limit, we recover the solution branch obtained by Spencer and Meiron. A
transition from a supercritical to a subcritical bifurcation is observed as
gravity is increased or the nonhydrostatic stress is decreased at fixed
gravity. At given values of the system parameters, we find a discrete, possibly
infinite, set of solution branches. This is reminiscent of dendritic or
eutectic growth, where similar solution sets exist, of which only a particular
one is linearly stable. Despite the high order of our expansions, the approach
is restricted to relatively small nondimensional amplitudes ($\lesssim 0.2$), a
disadvantage we can overcome by a variational approach that will be discussed
in a companion paper. At the critical point, we find that not only the first
Landau coefficient is negative but all of them up to the highest amplitude
order (15) we could compute so far.
