Effective and reliable clinical uses of dental ceramics necessitate an insightful analysis of the fracture behaviour under critical conditions. To better understand failure characteristics of porcelain veneered to zirconia core ceramic structures, thermally induced cracking during the cooling phase of fabrication is studied here by using the extended finite element method (XFEM). In this study, a transient thermal analysis of cooling is conducted first to determine the temperature distributions. The time-dependent temperature field is then imported to the XFEM model for viscoelastic thermomechanical analysis, which predicts thermally induced damage and cracking at different time steps. Temperature-dependent material properties are used in both transient thermal and thermomechanical analyses. Three typical ceramic structures are considered in this paper, namely bi-layered spheres, squat cylinders and dental crowns with thickness ratios of either 1:2 or 1:1. The XFEM fracture patterns exhibit good agreement with clinical observation and the in vitro experimental results obtained from scanning electron microscopy characterization. The study reveals that fast cooling can lead to thermal fracture of these different bi-layered ceramic structures, and cooling rate (in terms of heat transfer coefficient) plays a critical role in crack initiation and propagation. By exploring different cooling rates, the heat transfer coefficient thresholds of fracture are determined for different structures, which are of clear clinical implication