Main objective of this study is to investigate the thermal behavior of ceramic laminate veneer restorations of the maxillary central incisor with different incisal preparations such as butt joint and palatinal chamfer using fmite element method. In addition, it is also aimed to understand the effect of different thermal loads which simulates hot and cold liquid imbibing in the mouth. Three-dimensional solid models of the sound tooth and prepared veneer restorations were obtained using micro-computed tomography images. Each ceramic veneer restoration was made up of ceramic, luting resin cement and adhesive layer which were generated based on the scanned images using computer-aided design software. Our solid model also included the remaining dental tissues such as periodontal ligament and surrounding cortical and spongy bones. Time-dependent linear thermal analyses were carried out to compare temperature changes and stress distributions of the sound and restored tooth models. The liquid is firstly in contact with the crown area where the maximum stresses were obtained. For the restorations, stresses on palatinal surfaces were found larger than buccal surfaces. Through interior tissues, the effect of thermal load diminished and smaller stress distributions were obtained near pulp and root-dentin regions. We found that the palatinal chamfer restoration presents comparatively larger stresses than the butt joint preparation. In addition, cold thermal loading showed larger temperature changes and stress distributions than those of hot thermal loading independent from the restoration technique.