With consideration of a high-rise coupled building system, a flexible beams-based analytical model is setup to characterize the dynamic behavior of the system. The general motion equation for the two beams interconnected by multiple viscous/visco-elastic dampers is rewritten into a non-dimensional form to identify the minimal set of parameters goving the dynamic characteristics. The corresponding exact solution suitable for arbitrary boundary conditions is presented. Furthermore, the methodology for computing the coefficients of the modal shape function is proposed. As an example, the explicit expression of the modal shape function is derived, provided only one damper is adopted to connect the adjacent buildings. Finally, to validate the proposed methodologies, three case studies are performed, in which the existence of the overdamping and the optimal damping coefficient are revealed. In the case of using one damper in connecting two similar buildings, the estimating equations for the first modal damping ratio are formulated.