To achieve response analysis in a time region to a model to be analyzed having an energy transfer boundary as a wave motion boundary model.

In performing time history response analysis of the model to be analyzed, having an energy transfer boundary as the wave motion boundary model, a transfer boundary matrix of the energy transfer boundary including an external ground is calculated (48); impulse response in a time region of the energy transfer boundary which includes the external ground is calculated, based on the value of the transfer boundary matrix when vibration has each frequency (50), using as a mathematical expression for defining the impulse response, where the mathematical expression includes at least a rigidity term that depends on the displacement of the object and consists of a simultaneous component and a time delay component, and an attenuation term that depends on the speed of the object and consists of a simultaneous component and a time delay component; and the response analysis in the time region is applied to the model to be analyzed by using the calculated impulse response (54).