This paper studies the combinatorial Yamabe flow on hyperbolic surfaces with boundary. It is proved by applying a variational principle that the length of boundary components is uniquely determined by the combinatorial conformal factor. The combinatorial Yamabe flow is a gradient flow of a concave function. The long time behavior of the flow and the geometric meaning is investigated.