Cracking patterns in ceramics coating under thermal stress are analyzed and simulated by "Analysis/FEM Strain Accommodation Method". This method combines the strain of a coated structure obtained from thermal-elastic analysis with the strain of substrate from Finite Element Method (FEM). The computer simulation of cracking patterns consists of random distribution of coating strength and interfacial strength, discrimination of cracking patterns, FEM analysis of strain in substrate, calculation of crack size from theoretical equation, crack initiation, propagation and connection. The method will be very effective not only to estimate the strengths of coating and interface but also to predict where and when the cracking or delamination occurs in coated structures with complex shapes under thermal stress.