A self-excited vibration evaluation method for evaluating self-excited vibration of a tube bundle arranged in a fluid so as to be supported by a support member , taking account of the friction damping effect applied to the tube bundle , the method comprising: for each of at least one eigenmode of the tube bundle , a time history response analysis step of performing time history response analysis of obtaining the vibration amplitude of the tube bundle as each value of a negative damping ratio being an input, while gradually increasing the value of the negative damping ratio corresponding to an excitation force of the fluid; a critical flow velocity calculation step of obtaining a critical negative damping ratio corresponding to a critical point at which the vibration amplitude of the tube bundle diverges in the time history response analysis, while the negative damping ratio is gradually increased, and calculating a critical flow velocity from the critical negative damping ratio; an input step of inputting an expected flow velocity of the fluid; and an evaluation step of evaluating the self-excited vibration of the tube bundle for each eigenmode by comparing the expected flow velocity of the fluid with the critical flow velocity, wherein the time history response analysis includes calculation which includes time-series simulation of vibration amplitude which occurs when an excitation force corresponding to the negative damping ratio is applied as an external force term to a vibration analysis model of the tube bundle , and wherein the vibration analysis model determines a magnitude of a friction force between the tube bundle and the support member , by assuming a distribution of a contact load acting between the tube bundle and the support member . for each of at least one eigenmode of the tube bundle , a time history response analysis step of performing time history response analysis of obtaining the vibration amplitude of the tube bundle as each value of a negative damping ratio being an input, while gradually increasing the value of the negative damping ratio corresponding to an excitation force of the fluid; a critical flow velocity calculation step of obtaining a critical negative damping ratio corresponding to a critical point at which the vibration amplitude of the tube bundle diverges in the time history response analysis, while the negative damping ratio is gradually increased, and calculating a critical flow velocity from the critical negative damping ratio; an input step of inputting an expected flow velocity of the fluid; and an evaluation step of evaluating the self-excited vibration of the tube bundle for each eigenmode by comparing the expected flow velocity of the fluid with the critical flow velocity, wherein the time history response analysis includes calculation which includes time-series simulation of vibration amplitude which occurs when an excitation force corresponding to the negative damping ratio is applied as an external force term to a vibration analysis model of the tube bundle , and wherein the vibration analysis model determines a magnitude of a friction force between the tube bundle and the support member , by assuming a distribution of a contact load acting between the tube bundle and the support member . <EOS>