A computationally efficient numerical algorithm for the transient response of high-speed elastic linkages

A new numerical algorithm, easily adaptable for computer simulation, is developed to approximate the transient response of a single degree-of-freedom vibrating system; governing differential equation is linear and second order with time-dependent and periodic coefficients. This is accomplished by first solving the classical linear single degree-of-freedom problem with constant coefficients. The system is excited by a periodic forcing function possessing a certain degree of smoothness. The integration terms in the solution are systematically expanded into two groups of terms: one consists of nonintegral terms while the other contains only integral terms. TJie final integral terms are bounded. For certain combinations of frequency and damping, within the sub-resonant frequency range, the relative size of the integral terms are demonstrated to be small. The algebraic expansion (non-integral) terms then approximate the solution. The solution to a single degree-of-freedom system with time-dependent and periodic parameters is made possible by discretizing the forcing period into a number of intervals and assuming the system parameters as constant over each interval. The numerical algorithm is tlien employed to solve an elastic linkage problem via modal superposition. Convergence of the solution is verified by refining the number of intervals of discretization.