Proliferating epithelia are sheets of dividing cells that tightly adhere to each other in many metazoans. Regulation of cell division is important for both normal and diseased epithelia. Quantitative analysis of epithelia shows significant conservation in the dominance of hexagonal cells in both animal and plant. However, there exist species with significant varied amount of hexagonal cells. The mechanism controlling these distributions is poorly understood. In this study, we use a computational method called cellGPS (cell global pattern simulator) to study the dynamic process of proliferating epithelia based on a mechanical off-lattice cellular model. This method keeps track of full geometric and biomechanical properties of epithelia. It models cell shape, area, perimeters and surface tension. Two biologically relevant division rules, namely, the orthogonal orientation found in plants and the largest-side orientation found in the mitotic spindle in human cells, are also incorporated in the division process. With this method, we found that regardless of division rules, our mechanical model can produce the common epithelia pattern of the dominance of hexagonal cells. In addition, we found that division rule of largest-side orientation lead to the formation of cell pattern with around 45% hexagons, which is observed in many animal epithelia and some plant tissues. This suggests that reducing stress on the cell wall by dividing the largest side is widely used in animal epithelia. The percentage of hexagons from division following orthogonal orientation is found to be about 60%, consistent with observations found in certain plant tissues. We conclude that both mechanical force and division orientation play important roles during epithelial proliferation. Mechanical model itself leads to topological dominance of hexagonal cells, and different rules of division orientation lead to varying amount of hexagonal cells.