We study the effects of the MSSM contribution on B -> phi K decays using the perturbative QCD approach. In this approach, strong phases can be calculated, so that we can predict the values of CP asymmetries with the MSSM contribution. We predict a large relative strong phase between the penguin amplitude and the chromomagnetic penguin amplitude. If there is a new CP violating phase in the chromomagnetic penguin amplitude, then the CP asymmetries may change significantly from the SM prediction. We parametrize the new physics contributions that appear in the Wilson coefficients. We maximize the new physics parameters up to the point where it is limited by experimental constraints. In the case of the LR insertion, we find that the direct CP asymmetries can reach about 85% and the indirect CP asymmetry can reach about -30%.