An analytical method for determining the thermal stress distributions in a 3-layered beam is developed, which is focused on the end effects. This method gives the stress distributions which satisfy the stress-free boundary condition at the end completely. For verification of the method, a numerical example which was introduced by other researchers is treated. The stress distributions agree with the results of other researchers. The results show that the show and peeling stresses at the interfaces are significant near the edge and become negligible in the interior region.