The transfer matrix method for a distributed flexural vibrational system possesses some defects such as cancelling, attributable to the sum and difference of the hyperbolic and trigonometrical function appearing in the elements of the transfer matrices, and occurrence of numerical instability due to numerical imbalance among elements and multiplication of the matrices. The authors apply the concept of the transfer influence coefficient method, treated in the previous report, to a distributed flexural vibrational system, and formulate an algorithm for free vibrational analysis with a high speed and with high accuracy which succeeds in overcoming all these defects of the transfer matrix method. The validity of the present algorithm is demonstrated by a relatively simple example computing critical speeds of a rotating shaft, and is also compared with that of the transfer matrix method on a personal computer.
