In the cellular networks, to support the increasing data rate requirements, many base stations (BSs) with low transmit power and small coverage area are deployed in addition to classical macro cell BSs. Low power nodes, such as micro, pico, and femto nodes (indoor and outdoor), which complement the conventional macro networks, are placed primarily to increase capacity in hotspots (such as shopping malls and conference centers) and to enhance coverage of macro cells near the cell boundary. Combining macro and small cells results in heterogeneous networks (HetNets). An accurate node (BS or user equipment (UE)) model is important in the research, design, evaluation, and deployment of 5G HetNets. The distance between transmitter (TX), receiver (RX), and interferer determines the received signal power and interference signal power. Therefore, the spatial placement of BSs and UEs greatly impacts the performance of cellular networks. However, the investigation on the spatial distribution of UE is limited, though there is ample research on the topic of the spatial distribution of BS. In HetNets, UEs tend to cluster around BSs or social attractors (SAs). The spatial distribution of these UEs is non-uniform. Therefore, the analysis of the impact of non-uniformity of UE distribution on HetNets is essential for designing efficient HetNets. This thesis presents a non-uniform user distribution model based on the existing K-tier BS distribution. Our proposed non-uniform user distribution model is such that a Poisson cluster process with the cluster centers located at SAs in which SAs have a base station offset with their BSs. There are two parameters (cluster radius and base station offset) the combination of which can cover many possible non-uniformity. The heterogeneity analysis of the proposed nonuniform user distribution model is also given. The downlink performance analysis of the designed non-uniform user model is investigated. The numerical results show that our theoretical results closely match the simulation results. Moreover, the effect of BS parameters of small cells such as BS density, BS cell extension bias factor, and BS transmit power is included. At the same time, the uplink coverage probability by the theoretical derivation is also analyzed based on some simplifying assumptions as a result of the added complexity of the uplink analysis due to the UEs’ mobile position and the uplink power control. However, the numerical results show a small gap between the theoretical results and the simulation results, suggesting that our simplifying assumptions are acceptable if the system requirement is not very strict. In addition to the effect of BS density, BS cell extension bias factor, and BS transmit power, the effect of fractional power control factor in the uplink is also introduced. The comparison between the downlink and the uplink is discussed and summarized at the end. The main goal of this thesis is to develop a comprehensive framework of the non-uniform user distribution in order to produce a tractable analysis of HetNets in the downlink and the uplink using the tools of stochastic geometry