An efficient mesh model repair method is proposed in this paper, which can repair the holes with boundary continuity. It is mainly composed of the following steps: Firstly, the hole boundary is triangulated to generate the initial surface patch, of which the edges are distributed along the boundary averagely similar to an curtain covering at a window. Secondly, the initial surface patch is refined, and the interior edges are relaxed while refining to get a Delaunay-like refinement surface patch, which approximates the density of the surrounding mesh. Finally, the refinement surface patch is reshaped based on the minimization of the surface energy solved by the order discrete Euler-Lagrange equation in order to realize the continuity repair. The proposed method is proved to be efficient and robust by experiments of different kinds of boundary conditions.