The objective of this study is to develop three dimensional finite elements models to investigate the effect of many factors on the behavior of dowel Jointed Portland Cement Concrete Pavement (JPCCP). The 3D-FE model is applied to analyze the critical bending stresses and deflections for understanding some modes of distress. Moreover, developed equation for predicting the non-uniform distribution of the modulus of subgrade reaction beneath concrete slab subjected to concentrated load is achieved. The accuracy of the developed equation is verified by a comparison with winkler approach such as SAP program. Winkler's solutions are found to overestimate the maximum deflections by about 7.5% more than 3D-FE solutions up to a horizontal distance of 60 inch from the load center. In general, the results show that the developed 3D-FE is suitable for identifying the effect of different design features on the structural response of rigid pavements.