Summary
In this paper, a novel numerical integration scheme is proposed for fractional‐order viscoelastic analysis of hereditary‐aging structures. More precisely, the idea of aging is first introduced through a new phenomenological viscoelastic model characterized by variable‐order fractional operators. Then, the presented fractional‐order viscoelastic model is included in a variational formulation, conceived for any viscous kernel and discretized in time by employing a discontinuous Galerkin method. The accuracy of the resulting finite element (FE) scheme is analyzed through a model problem, whose exact solution is known; and the most significant variables affecting the solution quality, such as the number of Gaussian quadrature points and time subintervals, are then investigated in terms of error and computational cost. Moreover, the proposed FE integration scheme is applied to study the short‐ and long‐term behavior of concrete structures, which, due to the severe aging exhibited during their service life, represents one of the most challenging time‐dependent behavior to be investigated. Eventually, also the Euler implicit method, commonly used in commercial software, is compared.