A shape optimization method for volume minimization of thermoelastic body is presented. In this study, the volume minimization problem is formulated, while mean compliance is constrained to less than a desired value, in the thermoelastic body. The shape gradient of the shape optimization problems was derived theoretically using the adjoint variable method, the Lagrange multiplier method and the formulae of the material derivative. Reshaping was accomplished using a traction method that was proposed as a solution to the shape optimization problems. The validity of the proposed method was confirmed based on the results of 2D numerical analysis.