An H_∞ synthesis method for control of a flexible joint, with non-linear spring characteristic, is proposed. The first step of the synthesis method is to extend the joint model with an uncertainty description of the stiffness parameter. In the second step, a non-linear optimisation problem, based on nominal performance and robust stability requirements, has to be solved. Using the Lyapunov shaping paradigm and a change of variables, the non-linear optimisation problem can be rewritten as a convex, yet conservative, LMI problem. The method is motivated by the assumption that the joint operates in a specific stiffness region of the non-linear spring most of the time, hence the performance requirements are only valid in that region. However, the controller must stabilise the system in all stiffness regions. The method is validated in simulations on a non-linear flexible joint model originating from an industrial robot.