It is generally believed that the failure of TBCs is attributed to the spallation occurred in the ceramic coat. The spallation is closed linked with sinuate morphology factors, including its amplitude and period, at the TGO/bond coat interface. In this work, dependence of the residual stress distribution on the sinuate morphology in the TBCs has been studied by means of Finite Element Method (FEM) simulation for isothermally annealed specimens. The simulation results indicated that the maximum value of residual stress existed inside the TGO layer. It was also found that the maximum residual stress occurred at different points, near the TGO/bond coat interface at the peak of the sinuate interface, while near the TGO/ceramic coat interface at the valley, respectively. And the maximum residual stress increased with increasing the ratio of the amplitude to period in the sine morphology, which has been proved by the thermal cycle experimental results.