Martensitic phase transformation of zirconia (ZrO2), when realized from the superelastic effect, has immense potential in energy applications, such as micro‐actuation and high‐energy dissipation/damping. In this study, martensitic transformation in ZrO2 is studied using a newly developed microscale phase field model, which is scale‐independent and contains an averaged gradient term. The model is implemented using the FEniCS package in Python. With a thorough explanation of the model construction, the mechanics of phase transformation are detailed, and the material behaviors under different stress conditions are discussed. Mesh sensitivity, multivariant evolution, and the effect of particle size are analyzed to understand the martensite evolution in superelastic ZrO2. This work provides new insight into the martensitic phase transformation behaviors of micron‐sized spherical ZrO2 particles.