A numerical integration approach for fractional‐order viscoelastic analysis of hereditary‐aging structures