A two dimensional thermo elasto-plastic finite clement stress analysis was performed to study stress distributions in functionally gradient material. The upper $ZrO_2$ surface is heated at 1200K until a steady state is established and cooled at 300K. The influences on the thermal stress distributions due to the difference of compositional gradient exponent p were investigated. In this study, we obtained the thermal stresses are low for p=1.