The interactive dynamics of two relativistic laser beams copropagating in underdense plasmas is studied using a coupled model equation for the relativistic laser propagation. It is shown that the relative phase difference between the two laser pulses plays a significant role on their interaction processes. When the relative phase varies, the two laser beams display different features, such as attraction, repulsion, and energy shift. Particularly, energy flow from the phase-advanced beam to the spot domain of the phase-delayed beam is observed when the relative phase difference is between 0 and π. When the relative phase is larger than π/2, repulsion is dominant and the interaction gradually becomes weak. When relative phase difference is smaller than π/2, attraction becomes dominant and as the phase difference decreases, the phase-advanced beam shifts its most energy into the spot domain of the phase-delayed beam. These conclusions are verified by our three dimensional particle-in-cell simulations. This provides an efficient way to manipulate the energy distribution of relativistically intense laser pulses in plasmas by adjusting their relative phase.