An Efficient Numerical Method for an Approximate Solution of the Beam Equation

Abstract: In this paper, we propose a horizontal type method of lines numerical scheme for the unsteady Euler-Bernoulli beam equation. The problem is initially reformulated as a first order system of initial value problems and a suitable one-step difference scheme is used for the highest order temporal derivative which leads to a system of steady beam equations. Then resulted family of steady problems is solved iteratively by the finite element method with Hermite cubic basis functions. This iterative procedure leads to approximations for both the solution of the unsteady problem and its derivatives. All these approximations are compared with the exact ones to illustrate the performance of the proposed method. Moreover, the optimization of the mesh parameters is discussed for both steady and unsteady problems by logarithmic scale plot.