Based on the principle of complementary energy, an analytical method is developed which focused on the end effects for determining thermal stress distributions in the clad beam. This method gives the stress distributions which completely satisfy the stress-free boundary condition at the edge. Numerical results shows that shear and peeling stress at the interface between the substrate and clad are significant near the edge and become negligible in the interior region. Even thought the relative location where the maximum or minimum stresses take place moves to interior as the length of the beam becomes smaller, the absolute location from the free end and the value of these stresses are the same in spite of the variation of the length of beam.