In this work, ferroelastic domain switching and R‐curve behavior in lead zirconate titanate (Nb/Ce co‐doped Pb(Zr0.52Ti0.48)O3, ab. PZT‐NC)‐based ferroelectric ceramics were investigated, using the indentation‐strength‐in‐bending (ISB) method. Firstly, Vickers indentation test examined the notable fracture anisotropy of PZT‐NC ceramics between the poling direction and its perpendicular direction, and the crack open displacement (COD) profiles in the two directions were also theoretically calculated from the indentation fracture mechanics. And then two kinds of ferroelastic domain switching modes (in‐plane and out‐of‐plane) were used for explaining such anisotropic propagation behavior of indentation cracks. The subsequent three‐point bending test illustrated the dependence of fracture strength on indentation load and the rising crack growth resistance curves (R‐curves) in two directions. The resulted R‐curves were fitted by the Hill's type Growth Function successfully, giving the reasonable values of crack extension exponential (n), plateau fracture toughness (Kmax), and initial fracture toughness (Kini). The in‐plane ferroelastic domain switching was identified as a more significant toughening mechanism for PZT‐NC ceramics than the out‐of‐plane switching due to more switchable domains.