A coupled simulation of a structural-elastic phenomenon and a heat conduction phenomenon of a simulation target including plural particles is performed. Here, numerical calculation of a motion equation capable of being transformed into an equation of the same form as that of a heat conduction equation is performed with respect to a term of a spatial temperature distribution and a term of a derivative of temperature with respect to time, to perform a simulation of the heat conduction phenomenon of the simulation target.