Abstract
This paper presents a numerical analysis method for solving the shape identification problems of thermoelastic fields. The square error integral between the actual thermal deformation distributions and the prescribed thermal deformation distributions on the prescribed sub-boundaries is used as the objective functional. The shape gradient of the shape identification 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 domain optimization problems. A new numerical procedure for the shape identification was proposed. The validity of the proposed method was confirmed by the results of 3D numerical analysis.