This paper presents a numerical method for dynamic calculation of third order systems involving fractional operators. Using the Liouville-Rieman’s definition of fractional derivatives, a numerical algorithm is developed on base of the well-known Newmark integration method to calculate dynamic response of third order systems. Then, we apply this method to calculate linear vibrations of viscoelastic systems containing fractional derivatives. | Vietnam Journal of Mechanics, VAST, Vol. 34, No. 2 (2012), pp. 91 – 99 NUMERICAL CALCULATING LINEAR VIBRATIONS OF THIRD ORDER SYSTEMS INVOLVING FRACTIONAL OPERATORS Nguyen Van Khang1 , Tran Dinh Son2 , Bui Thi Thuy2 1 Hanoi University of Technology, Vietnam 2 Hanoi University of Mining and Geology, Vietnam Abstract. This paper presents a numerical method for dynamic calculation of third order systems involving fractional operators. Using the Liouville-Rieman’s definition of fractional derivatives, a numerical algorithm is developed on base of the well-known Newmark integration method to calculate dynamic response of third order systems. Then, we apply this method to calculate linear vibrations of viscoelastic systems containing fractional derivatives. Key words: Fractional order derivative, numerical method, vibration, third order system. 1. INTRODUCTION In 1959 Newmark presented a family of single-step integration methods for the solution of structural dynamic problems [1, 2]. During the past time Newmark’s method has been applied to the dynamic analysis of many practical engineering structures. It has been modified and improved by many other researchers such as Wilson, Hilber, Hughes and Taylor. However, these methods are only used for the system of second order equations. The concepts of fractional derivatives [3, 4, 5] appeared many years ago and are introduced by famous mathematicians like Riemann, Liouville, Gr¨ unwald, Letnikov, Caputo. The concept of fractional operators in engineering applications is now increasingly attractive in the formulations of the constitutive law for some viscoelastic materials. In [6, 7, 8] Shimizu and Zhang have used the Newmark integration method for calculating the vibrations of second order systems involving fractional derivatives. Many vibration problems in engineering lead the system of differential equations of third order. In this paper we present the using Newmark integration method for calculating vibrations of .
