Fundamental solution and its validation by numerical inverse Laplace transformation and FEM for a damped Timoshenko beam subjected to impact and moving loads

• Published in: Journal of vibration and control 2019-02, Vol.25 (3), p.593-611
• Main Authors: Zhang, Xiayang, Liang, Haoquan, Zhao, Meijuan
• Format: Article
• Language: English
• Subjects: Boundary conditions; Case studies; Convergence; Damping; Finite element method; Impact loads; Laplace transforms; Load; Low speed; Mathematical models; Mode superposition method; Moving loads; Timoshenko beams
• ISSN: 1077-5463; 1741-2986
• DOI: 10.1177/1077546318790867

This paper, taking the clamped boundary condition as an example, develops Su and Ma's fundamental solutions of the dynamic responses of a Timoshenko beam subjected to impact load. Based on that, a further extension regarding the general moving load case is also established. Kelvin–Voigt damping, whether proportionally or nonproportionally damped, is incorporated into the model, making it more comprehensive than the model of Su and Ma. Numerical inverse Laplace transformation is introduced to obtain the time-domain solution, where Durbin's formula and the corresponding convergence criteria are utilized in numerical experiments. Further, the real modal superposition method is applied at an analytical level to validate the numerical results by applying a proportionally damped condition. Total comparisons are made between the methods by sufficient case studies. The dynamic responses with and without damping effect are computed with wider slenderness to verify the correctness and effectiveness of the numerical results. Furthermore, parametric studies regarding the damping coefficients are performed to explore the nonproportional damping effect. The results show that the structural damping has significant influences on the dynamic behaviors and is especially stronger at small slender ratios. As the damping decreases the inherent frequencies and excites the low-frequency modal components more actively, a resonant phenomenon appears in high slenderness case when the beam experiences a low-speed moving load. Additionally, the computations in the moving load case indicate that the algorithm convergence is preferable when the number of grids exceeds 1000.