Simulation of fractionally damped mechanical systems by means of a Newmark-diffusive scheme

A Newmark-diffusive scheme is presented for the time-domain solution of dynamic systems containing fractional derivatives. This scheme combines a classical Newmark time-integration method used to solve second-order mechanical systems (obtained for example after finite element discretization), with a diffusive representation based on the transformation of the fractional operator into a diagonal system of linear differential equations, which can be seen as internal memory variables. The focus is given on the algorithm implementation into a finite element framework, the strategies for choosing diffusive parameters, and applications to beam structures with a fractional Zener model.