“…In 1921, Tymoshenko [1] optimized the Euler-Bernoulli beam model and the Rayleigh beam model and proposed the following hyperbolic system of two coupled wave equations ρ1ϕtt − k(ϕx + ψ)x = 0, ρ2ψtt − bψxx + k(ϕx + ψ) = 0, (1.1) which is called Timoshenko beam model, where ϕ and ψ are the deflection of the beam from its equilibrium position and the rotation of the neutral axis, respectively, ρ1 = ρA, ρ2 = ρI, b = EI and k = k GA are positive constants with ρ is the density, A is the cross-sectional area, I is the second moment of area of the cross-sectional area, E is the Young modulus of elasticity, G is the modulus of rigidity, k is the transverse shear factor. However, it was later discovered that the Timoshenko beam model admits two wave speeds k/ρ1 and b/ρ2, which contributes to a physical paradox called the second spectrum (see, for example, [2,3,4]). Based on these reasons, Elishakoff [5] proposed the following truncated version model by combining d'Alembert's principle for dynamic equilibrium from Timoshenko hypothesis,…”