Numerical analysis of viscoelastic pipe with variable fractional order model based on shifted Legendre polynomials algorithm