A computer-implemented elastic response performance prediction method that predicts an elastic response performance expressing deformation behavior of a rubber product, the elastic response performance prediction method comprising: generating a model of the rubber product, operating an input-output terminal to set a constitutive equation as a framework condition of the rubber product, applying the framework condition to the model of the rubber product, and reconstructing model data of the rubber product, deforming the model under the framework conditions and predicting the elastic response performance of the rubber product by employing the constitutive equation that expresses temperature and strain dependence of strain energy in the rubber product expressed using a parameter representing intermolecular interaction, wherein the constitutive equation is the following Equation: Δ ⁢ A = U 1 − TS 1 + p V 1 − V 0 − U 0 − TS 0 wherein A represents Helmholz free energy, U 0 represents internal energy in a non-deformed state, U 1 represents internal energy in a deformed state, p represents pressure, V 0 represents volume in a non-deformed state, V 1 represents volume in a deformed state, T represents absolute temperature, S 0 represents entropy in a non-deformed state, and S 1 represents entropy in a deformed state, with each of the terms of Equation expressed by the following Equations to: U 1 − T ⋅ S 1 = e β ′ ⋅ κ κ cosh 2 β ′ I 1 ′ − 3 + 2 I 1 ′ − 3 ⋅ sinh 2 β ′ I 1 ′ − 3 e β ′ ⋅ κ cosh 2 β ′ I 1 ′ − 3 + 1 − N β ⋅ 1 2 I 1 + 3 100 ⁢ n 3 I 1 2 − 4 I 2 + 99 12250 ⁢ n 5 I 1 3 − 12 I 1 I 2 p V 1 − V 0 = B ⋅ V 1 − V 0 = B I 3 1 2 − 1 2 − 1 β ′ ln 1 + e β ′ ⋅ κ cosh 2 β ′ I 1 ′ − 3 − ln 1 + e β ′ ⋅ κ U 0 − T ⋅ S 0 = κ ⋅ e β ′ ⋅ κ e β ′ ⋅ κ + 1 + N β 3 2 + 45 100 ⁢ n + 2673 12250 n 2 wherein: I 1, I 2, and I 3 are expressed as functions of three extension ratios of deformation λ 1, λ 2 and λ 3 in xyz directions in three dimensional axes of rubber by I 1 = λ 1 2 + λ 2 2 + λ 3 2, I 2 = λ 1 2 · λ 2 2 + λ 2 2 · λ 3 2 + λ 3 2 · λ 1 2, and I 3 = λ 1 2 · λ 2 2 · λ 3 2; n represents a number of links between cross-linked points in a statistical molecule chain; κ expresses an intermolecular interaction energy coefficient; β = 1/RT and β' = 1/R (T - T g) wherein R is a gas constant and T g is a glass transition temperature; and I 1 ' is expressed using a local interaction function λ micro as a parameter expressing the intermolecular interaction by the following Equation: I 1 ′ = λ micro 2 λ 1 2 + λ 2 2 + λ 3 2 = λ micro 2 ⋅ I 1 generating a model of the rubber product, operating an input-output terminal to set a constitutive equation as a framework condition of the rubber product, applying the framework condition to the model of the rubber product, and reconstructing model data of the rubber product, deforming the model under the framework conditions and predicting the elastic response performance of the rubber product by employing the constitutive equation that expresses temperature and strain dependence of strain energy in the rubber product expressed using a parameter representing intermolecular interaction, wherein the constitutive equation is the following Equation: Δ ⁢ A = U 1 − TS 1 + p V 1 − V 0 − U 0 − TS 0 wherein A represents Helmholz free energy, U 0 represents internal energy in a non-deformed state, U 1 represents internal energy in a deformed state, p represents pressure, V 0 represents volume in a non-deformed state, V 1 represents volume in a deformed state, T represents absolute temperature, S 0 represents entropy in a non-deformed state, and S 1 represents entropy in a deformed state, with each of the terms of Equation expressed by the following Equations to: U 1 − T ⋅ S 1 = e β ′ ⋅ κ κ cosh 2 β ′ I 1 ′ − 3 + 2 I 1 ′ − 3 ⋅ sinh 2 β ′ I 1 ′ − 3 e β ′ ⋅ κ cosh 2 β ′ I 1 ′ − 3 + 1 − N β ⋅ 1 2 I 1 + 3 100 ⁢ n 3 I 1 2 − 4 I 2 + 99 12250 ⁢ n 5 I 1 3 − 12 I 1 I 2 p V 1 − V 0 = B ⋅ V 1 − V 0 = B I 3 1 2 − 1 2 − 1 β ′ ln 1 + e β ′ ⋅ κ cosh 2 β ′ I 1 ′ − 3 − ln 1 + e β ′ ⋅ κ U 0 − T ⋅ S 0 = κ ⋅ e β ′ ⋅ κ e β ′ ⋅ κ + 1 + N β 3 2 + 45 100 ⁢ n + 2673 12250 n 2 wherein: I 1, I 2, and I 3 are expressed as functions of three extension ratios of deformation λ 1, λ 2 and λ 3 in xyz directions in three dimensional axes of rubber by I 1 = λ 1 2 + λ 2 2 + λ 3 2, I 2 = λ 1 2 · λ 2 2 + λ 2 2 · λ 3 2 + λ 3 2 · λ 1 2, and I 3 = λ 1 2 · λ 2 2 · λ 3 2; n represents a number of links between cross-linked points in a statistical molecule chain; κ expresses an intermolecular interaction energy coefficient; β = 1/RT and β' = 1/R (T - T g) wherein R is a gas constant and T g is a glass transition temperature; and I 1 ' is expressed using a local interaction function λ micro as a parameter expressing the intermolecular interaction by the following Equation: I 1 ′ = λ micro 2 λ 1 2 + λ 2 2 + λ 3 2 = λ micro 2 ⋅ I 1