This paper describes a numerical study of the punching shear resistance of unbonded post-tensioned slabs without shear reinforcement. This research aimed to develop a methodology for modeling unbonded tendons and numerically evaluate the prestressing effects on the punching shear capacity. To validate the modeling approach, a series of well documented experimental tests were simulated using the finite element software DIANA. The nonlinear analyses were performed using three-dimensional solid elements, incorporating the cracking behavior of concrete by the smeared crack approach. In addition, interface elements were included, providing bond-slip properties to simulate the friction between tendons and concrete. A good agreement was found between the numerical results and experimental data. Load capacity, cracks patterns, and the prestressing effects were accurately simulated. After the validation, a parametric study was conducted to analyze the influence of distribution, force and profile of prestressing tendons. Finally, the numerical results were compared with current design code provisions and the approach provided by the Critical Shear Crack Theory.