In this paper we propose a Hamiltonian approach to gapped topological phases on open surfaces. Our setting is an extension of the Levin-Wen model to a 2d graph on an open surface, whose boundary is part of the graph. We systematically construct a series of boundary Hamiltonians such that each of them, when combined with the usual Levin-Wen bulk Hamiltonian, gives rise to a gapped energy spectrum which is topologically protected. It is shown that the corresponding wave functions are robust under changes of the underlying graph that maintain the spatial topology of the system. We derive explicit ground-state wavefunctions of the system on a disk as well as on a cylinder. For boundary quasiparticle excitations, we are able to construct their creation, annihilation, measuring and hopping operators etc. Given a bulk string-net theory, our approach provides a classification scheme of possible types of gapped boundary conditions by Frobenius algebras (modulo Morita equivalence) of the bulk fusion category; the boundary quasiparticles are characterized by bimodules of the pertinent Frobenius algebras. Our approach also offers a set of concrete tools for computations. We illustrate our approach by a few examples.