Analytical model based on the Bernoulli beam theory and strain compatibility conditions at the interfaces between the two layers have been developed to predict the distribution of thermal stresses within the multi-layered structure due to the mismatch of thermal expansion. The closed-form solution of thermal stresses related to the material properties and geometry were obtained. It is useful to provide a simple and efficient analytical model, so that the stress level in the layers can be accurately estimated. The analytical results are compared with finite element results. Good agreement demonstrates that the proposed approach is able to provide an efficient way for the calculation of the thermal stresses.