A numerical method to analyze three-dimensional dynamic thermal stresses in a semi-infinite plate is presented, in which instantaneous heating is applied to part of the plate surface. For the numerical analysis, the finite difference method based on integration along bicharacteristics is employed. From the numerical results, it is found that the propagation and reflection of the stress waves under instantaneous heating can be examined in detail by this method.