Stationary solutions for two layer lubrication equations We investigate stationary solutions of flows of thin liquid bilayers in an energetic formulation which is motivated by the gradient flow structure of its lubrication approximation. The corresponding energy favors the liquid substrate to be only partially covered by the upper liquid. This is expressed by a negative spreading coefficient which arises from an intermolecular potential combining attractive and repulsive forces and leads to an ultra thin layer of thickness e. For the corresponding lubrication models existence of stationary solutions is proven. In the limit e to 0 matched asymptotic analysis is applied to derive sharp interface models and the corresponding contact angles i.e. the Neumann triangle. In addition we use G convergence and derive the equivalent sharp interface models rigorously in this limit. For the resulting model existence and uniqueness of energetic minimizers are proven. The minimizers agree with solutions obtained by matched asymptotics.