Finite element implementation of Huet-Sayegh and 2S2P1D models for analysis of asphalt pavement structures in time domain

• Published in: Road materials and pavement design Vol. 23; no. 1; pp. 22 - 46
• Main Authors: Tai Nguyen, H. T., Nguyen, Duy-Liem, Tran, Vu-Tu, Nguyen, Mai-Lan
• Format: Journal Article
• Language: English
• Published: Abingdon Taylor & Francis 02-01-2022; Lavoisier
• Subjects: 2S2P1D; Accuracy; Algorithms; Approximation; Asphalt pavements; Decay rate; decay rates estimation; Engineering Sciences; Exponential functions; Finite element method; Huet-Sayegh; Mathematical models; Numerical analysis; relaxation modulus; Rheological properties; Stress-strain relationships; sum of exponential functions; Time domain analysis; sum of exponential functions; Huet-Sayegh; 2S2P1D; relaxation modulus; finite element method; decay rates estimation
• ISSN: 1468-0629; 2164-7402
• DOI: 10.1080/14680629.2020.1809501

The Huet-Sayegh (HS) and 2S2P1D rheological models represent very well the linear viscoelastic behaviour of asphalt concrete and binder while using less parameters than other models. However, the finite element implementation in time domain of these two models encounters a great difficulty. In this paper, the author aims at developing a numerical implementation in time domain of the HS and 2S2P1D models. Firstly, the relaxation modulus is analytically approximated using Schapery and Park's method or exactly calculated using the numerical algorithm proposed by Heck. Secondly, the relaxation modulus is once again approximated by the sum of exponential functions for deriving a recursive formula of stress-strain relationship. Since the relaxation modulus is approximated twice, a special technique is used to estimate the appropriate decay rates and coefficients of the sum of exponential functions to ensure high accuracy of the obtained results. It is shown that the calculations based on the proposed formulation achieve a very good level of accuracy.