Numerically Stable form of Green’s Function for a Free-Free Uniform Timoshenko Beam

Beam models are widely applied in civil engineering, transport, and industry because the beams are basic structural elements. When dealing with the high-order modes of beam in the context of applying the modal analysis method, the numerical instability issue affects the numeric simulation accuracy in many boundary conditions. There are two solutions in literature to overcome this shortcoming, namely refinement of the asymptotic form for the high order modes and reshaping the terms within the equation of the modes to eliminate the source of the numerical instability. In this paper, the numerical instability issue is signalled when the standard form of Green’s function, which includes hyperbolic functions, is applied to a free-free Timoshenko length-long beam. A new way is proposed based on new set of eigenfunctions, including an exponential function, to construct a new form of Green’s function. To this end, it starts from a new general form of Green’s function and the characteristic equation is obtained; then, based on the boundary condition, the Green’s function associated to the differential operator of the free-free Timoshenko beam is distilled. The numerical stability of the new form of the Green’s function is verified in a numerical application and the results are compared with those obtained by using the standard form of the Green’s function.