Multiple crack simulations were conducted on the delamination process of thermal barrier coating applied to gas turbine blades using the polycrystal FEM model of thermal cycle test pieces. The simulation results of multiple cracks during the thermal cycle showed quite similar morphology to the observation on test pieces by modifying the Young's modulus of topcoat layer according to the fraction of pore. Moreover, statistical distribution of crack length also showed similar trend with the measured result.Therefore this polycrystal model could be an effective tool for the analysis of crack morphology and distribution in the delamination process of TBC.