In the core of the nuclear power plant of PWR, several cases of fuel failure by unknown causes have been experienced for various fuel types. From the common features of the failure pattern, failure lead time, flow conditions, and flow induced vibration characteristics in nuclear fuel bundles, it is deduced that the fretting wear failure of the fuel rod at the spacer grid position is due to the fluidelastic vibration. In the past, fluidelastic vibration was simulated by quasi -static semi-analytical model, so called the static model, which could not account for the interaction between the rods within a bundle. To overcome this defect and to provide for more flexibilities applicable to the fuel bundle, Tanaka's unsteady model was modified to accomodate the geometrical differences and governing parameter changes during the operations such as the number of rods, pitch to diameter ratio (P/D), spring force, damping coefficient, etc. The critical velocity was calculated by solving the governing equations with the MATLAB code. A comparison between the estimated critical velocity and the test result shows a good agreement. Finally, the level of decrease of the critical velocity due to the reduction in the spring force and reduced damping coefficient due to the radiation exposure is also estimated.