Hybrid beamforming via large antenna arrays has a great potential for increasing data rate in cellular networks by delivering multiple data streams simultaneously. In this paper, several beamforming design algorithms are proposed based on long-term channel information in macro-cellular environments where the base station is equipped with a massive phased array under per-antenna power constraint. Using an adaptive scheme, beamforming vectors are updated whenever the long-term channel information changes. First, the problem is studied when the base station has a single RF chain (single-beam scenario). Semi-definite relaxation (SDR) with randomization is used to solve the problem. As a second approach, a low-complexity heuristic beam composition algorithm is proposed which performs very close to the upper-bound obtained by SDR. Next, the problem is studied for a generic number of RF chains (multi-beam scenario) where the Gradient Projection method is used to obtain local solutions. Numerical results reveal that using massive antenna arrays with optimized beamforming vectors can lead to five-fold network throughput improvement over systems with conventional antennas.