The performances of piezoelectric PN junction photovoltaic cells are closely related to the internal potential barrier configurations and the distributions of carriers,and can be tuned through carrier transport characteristic changes by the piezopotentials under the piezo-effect.However,the classical PN junction model fails to describe the coupling effect between multiple physical fields and carriers in the potential barrier zone due to the depletion layer assumptions and others,and in turn gives severely distorted results.Herein a me-chanics-electricity-photonics-carrier global multi-field coupling model was developed to investigate the tuning mechanism for mechanical loadings on the output characteristics of ZnO photovoltaic cells.The numerical re-sults indicate that,the short-circuit current,the open-circuit voltage,and the maximum output power of the photovoltaic cell increase with the compressive stress under a fixed light intensity,while tensile stresses are not conducive to improving the performances of photovoltaic cells.In addition,a better tuning effect occurs with a loading region wider than the illuminated region,or with both these two external fields acting in the same side of the n/p-zone.